Showing 1–50 of 87 results for author: Fielding, D

Illuminating M82: Simulating X-ray Emission from Galactic Winds in a Starburst Galaxy

Authors: Allison Lin, Lachlan Lancaster, Viraj Pandya, Drummond B. Fielding, Greg L. Bryan, John A. ZuHone, Sebastian Lopez, Laura A. Lopez, Christopher Carr

Abstract: We generate mock X-ray observations from a suite of idealized high-resolution ($\sim 4$ pc), tall-box ($\sim 2 \times 2 \times 8$ kpc$^3$) simulations of star formation driven galactic winds in an M82-like system, varying the spatial resolution as well as the strength and distribution of supernova (SN) energy injection. We compare our mock X-ray observations with deep Chandra observations of the h… ▽ More We generate mock X-ray observations from a suite of idealized high-resolution ($\sim 4$ pc), tall-box ($\sim 2 \times 2 \times 8$ kpc$^3$) simulations of star formation driven galactic winds in an M82-like system, varying the spatial resolution as well as the strength and distribution of supernova (SN) energy injection. We compare our mock X-ray observations with deep Chandra observations of the hot plasma around M82. While the simulated total X-ray luminosity, $L_X$, increases with resolution and when SNe feedback is spatially distributed, even in the best case scenario, our simulated $L_X$ is a factor of $\sim 50-100$ lower than observed and the surface brightness profiles of X-ray emission, $S_X(z)$, fall off too quickly with distance from the galaxy. Past results were able to reproduce these observables and we discuss potential simulation differences that could explain this discrepancy. We make the first comparison of the X-ray spectrum of our simulations to observations and find that our simulated spectrum is too soft, with a deficit of hard X-ray photons at $\gtrsim 1$ keV. We discuss how physical processes missing from our simulations and prior work (e.g., thermal conduction and cosmic rays) could help resolve this discrepancy. △ Less

Submitted 13 July, 2026; originally announced July 2026.

Comments: 16 pages, 7 figures, 2 tables. Submitted to The Astrophysical Journal

Phase-dependent magnetic coherence in the turbulent interstellar medium

Authors: Iryna S. Butsky, Caleb Redshaw, Minjie Lei, Susan E. Clark, Drummond B. Fielding

Abstract: Magnetic fields permeate the multiphase interstellar medium (ISM), yet their phase-dependent structure remains poorly constrained by observations. Dust polarization and \ion{H}{1} emission together offer complementary probes of the plane-of-sky magnetic field and cold neutral medium (CNM) gas structure, respectively. Recent observational work has shown that in the diffuse ISM, the dust polarizatio… ▽ More Magnetic fields permeate the multiphase interstellar medium (ISM), yet their phase-dependent structure remains poorly constrained by observations. Dust polarization and \ion{H}{1} emission together offer complementary probes of the plane-of-sky magnetic field and cold neutral medium (CNM) gas structure, respectively. Recent observational work has shown that in the diffuse ISM, the dust polarization fraction correlates positively with the CNM mass fraction ($f_{\rm CNM}$) but not with total \ion{H}{1} column density, suggesting a phase-dependent magnetic field geometry. Here, we use extremely high-resolution ($2048^3$) simulations of the turbulent, magnetized, multiphase ISM to investigate the physical origin of this trend. By constructing synthetic \ion{H}{1} and dust polarization maps, we directly compare our simulations to the observational results of \citet{Lei:2024}. We recover a positive $f_{\rm CNM}$-polarization correlation most clearly for sightlines intersecting fewer than $\sim$20 discrete CNM clouds, while the trend becomes weak or intermittent for larger cloud counts, consistent with the expectation that high-Galactic-latitude sightlines contain relatively few independent cold structures. We show that this correlation reflects genuine phase-dependent magnetic structure: CNM clouds tend to be elongated along the local magnetic field and, when normalized by column density, exhibit lower magnetic disorder than the warm neutral medium (WNM). We further demonstrate that apparent discrepancies between simulation- and observation-based measures of magnetic disorder arise from whether disorder is quantified per unit path length or per unit mass. Our results support a picture in which CNM structures host relatively ordered magnetic fields, producing higher polarization fractions along CNM-dominated sightlines in the diffuse ISM. △ Less

Submitted 9 June, 2026; originally announced June 2026.

Comments: 10 pages, 6 figures. Submitted to ApJ

Clumps in a Cocoon: Geometry and Mixing Set the Universal X-ray to H$α$ Surface Brightness Ratio

Authors: Zirui Chen, S. Peng Oh, Drummond B. Fielding, Lachlan Lancaster, Yuan Li, Brent Tan

Abstract: Recent observations reveal a universal X-ray to H$α$ surface-brightness ratio, ${\rm SB}_{\rm X}/{\rm SB}_{\rm Hα}\sim 3$, in galactic winds, ram-pressure stripped tails, and cluster filaments. This is surprising because H$α$ traces cold ($\sim 10^4$ K) gas while X-rays trace much hotter ($\sim 10^{6}$--$10^{7}$ K) gas. Plane-parallel mixing-layer models do not recover this ratio, and can be off b… ▽ More Recent observations reveal a universal X-ray to H$α$ surface-brightness ratio, ${\rm SB}_{\rm X}/{\rm SB}_{\rm Hα}\sim 3$, in galactic winds, ram-pressure stripped tails, and cluster filaments. This is surprising because H$α$ traces cold ($\sim 10^4$ K) gas while X-rays trace much hotter ($\sim 10^{6}$--$10^{7}$ K) gas. Plane-parallel mixing-layer models do not recover this ratio, and can be off by orders of magnitude. Motivated by recent work showing that geometry controls the temperature PDF of multiphase gas (Chen & Oh 2026), we run 3D wind-tunnel simulations in the high density contrast ($χ\sim 10^3$) regime. In this limit, the cold phase shatters into many small H$α$-emitting clumps, while X-ray-emitting gas forms a volume-filling cocoon around them. After smoothing on the tail-width scale, the measured surface-brightness ratio converges to the observed value, which can be understood theoretically. The H$α$ luminosity fraction is set by atomic physics, whereas the X-ray luminosity fraction is set by the residence time of gas in the X-ray-emitting band. This residence time is much shorter than the cooling time at X-ray temperatures, but scales roughly inversely with pressure, suggesting that it is tied to the cooling time at a lower-temperature outlet of the mixing cascade. This framework naturally explains why the observed ratio is order unity, and robust to changes in gas pressure. △ Less

Submitted 5 June, 2026; originally announced June 2026.

Comments: 6 pages, 3 figures. Submitted to MNRAS Letters

The Origin of Da Scaling: Suppressed Cooling in Fast-Cooling Mixing Layers

Authors: Lachlan Lancaster, Drummond B. Fielding, Rajsekhar Mohapatra, Greg L. Bryan

Abstract: In numerical experiments simulating Turbulent Radiative Mixing Layers (TRMLs) it is observed that as the cooling time in the mixed gas, $t_{\rm cool}$, becomes very short compared to the dynamical time of the turbulence, $t_{\rm eddy}/t_{\rm cool} \gg 1$, there is a change in the scaling behavior of the total energy radiated in the TRML as a function of this ratio, also known as the Damköhler numb… ▽ More In numerical experiments simulating Turbulent Radiative Mixing Layers (TRMLs) it is observed that as the cooling time in the mixed gas, $t_{\rm cool}$, becomes very short compared to the dynamical time of the turbulence, $t_{\rm eddy}/t_{\rm cool} \gg 1$, there is a change in the scaling behavior of the total energy radiated in the TRML as a function of this ratio, also known as the Damköhler number, ${\rm Da} \equiv t_{\rm eddy}/t_{\rm cool}$, from $\dot{E}_{\rm cool} \propto {\rm Da}^{1/2}$ to $\dot{E}_{\rm cool} \propto {\rm Da}^{1/4}$. The latter, so-called "fast-cooling," regime is of particular interest as many astrophysical mixing layers lie in this regime. We demonstrate that the origin of this change is the suppression of turbulent folding of the surface by the ram-pressure of the inflowing gas, which becomes much greater than the turbulent pressure in this regime. We present an argument that reproduces the $\dot{E}_{\rm cool} \propto {\rm Da}^{1/4}$ behavior by appealing to the suppression of the fractal structure of the interface by the ram-pressure of the inflowing gas. △ Less

Submitted 2 June, 2026; originally announced June 2026.

Comments: 11 pages, 5 figures, submitted to ApJ

Ceci n'est pas une Couche de Mélange: The Meaning of Resolved Turbulent Radiative Mixing

Authors: Lachlan Lancaster, Rajsekhar Mohapatra, Drummond B. Fielding, Greg L. Bryan

Abstract: Turbulent Radiative Mixing Layers (TRMLs) are of fundamental importance to the transport of energy and momentum in multi-phase, astrophysical fluids. We use measurements of the "micro" and "macro" properties of these layers in high-resolution \texttt{AthenaK} simulations to investigate when their properties can be considered \textit{well}-resolved. In particular, we demonstrate that the previously… ▽ More Turbulent Radiative Mixing Layers (TRMLs) are of fundamental importance to the transport of energy and momentum in multi-phase, astrophysical fluids. We use measurements of the "micro" and "macro" properties of these layers in high-resolution \texttt{AthenaK} simulations to investigate when their properties can be considered \textit{well}-resolved. In particular, we demonstrate that the previously noticed resolution independence of total cooling, $\dot{E}_{\rm cool}$, in these simulations is due to a remarkable, and perhaps fortuitous, cancellation of the countervailing effects of numerical dissipation and numerical viscosity. This calls into question the degree to which we can trust the results of these experiments, as there is no physical picture that explains this cancellation. We also demonstrate that in order to correctly resolve the phase structure in these layers, important for accurate predictions of their observable properties, one must resolve the scale on which turbulent diffusion acts on time-scales comparable to the cooling time. This "turbulent Field length", $λ_{\rm F,turb}$, is where the eddy turnover time is equal to the cooling time ($t_{\rm eddy}(λ_{\rm F,turb}) = t_{\rm cool}$). We demonstrate that resolving this scale results in converged phase-structure and spatially resolved transitions in the gas phases. △ Less

Submitted 2 June, 2026; originally announced June 2026.

Comments: 19 pages, 10 figure, Accepted in ApJ

ArkenstoneBH. A model for high-specific energy black hole feedback in cosmological simulations

Authors: James M. Sullivan, Greg L. Bryan, Matthew C. Smith, Jake S. Bennett, Drummond B. Fielding, Bryan A. Terrazas, Sophie Koudmani, Rachel S. Somerville, Michaela Hirschmann

Abstract: AGN feedback is a key piece of galaxy evolution but is difficult to model due to its high specific energies, multiphase nature, and limited simulation resolutions. Arkenstone is a subgrid framework for representing multiphase flows in coarse resolution simulations that has been used to model stellar feedback driven galactic winds. It ensures the correct treatment of high specific energy feedback t… ▽ More AGN feedback is a key piece of galaxy evolution but is difficult to model due to its high specific energies, multiphase nature, and limited simulation resolutions. Arkenstone is a subgrid framework for representing multiphase flows in coarse resolution simulations that has been used to model stellar feedback driven galactic winds. It ensures the correct treatment of high specific energy feedback that would otherwise be challenging to model accurately in Lagrangian simulations. We introduce the new Arkenstone BH model, which extends the Arkenstone framework to model black hole feedback. We focus on describing the first piece of this framework, which follows the hot, high specific energy phase of these outflows. The second piece, which treats their multiphase structure with a scheme for modeling unresolved cold clouds, will be implemented and described in a later paper. We present Arkenstone BH in simulations of an isolated galaxy to demonstrate the framework and its ability to capture high specific energy feedback that interacts only weakly with cold, dense gas. We show how these energetic outflows suppress star formation in our isolated galaxy by counteracting the inflow of gas from the circumgalactic medium into the interstellar medium. This work is part of the "Learning the Universe" collaboration, which aims to understand the Universe's underlying physics and initial conditions. △ Less

Submitted 4 May, 2026; originally announced May 2026.

Comments: 23 pages, 21 figures

Resolving the Unresolved Galactic Winds in Multi-phase Models. I. Methodology and Application

Authors: Xinfeng Xu, Drummond Fielding, Timothy Heckman, Greg L. Bryan, Alaina Henry, Karla Z. Arellano-Cordova, Cody Carr, John Chisholm, Claude-Andre Faucher-Giguere, Matthew Hayes, Mason Huberty, Michael Jennings, Crystal L. Martin, Claudia Scarlata, Allison L. Strom

Abstract: Galactic winds shape galaxy evolution; however, the outflowing gas is complex: it consists of multiple ionization phases, and its properties vary spatially. Therefore, methods that combine high-fidelity observations with state-of-the-art galactic-wind models are limited. Here we investigate methods for fitting the column density profiles derived from high-quality outflow observations with the mult… ▽ More Galactic winds shape galaxy evolution; however, the outflowing gas is complex: it consists of multiple ionization phases, and its properties vary spatially. Therefore, methods that combine high-fidelity observations with state-of-the-art galactic-wind models are limited. Here we investigate methods for fitting the column density profiles derived from high-quality outflow observations with the multiphase, multiscale wind model from Fielding & Bryan 2022. We identify three key outflow parameters: the initial hot-phase mass-loading factor ($η_\text{ M,hot,0}$), the initial cool-phase mass-loading factor ($η_\text{ M,cool,0}$), and the initial cool-cloud mass. We obtain good fits for most galaxies, with tight constraints on $η_\text{ M,cool,0}$ and moderate constraints on the other two parameters. We find the inferred $η_\text{ M,cool,0}$ and $η_\text{ M,hot,0}$ are mostly of order unity, with significant scatter. The constraints on $η_\text{ M,hot,0}$ suggest that the interaction between the cool and hot phases allows us to indirectly constrain the properties of the hot wind from cool-outflow observations. The model also predicts various radial trends. First, for all galaxies, the cool-phase outflow velocity increases between $1-2$ times of the half-light radius, then reaches a plateau. Second, most galaxies exhibit increasing $η_\text{ M,cool}$ and decreasing $η_\text{ M,hot}$ with radius, with a few showing the reverse trends. These results are effective, model-conditional constraints, and are consistent with other recent multiphase simulations and observations. This highlights that the velocity-radius mapping encoded in UV absorption profiles enables recovery of outflow spatial structures from spatially integrated spectra. Our method paves the way for future broad parameter studies and guides updates of outflow simulations in future work. △ Less

Submitted 1 May, 2026; originally announced May 2026.

Comments: 18 pages, 8 figures, 1 table, accepted by the Astrophysical Journal

Introducing sapphire: Towards Hybrid Physics-Informed, Data-Driven Modeling of Galaxy Formation

Authors: Viraj Pandya, Greg L. Bryan, T. Lucas Makinen, Austen Gabrielpillai, Christopher Carr, Drummond B. Fielding, Lars Hernquist, Matthew Ho, Kartheik Iyer, Christian Kragh Jespersen, Sophie Koudmani, Marta Laska, Pablo Lemos, Christopher C. Lovell, Lucia A. Perez, William F. Robinson Jr., Rachel S. Somerville, Tjitske K. Starkenburg, Richard Stiskalek, Bryan Terrazas, G. Mark Voit

Abstract: Semi-analytic models (SAMs) have been treating galaxy populations as dynamical systems for $\gtrsim50$ years, but their evolution equations remain poorly constrained. We introduce sapphire, a modular, automatically differentiable, GPU-accelerated SAM written in JAX. For the first time, we compute exact Jacobian and Hessian matrices of a galaxy formation SAM, using the Pandya et al. (2023) nonlinea… ▽ More Semi-analytic models (SAMs) have been treating galaxy populations as dynamical systems for $\gtrsim50$ years, but their evolution equations remain poorly constrained. We introduce sapphire, a modular, automatically differentiable, GPU-accelerated SAM written in JAX. For the first time, we compute exact Jacobian and Hessian matrices of a galaxy formation SAM, using the Pandya et al. (2023) nonlinear differential equation system as an example. These allow efficient, interpretable local and global sensitivity analyses, which reveal that supernova energy loading is the key astrophysical parameter. We use gradient descent and Hamiltonian Monte Carlo (HMC) to perform comprehensive mock parameter recovery tests. These indicate that the $z=0$ stellar-to-halo-mass relation alone does not contain enough information to infer many astrophysical parameters. Using observations of star-forming galaxies from the MaNGA survey and the Behroozi et al. (2019) empirical model as one baseline, we derive multiple posteriors assuming different combinations of data, including $z=0$ interstellar medium gas fractions and metallicities. The inferred physical parameters suggest that galaxies self-regulate their star formation primarily through preventative rather than ejective feedback, though this remains uncertain due to the lack of satellite galaxies, black holes and multi-phase galactic atmosphere physics. Both Fisher and HMC forecasts demonstrate the potential of sapphire to enable precision inference for galaxy formation and cosmology in a hybrid physics-informed, data-driven way, but more work is needed to expand its library of models and methods. We make sapphire publicly available at https://github.com/virajpandya/sapphire. △ Less

Submitted 12 August, 2026; v1 submitted 7 April, 2026; originally announced April 2026.

Comments: Resubmitted to ApJ after minor text clarifications and streamlining, comments welcome, code to reproduce analysis and figures at https://github.com/virajpandya/sapphire

Zooming in on the circumgalactic medium with GIBLE: Cloud-scale simulations with cosmological initial conditions

Authors: Rahul Ramesh, Dylan Nelson, Drummond Fielding, Marcus Brüggen

Abstract: We conduct simulations of $\sim$kpc-scale cool clouds in the circumgalactic medium (CGM), using initial conditions sampled from a highly resolved cosmological magneto-hydrodynamical zoom-in of a Milky Way-like galaxy. We select ten distinct cold clouds with masses of $m_{\rm{cloud}}$ $\sim$ $10^{4.5-5}$ M$_\odot$, originally resolved at a mass resolution of $m_{\rm{gas}}$ $\sim$ $200$ M$_\odot$. T… ▽ More We conduct simulations of $\sim$kpc-scale cool clouds in the circumgalactic medium (CGM), using initial conditions sampled from a highly resolved cosmological magneto-hydrodynamical zoom-in of a Milky Way-like galaxy. We select ten distinct cold clouds with masses of $m_{\rm{cloud}}$ $\sim$ $10^{4.5-5}$ M$_\odot$, originally resolved at a mass resolution of $m_{\rm{gas}}$ $\sim$ $200$ M$_\odot$. To further resolve small-scale features and physics, we implement a targeted refinement scheme within spherical regions co-moving with each cloud, thereby boosting the local mass resolution by a factor of 1000, reaching $m_{\rm{gas}}$ $\sim$ $0.2$ M$_\odot$ (spatial resolution, $r_{\rm{gas,cloud}}$ $\sim$ $O(\rm{pc})$). The selected clouds have diverse properties, across a broad parameter space, resulting in heterogeneous evolution. For the clouds we study, radiative cooling is the dominant physical process enabling cloud survival, while magnetic fields play a comparatively smaller role. The motion of these clouds is governed not only by drag forces that decelerate them, but also by acceleration from momentum exchange with the complex background velocity field, which can cause them to move faster than ballistic projectiles set by their initial velocities. Our results suggest that the non-trivial details of realistic cosmological initial conditions -- specifically the complex density, temperature, and velocity fields -- may play an important role in subsequent cloud evolution, and that sampling the output of an existing large-scale simulation provides a self-consistent approach to capture these effects without ad hoc assumptions. △ Less

Submitted 26 February, 2026; originally announced February 2026.

Comments: Submitted to PASA; 13 pages, 7 figures (main draft)

Hot accretion onto spiral galaxies: the origin of extended and warped HI discs

Authors: Sriram Sankar, Jonathan Stern, Chris Power, Barbara Catinella, Drummond Fielding, Claude-André Faucher-Giguère, Imran Sultan, Michael Boylan-Kolchin, Joss Bland-Hawthorn

Abstract: Gas accretion, hot ($\sim 10^6$ K) atmospheres, and a tilt between the rotation axes of the disc and the atmosphere are all common predictions of standard galaxy evolution theory for massive star-forming galaxies at low redshift. Using idealised hydrodynamic simulations, we demonstrate that the central regions of hot galaxy atmospheres continuously condense into cool ($\sim 10^4$ K) discs, while b… ▽ More Gas accretion, hot ($\sim 10^6$ K) atmospheres, and a tilt between the rotation axes of the disc and the atmosphere are all common predictions of standard galaxy evolution theory for massive star-forming galaxies at low redshift. Using idealised hydrodynamic simulations, we demonstrate that the central regions of hot galaxy atmospheres continuously condense into cool ($\sim 10^4$ K) discs, while being replenished by an inflow from larger scales. The size and orientation of the condensed disc are determined by the angular momentum of the atmosphere, so the condensed disc is expected to often be tilted and more extended than the stellar disc. Continuous accretion from hot atmospheres can thus explain the ubiquity of extended and warped HI discs around local spirals, and also potentially provide the necessary fuel for star formation. This hot accretion scenario predicts the absence of significant HI from galaxy halos, consistent with recent 21 cm constraints on nearby spirals (the so-called `HI desert'). Moreover, our analysis indicates that observations of HI warps can be used to constrain the angular momentum, accretion rate, and gas metallicity of hot galaxy atmospheres, important parameters for disc galaxy evolution that are hard to determine by other means. △ Less

Submitted 3 June, 2026; v1 submitted 5 November, 2025; originally announced November 2025.

Comments: Accepted for publication in MNRAS. 17 pages, 12 figures

Journal ref: Mon Not R Astron Soc (2026)

Accretion-Driven Turbulence in the Circumgalactic Medium

Authors: Roy Goldner, Jonathan Stern, Drummond Fielding, Claude-André Faucher-Giguère, Yakov Faerman, Aharon Kakoly

Abstract: Simulations suggest that turbulence is ubiquitous in the circumgalactic medium (CGM), though the source and properties of CGM turbulence is uncertain. Using analytic considerations and hydrodynamic simulations we study how CGM turbulence is driven by gas accretion, thus providing a baseline for additional turbulence driving processes such as galaxy feedback. We demonstrate that in halos with mass… ▽ More Simulations suggest that turbulence is ubiquitous in the circumgalactic medium (CGM), though the source and properties of CGM turbulence is uncertain. Using analytic considerations and hydrodynamic simulations we study how CGM turbulence is driven by gas accretion, thus providing a baseline for additional turbulence driving processes such as galaxy feedback. We demonstrate that in halos with mass $\sim 10^{10}-10^{12} M_{\odot}$ at $0 < z < 2$, accretion amplifies mild turbulent velocities near the virial radius of $σ_t(R_{\rm vir}) \sim 10 \, {\rm km \, s^{-1}}$ to virial velocities at inner CGM radii, $σ_t(0.1 R_{\rm vir}) \approx v_{\rm vir} \sim 100 \, {\rm km \, s^{-1}}$. Rapid cooling at these inner radii further implies that thermal pressure support is small, and the gas is dominated by the cool and warm ($\sim 10^4-10^5 \, {\rm K}$) phases. Inner CGM energetics in these halos is thus dominated by turbulence, with gas density distributions and velocity structure functions similar to those seen in simulations of isothermal supersonic turbulence, rather than those seen in subsonically turbulent stratified media such as the ICM. The accretion rate in these systems is regulated by the turbulence dissipation rate, in contrast with being regulated by the cooling rate as in more massive halos. We argue that galaxy feedback is unlikely to qualitatively change our conclusions unless it significantly depletes the CGM or continuously injects high specific energy material ($\gg v^2_{\rm vir}$). Such `turbulence-dominated' CGM can be identified in observations via the predicted wide lognormal ionization distributions and large velocity dispersions in UV absorption spectra. △ Less

Submitted 14 April, 2026; v1 submitted 31 October, 2025; originally announced October 2025.

Comments: 21 pages, 16 figures, submitted to MNRAS, Accepted to MNRAS. Minor changes following a positive referee report

Journal ref: Mon Not R Astron Soc (2026)

Modeling Emission-Line Surface Brightness in a Multiphase Galactic Wind: An O VI Case Study

Authors: Zirui Chen, Zixuan Peng, Kate H. R. Rubin, Timothy M. Heckman, Matthew J. Hayes, Yakov Faerman, Crystal L. Martin, S. Peng Oh, Drummond B. Fielding

Abstract: We present a fast and robust analytic framework for predicting surface brightness (SB) of emission lines in galactic winds as a function of radius up to $\sim 100$ kpc out in the circum-galactic medium. We model multiphase structure in galactic winds by capturing emission from both the volume-filling hot phase (T $\sim 10^{6-7}$ K) and turbulent radiative mixing layers that host intermediate tempe… ▽ More We present a fast and robust analytic framework for predicting surface brightness (SB) of emission lines in galactic winds as a function of radius up to $\sim 100$ kpc out in the circum-galactic medium. We model multiphase structure in galactic winds by capturing emission from both the volume-filling hot phase (T $\sim 10^{6-7}$ K) and turbulent radiative mixing layers that host intermediate temperature gas at the boundaries of cold clouds (T $\sim 10^4$ K). Our multiphase framework makes significantly different predictions of emission signatures compared to traditional single-phase models and explains the paucity of OVI SB measurements in the literature. After accounting for ram pressure equilibrium between the cold clouds and hot wind in supersonic outflows, non-equilibrium ionization effects, and energy budgets other than mechanical energy from core-collapse supernovae, our OVI SB predictions qualitatively match observational results. Our framework provides constraints on the optimal galactic wind properties that facilitate OVI emission observations, including star formation rate surface density, hot phase mass loading factor, and thermalization efficiency factor. These constraints are consistent with existing observations and can help inform future target selections. △ Less

Submitted 27 March, 2026; v1 submitted 2 October, 2025; originally announced October 2025.

Comments: 19 pages, 15 figures. Accepted for publication in MNRAS, 27 March 2026. Zirui Chen and Zixuan Peng are co-first authors who made equal contributions to this work

Journal ref: Mon Not R Astron Soc (2026)

Self-Similar Cosmic-Ray Transport in High-Resolution Magnetohydrodynamic Turbulence

Authors: Philipp Kempski, Drummond B. Fielding, Eliot Quataert, Robert J. Ewart, Philipp Grete, Matthew W. Kunz, Alexander A. Philippov, James Stone

Abstract: We study the propagation of cosmic rays (CRs) through a simulation of magnetohydrodynamic (MHD) turbulence at unprecedented resolution of $10{,}240^3$. We drive turbulence that is subsonic and super-Alfvénic, characterized by $δB_{\rm rms}/B_0=2$. The high resolution enables an extended inertial range such that the Alfvén scale $l_A$, where $δB (l_A)\approx B_0$, is well resolved. This allows us t… ▽ More We study the propagation of cosmic rays (CRs) through a simulation of magnetohydrodynamic (MHD) turbulence at unprecedented resolution of $10{,}240^3$. We drive turbulence that is subsonic and super-Alfvénic, characterized by $δB_{\rm rms}/B_0=2$. The high resolution enables an extended inertial range such that the Alfvén scale $l_A$, where $δB (l_A)\approx B_0$, is well resolved. This allows us to properly capture how the cascade transitions from large amplitudes on large scales to small amplitudes on small scales. We find that sharp bends in the magnetic field are key mediators of particle transport even on small scales via resonant curvature scattering. We further find that particle scattering in the turbulence shows strong hints of self-similarity: (1) the diffusion has weak energy dependence over almost two decades in particle energy and (2) the particles' random walk exhibits a broad power-law distribution of collision times such that the diffusion is dominated by the rarest, long-distance excursions. Our results suggest that large-amplitude MHD turbulence can provide efficient scattering over a wide range of CR energies and may help explain many CR observations above a $\sim$TeV: the flattening of the B/C spectrum, the hardening of CR primary spectra and the weak dependence of arrival anisotropy on CR energy. △ Less

Submitted 14 July, 2025; originally announced July 2025.

Comments: Submitted to ApJL

Simulating High-Velocity Clouds in the Observational Plane: An Initial Study with the Smith Cloud

Authors: Lori E. Porter, Matthew Abruzzo, Greg L. Bryan, Mary Putman, Yong Zheng, Drummond Fielding

Abstract: High-velocity clouds (HVCs) may fuel future star formation in the Milky Way, but they must first survive their passage through the hot halo. While recent work has improved our understanding of the survival criterion for cloud-wind interactions, few observational comparisons exist that test this criterion. We therefore present an initial comparison of simulations with the Smith Cloud (SC; $d=$ 12.4… ▽ More High-velocity clouds (HVCs) may fuel future star formation in the Milky Way, but they must first survive their passage through the hot halo. While recent work has improved our understanding of the survival criterion for cloud-wind interactions, few observational comparisons exist that test this criterion. We therefore present an initial comparison of simulations with the Smith Cloud (SC; $d=$ 12.4 kpc, $l, b = 40^{\circ}, -13^{\circ}$) as mapped with the GALFA-HI survey. We use the Smith Cloud's observed properties to motivate simulations of comparable clouds in wind tunnel simulations with Enzo-E, an MHD code. For both observations and simulations, we generate moment maps, characterize turbulence through a projected first-order velocity structure function (VSF), and do the same for HI column density with a normalized autocovariance function. We explore how initial cloud conditions (such as radius, metallicity, thermal pressure, viewing angle, and distance) affect these statistics, demonstrating that the small-scale VSF is sensitive to cloud turbulence while large scales depend on cloud bulk velocity and viewing angle. We find that some simulations reproduce key observational features (particularly the correlation between column density and velocity dispersion) but none match all observational probes at the same time (the large scales of the column density autocovariance is particularly challenging). We find that the simulated cloud (cloud C) showing growth via a turbulent radiative mixing layer (TRML) is the best match, implying the importance of TRML-mediated cooling for Milky Way HVCs. We conclude by suggesting improvements for simulations to better match observed HVCs. △ Less

Submitted 30 May, 2025; originally announced June 2025.

Comments: 19 pages, 12 figures. Submitted to MNRAS

Journal ref: Mon Not R Astron Soc (2025) 1188-1207

Neutral gas phase distribution from HI morphology: phase separation with scattering spectra and variational autoencoders

Authors: Minjie Lei, S. E. Clark, Rudy Morel, E. Allys, Iryna S. Butsky, Caleb Redshaw, Drummond B. Fielding

Abstract: Unraveling the multi-phase structure of the diffuse interstellar medium (ISM) as traced by neutral hydrogen (HI) is essential to understanding the lifecycle of the Milky Way. However, HI phase separation is a challenging and under-constrained problem. The neutral gas phase distribution is often inferred from the spectral line structure of HI emission. In this work, we develop a data-driven phase s… ▽ More Unraveling the multi-phase structure of the diffuse interstellar medium (ISM) as traced by neutral hydrogen (HI) is essential to understanding the lifecycle of the Milky Way. However, HI phase separation is a challenging and under-constrained problem. The neutral gas phase distribution is often inferred from the spectral line structure of HI emission. In this work, we develop a data-driven phase separation method that extracts HI phase structure solely from the spatial morphology of HI emission intensity structures. We combine scattering spectra (SS) statistics with a Gaussian-mixture variational autoencoder (VAE) model to: 1. derive an interpretable statistical model of different HI phases from their multi-scale morphological structures; 2. use this model to decompose the 2D channel maps of GALFA-HI emission in diffuse high latitude ($|b|>30$\degree) regions over narrow velocity channels ($Δv=3$ km/s) into cold neutral medium (CNM), warm neutral medium (WNM), and noise components. We integrate our CNM map over velocity channels to compare it to an existing map produced by a spectrum-based method, and find that the two maps are highly correlated, while ours recovers more spatially coherent structures at small scales. Our work illustrates and quantifies a clear physical connection between the HI morphology and HI phase structure, and unlocks a new avenue for improving future phase separation techniques by making use of both HI spectral and spatial information to decompose HI in 3D position-position-velocity (PPV) space. These results are consistent with a physical picture where processes that drive HI phase transitions also shape the morphology of HI gas, imprinting a sparse, filamentary CNM that forms out of a diffuse, extended WNM. △ Less

Submitted 14 October, 2025; v1 submitted 26 May, 2025; originally announced May 2025.

Turbulence-dominated CGM: the origin of UV absorbers with equivalent widths of $\sim1$Å

Authors: Aharon Kakoly, Jonathan Stern, Claude-André Faucher-Giguère, Drummond B. Fielding, Roy Goldner, Guochao Sun, Cameron B. Hummels

Abstract: Theoretical arguments and observations suggest that in massive halos ($>10^{12}\,M_\odot$), the circumgalactic medium (CGM) is dominated by a 'hot' phase with gas temperature near the virial temperature ($T \approx T_{\rm vir}$) and a quasi-hydrostatic pressure profile. Lower-mass halos are however unlikely to be filled with a similar quasi-static hot phase, due to rapid radiative cooling. Using t… ▽ More Theoretical arguments and observations suggest that in massive halos ($>10^{12}\,M_\odot$), the circumgalactic medium (CGM) is dominated by a 'hot' phase with gas temperature near the virial temperature ($T \approx T_{\rm vir}$) and a quasi-hydrostatic pressure profile. Lower-mass halos are however unlikely to be filled with a similar quasi-static hot phase, due to rapid radiative cooling. Using the FIRE cosmological zoom simulations, we demonstrate that the hot phase is indeed sub-dominant at inner radii ($\lesssim 0.3\,R_{\rm vir}$) of $\lesssim 10^{12}\,M_\odot$ halos, and the inner CGM is instead filled with $T \ll T_{\rm vir}$ gas originating in outflows and inflows, with a turbulent velocity comparable to the halo virial velocity. The turbulent velocity thus exceeds the mass-weighted sound speed in the inner CGM, and the turbulence is supersonic. UV absorption features from such CGM trace the wide lognormal density distributions of the predominantly cool and turbulent volume-filling phase, in contrast with tracing localized cool 'clouds' embedded in a hot medium. We predict equivalent widths of $W_λ\sim 2λv_c/c \sim 1A$ for a broad range of strong UV and EUV transitions (Mg II, C II, C IV, Si II-IV, O III-V) in sightlines through inner CGM dominated by turbulent pressure of $\lesssim L^*$ galaxies at redshifts $0 \leq z \lesssim 2$, where $λ$ is the transition wavelength, $v_{\rm c}$ is the halo circular velocity and $c$ is the speed of light. Comparison of our predictions with observational constraints suggests that star-forming dwarf and $\lesssim L^*$ galaxies are generally dominated by turbulent pressure in their inner CGM, rather than by thermal pressure. The inner CGM surrounding these galaxies is thus qualitatively distinct from that around quenched galaxies and massive disks such as the Milky-Way and M31, in which thermal pressure likely dominates. △ Less

Submitted 24 August, 2025; v1 submitted 23 April, 2025; originally announced April 2025.

Comments: 21 pages, 12 figures (+2 in appendix), submitted to MNRAS. Resubmission following a positive referee report at MNRAS; Fig. 12 has been substantially updated

Applying a star formation model calibrated on high-resolution interstellar medium simulations to cosmological simulations of galaxy formation

Authors: Jan D. Burger, Volker Springel, Eve C. Ostriker, Chang-Goo Kim, Sarah M. R. Jeffreson, Matthew C. Smith, Rüdiger Pakmor, Sultan Hassan, Drummond Fielding, Lars Hernquist, Greg L. Bryan, Rachel S. Somerville, Jake S. Bennett, Rainer Weinberger

Abstract: Modern high-resolution simulations of the interstellar medium (ISM) have shown that key factors in governing star formation are the competing influences of radiative dissipation, pressure support driven by stellar feedback, and the relentless pull of gravity. Cosmological simulations of galaxy formation, such as IllustrisTNG or ASTRID, are however not able to resolve this physics in detail and the… ▽ More Modern high-resolution simulations of the interstellar medium (ISM) have shown that key factors in governing star formation are the competing influences of radiative dissipation, pressure support driven by stellar feedback, and the relentless pull of gravity. Cosmological simulations of galaxy formation, such as IllustrisTNG or ASTRID, are however not able to resolve this physics in detail and therefore need to rely on approximate treatments. These have often taken the form of empirical subgrid models of the ISM expressed in terms of an effective equation of state (EOS) that relates the mean ISM pressure to the mean gas density. Here we seek to improve these heuristic models by directly fitting their key ingredients to results of the high-resolution TIGRESS simulations, which have shown that the dynamical equilibrium of the ISM can be understood in terms of a pressure-regulated, feedback modulated (PRFM) model for star formation. Here we explore a simple subgrid model that draws on the PRFM concept but uses only local quantities. It accurately reproduces PRFM for pure gas disks, while it predicts slightly less star formation than PRFM in the presence of an additional thin stellar disk. We compare the properties of this model with the older Springel and Hernquist and TNG prescriptions, and apply all three to isolated simulations of disk galaxies as well as to a set of high-resolution zoom-in simulations carried out with a novel 'multi-zoom' technique that we introduce in this study. The softer EOS implied by TIGRESS produces substantially thinner disk galaxies, which has important ramifications for disk stability and galaxy morphology. The total stellar mass of galaxies is however hardly modified at low redshift, reflecting the dominating influence of large-scale gaseous inflows and outflows to galaxies, which are not sensitive to the EOS itself △ Less

Submitted 18 February, 2025; originally announced February 2025.

Comments: 22 pages, 21 figures, to be submitted to MNRAS. This is a Learning the Universe publication

Journal ref: Mon Not R Astron Soc (2025) 1390-1411

The Type Ia Supernova and Asymptotic Giant Branch Stellar Ejecta-regulated Interstellar Medium of Massive Galaxies

Authors: Rajsekhar Mohapatra, Eliot Quataert, Drummond Fielding, Minghao Guo

Abstract: Observations and theory suggest that Type Ia supernovae (SNIa) heating and mass loss from asymptotic giant branch (AGB) stars play a crucial role in the interstellar medium (ISM) of massive galaxies. We perform 3D hydrodynamic simulations of the central few kiloparsecs of massive galaxies, including radiative cooling and mass and energy injection from AGB winds and SNIa (resolving each SNIa remnan… ▽ More Observations and theory suggest that Type Ia supernovae (SNIa) heating and mass loss from asymptotic giant branch (AGB) stars play a crucial role in the interstellar medium (ISM) of massive galaxies. We perform 3D hydrodynamic simulations of the central few kiloparsecs of massive galaxies, including radiative cooling and mass and energy injection from AGB winds and SNIa (resolving each SNIa remnant, a few $\times10~\mathrm{pc}$ in size), excluding black hole feedback. We study systems with different initial core thermodynamic profiles, focusing on NGC 1399. Our simulations reproduce its observed density and entropy profiles well. Over $100~\mathrm{Myr}$, two steady-state profiles emerge, depending on the inner circumgalactic medium (CGM) pressure and the ratio of Ia heating to cooling: (i) if SNIa heating is less than cooling, a cooling flow develops; (ii) if SNIa heating is comparable to or exceeds cooling, SNIa heating drives a slow subsonic outflow of AGB ejecta, with black hole accretion at small radii. This outflow, pressure-confined by the CGM, adapts the ISM to the CGM properties: a low entropy CGM results in a dense, low entropy ISM with higher black hole accretion, while a high entropy CGM leads to a less dense, high entropy ISM with lower accretion. This suggests that the AGB-SNIa regulated ISM connects CGM and galaxy scales, potentially influencing black hole feedback in massive halos. Approximate methods of modeling Ia heating, such as clustered SNIa and smoothly distributed heating, produce unrealistic ISM profiles over $100~\mathrm{Myr}$, highlighting the importance of resolving SNIa in simulations. △ Less

Submitted 13 August, 2025; v1 submitted 7 February, 2025; originally announced February 2025.

Comments: Updated to APJ accepted version, simulation movies at: https://youtube.com/playlist?list=PLuaNgQ1v_KMZh_O4fznaKMyvOEmhUav_T&si=9HbIdFQMIkDhd9dt

Journal ref: ApJ 989 103 (2025)

Physical Origins of Outflowing Cold Clouds in Local Star-forming Dwarf Galaxies

Authors: Zixuan Peng, Crystal L. Martin, Zirui Chen, Drummond B. Fielding, Xinfeng Xu, Timothy Heckman, Lise Ramambason, Yuan Li, Cody Carr, Weida Hu, Zuyi Chen, Claudia Scarlata, Alaina Henry

Abstract: We study the physical origins of outflowing cold clouds in a sample of 14 low-redshift dwarf ($M_{\ast} \lesssim 10^{10}$ $M_{\odot}$) galaxies from the COS Legacy Archive Spectroscopic SurveY (CLASSY) using Keck/ESI data. Outflows are traced by broad (FWHM ~ 260 $\rm{km}$ $\rm{s^{-1}}$) and very-broad (VB; FWHM ~ 1200 $\rm{km}$ $\rm{s^{-1}}$) velocity components in strong emission lines like [O I… ▽ More We study the physical origins of outflowing cold clouds in a sample of 14 low-redshift dwarf ($M_{\ast} \lesssim 10^{10}$ $M_{\odot}$) galaxies from the COS Legacy Archive Spectroscopic SurveY (CLASSY) using Keck/ESI data. Outflows are traced by broad (FWHM ~ 260 $\rm{km}$ $\rm{s^{-1}}$) and very-broad (VB; FWHM ~ 1200 $\rm{km}$ $\rm{s^{-1}}$) velocity components in strong emission lines like [O III] $λ5007$ and $\rm{H}α$. The maximum velocities ($v_{\rm{max}}$) of broad components correlate positively with SFR, unlike the anti-correlation observed for VB components, and are consistent with superbubble models. In contrast, supernova-driven galactic wind models better reproduce the $v_{\rm{max}}$ of VB components. Direct radiative cooling from a hot wind significantly underestimates the luminosities of both broad and VB components. A multi-phase wind model with turbulent radiative mixing reduces this discrepancy to at least one dex for most VB components. Stellar photoionization likely provides additional energy since broad components lie in the starburst locus of excitation diagnostic diagrams. We propose a novel interpretation of outflow origins in star-forming dwarf galaxies$-$broad components trace expanding superbubble shells, while VB components originate from galactic winds. One-zone photoionization models fail to explain the low-ionization lines ([S II] and [O I]) of broad components near the maximal starburst regime, which two-zone photoionization models with density-bounded channels instead reproduce. These two-zone models indicate anisotropic leakage of Lyman continuum photons through low-density channels formed by expanding superbubbles. Our study highlights extreme outflows ($v_{\rm{max}} \gtrsim 1000$ $\rm{km}$ $\rm{s^{-1}}$) in 9 out of 14 star-forming dwarf galaxies, comparable to AGN-driven winds. △ Less

Submitted 6 December, 2024; originally announced December 2024.

Comments: 42 pages, 22 figures, 7 tables. This is the revised version submitted to ApJ, incorporating changes based on the first round of referee feedback

Journal ref: The Astrophysical Journal 981 (2025) 2

A Unified Model of Cosmic Ray Propagation and Radio Extreme Scattering Events from Intermittent Interstellar Structures

Authors: Philipp Kempski, Dongzi Li, Drummond B. Fielding, Eliot Quataert, E. Sterl Phinney, Matthew W. Kunz, Dylan L. Jow, Alexander A. Philippov

Abstract: Intermittent magnetic structures are a plausible candidate for explaining cosmic-ray (CR) diffusion rates derived from observed CR energy spectra. Independently, studies of extreme scattering events (ESEs) of radio quasars and pulsar scintillation have hinted that very straight, large-aspect-ratio, magnetic current sheets may be responsible for the localized large scattering of radio waves. The re… ▽ More Intermittent magnetic structures are a plausible candidate for explaining cosmic-ray (CR) diffusion rates derived from observed CR energy spectra. Independently, studies of extreme scattering events (ESEs) of radio quasars and pulsar scintillation have hinted that very straight, large-aspect-ratio, magnetic current sheets may be responsible for the localized large scattering of radio waves. The required shortest axis of the typical structures producing ESEs is of the same scale ($\sim$AU) as the gyroradii of $\sim$GeV CRs. In this paper, we propose that the same magnetic/density sheets can produce large scattering of both CRs and radio waves. We demonstrate that the geometry and volume filling factor of the sheets derived from quasar ESEs can explain the observed mean free path of GeV CRs without introducing free parameters. The model places constraints on the sheet geometry, such as straightness and large aspect ratio, and assumes the statistics of the sheets are similar throughout the Galactic volume. We, therefore, discuss observational tests of the sheet model, which includes observations of echoes in pulsars and fast radio bursts, gravitationally lensed quasars, the distribution of ESE durations, and spatial correlations between ESE events and rotation-measure fluctuations. Such tests will be enabled by upcoming wide-field radio instruments, including Canadian Hydrogen Observatory and Radio-transient Detector (CHORD) and Deep Synoptic Array 2000 Antennas (DSA-2000). △ Less

Submitted 4 December, 2024; originally announced December 2024.

Comments: 21 Pages, 9 Figures. Comments welcome

The Well: a Large-Scale Collection of Diverse Physics Simulations for Machine Learning

Authors: Ruben Ohana, Michael McCabe, Lucas Meyer, Rudy Morel, Fruzsina J. Agocs, Miguel Beneitez, Marsha Berger, Blakesley Burkhart, Keaton Burns, Stuart B. Dalziel, Drummond B. Fielding, Daniel Fortunato, Jared A. Goldberg, Keiya Hirashima, Yan-Fei Jiang, Rich R. Kerswell, Suryanarayana Maddu, Jonah Miller, Payel Mukhopadhyay, Stefan S. Nixon, Jeff Shen, Romain Watteaux, Bruno Régaldo-Saint Blancard, François Rozet, Liam H. Parker , et al. (2 additional authors not shown)

Abstract: Machine learning based surrogate models offer researchers powerful tools for accelerating simulation-based workflows. However, as standard datasets in this space often cover small classes of physical behavior, it can be difficult to evaluate the efficacy of new approaches. To address this gap, we introduce the Well: a large-scale collection of datasets containing numerical simulations of a wide va… ▽ More Machine learning based surrogate models offer researchers powerful tools for accelerating simulation-based workflows. However, as standard datasets in this space often cover small classes of physical behavior, it can be difficult to evaluate the efficacy of new approaches. To address this gap, we introduce the Well: a large-scale collection of datasets containing numerical simulations of a wide variety of spatiotemporal physical systems. The Well draws from domain experts and numerical software developers to provide 15TB of data across 16 datasets covering diverse domains such as biological systems, fluid dynamics, acoustic scattering, as well as magneto-hydrodynamic simulations of extra-galactic fluids or supernova explosions. These datasets can be used individually or as part of a broader benchmark suite. To facilitate usage of the Well, we provide a unified PyTorch interface for training and evaluating models. We demonstrate the function of this library by introducing example baselines that highlight the new challenges posed by the complex dynamics of the Well. The code and data is available at https://github.com/PolymathicAI/the_well. △ Less

Submitted 21 February, 2025; v1 submitted 30 November, 2024; originally announced December 2024.

Comments: 38th Conference on Neural Information Processing Systems (NeurIPS 2024) Track on Datasets and Benchmarks

Prevention is better than cure? Feedback from high specific energy winds in cosmological simulations with Arkenstone

Authors: Jake S. Bennett, Matthew C. Smith, Drummond B. Fielding, Greg L. Bryan, Chang-Goo Kim, Volker Springel, Lars Hernquist, Rachel S. Somerville, Laura Sommovigo

Abstract: We deploy the new Arkenstone galactic wind model in cosmological simulations for the first time, allowing us to robustly resolve the evolution and impact of high specific energy winds. In a (25 $h^{-1}$ Mpc)$^3$ box we perform a set of numerical experiments that systematically vary the mass and energy loadings of such winds, finding that their energy content is the key parameter controlling the st… ▽ More We deploy the new Arkenstone galactic wind model in cosmological simulations for the first time, allowing us to robustly resolve the evolution and impact of high specific energy winds. In a (25 $h^{-1}$ Mpc)$^3$ box we perform a set of numerical experiments that systematically vary the mass and energy loadings of such winds, finding that their energy content is the key parameter controlling the stellar to dark matter mass ratio. Increasing the mass loading, at fixed energy, actually results in mildly enhanced star formation, counter to prevailing wisdom, due to the wind becoming cooler. Of the simple parametrisations that we test, we find that an energy loading that scales inversely with halo mass best matches a wide range of observations and can do so with mass loadings drastically lower than those in most previous cosmological simulations. In this scenario, much less material is ejected from the interstellar medium. Instead, winds both heat gas in the circumgalactic medium, slowing infall onto the galaxy, and also drive shocks beyond the virial radius, decreasing the halo-scale accretion rate. We can also report that a much lower fraction of the available supernova energy is needed in preventative galaxy regulation than required by ejective wind feedback models such as IllustrisTNG. This is a Learning the Universe collaboration publication. △ Less

Submitted 8 September, 2025; v1 submitted 16 October, 2024; originally announced October 2024.

Comments: Accepted in MNRAS. 23 pages, 19 figures. Update includes additional plots investigating mass fluxes and metals in our simulations. This is a Learning the Universe publication

Journal ref: Mon Not R Astron Soc (2026)

The Metallicity Mapping of the Ionized Diffuse Gas at the Milky Way Disk-halo Interface

Authors: Bo-Eun Choi, Jessica K. Werk, Kirill Tchernyshyov, J. Xavier Prochaska, Yong Zheng, Mary E. Putman, Drummond B. Fielding, Jay Strader

Abstract: Metals in the diffuse, ionized gas at the boundary between the Milky Way's interstellar medium (ISM) and circumgalactic medium (CGM), known as the disk-halo interface (DHI), are valuable tracers of the feedback processes that drive the Galactic fountain. However, metallicity measurements in this region are challenging due to obscuration by the Milky Way ISM and uncertain ionization corrections tha… ▽ More Metals in the diffuse, ionized gas at the boundary between the Milky Way's interstellar medium (ISM) and circumgalactic medium (CGM), known as the disk-halo interface (DHI), are valuable tracers of the feedback processes that drive the Galactic fountain. However, metallicity measurements in this region are challenging due to obscuration by the Milky Way ISM and uncertain ionization corrections that affect the total hydrogen column density. In this work, we constrain the ionization corrections to neutral hydrogen column densities using precisely measured electron column densities from the dispersion measure of pulsars that lie in the same globular clusters as UV-bright targets with high-resolution absorption spectroscopy. We address the blending of absorption lines with the ISM by jointly fitting Voigt profiles to all absorption components. We present our metallicity estimates for the DHI of the Milky Way based on detailed photoionization modeling to the absorption from ionized metal lines and ionization-corrected total hydrogen columns. Generally, the gas clouds show a large scatter in metallicity, ranging between $0.04-3.2\ Z_{\odot}$, implying that the DHI consists of a mixture of gaseous structures having multiple origins. We estimate the inflow and outflow timescales of the DHI ionized clouds to be $6 - 35$ Myr. We report the detection of an infalling cloud with super-solar metallicity that suggests a Galactic fountain mechanism, whereas at least one low-metallicity outflowing cloud ($Z < 0.1\ Z_{\odot}$) poses a challenge for Galactic fountain and feedback models. △ Less

Submitted 8 October, 2024; originally announced October 2024.

Comments: 28 pages, 15 figures, accepted for publication in ApJ

AthenaK: A Performance-Portable Version of the Athena++ AMR Framework

Authors: James M. Stone, Patrick D. Mullen, Drummond Fielding, Philipp Grete, Minghao Guo, Philipp Kempski, Elias R. Most, Christopher J. White, George N. Wong

Abstract: We describe AthenaK: a new implementation of the Athena++ block-based adaptive mesh refinement (AMR) framework using the Kokkos programming model. Finite volume methods for Newtonian, special relativistic (SR), and general relativistic (GR) hydrodynamics and magnetohydrodynamics (MHD), and GR-radiation hydrodynamics and MHD, as well as a module for evolving Lagrangian tracer or charged test partic… ▽ More We describe AthenaK: a new implementation of the Athena++ block-based adaptive mesh refinement (AMR) framework using the Kokkos programming model. Finite volume methods for Newtonian, special relativistic (SR), and general relativistic (GR) hydrodynamics and magnetohydrodynamics (MHD), and GR-radiation hydrodynamics and MHD, as well as a module for evolving Lagrangian tracer or charged test particles (e.g., cosmic rays) are implemented using the framework. In two companion papers we describe (1) a new solver for the Einstein equations based on the Z4c formalism and (2) a GRMHD solver in dynamical spacetimes also implemented using the framework, enabling new applications in numerical relativity. By adopting Kokkos, the code can be run on virtually any hardware, including CPUs, GPUs from multiple vendors, and emerging ARM processors. AthenaK shows excellent performance and weak scaling, achieving over one billion cell updates per second for hydrodynamics in three-dimensions on a single NVIDIA Grace Hopper processor and with a typical parallel efficiency of 80% on 65536 AMD GPUs on the OLCF Frontier system. Such performance portability enables AthenaK to leverage modern exascale computing systems for challenging applications in astrophysical fluid dynamics, numerical relativity, and multimessenger astrophysics. △ Less

Submitted 24 September, 2024; originally announced September 2024.

Comments: Submitted to ApJS

The Outflowing [OII] Nebulae of Compact Starburst Galaxies at z $\sim$ 0.5

Authors: Serena Perrotta, Alison L. Coil, David S. N. Rupke, Wenmeng Ning, Brendan Duong, Aleksandar M. Diamond-Stanic, Drummond B. Fielding, James E. Geach, Ryan C. Hickox, John Moustakas, Gregory H. Rudnick, Paul H. Sell, Cameren N. Swiggum, Christy A. Tremonti

Abstract: High-velocity outflows are ubiquitous in compact, massive (M$_* \sim$ 10$^{11}$ M$_{\odot}$), z $\sim$ 0.5 galaxies with extreme star formation surface densities ($Σ_{SFR} \sim$ 2000 M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$). We have previously detected and characterized these outflows using MgII absorption lines. To probe their full extent, we present Keck/KCWI integral field spectroscopy of the [OII] an… ▽ More High-velocity outflows are ubiquitous in compact, massive (M$_* \sim$ 10$^{11}$ M$_{\odot}$), z $\sim$ 0.5 galaxies with extreme star formation surface densities ($Σ_{SFR} \sim$ 2000 M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$). We have previously detected and characterized these outflows using MgII absorption lines. To probe their full extent, we present Keck/KCWI integral field spectroscopy of the [OII] and MgII emission nebulae surrounding all of the 12 galaxies in this study. We find that [OII] is more effective than MgII in tracing low surface brightness, extended emission in these galaxies. The [OII] nebulae are spatially extended beyond the stars, with radial extent R$_{90}$ between 10 and 40 kpc. The nebulae exhibit non-gravitational motions, indicating galactic outflows with maximum blueshifted velocities ranging from -335 to -1920 km s$^{-1}$. The outflow kinematics correlate with the bursty star formation histories of these galaxies. Galaxies with the most recent bursts of star formation (within the last $<$ 3 Myr) exhibit the highest central velocity dispersions ($σ>$ 400 km s$^{-1}$), while the oldest bursts have the lowest-velocity outflows. Many galaxies exhibit both high-velocity cores and more extended, slower-moving gas indicative of multiple outflow episodes. The slower, larger outflows occurred earlier and have decelerated as they propagate into the CGM and mix on timescales $>$ 50 Myr. △ Less

Submitted 16 September, 2024; originally announced September 2024.

Comments: Accepted for publication in The Astrophysical Journal

Towards Implementation of the Pressure-Regulated, Feedback-Modulated Model of Star Formation in Cosmological Simulations: Methods and Application to TNG

Authors: Sultan Hassan, Eve C. Ostriker, Chang-Goo Kim, Greg L. Bryan, Jan D. Burger, Drummond B. Fielding, John C. Forbes, Shy Genel, Lars Hernquist, Sarah M. R. Jeffreson, Bhawna Motwani, Matthew C. Smith, Rachel S. Somerville, Ulrich P. Steinwandel, Romain Teyssier

Abstract: Traditional star formation subgrid models implemented in cosmological galaxy formation simulations, such as that of Springel & Hernquist (2003, hereafter SH03), employ adjustable parameters to satisfy constraints measured in the local Universe. In recent years, however, theory and spatially-resolved simulations of the turbulent, multiphase, star-forming ISM have begun to produce new first-principl… ▽ More Traditional star formation subgrid models implemented in cosmological galaxy formation simulations, such as that of Springel & Hernquist (2003, hereafter SH03), employ adjustable parameters to satisfy constraints measured in the local Universe. In recent years, however, theory and spatially-resolved simulations of the turbulent, multiphase, star-forming ISM have begun to produce new first-principles models, which when fully developed can replace traditional subgrid prescriptions. This approach has advantages of being physically motivated and predictive rather than empirically tuned, and allowing for varying environmental conditions rather than being tied to local Universe conditions. As a prototype of this new approach, by combining calibrations from the TIGRESS numerical framework with the Pressure-Regulated Feedback-Modulated (PRFM) theory, simple formulae can be obtained for both the gas depletion time and an effective equation of state. Considering galaxies in TNG50, we compare the "native" simulation outputs with post-processed predictions from PRFM. At TNG50 resolution, the total midplane pressure is nearly equal to the total ISM weight, indicating that galaxies in TNG50 are close to satisfying vertical equilibrium. The measured gas scale height is also close to theoretical equilibrium predictions. The slopes of the effective equations of states are similar, but with effective velocity dispersion normalization from SH03 slightly larger than that from current TIGRESS simulations. Because of this and the decrease in PRFM feedback yield at high pressure, the PRFM model predicts shorter gas depletion times than the SH03 model at high densities and redshift. Our results represent a first step towards implementing new, numerically calibrated subgrid algorithms in cosmological galaxy formation simulations. △ Less

Submitted 13 September, 2024; originally announced September 2024.

Comments: 33 pages, 12 figures, Accepted for publication in ApJ. This is a Learning the Universe Publication. All codes and data used to produce this work can be found at the following $\href{https://github.com/sultan-hassan/tng50-post-processing-prfm}{GitHub \,Link.}$

Arkenstone -- II. A model for unresolved cool clouds entrained in galactic winds in cosmological simulations

Authors: Matthew C. Smith, Drummond B. Fielding, Greg L. Bryan, Jake S. Bennett, Chang-Goo Kim, Eve C. Ostriker, Rachel S. Somerville

Abstract: Arkenstone is a new scheme that allows multiphase, stellar feedback-driven winds to be included in coarse resolution cosmological simulations. The evolution of galactic winds and their subsequent impact on the circumgalactic medium are altered by exchanges of mass, energy, momentum, and metals between their component phases. These exchanges are governed by complex, small-scale physical processes t… ▽ More Arkenstone is a new scheme that allows multiphase, stellar feedback-driven winds to be included in coarse resolution cosmological simulations. The evolution of galactic winds and their subsequent impact on the circumgalactic medium are altered by exchanges of mass, energy, momentum, and metals between their component phases. These exchanges are governed by complex, small-scale physical processes that cannot be resolved in cosmological simulations. In this second presentation paper, we describe Arkenstone's novel cloud particle approach for modelling unresolvable cool clouds entrained in hot, fast winds. This general framework allows models of the cloud-wind interaction, derived from state-of-the-art high-resolution simulations, to be applied in a large-scale context. In this work, we adopt a cloud evolution model that captures simultaneous cloud mass loss to and gain from the ambient hot phase via turbulent mixing and radiative cooling, respectively. We demonstrate the scheme using non-cosmological idealized simulations of a galaxy with a realistic circumgalactic medium component, using the Arepo code. We show that the ability of a high-specific energy wind component to perform preventative feedback may be limited by heavy loading of cool clouds coupled into it. We demonstrate that the diverging evolution of clouds of initially differing masses leads to a complex velocity field for the cool phase and a cloud mass function that varies both spatially and temporally in a non-trivial manner. These latter two phenomena can manifest in the simulation because of our choice of a Lagrangian discretisation of the cloud population, in contrast to other proposed schemes. This is a Learning the Universe publication. △ Less

Submitted 27 August, 2024; originally announced August 2024.

Comments: Submitted to MNRAS. 27 pages, 11 figures. This is a Learning the Universe publication

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.

Comments: 19 pages, 11 figures, submitted to ApJ, comments welcome!

Zooming in on the Circumgalactic Medium with GIBLE: Tracing the Origin and Evolution of Cold Clouds

Authors: Rahul Ramesh, Dylan Nelson, Drummond Fielding, Marcus Brüggen

Abstract: We use the GIBLE suite of cosmological zoom-in simulations of Milky Way-like galaxies with additional super-Lagrangian refinement in the circumgalactic medium (CGM) to quantify the origin and evolution of CGM cold gas clouds. The origin of $z$\,$=$\,$0$ clouds can be traced back to recent ($\lesssim$\,$2$\,Gyr) outflows from the central galaxy ($\sim$\,45\,$\%$), condensation out of the hot phase… ▽ More We use the GIBLE suite of cosmological zoom-in simulations of Milky Way-like galaxies with additional super-Lagrangian refinement in the circumgalactic medium (CGM) to quantify the origin and evolution of CGM cold gas clouds. The origin of $z$\,$=$\,$0$ clouds can be traced back to recent ($\lesssim$\,$2$\,Gyr) outflows from the central galaxy ($\sim$\,45\,$\%$), condensation out of the hot phase of the CGM in the same time frame ($\sim$\,45\,$\%$), and to a lesser degree to satellite galaxies ($\lesssim$\,5\,$\%$). We find that in-situ condensation results from rapid cooling around local over-densities primarily seeded by the dissolution of the previous generation of clouds into the hot halo. About $\lesssim$\,10\,$\%$ of the cloud population is long lived, with their progenitors having already assembled $\sim$\,$2$\,Gyr ago. Collective cloud-cloud dynamics are crucial to their evolution, with coalescence and fragmentation events occurring frequently ($\gtrsim$\,20\,Gyr$^{-1}$). These interactions are modulated by non-vanishing pressure imbalances between clouds and their interface layers. The gas content of clouds is in a constant state of flux, with clouds and their surroundings exchanging mass at a rate of \mbox{$\gtrsim$\,$10^3$\,M$_\odot$\,Myr$^{-1}$}, depending on cloud relative velocity and interface vorticity. Furthermore, we find that a net magnetic tension force acting against the density gradient is capable of inhibiting cloud-background mixing. Our results show that capturing the distinct origins of cool CGM clouds, together with their physical evolution, requires high-resolution, cosmological galaxy formation simulations with both stellar and supermassive black hole feedback-driven outflows. △ Less

Submitted 28 June, 2024; originally announced July 2024.

Comments: Submitted to A&A

Journal ref: A&A 699, A10 (2025)

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.

Comments: 15 pages, 3 figures, Submitted to ApJ

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.

Comments: 18 pages, 5 figures, submitted to ApJ

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.

Comments: Resubmitted to ApJ after minor revision

JWST Observations of Starbursts: Cold Clouds and Plumes Launching in the M82 Outflow

Authors: Deanne B. Fisher, Alberto D. Bolatto, John Chisholm, Drummond Fielding, Rebecca C. Levy, Elizabeth Tarantino, Martha L. Boyer, Serena A. Cronin, Laura A. Lopez, J. D. Smith, Danielle A. Berg, Sebastian Lopez, Sylvain Veilleux, Paul P. van der Werf, Torsten Böker, Leindert A. Boogaard, Laura Lenkić, Simon C. O. Glover, Vicente Villanueva, Divakara Mayya, Thomas S. -Y. Lai, Daniel A. Dale, Kimberly L. Emig, Fabian Walter, Monica Relaño , et al. (6 additional authors not shown)

Abstract: In this paper we study the filamentary substructure of 3.3 $μ$m PAH emission from JWST/NIRCam observations in the base of the M82 star-burst driven wind. We identify plume-like substructure within the PAH emission with widths of $\sim$50 pc. Several of the plumes extend to the edge of the field-of-view, and thus are at least 200-300 pc in length. In this region of the outflow, the vast majority (… ▽ More In this paper we study the filamentary substructure of 3.3 $μ$m PAH emission from JWST/NIRCam observations in the base of the M82 star-burst driven wind. We identify plume-like substructure within the PAH emission with widths of $\sim$50 pc. Several of the plumes extend to the edge of the field-of-view, and thus are at least 200-300 pc in length. In this region of the outflow, the vast majority ($\sim$70\%) of PAH emission is associated with the plumes. We show that those structures contain smaller scale "clouds" with widths that are $\sim$5-15 pc, and they are morphologically similar to the results of "cloud-crushing" simulations. We estimate the cloud-crushing time-scales of $\sim$0.5-3 Myr, depending on assumptions. We show this time scale is consistent with a picture in which these observed PAH clouds survived break-out from the disk rather than being destroyed by the hot wind. The PAH emission in both the midplane and the outflow is shown to tightly correlate with that of Pa$α$ emission (from HST/NICMOS data), at the scale of both plumes and clouds, though the ratio of PAH-to-Pa$α$ increases at further distances from the midplane. Finally, we show that the outflow PAH emission is suppressed in regions of the M82 wind that are bright in X-ray emission. Overall, our results are broadly consistent with a picture in which cold gas in galactic outflows is launched via hierarchically structured plumes, and those small scale clouds are more likely to survive the wind environment when collected into the larger plume structure. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: Submitted to MNRAS

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.

Comments: 15 pages, 11 figures, accepted by MNRAS

Zooming in on the Circumgalactic Medium with GIBLE: the Topology and Draping of Magnetic Fields around Cold Clouds

Authors: Rahul Ramesh, Dylan Nelson, Drummond Fielding, Marcus Brüggen

Abstract: We use a cosmological zoom-in simulation of a Milky Way-like galaxy to study and quantify the topology of magnetic field lines around cold gas clouds in the circumgalactic medium (CGM). This simulation is a new addition to Project GIBLE, a suite of cosmological magnetohydrodynamical simulations of galaxy formation with preferential super-Lagrangian refinement in the CGM, reaching an unprecedented… ▽ More We use a cosmological zoom-in simulation of a Milky Way-like galaxy to study and quantify the topology of magnetic field lines around cold gas clouds in the circumgalactic medium (CGM). This simulation is a new addition to Project GIBLE, a suite of cosmological magnetohydrodynamical simulations of galaxy formation with preferential super-Lagrangian refinement in the CGM, reaching an unprecedented (CGM) gas mass resolution of $\sim$ $225$ M$_\odot$. To maximize statistics and resolution, we focus on a sample of $\sim$ $200$ clouds with masses of $\sim$ $10^6$ M$_\odot$. The topology of magnetic field lines around clouds is diverse, from threading to draping, and there is large variation in the magnetic curvature ($κ$) within cloud-background interfaces. We typically find little variation of $κ$ between upstream and downstream cloud faces, implying that strongly draped configurations are rare. In addition, $κ$ correlates strongly with multiple properties of the interface and the ambient background, including cloud overdensity and relative velocity, suggesting that cloud properties impact the topology of interface magnetic fields. △ Less

Submitted 1 April, 2024; originally announced April 2024.

Comments: Accepted for publication (A&A Letters); 5 pages, 4 figures

Journal ref: A&A 684, L16 (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 (Deep near-UV observations of Entrained gas in Turbulent galaxies) sample, using Keck/KCWI observations of H$β$ and [OIII]~$λ$5007. We measure $\sim460$ lines-of-sight that contain outflows, and use these to study scaling relationships of outflow velocity ($v_{\rm out}$), mass-loading fact… ▽ More We measure resolved (kiloparsec-scale) outflow properties in a sample of 10 starburst galaxies from the DUVET (Deep near-UV observations of Entrained gas in Turbulent galaxies) sample, using Keck/KCWI observations of H$β$ and [OIII]~$λ$5007. We measure $\sim460$ 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.5}$. 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 18 November, 2024; v1 submitted 27 February, 2024; originally announced February 2024.

Comments: 23 pages, 9 figures, plus 5 figures in the Appendix. Supplementary table of results in Ancillary Files. Resubmitted to MNRAS after referee report

Observational Signatures of AGN Feedback in the Morphology and the Ionization States of Milky Way-like Galaxies

Authors: Nadia Qutob, Razieh Emami, Kung-Yi Su, Randall Smith, Lars Hernquist, Dian P. Triani, Cameron Hummels, Drummond Fielding, Philip F. Hopkins, Rachel S. Somerville, David R. Ballantyne, Mark Vogelsberger, Grant Tremblay, James F. Steiner, Douglas Finkbeiner, Ramesh Narayan, Minjung Park, Josh Grindlay, Priyamvada Natarajan, Christopher C. Hayward, Dušan Kereš, Sam B. Ponnada, Sirio Belli, Rebecca Davies, Gabriel Maheson , et al. (2 additional authors not shown)

Abstract: We make an in-depth analysis of different AGN jet models' signatures, inducing quiescence in galaxies with a halo mass of $10^{12} M_\odot$. Three jet models, including cosmic ray-dominant, hot thermal, and precessing kinetic jets, are studied at two energy flux levels each, compared to a jet-free, stellar feedback-only simulation. We examine the distribution of Mg II, O VI, and O VIII ions, along… ▽ More We make an in-depth analysis of different AGN jet models' signatures, inducing quiescence in galaxies with a halo mass of $10^{12} M_\odot$. Three jet models, including cosmic ray-dominant, hot thermal, and precessing kinetic jets, are studied at two energy flux levels each, compared to a jet-free, stellar feedback-only simulation. We examine the distribution of Mg II, O VI, and O VIII ions, alongside gas temperature and density profiles. Low-energy ions, like Mg II, concentrate in the ISM, while higher energy ions, e.g., O VIII, prevail at the AGN jet cocoon's edge. High-energy flux jets display an isotropic ion distribution with lower overall density. High-energy thermal or cosmic ray jets pressurize at smaller radii, significantly suppressing core density. The cosmic ray jet provides extra pressure support, extending cool and warm gas distribution. A break in the ion-to-mass ratio slope in O VI and O VIII is demonstrated in the ISM-to-CGM transition (between 10-30 kpc), growing smoothly towards the CGM at greater distances. △ Less

Submitted 22 December, 2023; originally announced December 2023.

Comments: 17 pages, 13 figures

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.

Comments: 22 pages, 7 figures. Submitted to AAS Journals, with minor modifications. Comments welcome. (1) Co-first authors who made equal contributions to this work

Unraveling Jet Quenching Criteria Across L* Galaxies and Massive Cluster Ellipticals

Authors: Kung-Yi Su, Greg L. Bryan, Christopher C. Hayward, Rachel S. Somerville, Philip F. Hopkins, Razieh Emami, Claude-André Faucher-Giguère, Eliot Quataert, Sam B. Ponnada, Drummond Fielding, Dušan Kereš

Abstract: In the absence of supplementary heat, the radiative cooling of halo gas around massive galaxies (Milky Way mass and above) leads to an excess of cold gas or stars beyond observed levels. AGN jet-induced heating is likely essential, but the specific properties of the jets remain unclear. Our previous work (Su et al. 2021) concludes from simulations of a halo with $10^{14} M_\odot$ that a successful… ▽ More In the absence of supplementary heat, the radiative cooling of halo gas around massive galaxies (Milky Way mass and above) leads to an excess of cold gas or stars beyond observed levels. AGN jet-induced heating is likely essential, but the specific properties of the jets remain unclear. Our previous work (Su et al. 2021) concludes from simulations of a halo with $10^{14} M_\odot$ that a successful jet model should have an energy flux comparable to the free-fall energy flux at the cooling radius and should inflate a sufficiently wide cocoon with a long enough cooling time. In this paper, we investigate three jet modes with constant fluxes satisfying the criteria, including high-temperature thermal jets, cosmic ray (CR)-dominant jets, and widely precessing kinetic jets in $10^{12}-10^{15}\,{\rm M}_{\odot}$ halos using high-resolution, non-cosmological MHD simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, conduction, and viscosity. We find that scaling the jet energy according to the free-fall energy at the cooling radius can successfully suppress the cooling flows and quench galaxies without obviously violating observational constraints. We investigate an alternative scaling method in which we adjust the energy flux based on the total cooling rate within the cooling radius. However, we observe that the strong interstellar medium (ISM) cooling dominates the total cooling rate in this scaling approach, resulting in a jet flux that exceeds the amount needed to suppress the cooling flows. With the same energy flux, the CR-dominant jet is most effective in suppressing the cooling flow across all the surveyed halo masses due to the enhanced CR pressure support. We confirm that the criteria for a successful jet model, which we proposed in Su et al. (2021), work across a much wider range, encompassing halo masses of $10^{12}-10^{15} {\rm M_\odot}$. △ Less

Submitted 2 November, 2023; v1 submitted 26 October, 2023; originally announced October 2023.

Comments: 15 pages, 10 figures, Submitted to MNRAS; comments welcome!

Ionized Gas Extended Over 40 kpc in an Odd Radio Circle Host Galaxy

Authors: Alison L. Coil, Serena Perrotta, David S. N. Rupke, Cassandra Lochhaas, Christy A. Tremonti, Aleks Diamond-Stanic, Drummond Fielding, James Geach, Ryan C. Hickox, John Moustakas, Gregory H. Rudnick, Paul Sell, Kelly E. Whalen

Abstract: A new class of extragalactic astronomical sources discovered in 2021, named Odd Radio Circles (ORCs, Norris et al. 2021), are large rings of faint, diffuse radio continuum emission spanning ~1 arcminute on the sky. Galaxies at the centers of several ORCs have photometric redshifts of z~0.3-0.6, implying physical scales of several 100 kiloparsecs in diameter for the radio emission, the origin of wh… ▽ More A new class of extragalactic astronomical sources discovered in 2021, named Odd Radio Circles (ORCs, Norris et al. 2021), are large rings of faint, diffuse radio continuum emission spanning ~1 arcminute on the sky. Galaxies at the centers of several ORCs have photometric redshifts of z~0.3-0.6, implying physical scales of several 100 kiloparsecs in diameter for the radio emission, the origin of which is unknown. Here we report spectroscopic data on an ORC including strong [OII] emission tracing ionized gas in the central galaxy of ORC4 at z=0.4512. The physical extent of the [OII] emission is ~40 kpc in diameter, larger than expected for a typical early-type galaxy (Pandya et al, 2017) but an order of magnitude smaller than the large-scale radio continuum emission. We detect a ~200 km/s velocity gradient across the [OII] nebula, as well as a high velocity dispersion of ~180 km/s. The [OII] equivalent width (EW, ~50 Ang) is extremely high for a quiescent galaxy. The morphology, kinematics, and strength of the [OII] emission are consistent with the infall of shock ionized gas near the galaxy, following a larger-scale, outward moving shock driven by a galactic wind. Both the extended optical and radio emission, while observed on very different scales, may therefore result from the same dramatic event. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 7 figures, accepted to Nature

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.

Comments: 17 pages, 12 figures, 2 tables Accepted for publication in MNRAS on Oct 6, 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.

Comments: 12 Pages, 6 Figures

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.

Comments: 6.5 pages, 4 figures, submitted to ApJL

Accretion onto disk galaxies via hot and rotating CGM inflows

Authors: Jonathan Stern, Drummond Fielding, Zachary Hafen, Kung-Yi Su, Nadav Naor, Claude-André Faucher-Giguère, Eliot Quataert, James Bullock

Abstract: Observed accretion rates onto the Milky-Way and other local spirals fall short of that required to sustain star formation for cosmological timescales. A potential avenue for this unseen accretion is an inflow in the volume-filling hot phase ($\sim10^6$ K) of the circumgalactic medium (CGM), as suggested by some cosmological simulations. Using hydrodynamic simulations and a new analytic solution va… ▽ More Observed accretion rates onto the Milky-Way and other local spirals fall short of that required to sustain star formation for cosmological timescales. A potential avenue for this unseen accretion is an inflow in the volume-filling hot phase ($\sim10^6$ K) of the circumgalactic medium (CGM), as suggested by some cosmological simulations. Using hydrodynamic simulations and a new analytic solution valid in the slow-rotation limit, we show that a hot inflow spins up as it approaches the galaxy, while remaining hot, subsonic and quasi-spherical. At the radius of angular momentum support ($\approx15$ kpc for the Milky-Way) the hot flow flattens into a disk geometry and then cools from $\sim10^6$ K to $\sim10^4$ K at the disk-halo interface. Cooling affects all hot gas, rather than just a subset of individual gas clouds, implying that accretion via hot inflows does not rely on local thermal instability in contrast with 'precipitation' models for galaxy accretion. Prior to cooling and accretion the inflow completes $\sim t_{\rm cool}/t_{\rm ff}$ radians of rotation, where $t_{\rm cool}/t_{\rm ff}$ is the cooling time to free-fall time ratio in hot gas immediately outside the galaxy. The ratio $t_{\rm cool}/t_{\rm ff}$ may thus govern the development of turbulence and enhancement of magnetic fields in gas accreting onto low-redshift spirals. We argue that accretion via hot inflows can explain the observed truncation of nearby thin stellar disks at $\approx4$ disk radii. We also show that if rotating hot inflows are common in Milky-Way size disk galaxies, as predicted, then signatures should be observable with X-ray telescopes and FRB surveys. △ Less

Submitted 15 February, 2024; v1 submitted 31 May, 2023; originally announced June 2023.

Comments: 21 pages, 12 figures, resubmitted to MNRAS following referee report

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.

Comments: 34 pages, 37 figures; Accepted for publication in MNRAS

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.

Comments: Published in MNRAS, 16 pages, 9 Figures

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.

Comments: Published MNRAS, 27 pages, 17 figures; updated for consistency with journal version

Journal ref: MNRAS, 527, 1216 (2024)

A framework for improving the accessibility of research papers on arXiv.org

Authors: Shamsi Brinn, Christopher Cameron, David Fielding, Charles Frankston, Alison Fromme, Peter Huang, Mark Nazzaro, Stephanie Orphan, Steinn Sigurdsson, Ryan Tay, Miranda Yang, Qianyu Zhou

Abstract: The research content hosted by arXiv is not fully accessible to everyone due to disabilities and other barriers. This matters because a significant proportion of people have reading and visual disabilities, it is important to our community that arXiv is as open as possible, and if science is to advance, we need wide and diverse participation. In addition, we have mandates to become accessible, and… ▽ More The research content hosted by arXiv is not fully accessible to everyone due to disabilities and other barriers. This matters because a significant proportion of people have reading and visual disabilities, it is important to our community that arXiv is as open as possible, and if science is to advance, we need wide and diverse participation. In addition, we have mandates to become accessible, and accessible content benefits everyone. In this paper, we will describe the accessibility problems with research, review current mitigations (and explain why they aren't sufficient), and share the results of our user research with scientists and accessibility experts. Finally, we will present arXiv's proposed next step towards more open science: offering HTML alongside existing PDF and TeX formats. An accessible HTML version of this paper is also available at https://info.arxiv.org/about/accessibility_research_report.html △ Less

Submitted 2 February, 2024; v1 submitted 14 December, 2022; originally announced December 2022.

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.

Comments: submitted to ApJ, 26 pages, 14 figures, comments welcome

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.

Comments: Updated to reflect version accepted by ApJ; minor text edits to clarify limitations and caveats of model