Researcher profile

Eliot Quataert

Eliot Quataert contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
42works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

42 published item(s)

preprint2022arXiv

Cosmic-Ray Driven Outflows to Mpc Scales from $L_{\ast}$ Galaxies

We study the effects of cosmic rays (CRs) on outflows from star-forming galaxies in the circum and inter-galactic medium (CGM/IGM), in high-resolution, fully-cosmological FIRE-2 simulations (accounting for mechanical and radiative stellar feedback, magnetic fields, anisotropic conduction/viscosity/CR diffusion and streaming, and CR losses). We showed previously that massive ($M_{\rm halo}\gtrsim 10^{11}\,M_{\odot}$), low-redshift ($z\lesssim 1-2$) halos can have CR pressure dominate over thermal CGM pressure and balance gravity, giving rise to a cooler CGM with an equilibrium density profile. This dramatically alters outflows. Absent CRs, high gas thermal pressure in massive halos "traps" galactic outflows near the disk, so they recycle. With CRs injected in supernovae as modeled here, the low-pressure halo allows "escape" and CR pressure gradients continuously accelerate this material well into the IGM in "fast" outflows, while lower-density gas at large radii is accelerated in-situ into "slow" outflows that extend to $>$Mpc scales. CGM/IGM outflow morphologies are radically altered: they become mostly volume-filling (with inflow in a thin mid-plane layer) and coherently biconical from the disk to $>$Mpc. The CR-driven outflows are primarily cool ($T\sim10^{5}\,$K) and low-velocity. All of these effects weaken and eventually vanish at lower halo masses ($\lesssim 10^{11}\,M_{\odot}$) or higher redshifts ($z\gtrsim 1-2$), reflecting the ratio of CR to thermal+gravitational pressure in the outer halo. We present a simple analytic model which explains all of the above phenomena.

preprint2022arXiv

Effects of Different Cosmic Ray Transport Models on Galaxy Formation

Cosmic rays (CRs) with ~GeV energies can contribute significantly to the energy and pressure budget in the interstellar, circumgalactic, and intergalactic medium (ISM, CGM, IGM). Recent cosmological simulations have begun to explore these effects, but almost all studies have been restricted to simplified models with constant CR diffusivity and/or streaming speeds. Physical models of CR propagation/scattering via extrinsic turbulence and self-excited waves predict transport coefficients which are complicated functions of local plasma properties. In a companion paper, we consider a wide range of observational constraints to identify proposed physically-motivated cosmic-ray propagation scalings which satisfy both detailed Milky Way (MW) and extra-galactic $γ$-ray constraints. Here, we compare the effects of these models relative to simpler 'diffusion+streaming' models on galaxy and CGM properties at dwarf through MW mass scales. The physical models predict large local variations in CR diffusivity, with median diffusivity increasing with galacto-centric radii and decreasing with galaxy mass and redshift. These effects lead to a more rapid dropoff of CR energy density in the CGM (compared to simpler models), in turn producing weaker effects of CRs on galaxy star formation rates (SFRs), CGM absorption profiles and galactic outflows. The predictions of the more physical CR models tend to lie 'in between' models which ignore CRs entirely and models which treat CRs with constant diffusivity.

preprint2022arXiv

Galaxies lacking dark matter produced by close encounters in a cosmological simulation

The standard cold dark matter plus cosmological constant model predicts that galaxies form within dark-matter haloes, and that low-mass galaxies are more dark-matter dominated than massive ones. The unexpected discovery of two low-mass galaxies lacking dark matter immediately provoked concerns about the standard cosmology and ignited explorations of alternatives, including self-interacting dark matter and modified gravity. Apprehension grew after several cosmological simulations using the conventional model failed to form adequate numerical analogues with comparable internal characteristics (stellar masses, sizes, velocity dispersions and morphologies). Here we show that the standard paradigm naturally produces galaxies lacking dark matter with internal characteristics in agreement with observations. Using a state-of-the-art cosmological simulation and a meticulous galaxy-identification technique, we find that extreme close encounters with massive neighbours can be responsible for this. We predict that approximately 30 percent of massive central galaxies (with at least 1e11 solar masses in stars) harbour at least one dark-matter-deficient satellite (with 1e8 - 1e9 solar masses in stars). This distinctive class of galaxies provides an additional layer in our understanding of the role of interactions in shaping galactic properties. Future observations surveying galaxies in the aforementioned regime will provide a crucial test of this scenario.

preprint2022arXiv

Gas infall and radial transport in cosmological simulations of Milky Way-mass disks

Observations indicate that a continuous supply of gas is needed to maintain observed star formation rates in large, disky galaxies. To fuel star formation, gas must reach the inner regions of such galaxies. Despite its crucial importance for galaxy evolution, how and where gas joins galaxies is poorly constrained observationally and is rarely explored in fully cosmological simulations. To investigate gas accretion in the vicinity of galaxies, we analyze the FIRE-2 cosmological zoom-in simulations for 4 Milky Way mass galaxies (M_halo ~ 10E12 solar masses), focusing on simulations with cosmic ray physics. We find that at z~0, gas approaches the disk with angular momentum similar to the gaseous disk edge and low radial velocities, piling-up near the edge and settling into full rotational support. Accreting gas moves predominantly parallel to the disk with small but nonzero vertical velocity components, and joins the disk largely in the outskirts as opposed to "raining" down onto the disk. Once in the disk, gas trajectories are complex, being dominated by spiral arm induced oscillations and feedback. However, time and azimuthal averages show clear but slow net radial infall with transport speeds of 1-3 km/s and net mass fluxes through the disk on the order of one solar mass per year, comparable to the star formation rates of the galaxies and decreasing towards galactic center as gas is sunk into star formation. These rates are slightly higher in simulations without cosmic rays (1-7 km/s, ~4-5 solar masses per year). We find overall consistency of our results with observational constraints and discuss prospects of future observations of gas flows in and around galaxies.

preprint2022arXiv

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

The fast solar wind's high speeds and nonthermal features require that significant 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 struggles to explain the observed dominance of ion over electron heating; and high-frequency ion-cyclotron waves (ICWs), which explain the nonthermal 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 in time, generating small parallel scales and high-frequency ICW heating from low-frequency turbulence. The resulting turbulence and ion distribution function also closely match in-situ measurements from Parker Solar Probe and other spacecraft, explaining, among other features, the decades-long puzzle of the steep "transition range" observed in magnetic fluctuation spectra. The theory predicts a causal link between plasma expansion and the ion-to-electron heating ratio. Given the observational association between wind speed and expansion, we argue that the helicity barrier could play a key role in regulating the bimodal speed distribution of the solar wind.

preprint2022arXiv

Hot-mode accretion and the physics of thin-disk galaxy formation

We use FIRE simulations to study disk formation in z~0, Milky Way-mass galaxies, and conclude that a key ingredient for the formation of thin stellar disks is the ability for accreting gas to develop an aligned angular momentum distribution via internal cancellation *prior* to joining the galaxy. Among galaxies with a high fraction (>70%) of their young stars in a thin disk (h/R~0.1) we find that: (i) hot, virial-temperature gas dominates the inflowing gas mass on halo scales (>~20 kpc), with radiative losses offset by compression heating; (ii) this hot accretion proceeds until angular momentum support slows inward motion, at which point the gas cools to T~10^4 K or less; (iii) prior to cooling, the accreting gas develops an angular momentum distribution that is aligned with the galaxy disk, and while cooling transitions from a quasi-spherical spatial configuration to a more flattened, disk-like configuration. We show that the existence of this "rotating cooling flow" accretion mode is strongly correlated with the fraction of stars forming in a thin disk among a sample of 17 z~0 galaxies spanning a halo mass range of 10^10.5 solar masses to 10^12 solar masses, or a stellar mass range 10^8 solar masses to 10^11 solar masses. Notably, galaxies with a thick disk or irregular morphology do not undergo significant angular momentum alignment of gas prior to accretion and show no correspondence between halo gas cooling and flattening. Our results suggest that rotating cooling flows (or, more generally, rotating subsonic flows) that become coherent and angular momentum-supported prior to accretion onto the galaxy are likely a necessary condition for the formation of thin, star-forming disk galaxies in a LambdaCDM universe.

preprint2022arXiv

Long-term evolution of supercritical black hole accretion with outflows: a subgrid feedback model for cosmological simulations

We study the long-term evolution of the global structure of axisymmetric accretion flows onto a black hole (BH) at rates substantially higher than the Eddington value ($\dot{M}_{\rm Edd}$), performing two-dimensional hydrodynamical simulations with and without radiative diffusion. In the high-accretion optically-thick limit, where the radiation energy is efficiently trapped within the inflow, the accretion flow becomes adiabatic and comprises of turbulent gas in the equatorial region and strong bipolar outflows. As a result, the mass inflow rate decreases toward the center as $\dot{M}_{\rm in}\propto r^{p}$ with $p\sim 0.5-0.7$ and a small fraction of the inflowing gas feeds the nuclear BH. Thus, super-Eddington accretion is sustained only when a larger amount of gas is supplied from larger radii at $> 100-1000~\dot{M}_{\rm Edd}$. The global structure of the flow settles down to a quasi-steady state in millions of the orbital timescale at the BH event horizon, which is $> 10-100$ times longer than that addressed in previous (magneto-)RHD simulation studies. Energy transport via radiative diffusion accelerates the outflow near the poles in the inner region but does not change the overall properties of the accretion flow compared to the cases without diffusion. Based on our simulation results, we provide a mechanical feedback model for super-Eddington accreting BHs. This can be applied as a sub-grid model in large-scale cosmological simulations that do not sufficiently resolve galactic nuclei, and to the formation of the heaviest gravitational-wave sources via accretion in dense environments.

preprint2022arXiv

Numerical Simulations of the Random Angular Momentum in Convection: Implications for Supergiant Collapse to Form Black Holes

During the core collapse of massive stars that do not undergo a canonical energetic explosion, some of the hydrogen envelope of a red supergiant (RSG) progenitor may infall onto the newborn black hole (BH). Within the Athena++ framework, we perform three-dimensional, hydrodynamical simulations of idealized models of supergiant convection and collapse in order to assess whether the infall of the convective envelope can give rise to rotationally-supported material, even if the star has zero angular momentum overall. Our dimensionless, polytropic models are applicable to the optically-thick hydrogen envelope of non-rotating RSGs and cover a factor of 20 in stellar radius. At all radii, the specific angular momentum due to random convective flows implies associated circularization radii of 10 - 1500 times the innermost stable circular orbit of the BH. During collapse, the angular momentum vector of the convective flows is approximately conserved and is slowly varying on the timescale relevant to forming disks at small radii. Our results indicate that otherwise failed explosions of RSGs lead to the formation of rotationally-supported flows that are capable of driving outflows to large radii and powering observable transients. When the BH is able to accrete most of the hydrogen envelope, the final BH spin parameter is $\sim$ 0.5, even though the star is non-rotating. For fractional accretion of the envelope, the spin parameter is generally lower and never exceeds 0.8. We discuss the implications of our results for transients produced by RSG collapse to a black hole.

preprint2022arXiv

Observational Signatures of Black Hole Accretion: Rotating vs. Spherical Flows with Tilted Magnetic Fields

We study the observational signatures of magnetically arrested black hole accretion with non-rotating inflow onto a rotating black hole; we consider a range of angles between the black hole spin and the initial magnetic field orientation. We compare the results of our General Relativistic Magneto-Hydrodynamic simulations to more commonly used rotating initial conditions and to the Event Horizon Telescope (EHT) observations of M87. We find that the mm intensity images, polarization images, and synchrotron emission spectra are very similar among the different simulations when post-processed with the same electron temperature model; observational differences due to different electron temperature models are significantly larger than those due to the different realizations of magnetically arrested accretion. The orientation of the mm synchrotron polarization is particularly insensitive to the initial magnetic field orientation, the electron temperature model, and the rotation of the inflowing plasma. The largest difference among the simulations with different initial rotation and magnetic tilt is in the strength and stability of the jet; spherical inflow leads to kink-unstable jets. We discuss the implications of our results for current and future EHT observations and for theoretical models of event-horizon-scale black hole accretion.

preprint2022arXiv

Optical to X-ray Signatures of Dense Circumstellar Interaction in Core-Collapse Supernovae

Progenitors of core-collapse supernovae (SNe) can shed significant mass to circumstellar material (CSM) in the months--years preceding core-collapse. The ensuing SN explosion launches ejecta that may subsequently collide with this CSM, producing shocks that can power emission across the electromagnetic spectrum. In this work we explore the thermal signatures of dense CSM interaction, when the CSM density profile is truncated at some outer radius. CSM with optical depth $>c/v$ (where $v$ is the shock velocity) will produce primarily $\sim$blackbody optical/UV emission whereas lower optical-depth CSM will power bremsstrahlung X-ray emission. Focusing on the latter, we derive light-curves and spectra of the resulting X-ray transients, that include a detailed treatment of Comptonization. Due to strong photoelectric absorption, the X-ray light-curve is dominated by the `post-interaction' phase that occurs after the shock reaches the CSM truncation radius. We treat this regime here for the first time. Using these results, we present the phase-space of optical, UV, and X-ray transients as a function of CSM properties, and discuss detectability prospects. We find that ROSAT would not have been sensitive to CSM X-ray transients but that eROSITA is expected to detect many such events. Future wide-field UV missions such as ULTRASAT will dramatically enhance sensitivity to large optical-depth CSM configurations. Finally, we present a framework within which CSM properties may be directly inferred from observable features of X-ray transients. This can serve as an important tool for studying stellar mass loss using SN X-ray detections.

preprint2022arXiv

Reconciling cosmic-ray transport theory with phenomenological models motivated by Milky-Way data

Phenomenological models of cosmic-ray (CR) transport in the Milky Way (MW) can reproduce a wide range of observations assuming that CRs scatter off of magnetic-field fluctuations with spectrum $\propto k^{-δ}$ and $δ\sim [1.4,1.67]$. We study the extent to which such models can be reconciled with current microphysical theories of CR transport, specifically self-confinement due to the streaming instability and/or extrinsic turbulence due to a cascade of MHD fast modes. We first review why it is that on their own neither theory is compatible with observations. We then highlight that CR transport is a strong function of local plasma conditions in the multi-phase interstellar medium (ISM), and may be diffusive due to turbulence in some regions and streaming due to self-confinement in others. A multi-phase combination of scattering mechanisms can in principle reproduce the main trends in the proton spectrum and the boron-to-carbon ratio (B/C). However, models with a combination of scattering by self-excited waves and fast-mode turbulence require significant fine-tuning due to fast-mode damping, unlike phenomenological models that assume undamped Kolmogorov turbulence. The assumption that fast modes follow a weak cascade is also not well justified theoretically, as the weak cascade is suppressed by wave steepening and weak-shock dissipation even in subsonic turbulence. These issues suggest that there may be a significant theoretical gap in our understanding of MHD turbulence. We discuss a few topics at the frontier of MHD turbulence theory that bear on this (possible) gap and that may be relevant for CR scattering.

preprint2022arXiv

Research Opportunities in Plasma Astrophysics

Major scientific questions and research opportunities are described on 10 unprioritized plasma astrophysics topics: (1) magnetic reconnection, (2) collisionless shocks and particle acceleration, (3) waves and turbulence, (4) magnetic dynamos, (5) interface and shear instabilities, (6) angular momentum transport, (7) dusty plasmas, (8) radiative hydrodynamics, (9) relativistic, pair-dominated and strongly magnetized plasmas, (10) jets and outflows. Note that this is a conference report from a Workshop on Opportunities in Plasma Astrophysics (WOPA, https://w3.pppl.gov/conferences/2010/WOPA/) in January 2010, that attracted broad representation from the community and was supported by the U.S. Department of Energy, National Aeronautics and Space Administration, National Science Foundation, American Physical Society's Topical Group for Plasma Astrophysics and Division of Plasma Physics, and Center for Magnetic Self-Organization in Laboratory and Astrophysical Plasmas. Although there has been much planning and many developments in both science and infrastructure since the report was written, most of the motivation, priorities, problems and technical challenges discussed therein remain unaddressed and are relevant at the time of posting.

preprint2022arXiv

Supercritical growth pathway to overmassive black holes at cosmic dawn: coevolution with massive quasar hosts

Observations of the most luminous quasars at high redshifts ($z > 6$) have revealed that the largest supermassive black holes (SMBHs) at those epochs tend to be substantially overmassive relative to their host galaxies compared to the local relations, suggesting they experienced rapid early growth phases. We propose an assembly model for the SMBHs that end up in rare massive $\sim10^{12}~M_{\odot}$ host halos at $z \sim 6-7$, applying a kinetic feedback prescription for BHs accreting above the Eddington rate, provided by radiation hydrodynamic simulations for the long-term evolution of the accretion-flow structure. The large inflow rates into these halos during their assembly enable the formation of $>10^9~M_{\odot}$ SMBHs by $z \sim 6$, even starting from stellar-mass seeds at $z \sim 30$, and even in the presence of outflows that reduce the BH feeding rate, especially at early times. This mechanism also naturally yields a high BH-to-galaxy mass ratio of $> 0.01$ before the SMBH mass reaches $M_{\rm BH} > 10^9~M_{\odot}$ by $z \sim 6$. These fast-growing SMBH progenitors are bright enough to be detected by upcoming observations with the James Webb Space Telescope over a wide range of redshift ($7 < z < 15$), regardless of how they were seeded.

preprint2022arXiv

Testing Physical Models for Cosmic Ray Transport Coefficients on Galactic Scales: Self-Confinement and Extrinsic Turbulence at GeV Energies

The microphysics of ~GeV cosmic ray (CR) transport on galactic scales remain deeply uncertain, with almost all studies adopting simple prescriptions (e.g. constant-diffusivity). We explore different physically-motivated, anisotropic, dynamical CR transport scalings in high-resolution cosmological FIRE simulations of dwarf and ~$L_{\ast}$ galaxies where scattering rates vary with local plasma properties motivated by extrinsic turbulence (ET) or self-confinement (SC) scenarios, with varying assumptions about e.g. turbulent power spectra on un-resolved scales, Alfven-wave damping, etc. We self-consistently predict observables including $γ$-rays ($L_γ$), grammage, residence times, and CR energy densities to constrain the models. We demonstrate many non-linear dynamical effects (not captured in simpler models) tend to enhance confinement. For example, in multi-phase media, even allowing arbitrary fast transport in neutral gas does not substantially reduce CR residence times (or $L_γ$), as transport is rate-limited by the ionized WIM and &#39;inner CGM&#39; gaseous halo ($10^{4}-10^{6}$ K gas within 10-30 kpc), and $L_γ$ can be dominated by trapping in small &#39;patches.&#39; Most physical ET models contribute negligible scattering of ~1-10 GeV CRs, but it is crucial to account for anisotropy and damping (especially of fast modes) or else scattering rates would violate observations. We show that the most widely-assumed scalings for SC models produce excessive confinement by factors >100 in the WIM and inner CGM, where turbulent and Landau damping dominate. This suggests either a breakdown of quasi-linear theory used to derive the CR transport parameters in SC, or that other novel damping mechanisms dominate in intermediate-density ionized gas.

preprint2022arXiv

The In-situ Origins of Dwarf Stellar Outskirts in FIRE-2

Extended, old, and round stellar halos appear to be ubiquitous around high-mass dwarf galaxies ($10^{8.5}<M_\star/M_\odot<10^{9.6}$) in the observed universe. However, it is unlikely that these dwarfs have undergone a sufficient number of minor mergers to form stellar halos that are composed of predominantly accreted stars. Here, we demonstrate that FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulations are capable of producing dwarf galaxies with realistic structure, including both a thick disk and round stellar halo. Crucially, these stellar halos are formed in-situ, largely via the outward migration of disk stars. However, there also exists a large population of &#34;non-disky&#34; dwarfs in FIRE-2 that lack a well-defined disk/halo and do not resemble the observed dwarf population. These non-disky dwarfs tend to be either more gas poor or to have burstier recent star formation histories than the disky dwarfs, suggesting that star formation feedback may be preventing disk formation. Both classes of dwarfs underscore the power of a galaxy&#39;s intrinsic shape -- which is a direct quantification of the distribution of the galaxy&#39;s stellar content -- to interrogate the feedback implementation in simulated galaxies.

preprint2021arXiv

A stripped-companion origin for Be stars: clues from the putative black holes HR 6819 and LB-1

HR 6819 is a bright ($V=5.36$), blue star recently proposed to be a triple containing a detached black hole (BH). We show that the system is a binary and does not contain a BH. Using spectral decomposition, we disentangle the observed composite spectra into two components: a rapidly rotating Be star and a slowly rotating B star with low surface gravity $(\log g \approx 2.75)$. Both stars show periodic radial velocity (RV) variability, but the RV semi-amplitude of the B star&#39;s orbit is $K_{\rm B}= (62.7 \pm 1)\,\rm km\,s^{-1}$, while that of the Be star is only $K_{\rm Be} = (4.5\pm 2)\,\rm km\,s^{-1}$. This implies that the B star is less massive by at least a factor of 10. The surface abundances of the B star bear imprints of CNO burning. We argue that the B star is a bloated, recently stripped helium star with mass $\approx 0.5\,M_{\odot}$ that is currently contracting to become a hot subdwarf. The orbital motion of the Be star obviates the need for a BH to explain the B star&#39;s motion. A stripped-star model reproduces the observed luminosity of the system, while a normal star with the B star&#39;s temperature and gravity would be more than 10 times too luminous. HR 6819 and the binary LB-1 probably formed through similar channels. We use MESA models to investigate their evolutionary history, finding that they likely formed from intermediate-mass ($3-7\,M_{\odot}$) primaries stripped by slightly lower-mass secondaries and are progenitors to Be + sdOB binaries such as $ϕ$ Persei. The lifetime of their current evolutionary phase is on average $2\times 10^5$ years, of order half a percent of the total lifetime of the Be phase. This implies that many Be stars have hot subdwarf and white dwarf companions, and that a substantial fraction ($20-100\%$) of field Be stars form through accretion of material from a binary companion.

preprint2021arXiv

Exploring the epoch of hydrogen reionization using FRBs

We describe three different methods for exploring the hydrogen reionization epoch using fast radio bursts (FRBs) and provide arguments for the existence of FRBs at high redshift (z). The simplest way, observationally, is to determine the maximum dispersion measure (DM$_{\rm max}$) of FRBs for an ensemble that includes bursts during the reionization. The DM$_{\rm max}$ provides information regarding reionization much like the optical depth of the CMB to Thomson scattering does, and it has the potential to be more accurate than constraints from Planck, if DM$_{\rm max}$ can be measured to a precision better than 500 $\mbox{pc cm}^{-3}$. Another method is to measure redshifts of about 40 FRBs between z of 6-10 with$\sim10\%$ accuracy to obtain the average electron density in 4 different z-bins with $\sim4\%$ accuracy. These two methods don&#39;t require knowledge of the FRB luminosity function and its possible redshift evolution. Finally, we show that the reionization history is reflected in the number of FRBs per unit DM, given a fluence limited survey of FRBs that includes bursts during the reionization epoch; we show using FIRE simulations that the contributions to DM from the FRB host galaxy $\&$ CGM during the reionization era is a small fraction of the observed DM. This third method requires no redshift information but does require knowledge of the FRB luminosity function.

preprint2021arXiv

Magnetically Modified Spherical Accretion in GRMHD: Reconnection-Driven Convection and Jet Propagation

We present 3D general relativistic magnetohydrodynamic(GRMHD) simulations of zero angular momentum accretion around a rapidly rotating black hole, modified by the presence of initially uniform magnetic fields. We consider serveral angles between the magnetic field direction and the black hole spin. In the resulting flows, the midplane dynamics are governed by magnetic reconnection-driven turbulence in a magnetically arrested (or a nearly arrested) state. Electromagnetic jets with outflow efficiencies ~10-200% occupy the polar regions, reaching several hundred gravitational radii before they dissipate due to the kink instability. The jet directions fluctuate in time and can be tilted by as much as ~30 degrees with respect to black hole spin, but this tilt does not depend strongly on the tilt of the initial magnetic field. A jet forms even when there is no initial net vertical magnetic flux since turbulent, horizon-scale fluctuations can generate a net vertical field locally. Peak jet power is obtained for an initial magnetic field tilted by 40-80 degrees with respect to the black hole spin because this maximizes the amount of magnetic flux that can reach the black hole. These simulations may be a reasonable model for low luminosity black hole accretion flows such as Sgr A* or M87.

preprint2021arXiv

The Effects of Tilt on the Time Variability of Millimeter and Infrared Emission from Sagittarius A*

Using a combination of general-relativistic magnetohydrodynamics simulations and ray tracing of synchrotron emission, we study the effect of modest (24 degrees) misalignment between the black hole spin and plasma angular momentum, focusing on the variability of total flux, image centroids, and image sizes. We consider both millimeter and infrared (IR) observables motivated by Sagittarius A* (Sgr A*), though our results apply more generally to optically thin flows. For most quantities, tilted accretion is more variable, primarily due to a significantly hotter and denser &#34;coronal&#34; region well off the disk midplane. We find (1) a 150% increase in millimeter light curve variability when adding tilt to the flow; (2) the tilted image centroid in the millimeter shifts on a scale of 3.7 microarcseconds over 28 hours (5000 gravitational times) for some electron temperature models; (3) tilted disk image diameters in the millimeter can be 10% larger (52 versus 47 microarcseconds) than those of aligned disks at certain viewing angles; (4) the tilted models produce significant IR flux, similar to that seen in Sgr A*, with comparable or even greater variability than observed; (5) for some electron models, the tilted IR centroid moves by more than 50 microarcseconds over several hours, in a similar fashion to the centroid motion detected by the GRAVITY interferometer.

preprint2021arXiv

The impact of astrophysical dust grains on the confinement of cosmic rays

We argue that charged dust grains could significantly impact the confinement and transport of galactic cosmic rays. For sub-GeV to ~1000GeV cosmic rays, small-scale parallel Alfvén waves, which isotropize cosmic rays through gyro-resonant interactions, are also gyro-resonant with charged grains. If the dust is nearly stationary, as in the bulk of the interstellar medium, Alfvén waves are damped by dust. This will reduce the amplitude of Alfvén waves produced by the cosmic rays through the streaming instability, thus enhancing cosmic-ray transport. In well-ionized regions, the dust damping rate is larger by a factor of ~10 than other mechanisms that damp parallel Alfvén waves at the scales relevant for ~GeV cosmic rays, suggesting that dust could play a key role in regulating cosmic-ray transport. In astrophysical situations in which the dust moves through the gas with super-Alfvénic velocities, Alfvén waves are rendered unstable, which could directly scatter cosmic rays. This interaction has the potential to create a strong feedback mechanism where dust, driven through the gas by radiation pressure, then strongly enhances the confinement of cosmic rays, increasing their capacity to drive outflows. This mechanism may act in the circumgalactic medium around star-forming galaxies and active galactic nuclei.

preprint2020arXiv

Ab Initio Horizon-Scale Simulations of Magnetically Arrested Accretion in Sagittarius A* Fed by Stellar Winds

We present 3D general relativistic magnetohydrodynamic (GRMHD) simulations of the accretion flow surrounding Sagittarius A* that are initialized using larger-scale MHD simulations of the $\sim$ 30 Wolf--Rayet (WR) stellar winds in the Galactic center. The properties of the resulting accretion flow on horizon scales are set not by ad hoc initial conditions but by the observationally constrained properties of the WR winds with limited free parameters. For this initial study we assume a non-spinning black hole. Our simulations naturally produce a $\sim 10^{-8} M_\odot$ yr$^{-1}$ accretion rate, consistent with previous phenomenological estimates. We find that a magnetically arrested flow is formed by the continuous accretion of coherent magnetic field being fed from large radii. Near the event horizon, the magnetic field is so strong that it tilts the gas with respect to the initial angular momentum and concentrates the originally quasi-spherical flow to a narrow disk-like structure. We also present 230 GHz images calculated from our simulations where the inclination angle and physical accretion rate are not free parameters but are determined by the properties of the WR stellar winds. The image morphology is highly time variable. Linear polarization on horizon scales is coherent with weak internal Faraday rotation.

preprint2020arXiv

Black widow evolution: magnetic braking by an ablated wind

Black widows are close binary systems in which a millisecond pulsar is orbited by a companion a few per cent the mass of the sun. It has been suggested that the pulsar&#39;s rotationally powered $γ$-ray luminosity gradually evaporates the companion, eventually leaving behind an isolated millisecond pulsar. The evaporation efficiency is determined by the temperature $T_{\rm ch}\propto F^{2/3}$ to which the outflow is heated by the flux $F$ on a dynamical time-scale. Evaporation is most efficient for companions that fill their Roche lobes. In this case, the outflow is dominated by a cap around the L1 point with an angle $θ_g\sim (T_{\rm ch}/T_g)^{1/2}$, and the evaporation time is $t_{\rm evap}=0.46(T_{\rm ch}/T_g)^{-2}\textrm{ Gyr}$, where $T_g>T_{\rm ch}$ is the companion&#39;s virial temperature. We apply our model to the observed black widow population, which has increased substantially over the last decade, considering each system&#39;s orbital period, companion mass, and pulsar spin-down power. While the original (Fruchter et al. 1988) black widow evaporates its companion on a few Gyr time-scale, direct evaporation on its own is too weak to explain the overall population. We propose instead that the evaporative wind couples to the companion&#39;s magnetic field, removes angular momentum from the binary, and maintains stable Roche-lobe overflow. While a stronger wind carries more mass, it also reduces the Alfvén radius, making this indirect magnetic braking mechanism less dependent on the flux $t_{\rm mag}\propto t_{\rm evap}^{1/3}$. This reduces the scatter in evolution times of observed systems, thus better explaining the combined black widow and isolated millisecond pulsar populations.

preprint2020arXiv

But What About... Cosmic Rays, Magnetic Fields, Conduction, & Viscosity in Galaxy Formation

We present a suite of high-resolution cosmological simulations, using the FIRE-2 feedback physics together with explicit treatment of magnetic fields, anisotropic conduction and viscosity, and cosmic rays (CRs) injected by supernovae (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultra-faint dwarf ($M_{\ast}\sim 10^{4}\,M_{\odot}$, $M_{\rm halo}\sim 10^{9}\,M_{\odot}$) through Milky Way masses, systematically vary CR parameters (e.g. the diffusion coefficient $κ$ and streaming velocity), and study an ensemble of galaxy properties (masses, star formation histories, mass profiles, phase structure, morphologies). We confirm previous conclusions that magnetic fields, conduction, and viscosity on resolved ($\gtrsim 1\,$pc) scales have small effects on bulk galaxy properties. CRs have relatively weak effects on all galaxy properties studied in dwarfs ($M_{\ast} \ll 10^{10}\,M_{\odot}$, $M_{\rm halo} \lesssim 10^{11}\,M_{\odot}$), or at high redshifts ($z\gtrsim 1-2$), for any physically-reasonable parameters. However at higher masses ($M_{\rm halo} \gtrsim 10^{11}\,M_{\odot}$) and $z\lesssim 1-2$, CRs can suppress star formation by factors $\sim 2-4$, given relatively high effective diffusion coefficients $κ\gtrsim 3\times10^{29}\,{\rm cm^{2}\,s^{-1}}$. At lower $κ$, CRs take too long to escape dense star-forming gas and lose energy to hadronic collisions, producing negligible effects on galaxies and violating empirical constraints from $γ$-ray emission. But around $κ\sim 3\times10^{29}\,{\rm cm^{2}\,s^{-1}}$, CRs escape the galaxy and build up a CR-pressure-dominated halo which supports dense, cool ($T\ll 10^{6}$ K) gas that would otherwise rain onto the galaxy. CR heating (from collisional and streaming losses) is never dominant.

preprint2020arXiv

No missing photons for reionization: moderate ionizing photon escape fractions from the FIRE-2 simulations

We present the escape fraction of hydrogen ionizing photons (f_esc) from a sample of 34 high-resolution cosmological zoom-in simulations of galaxies at z>5 in the Feedback in Realistic Environments project, post-processed with a Monte Carlo radiative transfer code for ionizing radiation. Our sample consists of 8500 halos in M_vir~10^8--10^{12} M_sun (M_star~10^4--10^{10} M_sun) at z=5--12. We find the sample average <f_esc> increases with halo mass for M_vir~10^8--10^{9.5} M_sun, becomes nearly constant for M_vir~10^{9.5}--10^{11} M_sun, and decreases at M_vir>10^{11} M_sun. Equivalently, <f_esc> increases with stellar mass up to M_star~10^8 M_sun and decreases at higher masses. Even applying single-star stellar population synthesis models, we find a moderate <f_esc>~0.2 for galaxies at M_star~10^8 M_sun. Nearly half of the escaped ionizing photons come from stars 1--3 Myr old and the rest from stars 3--10 Myr old. Binaries only have a modest effect, boosting <f_esc> by ~25--35% and the number of escaped photons by 60--80%. Most leaked ionizing photons are from vigorously star-forming regions that usually contain a feedback-driven kpc-scale superbubble surrounded by a dense shell. The shell is forming stars while accelerated, so new stars formed earlier in the shell are already inside the shell. Young stars in the bubble and near the edge of the shell can fully ionize some low-column-density paths pre-cleared by feedback, allowing a large fraction of their ionizing photons to escape. The decrease of <f_esc> at the high-mass end is due to dust attenuation, while at the low-mass end, <f_esc> decreases owing to inefficient star formation (and hence feedback). At fixed mass, <f_esc> tends to increase with redshift. Our simulations produce sufficient ionizing photons for cosmic reionization.

preprint2020arXiv

Not so fast: LB-1 is unlikely to contain a 70 $M_{\odot}$ black hole

The recently discovered binary LB-1 has been reported to contain a $\sim$\,$70\,M_{\odot}$ black hole (BH). The evidence for the unprecedentedly high mass of the unseen companion comes from reported radial velocity (RV) variability of the H$α$ emission line, which has been proposed to originate from an accretion disk around a BH. We show that there is in fact no evidence for RV variability of the H$α$ emission line, and that its apparent shifts instead originate from shifts in the luminous star&#39;s H$α$ absorption line. If not accounted for, such shifts will cause a stationary emission line to appear to shift in anti-phase with the luminous star. We show that once the template spectrum of a B star is subtracted from the observed Keck/HIRES spectra of LB-1, evidence for RV variability vanishes. Indeed, the data rule out periodic variability of the line with velocity semi-amplitude $K_{\rm Hα} > 1.3\,\rm km\,s^{-1}$. This strongly suggests that the observed H$α$ emission does not originate primarily from an accretion disk around a BH, and thus that the mass ratio cannot be constrained from the relative velocity amplitudes of the emission and absorption lines. The nature of the unseen companion remains uncertain, but a &#34;normal&#34; stellar-mass BH with mass $5\lesssim M/M_{\odot}\lesssim 20 $ seems most plausible. The H$α$ emission likely originates primarily from circumbinary material, not from either component of the binary.

preprint2020arXiv

On the Comparison of AGN with GRMHD Simulations: I. Sgr A*

We present models of Galactic Center emission in the vicinity of Sagittarius A* that use parametrizations of the electron temperature or energy density. These models include those inspired by two-temperature general relativistic magnetohydrodynamic (GRMHD) simulations as well as jet-motivated prescriptions generalizing equipartition of particle and magnetic energies. From these models, we calculate spectra and images and classify them according to their distinct observational features. Some models produce morphological and spectral features, e.g., image sizes, the sub-mm bump and low frequency spectral slope compatible with observations. Models with spectra consistent with observations produce the most compact images, with the most prominent, asymmetric photon rings. Limb brightened outflows are also visible in many models. Of all the models we consider, that which represents the current data the best is one in which electrons are relativistically hot when magnetic pressure is larger than the thermal pressure, but cold (i.e., negligibly contributing to the emission) otherwise. This work is part of a series also applying the &#34;observing&#34; simulations methodology to near-horizon regions of supermassive black holes in M87 and 3C 279.

preprint2020arXiv

Properties of the Circumgalactic Medium in Cosmic Ray-Dominated Galaxy Halos

We investigate the impact of cosmic rays (CRs) on the circumgalactic medium (CGM) in FIRE-2 simulations, for ultra-faint dwarf through Milky Way (MW)-mass halos hosting star-forming (SF) galaxies. Our CR treatment includes injection by supernovae, anisotropic streaming and diffusion along magnetic field lines, collisional and streaming losses, with constant parallel diffusivity $κ\sim3\times10^{29}\,\mathrm{cm^2\ s^{-1}}$ chosen to match $γ$-ray observations. With this, CRs become more important at larger halo masses and lower redshifts, and dominate the pressure in the CGM in MW-mass halos at $z\lesssim 1-2$. The gas in these &#34;CR-dominated&#34; halos differs significantly from runs without CRs: the gas is primarily cool (a few $\sim10^{4}\,$K), and the cool phase is volume-filling and has a thermal pressure below that needed for virial or local thermal pressure balance. Ionization of the &#34;low&#34; and &#34;mid&#34; ions in this diffuse cool gas is dominated by photo-ionization, with O VI columns $\gtrsim 10^{14.5}\,\mathrm{cm^{-2}}$ at distances $\gtrsim 150\,\mathrm{kpc}$. CR and thermal gas pressure are locally anti-correlated, maintaining total pressure balance, and the CGM gas density profile is determined by the balance of CR pressure gradients and gravity. Neglecting CRs, the same halos are primarily warm/hot ($T\gtrsim 10^{5}\,$K) with thermal pressure balancing gravity, collisional ionization dominates, O VI columns are lower and Ne VIII higher, and the cool phase is confined to dense filaments in local thermal pressure equilibrium with the hot phase.

preprint2020arXiv

Self-consistent proto-globular cluster formation in cosmological simulations of high-redshift galaxies

We report the formation of bound star clusters in a sample of high-resolution cosmological zoom-in simulations of z>5 galaxies from the FIRE project. We find that bound clusters preferentially form in high-pressure clouds with gas surface densities over 10^4 Msun pc^-2, where the cloud-scale star formation efficiency is near unity and young stars born in these regions are gravitationally bound at birth. These high-pressure clouds are compressed by feedback-driven winds and/or collisions of smaller clouds/gas streams in highly gas-rich, turbulent environments. The newly formed clusters follow a power-law mass function of dN/dM~M^-2. The cluster formation efficiency is similar across galaxies with stellar masses of ~10^7-10^10 Msun at z>5. The age spread of cluster stars is typically a few Myrs and increases with cluster mass. The metallicity dispersion of cluster members is ~0.08 dex in [Z/H] and does not depend on cluster mass significantly. Our findings support the scenario that present-day old globular clusters (GCs) were formed during relatively normal star formation in high-redshift galaxies. Simulations with a stricter/looser star formation model form a factor of a few more/fewer bound clusters per stellar mass formed, while the shape of the mass function is unchanged. Simulations with a lower local star formation efficiency form more stars in bound clusters. The simulated clusters are larger than observed GCs due to finite resolution. Our simulations are among the first cosmological simulations that form bound clusters self-consistently in a wide range of high-redshift galaxies.

preprint2020arXiv

Sound-Wave Instabilities in Dilute Plasmas with Cosmic Rays: Implications for Cosmic-Ray Confinement and the Perseus X-ray Ripples

Weakly collisional, magnetised plasmas characterised by anisotropic viscosity and conduction are ubiquitous in galaxies, halos and the intracluster medium (ICM). Cosmic rays (CRs) play an important role in these environments as well, by providing additional pressure and heating to the thermal plasma. We carry out a linear stability analysis of weakly collisional plasmas with cosmic rays using Braginskii MHD for the thermal gas. We assume that the CRs stream at the Alfvén speed, which in a weakly collisional plasma depends on the pressure anisotropy ($Δp$) of the thermal plasma. We find that this $Δp$-dependence introduces a phase shift between the CR-pressure and gas-density fluctuations. This drives a fast-growing acoustic instability: CRs offset the damping of acoustic waves by anisotropic viscosity and give rise to wave growth when the ratio of CR pressure to gas pressure is $\gtrsim αβ^{-1/2}$, where $β$ is the ratio of thermal to magnetic pressure, and $α$, typically $\lesssim 1$, depends on other dimensionless parameters. In high-$β$ environments like the ICM, this condition is satisfied for small CR pressures. We speculate that the instability studied here may contribute to the scattering of high-energy CRs and to the excitation of sound waves in galaxy-halo, group and cluster plasmas, including the long-wavelength X-ray fluctuations in \textit{Chandra} observations of the Perseus cluster. It may also be important in the vicinity of shocks in dilute plasmas (e.g., cluster virial shocks or galactic wind termination shocks), where the CR pressure is locally enhanced.

preprint2020arXiv

Synthetic Gaia surveys from the FIRE cosmological simulations of Milky Way-mass galaxies

With Gaia Data Release 2, the astronomical community is entering a new era of multidimensional surveys of the Milky Way. This new phase-space view of our Galaxy demands new tools for comparing observations to simulations of Milky-Way-mass galaxies in a cosmological context, to test the physics of both dark matter and galaxy formation. We present ananke, a framework for generating synthetic phase-space surveys from high-resolution baryonic simulations, and use it to generate a suite of synthetic surveys resembling Gaia DR2 in data structure, magnitude limits, and observational errors. We use three cosmological simulations of Milky-Way-mass galaxies from the Latte suite of the Feedback In Realistic Environments (FIRE) project, which feature self-consistent clustering of star formation in dense molecular clouds and thin stellar/gaseous disks in live cosmological halos with satellite dwarf galaxies and stellar halos. We select three solar viewpoints from each simulation to generate nine synthetic Gaia-like surveys. We sample synthetic stars by assuming each star particle (of mass 7070 $M_{\odot}$) represents a single stellar population. At each viewpoint, we compute dust extinction from the simulated gas metallicity distribution and apply a simple error model to produce a synthetic Gaia-like survey that includes both observational properties and a pointer to the generating star particle. We provide the complete simulation snapshot at $z = 0$ for each simulated galaxy. We describe data access points, the data model, and plans for future upgrades. These synthetic surveys provide a tool for the scientific community to test analysis methods and interpret Gaia data.

preprint2020arXiv

The Effect of Jet-Ejecta Interaction on the Viewing Angle Dependence of Kilonova Light Curves

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 multi-dimensional Monte Carlo radiation transport simulations of the resulting viewing angle-dependent light curves and spectra starting at 1.5h after merger. We find that on this timescale, jet shock heating does not affect the kilonova emission. 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 ($\lesssim 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.

preprint2020arXiv

The Effects of Tilt on the Images of Black Hole Accretion Flows

We analyze two 3D general-relativistic magnetohydrodynamic accretion simulations in the context of how they would manifest in Event Horizon Telescope (EHT) observations of supermassive black holes. The two simulations differ only in whether the initial angular momentum of the plasma is aligned with the rapid (a = 0.9) spin of the black hole. Both have low net magnetic flux. Ray tracing is employed to generate resolved images of the synchrotron emission. When using parameters appropriate for Sgr A* and assuming a viewing angle aligned with the black hole spin, we find the most prominent difference is that the central shadow in the image is noticeably eccentric in tilted models, with the ring of emission relatively unchanged. Applying this procedure to M87 with a viewing angle based on the large-scale jet, we find that adding tilt increases the angular size of the ring for fixed black hole mass and distance, while at the same time increasing the number of bright spots in the image. Our findings illustrate observable features that can distinguish tilted from aligned flows. They also show that tilted models can be viable for M87, and that not accounting for tilt can bias inferences of physical parameters. Future modeling of horizon-scale observations should account for potential angular momentum misalignment, which is likely generic at the low accretion rates appropriate for EHT targets.

preprint2020arXiv

The Impact of Type Ia Supernovae in Quiescent Galaxies: I. Formation of the Multiphase Interstellar medium

A cool phase of the interstellar medium has been observed in many giant elliptical galaxies, but its origin remains unclear. We propose that uneven heating from Type Ia supernovae (SNe Ia), together with radiative cooling, can lead to the formation of the cool phase. The basic idea is that since SNe Ia explode randomly, gas parcels which are not directly heated by SN shocks will cool, forming multiphase gas. We run a series of idealized high-resolution numerical simulations, and find that cool gas develops even when the overall SNe heating rate $H$ exceeds the cooling rate $C$ by a factor as large as 1.4. We also find that the time for multiphase gas development depends on the gas temperature. When the medium has a temperature $T = 3\times 10^6$ K, the cool phase forms within one cooling time \tc; however, the cool phase formation is delayed to a few times \tc\ for higher temperatures. The main reason for the delay is turbulent mixing. Cool gas formed this way would naturally have a metallicity lower than that of the hot medium. For constant $H/C$, there is more turbulent mixing for higher temperature gas. We note that this mechanism of producing cool gas cannot be captured in cosmological simulations, which usually fail to resolve individual SN remnants.

preprint2020arXiv

The Impact of Type Ia Supernovae in Quiescent Galaxies: II. Energetics and Turbulence

Type Ia supernovae (SNe Ia) provide unique and important feedback in quiescent galaxies, but their impact has been underappreciated. In this paper, we analyze a series of high-resolution simulations to examine the energetics and turbulence of the medium under SNe Ia. We find that when SN remnants are resolved, their effects differ distinctly from a volumetric heating term, as is commonly assumed in unresolved simulations. First, the net heating is significantly higher than expected, by 30$\pm$10\% per cooling time. This is because a large fraction of the medium is pushed into lower densities which cool inefficiently. Second, the medium is turbulent; the root-mean-squared (RMS) velocity of the gas to 20-50 km s$^{-1}$ on a driving scale of tens of parsec. The velocity field of the medium is dominated by compressional modes, which are larger than the solenoidal components by a factor of 3-7. Third, the hot gas has a very broad density distribution. The ratio between the density fluctuations and the RMS Mach number, parameterized as $b$, is 2-20. This is in contrast to previous simulations of turbulent media, which have found $b\lesssim$ 1. The reason for the difference is mainly caused by the \textit{localized} heating of SNe Ia, which creates a large density contrast. Last, the typical length scale of a density fluctuation grows with time, forming increasingly larger bubbles and filamentary ridges. These underlying density fluctuations need to be included when X-ray observations are interpreted.

preprint2020arXiv

The maximum accretion rate of hot gas in dark matter halos

We revisit the question of &#39;hot mode&#39; versus &#39;cold mode&#39; accretion onto galaxies using steady-state cooling flow solutions and idealized 3D hydrodynamic simulations. We demonstrate that for the hot accretion mode to exist, the cooling time is required to be longer than the free-fall time near the radius where the gas is rotationally-supported, R_circ, i.e. the existence of the hot mode depends on physical conditions at the galaxy scale rather than on physical conditions at the halo scale. When allowing for the depletion of the halo baryon fraction relative to the cosmic mean, the longer cooling times imply that a virialized gaseous halo may form in halo masses below the threshold of ~10^12 M_sun derived for baryon-complete halos. We show that for any halo mass there is a maximum accretion rate for which the gas is virialized throughout the halo and can accrete via the hot mode of Mdot_crit ~ 0.7(v_c/100 km/s)^5.4 (R_circ / 10 kpc) (Z/Z_sun)^-0.9 M_sun / yr, where Z and v_c are the metallicity and circular velocity measured at R_circ. For accretion rates >~Mdot_crit the volume-filling gas phase can in principle be `transonic&#39; -- virialized in the outer halo but cool and free-falling near the galaxy. We compare Mdot_crit to the average star formation rate (SFR) in halos at 0<z<10 implied by the stellar-mass halo-mass relation. For a plausible metallicity evolution with redshift, we find that SFR <~ Mdot_crit at most masses and redshifts, suggesting that the SFR of galaxies could be primarily sustained by the hot mode in halo masses well below the classic threshold of ~10^12 M_sun.

preprint2020arXiv

The Structure of Radiatively Inefficient Black Hole Accretion Flows

We run three long-timescale general-relativistic magnetohydrodynamic simulations of radiatively inefficient accretion flows onto non-rotating black holes. Our aim is to achieve steady-state behavior out to large radii and understand the resulting flow structure. A simulation with adiabatic index Gamma = 4/3 and small initial alternating poloidal magnetic field loops is run to a time of 440,000 GM/c^3, reaching inflow equilibrium inside a radius of 370 GM/c^2. Variations with larger alternating field loops and with Gamma = 5/3 are run to 220,000 GM/c^3, attaining equilibrium out to 170 GM/c^2 and 440 GM/c^2. There is no universal self-similar behavior obtained at radii in inflow equilibrium: the Gamma = 5/3 simulation shows a radial density profile with power law index ranging from -1 in the inner regions to -1/2 in the outer regions, while the others have a power-law slope ranging from -1/2 to close to -2. Both simulations with small field loops reach a state with polar inflow of matter, while the more ordered initial field has polar outflows. However, unbound outflows remove only a factor of order unity of the inflowing material over a factor of ~300 in radius. Our results suggest that the dynamics of radiatively inefficient accretion flows are sensitive to how the flow is fed from larger radii, and may differ appreciably in different astrophysical systems. Millimeter images appropriate for Sgr A* are qualitatively (but not quantitatively) similar in all simulations, with a prominent asymmetric image due to Doppler boosting.

preprint2020arXiv

The Surprisingly Small Impact of Magnetic Fields On The Inner Accretion Flow of Sagittarius A* Fueled By Stellar Winds

We study the flow structure in 3D magnetohydrodynamic (MHD) simulations of accretion onto Sagittarius A* via the magnetized winds of the orbiting Wolf-Rayet stars. These simulations cover over 3 orders of magnitude in radius to reach $\approx$ 300 gravitational radii, with only one poorly constrained parameter (the magnetic field in the stellar winds). Even for winds with relatively weak magnetic fields (e.g., plasma $β$ $\sim$ $10^6$), flux freezing/compression in the inflowing gas amplifies the field to $β$ $\sim$ few well before it reaches the event horizon. Overall, the dynamics, accretion rate, and spherically averaged flow profiles (e.g., density, velocity) in our MHD simulations are remarkably similar to analogous hydrodynamic simulations. We attribute this to the broad distribution of angular momentum provided by the stellar winds, which sources accretion even absent much angular momentum transport. We find that the magneto-rotational instability is not important because of i) strong magnetic fields that are amplified by flux freezing/compression, and ii) the rapid inflow/outflow times of the gas and inefficient radiative cooling preclude circularization. The primary effect of magnetic fields is that they drive a polar outflow that is absent in hydrodynamics. The dynamical state of the accretion flow found in our simulations is unlike the rotationally supported tori used as initial conditions in horizon scale simulations, which could have implications for models being used to interpret Event Horizon Telescope and GRAVITY observations of Sgr A*.

preprint2020arXiv

Thermal Instability of Halo Gas Heated by Streaming Cosmic Rays

Heating of virialized gas by streaming cosmic rays (CRs) may be energetically important in galaxy halos, groups and clusters. We present a linear thermal stability analysis of plasmas heated by streaming CRs. We separately treat equilibria with and without background gradients, and with and without gravity. We include both CR streaming and diffusion along the magnetic-field direction. Thermal stability depends strongly on the ratio of CR pressure to gas pressure, which determines whether modes are isobaric or isochoric. Modes with $\mathbf{k \cdot B }\neq 0$ are strongly affected by CR diffusion. When the streaming time is shorter than the CR diffusion time, thermally unstable modes (with $\mathbf{k \cdot B }\neq 0$) are waves propagating at a speed $\propto$ the Alfvén speed. Halo gas in photoionization equilibrium is thermally stable independent of CR pressure, while gas in collisional ionization equilibrium is unstable for physically realistic parameters. In gravitationally stratified plasmas, the oscillation frequency of thermally overstable modes can be higher in the presence of CR streaming than the buoyancy/free-fall frequency. This may modify the critical $t_{\rm cool}/t_{\rm ff}$ at which multiphase gas is present. The criterion for convective instability of a stratified, CR-heated medium can be written in the familiar Schwarzschild form $d s_{\rm eff} / d z < 0$, where $s_{\rm eff}$ is an effective entropy involving the gas and CR pressures. We discuss the implications of our results for the thermal evolution and multiphase structure of galaxy halos, groups and clusters.

preprint2019arXiv

Direct Detection of Black Hole-Driven Turbulence in the Centers of Galaxy Clusters

Supermassive black holes (SMBHs) are thought to provide energy that prevents catastrophic cooling in the centers of massive galaxies and galaxy clusters. However, it remains unclear how this &#34;feedback&#34; process operates. We use high-resolution optical data to study the kinematics of multi-phase filamentary structures by measuring the velocity structure function (VSF) of the filaments over a wide range of scales in the centers of three nearby galaxy clusters: Perseus, Abell 2597 and Virgo. We find that the motions of the filaments are turbulent in all three clusters studied. There is a clear correlation between features of the VSFs and the sizes of bubbles inflated by SMBH driven jets. Our study demonstrates that SMBHs are the main driver of turbulent gas motions in the centers of galaxy clusters and suggests that this turbulence is an important channel for coupling feedback to the environment. Our measured amplitude of turbulence is in good agreement with Hitomi Doppler line broadening measurement and X-ray surface brightness fluctuation analysis, suggesting that the motion of the cold filaments is well-coupled to that of the hot gas. The smallest scales we probe are comparable to the mean free path in the intracluster medium (ICM). Our direct detection of turbulence on these scales provides the clearest evidence to date that isotropic viscosity is suppressed in the weakly-collisional, magnetized intracluster plasma.

preprint2019arXiv

Low-frequency variability in massive stars: Core generation or surface phenomenon?

Bowman et al. (2019) reported low-frequency photometric variability in 164 O- and B-type stars observed with K2 and TESS. They interpret these motions as internal gravity waves, which could be excited stochastically by convection in the cores of these stars. The detection of internal gravity waves in massive stars would help distinguish between massive stars with convective or radiative cores, determine core size, and would provide important constraints on massive star structure and evolution. In this work, we study the observational signature of internal gravity waves generated by core convection. We calculate the \textit{wave transfer function}, which links the internal gravity wave amplitude at the base of the radiative zone to the surface luminosity variation. This transfer function varies by many orders of magnitude for frequencies $\lesssim 1 \, {\rm d}^{-1}$, and has regularly-spaced peaks near $1 \, {\rm d}^{-1}$ due to standing modes. This is inconsistent with the observed spectra which have smooth ``red noise&#39;&#39; profiles, without the predicted regularly-spaced peaks. The wave transfer function is only meaningful if the waves stay predominately linear. We next show that this is the case: low frequency traveling waves do not break unless their luminosity exceeds the radiative luminosity of the star; and, the observed luminosity fluctuations at high frequencies are so small that standing modes would be stable to nonlinear instability. These simple calculations suggest that the observed low-frequency photometric variability in massive stars is not due to internal gravity waves generated in the core of these stars. We finish with a discussion of (sub)surface convection, which produces low-frequency variability in low-mass stars, very similar to that observed in Bowman et al. (2019) in higher mass stars.

preprint2019arXiv

The Fate of AGB Wind in Massive Galaxies and the ICM

Asymptotic Giant Branch (AGB) winds from evolved stars not only provide a non-trivial amount of mass and energy return, but also produce dust grains in massive elliptical galaxies. Due to the fast stellar velocity and the high ambient temperature, the wind is thought to form a comet-like tail, similar to Mira in the Local Bubble. Many massive elliptical galaxies and cluster central galaxies host extended dusty cold filaments. The fate of the cold dusty stellar wind and its relation to cold filaments are not well understood. In this work, we carry out both analytical and numerical studies of the interaction between an AGB wind and the surrounding hot gas. We find that the cooling time of the tail is inversely proportional to the ambient pressure. In the absence of cooling, or in low pressure environments (e.g., the outskirts of elliptical galaxies), AGB winds are quickly mixed into the hot gas, and all the AGB winds have similar appearance and head-to-tail ratio. In high pressure environments, such as the Local Bubble and the central regions of massive elliptical galaxies, some of the gas in the mixing layer between the stellar wind and the surrounding hot gas can cool efficiently and cause the tail to become longer. Our simulated tail of Mira itself has similar length and velocity to that observed, and appears similar to the simulated AGB tail in the central regions of massive galaxies. We speculate that instead of thermal instability, the induced condensation at the mixing layer of AGB winds may be the origin of cold filaments in massive galaxies and galaxy clusters. This naturally explains the existence of dust and PAH in the filaments.

preprint2019arXiv

The Progenitors of Calcium-Strong Transients

A new class of faint, spectroscopically peculiar transients has emerged in the last decade. We term these events &#34;calcium-strong transients&#34; (CaSTs) because of their atypically high calcium-to-oxygen nebular line ratios. Previous studies have struggled to deduce the identity of their progenitors due to a combination of their extremely extended radial distributions with respect to their host galaxies and their relatively high rate of occurrence. In this work, we find that the CaST radial distribution is consistent with the radial distribution of two populations of stars: old (ages > 5 Gyr), low-metallicity (Z/Zsol < 0.3) stars and globular clusters. While no obvious progenitor scenario arises from considering old, metal-poor stars, the alternative production site of globular clusters leads us to narrow down the list of possible candidates to three binary scenarios: mergers of helium and oxygen/neon white dwarfs; tidal disruptions of helium white dwarfs by neutron stars; and stable accretion from low-mass helium-burning stars onto white dwarfs. While rare in the field, these binary systems can be formed dynamically at much higher rates in globular clusters. Subsequent binary hardening both increases their interaction rate and ejects them from their parent globular clusters prior to mass transfer contact. Their production in, and ejection from, globular clusters may explain their radial distribution and the absence of globular clusters at their explosion site. This model predicts a currently undiscovered high rate of CaSTs in nuclear star clusters. Alternatively, an undetermined progenitor scenario involving old, low-metallicity stars may instead hold the key to understanding CaSTs.