Showing 1–50 of 144 results for author: Scannapieco, E

The hydrodynamic response of small-scale structure to reionization drives large IGM temperature fluctuations that persist to z = 4

Authors: Christopher Cain, Evan Scannapieco, Matthew McQuinn, Anson D'Aloisio, Hy Trac

Abstract: The thermal history and structure of the intergalactic medium (IGM) at $z \geq 4$ is a key boundary condition for reionization, which has been measured using the Ly$α$ forest of high-redshift quasars. It is also a key input for studies that use the forest to constrain the particle masses of alternative dark matter candidates. Most such inferences rely on simulations that lack the spatial resolutio… ▽ More The thermal history and structure of the intergalactic medium (IGM) at $z \geq 4$ is a key boundary condition for reionization, which has been measured using the Ly$α$ forest of high-redshift quasars. It is also a key input for studies that use the forest to constrain the particle masses of alternative dark matter candidates. Most such inferences rely on simulations that lack the spatial resolution to fully resolve the hydrodynamic response of IGM filaments and minihalos caused by HI reionization heating. In this letter, we use high-resolution hydrodynamic+radiative transfer simulations to study the impact of these on the IGM thermal structure. We find that the adiabatic heating and cooling driven by the expansion of initially cold gas filaments and minihalos drives significant temperature fluctuations on small scales. These likely persist in much of the IGM until at least $z = 4$. Capturing this effect requires resolving the characteristic clumping scale of cold, pre-ionized gas, which demands spatial resolutions of at least $2$ $h^{-1}$kpc. Pre-heating of the IGM by X-Ray sources can slightly reduce the effect. Our preliminary estimate of the effect on the Ly$α$ forest finds that, at $\log(k /[{\rm km^{-1} s}]) = -1.0$, the forest flux power (at fixed mean flux) can increase by $10-20\%$ when going from $8$ and $2$ $h^{-1}$kpc resolution at $z = 4-5$ for gas ionized at $z < 7$. These findings motivate a more careful analysis of how temperature fluctuations driven by pressure smoothing from reionization affect the Ly$α$ forest. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: 6+1 pages, 2+1 figures, to be submitted to MNRAS letters. Comments welcome

The imprint of magnetic fields on absorption spectra from circumgalactic wind-cloud systems

Authors: Benedetta Casavecchia, Wladimir Banda-Barragan, Marcus Brueggen, Fabrizio Brighenti, Evan Scannapieco

Abstract: Galactic winds probe how stellar feedback regulates the mass and metallicity of galaxies. Observations show that galactic winds are multiphase and magnetised. In the local Universe, the dense phase is traced by emission and absorption lines, which reveal the presence of fast-moving clouds embedded in hot streams. Simulations tell us that magnetic fields can shield such clouds and delay their disru… ▽ More Galactic winds probe how stellar feedback regulates the mass and metallicity of galaxies. Observations show that galactic winds are multiphase and magnetised. In the local Universe, the dense phase is traced by emission and absorption lines, which reveal the presence of fast-moving clouds embedded in hot streams. Simulations tell us that magnetic fields can shield such clouds and delay their disruption, but there is little discussed on their observational effects. Using 3D MHD simulations, we study the influence of two orientations of the magnetic field (aligned and transverse) on the cloud morphology, temperature and density structure, mixing fraction, ion kinematics, column densities, and absorption spectra. We study supersonic wind-cloud systems with radiative processes, and develop a framework to extract ion column density maps and synthetic absorption spectra. The framework relies on studying ion populations and creating down-the-barrel spectra via an interface that links our PLUTO simulations to TRIDENT using YT, CLOUDY, and STARBURST99. We find that the transverse magnetic field makes the cloud asymmetric, shields and protects dense cold gas, and reduces mixing fractions compared to the aligned case. Ions can reach higher velocities in the transverse field case. The imprints of the initial orientation of the field on the synthetic spectra are: in the cold phase we find no signature of C ii and Si ii when the field is aligned, in the intermediate phase traced by C iv and Si iv we find broader lines in the transverse case, and in the warm phase we find deeper lines for O vi and N v in the aligned case, but they are less sensitive to the field orientation. Magnetic fields significantly affect the absorption spectra of cold clouds. Intermediate ions are the most sensitive to the magnetic field orientation and can potentially yield information about magnetic field topology. △ Less

Submitted 8 July, 2024; v1 submitted 2 February, 2024; originally announced February 2024.

Comments: 18 pages, 13 figures

Distinguishing AGN Feedback Models with the Thermal Sunyaev-Zel'dovich Effect

Authors: Skylar Grayson, Evan Scannapieco, Romeel Davé

Abstract: Current models of galaxy formation require strong feedback from active galactic nuclei (AGN) to explain the observed lack of star formation in massive galaxies since z~2 but direct evidence of this energy input is limited. We use the SIMBA cosmological galaxy formation simulations to assess the ability of thermal Sunyaev-Zel'dovich (tSZ) measurements to provide such evidence, by mapping the pressu… ▽ More Current models of galaxy formation require strong feedback from active galactic nuclei (AGN) to explain the observed lack of star formation in massive galaxies since z~2 but direct evidence of this energy input is limited. We use the SIMBA cosmological galaxy formation simulations to assess the ability of thermal Sunyaev-Zel'dovich (tSZ) measurements to provide such evidence, by mapping the pressure structure of the circumgalactic medium around massive z~0.2-1.5 galaxies. We undertake a stacking approach to calculate the total tSZ signal and its radial profile in simulations with varying assumptions of AGN feedback, and we assess its observability with current and future telescopes. By convolving our predictions with the 2.1' beam of the Atacama Cosmology Telescope (ACT), we show that current observations at z~1 are consistent with SIMBA's fiducial treatment of AGN feedback, and inconsistent with SIMBA models without feedback. At z~0.5, observational signals lie between SIMBA run with and without AGN feedback, suggesting AGN in SIMBA may inject too much energy at late times. By convolving our data with a 9.5'' beam corresponding to the TolTEC camera on the Large Millimeter Telescope Alfonso Serrano (LMT), we predict a unique profile for AGN feedback that can be distinguished with future higher-resolution measurements. Finally, we explore a novel approach to quantify the non-spherically symmetric features surrounding our galaxies by plotting radial profiles representing the component of the stack with m-fold symmetry. △ Less

Submitted 2 October, 2023; originally announced October 2023.

Comments: 16 pages, 13 figures, 3 tables. Accepted by ApJ

Constraining Circumgalactic Turbulence with QSO Absorption-line Measurements

Authors: Brad Koplitz, Edward Buie II, Evan Scannapieco

Abstract: Our knowledge of the circumgalactic medium (CGM) is mostly based on quasar absorption-line measurements. These have uncovered a multiphase medium that is likely highly turbulent, but constraints of this turbulence are limited to measurements of the non-thermal width of absorption-line components ($b_{turb}$) and the line-of-sight velocity dispersion between components ($σ_{LOS}$). Here we analyze… ▽ More Our knowledge of the circumgalactic medium (CGM) is mostly based on quasar absorption-line measurements. These have uncovered a multiphase medium that is likely highly turbulent, but constraints of this turbulence are limited to measurements of the non-thermal width of absorption-line components ($b_{turb}$) and the line-of-sight velocity dispersion between components ($σ_{LOS}$). Here we analyze a suite of CGM simulations to determine how well these indirect measures are related to the underlying CGM. Our simulations track the non-equilibrium evolution of all commonly observed ions, and consist of two main types: small-scale simulations of regions of homogeneous CGM turbulence and global simulations of inhomogeneous turbulence throughout a galactic halo. From each simulation, we generate mock spectra of Si II, Si IV, C IV, and O VI, which allow us to directly compare $b_{turb}$ and $σ_{LOS}$ to the true line-of-sight turbulence ($σ_{1D}$). In the small-scale simulations, $b_{turb}$ is only weakly correlated with $σ_{1D}$, likely because it measures random motions within individual warm CGM clouds, which do not sample the overall random motions. Meanwhile, $σ_{LOS}$ and $σ_{1D}$ are strongly correlated, with $σ_{1D}\approxσ_{LOS}+10$ km s$^{-1}$ in the densest regions we simulated, though, the strength of this correlation depended weakly on the gas phase being probed. Our large-scale simulations also indicate that $b_{turb}$ and $σ_{1D}$ are largely uncorrelated, and that $σ_{1D}\approxσ_{LOS}+10$ kms$^{-1}$ on average, although it varies along individual sightlines. Moreover, the $σ_\mathrm{LOS}$ distributions from our global simulations are similar to recent observations, suggesting that this quantity may provide useful constraints on circumgalactic turbulence regardless of the axis probed. △ Less

Submitted 23 August, 2023; originally announced August 2023.

Comments: 16 pages, 9 figures, 2 tables, ApJ, in press

Evidence of Extended Dust and Feedback around z$\approx$1 Quiescent Galaxies via Millimeter Observations

Authors: Jeremy Meinke, Seth Cohen, Jenna Moore, Kathrin Böckmann, Philip Mauskopf, Evan Scannapieco

Abstract: We use public data from the South Pole Telescope (SPT) and Atacama Cosmology Telescope (ACT) to measure radial profiles of the thermal Sunyaev-Zel'dovich (tSZ) effect and dust emission around massive quiescent galaxies at $z\approx1.$ Using survey data from the Dark Energy Survey (DES) and Wide-Field infrared Survey Explorer (WISE), we selected $387,627$ quiescent galaxies within the ACT field, wi… ▽ More We use public data from the South Pole Telescope (SPT) and Atacama Cosmology Telescope (ACT) to measure radial profiles of the thermal Sunyaev-Zel'dovich (tSZ) effect and dust emission around massive quiescent galaxies at $z\approx1.$ Using survey data from the Dark Energy Survey (DES) and Wide-Field infrared Survey Explorer (WISE), we selected $387,627$ quiescent galaxies within the ACT field, with a mean stellar $\log_{10}(M_{\star}/\rm{M_{\odot}})$ of $11.40$. A subset of $94,452$ galaxies, with a mean stellar $\log_{10}(M_{\star}/\rm{M_{\odot}})$ of $11.36,$ are also covered by SPT. In $0.5$ arcminute radial bins around these galaxies, we detect the tSZ profile at levels up to $11σ$, and dust profile up to $20σ.$ Both profiles are extended, and the dust profile slope at large radii is consistent with galaxy clustering. We analyze the thermal energy and dust mass versus stellar mass via integration within $R=2.0$ arcminute circular apertures and fit them with a forward-modeled power-law to correct for our photometric stellar mass uncertainty. At the mean log stellar mass of our overlap and wide-area samples, respectively, we extract thermal energies from the tSZ of $E_{\rm{pk}}=6.45_{-1.52}^{+1.67}\times10^{60}{\rm{ erg}}$ and $8.20_{-0.52}^{+0.52}\times10^{60}{\rm{ erg}},$ most consistent with moderate to high levels of active galactic nuclei feedback acting upon the circumgalactic medium. Dust masses at the mean log stellar mass are $M_{\rm{d,pk}}=6.23_{-0.67}^{+0.67}\times10^{8}\rm{ M_{\odot}}$ and $6.76_{-0.56}^{+0.56}\times10^{8}\rm{ M_{\odot}},$ and we find a greater than linear dust-to-stellar mass relation, which indicates that the more massive galaxies in our study retain more dust. Our work highlights current capabilities of stacking millimeter data around individual galaxies and potential for future use. △ Less

Submitted 7 June, 2023; originally announced June 2023.

Comments: 23 pages, 12 figures, accepted for publication in ApJ

The Launching of Cold Clouds by Galaxy Outflows V: The Role of Anisotropic Thermal Conduction

Authors: M. Brüggen, E. Scannapieco, P. Grete

Abstract: Motivated by observations of multiphase galaxy outflows, we explore the impact of isotropic and anisotropic electron thermal conduction on the evolution of radiatively-cooled, cold clouds embedded in hot, magnetized winds. Using the adaptive mesh refinement code AthenaPK, we conduct simulations of clouds impacted by supersonic and transonic flows with magnetic fields initially aligned parallel and… ▽ More Motivated by observations of multiphase galaxy outflows, we explore the impact of isotropic and anisotropic electron thermal conduction on the evolution of radiatively-cooled, cold clouds embedded in hot, magnetized winds. Using the adaptive mesh refinement code AthenaPK, we conduct simulations of clouds impacted by supersonic and transonic flows with magnetic fields initially aligned parallel and perpendicular to the flow direction. In cases with isotropic thermal conduction, an evaporative wind forms, stabilizing against instabilities and leading to a mass loss rate that matches the hydrodynamic case. In anisotropic cases, the impact of conduction is more limited and strongly dependent on the field orientation. In runs with initially perpendicular fields, the field lines are folded back into the tail, strongly limiting conduction, but magnetic fields act to dampen instabilities and slow the stretching of the cloud in the flow direction. In the parallel case, anisotropic conduction aids cloud survival by forming a radiative wind near the front of the cloud, which suppresses instabilities and reduces mass loss. In all cases, anisotropic conduction has a minimal impact on the acceleration of the cloud. △ Less

Submitted 16 May, 2023; v1 submitted 19 April, 2023; originally announced April 2023.

Comments: accepted by ApJ

The Effects of Radiative Feedback and Supernova Induced Turbulence on Early Galaxies

Authors: Richard Sarmento, Evan Scannapieco

Abstract: The recently launched James Webb Space Telescope promises unparalleled advances in our understanding of the first stars and galaxies, but realizing this potential requires cosmological simulations that capture the key physical processes that affected these objects. Here we show that radiative transfer and subgrid turbulent mixing are two such processes. By comparing simulations with and without ra… ▽ More The recently launched James Webb Space Telescope promises unparalleled advances in our understanding of the first stars and galaxies, but realizing this potential requires cosmological simulations that capture the key physical processes that affected these objects. Here we show that radiative transfer and subgrid turbulent mixing are two such processes. By comparing simulations with and without radiative transfer but with exactly the same physical parameters and subgrid turbulent mixing model, we show that tracking radiative transfer suppresses the Population III (Pop III) star formation density by a factor of approximately 4. In both simulations, $\gtrsim 90\%$ of Pop III stars are found in the unresolved pristine regions tracked by our subgrid model, which does a better job at modeling the regions surrounding proto-galaxy cores where metals from supernovae take tens of Myrs to mix thoroughly. At the same time, radiative transfer suppresses Pop III star formation, via the development of ionized bubbles that slows gas accretion in these regions, and it results in compact high-redshift galaxies that are surrounded by isolated low mass satellites. Thus turbulent mixing and radiative transfer are both essential processes that must be included to accurately model the morphology, composition, and growth of primordial galaxies. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: Accepted by ApJ, Jul 2022, Comments welcome

Modeling Photoionized Turbulent Material in the Circumgalactic Medium III: Effects of Co-rotation and Magnetic Fields

Authors: Edward Buie II, Evan Scannapieco, G. Mark Voit

Abstract: Absorption-line measurements of the circumgalactic medium (CGM) display a highly non-uniform distribution of lower ionization state species accompanied by more widespread higher ionization state material. This suggests that the CGM is a dynamic, multiphase medium, such as arises in the presence of turbulence. To better understand this evolution, we perform hydrodynamic and magneto-hydrodynamic (MH… ▽ More Absorption-line measurements of the circumgalactic medium (CGM) display a highly non-uniform distribution of lower ionization state species accompanied by more widespread higher ionization state material. This suggests that the CGM is a dynamic, multiphase medium, such as arises in the presence of turbulence. To better understand this evolution, we perform hydrodynamic and magneto-hydrodynamic (MHD) simulations of the CGM surrounding Milky Way-like galaxies. In both cases, the CGM is initially in hydrostatic balance in a $10^{12}$ solar masses dark matter gravitational potential, and the simulations include rotation in the inner halo and turbulence that decreases radially. They also track ionizations, recombinations, and species-by-species radiative cooling in the presence of the redshift-zero UV background, employing the MAIHEM non-equilibrium chemistry package. We find that after 9 Gyrs of evolution, the presence of a magnetic field leads to an overall hotter CGM, with cool gas in the center where magnetic pressure dominates. While the non-MHD run produces more cold clouds overall, we find similar Si IV/O VI and N V/O VI ratios between the MHD and non-MHD runs, which are both very different from their equilibrium values. The non-MHD halo develops cool, low angular momentum filaments above the central disk, in comparison to the MHD run that has more efficient angular momentum transport, especially for the cold gas which forms a more ordered and extended disk late into its evolution. △ Less

Submitted 23 January, 2022; originally announced January 2022.

Comments: 24 pages, 16 figures

A New Model For Including Galactic Winds in Simulations of Galaxy Formation II: Implementation of PhEW in Cosmological Simulations

Authors: Shuiyao Huang, Neal Katz, J'Neil Cottle, Evan Scannapieco, Romeel Davé, David H. Weinberg

Abstract: Although galactic winds play a critical role in regulating galaxy formation, hydrodynamic cosmological simulations do not resolve the scales that govern the interaction between winds and the ambient circumgalactic medium (CGM). We implement the Physically Evolved Wind (PhEW) model of Huang et al. (2020) in the GIZMO hydrodynamics code and perform test cosmological simulations with different choice… ▽ More Although galactic winds play a critical role in regulating galaxy formation, hydrodynamic cosmological simulations do not resolve the scales that govern the interaction between winds and the ambient circumgalactic medium (CGM). We implement the Physically Evolved Wind (PhEW) model of Huang et al. (2020) in the GIZMO hydrodynamics code and perform test cosmological simulations with different choices of model parameters and numerical resolution. PhEW adopts an explicit subgrid model that treats each wind particle as a collection of clouds that exchange mass, metals, and momentum with their surroundings and evaporate by conduction and hydrodynamic instabilities as calibrated on much higher resolution cloud scale simulations. In contrast to a conventional wind algorithm, we find that PhEW results are robust to numerical resolution and implementation details because the small scale interactions are defined by the model itself. Compared to conventional wind simulations with the same resolution, our PhEW simulations produce similar galaxy stellar mass functions at $z\geq 1$ but are in better agreement with low-redshift observations at $M_* < 10^{11}M_\odot$ because PhEW particles shed mass to the CGM before escaping low mass halos. PhEW radically alters the CGM metal distribution because PhEW particles disperse metals to the ambient medium as their clouds dissipate, producing a CGM metallicity distribution that is skewed but unimodal and is similar between cold and hot gas. While the temperature distributions and radial profiles of gaseous halos are similar in simulations with PhEW and conventional winds, these changes in metal distribution will affect their predicted UV/X-ray properties in absorption and emission. △ Less

Submitted 22 February, 2022; v1 submitted 2 June, 2021; originally announced June 2021.

Comments: 23 pages, 17 figures, MNRAS accepted

The Thermal Sunyaev-Zel'dovich Effect from Massive, Quiescent 0.5 $\leq$ z $\leq$ 1.5 Galaxies

Authors: Jeremy Meinke, Kathrin Böckmann, Seth Cohen, Philip Mauskopf, Evan Scannapieco, Richard Sarmento, Emily Lunde, J'Neil Cottle

Abstract: We use combined South Pole Telescope (SPT)+Planck temperature maps to analyze the circumgalactic medium (CGM) encompassing 138,235 massive, quiescent 0.5 $\leq$ z $\leq$ 1.5 galaxies selected from data from the Dark Energy Survey (DES) and Wide-Field Infrared Survey Explorer (WISE). Images centered on these galaxies were cut from the 1.85 arcmin resolution maps with frequency bands at 95, 150, and… ▽ More We use combined South Pole Telescope (SPT)+Planck temperature maps to analyze the circumgalactic medium (CGM) encompassing 138,235 massive, quiescent 0.5 $\leq$ z $\leq$ 1.5 galaxies selected from data from the Dark Energy Survey (DES) and Wide-Field Infrared Survey Explorer (WISE). Images centered on these galaxies were cut from the 1.85 arcmin resolution maps with frequency bands at 95, 150, and 220 GHz. The images were stacked, filtered, and fit with a gray-body dust model to isolate the thermal Sunyaev-Zel'dovich (tSZ) signal, which is proportional to the total energy contained in the CGM of the galaxies. We separate these $M_{\star} = 10^{10.9} M_\odot$ - $10^{12} M_\odot$ galaxies into 0.1 dex stellar mass bins, detecting tSZ per bin up to $5.6σ$ and a total signal-to-noise ratio of $10.1σ$. We also detect dust with an overall signal-to-noise ratio of $9.8σ$, which overwhelms the tSZ at 150GHz more than in other lower-redshift studies. We correct for the $0.16$ dex uncertainty in the stellar mass measurements by parameter fitting for an unconvolved power-law energy-mass relation, $E_{\rm therm} = E_{\rm therm,peak} \left(M_\star/M_{\star,{\rm peak}} \right)^α$, with the peak stellar mass distribution of our selected galaxies defined as $M_{\star,{\rm peak}}= 2.3 \times 10^{11} M_\odot$. This yields an $E_{\rm therm,peak}= 5.98_{-1.00}^{+1.02} \times 10^{60}$ erg and $α=3.77_{-0.74}^{+0.60}$. These are consistent with $z \approx 0$ observations and within the limits of moderate models of active galactic nuclei (AGN) feedback. We also compute the radial profile of our full sample, which is similar to that recently measured at lower-redshift by Schaan et al. (2021). △ Less

Submitted 20 June, 2023; v1 submitted 1 March, 2021; originally announced March 2021.

Comments: 23 pages, 9 figures, Published in ApJ. Corrected y-axis label/units of Fig. 5 and fixed minor syntax/typos

Shock-multicloud interactions in galactic outflows -- II. Radiative fractal clouds and cold gas thermodynamics

Authors: Wladimir Banda-Barragán, Marcus Brüggen, Volker Heesen, Evan Scannapieco, J'Neil Cottle, Christoph Federrath, Alexander Y. Wagner

Abstract: Galactic winds are crucial to the cosmic cycle of matter, transporting material out of the dense regions of galaxies. Observations show the coexistence of different temperature phases in such winds, which is not easy to explain. We present a set of 3D shock-multicloud simulations that account for radiative heating and cooling at temperatures between $10^2\,\rm K$ and $10^7\,\rm K$. The interplay b… ▽ More Galactic winds are crucial to the cosmic cycle of matter, transporting material out of the dense regions of galaxies. Observations show the coexistence of different temperature phases in such winds, which is not easy to explain. We present a set of 3D shock-multicloud simulations that account for radiative heating and cooling at temperatures between $10^2\,\rm K$ and $10^7\,\rm K$. The interplay between shock heating, dynamical instabilities, turbulence, and radiative heating and cooling creates a complex multi-phase flow with a rain-like morphology. Cloud gas fragments and is continuously eroded, becoming efficiently mixed and mass loaded. The resulting warm mixed gas then cools down and precipitates into new dense cloudlets, which repeat the process. Thus, radiative cooling is able to sustain fast-moving dense gas by aiding condensation of gas from warm clouds and the hot wind. In the ensuing outflow, hot gas with temperatures $\gtrsim 10^6\,\rm K$ outruns the warm and cold phases, which reach thermal equilibrium near $\approx 10^4\,\rm K$ and $\approx 10^2\,\rm K$, respectively. Although the volume filling factor of hot gas is higher in the outflow, most of the mass is concentrated in dense gas cloudlets and filaments with these temperatures. More porous multicloud layers result in more vertically extended outflows, and dense gas is more efficiently produced in more compact layers. The cold phase is not accelerated by ram-pressure, but, instead, precipitates from warm and mixed gas out of thermal equilibrium. This cycle can explain the presence of high-velocity H\,{\sc i} gas with $N_{\rm H\,{\scriptstyle I}}=10^{19-21}\,\rm cm^{-2}$ and $Δv_{\rm FWHM}\lesssim37\,\rm km\,s^{-1}$ in the Galactic centre outflow. △ Less

Submitted 29 June, 2021; v1 submitted 10 November, 2020; originally announced November 2020.

Comments: Accepted for publication in MNRAS, 20 pages, 13 figures. This version includes an additional clump-finding analysis in the Appendix

The Launching of Cold Clouds by Galaxy Outflows. IV. Cosmic-Ray-Driven Acceleration

Authors: Marcus Brüggen, Evan Scannapieco

Abstract: We carry out a suite of simulations of the evolution of cosmic-ray (CR) driven, radiatively-cooled cold clouds embedded in hot material, as found in galactic outflows. In such interactions, CRs stream towards the cloud at the Alfvén speed, which decreases dramatically at the cloud boundary, leading to a bottleneck in which pressure builds up in front of the cloud. At the same time, CRs stream alon… ▽ More We carry out a suite of simulations of the evolution of cosmic-ray (CR) driven, radiatively-cooled cold clouds embedded in hot material, as found in galactic outflows. In such interactions, CRs stream towards the cloud at the Alfvén speed, which decreases dramatically at the cloud boundary, leading to a bottleneck in which pressure builds up in front of the cloud. At the same time, CRs stream along the sides of the cloud, forming a boundary layer where large filaments develop. Shear in this boundary layer is the primary mode of cloud destruction, which is relatively slow in all cases, but slowest in the cases with the lowest Alfvén speeds. Thus the CR ray pressure in the bottleneck region has sufficient time to accelerate the cold clouds efficiently. Furthermore, radiative cooling has relatively little impact on these interactions. Our simulations are two-dimensional and limited by a simplified treatment of cosmic-ray dynamics, the neglect of CR heating, and an idealized magnetic field geometry. Nevertheless, our results suggest that cosmic rays, when acting as the primary source of momentum input, are capable of accelerating clouds to velocities comparable to those observed in galaxy outflows. △ Less

Submitted 14 October, 2020; originally announced October 2020.

Comments: 14 pages, 10 figures, ApJ in press

Shock-multicloud interactions in galactic outflows -- I. Cloud layers with log-normal density distributions

Authors: Wladimir Banda-Barragán, Marcus Brüggen, Christoph Federrath, Alexander Y. Wagner, Evan Scannapieco, J'Neil Cottle

Abstract: We report three-dimensional hydrodynamical simulations of shocks (${\cal M_{\rm shock}}\geq 4$) interacting with fractal multicloud layers. The evolution of shock-multicloud systems consists of four stages: a shock-splitting phase in which reflected and refracted shocks are generated, a compression phase in which the forward shock compresses cloud material, an expansion phase triggered by internal… ▽ More We report three-dimensional hydrodynamical simulations of shocks (${\cal M_{\rm shock}}\geq 4$) interacting with fractal multicloud layers. The evolution of shock-multicloud systems consists of four stages: a shock-splitting phase in which reflected and refracted shocks are generated, a compression phase in which the forward shock compresses cloud material, an expansion phase triggered by internal heating and shock re-acceleration, and a mixing phase in which shear instabilities generate turbulence. We compare multicloud layers with narrow ($σ_ρ=1.9\barρ$) and wide ($σ_ρ=5.9\barρ$) log-normal density distributions characteristic of Mach $\approx 5$ supersonic turbulence driven by solenoidal and compressive modes. Our simulations show that outflowing cloud material contains imprints of the density structure of their native environments. The dynamics and disruption of multicloud systems depend on the porosity and the number of cloudlets in the layers. `Solenoidal' layers mix less, generate less turbulence, accelerate faster, and form a more coherent mixed-gas shell than the more porous `compressive' layers. Similarly, multicloud systems with more cloudlets quench mixing via a shielding effect and enhance momentum transfer. Mass loading of diffuse mixed gas is efficient in all models, but direct dense gas entrainment is highly inefficient. Dense gas only survives in compressive clouds, but has low speeds. If normalised with respect to the shock-passage time, the evolution shows invariance for shock Mach numbers $\geq10$ and different cloud-generating seeds, and slightly weaker scaling for lower Mach numbers and thinner cloud layers. Multicloud systems also have better convergence properties than single-cloud systems, with a resolution of $8$ cells per cloud radius being sufficient to capture their overall dynamics. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Comments: Accepted for publication in MNRAS; 19 pages plus an appendix; 14 figures

Limits to Rest-Frame Ultraviolet Emission From Far-Infrared-Luminous z~6 Quasar Hosts

Authors: Madeline A. Marshall, Mira Mechtley, Rogier A. Windhorst, Seth H. Cohen, Rolf A. Jansen, Linhua Jiang, Victoria R. Jones, J. Stuart B. Wyithe, Xiaohui Fan, Nimish P. Hathi, Knud Jahnke, William C. Keel, Anton M. Koekemoer, Victor Marian, Keven Ren, Jenna Robinson, Huub J. A. Röttgering, Russell E. Ryan Jr., Evan Scannapieco, Donald P. Schneider, Glenn Schneider, Brent M. Smith, Haojing Yan

Abstract: We report on a Hubble Space Telescope search for rest-frame ultraviolet emission from the host galaxies of five far-infrared-luminous $z\simeq{}6$ quasars and the $z=5.85$ hot-dust free quasar SDSS J0005-0006. We perform 2D surface brightness modeling for each quasar using a Markov-Chain Monte-Carlo estimator, to simultaneously fit and subtract the quasar point source in order to constrain the und… ▽ More We report on a Hubble Space Telescope search for rest-frame ultraviolet emission from the host galaxies of five far-infrared-luminous $z\simeq{}6$ quasars and the $z=5.85$ hot-dust free quasar SDSS J0005-0006. We perform 2D surface brightness modeling for each quasar using a Markov-Chain Monte-Carlo estimator, to simultaneously fit and subtract the quasar point source in order to constrain the underlying host galaxy emission. We measure upper limits for the quasar host galaxies of $m_J>22.7$ mag and $m_H>22.4$ mag, corresponding to stellar masses of $M_\ast<2\times10^{11}M_\odot$. These stellar mass limits are consistent with the local $M_{\textrm{BH}}$-$M_\ast$ relation. Our flux limits are consistent with those predicted for the UV stellar populations of $z\simeq6$ host galaxies, but likely in the presence of significant dust ($\langle A_{\mathrm{UV}}\rangle\simeq 2.6$ mag). We also detect a total of up to 9 potential $z\simeq6$ quasar companion galaxies surrounding five of the six quasars, separated from the quasars by 1.4''-3.2'', or 8.4-19.4 kpc, which may be interacting with the quasar hosts. These nearby companion galaxies have UV absolute magnitudes of -22.1 to -19.9 mag, and UV spectral slopes $β$ of -2.0 to -0.2, consistent with luminous star-forming galaxies at $z\simeq6$. These results suggest that the quasars are in dense environments typical of luminous $z\simeq6$ galaxies. However, we cannot rule out the possibility that some of these companions are foreground interlopers. Infrared observations with the James Webb Space Telescope will be needed to detect the $z\simeq6$ quasar host galaxies and better constrain their stellar mass and dust content. △ Less

Submitted 27 July, 2020; originally announced July 2020.

Comments: 22 pages, 13 figures. Accepted for publication in ApJ

Modeling Photoionized Turbulent Material in the Circumgalactic Medium II: Effect of Turbulence within a Stratified Medium

Authors: Edward Buie II, William J. Gray, Evan Scannapieco, Mohammadtaher Safarzadeh

Abstract: The circumgalactic medium (CGM) of nearby star-forming galaxies shows clear indications of OVI absorption accompanied by little to no detectable NV absorption. This unusual spectral signature, accompanied by highly non-uniform absorption from lower ionization state species, indicates that the CGM must be viewed as a dynamic, multiphase medium, such as occurs in the presence of turbulence. Motivate… ▽ More The circumgalactic medium (CGM) of nearby star-forming galaxies shows clear indications of OVI absorption accompanied by little to no detectable NV absorption. This unusual spectral signature, accompanied by highly non-uniform absorption from lower ionization state species, indicates that the CGM must be viewed as a dynamic, multiphase medium, such as occurs in the presence of turbulence. Motivated by previous isotropic turbulent simulations, we carry out chemodynamical simulations of stratified media in a Navarro-Frenk-White (NFW) gravitational potential with a total mass of $10^{12}$ solar masses and turbulence that decreases radially. The simulations assume a metallicity of 0.3 solar, a redshift zero metagalatic UV background, and they track ionizations, recombinations, and species-by-species radiative cooling using the MAIHEM package. We compare a suite of ionic column densities with the COS-Halos sample of low-redshift star-forming galaxies. Turbulence with an average one-dimensional velocity dispersion approximately 40 km/s, corresponding to an energy injection rate of approximately $10^{49}$ erg/yr, produces a CGM that matches many of the observed ionic column densities and ratios. In this simulation, the NVI to OVI ratio is suppressed from its equilibrium value due to a combination of radiative cooling and cooling from turbulent mixing. This level of turbulence is consistent with expectations from observations of better constrained, higher-mass systems, and could be sustained by energy input from supernovae, gas inflows, and dynamical friction from dark matter subhalos. We also conduct a higher resolution run which yields smaller-scale structures, but remains in agreement with observations. △ Less

Submitted 4 June, 2020; originally announced June 2020.

Comments: 24 pages, 19 Figures, ApJ, in press

A New Model For Including Galactic Winds in Simulations of Galaxy Formation I: Introducing the Physically Evolved Winds (PhEW) Model

Authors: Shuiyao Huang, Neal Katz, Evan Scannapieco, J'Neil Cottle, Romeel Davé, David H. Weinberg, Molly S. Peeples, Marcus Brüggen

Abstract: The propagation and evolution of cold galactic winds in galactic haloes is crucial to galaxy formation models. However, modelling of this process in hydrodynamic simulations of galaxy formation is over-simplified owing to a lack of numerical resolution and often neglects critical physical processes such as hydrodynamic instabilities and thermal conduction. We propose an analytic model, Physically… ▽ More The propagation and evolution of cold galactic winds in galactic haloes is crucial to galaxy formation models. However, modelling of this process in hydrodynamic simulations of galaxy formation is over-simplified owing to a lack of numerical resolution and often neglects critical physical processes such as hydrodynamic instabilities and thermal conduction. We propose an analytic model, Physically Evolved Winds (PhEW), that calculates the evolution of individual clouds moving supersonically through a uniform ambient medium. Our model reproduces predictions from very high resolution cloud-crushing simulations that include isotropic thermal conduction over a wide range of physical conditions. We discuss the implementation of this model into cosmological hydrodynamic simulations of galaxy formation as a sub-grid prescription to model galactic winds more robustly both physically and numerically. △ Less

Submitted 22 February, 2022; v1 submitted 27 May, 2020; originally announced May 2020.

Comments: 21 pages, 9 figures, accepted for publication in MNRAS

The Launching of Cold Clouds by Galaxy Outflows III: The Influence of Magnetic Fields

Authors: J'Neil Cottle, Evan Scannapieco, Marcus Bruggen, Wladimir Banda-Barragan, Christoph Federrath

Abstract: Motivated by observations of outflowing galaxies, we investigate the combined impact of magnetic fields and radiative cooling on the evolution of cold clouds embedded in a hot wind. We perform a collection of three-dimensional adaptive mesh refinement, magnetohydrodynamical simulations that span two resolutions, and include fields that are aligned and transverse to the oncoming, super-Alfvénic mat… ▽ More Motivated by observations of outflowing galaxies, we investigate the combined impact of magnetic fields and radiative cooling on the evolution of cold clouds embedded in a hot wind. We perform a collection of three-dimensional adaptive mesh refinement, magnetohydrodynamical simulations that span two resolutions, and include fields that are aligned and transverse to the oncoming, super-Alfvénic material. Aligned fields have little impact on the overall lifetime of the clouds over the non-magnetized case, although they do increase the mixing between the wind and cloud material by a factor of $\approx 3.$ Transverse fields lead to magnetic draping, which isolates the clouds, but they also squeeze material in the direction perpendicular to the field lines, which leads to rapid mass loss. A resolution study suggests that the magnetized simulations have somewhat better convergence properties than non-magnetized simulations, and that a resolution of 64 zones per cloud radius is sufficient to accurately describe these interactions. We conclude that the combined effects of radiative cooling and magnetic fields are dependent on field orientation, but are unlikely to enhance cloud lifetimes beyond the effect of radiative cooling alone. △ Less

Submitted 18 February, 2020; originally announced February 2020.

Comments: 15 pages, 14 figures, accepted to ApJ

Interpreting Observations of Absorption Lines in the Circumgalactic Medium with a Turbulent Medium

Authors: Edward Buie II, Michele Fumagalli, Evan Scannapieco

Abstract: Single-phase photoionization equilibrium (PIE) models are often used to infer the underlying physical properties of galaxy halos probed in absorption with ions at different ionization potentials. To incorporate the effects of turbulence, we use the MAIHEM code to model an isotropic turbulent medium exposed to a redshift zero metagalactic UV background, while tracking the ionizations, recombination… ▽ More Single-phase photoionization equilibrium (PIE) models are often used to infer the underlying physical properties of galaxy halos probed in absorption with ions at different ionization potentials. To incorporate the effects of turbulence, we use the MAIHEM code to model an isotropic turbulent medium exposed to a redshift zero metagalactic UV background, while tracking the ionizations, recombinations, and species-by-species radiative cooling for a wide range of ions. By comparing observations and simulations over a wide range of turbulent velocities, densities, and metallicity with a Markov chain Monte Carlo technique, we find that MAIHEM models provide an equally good fit to the observed low-ionization species compared to PIE models, while reproducing at the same time high-ionization species such as \ion{Si}{4} and \ion{O}{6}. By including multiple phases, MAIHEM models favor a higher metallicity ($Z/Z_\odot \approx 40\%$) for the circumgalactic medium compared to PIE models. Furthermore, all of the solutions require some amount of turbulence ($σ_{\rm 3D} \geqslant 26\ {\rm km}\ {\rm s}^{-1}$). Correlations between turbulence, metallicity, column density, and impact parameter are discussed alongside mechanisms that drive turbulence within the halo. △ Less

Submitted 14 January, 2020; originally announced January 2020.

Comments: 19 pages, 13 figures, accepted in ApJ

Catastrophic Cooling in Superwinds: Line Emission and Non-equilibrium Ionization

Authors: William J Gray, M. S. Oey, Sergiy Silich, Evan Scannapieco

Abstract: Outflows are a pervasive feature of mechanical feedback from super star clusters (SSC) in starburst galaxies, playing a fundamental role in galaxy evolution. Observations are now starting to confirm that outflows can undergo catastrophic cooling, suppressing adiabatic superwinds. Here we present a suite of one-dimensional, hydrodynamic simulations that study the ionization structure of these outfl… ▽ More Outflows are a pervasive feature of mechanical feedback from super star clusters (SSC) in starburst galaxies, playing a fundamental role in galaxy evolution. Observations are now starting to confirm that outflows can undergo catastrophic cooling, suppressing adiabatic superwinds. Here we present a suite of one-dimensional, hydrodynamic simulations that study the ionization structure of these outflows and the resulting line emission generated by the cooling gas. We use the non-equilibrium atomic chemistry package within MAIHEM, our modified version of FLASH, which evolves the ionization state of the gas and computes the total cooling rate on an ion-by-ion basis. We find that catastrophically cooling models produce strong nebular line emission compared to adiabatic outflows. We also show that such models exhibit non-equilibrium conditions, thereby generating more highly ionized states than equivalent equilibrium models. When including photoionization from the parent SSC, catastrophically cooling models show strong C IV λ1549 and O VI λ1037 emission. For density bounded photoionization, He II λ1640, λ4686, C III] λ1908, Si IV λ1206, and Si III λ1400 are also strongly enhanced. These lines are seen in extreme starbursts where catastrophic cooling is likely to occur, suggesting that they may serve as diagnostics of such conditions. The higher ionization generated by these flows may help to explain line emission that cannot be attributed to SSC photoionization alone. △ Less

Submitted 28 October, 2019; originally announced October 2019.

Comments: 29 pages, 13 figures, accepted to ApJ

Magnetic helicity dissipation and production in an ideal MHD code

Authors: Axel Brandenburg, Evan Scannapieco

Abstract: We study a turbulent helical dynamo in a periodic domain by solving the ideal magnetohydrodynamic (MHD) equations with the FLASH code using the divergence-cleaning eight-wave method and compare our results with direct numerical simulations (DNS) using the Pencil Code. At low resolution, FLASH reproduces the DNS results qualitatively by developing the large-scale magnetic field expected from DNS, b… ▽ More We study a turbulent helical dynamo in a periodic domain by solving the ideal magnetohydrodynamic (MHD) equations with the FLASH code using the divergence-cleaning eight-wave method and compare our results with direct numerical simulations (DNS) using the Pencil Code. At low resolution, FLASH reproduces the DNS results qualitatively by developing the large-scale magnetic field expected from DNS, but at higher resolution, no large-scale magnetic field is obtained. In all those cases in which a large-scale magnetic field is generated, the ideal MHD results yield too little power at small scales. As a consequence, the small-scale current helicity is too small compared with that of the DNS. The resulting net current helicity has then always the wrong sign, and its statistical average also does not approach zero at late times, as expected from the DNS. Our results have implications for astrophysical dynamo simulations of stellar and galactic magnetism using ideal MHD codes. △ Less

Submitted 25 January, 2020; v1 submitted 14 October, 2019; originally announced October 2019.

Comments: 9 pages, 14 figures, 2 tables, ApJ 889, 55 (2020)

Report number: NORDITA-2019-091

Journal ref: Astrophys. J. 889, 55 (2020)

Warped diffusive radio halo around the quiescent spiral edge-on galaxy NGC 4565

Authors: V. Heesen, L. Whitler, P. Schmidt, A. Miskolczi, S. S. Sridhar, C. Horellou, R. Beck, G. Gürkan, E. Scannapieco, M. Brüggen, G. H. Heald, M. Krause, R. Paladino, B. Nikiel-Wroczyński, R. -J. Dettmar

Abstract: Cosmic rays play a pivotal role in launching galactic winds, particularly in quiescently star-forming galaxies where the hot gas alone is not sufficient to drive a wind. Except for the Milky Way, not much is known about the transport of cosmic rays in galaxies. In this Letter, we present low-frequency observations of the nearby edge-on spiral galaxy NGC 4565 using the LOw-Frequency ARray (LOFAR).… ▽ More Cosmic rays play a pivotal role in launching galactic winds, particularly in quiescently star-forming galaxies where the hot gas alone is not sufficient to drive a wind. Except for the Milky Way, not much is known about the transport of cosmic rays in galaxies. In this Letter, we present low-frequency observations of the nearby edge-on spiral galaxy NGC 4565 using the LOw-Frequency ARray (LOFAR). With our deep 144-MHz observations, we obtain a clean estimate of the emission originating from old cosmic-ray electrons (CRe), which is almost free from contamination by thermal emission. We measured vertical profiles of the non-thermal radio continuum emission that we fitted with Gaussian and exponential functions. The different profile shapes correspond to 1D cosmic-ray transport models of pure diffusion and advection, respectively. We detect a warp in the radio continuum that is reminiscent of the previously known HI warp. Because the warp is not seen at GHz-frequencies in the radio continuum, its minimum age must be about 100 Myr. The warp also explains the slight flaring of the thick radio disc that can otherwise be well described by a Gaussian profile with an FWHM of 65 arcsec (3.7 kpc). The diffusive radio halo together with the extra-planar X-ray emission may be remnants of enhanced star-forming activity in the past where the galaxy had a galactic wind, as GHz-observations indicate only a weak outflow in the last 40 Myr. NGC 4565 could be in transition from an outflow- to an inflow-dominated phase. △ Less

Submitted 16 July, 2019; originally announced July 2019.

Comments: 5 pages, 4 figures, 2 tables; A\&A Letters accepted; appendix with cosmic-ray transport models

Journal ref: A&A 628, L3 (2019)

Measuring the delay time distribution of binary neutron stars. II. Using the redshift distribution from third-generation gravitational wave detectors network

Authors: Mohammadtaher Safarzadeh, Edo Berger, Ken K. -Y. Ng, Hsin-Yu Chen, Salvatore Vitale, Chris Whittle, Evan Scannapieco

Abstract: We investigate the ability of current and third-generation gravitational wave (GW) detectors to determine the delay time distribution (DTD) of binary neutron stars (BNS) through a direct measurement of the BNS merger rate as a function of redshift. We assume that the DTD follows a power law distribution with a slope $Γ$ and a minimum merger time $t_{\rm min}$, and also allow the overall BNS format… ▽ More We investigate the ability of current and third-generation gravitational wave (GW) detectors to determine the delay time distribution (DTD) of binary neutron stars (BNS) through a direct measurement of the BNS merger rate as a function of redshift. We assume that the DTD follows a power law distribution with a slope $Γ$ and a minimum merger time $t_{\rm min}$, and also allow the overall BNS formation efficiency per unit stellar mass to vary. By convolving the DTD and mass efficiency with the cosmic star formation history, and then with the GW detector capabilities, we explore two relevant regimes. First, for the current generation of GW detectors, which are only sensitive to the local universe, but can lead to precise redshift determinations via the identification of electromagnetic counterparts and host galaxies, we show that the DTD parameters are strongly degenerate with the unknown mass efficiency and therefore cannot be determined uniquely. Second, for third-generation detectors such as Einstein Telescope (ET) and Cosmic Explorer (CE), which will detect BNS mergers at cosmological distances, but with a redshift uncertainty inherent to GW-only detections ($δ(z)/z\approx 0.1z$), we show that the DTD and mass efficiency can be well-constrained to better than 10\% with a year of observations. This long-term approach to determining the DTD through a direct mapping of the BNS merger redshift distribution will be supplemented by more near term studies of the DTD through the properties of BNS merger host galaxies at $z\approx 0$ (Safarzadeh & Berger 2019). △ Less

Submitted 20 May, 2019; v1 submitted 24 April, 2019; originally announced April 2019.

Comments: 10 pages, Accepted to ApJ Letters

Non-equilibrium Ionization States Within Galactic Outflows: Explaining Their O VI and N V Column Densities

Authors: William J Gray, Evan Scannapieco, Matthew D. Lehnert

Abstract: We present a suite of one-dimensional spherically-symmetric hydrodynamic simulations that study the atomic ionization structure of galactic outflows. We track the ionization state of the outflowing gas with a non-equilibrium atomic chemistry network that includes photoionization, photo-heating, and ion-by-ion cooling. Each simulation describes a steady-state outflow that is defined by its mass and… ▽ More We present a suite of one-dimensional spherically-symmetric hydrodynamic simulations that study the atomic ionization structure of galactic outflows. We track the ionization state of the outflowing gas with a non-equilibrium atomic chemistry network that includes photoionization, photo-heating, and ion-by-ion cooling. Each simulation describes a steady-state outflow that is defined by its mass and energy input rates, sonic radius, metallicity, and UV flux from both the host galaxy and meta-galactic background. We find that for a large range of parameter choices, the ionization state of the material departs strongly from what it would be in photo-ionization equilibrium, in conflict with what is commonly assumed in the analysis of observations. In addition, nearly all the models reproduce the low N V to O VI column density ratios and the relatively high O VI column densities that are observed. △ Less

Submitted 14 March, 2019; originally announced March 2019.

Comments: 16 pages, 11 figures, accepted to ApJ. Comments welcome

Following the Cosmic Evolution of Pristine Gas III: The Observational Consequences of the Unknown Properties of Population III Stars

Authors: Richard Sarmento, Evan Scannapieco, Benoit Côté

Abstract: We study the observational consequences of several unknown properties of Population III (Pop III) stars using large-scale cosmological simulations that include a subgrid model to track the unresolved mixing of pollutants. Varying the value of the critical metallicity that marks the boundary between Pop III and Population II (Pop II) star formation across 2 dex has a negligible effect on the fracti… ▽ More We study the observational consequences of several unknown properties of Population III (Pop III) stars using large-scale cosmological simulations that include a subgrid model to track the unresolved mixing of pollutants. Varying the value of the critical metallicity that marks the boundary between Pop III and Population II (Pop II) star formation across 2 dex has a negligible effect on the fraction of Pop III stars formed and the subsequent fraction of Pop III flux from high-redshift galaxies. However, adopting a log normal initial mass function (IMF) for Pop III stars, in place of a baseline Salpeter IMF, results in a Pop III star formation rate density (SFRD) that is 1/4 of the baseline rate. The flux from high-redshift galaxies modeled with this IMF is highly bimodal, resulting in a tiny fraction of $z \leq 8$ galaxies with more than 75\% of their flux coming from Pop III stars. However, at $z=9$, right before reionization in our simulations, $\approx$ 20\% of galaxies are Pop III-bright with $m_{\rm UV} \le 31.4$ mag and at least 75\% of their flux generated by Pop III stars . Additionally, the log normal Pop III IMF results in a population of carbon enhanced, metal poor stars in reasonable agreement with MW halo observations. Our analysis supports the conclusion that the Pop III IMF was dominated by stars in the 20-120$M_{\odot}$ range that generate SN with carbon-enhanced ejecta. △ Less

Submitted 11 January, 2019; originally announced January 2019.

Comments: Accepted by ApJ

On Neutron Star Mergers as the Source of r-process Enhanced Metal Poor Stars in the Milky Way

Authors: Mohammadtaher Safarzadeh, Richard Sarmento, Evan Scannapieco

Abstract: We model the history of Galactic r-process enrichment using high-redshift, high-resolution zoom cosmological simulations of a Milky Way (MW) type halo. We assume that all r-process sources are neutron star mergers (NSMs) with a power law delay time distribution. We model the time to mix pollutants at subgrid scales, which allows us to to better compute the properties of metal poor (MP) and carbon… ▽ More We model the history of Galactic r-process enrichment using high-redshift, high-resolution zoom cosmological simulations of a Milky Way (MW) type halo. We assume that all r-process sources are neutron star mergers (NSMs) with a power law delay time distribution. We model the time to mix pollutants at subgrid scales, which allows us to to better compute the properties of metal poor (MP) and carbon enhanced metal poor (CEMP) stars, along with statistics of their r-process enhanced subclasses. Our simulations underpredict the cumulative ratios of r-process enhanced MP and CEMP stars (MP-r, CEMP-r) over MP and CEMP stars by about one order of magnitude, even when the minimum coalescence time of the double neutron stars ($t_{\rm min}$) is set to 1 Myr. No r-process enhanced stars form if $t_{\rm min}=100$ Myr. Our results show that even when we adopt the r-process yield estimates observed in GW170817, NSMs by themselves can only explain the observed frequency of r-process enhanced stars if either the birth rate of double neutron stars per unit mass of stars is boosted to $\approx10^{-4} M_\odot^{-1}$ or the Europium yield of each NSM event is boosted to $\approx 10^{-4} M_{\odot}$. △ Less

Submitted 25 March, 2019; v1 submitted 6 December, 2018; originally announced December 2018.

Comments: 11 pages, Accepted for publication at ApJ

Calibrating the relation of low-frequency radio continuum to star formation rate at 1 kpc scale with LOFAR

Authors: V. Heesen, E. Buie II, C. J. Huff, L. A. Perez, J. G. Woolsey, D. A. Rafferty, A. Basu, R. Beck, E. Brinks, C. Horellou, E. Scannapieco, M. Brüggen, R. -J. Dettmar, K. Sendlinger, B. Nikiel-Wroczyński, K. T. Chyży, P. N. Best, G. H. Heald, R. Paladino

Abstract: Radio continuum (RC) emission in galaxies allows us to measure star formation rates (SFRs) unaffected by extinction due to dust, of which the low-frequency part is uncontaminated from thermal (free-free) emission. We calibrate the conversion from the spatially resolved 140 MHz RC emission to the SFR surface density ($Σ_{\rm SFR}$) at 1 kpc scale. We used recent observations of three galaxies (NGC… ▽ More Radio continuum (RC) emission in galaxies allows us to measure star formation rates (SFRs) unaffected by extinction due to dust, of which the low-frequency part is uncontaminated from thermal (free-free) emission. We calibrate the conversion from the spatially resolved 140 MHz RC emission to the SFR surface density ($Σ_{\rm SFR}$) at 1 kpc scale. We used recent observations of three galaxies (NGC 3184, 4736, and 5055) from the LOFAR Two-metre Sky Survey (LoTSS), and archival LOw-Frequency ARray (LOFAR) data of NGC 5194. Maps were created with the facet calibration technique and converted to radio $Σ_{\rm SFR}$ maps using the Condon relation. We compared these maps with hybrid $Σ_{\rm SFR}$ maps from a combination of GALEX far-ultraviolet and Spitzer 24 $μ\rm m$ data using plots tracing the relation at $1.2\times 1.2$-kpc$^2$ resolution. The RC emission is smoothed with respect to the hybrid $Σ_{\rm SFR}$ owing to the transport of cosmic-ray electrons (CREs). This results in a sublinear relation $(Σ_{\rm SFR})_{\rm RC} \propto [(Σ_{\rm SFR})_{\rm hyb}]^{a}$, where $a=0.59\pm 0.13$ (140 MHz) and $a=0.75\pm 0.10$ (1365 MHz). Both relations have a scatter of $σ= 0.3~\rm dex$. If we restrict ourselves to areas of young CREs ($α> -0.65$; $I_ν\propto ν^α$), the relation becomes almost linear at both frequencies with $a\approx 0.9$ and a reduced scatter of $σ= 0.2~\rm dex$. We then simulate the effect of CRE transport by convolving the hybrid $Σ_{\rm SFR}$ maps with a Gaussian kernel until the RC-SFR relation is linearised; CRE transport lengths are $l=1$-5 kpc. Solving the CRE diffusion equation, we find diffusion coefficients of $D=(0.13$-$1.5) \times 10^{28} \rm cm^2\,s^{-1}$ at 1 GeV. A RC-SFR relation at $1.4$ GHz can be exploited to measure SFRs at redshift $z \approx 10$ using $140$ MHz observations. △ Less

Submitted 19 November, 2018; originally announced November 2018.

Comments: 9 figures, 6 tables and 17 pages. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available on lofar-surveys.org upon publication of the journal

Hot Atmospheres, Cold Gas, AGN Feedback and the Evolution of Early Type Galaxies: a Topical Perspective

Authors: N. Werner, B. R. McNamara, E. Churazov, E. Scannapieco

Abstract: Most galaxies comparable to or larger than the mass of the Milky Way host hot, X-ray emitting atmospheres, and many such galaxies are radio sources. Hot atmospheres and radio jets and lobes are the ingredients of radio-mechanical active galactic nucleus (AGN) feedback. While a consensus has emerged that such feedback suppresses cooling of hot cluster atmospheres, less attention has been paid to ma… ▽ More Most galaxies comparable to or larger than the mass of the Milky Way host hot, X-ray emitting atmospheres, and many such galaxies are radio sources. Hot atmospheres and radio jets and lobes are the ingredients of radio-mechanical active galactic nucleus (AGN) feedback. While a consensus has emerged that such feedback suppresses cooling of hot cluster atmospheres, less attention has been paid to massive galaxies where similar mechanisms are at play. Observation indicates that the atmospheres of elliptical and S0 galaxies were accreted externally during the process of galaxy assembly and augmented significantly by stellar mass loss. Their atmospheres have entropy and cooling time profiles that are remarkably similar to those of central cluster galaxies. About half display filamentary or disky nebulae of cool and cold gas, much of which has likely cooled from the hot atmospheres. We review the observational and theoretical perspectives on thermal instabilities in galactic atmospheres and the evidence that AGN heating is able to roughly balance the atmospheric cooling. Such heating and cooling may be regulating star formation in all massive spheroids at late times. △ Less

Submitted 19 December, 2018; v1 submitted 12 November, 2018; originally announced November 2018.

Comments: Final version

Journal ref: Space Sci Rev (2019) 215: 5

Astrophysics with the Spatially and Spectrally Resolved Sunyaev-Zeldovich Effects: A Millimetre/Submillimetre Probe of the Warm and Hot Universe

Authors: Tony Mroczkowski, Daisuke Nagai, Kaustuv Basu, Jens Chluba, Jack Sayers, Rémi Adam, Eugene Churazov, Abigail Crites, Luca Di Mascolo, Dominique Eckert, Juan Macias-Perez, Frédéric Mayet, Laurence Perotto, Etienne Pointecouteau, Charles Romero, Florian Ruppin, Evan Scannapieco, John ZuHone

Abstract: In recent years, observations of the Sunyaev-Zeldovich (SZ) effect have had significant cosmological implications and have begun to serve as a powerful and independent probe of the warm and hot gas that pervades the Universe. As a few pioneering studies have already shown, SZ observations both complement X-ray observations -- the traditional tool for studying the intra-cluster medium -- and bring… ▽ More In recent years, observations of the Sunyaev-Zeldovich (SZ) effect have had significant cosmological implications and have begun to serve as a powerful and independent probe of the warm and hot gas that pervades the Universe. As a few pioneering studies have already shown, SZ observations both complement X-ray observations -- the traditional tool for studying the intra-cluster medium -- and bring unique capabilities for probing astrophysical processes at high redshifts and out to the low-density regions in the outskirts of galaxy clusters. Advances in SZ observations have largely been driven by developments in centimetre-, millimetre-, and submillimetre-wave instrumentation on ground-based facilities, with notable exceptions including results from the Planck satellite. Here we review the utility of the thermal, kinematic, relativistic, non-thermal, and polarised SZ effects for studies of galaxy clusters and other large scale structures, incorporating the many advances over the past two decades that have impacted SZ theory, simulations, and observations. We also discuss observational results, techniques, and challenges, and aim to give an overview and perspective on emerging opportunities, with the goal of highlighting some of the exciting new directions in this field. △ Less

Submitted 21 January, 2020; v1 submitted 6 November, 2018; originally announced November 2018.

Comments: This version corrects a few typos. 72 pages

Journal ref: Space Sci Rev (2019) 215:17

r-process enrichment of ultra-faint dwarf galaxies by fast merging double neutron stars

Authors: Mohammadtaher Safarzadeh, Enrico Ramirez-Ruiz, Jeff. J. Andrews, Tassos Fragos, Phillip Macias, Evan Scannapieco

Abstract: The recent aLIGO/aVirgo discovery of gravitational waves from the neutron star merger (NSM) GW170817 and the follow up kilonova observations have shown that NSMs produce copious amount of r-process material. However, it is difficult to reconcile the large natal kicks and long average merging times of Double Neutron Stars (DNSs), with the levels of r-process enrichment seen in ultra faint dwarf (UF… ▽ More The recent aLIGO/aVirgo discovery of gravitational waves from the neutron star merger (NSM) GW170817 and the follow up kilonova observations have shown that NSMs produce copious amount of r-process material. However, it is difficult to reconcile the large natal kicks and long average merging times of Double Neutron Stars (DNSs), with the levels of r-process enrichment seen in ultra faint dwarf (UFD) galaxies such as Reticulum II and Tucana III. Assuming that such dwarf systems have lost a significant fraction of their stellar mass through tidal stripping, we conclude that contrary to most current models, it is the DNSs with rather large natal kicks but very short merging timescales that can enrich UFD-type galaxies. These binaries are either on highly eccentric orbits, or form with very short separations due to an additional mass-transfer between the first-born neutron star and a naked helium star, progenitor of the second-born neutron star. These DNSs are born with a frequency that agrees with the statistics of the r-process UFDs, and merge well within the virial radius of their host halos, therefore contributing significantly to their r-process enrichment. △ Less

Submitted 12 January, 2019; v1 submitted 9 October, 2018; originally announced October 2018.

Comments: Accepted for publication in ApJ

Modeling Star Formation as a Markov Process in a Supersonic Gravoturbulent Medium

Authors: Evan Scannapieco, Mohammadtaher Safarzadeh

Abstract: Molecular clouds exhibit lognormal probability density functions (PDF) of mass densities, which are thought to arise as a consequence of isothermal, supersonic turbulence. Star formation is then widely assumed to occur in perturbations in which gravitational collapse is faster than the rate of change due to turbulent motions. Here we use direct numerical simulations to measure this rate as a funct… ▽ More Molecular clouds exhibit lognormal probability density functions (PDF) of mass densities, which are thought to arise as a consequence of isothermal, supersonic turbulence. Star formation is then widely assumed to occur in perturbations in which gravitational collapse is faster than the rate of change due to turbulent motions. Here we use direct numerical simulations to measure this rate as a function of density for a range of turbulent Mach numbers, and show it is faster at high densities than at low densities. Furthermore, we show that both the density PDF and rate of change arise naturally in a simple model of turbulence as a continuous Markov process. The one-dimensional Langevin equation that describes this evolution depends on only two parameters, yet it captures the full evolution seen in direct three-dimensional simulations. If it is modified to include gravity, the Langevin equation also reproduces the rate of material collapsing to high densities seen in turbulent simulations including self-gravity. When generalized to include both temperature and density, similar analyses are likely applicable throughout astrophysics. △ Less

Submitted 20 September, 2018; originally announced September 2018.

Comments: 6 Pages, ApJ Letters, in Press

Modeling Photoionized Turbulent Material in the Circumgalactic Medium

Authors: Edward Buie II, William Gray, Evan Scannapieco

Abstract: The circumgalactic medium (CGM) of nearby star-forming galaxies show clear indications of \ion{O}{6} absorption accompanied by little to no \ion{N}{5} absorption. This unusual spectral signature, accompanied by absorption from lower ionization state species whose columns vary by orders of magnitude along \st{difference} \textbf{different} sightlines, indicates that the CGM must be viewed as a dyna… ▽ More The circumgalactic medium (CGM) of nearby star-forming galaxies show clear indications of \ion{O}{6} absorption accompanied by little to no \ion{N}{5} absorption. This unusual spectral signature, accompanied by absorption from lower ionization state species whose columns vary by orders of magnitude along \st{difference} \textbf{different} sightlines, indicates that the CGM must be viewed as a dynamic, multiphase medium, such as occurs in the presence of turbulence. To explore this possibility, we carry out a series of chemodynamical simulations of a isotropic turbulent media, using the MAIHEM package. The simulations assume a metallicity of $0.3\ Z_{\odot}$ and a redshift zero metagalatic UV background, and they track ionizations, recombinations, and species-by-species radiative cooling for a wide range of elements. We find that turbulence with a one-dimensional velocity dispersion of $σ_{1D} \approx 60$ km/s replicates many of the observed features within the CGM, such as clumping of low ionization-state ions and the existence of \ion{O}{6} at moderate ionization parameters. However, unlike observations, \ion{N}{5} often arises in our simulations with derived column densities of a similar magnitude to those of \ion{O}{6}. While higher values of $σ_{1D}$ lead to a thermal runaway in our isotropic simulations, this would not be the case in stratified media, and thus we speculate that more complex models of the turbulence may well match the absence of \ion{N}{5} in the CGM of star-forming galaxies. △ Less

Submitted 14 September, 2018; originally announced September 2018.

Comments: 10 pages, 8 figures, accepted to ApJ

Journal ref: Edward Buie II et al. 2018 ApJ 864 114

Column Density Profiles of Cold Clouds Driven by Galactic Outflows

Authors: J'Neil Cottle, Evan Scannapieco, Marcus Bruggen

Abstract: Absorption line studies are essential to understanding the origin, nature, and impact of starburst-driven galactic outflows. Such studies have revealed a multiphase medium with a number of poorly-understood features leading to a need to study the ionization mechanism of this gas. To better interpret these observations, we make use of a suite of adaptive mesh refinement hydrodynamic simulations of… ▽ More Absorption line studies are essential to understanding the origin, nature, and impact of starburst-driven galactic outflows. Such studies have revealed a multiphase medium with a number of poorly-understood features leading to a need to study the ionization mechanism of this gas. To better interpret these observations, we make use of a suite of adaptive mesh refinement hydrodynamic simulations of cold, atomic clouds driven by hot, supersonic outflows, including the effect of radiative cooling, thermal conduction, and an ionizing background characteristic of a starbursting galaxy. Using a new analysis tool, Trident, we estimate the equilibrium column density distributions for ten different ions: H I, Mg II, C II, C III, C IV, Si III, Si IV, N V, O VI, and Ne VIII. These are fit to model profiles with two parameters describing the maximum column density and coverage, and for each ion we provide a table of these fit parameters, along with average velocities and line widths. Our results are most sensitive to Mach number and conduction efficiency, with higher Mach numbers and more efficient conduction leading to more compact, high column density clouds. We use our results to interpret down-the-barrel observations of outflows and find that the adopted ionization equilibrium model overpredicts column densities of ions such as Si IV and does not adequately capture the observed trends for N V and O VI, implying the presence of strong non equilibrium ionization effects. △ Less

Submitted 24 July, 2018; originally announced July 2018.

Comments: 16 pages, 6 figures, accepted for publication in ApJ

A limit on the warm dark matter particle mass from the redshifted 21 cm absorption line

Authors: Mohammadtaher Safarzadeh, Evan Scannapieco, Arif Babul

Abstract: The recent EDGES collaboration detection of an absorption signal at a central frequency of $ν= 78 \pm 1$ MHz points to the presence of a significant Lyman-$α$ background by a redshift of $z=18$. The timing of this signal constrains the dark matter particle mass ($m_χ$) in the warm dark matter (WDM) cosmological model. WDM delays the formation of small-scale structures, and therefore a stringent lo… ▽ More The recent EDGES collaboration detection of an absorption signal at a central frequency of $ν= 78 \pm 1$ MHz points to the presence of a significant Lyman-$α$ background by a redshift of $z=18$. The timing of this signal constrains the dark matter particle mass ($m_χ$) in the warm dark matter (WDM) cosmological model. WDM delays the formation of small-scale structures, and therefore a stringent lower limit can be placed on $m_χ,$ based on the presence of a sufficiently strong Ly-$α$ background due to star formation at $z=18$. Our results show that the coupling the spin temperature to the gas through Ly-$α$ pumping requires a minimum mass of $m_χ>3$ keV if atomic cooling halos dominate the star formation rate at $z=18,$ and $m_χ>2$ keV if ${\rm\,H_2}$ cooling halos also form stars efficiently at this redshift. These limits match or exceed the most stringent limits cited to date in the literature, even in the face of the many uncertainties regarding star-formation at high redshift. △ Less

Submitted 16 May, 2018; v1 submitted 19 March, 2018; originally announced March 2018.

Comments: Accepted for publication in ApJ Letters

Selecting ultra-faint dwarf candidate progenitors in cosmological N-body simulations at high redshifts

Authors: Mohammadtaher Safarzadeh, Alexander P. Ji, Gregory A. Dooley, Anna Frebel, Evan Scannapieco, Facundo A. Gómez, Brian W. O'Shea

Abstract: The smallest satellites of the Milky Way ceased forming stars during the epoch of reionization and thus provide archaeological access to galaxy formation at $z>6$. Numerical studies of these ultra-faint dwarf galaxies (UFDs) require expensive cosmological simulations with high mass resolution that are carried out down to $z=0$. However, if we are able to statistically identify UFD host progenitors… ▽ More The smallest satellites of the Milky Way ceased forming stars during the epoch of reionization and thus provide archaeological access to galaxy formation at $z>6$. Numerical studies of these ultra-faint dwarf galaxies (UFDs) require expensive cosmological simulations with high mass resolution that are carried out down to $z=0$. However, if we are able to statistically identify UFD host progenitors at high redshifts \emph{with relatively high probabilities}, we can avoid this high computational cost. To find such candidates, we analyze the merger trees of Milky Way type halos from the high-resolution ${\it Caterpillar}$ suite of dark matter only simulations. Satellite UFD hosts at $z=0 $ are identified based on four different abundance matching (AM) techniques. All the halos at high redshifts are traced forward in time in order to compute the probability of surviving as satellite UFDs today. Our results show that selecting potential UFD progenitors based solely on their mass at z=12 (8) results in a 10\% (20\%) chance of obtaining a surviving UFD at $z=0$ in three of the AM techniques we adopted. We find that the progenitors of surviving satellite UFDs have lower virial ratios ($η$), and are preferentially located at large distances from the main MW progenitor, while they show no correlation with concentration parameter. Halos with favorable locations and virial ratios are $\approx 3$ times more likely to survive as satellite UFD candidates at $z=0.$ △ Less

Submitted 1 March, 2018; v1 submitted 11 December, 2017; originally announced December 2017.

Comments: 12 pages, 7 figures, accepted for publication at MNRAS after minor revisions

Using Real and Simulated Measurements of the Thermal Sunyaev-Zel'dovich Effect to Constrain Models of AGN Feedback

Authors: Alexander Spacek, Mark Richardson, Evan Scannapieco, Julien Devriendt, Yohan Dubois, Sebastien Peirani, Christophe Pichon

Abstract: Energetic feedback from active galactic nuclei (AGNs) is often used in simulations to resolve several outstanding issues in galaxy formation, but its impact is still not fully understood. Here we derive new constraints on AGN feedback by comparing observations and simulations of the thermal Sunyaev-Zel'dovich (tSZ) effect. We draw on observational results presented in Spacek et al. (2016, 2017) wh… ▽ More Energetic feedback from active galactic nuclei (AGNs) is often used in simulations to resolve several outstanding issues in galaxy formation, but its impact is still not fully understood. Here we derive new constraints on AGN feedback by comparing observations and simulations of the thermal Sunyaev-Zel'dovich (tSZ) effect. We draw on observational results presented in Spacek et al. (2016, 2017) who used data from the South Pole Telescope (SPT) and Atacama Cosmology Telescope (ACT) to measure the tSZ signal from >= 10^11 M_Sun and >= 1 Gyr galaxies at z=0.5-1.0 (low-z) and z=1.0-1.5 (high-z). Using the large-scale cosmological hydrodynamical simulations Horizon-AGN and Horizon-NoAGN, which include and omit AGN feedback, we extract simulated tSZ measurements around galaxies equivalent to the observational work. We find that the Horizon-AGN results only differ from the SPT measurements at levels of 0.4 sigma (low-z) and 0.6 sigma (high-z), but differ from the ACT measurements by 3.4 sigma (low-z) and 2.3 sigma (high-z). The Horizon-NoAGN results provide a slightly better fit to the SPT measurements by differing by 0.2 sigma (low-z) and 0.4 sigma (high-z), but a significantly better match to the ACT measurements by differing by only 0.5 sigma (low-z) and 1.4 sigma (high-z). We conclude that, while the lower-mass (<~ 5 x 10^11 M_Sun) SPT results allow for the presence AGN feedback energy, the higher-mass (>~ 5 x 10^11 M_Sun) ACT results show significantly less energy than predicted in the simulation including AGN feedback, while more closely matching the simulation without AGN feedback, indicating that AGN feedback may be milder than often predicted in simulations. △ Less

Submitted 11 August, 2018; v1 submitted 14 November, 2017; originally announced November 2017.

Comments: 12 pages, 8 figures, accepted for publication in ApJ

Report number: LA-UR-17-30403

Following the Cosmic Evolution of Pristine Gas II: The search for Pop III-Bright Galaxies

Authors: Richard Sarmento, Evan Scannapieco, Seth Cohen

Abstract: Direct observational searches for Population III (Pop III) stars at high redshift are faced with the question of how to select the most promising targets for spectroscopic follow-up. To help answer this, we use a large-scale cosmological simulation, augmented with a new subgrid model that tracks the fraction of pristine gas, to follow the evolution of high-redshift galaxies and the Pop III stars t… ▽ More Direct observational searches for Population III (Pop III) stars at high redshift are faced with the question of how to select the most promising targets for spectroscopic follow-up. To help answer this, we use a large-scale cosmological simulation, augmented with a new subgrid model that tracks the fraction of pristine gas, to follow the evolution of high-redshift galaxies and the Pop III stars they contain. We generate rest-frame ultraviolet (UV) luminosity functions for our galaxies and find that they are consistent with current $z \ge 7 $ observations. Throughout the redshift range $7 \le z \le 15$ we identify "Pop III-bright" galaxies as those with at least 75% of their flux coming from Pop III stars. While less than 1% of galaxies brighter than $m_{\rm UV, AB} = 31.4$ mag are Pop III--bright in the range $7\leq z \leq8$, roughly 17% of such galaxies are Pop III--bright at $z=9$, immediately before reionization occurs in our simulation. Moving to $z=10$, $m_{\rm UV, AB} = 31.4$ mag corresponds to larger, more luminous galaxies and the Pop III-bright fraction falls off to 5%. Finally, at the highest redshifts, a large fraction (29% at $z=14$ and 41% at $z=15)$ of all galaxies are Pop III-bright regardless of magnitude. While $m_{\rm UV, AB} = 31.4$ mag galaxies are extremely rare during this epoch, we find that 13% of galaxies at $z = 14$ are Pop III-bright with $m_{\rm UV, AB} \le 33$ mag, an intrisic magnitude within reach of the James Webb Space Telescope using lensing. Thus, we predict that the best redshift to search for luminous Pop III--bright galaxies is just before reionization, while lensing surveys for fainter galaxies should push to the highest redshifts possible. △ Less

Submitted 21 February, 2018; v1 submitted 26 October, 2017; originally announced October 2017.

Comments: Accepted for publication by ApJ, Updated version with correct Figure 6

The Fate of Gas rich Satellites in Clusters

Authors: Mohammadtaher Safarzadeh, Evan Scannapieco

Abstract: We investigate the stellar mass loss of gas rich galaxies falling into clusters due to the change in the gravitational potential caused by the ram pressure stripping of their gas. We model the satellites with exponential stellar and gas disk profiles, assume rapid ram pressure stripping, and follow the stellar orbits in the shocked potential. Due to the change of the potential, the stars move from… ▽ More We investigate the stellar mass loss of gas rich galaxies falling into clusters due to the change in the gravitational potential caused by the ram pressure stripping of their gas. We model the satellites with exponential stellar and gas disk profiles, assume rapid ram pressure stripping, and follow the stellar orbits in the shocked potential. Due to the change of the potential, the stars move from circular orbits to elliptical orbits with apocenters that are often outside the tidal radius, causing those stars to be stripped. We explore the impact of the redshift of infall, gas fraction, satellite halo mass and cluster mass on this process. The puffing of the satellites makes them appear as ultra diffuse galaxies, and the stripped stars contribute to the intracluster light. Our results show these effects are most significant for less-massive satellites, which have larger gas fractions when they are accreted into clusters. The preferential destruction of low mass systems causes the red fraction of cluster galaxies to be smaller at lower masses, an observation that is otherwise difficult to explain. △ Less

Submitted 17 October, 2017; v1 submitted 3 October, 2017; originally announced October 2017.

Comments: 9 pages, 8 figures, Accepted for publication in ApJ

The Effect of Turbulence on Nebular Emission Line Ratios

Authors: William J. Gray, Evan Scannapieco

Abstract: Motivated by the observed differences in the nebular emission of nearby and high-redshift galaxies, we carry out a set of direct numerical simulations of turbulent astrophysical media exposed to a UV background. The simulations assume a metallicity of $Z/Z_{\odot}$=0.5 and explicitly track ionization, recombination, charge transfer, and ion-by-ion radiative cooling for several astrophysically impo… ▽ More Motivated by the observed differences in the nebular emission of nearby and high-redshift galaxies, we carry out a set of direct numerical simulations of turbulent astrophysical media exposed to a UV background. The simulations assume a metallicity of $Z/Z_{\odot}$=0.5 and explicitly track ionization, recombination, charge transfer, and ion-by-ion radiative cooling for several astrophysically important elements. Each model is run to a global steady state that depends on the ionization parameter $U$, and the one-dimensional turbulent velocity dispersion, $σ_{\rm 1D}$, and the turbulent driving scale. We carry out a suite of models with a T=42,000K blackbody spectrum, $n_e$ = 100 cm$^{-3}$ and $σ_{\rm 1D}$ ranging between 0.7 to 42 km s$^{-1},$ corresponding to turbulent Mach numbers varying between 0.05 and 2.6. We report our results as several nebular diagnostic diagrams and compare them to observations of star-forming galaxies at a redshift of $z\approx$2.5, whose higher surface densities may also lead to more turbulent interstellar media. We find that subsonic, transsonic turbulence, and turbulence driven on scales of 1 parsec or greater, have little or no effect on the line ratios. Supersonic, small-scale turbulence, on the other hand, generally increases the computed line emission. In fact with a driving scale $\approx 0.1$ pc, a moderate amount of turbulence, $σ_{\rm 1D}$=21-28 km s$^{-1},$ can reproduce many of the differences between high and low redshift observations without resorting to harder spectral shapes. △ Less

Submitted 3 October, 2017; originally announced October 2017.

Comments: 18 pages, 9 figures, accepted to ApJ, comments welcome

Simulating Neutron Star Mergers as r-process Sources in Ultra Faint Dwarf Galaxies

Authors: Mohammadtaher Safarzadeh, Evan Scannapieco

Abstract: To explain the high observed abundances of r-process elements in local ultra- faint dwarf (UFD) galaxies, we perform cosmological zoom simulations that include r-process production from neutron star mergers (NSMs). We model star-formation stochastically and simulate two different halos with total masses $\approx 10^8 M_{\odot}$ at z = 6. We find that the final distribution of [Eu/H] vs. [Fe/H] is… ▽ More To explain the high observed abundances of r-process elements in local ultra- faint dwarf (UFD) galaxies, we perform cosmological zoom simulations that include r-process production from neutron star mergers (NSMs). We model star-formation stochastically and simulate two different halos with total masses $\approx 10^8 M_{\odot}$ at z = 6. We find that the final distribution of [Eu/H] vs. [Fe/H] is relatively insensitive to the energy by which the r-process material is ejected into the interstellar medium, but strongly sensitive to the environment in which the NSM event occurs. In one halo the NSM event takes place at the center of the stellar distribution, leading to high-levels of r-process enrichment such as seen in a local UFD, Reticulum II (Ret II). In a second halo, the NSM event takes place outside of the densest part of the galaxy, leading to a more extended r-process distribution. The subsequent star formation occurs in an interstellar medium with shallow levels of r-process enrichment which results in stars with low levels of [Eu/H] compared to Ret II stars even when the maximum possible r-process mass is assumed to be ejected. This suggests that the natal kicks of neutron stars may also play an important role in determining the r-process abundances in UFD galaxies, a topic that warrants further theoretical investigation. △ Less

Submitted 6 July, 2017; originally announced July 2017.

Comments: Accepted for publication in MNRAS,10 pages, 6 Figures

The Production of Cold Gas Within Galaxy Outflows

Authors: Evan Scannapieco

Abstract: I present a suite of three-dimensional simulations of the evolution of initially-hot material ejected by starburst-driven galaxy outflows. The simulations are conducted in a comoving frame that moves with the material, tracking atomic/ionic cooling, Compton cooling, and dust cooling and destruction. Compton cooling is most efficient of these processes, while the main role of atomic/ionic cooling i… ▽ More I present a suite of three-dimensional simulations of the evolution of initially-hot material ejected by starburst-driven galaxy outflows. The simulations are conducted in a comoving frame that moves with the material, tracking atomic/ionic cooling, Compton cooling, and dust cooling and destruction. Compton cooling is most efficient of these processes, while the main role of atomic/ionic cooling is to enhance density inhomogeneities. Dust, on the other hand, has little effect on the outflow evolution, and is rapidly destroyed in all the simulations except the case with the smallest mass flux. I use the results to construct a simple steady-state model of the observed UV/optical emission from each outflow. The velocity profiles in this case are dominated by geometric effects, and the overall luminosities are extremely strong functions of the properties of the host system, as observed in ultra-luminous infrared galaxies (ULIRGs). Furthermore the luminosities and maximum velocities in several models are consistent with emission-line observations of ULIRGs, although the velocities are significantly greater than observed in absorption-line studies. It may be that absorption line observations of galaxy outflows probe entrained cold material at small radii, while emission-line observations probe cold material condensing from the initially hot medium at larger distances. △ Less

Submitted 2 February, 2017; originally announced February 2017.

Comments: 17 Pages, 10 Figures, ApJ, in Press

The Impact of Unresolved Turbulence on the Escape Fraction of Lyman Continuum Photons

Authors: Mohammadtaher Safarzadeh, Evan Scannapieco

Abstract: We investigate the relation between the turbulent Mach number (\mach) and the escape fraction of Lyman continuum photons ($f_{\rm esc}$) in high-redshift galaxies. Approximating the turbulence as isothermal and isotropic, we show that the increase in the variance in column densities from $\mathcal{M}=1$ to $\mathcal{M}=10$ causes $f_{\rm esc}$ to increase by $\approx 25$\%, and the increase from… ▽ More We investigate the relation between the turbulent Mach number (\mach) and the escape fraction of Lyman continuum photons ($f_{\rm esc}$) in high-redshift galaxies. Approximating the turbulence as isothermal and isotropic, we show that the increase in the variance in column densities from $\mathcal{M}=1$ to $\mathcal{M}=10$ causes $f_{\rm esc}$ to increase by $\approx 25$\%, and the increase from $\mathcal{M}=1$ to $\mathcal{M}=20$ causes $f_{\rm esc}$ to increases by $\approx 50$\% for a medium with opacity $τ\approx1$. At a fixed Mach number, the correction factor for escape fraction relative to a constant column density case scales exponentially with the opacity in the cell, which has a large impact for simulated star forming regions. Furthermore, in simulations of isotropic turbulence with full atomic/ionic cooling and chemistry, the fraction of HI drops by a factor of $\approx 2.5$ at $\mathcal{M}\approx10$ even when the mean temperature is $\approx5\times10^3 K$. If turbulence is unresolved, these effects together enhance $f_{\rm esc}$ by a factor $>3$ at Mach numbers above 10. Such Mach numbers are common at high-redshifts where vigorous turbulence is driven by supernovae, gravitational instabilities, and merger activity, as shown both by numerical simulations and observations. These results, if implemented in the current hydrodynamical cosmological simulations to account for unresolved turbulence, can boost the theoretical predictions of the Lyman Continuum photon escape fraction and further constrain the sources of reionization. △ Less

Submitted 1 November, 2016; originally announced November 2016.

Comments: 5 pages, 3 figures, accepted for publication in ApJ Letters

Following The Cosmic Evolution Of Pristine Gas I: Implications For Milky Way Halo Stars

Authors: Richard Sarmento, Evan Scannapieco, Liubin Pan

Abstract: We make use of new subgrid model of turbulent mixing to accurately follow the cosmological evolution of the first stars, the mixing of their supernova ejecta, and the impact on the chemical composition of the Galactic Halo. Using the cosmological adaptive mesh refinement code RAMSES, we implement a model for the pollution of pristine gas as described in Pan et al. Tracking the metallicity of Pop I… ▽ More We make use of new subgrid model of turbulent mixing to accurately follow the cosmological evolution of the first stars, the mixing of their supernova ejecta, and the impact on the chemical composition of the Galactic Halo. Using the cosmological adaptive mesh refinement code RAMSES, we implement a model for the pollution of pristine gas as described in Pan et al. Tracking the metallicity of Pop III stars with metallicities below a critical value allows us to account for the fraction of Z < Zcrit stars formed even in regions in which the gas' average metallicity is well above Zcrit. We demonstrate that such partially-mixed regions account for 0.5 to 0.7 of all Pop III stars formed up to z = 5. Additionally, we track the creation and transport of "primordial metals" (PM) generated by Pop III supernovae (SNe). These neutron-capture deficient metals are taken up by second-generation stars and likely lead to unique abundance signatures characteristic of carbon-enhanced, metal-poor (CEMP-no) stars. As an illustrative example, we associate primordial metals with abundance ratios used by Keller et al. to explain the source of metals in the star SMSS J031300.36-670839.3, finding good agreement with the observed [Fe/H], [C/H], [O/H], and [Mg/Ca] ratios in CEMP-no Milky Way halo stars. Similar future simulations will aid in further constraining the properties of Pop III stars using CEMP observations, as well as improve predictions of the spatial distribution of Pop III stars, as will be explored by the next generation of ground- and space-based telescopes. △ Less

Submitted 1 December, 2016; v1 submitted 31 October, 2016; originally announced November 2016.

Comments: 19 pages, 13 figures, accepted by ApJ Oct 2016

Numerical Simulation of Star Formation by the Bow Shock of the Centaurus A Jet

Authors: Carl L. Gardner, Jeremiah R. Jones, Evan Scannapieco, Rogier A. Windhorst

Abstract: Recent Hubble Space Telescope (HST) observations of the extragalactic radio source Centaurus A (Cen A) display a young stellar population around the southwest tip of the inner filament 8.5 kpc from the Cen A galactic center, with ages in the range of 1-3 Myr. Crockett et al. (2012) argue that the transverse bow shock of the Cen A jet triggered this star formation as it impacted dense molecular cor… ▽ More Recent Hubble Space Telescope (HST) observations of the extragalactic radio source Centaurus A (Cen A) display a young stellar population around the southwest tip of the inner filament 8.5 kpc from the Cen A galactic center, with ages in the range of 1-3 Myr. Crockett et al. (2012) argue that the transverse bow shock of the Cen A jet triggered this star formation as it impacted dense molecular cores of clouds in the filament. To test this hypothesis, we perform three-dimensional numerical simulations of induced star formation by the jet bow shock in the inner filament of Cen A, using a positivity preserving WENO method to solve the equations of gas dynamics with radiative cooling. We find that star clusters form inside a bow-shocked molecular cloud when the maximum initial density of the cloud is > 40 H2 molecules/cm^3. In a typical molecular cloud of mass 10^6 M_sun and diameter 200 pc, approximately 20 star clusters of mass 10^4 M_sun are formed, matching the HST images. △ Less

Submitted 9 December, 2016; v1 submitted 6 October, 2016; originally announced October 2016.

Comments: 16 pages, 9 figures

Searching For Fossil Evidence of AGN Feedback in WISE-Selected Stripe-82 Galaxies By Measuring the Thermal Sunyaev-Zel'dovich Effect With the Atacama Cosmology Telescope

Authors: Alexander Spacek, Evan Scannapieco, Seth Cohen, Bhavin Joshi, Philip Mauskopf

Abstract: We directly measure the thermal energy of the gas surrounding galaxies through the thermal Sunyaev-Zel'dovich (tSZ) effect. We perform a stacking analysis of microwave background images from the Atacama Cosmology Telescope, around 1179 massive quiescent elliptical galaxies at 0.5 <= z <= 1.0 ('low-z') and 3274 galaxies at 1.0 <= z <= 1.5 ('high-z'), selected using data from the Wide-Field Infrared… ▽ More We directly measure the thermal energy of the gas surrounding galaxies through the thermal Sunyaev-Zel'dovich (tSZ) effect. We perform a stacking analysis of microwave background images from the Atacama Cosmology Telescope, around 1179 massive quiescent elliptical galaxies at 0.5 <= z <= 1.0 ('low-z') and 3274 galaxies at 1.0 <= z <= 1.5 ('high-z'), selected using data from the Wide-Field Infrared Survey Explorer All-Sky Survey and the Sloan Digital Sky Survey (SDSS) within the SDSS Stripe-82 field. The gas surrounding these galaxies is expected to contain energy from past episodes of active galactic nucleus (AGN) feedback, and after using modeling to subtract undetected contaminants, we detect a tSZ signal at a significance of 0.9-sigma for our low-z galaxies and 1.8-sigma for our high-z galaxies. We then include data from the high-frequency Planck bands for a subset of 227 low-z galaxies and 529 high-z galaxies and find low-z and high-z tSZ detections of 1.0-sigma and 1.5-sigma, respectively. These results indicate an average thermal heating around these galaxies of 5.6(+5.9/-5.6) x 10^60 erg for our low-z galaxies and 7.0(+4.7/-4.4) x 10^60 erg for our high-z galaxies. Based on simple heating models, these results are consistent with gravitational heating without additional heating due to AGN feedback. △ Less

Submitted 8 December, 2016; v1 submitted 6 October, 2016; originally announced October 2016.

Comments: 17 pages, 9 figures, accepted for publication in ApJ. Updated to reflect the final changes made for the published paper. arXiv admin note: substantial text overlap with arXiv:1601.01330

Journal ref: 2017, ApJ 834 102

On the Formation of Molecular Clumps in QSO Outflows

Authors: Andrea Ferrara, Evan Scannapieco

Abstract: We study the origin of the cold molecular clumps in quasar outflows, recently detected in CO and HCN emission. We first describe the physical properties of such radiation-driven outflows and show that a transition from a momentum- to an energy-driven flow must occur at a radial distance of R ~ 0.25 kpc. During this transition, the shell of swept up material fragments due to Rayleigh-Taylor instabi… ▽ More We study the origin of the cold molecular clumps in quasar outflows, recently detected in CO and HCN emission. We first describe the physical properties of such radiation-driven outflows and show that a transition from a momentum- to an energy-driven flow must occur at a radial distance of R ~ 0.25 kpc. During this transition, the shell of swept up material fragments due to Rayleigh-Taylor instabilities, but these clumps contain little mass and are likely to be rapidly ablated by the hot gas in which they are immersed. We then explore an alternative scenario in which clumps form from thermal instabilities at R >~ 1 kpc, possibly containing enough dust to catalyze molecule formation. We investigate this processes with 3D two-fluid (gas+dust) numerical simulations of a kpc^3 patch of the outflow, including atomic and dust cooling, thermal conduction, dust sputtering, and photoionization from the QSO radiation field. In all cases, dust grains are rapidly destroyed in ~10,000 years; and while some cold clumps form at later times, they are present only as transient features, which disappear as cooling becomes more widespread. In fact, we only find a stable two-phase medium with dense clumps if we artificially enhance the QSO radiation field by a factor 100. This result, together with the complete destruction of dust grains, renders the interpretation of molecular outflows a very challenging problem. △ Less

Submitted 8 September, 2016; originally announced September 2016.

Comments: 17 pages, 12 figures, ApJ, In press

Comparing Simulations of AGN Feedback

Authors: Mark L. A. Richardson, Evan Scannapieco, Julien Devriendt, Adrianne Slyz, Robert J. Thacker, Yohan Dubois, James Wurster, Joseph Silk

Abstract: We perform adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) cosmological zoom simulations of a region around a forming galaxy cluster, comparing the ability of the methods to handle successively more complex baryonic physics. In the simplest, non-radiative case, the two methods are in good agreement with each other, but the SPH simulations generate central cores with slight… ▽ More We perform adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) cosmological zoom simulations of a region around a forming galaxy cluster, comparing the ability of the methods to handle successively more complex baryonic physics. In the simplest, non-radiative case, the two methods are in good agreement with each other, but the SPH simulations generate central cores with slightly lower entropies and virial shocks at slightly larger radii, consistent with what has been seen in previous studies. The inclusion of radiative cooling, star formation, and stellar feedback leads to much larger differences between the two methods. Most dramatically, at z=5, rapid cooling in the AMR case moves the accretion shock well within the virial radius, while this shock remains near the virial radius in the SPH case, due to excess heating, coupled with poorer capturing of the shock width. On the other hand, the addition of feedback from active galactic nuclei (AGN) to the simulations results in much better agreement between the methods. In this case both simulations display halo gas entropies of 100 keV cm^2, similar decrements in the star-formation rate, and a drop in the halo baryon content of roughly 30%. This is consistent with AGN growth being self-regulated, regardless of the numerical method. However, the simulations with AGN feedback continue to differ in aspects that are not self-regulated, such that in SPH a larger volume of gas is impacted by feedback, and the cluster still has a lower entropy central core. △ Less

Submitted 11 May, 2016; originally announced May 2016.

Comments: 22 pages, 20 figures, 3 tables, Accepted to ApJ, comments welcome

The Launching of Cold Clouds by Galaxy Outflows II: The Role of Thermal Conduction

Authors: Marcus Brüggen, Evan Scannapieco

Abstract: We explore the impact of electron thermal conduction on the evolution of radiatively-cooled cold clouds embedded in flows of hot and fast material, as occur in outflowing galaxies. Performing a parameter study of three-dimensional adaptive mesh refinement hydrodynamical simulations, we show that electron thermal conduction causes cold clouds to evaporate, but it can also extend their lifetimes by… ▽ More We explore the impact of electron thermal conduction on the evolution of radiatively-cooled cold clouds embedded in flows of hot and fast material, as occur in outflowing galaxies. Performing a parameter study of three-dimensional adaptive mesh refinement hydrodynamical simulations, we show that electron thermal conduction causes cold clouds to evaporate, but it can also extend their lifetimes by compressing them into dense filaments. We distinguish between low column-density clouds, which are disrupted on very short times, and high-column density clouds with much-longer disruption times that are set by a balance between impinging thermal energy and evaporation. We provide fits to the cloud lifetimes and velocities that can be used in galaxy-scale simulations of outflows, in which the evolution of individual clouds cannot be modeled with the required resolution. Moreover, we show that the clouds are only accelerated to a small fraction of the ambient velocity because compression by evaporation causes the clouds to present a small cross-section to the ambient flow. This means that either magnetic fields must suppress thermal conduction, or that the cold clouds observed in galaxy outflows are not formed of cold material carried out from the galaxy. △ Less

Submitted 4 February, 2016; originally announced February 2016.

Comments: accepted by ApJ

Constraining AGN Feedback in Massive Ellipticals with South Pole Telescope Measurements of the Thermal Sunyaev-Zel'dovich Effect

Authors: Alexander Spacek, Evan Scannapieco, Seth Cohen, Bhavin Joshi, Philip Mauskopf

Abstract: Energetic feedback due to active galactic nuclei (AGN) is likely to play an important role in the observed anti-hierarchical trend in the evolution of galaxies, and yet the energy injected into the circumgalactic medium by this process is largely unknown. One promising approach to constrain this feedback is through measurements of CMB spectral distortions due to the thermal Sunyaev-Zel'dovich (tSZ… ▽ More Energetic feedback due to active galactic nuclei (AGN) is likely to play an important role in the observed anti-hierarchical trend in the evolution of galaxies, and yet the energy injected into the circumgalactic medium by this process is largely unknown. One promising approach to constrain this feedback is through measurements of CMB spectral distortions due to the thermal Sunyaev-Zel'dovich (tSZ) effect, whose magnitude is directly proportional to the energy input by AGN. Here we co-add South Pole Telescope SZ (SPT-SZ) survey data around a large set of massive quiescent elliptical galaxies at z >= 0.5. We use data from the Blanco Cosmology Survey and VISTA Hemisphere Survey to create a large catalog of galaxies split up into two redshift bins, with 3394 galaxies at 0.5 <= z <= 1.0 and 924 galaxies at 1.0 <= z <= 1.5, with typical stellar masses of 1.5 x 10^11 M_Sun. We then co-add the emission around these galaxies, resulting in a measured tSZ signal at 2.2 sigma significance for the lower redshift bin and a contaminating signal at 1.1 sigma for the higher redshift bin. To remove contamination due to dust emission, we use SPT-SZ source counts to model a contaminant source population in both the SPT-SZ bands and Planck high-frequency bands for a subset of 937 low-redshift galaxies and 240 high-redshift galaxies. This increases our detection to 3.6 sigma for low redshifts and 0.9 sigma for high redshifts. We find the mean angularly-integrated Compton-y values to be 2.2 (-0.7+0.9) x 10^-7 Mpc^2 for low redshifts and 1.7 (-1.8+2.2) x 10^-7 Mpc^2 for high redshifts, corresponding to total thermal energies of 7.6 (-2.3+3.0) x 10^60 ergs and 6.0 (-6.3+7.7) x 10^60 ergs, respectively. These numbers are higher than expected from simple theoretical models that do not include AGN feedback, and serve as constraints that can be applied to current simulations of massive galaxy formation. (abridged) △ Less

Submitted 8 June, 2016; v1 submitted 6 January, 2016; originally announced January 2016.

Comments: 23 pages, 13 figures, published in ApJ. Updated to reflect minor changes made for the published paper, plus fixed feedback efficiency results which were off by a factor of 5. The overall results weren't significantly affected

Journal ref: 2016, ApJ 819 128

Galaxy Outflows Without Supernovae

Authors: Sharanya Sur, Evan Scannapieco, Eve C. Ostriker

Abstract: High surface density, rapidly star-forming galaxies are observed to have $\approx 50-100\,{\rm km\,s^{-1}}$ line-of-sight velocity dispersions, which are much higher than expected from supernova driving alone, but may arise from large-scale gravitational instabilities. Using three-dimensional simulations of local regions of the interstellar medium, we explore the impact of high velocity dispersion… ▽ More High surface density, rapidly star-forming galaxies are observed to have $\approx 50-100\,{\rm km\,s^{-1}}$ line-of-sight velocity dispersions, which are much higher than expected from supernova driving alone, but may arise from large-scale gravitational instabilities. Using three-dimensional simulations of local regions of the interstellar medium, we explore the impact of high velocity dispersions that arise from these disk instabilities. Parametrizing disks by their surface densities and epicyclic frequencies, we conduct a series of simulations that probe a broad range of conditions. Turbulence is driven purely horizontally and on large scales, neglecting any energy input from supernovae. We find that such motions lead to strong global outflows in the highly-compact disks that were common at high redshifts, but weak or negligible mass loss in the more diffuse disks that are prevalent today. Substantial outflows are generated if the one-dimensional horizontal velocity dispersion exceeds $\approx 35\,{\rm km\,s^{-1}},$ as occurs in the dense disks that have star formation rate densities above $\approx 0.1\,{\rm M}_\odot\,{\rm yr}^{-1}\,{\rm kpc}^{-2}.$ These outflows are triggered by a thermal runaway, arising from the inefficient cooling of hot material coupled with successive heating from turbulent driving. Thus, even in the absence of stellar feedback, a critical value of the star-formation rate density for outflow generation can arise due to a turbulent heating instability. This suggests that in strongly self-gravitating disks, outflows may be enhanced by, but need not caused by, energy input from supernovae. △ Less

Submitted 4 January, 2016; originally announced January 2016.

Comments: Accepted in ApJ

Report number: NSF-KITP-16-002

Atomic Chemistry in Turbulent Astrophysical Media II: Effect of the Redshift Zero Metagalactic Background

Authors: William J Gray, Evan Scannapieco

Abstract: We carry out direct numerical simulations of turbulent astrophysical media exposed to the redshift zero metagalactic background. The simulations assume solar composition and explicitly track ionizations, recombinations, and ion-by-ion radiative cooling for hydrogen, helium, carbon, nitrogen, oxygen, neon, sodium, magnesium, silicon, sulfur, calcium, and iron. Each run reaches a global steady state… ▽ More We carry out direct numerical simulations of turbulent astrophysical media exposed to the redshift zero metagalactic background. The simulations assume solar composition and explicitly track ionizations, recombinations, and ion-by-ion radiative cooling for hydrogen, helium, carbon, nitrogen, oxygen, neon, sodium, magnesium, silicon, sulfur, calcium, and iron. Each run reaches a global steady state that not only depends on the ionization parameter, $U,$ and mass-weighted average temperature, $T_{\rm MW},$ but also on the the one-dimensional turbulent velocity dispersion, \soned. We carry out runs that span a grid of models with $U$ ranging from 0 to 10$^{-1}$ and \soned\ ranging from 3.5 to 58 km s$^{-1}$, and we vary the product of the mean density and the driving scale of the turbulence, $nL,$ which determines the average temperature of the medium, from $nL =10^{16}$ to $nL =10^{20}$ cm$^{-2}$. The turbulent Mach numbers of our simulations vary from $M \approx 0.5$ for the lowest velocity dispersions cases to $M \approx 20$ for the largest velocity dispersion cases. When $M \lesssim1,$ turbulent effects are minimal, and the species abundances are reasonably described as those of a uniform photoionized medium at a fixed temperature. On the other hand, when $M \gtrsim 1,$ dynamical simulations such as the ones carried out here are required to accurately predict the species abundances. We gather our results into a set of tables, to allow future redshift zero studies of the intergalactic medium to account for turbulent effects. △ Less

Submitted 16 November, 2015; originally announced November 2015.

Comments: 22 pages, 12 figures. Accepted to ApJ. Comments welcome