-
Pumping Iron: How turbulent metal diffusion impacts multiphase galactic outflows
Authors:
Ulrich P. Steinwandel,
Douglas Rennehan,
Matthew E. Orr,
Drummond B. Fielding,
Chang-Goo Kim
Abstract:
Most numerical simulations of galaxy formation and evolution are unable to properly resolve the turbulent cascade at or below the resolution scale and turbulence models are required to capture the motion of eddies on those unresolved scales. In this study, we investigate the impact of turbulent metal diffusion models on multiphase outflows originating from dwarf galaxies (…
▽ More
Most numerical simulations of galaxy formation and evolution are unable to properly resolve the turbulent cascade at or below the resolution scale and turbulence models are required to capture the motion of eddies on those unresolved scales. In this study, we investigate the impact of turbulent metal diffusion models on multiphase outflows originating from dwarf galaxies ($M_{\rm halo} \sim 10^{10} - 10^{11}$ M$_\odot$). We use our state-of-the-art numerical model for the formation of single stars and non-equilibrium cooling and hydrogen chemistry. Our simulations are carried out at a mass resolution of $\sim$1 M$_{\odot}$, where the individual supernova explosions are resolved in terms of hot-phase generation and momentum input. We find that mass, energy, and metal loading factors are only weakly affected by the inclusion of a metal diffusion model. The metal enrichment factor at low altitude above the galactic disk is higher by around 20 per cent when the metal diffusion model is included. Specifically, we find more efficient cooling in the cold interstellar medium, as higher amounts of metals are kept in the cold dense phase. The most striking effect of the metal diffusion model is that, without metal diffusion, there is more rapid cooling in the hot phase and a reduced sound speed by a factor of two. Specifically, we find that the hot phase is more metal enriched in the case without metal diffusion leading to more rapid (over) cooling of that phase which is consistent with the higher sound speed we find in the runs with metal diffusion.
△ Less
Submitted 19 July, 2024;
originally announced July 2024.
-
Equilibrium States of Galactic Atmospheres II: Interpretation and Implications
Authors:
G. M. Voit,
C. Carr,
D. B. Fielding,
V. Pandya,
G. L. Bryan,
M. Donahue,
B. D. Oppenheimer,
R. S. Somerville
Abstract:
The scaling of galaxy properties with halo mass suggests that feedback loops regulate star formation, but there is no consensus yet about how those feedback loops work. To help clarify discussions of galaxy-scale feedback, Paper I presented a very simple model for supernova feedback that it called the minimalist regulator model. This followup paper interprets that model and discusses its implicati…
▽ More
The scaling of galaxy properties with halo mass suggests that feedback loops regulate star formation, but there is no consensus yet about how those feedback loops work. To help clarify discussions of galaxy-scale feedback, Paper I presented a very simple model for supernova feedback that it called the minimalist regulator model. This followup paper interprets that model and discusses its implications. The model itself is an accounting system that tracks all of the mass and energy associated with a halo's circumgalactic baryons--the central galaxy's atmosphere. Algebraic solutions for the equilibrium states of that model reveal that star formation in low-mass halos self-regulates primarily by expanding the atmospheres of those halos, ultimately resulting in stellar masses that are insensitive to the mass-loading properties of galactic winds. What matters most is the proportion of supernova energy that couples with circumgalactic gas. However, supernova feedback alone fails to expand galactic atmospheres in higher-mass halos. According to the minimalist regulator model, an atmospheric contraction crisis ensues, which may be what triggers strong black-hole feedback. The model also predicts that circumgalactic medium properties emerging from cosmological simulations should depend largely on the specific energy of the outflows they produce, and we interpret the qualitative properties of several numerical simulations in light of that prediction.
△ Less
Submitted 11 June, 2024;
originally announced June 2024.
-
Equilibrium States of Galactic Atmospheres I: The Flip Side of Mass Loading
Authors:
G. M. Voit,
V. Pandya,
D. B. Fielding,
G. L. Bryan,
C. Carr,
M. Donahue,
B. D. Oppenheimer,
R. S. Somerville
Abstract:
This paper presents a new framework for understanding the relationship between a galaxy and its circumgalactic medium (CGM). It focuses on how imbalances between heating and cooling cause either expansion or contraction of the CGM. It does this by tracking \textit{all} of the mass and energy associated with a halo's baryons, including their gravitational potential energy, even if feedback has push…
▽ More
This paper presents a new framework for understanding the relationship between a galaxy and its circumgalactic medium (CGM). It focuses on how imbalances between heating and cooling cause either expansion or contraction of the CGM. It does this by tracking \textit{all} of the mass and energy associated with a halo's baryons, including their gravitational potential energy, even if feedback has pushed some of those baryons beyond the halo's virial radius. We show how a star-forming galaxy's equilibrium state can be algebraically derived within the context of this framework, and we analyze how the equilibrium star formation rate depends on supernova feedback. We consider the consequences of varying the mass loading parameter etaM = Mdot_wind / Mdot_* relating a galaxy's gas mass outflow rate (Mdot_wind) to its star formation rate (Mdot_*) and obtain results that challenge common assumptions. In particular, we find that equilibrium star formation rates in low-mass galaxies are generally insensitive to mass loading, and when mass loading does matter, increasing it actually results in \textit{more} star formation because more supernova energy is needed to resist atmospheric contraction.
△ Less
Submitted 11 June, 2024;
originally announced June 2024.
-
Metallicity Dependence of Pressure-Regulated Feedback-Modulated Star Formation in the TIGRESS-NCR Simulation Suite
Authors:
Chang-Goo Kim,
Eve C. Ostriker,
Jeong-Gyu Kim,
Munan Gong,
Greg L. Bryan,
Drummond B. Fielding,
Sultan Hassan,
Matthew Ho,
Sarah M. R. Jeffreson,
Rachel S. Somerville,
Ulrich P. Steinwandel
Abstract:
We present a new simulation suite for the star-forming interstellar medium (ISM) in galactic disks using the TIGRESS-NCR framework. Distinctive aspects of our simulation suite are: (1) sophisticated and comprehensive numerical treatments of essential physical processes including magnetohydrodynamics, self-gravity, and galactic differential rotation, as well as photochemistry, cooling, and heating…
▽ More
We present a new simulation suite for the star-forming interstellar medium (ISM) in galactic disks using the TIGRESS-NCR framework. Distinctive aspects of our simulation suite are: (1) sophisticated and comprehensive numerical treatments of essential physical processes including magnetohydrodynamics, self-gravity, and galactic differential rotation, as well as photochemistry, cooling, and heating coupled with ray-tracing UV radiation transfer and resolved supernova feedback and (2) wide parameter coverage including metallicity over $Z'\equiv Z/Z_\odot\sim0.1-3$, gas surface density $Σ_{\rm gas}\sim5-150 M_{\odot}{\rm pc^{-2}}$, and stellar surface density $Σ_{\rm star}\sim 1-50 M_{\odot}{\rm pc^{-2}}$. The range of emergent star formation rate surface density is $Σ_{\rm SFR}\sim 10^{-4}-0.5 M_{\odot}{\rm kpc^{-2}yr^{-1}}$ and ISM total midplane pressure is $P_{\rm tot}/k_B=10^3-10^6{\rm cm^{-3}K}$, with $P_{\rm tot}$ equal to the ISM weight $W$. For given $Σ_{\rm gas}$ and $Σ_{\rm star}$, we find $Σ_{\rm SFR} \propto Z'^{0.3}$. We provide an interpretation based on the pressure-regulated feedback-modulated (PRFM) star formation theory. We characterize feedback modulation in terms of the yield $Υ$, defined as the ratio of each stress to $Σ_{\rm SFR}$. The thermal feedback yield varies sensitively with both weight and metallicity as $Υ_{\rm th}\propto W^{-0.46}Z'^{-0.53}$, while the combined turbulent and magnetic feedback yield shows weaker dependence $Υ_{\rm turb+mag}\propto W^{-0.22}Z'^{-0.18}$. The reduction in $Σ_{\rm SFR}$ at low metallicity is due mainly to enhanced thermal feedback yield, resulting from reduced attenuation of UV radiation. With the metallicity-dependent calibrations we provide, PRFM theory can be used for a new subgrid star formation prescription in cosmological simulations where the ISM is unresolved.
△ Less
Submitted 6 June, 2024; v1 submitted 29 May, 2024;
originally announced May 2024.
-
DUVET: Resolved direct metallicity measurements in the outflow of starburst galaxy NGC 1569
Authors:
Magdalena J. Hamel-Bravo,
Deanne B. Fisher,
Danielle Berg,
Bjarki Björgvinsson,
Alberto D. Bolatto,
Alex J. Cameron,
John Chisholm,
Drummond B. Fielding,
Rodrigo Herrera-Camus,
Glenn G. Kacprzak,
Miao Li,
Barbara Mazzilli Ciraulo,
Anna F. McLeod,
Daniel K. McPherson,
Nikole M. Nielsen,
Bronwyn Reichardt Chu,
Ryan J. Rickards Vaught,
Karin Sandstrom
Abstract:
We present the results of direct-method metallicity measurements in the disk and outflow of the low-metallicity starburst galaxy NGC 1569. We use Keck Cosmic Web Imager observations to map the galaxy across 54$\arcsec$ (800 pc) along the major axis and 48$\arcsec$ (700 pc) along the minor axis with a spatial resolution of 1$\arcsec$ ($\sim$15 pc). We detect common strong emission lines ([\ion{O}{I…
▽ More
We present the results of direct-method metallicity measurements in the disk and outflow of the low-metallicity starburst galaxy NGC 1569. We use Keck Cosmic Web Imager observations to map the galaxy across 54$\arcsec$ (800 pc) along the major axis and 48$\arcsec$ (700 pc) along the minor axis with a spatial resolution of 1$\arcsec$ ($\sim$15 pc). We detect common strong emission lines ([\ion{O}{III}] $λ$5007, H$β$, [\ion{O}{II}] $λ$3727) and the fainter [\ion{O}{III}] $λ$4363 auroral line, which allows us to measure electron temperature ($T_e$) and metallicity. Theory suggests that outflows drive metals out of the disk driving observed trends between stellar mass and gas-phase metallicity. Our main result is that the metallicity in the outflow is similar to that of the disk, $Z_{\rm out} / Z_{\rm ISM} \approx 1$. This is consistent with previous absorption line studies in higher mass galaxies. Assumption of a mass-loading factor of $\dot{M}_{\rm out}/{\rm SFR}\sim3$ makes the metal-loading of NGC 1569 consistent with expectations derived from the mass-metallicity relationship. Our high spatial resolution metallicity maps reveal a region around a supermassive star cluster (SSC-B) with distinctly higher metallicity and higher electron density, compared to the disk. Given the known properties of SSC-B the higher metallicity and density of this region are likely the result of star formation-driven feedback acting on the local scale. Overall, our results are consistent with the picture in which metal-enriched winds pollute the circumgalactic medium surrounding galaxies, and thus connect the small-scale feedback processes to large-scale properties of galaxy halos.
△ Less
Submitted 6 April, 2024;
originally announced April 2024.
-
DUVET: sub-kiloparsec resolved star formation driven outflows in a sample of local starbursting disk galaxies
Authors:
Bronwyn Reichardt Chu,
Deanne B. Fisher,
John Chisholm,
Danielle Berg,
Alberto Bolatto,
Alex J. Cameron,
Drummond B. Fielding,
Rodrigo Herrera-Camus,
Glenn G. Kacprzak,
Miao Li,
Anna F. McLeod,
Daniel K. McPherson,
Nikole M. Nielsen,
Ryan Rickards Vaught,
Sophia G. Ridolfo,
Karin Sandstrom
Abstract:
We measure resolved (kiloparsec-scale) outflow properties in a sample of 10 starburst galaxies from the DUVET sample, using Keck/KCWI observations of H$β$ and [OIII]~$λ$5007. We measure $\sim450$ lines-of-sight that contain outflows, and use these to study scaling relationships of outflow velocity ($v_{\rm out}$), mass-loading factor ($η$; mass outflow rate per SFR) and mass flux (…
▽ More
We measure resolved (kiloparsec-scale) outflow properties in a sample of 10 starburst galaxies from the DUVET sample, using Keck/KCWI observations of H$β$ and [OIII]~$λ$5007. We measure $\sim450$ lines-of-sight that contain outflows, and use these to study scaling relationships of outflow velocity ($v_{\rm out}$), mass-loading factor ($η$; mass outflow rate per SFR) and mass flux ($\dotΣ_{\rm out}$; mass outflow rate per area) with co-located SFR surface density ($Σ_{\rm SFR}$) and stellar mass surface density ($Σ_{\ast}$). We find strong, positive correlations of $\dotΣ_{\rm out} \propto Σ_{\rm SFR}^{1.2}$ and $\dotΣ_{\rm out} \propto Σ_{\ast}^{1.7}$. We also find shallow correlations between $v_{\rm out}$ and both $Σ_{\rm SFR}$ and $Σ_{\ast}$. Our resolved observations do not suggest a threshold in outflows with $Σ_{\rm SFR}$, but rather we find that the local specific SFR ($Σ_{\rm SFR}/Σ_\ast$) is a better predictor of where outflows are detected. We find that outflows are very common above $Σ_{\rm SFR}/Σ_\ast\gtrsim 0.1$~Gyr$^{-1}$ and rare below this value. We argue that our results are consistent with a picture in which outflows are driven by supernovae, and require more significant injected energy in higher mass surface density environments to overcome local gravity. The correlations we present here provide a statistically robust, direct comparison for simulations and higher redshift results from JWST.
△ Less
Submitted 27 February, 2024;
originally announced February 2024.
-
CloudFlex: A Flexible Parametric Model for the Small-Scale Structure of the Circumgalactic Medium
Authors:
Cameron B. Hummels,
Kate H. R. Rubin,
Evan E. Schneider,
Drummond B. Fielding
Abstract:
We present CloudFlex, a new open-source tool for predicting the absorption-line signatures of cool gas in galaxy halos with complex small-scale structure. Motivated by analyses of cool material in hydrodynamical simulations of turbulent, multiphase media, we model individual cool gas structures as assemblies of cloudlets with a power-law distribution of cloudlet mass $\propto m_{\rm cl}^{-α}$ and…
▽ More
We present CloudFlex, a new open-source tool for predicting the absorption-line signatures of cool gas in galaxy halos with complex small-scale structure. Motivated by analyses of cool material in hydrodynamical simulations of turbulent, multiphase media, we model individual cool gas structures as assemblies of cloudlets with a power-law distribution of cloudlet mass $\propto m_{\rm cl}^{-α}$ and relative velocities drawn from a turbulent velocity field. The user may specify $α$, the lower limit of the cloudlet mass distribution ($m_{\rm cl,min}$), and several other parameters that set the total mass, size, and velocity distribution of the complex. We then calculate the MgII 2796 absorption profiles induced by the cloudlets along pencil-beam lines of sight. We demonstrate that at fixed metallicity, the covering fraction of sightlines with equivalent widths $W_{2796} < 0.3$ Ang increases significantly with decreasing $m_{\rm cl,min}$, cool cloudlet number density ($n_{\rm cl}$), and cloudlet complex size. We then present a first application, using this framework to predict the projected $W_{2796}$ distribution around ${\sim}L^*$ galaxies. We show that the observed incidences of $W_{2796}>0.3$ Ang sightlines within 10 kpc < $R_{\perp}$ < 50 kpc are consistent with our model over much of parameter space. However, they are underpredicted by models with $m_{\rm cl,min}\ge100M_{\odot}$ and $n_{\rm cl}\ge0.03$ $\rm cm^{-3}$, in keeping with a picture in which the inner cool circumgalactic medium (CGM) is dominated by numerous low-mass cloudlets ($m_{\rm cl}\lesssim100M_{\odot}$) with a volume filling factor ${\lesssim}1\%$. When used to simultaneously model absorption-line datasets built from multi-sightline and/or spatially-extended background probes, CloudFlex will enable detailed constraints on the size and velocity distributions of structures comprising the photoionized CGM.
△ Less
Submitted 9 November, 2023;
originally announced November 2023.
-
Seeding the CGM: How Satellites Populate the Cold Phase of Milky Way Halos
Authors:
Manami Roy,
Kung-Yi Su,
Stephanie Tonnesen,
Drummond B. Fielding,
Claude-André Faucher-Giguère
Abstract:
The origin of the cold phase in the CGM is a highly debated question. We investigate the contribution of satellite galaxies to the cold gas budget in the circumgalactic medium (CGM)of a Milky Way-like host galaxy. We perform controlled experiments with three different satellite mass distributions and identify several mechanisms by which satellites can add cold gas to the CGM, including ram pressur…
▽ More
The origin of the cold phase in the CGM is a highly debated question. We investigate the contribution of satellite galaxies to the cold gas budget in the circumgalactic medium (CGM)of a Milky Way-like host galaxy. We perform controlled experiments with three different satellite mass distributions and identify several mechanisms by which satellites can add cold gas to the CGM, including ram pressure stripping and induced cooling in the mixing layer of the stripped cold gas. These two mechanisms contribute a comparable amount of cold gas to the host CGM. We find that the less massive satellites ($\leq 10^9 M_\odot$) not only lose all of their cold gas in a short period ($\sim$ 0.5-1 Gyr), but their stripped cold clouds also mix with the hot CGM gas and get heated up quickly. However, stellar feedback from these less massive satellites can hugely alter the fate of their stripped gas. Feedback speeds up the destruction of the stripped cold clouds from these satellites by making them more diffuse with more surface area. On the other hand, the more massive satellites (LMC or SMC-like $\sim 10^{10} M_\odot$) can add cold gas to the total gas budget of the host CGM for several Gyrs.
△ Less
Submitted 6 October, 2023;
originally announced October 2023.
-
DUVET Survey: Mapping Outflows in the Metal-Poor Starburst Mrk 1486
Authors:
Daniel K. McPherson,
Deanne B. Fisher,
Nikole M. Nielsen,
Glenn G. Kacprzak,
Bronwyn Reichardt Chu,
Alex J. Cameron,
Alberto D. Bolatto,
John Chisholm,
Drummond B. Fielding,
Danielle Berg,
Rodrigo Herrera-Camus,
Miao Li,
Ryan J. Rickards Vaught,
Karin Sandstrom
Abstract:
We present a method to characterize star-formation driven outflows from edge-on galaxies and apply this method to the metal-poor starburst galaxy, Mrk 1486. Our method uses the distribution of emission line flux (from H$β$ and [OIII] 5007) to identify the location of the outflow and measure the extent above the disk, the opening angle, and the transverse kinematics. We show that this simple techni…
▽ More
We present a method to characterize star-formation driven outflows from edge-on galaxies and apply this method to the metal-poor starburst galaxy, Mrk 1486. Our method uses the distribution of emission line flux (from H$β$ and [OIII] 5007) to identify the location of the outflow and measure the extent above the disk, the opening angle, and the transverse kinematics. We show that this simple technique recovers a similar distribution of the outflow without requiring complex modelling of line-splitting or multi-Gaussian components, and is therefore applicable to lower spectral resolution data. In Mrk 1486 we observe an asymmetric outflow in both the location of the peak flux and total flux from each lobe. We estimate an opening angle of $17-37^{\circ}$ depending on the method and assumptions adopted. Within the minor axis outflows, we estimate a total mass outflow rate of $\sim2.5$ M$_{\odot}$ yr$^{-1}$, which corresponds to a mass loading factor of $η=0.7$. We observe a non-negligible amount of flux from ionized gas outflowing along the edge of the disk (perpendicular to the biconical components), with a mass outflow rate $\sim0.9$ M$_{\odot}$ yr$^{-1}$. Our results are intended to demonstrate a method that can be applied to high-throughput, low spectral resolution observations, such as narrow band filters or low spectral resolution IFS that may be more able to recover the faint emission from outflows.
△ Less
Submitted 13 August, 2023;
originally announced August 2023.
-
TuRMoiL of Survival: A Unified Survival Criterion for Cloud-Wind Interactions
Authors:
Matthew W. Abruzzo,
Drummond B. Fielding,
Greg L. Bryan
Abstract:
Cloud-wind interactions play an important role in long-lived multiphase flows in many astrophysical contexts. When this interaction is primarily mediated by hydrodynamics and radiative cooling, the survival of clouds can be phrased in terms of the comparison between a timescale that dictates the evolution of the cloud-wind interaction, (the dynamical time-scale $τ_{\rm dyn}$) and the relevant cool…
▽ More
Cloud-wind interactions play an important role in long-lived multiphase flows in many astrophysical contexts. When this interaction is primarily mediated by hydrodynamics and radiative cooling, the survival of clouds can be phrased in terms of the comparison between a timescale that dictates the evolution of the cloud-wind interaction, (the dynamical time-scale $τ_{\rm dyn}$) and the relevant cooling timescale $τ_{\rm cool}$. Previously proposed survival criteria, which can disagree by large factors about the size of the smallest surviving clouds, differ in both their choice of $τ_{\rm cool}$ and (to a lesser extent) $τ_{\rm dyn}$. Here we present a new criterion which agrees with a previously proposed empirical formulae but is based on simple physical principles. The key insight is that clouds can grow if they are able to mix and cool gas from the hot wind faster than it advects by the cloud. Whereas prior criteria associate $τ_{\rm dyn}$ with the cloud crushing timescale, our new criterion links it to the characteristic cloud-crossing timescale of a hot-phase fluid element, making it more physically consistent with shear-layer studies. We develop this insight into a predictive expression and validate it with hydrodynamic ENZO-E simulations of ${\sim}10^4\, {\rm K}$, pressure-confined clouds in hot supersonic winds, exploring, in particular, high wind/cloud density contrasts, where disagreements are most pronounced. Finally, we illustrate how discrepancies among previous criteria primarily emerged due to different choices of simulation conditions and cooling properties, and discuss how they can be reconciled.
△ Less
Submitted 6 July, 2023;
originally announced July 2023.
-
Cloud Atlas: Navigating the Multiphase Landscape of Tempestuous Galactic Winds
Authors:
Brent Tan,
Drummond B. Fielding
Abstract:
Galaxies comprise intricate networks of interdependent processes which together govern their evolution. Central among these are the multiplicity of feedback channels, which remain incompletely understood. One outstanding problem is the understanding and modeling of the multiphase nature of galactic winds, which play a crucial role in galaxy formation and evolution. We present the results of three…
▽ More
Galaxies comprise intricate networks of interdependent processes which together govern their evolution. Central among these are the multiplicity of feedback channels, which remain incompletely understood. One outstanding problem is the understanding and modeling of the multiphase nature of galactic winds, which play a crucial role in galaxy formation and evolution. We present the results of three dimensional magnetohydrodynamical tall box interstellar medium patch simulations with clustered supernova driven outflows. Fragmentation of the interstellar medium during superbubble breakout seeds the resulting hot outflow with a population of cool clouds. We focus on analyzing and modeling the origin and properties of these clouds. Their presence induces large scale turbulence, which in turn leads to complex cloud morphologies. Cloud sizes are well described by a power law distribution and mass growth rates can be modelled using turbulent radiative mixing layer theory. Turbulence provides significant pressure support in the clouds, while magnetic fields only play a minor role. We conclude that many of the physical insights and analytic scalings derived from idealized small scale simulations translate well to larger scale, more realistic turbulent magnetized winds, thus paving a path towards their necessary yet challenging inclusion in global-scale galaxy models.
△ Less
Submitted 7 December, 2023; v1 submitted 23 May, 2023;
originally announced May 2023.
-
Cosmic ray transport in large-amplitude turbulence with small-scale field reversals
Authors:
Philipp Kempski,
Drummond B. Fielding,
Eliot Quataert,
Alisa K. Galishnikova,
Matthew W. Kunz,
Alexander A. Philippov,
Bart Ripperda
Abstract:
The nature of cosmic ray (CR) transport in the Milky Way remains elusive. The predictions of current micro-physical CR transport models in magneto-hydrodynamic (MHD) turbulence are drastically different from what is observed. These models usually focus on MHD turbulence with a strong guide field and ignore the impact of turbulent intermittency on particle propagation. This motivates our studying t…
▽ More
The nature of cosmic ray (CR) transport in the Milky Way remains elusive. The predictions of current micro-physical CR transport models in magneto-hydrodynamic (MHD) turbulence are drastically different from what is observed. These models usually focus on MHD turbulence with a strong guide field and ignore the impact of turbulent intermittency on particle propagation. This motivates our studying the alternative regime of large-amplitude turbulence with $δB/B_0 \gg 1$, in which intermittent small-scale magnetic field reversals are ubiquitous. We study particle transport in such turbulence by integrating trajectories in stationary snapshots. To quantify spatial diffusion, we use a setup with continuous particle injection and escape, which we term the turbulent leaky box. We find that particle transport is very different from the strong-guide-field case. Low-energy particles are better confined than high-energy particles, despite less efficient pitch-angle isotropization at small energies. In the limit of weak guide field, energy-dependent confinement is driven by the energy-dependent (in)ability to follow reversing magnetic field lines exactly and by the scattering in regions of ``resonant curvature", where the field line bends on a scale that is of order the local particle gyro-radius. We derive a heuristic model of particle transport in magnetic folds that approximately reproduces the energy dependence of transport found numerically. We speculate that CR propagation in the Galaxy is regulated by the intermittent field reversals highlighted here and discuss the implications of our findings for CR transport in the Milky Way.
△ Less
Submitted 19 December, 2023; v1 submitted 24 April, 2023;
originally announced April 2023.
-
Arkenstone I: A Novel Method for Robustly Capturing High Specific Energy Outflows In Cosmological Simulations
Authors:
Matthew C. Smith,
Drummond B. Fielding,
Greg L. Bryan,
Chang-Goo Kim,
Eve C. Ostriker,
Rachel S. Somerville,
Jonathan Stern,
Kung-Yi Su,
Rainer Weinberger,
Chia-Yu Hu,
John C. Forbes,
Lars Hernquist,
Blakesley Burkhart,
Yuan Li
Abstract:
Arkenstone is a new model for multiphase, stellar feedback driven galactic winds designed for inclusion in coarse resolution cosmological simulations. In this first paper of a series, we describe the features that allow Arkenstone to properly treat high specific energy wind components and demonstrate them using idealised non-cosmological simulations of a galaxy with a realistic CGM, using the Arep…
▽ More
Arkenstone is a new model for multiphase, stellar feedback driven galactic winds designed for inclusion in coarse resolution cosmological simulations. In this first paper of a series, we describe the features that allow Arkenstone to properly treat high specific energy wind components and demonstrate them using idealised non-cosmological simulations of a galaxy with a realistic CGM, using the Arepo code. Hot, fast gas phases with low mass loadings are predicted to dominate the energy content of multiphase outflows. In order to treat the huge dynamic range of spatial scales involved in cosmological galaxy formation at feasible computational expense, cosmological volume simulations typically employ a Lagrangian code or else use adaptive mesh refinement with a quasi-Lagrangian refinement strategy. However, it is difficult to inject a high specific energy wind in a Lagrangian scheme without incurring artificial burstiness. Additionally, the low densities inherent to this type of flow result in poor spatial resolution. Arkenstone addresses these issues with a novel scheme for coupling energy into the ISM/CGM transition region which also provides the necessary level of refinement at the base of the wind. In the absence of our improvements, we show that poor spatial resolution near the sonic point of a hot, fast outflow leads to an underestimation of gas acceleration as the wind propagates. We explore the different mechanisms by which low and high specific energy winds can regulate the SFR of galaxies. In future work, we will demonstrate other aspects of the Arkenstone model.
△ Less
Submitted 7 March, 2024; v1 submitted 17 January, 2023;
originally announced January 2023.
-
The structure and composition of multiphase galactic winds in a Large Magellanic Cloud mass simulated galaxy
Authors:
Ulrich P. Steinwandel,
Chang-Goo Kim,
Greg L. Bryan,
Eve C. Ostriker,
Rachel S. Somerville,
Drummond B. Fielding
Abstract:
We present the first results from a high resolution simulation with a focus on galactic wind driving for an isolated galaxy with a halo mass of $\sim 10^{11}$ M$_{\odot}$ (similar to the Large Magellanic Cloud) and a total gas mass of $\sim 6 \times 10^{8}$ M$_{\odot}$, resulting in $\sim 10^{8}$ gas cells at $\sim 4$ M$_{\odot}$ mass resolution. We adopt a resolved stellar feedback model with non…
▽ More
We present the first results from a high resolution simulation with a focus on galactic wind driving for an isolated galaxy with a halo mass of $\sim 10^{11}$ M$_{\odot}$ (similar to the Large Magellanic Cloud) and a total gas mass of $\sim 6 \times 10^{8}$ M$_{\odot}$, resulting in $\sim 10^{8}$ gas cells at $\sim 4$ M$_{\odot}$ mass resolution. We adopt a resolved stellar feedback model with non-equilibrium cooling and heating, including photoelectric heating and photo-ionizing radiation, as well as supernovae (SNe), coupled to the second order meshless finite mass (MFM) method for hydrodynamics. These features make this the largest resolved-ISM galaxy model run to date. We find mean star formation rates around $0.05$ M$_{\odot}$ yr$^{-1}$ and evaluate typical time averaged loading factors for mass ($η_\mathrm{M}$ $\sim$ 1.0, in good agreement with recent observations) and energy ($η_\mathrm{E}$ $\sim$ 0.01). The bulk of the mass of the wind is transported by the warm ($T < 5 \times 10^5$K) phase, while there is a similar amount of energy transported in the warm and the hot phases ($T > 5 \times 10^5$K). We find an average opening angle of 30 degrees for the wind, decreasing with higher altitude above the midplane. The wind mass loading is decreasing (flat) for the warm (hot) phase as a function of the star formation surface rate density $Σ_{\rm SFR}$, while the energy loading shows inverted trends with $Σ_{\rm SFR}$, decreasing for the warm wind and increasing for the hot wind, although with very shallow slopes. These scalings are in good agreement with previous simulations of resolved wind driving in the multi-phase ISM.
△ Less
Submitted 7 December, 2022;
originally announced December 2022.
-
A unified model for the co-evolution of galaxies and their circumgalactic medium: the relative roles of turbulence and atomic cooling physics
Authors:
Viraj Pandya,
Drummond B. Fielding,
Greg L. Bryan,
Christopher Carr,
Rachel S. Somerville,
Jonathan Stern,
Claude-Andre Faucher-Giguere,
Zachary Hafen,
Daniel Angles-Alcazar,
John C. Forbes
Abstract:
The circumgalactic medium (CGM) plays a pivotal role in regulating gas flows around galaxies and thus shapes their evolution. However, the details of how galaxies and their CGM co-evolve remain poorly understood. We present a new time-dependent two-zone model that self-consistently tracks not just mass and metal flows between galaxies and their CGM but also the evolution of the global thermal and…
▽ More
The circumgalactic medium (CGM) plays a pivotal role in regulating gas flows around galaxies and thus shapes their evolution. However, the details of how galaxies and their CGM co-evolve remain poorly understood. We present a new time-dependent two-zone model that self-consistently tracks not just mass and metal flows between galaxies and their CGM but also the evolution of the global thermal and turbulent kinetic energy of the CGM. Our model accounts for heating and turbulence driven by both supernova winds and cosmic accretion as well as radiative cooling, turbulence dissipation, and halo outflows due to CGM overpressurization. We demonstrate that, depending on parameters, the CGM can undergo a phase transition (``thermalization'') from a cool, turbulence-supported phase to a virial-temperature, thermally-supported phase. This CGM phase transition is largely determined by the ability of radiative cooling to balance heating from supernova winds and turbulence dissipation. We perform an initial calibration of our model to the FIRE-2 cosmological hydrodynamical simulations and show that it can approximately reproduce the baryon cycles of the simulated halos. In particular, we find that, for these parameters, the phase transition occurs at high-redshift in ultrafaint progenitors and at low redshift in classical $M_{\rm vir}\sim10^{11}M_{\odot}$ dwarfs, while Milky Way-mass halos undergo the transition at $z\approx0.5$. We see a similar transition in the simulations though it is more gradual, likely reflecting radial dependence and multi-phase gas not captured by our model. We discuss these and other limitations of the model and possible future extensions.
△ Less
Submitted 27 August, 2023; v1 submitted 17 November, 2022;
originally announced November 2022.
-
Plasmoid Instability in the Multiphase Interstellar Medium
Authors:
Drummond B. Fielding,
Bart Ripperda,
Alexander A. Philippov
Abstract:
The processes controlling the complex clump structure, phase distribution, and magnetic field geometry that develops across a broad range of scales in the turbulent interstellar medium remains unclear. Using unprecedentedly high resolution three-dimensional magnetohydrodynamic simulations of thermally unstable turbulent systems, we show that large current sheets unstable to plasmoid-mediated recon…
▽ More
The processes controlling the complex clump structure, phase distribution, and magnetic field geometry that develops across a broad range of scales in the turbulent interstellar medium remains unclear. Using unprecedentedly high resolution three-dimensional magnetohydrodynamic simulations of thermally unstable turbulent systems, we show that large current sheets unstable to plasmoid-mediated reconnection form regularly throughout the volume. The plasmoids form in three distinct environments: (i) within cold clumps, (ii) at the asymmetric interface of the cold and warm phases, and (iii) within the warm, volume-filling phase. We then show that the complex magneto-thermal phase structure is characterized by a predominantly highly magnetized cold phase, but that regions of high magnetic curvature, which are the sites of reconnection, span a broad range in temperature. Furthermore, we show that thermal instabilities change the scale dependent anisotropy of the turbulent magnetic field, reducing the increase in eddy elongation at smaller scales. Finally, we show that most of the mass is contained in one contiguous cold structure surrounded by smaller clumps that follow a scale free mass distribution. These clumps tend to be highly elongated and exhibit a size versus internal velocity relation consistent with supersonic turbulence, and a relative clump distance-velocity scaling consistent with subsonic motion. We discuss the striking similarity of cold plasmoids to observed tiny scale atomic and ionized structures and HI fibers, and consider how the prevalence of plasmoids will modify the motion of charged particles thereby impacting cosmic ray transport and thermal conduction in the ISM and other similar systems.
△ Less
Submitted 29 April, 2024; v1 submitted 11 November, 2022;
originally announced November 2022.
-
Regulation of Star Formation by a Hot Circumgalactic Medium
Authors:
Christopher Carr,
Greg L. Bryan,
Drummond B. Fielding,
Viraj Pandya,
Rachel S. Somerville
Abstract:
Galactic outflows driven by supernovae (SNe) are thought to be a powerful regulator of a galaxy's star-forming efficiency. Mass, energy, and metal outflows ($η_M$, $η_E$, and $η_Z$, here normalized by the star formation rate, the SNe energy and metal production rates, respectively) shape galaxy properties by both ejecting gas and metals out of the galaxy and by heating the circumgalactic medium (C…
▽ More
Galactic outflows driven by supernovae (SNe) are thought to be a powerful regulator of a galaxy's star-forming efficiency. Mass, energy, and metal outflows ($η_M$, $η_E$, and $η_Z$, here normalized by the star formation rate, the SNe energy and metal production rates, respectively) shape galaxy properties by both ejecting gas and metals out of the galaxy and by heating the circumgalactic medium (CGM), preventing future accretion. Traditionally, models have assumed that galaxies self-regulate by ejecting a large fraction of the gas which enters the interstellar medium (ISM), even though such high mass-loadings are in growing tension with observations. To better understand how the relative importance of ejective (i.e. high mass-loading) vs preventative (i.e. high energy-loading) feedback affects the present-day properties of galaxies, we develop a simple gas-regulator model of galaxy evolution, where the stellar mass, ISM, and CGM are modeled as distinct reservoirs which exchange mass, metals, and energy at different rates within a growing halo. Focusing on the halo mass range from $10^{10}$ to $10^{12} M_{\odot}$, we demonstrate that, with reasonable parameter choices, we can reproduce the stellar-to-halo mass relation and the ISM-to-stellar mass relation with low mass-loaded ($η_M \sim 0.1-10$) but high energy-loaded ($η_E \sim 0.1-1$) winds, with self-regulation occurring primarily through heating and cooling of the CGM. We show that the model predictions are robust against changes to the mass-loading of outflows but are quite sensitive to our choice of the energy-loading, preferring $η_E \sim 1$ for the lowest mass halos and $\sim 0.1$ for Milky Way-like halos.
△ Less
Submitted 9 November, 2022;
originally announced November 2022.
-
The Anatomy of a Turbulent Radiative Mixing Layer: Insights from an Analytic Model with Turbulent Conduction and Viscosity
Authors:
Zirui Chen,
Drummond B. Fielding,
Greg L. Bryan
Abstract:
Turbulent Radiative Mixing Layers (TRMLs) form at the interface of cold, dense gas and hot, diffuse gas in motion with each other. TRMLs are ubiquitous in and around galaxies on a variety of scales, including galactic winds and the circumgalactic medium. They host the intermediate temperature gases that are efficient in radiative cooling, thus play a crucial role in controlling the cold gas supply…
▽ More
Turbulent Radiative Mixing Layers (TRMLs) form at the interface of cold, dense gas and hot, diffuse gas in motion with each other. TRMLs are ubiquitous in and around galaxies on a variety of scales, including galactic winds and the circumgalactic medium. They host the intermediate temperature gases that are efficient in radiative cooling, thus play a crucial role in controlling the cold gas supply, phase structure, and spectral features of galaxies. In this work, we develop an intuitive analytic 1.5 dimensional model for TRMLs that includes a simple parameterization of the effective turbulent conductivity and viscosity and a piece-wise power-law cooling curve. Our analytic model reproduces the mass flux, total cooling, and phase structure of 3D simulations of TRMLs at a fraction of the computational cost. It also reveals essential insights into the physics of TRMLs, particularly the importance of the viscous dissipation of relative kinetic energy in balancing radiative cooling as the shear Mach number approaches unity. This dissipation takes place both in the intermediate temperature phase, which reduces the enthalpy flux from the hot phase, and in the cold phase, which enhances radiative cooling. Additionally, our model provides a fast and easy way of computing the column density and surface brightness of TRMLs, which can be directly linked to observations.
△ Less
Submitted 29 June, 2023; v1 submitted 2 November, 2022;
originally announced November 2022.
-
Taming the TuRMoiL: The Temperature Dependence of Turbulence in Cloud-Wind Interactions
Authors:
Matthew W. Abruzzo,
Drummond B. Fielding,
Greg L. Bryan
Abstract:
Turbulent radiative mixing layers (TRMLs) play an important role in many astrophysical contexts where cool ($\lesssim 10^4$ K) clouds interact with hot flows (e.g., galactic winds, high velocity clouds, infalling satellites in halos and clusters). The fate of these clouds (as well as many of their observable properties) is dictated by the competition between turbulence and radiative cooling; howev…
▽ More
Turbulent radiative mixing layers (TRMLs) play an important role in many astrophysical contexts where cool ($\lesssim 10^4$ K) clouds interact with hot flows (e.g., galactic winds, high velocity clouds, infalling satellites in halos and clusters). The fate of these clouds (as well as many of their observable properties) is dictated by the competition between turbulence and radiative cooling; however, turbulence in these multiphase flows remains poorly understood. We have investigated the emergent turbulence arising in the interaction between clouds and supersonic winds in hydrodynamic ENZO-E simulations. In order to obtain robust results, we employed multiple metrics to characterize the turbulent velocity, $v_{\rm turb}$. We find four primary results, when cooling is sufficient for cloud survival. First, $v_{\rm turb}$ manifests clear temperature dependence. Initially, $v_{\rm turb}$ roughly matches the scaling of sound speed on temperature. In gas hotter than the temperature where cooling peaks, this dependence weakens with time until $v_{\rm turb}$ is constant. Second, the relative velocity between the cloud and wind initially drives rapid growth of $v_{\rm turb}$. As it drops (from entrainment), $v_{\rm turb}$ starts to decay before it stabilizes at roughly half its maximum. At late times cooling flows appear to support turbulence. Third, the magnitude of $v_{\rm turb}$ scales with the ratio between the hot phase sound crossing time and the minimum cooling time. Finally, we find tentative evidence for a length-scale associated with resolving turbulence. Under-resolving this scale may cause violent shattering and affect the cloud's large-scale morphological properties.
△ Less
Submitted 20 June, 2023; v1 submitted 27 October, 2022;
originally announced October 2022.
-
Hot-mode accretion and the physics of thin-disk galaxy formation
Authors:
Zachary Hafen,
Jonathan Stern,
James Bullock,
Alex B. Gurvich,
Sijie Yu,
Claude-Andre Faucher-Giguere,
Drummond B. Fielding,
Daniel Angles-Alcazar,
Eliot Quataert,
Andrew Wetzel,
Tjitske Starkenburg,
Michael Boylan-Kolchin,
Jorge Moreno,
Robert Feldmann,
Kareem El-Badry,
T. K. Chan,
Cameron Trapp,
Dusan Keres,
Philip F. Hopkins
Abstract:
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:…
▽ More
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.
△ Less
Submitted 6 June, 2022; v1 submitted 18 January, 2022;
originally announced January 2022.
-
Bursting Bubbles: Feedback from Clustered SNe and the Trade-off Between Turbulence and Outflows
Authors:
Matthew E. Orr,
Drummond B. Fielding,
Christopher C. Hayward,
Blakesley Burkhart
Abstract:
We present an analytic model for clustered supernovae (SNe) feedback in galaxy disks, incorporating the dynamical evolution of superbubbles formed from spatially overlapping SNe remnants. We propose two realistic outcomes for the evolution of superbubbles in galactic disks: (1) the expansion velocity of the shock front falls below the turbulent velocity dispersion of the ISM in the galaxy disk, wh…
▽ More
We present an analytic model for clustered supernovae (SNe) feedback in galaxy disks, incorporating the dynamical evolution of superbubbles formed from spatially overlapping SNe remnants. We propose two realistic outcomes for the evolution of superbubbles in galactic disks: (1) the expansion velocity of the shock front falls below the turbulent velocity dispersion of the ISM in the galaxy disk, whereupon the superbubble stalls and fragments, depositing its momentum entirely within the galaxy disk, or (2) the superbubble grows in size to reach the gas scale height, breaking out of the galaxy disk and driving galactic outflows/fountains. In either case, we find that superbubble breakup/breakout almost always occurs before the last Type-II SN ($\lesssim$40 Myr) in the recently formed star cluster, assuming a standard high-end IMF slope, and scalings between stellar lifetimes and masses. The threshold between these two cases implies a break in the effective strength of feedback in driving turbulence within galaxies, and a resulting change in the scalings of, e.g., star formation rates with gas surface density (the Kennicutt-Schmidt relation) and the star formation efficiency in galaxy disks.
△ Less
Submitted 2 May, 2022; v1 submitted 29 September, 2021;
originally announced September 2021.
-
Bursting Bubbles: Clustered Supernova Feedback in Local and High-redshift Galaxies
Authors:
Matthew E. Orr,
Drummond B. Fielding,
Christopher C. Hayward,
Blakesley Burkhart
Abstract:
We compare an analytic model for the evolution of supernova-driven superbubbles with observations of local and high-redshift galaxies, and the properties of intact HI shells in local star-forming galaxies. Our model correctly predicts the presence of superwinds in local star-forming galaxies (e.g., NGC 253) and the ubiquity of outflows near $z \sim 2$. We find that high-redshift galaxies may `capt…
▽ More
We compare an analytic model for the evolution of supernova-driven superbubbles with observations of local and high-redshift galaxies, and the properties of intact HI shells in local star-forming galaxies. Our model correctly predicts the presence of superwinds in local star-forming galaxies (e.g., NGC 253) and the ubiquity of outflows near $z \sim 2$. We find that high-redshift galaxies may `capture' 20-50\% of their feedback momentum in the dense ISM (with the remainder escaping into the nearby CGM), whereas local galaxies may contain $\lesssim$10\% of their feedback momentum from the central starburst. Using azimuthally averaged galaxy properties, we predict that most superbubbles stall and fragment \emph{within} the ISM, and that this occurs at, or near, the gas scale height. We find a consistent interpretation in the observed HI bubble radii and velocities, and predict that most will fragment within the ISM, and that those able to break-out originate from short dynamical time regions (where the dynamical time is shorter than feedback timescales). Additionally, we demonstrate that models with constant star cluster formation efficiency per Toomre mass are inconsistent with the occurrence of outflows from high-$z$ starbursts and local circumnuclear regions.
△ Less
Submitted 21 January, 2022; v1 submitted 29 September, 2021;
originally announced September 2021.
-
The Structure of Multiphase Galactic Winds
Authors:
Drummond B. Fielding,
Greg L. Bryan
Abstract:
We present a novel analytic framework to model the steady-state structure of multiphase galactic winds comprised of a hot, volume-filling component and a cold, clumpy component. We first derive general expressions for the structure of the hot phase for arbitrary mass, momentum, and energy sources terms. Next, informed by recent simulations, we parameterize the cloud-wind mass transfer rates, which…
▽ More
We present a novel analytic framework to model the steady-state structure of multiphase galactic winds comprised of a hot, volume-filling component and a cold, clumpy component. We first derive general expressions for the structure of the hot phase for arbitrary mass, momentum, and energy sources terms. Next, informed by recent simulations, we parameterize the cloud-wind mass transfer rates, which are set by the competition between turbulent mixing and radiative cooling. This enables us to cast the cloud-wind interaction as a source term for the hot phase and thereby simultaneously solve for the evolution of both phases fully accounting for their bidirectional influence. With this model, we explore the nature of galactic winds over a broad range of conditions. We find that: (i) with realistic parameter choices, we naturally produce a hot, low-density wind that transports energy while entraining a significant flux of cold clouds, (ii) mixing dominates the cold cloud acceleration and decelerates the hot wind, (iii) during mixing thermalization of relative kinetic energy provides significant heating, (iv) systems with low hot-phase mass loading factors and/or star formation rates can sustain higher initial cold phase mass loading factors, but the clouds are quickly shredded, and (v) systems with large hot-phase mass loading factors and/or star formation rates cannot sustain large initial cold-phase mass loading factors, but the clouds tend to grow with radius. Our results highlight the necessity of accounting for the multiphase structure of galactic winds, both physically and observationally, and have important implications for feedback in galactic systems.
△ Less
Submitted 11 August, 2021;
originally announced August 2021.
-
The Impact of Cosmic Rays on the Kinematics of the Circumgalactic Medium
Authors:
Iryna S. Butsky,
Jessica K. Werk,
Kirill Tchernyshyov,
Drummond B. Fielding,
Joseph Breneman,
Daniel Piacitelli,
Thomas R. Quinn,
N. Nicole Sanchez,
Akaxia Cruz,
Cameron B. Hummels,
Joseph N. Burchett,
Michael Tremmel
Abstract:
We use hydrodynamical simulations of two Milky Way-mass galaxies to demonstrate the impact of cosmic-ray pressure on the kinematics of cool and warm circumgalactic gas. Consistent with previous studies, we find that cosmic-ray pressure can dominate over thermal pressure in the inner 50 kpc of the circumgalactic medium (CGM), creating an overall cooler CGM than that of similar galaxy simulations ru…
▽ More
We use hydrodynamical simulations of two Milky Way-mass galaxies to demonstrate the impact of cosmic-ray pressure on the kinematics of cool and warm circumgalactic gas. Consistent with previous studies, we find that cosmic-ray pressure can dominate over thermal pressure in the inner 50 kpc of the circumgalactic medium (CGM), creating an overall cooler CGM than that of similar galaxy simulations run without cosmic rays. We generate synthetic sightlines of the simulated galaxies' CGM and use Voigt profile fitting methods to extract ion column densities, Doppler-b parameters, and velocity centroids of individual absorbers. We directly compare these synthetic spectral line fits with HST/COS CGM absorption-line data analyses, which tend to show that metallic species with a wide range of ionization potential energies are often kinematically aligned. Compared to the Milky-Way simulation run without cosmic rays, the presence of cosmic-ray pressure in the inner CGM creates narrower OVI absorption features and broader SiIII absorption features, a quality which is more consistent with observational data. Additionally, because the cool gas is buoyant due to nonthermal cosmic-ray pressure support, the velocity centroids of both cool and warm gas tend to align in the simulated Milky Way with feedback from cosmic rays. Our study demonstrates that detailed, direct comparisons between simulations and observations, focused on gas kinematics, have the potential to reveal the dominant physical mechanisms that shape the CGM.
△ Less
Submitted 17 August, 2022; v1 submitted 28 June, 2021;
originally announced June 2021.
-
Characterizing mass, momentum, energy and metal outflow rates of multi-phase galactic winds in the FIRE-2 cosmological simulations
Authors:
Viraj Pandya,
Drummond B. Fielding,
Daniel Anglés-Alcázar,
Rachel S. Somerville,
Greg L. Bryan,
Christopher C. Hayward,
Jonathan Stern,
Chang-Goo Kim,
Eliot Quataert,
John C. Forbes,
Claude-André Faucher-Giguère,
Robert Feldmann,
Zachary Hafen,
Philip F. Hopkins,
Dušan Kereš,
Norman Murray,
Andrew Wetzel
Abstract:
We characterize mass, momentum, energy and metal outflow rates of multi-phase galactic winds in a suite of FIRE-2 cosmological "zoom-in" simulations from the Feedback in Realistic Environments (FIRE) project. We analyze simulations of low-mass dwarfs, intermediate-mass dwarfs, Milky Way-mass halos, and high-redshift massive halos. Consistent with previous work, we find that dwarfs eject about 100…
▽ More
We characterize mass, momentum, energy and metal outflow rates of multi-phase galactic winds in a suite of FIRE-2 cosmological "zoom-in" simulations from the Feedback in Realistic Environments (FIRE) project. We analyze simulations of low-mass dwarfs, intermediate-mass dwarfs, Milky Way-mass halos, and high-redshift massive halos. Consistent with previous work, we find that dwarfs eject about 100 times more gas from their interstellar medium (ISM) than they form in stars, while this mass "loading factor" drops below one in massive galaxies. Most of the mass is carried by the hot phase ($>10^5$ K) in massive halos and the warm phase ($10^3-10^5$ K) in dwarfs; cold outflows ($<10^3$ K) are negligible except in high-redshift dwarfs. Energy, momentum and metal loading factors from the ISM are of order unity in dwarfs and significantly lower in more massive halos. Hot outflows have $2-5\times$ higher specific energy than needed to escape from the gravitational potential of dwarf halos; indeed, in dwarfs, the mass, momentum, and metal outflow rates increase with radius whereas energy is roughly conserved, indicating swept up halo gas. Burst-averaged mass loading factors tend to be larger during more powerful star formation episodes and when the inner halo is not virialized, but we see effectively no trend with the dense ISM gas fraction. We discuss how our results can guide future controlled numerical experiments that aim to elucidate the key parameters governing galactic winds and the resulting associated preventative feedback.
△ Less
Submitted 17 September, 2021; v1 submitted 11 March, 2021;
originally announced March 2021.
-
A simple model for mixing and cooling in cloud-wind interactions
Authors:
Matthew W. Abruzzo,
Greg L. Bryan,
Drummond B. Fielding
Abstract:
We introduce a simple entropy-based formalism to characterize the role of mixing in pressure-balanced multiphase clouds, and demonstrate example applications using Enzo-E (magneto)hydrodynamic simulations. Under this formalism, the high-dimensional description of the system's state at a given time is simplified to the joint distribution of mass over pressure ($P$) and entropy ($K=P/ρ^γ$). As a res…
▽ More
We introduce a simple entropy-based formalism to characterize the role of mixing in pressure-balanced multiphase clouds, and demonstrate example applications using Enzo-E (magneto)hydrodynamic simulations. Under this formalism, the high-dimensional description of the system's state at a given time is simplified to the joint distribution of mass over pressure ($P$) and entropy ($K=P/ρ^γ$). As a result, this approach provides a way for (empirically and analytically) quantifying the impact of different initial conditions and sets of physics on the system evolution. We find that mixing predominantly alters the distribution along the $K$ direction and illustrate how the formalism can be used to model mixing and cooling for fluid elements originating in the cloud. We further confirm and generalize a previously suggested criterion for cloud growth in the presence of radiative cooling, and demonstrate that the shape of the cooling curve, particularly at the low temperature end, can play an important role in controlling condensation. Moreover, we discuss the capacity of our approach to generalize such a criterion to apply to additional sets of physics, and to build intuition for the impact of subtle higher order effects not directly addressed by the criterion.
△ Less
Submitted 25 January, 2021;
originally announced January 2021.
-
A Framework for Multiphase Galactic Wind Launching using TIGRESS
Authors:
Chang-Goo Kim,
Eve C. Ostriker,
Drummond B. Fielding,
Matthew C. Smith,
Greg L. Bryan,
Rachel S. Somerville,
John C. Forbes,
Shy Genel,
Lars Hernquist
Abstract:
Galactic outflows have density, temperature, and velocity variations at least as large as that of the multiphase, turbulent interstellar medium (ISM) from which they originate. We have conducted a suite of parsec-resolution numerical simulations using the TIGRESS framework, in which outflows emerge as a consequence of interaction between supernovae (SNe) and the star-forming ISM. The outflowing ga…
▽ More
Galactic outflows have density, temperature, and velocity variations at least as large as that of the multiphase, turbulent interstellar medium (ISM) from which they originate. We have conducted a suite of parsec-resolution numerical simulations using the TIGRESS framework, in which outflows emerge as a consequence of interaction between supernovae (SNe) and the star-forming ISM. The outflowing gas is characterized by two distinct thermal phases, cool (T<10^4 K) and hot (T>10^6 K), with most mass carried by the cool phase and most energy and newly-injected metals carried by the hot phase. Both components have a broad distribution of outflow velocity, and especially for cool gas this implies a varying fraction of escaping material depending on the halo potential. Informed by the TIGRESS results, we develop straightforward analytic formulae for the joint probability density functions (PDFs) of mass, momentum, energy, and metal loading as distributions in outflow velocity and sound speed. The model PDFs have only two parameters, SFR surface density Σ_SFR and the metallicity of the ISM, and fully capture the behavior of the original TIGRESS simulation PDFs over Σ_SFR~(10^{-4},1)M_sun/kpc^2/yr. Employing PDFs from resolved simulations will enable galaxy formation subgrid model implementations with wind velocity and temperature (as well as total loading factors) that are based on theoretical predictions rather than empirical tuning. This is a critical step to incorporate advances from TIGRESS and other high-resolution simulations in future cosmological hydrodynamics and semi-analytic galaxy formation models. We release a python package to prototype our model and to ease its implementation.
△ Less
Submitted 18 October, 2020;
originally announced October 2020.
-
The Impact of Cosmic Rays on Thermal Instability in the Circumgalactic Medium
Authors:
Iryna S. Butsky,
Drummond B. Fielding,
Christopher C. Hayward,
Cameron B. Hummels,
Thomas R. Quinn,
Jessica K. Werk
Abstract:
Large reservoirs of cold (~ 10^4 K) gas exist out to and beyond the virial radius in the circumgalactic medium (CGM) of all types of galaxies. Photoionization modeling suggests that cold CGM gas has significantly lower densities than expected by theoretical predictions based on thermal pressure equilibrium with hot CGM gas. In this work, we investigate the impact of cosmic ray physics on the forma…
▽ More
Large reservoirs of cold (~ 10^4 K) gas exist out to and beyond the virial radius in the circumgalactic medium (CGM) of all types of galaxies. Photoionization modeling suggests that cold CGM gas has significantly lower densities than expected by theoretical predictions based on thermal pressure equilibrium with hot CGM gas. In this work, we investigate the impact of cosmic ray physics on the formation of cold gas via thermal instability. We use idealized three-dimensional magnetohydrodynamic simulations to follow the evolution of thermally unstable gas in a gravitationally stratified medium. We find that cosmic ray pressure lowers the density and increases the size of cold gas clouds formed through thermal instability. We develop a simple model for how the cold cloud sizes and the relative densities of cold and hot gas depend on cosmic ray pressure. Cosmic ray pressure can help counteract gravity to keep cold gas in the CGM for longer, thereby increasing the predicted cold mass fraction and decreasing the predicted cold gas inflow rates. Efficient cosmic ray transport, by streaming or diffusion, redistributes cosmic ray pressure from the cold gas to the background medium, resulting in cold gas properties that are in-between those predicted by simulations with inefficient transport and simulations without cosmic rays. We show that cosmic rays can significantly reduce galactic accretion rates and resolve the tension between theoretical models and observational constraints on the properties of cold CGM gas.
△ Less
Submitted 23 September, 2020; v1 submitted 11 August, 2020;
originally announced August 2020.
-
First results from SMAUG: The need for preventative stellar feedback and improved baryon cycling in semi-analytic models of galaxy formation
Authors:
Viraj Pandya,
Rachel S. Somerville,
Daniel Anglés-Alcázar,
Christopher C. Hayward,
Greg L. Bryan,
Drummond B. Fielding,
John C. Forbes,
Blakesley Burkhart,
Shy Genel,
Lars Hernquist,
Chang-Goo Kim,
Stephanie Tonnesen,
Tjitske Starkenburg
Abstract:
Semi-analytic models (SAMs) are a promising means of tracking the physical processes associated with galaxy formation, but many of their approximations have not been rigorously tested. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, we compare predictions from the FIRE-2 hydrodynamical "zoom-in" simulations to those from the Santa Cruz SAM run on the same…
▽ More
Semi-analytic models (SAMs) are a promising means of tracking the physical processes associated with galaxy formation, but many of their approximations have not been rigorously tested. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, we compare predictions from the FIRE-2 hydrodynamical "zoom-in" simulations to those from the Santa Cruz SAM run on the same halo merger trees, with an emphasis on the global mass flow cycle. Our study includes 13 halos spanning low-mass dwarfs (M_vir~10^10 M_sun at z=0), intermediate-mass dwarfs (M_vir~10^11 M_sun) and Milky Way-mass galaxies (M_vir~10^12 M_sun). The SAM and FIRE-2 predictions agree relatively well with each other in terms of stellar and interstellar mass, but differ dramatically on circumgalactic mass (the SAM is lower than FIRE-2 by ~3 orders of magnitude for dwarfs). Strikingly, the SAM predicts higher gas accretion rates for dwarfs compared to FIRE-2 by factors of ~10-100, and this is compensated for with higher mass outflow rates in the SAM. We argue that the most severe model discrepancies are caused by the lack of preventative stellar feedback and the assumptions for halo gas cooling and recycling in the SAM. As a first step towards resolving these model tensions, we present a simple yet promising new preventative stellar feedback model in which the energy carried by supernova-driven winds is allowed to heat some fraction of gas outside of halos to at least the virial temperature such that accretion is suppressed.
△ Less
Submitted 30 October, 2020; v1 submitted 29 June, 2020;
originally announced June 2020.
-
First results from SMAUG: Uncovering the Origin of the Multiphase Circumgalactic Medium with a Comparative Analysis of Idealized and Cosmological Simulations
Authors:
Drummond B. Fielding,
Stephanie Tonnesen,
Daniel DeFelippis,
Miao Li,
Kung-Yi Su,
Greg L. Bryan,
Chang-Goo Kim,
John C. Forbes,
Rachel S. Somerville,
Nicholas Battaglia,
Evan E. Schneider,
Yuan Li,
Ena Choi,
Christopher C. Hayward,
Lars Hernquist
Abstract:
We examine the properties of the circumgalactic medium (CGM) at low redshift in a range of simulated Milky Way mass halos. The sample is comprised of seven idealized simulations, an adaptive mesh refinement cosmological zoom-in simulation, and two groups of 50 halos with star forming or quiescent galaxies taken from the IllustrisTNG100 simulation. The simulations have very different setups, resolu…
▽ More
We examine the properties of the circumgalactic medium (CGM) at low redshift in a range of simulated Milky Way mass halos. The sample is comprised of seven idealized simulations, an adaptive mesh refinement cosmological zoom-in simulation, and two groups of 50 halos with star forming or quiescent galaxies taken from the IllustrisTNG100 simulation. The simulations have very different setups, resolution, and feedback models, but are analyzed in a uniform manner. By comparing median radial profiles and mass distributions of CGM properties, we isolate key similarities and differences. In doing so, we advance the efforts of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project that aims to understand the inherently multiscale galaxy formation process. In the cosmological simulations, the CGM exhibits nearly flat temperature distributions, and broad pressure and radial velocity distributions. In the idealized simulations, similar distributions are found in the inner CGM ($\lesssim 0.5 \, r_{\rm 200c}$) when strong galactic feedback models are employed, but the outer CGM ($\gtrsim 0.5 \, r_{\rm 200c}$) has a much less prominent cold phase, and narrower pressure and velocity distributions even in models with strong feedback. This comparative analysis demonstrates the dominant role feedback plays in shaping the inner CGM and the increased importance of cosmological effects, such as nonspherical accretion and satellite galaxies, in the outer CGM. Furthermore, our findings highlight that while cosmological simulations are required to capture the multiphase structure of the CGM at large radii, idealized simulations provide a robust framework to study how galactic feedback interacts with the inner CGM and thereby provide a reliable avenue to constrain feedback prescriptions.
△ Less
Submitted 29 June, 2020;
originally announced June 2020.
-
First results from SMAUG: Characterization of Multiphase Galactic Outflows from a Suite of Local Star-Forming Galactic Disk Simulations
Authors:
Chang-Goo Kim,
Eve C. Ostriker,
Rachel S. Somerville,
Greg L. Bryan,
Drummond B. Fielding,
John C. Forbes,
Christopher C. Hayward,
Lars Hernquist,
Viraj Pandya
Abstract:
Large scale outflows in star-forming galaxies are observed to be ubiquitous, and are a key aspect of theoretical modeling of galactic evolution in a cosmological context, the focus of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project. Gas blown out from galactic disks, similar to gas within galaxies, consists of multiple phases with large contrasts of density, temperatu…
▽ More
Large scale outflows in star-forming galaxies are observed to be ubiquitous, and are a key aspect of theoretical modeling of galactic evolution in a cosmological context, the focus of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project. Gas blown out from galactic disks, similar to gas within galaxies, consists of multiple phases with large contrasts of density, temperature, and other properties. To study multiphase outflows as emergent phenomena, we run a suite of ~pc-resolution local galactic disk simulations using the TIGRESS framework. Explicit modeling of the interstellar medium (ISM), including star formation and self-consistent radiative heating plus supernova feedback, regulates ISM properties and drives the outflow. We investigate the scaling of outflow mass, momentum, energy, and metal loading factors with galactic disk properties, including star formation rate (SFR) surface density (Σ_SFR~10^{-4}-1 M_sun/kpc^2/yr), gas surface density (~1-100 M_sun/pc^2), and total midplane pressure (or weight) (~10^3-10^6 k_B cm^{-3} K). The main components of outflowing gas are mass-delivering cool gas (T~10^4 K) and energy/metal-delivering hot gas (T~10^6 K). Cool mass outflow rates measured at outflow launch points (one or two scale heights) are 1-100 times the SFR (decreasing with Σ_SFR), although in massive galaxies most mass falls back due to insufficient outflow velocity. The hot galactic outflow carries mass comparable to 10% of the SFR, together with 10-20% of the energy and 30-60% of the metal mass injected by SN feedback. The characteristic outflow velocities of both phases scale very weakly with SFR, as v_out \propto Σ_SFR^{0.1~0.2}, consistent with observations. Importantly, our analysis demonstrates that in any physically-motivated cosmological wind model, it is crucial to include at least two distinct thermal wind components.
△ Less
Submitted 27 July, 2020; v1 submitted 29 June, 2020;
originally announced June 2020.
-
Multiphase Gas and the Fractal Nature of Radiative Turbulent Mixing Layers
Authors:
Drummond B. Fielding,
Eve C. Ostriker,
Greg L. Bryan,
Adam S. Jermyn
Abstract:
A common situation in galactic and intergalactic gas involves cold dense gas in motion relative to hot diffuse gas. Kelvin-Helmholtz instability creates a turbulent mixing layer and populates the intermediate-temperature phase, which often cools rapidly. The energy lost to cooling is balanced by the advection of hot high enthalpy gas into the mixing layer, resulting in growth and acceleration of t…
▽ More
A common situation in galactic and intergalactic gas involves cold dense gas in motion relative to hot diffuse gas. Kelvin-Helmholtz instability creates a turbulent mixing layer and populates the intermediate-temperature phase, which often cools rapidly. The energy lost to cooling is balanced by the advection of hot high enthalpy gas into the mixing layer, resulting in growth and acceleration of the cold phase. This process may play a major role in determining the interstellar medium and circumgalactic medium phase structure, and accelerating cold gas in galactic winds and cosmic filaments. Cooling in these mixing layers occurs in a thin corrugated sheet, which we argue has an area with fractal dimension $D=5/2$ and a thickness that adjusts to match the hot phase mixing time to the cooling time. These cooling sheet properties form the basis of a new model for how the cooling rate and hot gas inflow velocity depend on the size $L$, cooling time $t_{\rm cool}$, relative velocity $v_{\rm rel}$, and density contrast $ρ_{\rm cold}/ρ_{\rm hot}$ of the system. Entrainment is expected to be enhanced in environments with short $t_{\rm cool}$, large $v_{\rm rel}$, and large $ρ_{\rm cold}/ρ_{\rm hot}$. Using a large suite of three dimensional hydrodynamic simulations, we demonstrate that this fractal cooling layer model accurately captures the energetics and evolution of turbulent interfaces and can therefore be used as a foundation for understanding multiphase mixing with strong radiative cooling.
△ Less
Submitted 18 March, 2020;
originally announced March 2020.
-
The Turbulent Origin of Outflow and Spin Misalignment in Multiple Star Systems
Authors:
Stella S. R. Offner,
Michael M. Dunham,
Katherine I. Lee,
Hector G. Arce,
Drummond B. Fielding
Abstract:
The protostellar outflows of wide-separation forming binaries frequently appear misaligned. We use magneto-hydrodynamic simulations to investigate the alignment of protostellar spin and molecular outflows for forming binary pairs. We show that the protostellar pairs, which form from turbulent fragmentation within a single parent core, have randomly oriented angular momenta. Although the pairs migr…
▽ More
The protostellar outflows of wide-separation forming binaries frequently appear misaligned. We use magneto-hydrodynamic simulations to investigate the alignment of protostellar spin and molecular outflows for forming binary pairs. We show that the protostellar pairs, which form from turbulent fragmentation within a single parent core, have randomly oriented angular momenta. Although the pairs migrate to closer separations, their spins remain partially misaligned. We produce $^{12}$CO(2-1) synthetic observations of the simulations and characterize the outflow orientation in the emission maps. The CO-identified outflows exhibit a similar random distribution and are also statistically consistent with the observed distribution of molecular outflows. We conclude that observed misalignment provides a clear signature of binary formation via turbulent fragmentation. The persistence of misaligned outflows and stellar spins following dynamical evolution may provide a signature of binary origins for more evolved multiple star systems.
△ Less
Submitted 27 June, 2016;
originally announced June 2016.
-
The Turbulent Origin of Spin-Orbit Misalignment in Planetary Systems
Authors:
Drummond B. Fielding,
Christopher F. McKee,
Aristotle Socrates,
Andrew J. Cunningham,
Richard I. Klein
Abstract:
The turbulent environment from which stars form may lead to misalignment between the stellar spin and the remnant protoplanetary disk. By using hydrodynamic and magnetohydrodynamic simulations, we demonstrate that a wide range of stellar obliquities may be produced as a by-product of forming a star within a turbulent environment. We present a simple semi-analytic model that reveals this connection…
▽ More
The turbulent environment from which stars form may lead to misalignment between the stellar spin and the remnant protoplanetary disk. By using hydrodynamic and magnetohydrodynamic simulations, we demonstrate that a wide range of stellar obliquities may be produced as a by-product of forming a star within a turbulent environment. We present a simple semi-analytic model that reveals this connection between the turbulent motions and the orientation of a star and its disk. Our results are consistent with the observed obliquity distribution of hot Jupiters. Migration of misaligned hot Jupiters may, therefore, be due to tidal dissipation in the disk, rather than tidal dissipation of the star-planet interaction.
△ Less
Submitted 8 June, 2015; v1 submitted 17 September, 2014;
originally announced September 2014.