-
Cold Quasar Investigation: Comparing Star Formation Rates to Black Hole Growth
Authors:
Sasha Mintz,
Brandon Coleman,
Allison Kirkpatrick
Abstract:
Cold quasars are a rare population of luminous, unobscured quasars associated with host galaxies that have a high star formation rate. We aimed to study the host galaxies of sixty four of these cold quasars in order to probe how the supermassive black holes and host galaxies were coevolving. We compiled data from the XXL survey and crossmatched with the VHS, WISE, and HerMES surveys to obtain mult…
▽ More
Cold quasars are a rare population of luminous, unobscured quasars associated with host galaxies that have a high star formation rate. We aimed to study the host galaxies of sixty four of these cold quasars in order to probe how the supermassive black holes and host galaxies were coevolving. We compiled data from the XXL survey and crossmatched with the VHS, WISE, and HerMES surveys to obtain multiwavelength photometry spanning the Xray to the infrared and including optical spectroscopy. From the data, we calculated the supermassive black hole masses using broad emission from the magnesium II and hydrogen beta lines. We compared this with the stellar mass of the entire galaxy and find that the black holes are significantly more massive than would be predicted by local relations, indicating that the majority of black hole growth precedes the bulk of the the stellar mass formation. In addition to this, we created a spectral energy distribution for each galaxy to calculate the star formation rate. We compared the star formation rate with the black hole accretion rate and find that the stellar mass is rapidly increasing at a relative rate faster than the black hole growth, supporting the picture where the black hole grows first.
△ Less
Submitted 11 March, 2024;
originally announced March 2024.
-
Favorable conditions for heavy element nucleosynthesis in rotating proto-magnetar winds
Authors:
Tejas Prasanna,
Matthew S. B. Coleman,
Todd A. Thompson
Abstract:
The neutrino-driven wind cooling phase of proto-neutron stars (PNSs) follows successful supernovae. Wind models without magnetic fields or rotation fail to achieve the necessary conditions for production of the third $r$-process peak, but robustly produce a weak $r$-process in neutron-rich winds. Using 2D magnetohydrodynamic simulations with magnetar-strength magnetic fields and rotation, we show…
▽ More
The neutrino-driven wind cooling phase of proto-neutron stars (PNSs) follows successful supernovae. Wind models without magnetic fields or rotation fail to achieve the necessary conditions for production of the third $r$-process peak, but robustly produce a weak $r$-process in neutron-rich winds. Using 2D magnetohydrodynamic simulations with magnetar-strength magnetic fields and rotation, we show that the PNS rotation rate significantly affects the thermodynamic conditions of the wind. We show that high entropy material is quasi-periodically ejected from the closed zone of the PNS magnetosphere with the required thermodynamic conditions to produce heavy elements. We show that maximum entropy $S$ of the material ejected depends systematically on the magnetar spin period $P_{\star}$ and scales as $S \propto P_{\star}^{-5/6}$ for sufficiently rapid rotation. We present results from simulations at a constant neutrino luminosity representative of $\sim 1-2$ s after the onset of cooling for $P_{\star}$ ranging from 5 ms to 200 ms and a few simulations with evolving neutrino luminosity where we follow the evolution of the magnetar wind until $10-14$ s after the onset of cooling. We estimate at magnetar polar magnetic field strength $B_0=3\times 10^{15}$ G, $10^{15}$ G, and $5\times 10^{14}$ G that neutron-rich magnetar winds can respectively produce at least $\sim 2-5\times 10^{-5}$ M$_{\odot}$, $\sim 3-4\times 10^{-6}$ M$_{\odot}$, and $\sim 2.5\times 10^{-8}$ M$_{\odot}$ of material with the required parameters for synthesis of the third $r-$process peak, within $1-2$ s, 10 s, and 14 s in that order after the onset of cooling. We show that proton-rich magnetar winds can have favorable conditions for production of $p-$nuclei, even at a modest $B_0=5\times 10^{14}$ G.
△ Less
Submitted 8 February, 2024;
originally announced February 2024.
-
A Theory for Neutron Star and Black Hole Kicks and Induced Spins
Authors:
Adam Burrows,
Tianshu Wang,
David Vartanyan,
Matthew S. B. Coleman
Abstract:
Using twenty long-term 3D core-collapse supernova simulations, we find that lower compactness progenitors that explode quasi-spherically due to the short delay to explosion experience smaller neutron star recoil kicks in the $\sim$100$-$200 km s$^{-1}$ range, while higher compactness progenitors that explode later and more aspherically leave neutron stars with kicks in the $\sim$300$-$1000 km s…
▽ More
Using twenty long-term 3D core-collapse supernova simulations, we find that lower compactness progenitors that explode quasi-spherically due to the short delay to explosion experience smaller neutron star recoil kicks in the $\sim$100$-$200 km s$^{-1}$ range, while higher compactness progenitors that explode later and more aspherically leave neutron stars with kicks in the $\sim$300$-$1000 km s$^{-1}$ range. In addition, we find that these two classes are correlated with the gravitational mass of the neutron star. This correlation suggests that the survival of binary neutron star systems may in part be due to their lower kick speeds. We also find a correlation of the kick with both the mass dipole of the ejecta and the explosion energy. Furthermore, one channel of black hole birth leaves masses of $\sim$10 $M_{\odot}$, is not accompanied by a neutrino-driven explosion, and experiences small kicks. A second is through a vigorous explosion that leaves behind a black hole with a mass of $\sim$3.0 $M_{\odot}$ kicked to high speeds. We find that the induced spins of nascent neutron stars range from seconds to $\sim$10 milliseconds, {but do not yet see a significant spin/kick correlation for pulsars.} We suggest that if an initial spin biases the explosion direction, a spin/kick correlation {would be} a common byproduct of the neutrino mechanism of core-collapse supernovae. Finally, the induced spin in explosive black hole formation is likely large and in the collapsar range. This new 3D model suite provides a greatly expanded perspective and appears to explain some observed pulsar properties by default.
△ Less
Submitted 31 January, 2024; v1 submitted 20 November, 2023;
originally announced November 2023.
-
The Accretion History of AGN: The Spectral Energy Distributions of X-ray Luminous AGN
Authors:
Connor Auge,
David Sanders,
Ezequiel Treister,
C. Megan Urry,
Allison Kirkpatrick,
Nico Cappelluti,
Tonima Tasnim Ananna,
Médéric Boquien,
Mislav Baloković,
Francesca Civano,
Brandon Coleman,
Aritra Ghosh,
Jeyhan Kartaltepe,
Michael Koss,
Stephanie LaMassa,
Stefano Marchesi,
Alessandro Peca,
Meredith Powell,
Benny Trakhtenbrot,
Tracey Jane Turner
Abstract:
Spectral energy distributions (SEDs) from X-ray to far-infrared (FIR) wavelengths are presented for a sample of 1246 X-ray luminous active galactic nuclei (AGN; $L_{0.5-10\rm{keV}}>10^{43}$ erg s$^{-1}$), with $z_{\rm{spec}}<1.2$, selected from Stripe 82X, COSMOS, and GOODS-N/S. The rest-frame SEDs show a wide spread ($\sim2.5$ dex) in the relative strengths of broad continuum features at X-ray, u…
▽ More
Spectral energy distributions (SEDs) from X-ray to far-infrared (FIR) wavelengths are presented for a sample of 1246 X-ray luminous active galactic nuclei (AGN; $L_{0.5-10\rm{keV}}>10^{43}$ erg s$^{-1}$), with $z_{\rm{spec}}<1.2$, selected from Stripe 82X, COSMOS, and GOODS-N/S. The rest-frame SEDs show a wide spread ($\sim2.5$ dex) in the relative strengths of broad continuum features at X-ray, ultraviolet (UV), mid-infrared (MIR), and FIR wavelengths. A linear correlation (log-log slope of 0.7$\pm0.04$) is found between $L_{\rm{MIR}}$ and $L_{\rm{X}}$. There is significant scatter in the relation between the $L_{\rm{UV}}$ and $L_{\rm{X}}$ due to heavy obscuration, however the most luminous and unobscured AGN show a linear correlation (log-log slope of 0.8$\pm0.06$) in the relation above this scatter. The relation between $L_{\rm{FIR}}$ and $L_{\rm{X}}$ is predominantly flat, but with decreasing dispersion at $L_{\rm{X}}>10^{44}$ erg s$^{-1}$. The ratio between the "galaxy subtracted" bolometric luminosity and the intrinsic $L_{\rm{X}}$ increases from a factor of $\sim$$10-70$ from log $L_{\rm{bol}}/{\rm(erg\; s}^{-1})=44.5-46.5$. Characteristic SED shapes have been determined by grouping AGN based on relative strengths of the UV and MIR emission. The average $L_{1μ\rm{m}}$ is constant for the majority of these SED shapes, while AGN with the strongest UV and MIR emission have elevated $L_{1μ\rm{m}}$, consistent with the AGN emission dominating their SEDs at optical and NIR wavelengths. A strong correlation is found between the SED shape and both the $L_{\rm{X}}$ and $L_{\rm{bol}}$, such that $L_{\rm{bol}}/L_{\rm{X}}=20.4\pm1.8$, independent of the SED shape. This is consistent with an evolutionary scenario of increasing $L_{\rm{bol}}$ with decreasing obscuration as the AGN blows away circumnuclear gas.
△ Less
Submitted 18 August, 2023;
originally announced August 2023.
-
AGN's Deadness Over Cosmic Time: UVJ Diagrams of X-Ray-Selected AGN
Authors:
Rachel Cionitti,
Brandon Coleman,
Allison Kirkpatrick,
Greg Troiani
Abstract:
Active Galactic Nuclei (AGN) are intensely accreting supermassive black holes at the centers of massive galaxies. Though these objects occupy little spatial extent of the galaxy itself, they are thought to have far reaching affects, impacting the galaxy's star formation, and possibly it's lifespan until it becomes 'red and dead'. Typical galaxies demonstrate that, over cosmic time, they tend to se…
▽ More
Active Galactic Nuclei (AGN) are intensely accreting supermassive black holes at the centers of massive galaxies. Though these objects occupy little spatial extent of the galaxy itself, they are thought to have far reaching affects, impacting the galaxy's star formation, and possibly it's lifespan until it becomes 'red and dead'. Typical galaxies demonstrate that, over cosmic time, they tend to separate into a bimodal distribution of 'red and dead' or blue and star forming. We examine whether active galaxies evolve over cosmic time in a similar way, and whether this can reveal anything about the complexities of the relationship between an AGN and the host galaxy. We use the Stripe82X survey to identify 3940 X-ray AGN spanning z=0-2.5, and we measure the rest-frame UVJ colors of each galaxy. We classify AGN as star-forming or quiescent based on their location in a UVJ color diagram. We find that there is not a clear bimodal distribution between AGN in star forming and quiescent galaxies. Furthermore, the most luminous X-ray sources tend to lie in the star forming region, which may indicate a correlation between central engine activity and increased rates of star formation.
△ Less
Submitted 21 July, 2023;
originally announced July 2023.
-
The early evolution of magnetar rotation -- II. Rapidly rotating magnetars: Implications for Gamma-Ray Bursts and Super Luminous Supernovae
Authors:
Tejas Prasanna,
Matthew S. B. Coleman,
Matthias J. Raives,
Todd A. Thompson
Abstract:
Rapidly rotating magnetars have been associated with gamma-ray bursts (GRBs) and super-luminous supernovae (SLSNe). Using a suite of 2D magnetohydrodynamic simulations at fixed neutrino luminosity and a couple of evolutionary models with evolving neutrino luminosity and magnetar spin period, we show that magnetars are viable central engines for powering GRBs and SLSNe. We also present analytic est…
▽ More
Rapidly rotating magnetars have been associated with gamma-ray bursts (GRBs) and super-luminous supernovae (SLSNe). Using a suite of 2D magnetohydrodynamic simulations at fixed neutrino luminosity and a couple of evolutionary models with evolving neutrino luminosity and magnetar spin period, we show that magnetars are viable central engines for powering GRBs and SLSNe. We also present analytic estimates of the energy outflow rate from the proto-neutron star (PNS) as a function of polar magnetic field strength $B_0$, PNS angular velocity $Ω_{\star}$, PNS radius $R_{\star}$ and mass outflow rate $\dot{M}$. We show that rapidly rotating magnetars with spin periods $P_{\star}\lesssim 4$ ms and polar magnetic field strength $B_0\gtrsim 10^{15}$ G can release $10^{50}-5\times 10^{51}$ ergs of energy during the first $\sim2$ s of the cooling phase. Based on this result, it is plausible that sustained energy injection by magnetars through the relativistic wind phase can power GRBs. We also show that magnetars with moderate field strengths of $B_0\lesssim 5\times 10^{14}$ G do not release a large fraction of their rotational kinetic energy during the cooling phase and hence, are not likely to power GRBs. Although we cannot simulate to times greater than $\sim 3-5$ s after a supernova, we can hypothesize that moderate field strength magnetars can brighten the supernova light curves by releasing their rotational kinetic energy via magnetic dipole radiation on timescales of days to weeks, since these do not expend most of their rotational kinetic energy during the early cooling phase.
△ Less
Submitted 21 October, 2023; v1 submitted 25 May, 2023;
originally announced May 2023.
-
The James Webb Space Telescope Mission
Authors:
Jonathan P. Gardner,
John C. Mather,
Randy Abbott,
James S. Abell,
Mark Abernathy,
Faith E. Abney,
John G. Abraham,
Roberto Abraham,
Yasin M. Abul-Huda,
Scott Acton,
Cynthia K. Adams,
Evan Adams,
David S. Adler,
Maarten Adriaensen,
Jonathan Albert Aguilar,
Mansoor Ahmed,
Nasif S. Ahmed,
Tanjira Ahmed,
Rüdeger Albat,
Loïc Albert,
Stacey Alberts,
David Aldridge,
Mary Marsha Allen,
Shaune S. Allen,
Martin Altenburg
, et al. (983 additional authors not shown)
Abstract:
Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astrono…
▽ More
Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 years, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit.
△ Less
Submitted 10 April, 2023;
originally announced April 2023.
-
The Gravitational-Wave Signature of Core-Collapse Supernovae
Authors:
David Vartanyan,
Adam Burrows,
Tianshu Wang,
Matthew S. B. Coleman,
Christopher J. White
Abstract:
We calculate the gravitational-wave (GW) signatures of detailed 3D core-collapse supernova simulations spanning a range of massive stars. Most of the simulations are carried out to times late enough to capture more than 95% of the total GW emission. We find that the f/g-mode and f-mode of proto-neutron star oscillations carry away most of the GW power. The f-mode frequency inexorably rises as the…
▽ More
We calculate the gravitational-wave (GW) signatures of detailed 3D core-collapse supernova simulations spanning a range of massive stars. Most of the simulations are carried out to times late enough to capture more than 95% of the total GW emission. We find that the f/g-mode and f-mode of proto-neutron star oscillations carry away most of the GW power. The f-mode frequency inexorably rises as the proto-neutron star (PNS) core shrinks. We demonstrate that the GW emission is excited mostly by accretion plumes onto the PNS that energize modal oscillations and also high-frequency (``haze") emission correlated with the phase of violent accretion. The duration of the major phase of emission varies with exploding progenitor and there is a strong correlation between the total GW energy radiated and the compactness of the progenitor. Moreover, the total GW emissions vary by as much as three orders of magnitude from star to star. For black-hole formation, the GW signal tapers off slowly and does not manifest the haze seen for the exploding models. For such failed models, we also witness the emergence of a spiral shock motion that modulates the GW emission at a frequency near $\sim$100 Hertz that slowly increases as the stalled shock sinks. We find significant angular anisotropy of both the high- and low-frequency (memory) GW emissions, though the latter have very little power.
△ Less
Submitted 30 May, 2023; v1 submitted 6 February, 2023;
originally announced February 2023.
-
Magnetized Rotating Isothermal Winds
Authors:
Matthias J. Raives,
Matthew S. B. Coleman,
Todd A. Thompson
Abstract:
We consider the general problem of a Parker-type non-relativistic isothermal wind from a rotating and magnetic star. Using the magnetohydrodynamics (MHD) code athena++, we construct an array of simulations in the stellar rotation rate $Ω_\ast$ and the isothermal sound speed $c_T$, and calculate the mass, angular momentum, and energy loss rates across this parameter space. We also briefly consider…
▽ More
We consider the general problem of a Parker-type non-relativistic isothermal wind from a rotating and magnetic star. Using the magnetohydrodynamics (MHD) code athena++, we construct an array of simulations in the stellar rotation rate $Ω_\ast$ and the isothermal sound speed $c_T$, and calculate the mass, angular momentum, and energy loss rates across this parameter space. We also briefly consider the three dimensional case, with misaligned magnetic and rotation axes.
We discuss applications of our results to the spindown of normal stars, highly-irradiated exoplanets, and to nascent highly-magnetic and rapidly-rotating neutron stars born in massive star core collapse.8
△ Less
Submitted 30 October, 2023; v1 submitted 10 February, 2023;
originally announced February 2023.
-
Kicks and Induced Spins of Neutron Stars at Birth
Authors:
Matthew S. B. Coleman,
Adam Burrows
Abstract:
Using simulations of non-rotating supernova progenitors, we explore the kicks imparted to and the spins induced in the compact objects birthed in core collapse. We find that the recoil due to neutrino emissions can be a factor affecting core recoil, comparable to and at times larger than the corresponding kick due to matter recoil. This result would necessitate a revision of the general model of t…
▽ More
Using simulations of non-rotating supernova progenitors, we explore the kicks imparted to and the spins induced in the compact objects birthed in core collapse. We find that the recoil due to neutrino emissions can be a factor affecting core recoil, comparable to and at times larger than the corresponding kick due to matter recoil. This result would necessitate a revision of the general model of the origin of pulsar proper motions. In addition, we find that the sign of the net neutrino momentum can be opposite to the sign of the corresponding matter recoil. As a result, at times the pulsar recoil and ejecta can be in the same direction. Moreover, our results suggest that the duration of the dipole in the neutrino emissions can be shorter than the duration of the radiation of the neutron-star binding energy. This allows a larger dipole asymmetry to arise, but for a shorter time, resulting in kicks in the observed pulsar range. Furthermore, we find that the spin induced by the aspherical accretion of matter can leave the residues of collapse with spin periods comparable to those inferred for radio pulsars and that there seems to be a slight anti-correlation between the direction of the induced spin and the net kick direction. This could explain such a correlation among observed radio pulsars. Finally, we find that the kicks imparted to black holes are due to the neutrino recoil alone, resulting in birth kicks $\le$100 km s$^{-1}$ most of the time.
△ Less
Submitted 12 September, 2022; v1 submitted 6 September, 2022;
originally announced September 2022.
-
The Early Evolution of Magnetar Rotation I: Slowly Rotating "Normal" Magnetars
Authors:
Tejas Prasanna,
Matthew S. B. Coleman,
Matthias J. Raives,
Todd A. Thompson
Abstract:
In the seconds following their formation in core-collapse supernovae, "proto"-magnetars drive neutrino-heated magneto-centrifugal winds. Using a suite of two-dimensional axisymmetric MHD simulations, we show that relatively slowly rotating magnetars with initial spin periods of $P_{\star0}=50-500$ ms spin down rapidly during the neutrino Kelvin-Helmholtz cooling epoch. These initial spin periods a…
▽ More
In the seconds following their formation in core-collapse supernovae, "proto"-magnetars drive neutrino-heated magneto-centrifugal winds. Using a suite of two-dimensional axisymmetric MHD simulations, we show that relatively slowly rotating magnetars with initial spin periods of $P_{\star0}=50-500$ ms spin down rapidly during the neutrino Kelvin-Helmholtz cooling epoch. These initial spin periods are representative of those inferred for normal Galactic pulsars, and much slower than those invoked for gamma-ray bursts and super-luminous supernovae. Since the flow is non-relativistic at early times, and because the Alfvén radius is much larger than the proto-magnetar radius, spindown is millions of times more efficient than the typically-used dipole formula. Quasi-periodic plasmoid ejections from the closed zone enhance spindown. For polar magnetic field strengths $B_0\gtrsim5\times10^{14}$ G, the spindown timescale can be shorter than than the Kelvin-Helmholtz timescale. For $B_0\gtrsim10^{15}$ G, it is of order seconds in early phases. We compute the spin evolution for cooling proto-magnetars as a function of $B_0$, $P_{\star0}$, and mass ($M$). Proto-magnetars born with $B_0$ greater than $\simeq1.3\times10^{15}\,{\rm\,G}\,(P_{\star0}/{400\,\rm\,ms})^{-1.4}(M/1.4\,{\rm M}_\odot)^{2.2}$ spin down to periods $> 1$ s in just the first few seconds of evolution, well before the end of the cooling epoch and the onset of classic dipole spindown. Spindown is more efficient for lower $M$ and for larger $P_{\star0}$. We discuss the implications for observed magnetars, including the discrepancy between their characteristic ages and supernova remnant ages. Finally, we speculate on the origin of 1E 161348-5055 in the remnant RCW 103, and the potential for other ultra-slowly rotating magnetars.
△ Less
Submitted 18 August, 2022;
originally announced August 2022.
-
The Science Performance of JWST as Characterized in Commissioning
Authors:
Jane Rigby,
Marshall Perrin,
Michael McElwain,
Randy Kimble,
Scott Friedman,
Matt Lallo,
René Doyon,
Lee Feinberg,
Pierre Ferruit,
Alistair Glasse,
Marcia Rieke,
George Rieke,
Gillian Wright,
Chris Willott,
Knicole Colon,
Stefanie Milam,
Susan Neff,
Christopher Stark,
Jeff Valenti,
Jim Abell,
Faith Abney,
Yasin Abul-Huda,
D. Scott Acton,
Evan Adams,
David Adler
, et al. (601 additional authors not shown)
Abstract:
This paper characterizes the actual science performance of the James Webb Space Telescope (JWST), as determined from the six month commissioning period. We summarize the performance of the spacecraft, telescope, science instruments, and ground system, with an emphasis on differences from pre-launch expectations. Commissioning has made clear that JWST is fully capable of achieving the discoveries f…
▽ More
This paper characterizes the actual science performance of the James Webb Space Telescope (JWST), as determined from the six month commissioning period. We summarize the performance of the spacecraft, telescope, science instruments, and ground system, with an emphasis on differences from pre-launch expectations. Commissioning has made clear that JWST is fully capable of achieving the discoveries for which it was built. Moreover, almost across the board, the science performance of JWST is better than expected; in most cases, JWST will go deeper faster than expected. The telescope and instrument suite have demonstrated the sensitivity, stability, image quality, and spectral range that are necessary to transform our understanding of the cosmos through observations spanning from near-earth asteroids to the most distant galaxies.
△ Less
Submitted 10 April, 2023; v1 submitted 12 July, 2022;
originally announced July 2022.
-
The Essential Character of the Neutrino Mechanism of Core-Collapse Supernova Explosions
Authors:
Tianshu Wang,
David Vartanyan,
Adam Burrows,
Matthew S. B. Coleman
Abstract:
Calibrating with detailed 2D core-collapse supernova simulations, we derive a simple core-collapse supernova explosion condition based solely upon the terminal density profiles of state-of-the-art stellar evolution calculations of the progenitor massive stars. This condition captures the vast majority of the behavior of the one hundred 2D state-of-the-art models we performed to gauge its usefulnes…
▽ More
Calibrating with detailed 2D core-collapse supernova simulations, we derive a simple core-collapse supernova explosion condition based solely upon the terminal density profiles of state-of-the-art stellar evolution calculations of the progenitor massive stars. This condition captures the vast majority of the behavior of the one hundred 2D state-of-the-art models we performed to gauge its usefulness. The goal is to predict, without resort to detailed simulation, the explodability of a given massive star. We find that the simple maximum fractional ram pressure jump discriminant we define works well ~90% of the time and we speculate on the origin of the few false positives and false negatives we witness. The maximum ram pressure jump generally occurs at the time of accretion of the silicon/oxygen interface, but not always. Our results depend upon the fidelity with which the current implementation of our code Fornax adheres to Nature and issues concerning the neutrino-matter interaction, the nuclear equation of state, the possible effects of neutrino oscillations, grid resolution, the possible role of rotation and magnetic fields, and the accuracy of the numerical algorithms employed remain to be resolved. Nevertheless, the explodability condition we obtain is simple to implement, shows promise that it might be further generalized while still employing data from only the unstable Chandrasekhar progenitors, and is a more credible and robust simple explosion predictor than can currently be found in the literature.
△ Less
Submitted 23 September, 2022; v1 submitted 5 July, 2022;
originally announced July 2022.
-
Accretion History of AGN: Estimating the Host Galaxy Properties in X-ray Luminous AGN from z=0-3
Authors:
Brandon Coleman,
Allison Kirkpatrick,
Kevin C. Cooke,
Eilat Glikman,
Stephanie La Massa,
Stefano Marchesi,
Alessandro Peca,
Ezequiel Treister,
Connor Auge,
C. Megan Urry,
Dave Sanders,
Tracey Jane Turner,
Tonima Tasnim Ananna
Abstract:
We aim to determine the intrinsic far-Infrared (far-IR) emission of X-ray-luminous quasars over cosmic time. Using a 16 deg^2 region of the Stripe 82 field surveyed by XMM-Newton and Herschel Space Observatory, we identify 2905 X-ray luminous (LX > 10^42 erg/s) Active Galactic Nuclei (AGN) in the range z ~ 0-3. The IR is necessary to constrain host galaxy properties such as star formation rate (SF…
▽ More
We aim to determine the intrinsic far-Infrared (far-IR) emission of X-ray-luminous quasars over cosmic time. Using a 16 deg^2 region of the Stripe 82 field surveyed by XMM-Newton and Herschel Space Observatory, we identify 2905 X-ray luminous (LX > 10^42 erg/s) Active Galactic Nuclei (AGN) in the range z ~ 0-3. The IR is necessary to constrain host galaxy properties such as star formation rate (SFR) and gas mass. However, only 10% of our AGN are detected both in the X-ray and IR. Because 90% of the sample is undetected in the far-IR by Herschel, we explore the mean IR emission of these undetected sources by stacking their Herschel/SPIRE images in bins of X-ray luminosity and redshift. We create stacked spectral energy distributions from the optical to the far-IR, and estimate the median star formation rate, dust mass, stellar mass, and infrared luminosity using a fitting routine. We find that the stacked sources on average have similar SFR/L_bol ratios as IR detected sources. The majority of our sources fall on or above the main sequence line suggesting that X-ray selection alone does not predict the location of a galaxy on the main sequence. We also find that the gas depletion timescales of our AGN are similar to those of dusty star forming galaxies. This suggests that X-ray selected AGN host high star formation and that there are no signs of declining star formation.
△ Less
Submitted 13 June, 2022;
originally announced June 2022.
-
Boundary Layers of Accretion Disks: Wave-Driven Transport and Disk Evolution
Authors:
Matthew S. B. Coleman,
Roman R. Rafikov,
Alexander A. Philippov
Abstract:
Astrophysical objects possessing a material surface (white dwarfs, young stars, etc.) may accrete gas from the disc through the so-called surface boundary layer (BL), in which the angular velocity of the accreting gas experiences a sharp drop. Acoustic waves excited by the supersonic shear in the BL play an important role in mediating the angular momentum and mass transport through that region. He…
▽ More
Astrophysical objects possessing a material surface (white dwarfs, young stars, etc.) may accrete gas from the disc through the so-called surface boundary layer (BL), in which the angular velocity of the accreting gas experiences a sharp drop. Acoustic waves excited by the supersonic shear in the BL play an important role in mediating the angular momentum and mass transport through that region. Here we examine the characteristics of the angular momentum transport produced by the different types of wave modes emerging in the inner disc, using the results of a large suite of hydrodynamic simulations of the BLs. We provide a comparative analysis of the transport properties of different modes across the range of relevant disc parameters. In particular, we identify the types of modes which are responsible for the mass accretion onto the central object. We find the correlated perturbations of surface density and radial velocity to provide an important contribution to the mass accretion rate. Although the wave-driven transport is intrinsically non-local, we do observe a clear correlation between the angular momentum flux injected into the disc by the waves and the mass accretion rate through the BL. We find the efficiency of angular momentum transport (normalized by thermal pressure) to be a weak function of the flow Mach number. We also quantify the wave-driven evolution of the inner disc, in particular the modification of the angular frequency profile in the disc. Our results pave the way for understanding wave-mediated transport in future three-dimensional, magnetohydrodynamic studies of the BLs.
△ Less
Submitted 10 March, 2022; v1 submitted 4 November, 2021;
originally announced November 2021.
-
On the Origin of Pulsar and Magnetar Magnetic Fields
Authors:
Christopher J. White,
Adam Burrows,
Matthew S. B. Coleman,
David Vartanyan
Abstract:
In order to address the generation of neutron star magnetic fields, with particular focus on the dichotomy between magnetars and radio pulsars, we consider the properties of dynamos as inferred from other astrophysical systems. With sufficiently low (modified) Rossby number, convective dynamos are known to produce dipole-dominated fields whose strength scales with convective flux, and we argue tha…
▽ More
In order to address the generation of neutron star magnetic fields, with particular focus on the dichotomy between magnetars and radio pulsars, we consider the properties of dynamos as inferred from other astrophysical systems. With sufficiently low (modified) Rossby number, convective dynamos are known to produce dipole-dominated fields whose strength scales with convective flux, and we argue that these expectations should apply to the convective proto-neutron stars at the centers of core-collapse supernovae. We analyze a suite of three-dimensional simulations of core collapse, featuring a realistic equation of state and full neutrino transport, in this context. All our progenitor models, ranging from 9 solar masses to 25 solar masses, including one with initial rotation, have sufficiently vigorous proto-neutron-star convection to generate dipole fields of order ~10^15 gauss, if the modified Rossby number resides in the critical range. Thus, the magnetar/radio pulsar dichotomy may arise naturally in part from the distribution of core rotation rates in massive stars.
△ Less
Submitted 19 December, 2021; v1 submitted 2 November, 2021;
originally announced November 2021.
-
The Physical Effects of Progenitor Rotation: Comparing Two Long-Duration 3D Core-Collapse Supernova Simulations
Authors:
Matthew S. B. Coleman,
Adam Burrows,
Christopher J. White
Abstract:
We analyse and determine the effects of modest progenitor rotation in the context of core-collapse supernovae by comparing two separate long-duration three-dimensional simulations of 9 M$_{\odot}$ progenitors, one rotating with an initial spin period of $\sim$60 seconds and the other non-rotating. We determine that both models explode early, though the rotating model explodes a bit earlier. Despit…
▽ More
We analyse and determine the effects of modest progenitor rotation in the context of core-collapse supernovae by comparing two separate long-duration three-dimensional simulations of 9 M$_{\odot}$ progenitors, one rotating with an initial spin period of $\sim$60 seconds and the other non-rotating. We determine that both models explode early, though the rotating model explodes a bit earlier. Despite this difference, the asymptotic explosion energies ($\sim$10$^{50}$ ergs) and residual neutron star baryon masses ($\sim$1.3 M$_{\odot}$) are similar. We find that the proto-neutron star (PNS) core can deleptonize and cool significantly more quickly. Soon into the evolution of the rotating model, we witness more vigorous and extended PNS core convection that early in its evolution envelopes the entire inner sphere, not just a shell. Moreover, we see a corresponding excursion in both the $ν_e$ luminosity and gravitational-wave strain that may be diagnostic of this observed dramatic phenomenon. In addition, after bounce the innermost region of the rotating model seems to execute meridional circulation. The rotationally-induced growth of the convective PNS region may facilitate the growth of core B-fields by the dynamo mechanism by facilitating the achievement of the critical Rossby number condition for substantial growth of a dipole field, obviating the need for rapid rotation rates to create dipole fields of significance. The next step is to explore the progenitor-mass and spin dependencies across the progenitor continuum of the supernova explosion, dynamics, and evolution of PNS convection and its potential role in the generation of magnetar and pulsar magnetic fields.
△ Less
Submitted 3 May, 2022; v1 submitted 29 October, 2021;
originally announced November 2021.
-
The Collapse and Three-Dimensional Explosion of Three-Dimensional, vis à vis One-Dimensional, Massive-star Supernova Progenitor Models
Authors:
David Vartanyan,
Matthew S. B. Coleman,
Adam Burrows
Abstract:
The explosion outcome and diagnostics of core-collapse supernovae depend sensitively on the nature of the stellar progenitor, but most studies to date have focused exclusively on one-dimensional, spherically-symmetric massive star progenitors. We present some of the first core-collapse supernovae simulations of three-dimensional massive star supernovae progenitors, a 12.5- and a 15-M$_{\odot}$ mod…
▽ More
The explosion outcome and diagnostics of core-collapse supernovae depend sensitively on the nature of the stellar progenitor, but most studies to date have focused exclusively on one-dimensional, spherically-symmetric massive star progenitors. We present some of the first core-collapse supernovae simulations of three-dimensional massive star supernovae progenitors, a 12.5- and a 15-M$_{\odot}$ model, evolved in three-dimensions from collapse to bounce through explosion with the radiation-hydrodynamic code F{\sc{ornax}}. We compare the results using those starting from three-dimensional progenitors to three-dimensional simulations of spherically-symmetric, one-dimensional progenitors of the same mass. We find that the models evolved in three dimensions during the final stages of massive star evolution are more prone to explosion. The turbulence arising in these multi-dimensional initial models serve as seed turbulence that promotes shock revival. Detection of gravitational waves and neutrinos signals could reveal signatures of pre-bounce turbulence.
△ Less
Submitted 12 January, 2022; v1 submitted 22 September, 2021;
originally announced September 2021.
-
Boundary Layers of Accretion Disks: Discovery of Vortex-Driven Modes and Other Waves
Authors:
Matthew S. B. Coleman,
Roman R. Rafikov,
Alexander A. Philippov
Abstract:
Disk accretion onto weakly magnetized objects possessing a material surface must proceed via the so-called boundary layer (BL) - a region at the inner edge of the disk, in which the velocity of accreting material abruptly decreases from its Keplerian value. Supersonic shear arising in the BL is known to be conducive to excitation of acoustic waves that propagate into both the accretor and the disk…
▽ More
Disk accretion onto weakly magnetized objects possessing a material surface must proceed via the so-called boundary layer (BL) - a region at the inner edge of the disk, in which the velocity of accreting material abruptly decreases from its Keplerian value. Supersonic shear arising in the BL is known to be conducive to excitation of acoustic waves that propagate into both the accretor and the disk, enabling angular momentum and mass transport across the BL. We carry out a numerical exploration of different wave modes that operate near the BL, focusing on their morphological characteristics in the innermost parts of accretion disk. Using a large suite of simulations covering a broad range of Mach numbers (of the supersonic shear flow in the BL), we provide accurate characterization of the different types of modes, verifying their properties against analytical results, when available. We discover new types of modes, in particular, global spiral density waves launched by vortices forming in the disk near the BL as a result of the Rossby wave instability; this instability is triggered by the vortensity production in that region caused by the nonlinear damping of acoustic waves. Azimuthal wavenumbers of the dominant modes that we observe appear to increase monotonically with the Mach number of the runs, but a particular mix of modes found in a simulation is mildly stochastic. Our results provide a basis for better understanding of the angular momentum and mass transport across the BL as well as the emission variability in accreting objects.
△ Less
Submitted 11 October, 2021; v1 submitted 22 March, 2021;
originally announced March 2021.
-
Non-Parametric Density Reconstruction of the Galactic Bulge Area using Red Clump Stars in the VVV Survey
Authors:
Dylan Paterson,
Brenna Coleman,
Chris Gordon
Abstract:
Studies of the red clump giant population in the inner Milky Way suggest the Galactic bulge/bar has a boxy/peanut/X-shaped structure as predicted by its formation via a disc buckling instability. We used a non-parametric method of estimating the Galactic bulge morphology that is based on maximum entropy regularisation. This enabled us to extract the three-dimensional distribution of the red gian…
▽ More
Studies of the red clump giant population in the inner Milky Way suggest the Galactic bulge/bar has a boxy/peanut/X-shaped structure as predicted by its formation via a disc buckling instability. We used a non-parametric method of estimating the Galactic bulge morphology that is based on maximum entropy regularisation. This enabled us to extract the three-dimensional distribution of the red giant stars in the bulge from deep photometric catalogues of the VISTA Variables in the Via Lactea (VVV) survey. Our high-resolution reconstruction confirms the well-known boxy/peanut/X-shaped structure of the bulge. We also find spiral arm structures that extend to around three kpc in front of and behind the bulge and are on different sides of the bulge major axis. We show that the detection of these structures is robust to the uncertainties in the luminosity function.
△ Less
Submitted 13 September, 2020; v1 submitted 12 November, 2019;
originally announced November 2019.
-
Maximum Entropy Estimation of the Galactic Bulge Morphology via the VVV Red Clump
Authors:
Brenna Coleman,
Dylan Paterson,
Chris Gordon,
Oscar Macias,
Harrison Ploeg
Abstract:
The abundance and narrow magnitude dispersion of Red Clump (RC) stars make them a popular candidate for mapping the morphology of the bulge region of the Milky Way. Using an estimate of the RC's intrinsic luminosity function, we extracted the three-dimensional density distribution of the RC from deep photometric catalogues of the VISTA Variables in the Via Lactea (VVV) survey. We used maximum en…
▽ More
The abundance and narrow magnitude dispersion of Red Clump (RC) stars make them a popular candidate for mapping the morphology of the bulge region of the Milky Way. Using an estimate of the RC's intrinsic luminosity function, we extracted the three-dimensional density distribution of the RC from deep photometric catalogues of the VISTA Variables in the Via Lactea (VVV) survey. We used maximum entropy based deconvolution to extract the spatial distribution of the bulge from Ks-band star counts. We obtained our extrapolated non-parametric model of the bulge over the inner 40 by 40 degrees squared region of the Galactic centre. Our reconstruction also naturally matches onto a parametric fit to the bulge outside the VVV region and inpaints overcrowded and high extinction regions. We found a range of bulge properties consistent with other recent investigations based on the VVV data. In particular, we estimated the bulge mass to be in the range 13 to 17 billion solar masses, the X-component to be between 18% and 25% of the bulge mass, and the bulge angle with respect to the Sun-Galactic centre line to be between 18 and 32 degrees. Studies of the Fermi Large Area Telescope (LAT) gamma-ray Galactic centre excess suggests that the excess may be traced by Galactic bulge distributed sources. We applied our deconvolved density in a template fitting analysis of this Fermi-LAT GeV excess and found an improvement in the fit compared to previous parametric based templates.
△ Less
Submitted 10 December, 2020; v1 submitted 12 November, 2019;
originally announced November 2019.
-
An Extension of the Athena++ Framework for General Equations of State
Authors:
Matthew S. B. Coleman
Abstract:
We present modifications to the Athena++ framework to enable use of general equations of state (EOS). Part of our motivation for doing so is to model transient astrophysics phenomena, as these types of events are often not well approximated by an ideal gas. This necessitated changes to the Riemann solvers implemented in Athena++. We discuss the adjustments made to the HLLC, and HLLD solvers and EO…
▽ More
We present modifications to the Athena++ framework to enable use of general equations of state (EOS). Part of our motivation for doing so is to model transient astrophysics phenomena, as these types of events are often not well approximated by an ideal gas. This necessitated changes to the Riemann solvers implemented in Athena++. We discuss the adjustments made to the HLLC, and HLLD solvers and EOS calls required for arbitrary EOS. We demonstrate the reliability of our code in a number of tests which utilize a relatively simple, but non-trivial EOS based on hydrogen ionization, appropriate for the transition from atomic to ionized hydrogen. Additionally, we perform tests using an electron-positron Helmholtz EOS, appropriate for regimes where nuclear statistical equilibrium is a good approximation. These new complex EOS tests overall show that our modifications to Athena++ accurately solve the Riemann problem with linear convergence and linear-wave tests with quadratic convergence. We provide our test solutions as a means to check the accuracy of other hydrodynamic codes. Our tests and additions to Athena++ will enable further research into (magneto)hydrodynamic problems where realistic treatments of the EOS are required.
△ Less
Submitted 8 April, 2020; v1 submitted 11 September, 2019;
originally announced September 2019.
-
The Accretion History of AGN: A Newly Defined Population of Cold Quasars
Authors:
Allison Kirkpatrick,
C. Megan Urry,
Jason Brewster,
Kevin C. Cooke,
Michael Estrada,
Eilat Glikman,
Kurt Hamblin,
Tonima Tasnim Ananna,
Casey Carlile,
Brandon Coleman,
Jordan Johnson,
Jeyhan S. Kartaltepe,
Stephanie M. LaMassa,
Stefano Marchesi,
Meredith Powell,
Dave Sanders,
Ezequiel Treister,
Tracey Jane Turner
Abstract:
Quasars are the most luminous of active galactic nuclei (AGN), and are perhaps responsible for quenching star formation in their hosts. The Stripe 82X catalog covers 31.3 deg$^2$ of the Stripe 82 field, of which the 15.6 deg$^2$ covered with XMM-Newton is also covered by Herschel/SPIRE. We have 2500 X-ray detected sources with multi-wavelength counterparts, and 30% of these are unobscured quasars,…
▽ More
Quasars are the most luminous of active galactic nuclei (AGN), and are perhaps responsible for quenching star formation in their hosts. The Stripe 82X catalog covers 31.3 deg$^2$ of the Stripe 82 field, of which the 15.6 deg$^2$ covered with XMM-Newton is also covered by Herschel/SPIRE. We have 2500 X-ray detected sources with multi-wavelength counterparts, and 30% of these are unobscured quasars, with $L_X > 10^{44}\,$erg/s and $M_B < -23$. We define a new population of quasars which are unobscured, have X-ray luminosities in excess of $10^{44}\,$erg/s, have broad emission lines, and yet are also bright in the far-infrared, with a 250$μ$m flux density of $S_{\rm 250}>30$mJy. We refer to these Herschel-detected, unobscured quasars as "Cold Quasars". A mere 4% (21) of the X-ray- and optically-selected unobscured quasars in Stripe 82X are detected at 250$μ$m. These Cold Quasars lie at $z\sim1-3$, have $L_{\rm IR}>10^{12}\,L_\odot$, and have star formation rates of $\sim200-1400\,M_\odot$/yr. Cold Quasars are bluer in the mid-IR than the full quasar population, and 72% of our Cold Quasars have WISE W3 $<$ 11.5 [Vega], while only 19% of the full quasar sample meets this criteria. Crucially, Cold Quasars have on average $\sim9\times$ as much star formation as the main sequence of star forming galaxies at similar redshifts. Although dust-rich, unobscured quasars have occasionally been noted in the literature before, we argue that they should be considered as a separate class of quasars due to their high star formation rates. This phase is likely short-lived, as the central engine and immense star formation consume the gas reservoir. Cold Quasars are type-1 blue quasars that reside in starburst galaxies.
△ Less
Submitted 6 July, 2020; v1 submitted 13 August, 2019;
originally announced August 2019.
-
Gravitational-wave-moderated Accretion: The Case of ES Ceti
Authors:
Matthew S. B. Coleman,
Tejaswi Venumadhav,
Barak Zackay
Abstract:
We show that recent observations of the compact binary, AM CVn type system, ES Ceti are fully consistent with theoretical predictions of stable mass transfer moderated by angular momentum loss due to gravitational-wave radiation. One of the main predictions of this model (for degenerate donors) is a widening of the binary. The mass transfer rate inferred from the observed rate of change in the orb…
▽ More
We show that recent observations of the compact binary, AM CVn type system, ES Ceti are fully consistent with theoretical predictions of stable mass transfer moderated by angular momentum loss due to gravitational-wave radiation. One of the main predictions of this model (for degenerate donors) is a widening of the binary. The mass transfer rate inferred from the observed rate of change in the orbital frequency is consistent with that inferred from the observed flux using the recent Gaia DR2 parallax
△ Less
Submitted 12 March, 2019;
originally announced March 2019.
-
Solving the Riemann Problem for Realistic Astrophysical Fluids
Authors:
Zhuo Chen,
Matthew S. B. Coleman,
Eric G. Blackman,
Adam Frank
Abstract:
We present new methods to solve the Riemann problem both exactly and approximately for general equations of state (EoS) to facilitate realistic modeling and understanding of astrophysical flows. The existence and uniqueness of the new exact general EoS Riemann solution can be guaranteed if the EoS is monotone regardless of the physical validity of the EoS. We confirm that: (1) the solution of the…
▽ More
We present new methods to solve the Riemann problem both exactly and approximately for general equations of state (EoS) to facilitate realistic modeling and understanding of astrophysical flows. The existence and uniqueness of the new exact general EoS Riemann solution can be guaranteed if the EoS is monotone regardless of the physical validity of the EoS. We confirm that: (1) the solution of the new exact general EoS Riemann solver and the solution of the original exact Riemann solver match when calculating perfect gas Euler equations; (2) the solution of the new Harten-Lax-van Leer-Contact (HLLC) general EoS Riemann solver and the solution of the original HLLC Riemann solver match when working with perfect gas EoS; and (3) the solution of the new HLLC general EoS Riemann solver approaches the new exact solution. We solve the EoS with two methods, one is to interpolate 2D EoS tables by the bi-linear interpolation method, and the other is to analytically calculate thermodynamic variables at run-time. The interpolation method is more general as it can work with other monotone and realistic EoS while the analytic EoS solver introduced here works with a relatively idealized EoS. Numerical results confirm that the accuracy of the two EoS solvers is similar. We study the efficiency of these two methods with the HLLC general EoS Riemann solver and find that analytic EoS solver is faster in the test problems. However, we point out that a combination of the two EoS solvers may become favorable in some specific problems. Throughout this research, we assume local thermal equilibrium.
△ Less
Submitted 2 November, 2019; v1 submitted 11 March, 2019;
originally announced March 2019.
-
Convection Enhances Magnetic Turbulence in AM CVn Accretion Disks
Authors:
Matthew S. B. Coleman,
Omer Blaes,
Shigenobu Hirose,
Peter H. Hauschildt
Abstract:
We present the results of local, vertically stratified, radiation magnetohydrodynamic shearing box simulations of magnetorotational instability (MRI) turbulence for a (hydrogen poor) composition applicable to accretion disks in AM CVn type systems. Many of these accreting white dwarf systems are helium analogues of dwarf novae (DNe). We utilize frequency-integrated opacity and equation of state ta…
▽ More
We present the results of local, vertically stratified, radiation magnetohydrodynamic shearing box simulations of magnetorotational instability (MRI) turbulence for a (hydrogen poor) composition applicable to accretion disks in AM CVn type systems. Many of these accreting white dwarf systems are helium analogues of dwarf novae (DNe). We utilize frequency-integrated opacity and equation of state tables appropriate for this regime to accurately portray the relevant thermodynamics. We find bistability of thermal equilibria in the effective temperature, surface mass density plane typically associated with disk instabilities. Along this equilibrium curve (i.e. the S-curve) we find that the stress to thermal pressure ratio $α$ varied with peak values of $\sim 0.15$ near the tip of the upper branch. Similar to DNe, we found enhancement of $α$ near the tip of the upper branch caused by convection; this increase in $α$ occurred despite our choice of zero net vertical magnetic flux. Two notable differences we find between DN and AM CVn accretion disk simulations are that AM CVn disks are capable of exhibiting persistent convection in outburst, and ideal MHD is valid throughout quiescence for AM CVns. In contrast, DNe simulations only show intermittent convection, and non-ideal MHD effects are likely important in quiescence. By combining our previous work with these new results, we also find that convective enhancement of the MRI is anticorrelated with mean molecular weight.
△ Less
Submitted 12 March, 2018;
originally announced March 2018.
-
Convective Quenching of Field Reversals in Accretion Disc Dynamos
Authors:
Matthew S. B. Coleman,
Evan Yerger,
Omer Blaes,
Greg Salvesen,
Shigenobu Hirose
Abstract:
Vertically stratified shearing box simulations of magnetorotational turbulence commonly exhibit a so-called butterfly diagram of quasi-periodic azimuthal field reversals. However, in the presence of hydrodynamic convection, field reversals no longer occur. Instead, the azimuthal field strength fluctuates quasi-periodically while maintaining the same polarity, which can either be symmetric or antis…
▽ More
Vertically stratified shearing box simulations of magnetorotational turbulence commonly exhibit a so-called butterfly diagram of quasi-periodic azimuthal field reversals. However, in the presence of hydrodynamic convection, field reversals no longer occur. Instead, the azimuthal field strength fluctuates quasi-periodically while maintaining the same polarity, which can either be symmetric or antisymmetric about the disc midplane. Using data from the simulations of Hirose et al. (2014), we demonstrate that the lack of field reversals in the presence of convection is due to hydrodynamic mixing of magnetic field from the more strongly magnetized upper layers into the midplane, which then annihilate field reversals that are starting there. Our convective simulations differ in several respects from those reported in previous work by others, in which stronger magnetization likely plays a more important role than convection.
△ Less
Submitted 27 January, 2017;
originally announced January 2017.
-
Galactic Bulge Preferred Over Dark Matter for the Galactic Center Gamma-Ray Excess
Authors:
Oscar Macias,
Chris Gordon,
Roland M. Crocker,
Brenna Coleman,
Dylan Paterson,
Shunsaku Horiuchi,
Martin Pohl
Abstract:
An anomalous gamma-ray excess emission has been found in Fermi Large Area Telescope data covering the centre of the Galaxy. Several theories have been proposed for this `Galactic Centre Excess'. They include self-annihilation of dark matter particles, an unresolved population of millisecond pulsars, an unresolved population of young pulsars, or a series of burst events. Here we report on a new a…
▽ More
An anomalous gamma-ray excess emission has been found in Fermi Large Area Telescope data covering the centre of the Galaxy. Several theories have been proposed for this `Galactic Centre Excess'. They include self-annihilation of dark matter particles, an unresolved population of millisecond pulsars, an unresolved population of young pulsars, or a series of burst events. Here we report on a new analysis that exploits hydrodynamical modelling to register the position of interstellar gas associated with diffuse Galactic gamma-ray emission. We find evidence that the Galactic Centre Excess gamma rays are statistically better described by the stellar over-density in the Galactic bulge and the nuclear stellar bulge, rather than a spherical excess. Given its non-spherical nature, we argue that the Galactic Centre Excess is not a dark matter phenomenon but rather associated with the stellar population of the Galactic bulge and nuclear bulge.
△ Less
Submitted 25 September, 2018; v1 submitted 20 November, 2016;
originally announced November 2016.
-
Dwarf Nova Outbursts with Magnetorotational Turbulence
Authors:
M. S. B. Coleman,
I. Kotko,
O. Blaes,
J. -P. Lasota,
S. Hirose
Abstract:
The phenomenological Disc Instability Model has been successful in reproducing the observed light curves of dwarf nova outbursts by invoking an enhanced Shakura-Sunyaev $α$ parameter $\sim0.1-0.2$ in outburst compared to a low value $\sim0.01$ in quiescence. Recent thermodynamically consistent simulations of magnetorotational (MRI) turbulence with appropriate opacities and equation of state for dw…
▽ More
The phenomenological Disc Instability Model has been successful in reproducing the observed light curves of dwarf nova outbursts by invoking an enhanced Shakura-Sunyaev $α$ parameter $\sim0.1-0.2$ in outburst compared to a low value $\sim0.01$ in quiescence. Recent thermodynamically consistent simulations of magnetorotational (MRI) turbulence with appropriate opacities and equation of state for dwarf nova accretion discs have found that thermal convection enhances $α$ in discs in outburst, but only near the hydrogen ionization transition. At higher temperatures, convection no longer exists and $α$ returns to the low value comparable to that in quiescence. In order to check whether this enhancement near the hydrogen ionization transition is sufficient to reproduce observed light curves, we incorporate this MRI-based variation in $α$ into the Disc Instability Model, as well as simulation-based models of turbulent dissipation and convective transport. These MRI-based models can successfully reproduce observed outburst and quiescence durations, as well as outburst amplitudes, albeit with different parameters from the standard Disc Instability Models. The MRI-based model lightcurves exhibit reflares in the decay from outburst, which are not generally observed in dwarf novae. However, we highlight the problematic aspects of the quiescence physics in the Disc Instability Model and MRI simulations that are responsible for this behavior.
△ Less
Submitted 3 August, 2016;
originally announced August 2016.
-
Convection Causes Enhanced Magnetic Turbulence in Accretion Disks in Outburst
Authors:
Shigenobu Hirose,
Omer Blaes,
Julian H. Krolik,
Matthew S. B. Coleman,
Takayoshi Sano
Abstract:
We present the results of local, vertically stratified, radiation MHD shearing box simulations of MRI turbulence appropriate for the hydrogen ionizing regime of dwarf nova and soft X-ray transient outbursts. We incorporate the frequency-integrated opacities and equation of state for this regime, but neglect non-ideal MHD effects and surface irradiation, and do not impose net vertical magnetic flux…
▽ More
We present the results of local, vertically stratified, radiation MHD shearing box simulations of MRI turbulence appropriate for the hydrogen ionizing regime of dwarf nova and soft X-ray transient outbursts. We incorporate the frequency-integrated opacities and equation of state for this regime, but neglect non-ideal MHD effects and surface irradiation, and do not impose net vertical magnetic flux. We find two stable thermal equilibrium tracks in the effective temperature versus surface mass density plane, in qualitative agreement with the S-curve picture of the standard disk instability model. We find that the large opacity at temperatures near $10^4$K, a corollary of the hydrogen ionization transition, triggers strong, intermittent thermal convection on the upper stable branch. This convection strengthens the magnetic turbulent dynamo and greatly enhances the time-averaged value of the stress to thermal pressure ratio $α$, possibly by generating vertical magnetic field that may seed the axisymmetric magnetorotational instability, and by increasing cooling so that the pressure does not rise in proportion to the turbulent dissipation. These enhanced stress to pressure ratios may alleviate the order of magnitude discrepancy between the $α$-values observationally inferred in the outburst state and those that have been measured from previous local numerical simulations of magnetorotational turbulence that lack net vertical magnetic flux.
△ Less
Submitted 11 April, 2014; v1 submitted 12 March, 2014;
originally announced March 2014.