Engineering Papers⌕ Search

Engineering topics

Quataert, Eliot

Publications and source records attributed to Quataert, Eliot.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Local models of two-temperature accretion disc coronae – II. Ion thermal conduction and the absence of disc evaporation

We use local stratified shearing-box simulations with magnetic field-aligned thermal conduction to study an idealized model of the coupling between a cold, radiatively efficient accretion disc, and an overlying, hot, two-temperature corona. Evaporation of a cold disc by conduction from the hot corona has been proposed as a means of mediating the soft-to-hard state transitions observed in X-ray binary systems. We model the coronal plasma in our local disc patch as an MHD fluid subject to both free-streaming ion conduction and a parametrized cooling function that captures the collisional transfer of energy from hot ions to colder, rapidly cooling leptons. In all of our models, independent of the initial net vertical magnetic flux (NF) threading the disc, we find no evidence of disc evaporation. The ion heat flux into the disc is radiated away before conduction can heat the disc’s surface layers. When an initial NF is present, steady-state temperature, density, and outflow velocities in our model coronae are unaffected by conduction. Instead of facilitating disc evaporation, thermal conduction is more likely to feed the disc with plasma condensing out of the corona, particularly in flows without NF. Our work indicates that uncertainties in the amount of NF threading the disc hold far greater influence over whether or not the disc will evaporate into a radiatively inefficient accretion flow compared to thermal conduction. We speculate that a change in net flux mediates disc truncation/evaporation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Pathways to Discovery in Astronomy and Astrophysics for the 2020s

We live in an extraordinary period of discovery in astronomy and astrophysics. Six Nobel Prizes have been awarded over the past decade alone for discoveries based on astronomical data (dark energy, gravitational waves, neutrino oscillations, the discovery of exoplanets, cosmology, and supermassive black holes). Many of the ambitious scientific visions of the 2010 New Worlds, New Horizons (Astro2010)1 decadal survey are being fulfilled, but momentum has only grown. We stand on the threshold of new endeavors that will transform not only our understanding of the universe and the processes and physical paradigms that govern it but also humanity’s place in it. This report of the Committee for a Decadal Survey on Astronomy and Astrophysics 2020 (Astro2020) proposes a broad, integrated plan for space- and ground-based astronomy and astrophysics for the decade 2023–2032.2 It also lays the foundations for further advances in the following decade. This is the seventh in a sequence of decadal survey studies in this field from the National Academies of Sciences, Engineering, and Medicine. This survey examines the program of record, providing advice on the major projects from prior surveys that are yet to be completed. It also lays out priorities for future investments driven by scientific opportunities. The recommendations in this report advance foundational activities that support the people who drive innovation and discovery, and that promote the technologies and tools needed to carry out the science. The report also recommends sustaining activities on a broad range of cost and time scales, as well as activities that enable future visionary projects by maturing them scientifically and technically. Last, the recommendations set in motion the construction of frontier facilities that will change the view and understanding of the cosmos. The survey is bounded by plausible budget scenarios based on briefings from the sponsoring agencies— the Department of Energy (DOE), the National Aeronautics and Space Administration (NASA), and the National Science Foundation (NSF). Within these bounds, the survey aims high, reflecting this time of great scientific promise and progress, with opportunities to pursue some of the most compelling scientific quests of our times.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Millimeter observational signatures of flares in magnetically arrested black hole accretion models

In general relativistic magnetohydrodynamic (GRMHD) simulations, accreted magnetic flux on the black hole horizon episodically decays, during which magnetic reconnection heats up the plasma near the horizon, potentially powering high-energy flares like those observed in M87* and Sgr A*. We study the mm observational counterparts of such flaring episodes in very high resolution GRMHD simulations. The change in 230 GHz flux during the expected high energy flares depends primarily on the efficiency of accelerating γ ≳ 100 (T e ≳ 10 11 K) electrons. For models in which the electrons are heated to T e ∼ 10 11 K during flares, the hot plasma produced by reconnection significantly enhances 230 GHz emission and increases the size of the 230 GHz image. By contrast, for models in which the electrons are heated to higher temperatures (which we argue are better motivated), the reconnection-heated plasma is too hot to produce significant 230 GHz synchrotron emission, and the 230 GHz flux decreases during high energy flares. We do not find a significant change in the mm polarization during flares as long as the emission is Faraday thin. We also present expectations for the ring-shaped image as observed by the Event Horizon Telescope during flares, as well as multiwavelength synchrotron spectra. Our results highlight several limitations of standard post-processing prescriptions for the electron temperature in GRMHD simulations. We also discuss the implications of our results for current and future observations of flares in Sgr A*, M87*, and related systems. Appendices contain detailed convergence studies with respect to resolution and plasma magnetization.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Large-scale Evolution of Seconds-long Relativistic Jets from Black Hole–Neutron Star Mergers

We present the first numerical simulations that track the evolution of a black hole–neutron star (BH–NS) merger from premerger to r ≳ 10 11 cm. The disk that forms after a merger of mass ratio q = 2 ejects massive disk winds (3–5 × 10 –2 M ⊙ ). We introduce various postmerger magnetic configurations and find that initial poloidal fields lead to jet launching shortly after the merger. The jet maintains a constant power due to the constancy of the large-scale BH magnetic flux until the disk becomes magnetically arrested (MAD), where the jet power falls off as L j ~ t –2 . All jets inevitably exhibit either excessive luminosity due to rapid MAD activation when the accretion rate is high or excessive duration due to delayed MAD activation compared to typical short gamma-ray bursts (sGRBs). This provides a natural explanation for long sGRBs such as GRB 211211A but also raises a fundamental challenge to our understanding of jet formation in binary mergers. One possible implication is the necessity of higher binary mass ratios or moderate BH spins to launch typical sGRB jets. For postmerger disks with a toroidal magnetic field, dynamo processes delay jet launching such that the jets break out of the disk winds after several seconds. We show for the first time that sGRB jets with initial magnetization σ0 > 100 retain significant magnetization (σ $\gg$ 1) at r > 10 10 cm, emphasizing the importance of magnetic processes in the prompt emission. The jet–wind interaction leads to a power-law angular energy distribution by inflating an energetic cocoon whose emission is studied in a companion paper.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Hours-long Near-UV/Optical Emission from Mildly Relativistic Outflows in Black Hole–Neutron Star Mergers

The ongoing LIGO–Virgo–KAGRA observing run O4 provides an opportunity to discover new multimessenger events, including binary neutron star (BNS) mergers such as GW170817 and the highly anticipated first detection of a multimessenger black hole–neutron star (BH–NS) merger. While BNS mergers were predicted to exhibit early optical emission from mildly relativistic outflows, it has remained uncertain whether the BH–NS merger ejecta provides the conditions for similar signals to emerge. We present the first modeling of early near-ultraviolet/optical emission from mildly relativistic outflows in BH–NS mergers. Adopting optimal binary properties, a mass ratio of q = 2, and a rapidly rotating BH, we utilize numerical relativity and general relativistic magnetohydrodynamic (GRMHD) simulations to follow the binary's evolution from premerger to homologous expansion. We use an M1 neutrino transport GRMHD simulation to self-consistently estimate the opacity distribution in the outflows and find a bright near-ultraviolet/optical signal that emerges due to jet-powered cocoon cooling emission, outshining the kilonova emission at early time. The signal peaks at an absolute magnitude of ~–15 a few hours after the merger, longer than previous estimates, which did not consider the first principles–based jet launching. By late 2024, the Rubin Observatory will have the capability to track the entire signal evolution or detect its peak up to distances of ≳1 Gpc. In 2026, ULTRASAT will conduct all-sky surveys within minutes, detecting some of these events within ~200 Mpc. The BH–NS mergers with higher mass ratios or lower BH spins would produce shorter and fainter signals.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Multi-scale Dynamics of Kinetic Turbulence in Weakly Collisional, High-Beta Plasmas (Final Report for DOE Grant DE-SC0019046)

This grant was a collaborative grant between Prof Matt Kunz (Princeton) and Prof. Eliot Quataert. Quataert was initially a faculty member at UC Berkeley when the grant was funded but moved to Princeton during the timeframe of this grant. Over the course of this grant we made major progress on understanding turbulence and multi-scale dynamics in weakly collisional high-beta plasmas. The most important contributions included: We developed one of the most compelling theoretical explanations for a decades-old puzzle at the heart of our understanding of the origin of the solar wind. The puzzle is that many of the observations of ion temperatures and distribution functions in the solar wind are consistent with cyclotron resonant heating. However, theoretical models of MHD turbulence in the solar wind show that most of the turbulent energy remains at low frequencies below the cyclotron frequency. In Squire et al. (2022), we showed, however, that in imbalanced turbulence (in which there is an asymmetry in the Alfven-wave flux in opposite directions along the magnetic field), which is the norm in the fast solar wind, the turbulent energy reaches a ‘bottleneck’ near the ion Larmor radius, and the amplitude and characteristic frequency of the fluctuations grow until the cyclotron frequency is reached.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Angular-momentum Transport in Proto-neutron Stars and the Fate of Neutron Star Merger Remnants

Both the core collapse of rotating massive stars, and the coalescence of neutron star (NS) binaries result in the formation of a hot, differentially rotating NS remnant. The timescales over which differential rotation is removed by internal angular-momentum transport processes (viscosity) have key implications for the remnant’s long-term stability and the NS equation of state (EOS). Guided by a nonrotating model of a cooling proto-NS, we estimate the dominant sources of viscosity using an externally imposed angular-velocity profile Ω(r). Although the magneto-rotational instability provides the dominant source of effective viscosity at large radii, convection and/or the Tayler–Spruit dynamo dominate in the core of merger remnants where dΩ/dr ≥ 0. Furthermore, the viscous timescale in the remnant core is sufficiently short that solid-body rotation will be enforced faster than matter is accreted from rotationally supported outer layers. Guided by these results, we develop a toy model for how the merger remnant core grows in mass and angular momentum due to accretion. We find that merger remnants with sufficiently massive and slowly rotating initial cores may collapse to black holes via envelope accretion, even when the total remnant mass is less than the usually considered threshold ≈1.2 M TOV for forming a stable solid-body rotating NS remnant (where M TOV is the maximum nonrotating NS mass supported by the EOS). This qualitatively new picture of the post-merger remnant evolution and stability criterion has important implications for the expected electromagnetic counterparts from binary NS mergers and for multimessenger constraints on the NS EOS

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Birth of a Be star: an APOGEE search for Be stars forming through binary mass transfer

ABSTRACT Motivated by recent suggestions that many Be stars form through binary mass transfer, we searched the APOGEE survey for Be stars with bloated, stripped companions. From a well-defined parent sample of 297 Be stars, we identified one mass-transfer binary, HD 15124. The object consists of a main-sequence Be star ($M_{\rm Be}=5.3\pm 0.6\, {\rm M}_{\odot }$) with a low-mass ($M_{\rm donor}=0.92\pm 0.22\, {\rm M}_{\odot }$), subgiant companion on a 5.47-d orbit. The emission lines originate in an accretion disc caused by ongoing mass transfer, not from a decretion disc as in classical Be stars. Both stars have surface abundances bearing imprint of CNO processing in the donor’s core: the surface helium fraction is YHe ≈ 0.6, and the nitrogen-to-carbon ratio is 1000 times the solar value. The system’s properties are well-matched by binary evolution models in which mass transfer begins while a $3-5\, {\rm M}_{\odot }$ donor leaves the main sequence, with the originally less massive component becoming the Be star. These models predict that the system will soon become a detached Be + stripped star binary like HR 6819 and LB-1, with the stripped donor eventually contracting to become a core helium-burning sdO/B star. Discovery of one object in this short-lived (∼1 Myr) evolutionary phase implies the existence of many more that have already passed through it and are now Be + sdO/B binaries. We infer that $(10-60)\, {{\ \rm per\ cent}}$ of Be stars have stripped companions, most of which are $\sim 100\, \times$ fainter than the Be stars in the optical. Together with the dearth of main-sequence companions to Be stars and recent discovery of numerous Be + sdO/B binaries in the UV, our results imply that binarity plays an important role in the formation of Be stars.

El-Badry, Kareem (ORCID:0000000268711752)↗
DOE OSTI2022

High-frequency heating of the solar wind triggered by low-frequency turbulence

The fast solar wind’s high speeds and non-thermal features require that considerable heating occurs well above the Sun’s surface. Two leading theories seem incompatible: low-frequency ‘Alfvénic’ turbulence, which transports energy outwards and is observed ubiquitously by spacecraft but seems insufficient to explain the observed dominance of ion over electron heating; and high-frequency ion-cyclotron waves, which explain the non-thermal heating of ions but lack an obvious source. Here we argue that the recently proposed ‘helicity barrier’ effect, which limits electron heating by inhibiting the turbulent cascade of energy to the smallest scales, can unify these two paradigms. Our six-dimensional simulations show how the helicity barrier causes the large-scale energy to grow through time, generating small parallel scales and high-frequency ion-cyclotron-wave heating from low-frequency turbulence, while simultaneously explaining various other long-standing observational puzzles. Furthermore, the predicted causal link between plasma expansion and the ion-to-electron heating ratio suggests that the helicity barrier could contribute to key observed differences between fast and slow wind streams.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The impact of r -process heating on the dynamics of neutron star merger accretion disc winds and their electromagnetic radiation

Neutron star merger accretion discs can launch neutron-rich winds of >10 -2 M⊙. This ejecta is a prime site for r-process nucleosynthesis, which will produce a range of radioactive heavy nuclei. The decay of these nuclei releases enough energy to accelerate portions of the wind by ~0.1c. Here, we investigate the effect of r-process heating on the dynamical evolution of disc winds. We extract the wind from a 3D general relativistic magnetohydrodynamic simulation of a disc from a post-merger system. This is used to create inner boundary conditions for 2D hydrodynamic simulations that continue the original 3D simulation. We perform two such simulations: one that includes the r-process heating, and another one that does not. We follow the hydrodynamic simulations until the winds reach homology (60 s). Using time-dependent multifrequency multidimensional Monte Carlo radiation transport simulations, we then calculate the kilonova light curves from the winds with and without dynamical r-process heating. We find that the r-process heating can substantially alter the velocity distribution of the wind, shifting the mass-weighted median velocity from 0.06c to 0.12c. The inclusion of the dynamical r-process heating makes the light curve brighter and bluer at $\sim 1\, \mathrm{d}$ post-merger. However, the high-velocity tail of the ejecta distribution and the early ($\lesssim 1\, \mathrm{d}$) light curves are largely unaffected.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The effect of jet–ejecta interaction on the viewing angle dependence of kilonova light curves

ABSTRACT The merger of two neutron stars produces an outflow of radioactive heavy nuclei. Within a second of merger, the central remnant is expected to also launch a relativistic jet, which shock-heats and disrupts a portion of the radioactive ejecta. Within a few hours, emission from the radioactive material gives rise to an ultraviolet, optical, and infrared transient (a kilonova). We use the endstates of a suite of 2D relativistic hydrodynamic simulations of jet–ejecta interaction as initial conditions for multidimensional Monte Carlo radiation transport simulations of the resulting viewing angle-dependent light curves and spectra starting at $1.5\, \mathrm{h}$ after merger. We find that on this time-scale, jet shock heating does not affect the kilonova emission for the jet parameters we survey. However, the jet disruption to the density structure of the ejecta does change the light curves. The jet carves a channel into the otherwise spheroidal ejecta, revealing the hot, inner regions. As seen from near (≲30°) the jet axis, the kilonova is brighter by a factor of a few and bluer. The strength of this effect depends on the jet parameters, since the light curves of more heavily disrupted ejecta are more strongly affected. The light curves and spectra are also more heavily modified in the ultraviolet than in the optical.

Klion, Hannah↗
NASA NTRS2019

Lynx X-ray Observatory

X-ray observations are indispensable for understanding the cosmos. Their power is immense because much of the baryonic matter and the sites for the most active energy releases in the Universe are primarily observable in X-rays. For the 2030s and beyond, an X-ray observatory with power matching the capabilities in other wavebands is a necessary discovery engine for full exploration of the Universe. JWST and other upcoming major space- and ground-based facilities are expected to greatly expand science frontiers in the coming decades. is presents both a great opportunity and a challenge for a next-generation X-ray observatory. In many areas, such as tracing black holes during the CosmicDawn and understanding the formation and evolution of galaxies, an X-ray observatory is the logical next step. e challenge is that the X-ray science at these new frontiers requires expansion of capabilities by orders of magnitude beyond the current state of the art or anything already planned. Until recently, such gains were not technologically possible. is has changed thanks to recent breakthroughs and sustained maturation of key technologies for X-ray mirrors and detectors. We are reaping the fruits of U.S. investments in these areas over the past 10–15 years. An X-ray observatory that can extend the science frontiers of the post-JWST era is now entirely feasible. Lynx is the mission concept that realizes this vision. It will y revolutionary optics and instrumentation onboard a simple, proven spacecraft. In all aspects, Lynx will be a next-generation Great Observatory that is certain to make a profound impact across the astrophysical landscape. It will provide the depth and breadth to answer the fundamental questions that confront us today; just as importantly, it will have capabilities to address questions we have yet to even ask.

Gaskin, Jessica A.↗
DOE OSTI2018

Long-term GRMHD simulations of neutron star merger accretion discs: implications for electromagnetic counterparts

We investigate the long-term evolution of black hole accretion disks formed in neutron star mergers. These disks expel matter that contributes to an r-process kilonova, and can produce relativistic jets powering short gamma-ray bursts. Here we report the results of a three-dimensional, general-relativistic magnetohydrodynamic (GRMHD) simulation of such a disk which is evolved for long enough (~9s, or ~6×10 5 r$_g$/c) to achieve completion of mass ejection far from the disk. Our model starts with a poloidal field, and fully resolves the most unstable mode of the magnetorotational instability. We parameterize the dominant microphysics and neutrino cooling effects, and compare with axisymmetric hydrodynamic models with shear viscosity. The GRMHD model ejects mass in two ways: a prompt MHD-mediated outflow and a late-time, thermally-driven wind once the disk becomes advective. The total amount of unbound mass ejected (0.013M ⊙ , or ≃40% of the initial torus mass) is twice as much as in hydrodynamic models, with higher average velocity (0.1c) and a broad electron fraction distribution with a lower average value (0.16). Scaling the ejected fractions to a disk mass of ~0.1M ⊙ can account for the red kilonova from GW170817 but underpredicts the blue component. About ~10-3M ⊙ of material should undergo neutron freezout and could produce a bright kilonova precursor in the first few hours after the merger. With our idealized initial magnetic field configuration, we obtain a robust jet and sufficient ejecta with Lorentz factor ~1-10 to (over)produce the non-thermal emission from GW1708107.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Jet Dynamics in Compact Object Mergers: GW170817 Likely Had a Successful Jet

Here, we use relativistic hydrodynamic numerical calculations to study the interaction between a jet and a homologous outflow produced dynamically during binary neutron star mergers. We quantify how the thermal energy imparted by the jet and the ability of the jet to escape the ejecta depend on the parameters of the jet engine and the ejecta. Under our assumptions, a collimated jet initiated at early times compared to the engine duration, we show that successful breakout of the forward cocoon shock necessitates a jet that successfully escapes the ejecta. This is because the ejecta is expanding, and the forward shock from a failed jet stalls before it reaches the edge of the ejecta. This conclusion can be circumvented only for very energetic wide angle jets, with parameters that are uncomfortable given short-duration GRB observations. For successful jets, we find two regimes of jet breakout from the ejecta: early breakout on timescales shorter than the engine duration, and late breakout well after the engine shuts off. A late breakout can explain the observed delay between gravitational waves and gamma rays in GW170817. We show that for the entire parameter space of jet parameters surveyed here (covering energies ~10 48 –10 51 erg and opening angles θ j ~ 0.07–0.4) the thermal energy deposited by the jet is less than that produced by r-process heating on second timescales by at least an order of magnitude. Shock heating is thus energetically subdominant in setting the luminosity of thermally powered transients coincident with neutron star mergers (kilonovae).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Two-temperature GRRMHD Simulations of M87

Here, we present axisymmetric two-temperature general relativistic radiation magnetohydrodynamic simulations of the inner region of the accretion flow onto the supermassive black hole M87. We address uncertainties from previous modeling efforts through inclusion of models for (1) self-consistent dissipative and Coulomb electron heating (2) radiation transport (3) frequency-dependent synchrotron emission, self-absorption, and Compton scattering. We adopt a distance D = 16.7 Mpc, an observer angle θ = 20°, and consider black hole masses M/M ⊙ = (3.3 x 10 9 , 6.2 x 10 9 ) and spins a star = (0.5, 0.9375) in a four-simulation suite. For each (M, a star ), we identify the accretion rate that recovers the 230 GHz flux from very long baseline interferometry measurements. We report on disk thermodynamics at these accretion rates (M · /M · Edd ~ 10 -5 ). The disk remains geometrically thick; cooling does not lead to a thin disk component. While electron heating is dominated by Coulomb rather than dissipation for r ≳ 10GM/c 2 , the accretion disk remains two-temperature. Radiative cooling of electrons is not negligible, especially for r ≲ 10GM/c 2 . The Compton y parameter is of order unity. We then compare derived and observed or inferred spectra, millimeter images, and jet powers. Simulations with M/M ⊙ = 3.3 × 10 9 are in conflict with observations. These simulations produce millimeter images that are too small, while the low-spin simulation also overproduces X-rays. For M/M ⊙ = 6.2 x 10 9 , both simulations agree with constraints on radio/IR/X-ray fluxes and millimeter image sizes. Simulation jet power is a factor 10 2 –10 3 below inferred values, a possible consequence of the modest net magnetic flux in our models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Outbursts of luminous blue variable stars from variations in the helium opacity

Luminous blue variables are massive, evolved stars that exhibit large variations in luminosity and size on timescales from months to years, with high associated rates of mass loss1,2,3,4,5. In addition to this on-going variability, these stars exhibit outburst phases, during which their size increases and as a result their effective temperature decreases, typically to about 9,000 kelvin3,6. Outbursts are believed to be caused by the radiation force on the cooler, more opaque, outer layers of the star balancing or even exceeding the force of gravity, although the exact mechanisms are unknown and cannot be determined using one-dimensional, spherically symmetric models of stars because such models cannot determine the physical processes that occur in this regime7. Here we report three-dimensional simulations of massive, radiation-dominated stars, which show that helium opacity has an important role in triggering outbursts and setting the observed effective temperature during outbursts of about 9,000 kelvin. It probably also triggers the episodic mass loss at rates of 10 -7 to 10 -5 solar masses per year. The peak in helium opacity is evident in our three-dimensional simulations only because the density and temperature of the stellar envelope (the outer part of the star near the photosphere) need to be determined self-consistently with convection, which cannot be done in one-dimensional models that assume spherical symmetry. The simulations reproduce observations of long-timescale variability, and predict that convection causes irregular oscillations in the radii of the stars and variations in brightness of 10–30 per cent on a typical timescale of a few days. The amplitudes of these short-timescale variations are predicted to be even larger for cooler stars (in the outburst phase). This short-timescale variability should be observable with high-cadence observations.

Science & Technology - Other Topics↗