-
Numerical Modeling of Galactic Superwinds with Time-evolving Stellar Feedback
Authors:
A. Danehkar,
M. S. Oey,
W. J. Gray
Abstract:
Mass-loss and radiation feedback from evolving massive stars produce galactic-scale superwinds, sometimes surrounded by pressure-driven bubbles. Using the time-dependent stellar population typically seen in star-forming regions, we conduct hydrodynamic simulations of a starburst-driven superwind model coupled with radiative efficiency rates to investigate the formation of radiative cooling superwi…
▽ More
Mass-loss and radiation feedback from evolving massive stars produce galactic-scale superwinds, sometimes surrounded by pressure-driven bubbles. Using the time-dependent stellar population typically seen in star-forming regions, we conduct hydrodynamic simulations of a starburst-driven superwind model coupled with radiative efficiency rates to investigate the formation of radiative cooling superwinds and bubbles. Our numerical simulations depict the parameter space where radiative cooling superwinds with or without bubbles occur. Moreover, we employ the physical properties and time-dependent ionization states to predict emission line profiles under the assumption of collisional ionization and non-equilibrium ionization caused by wind thermal feedback in addition to photoionization created by the radiation background. We see the dependence of non-equilibrium ionization structures on the time-evolving ionizing source, leading to a deviation from collisional ionization in radiative cooling wind regions over time.
△ Less
Submitted 25 June, 2024;
originally announced June 2024.
-
Hydrodynamic Simulations and Time-dependent Photoionization Modeling of Starburst-driven Superwinds
Authors:
A. Danehkar,
M. S. Oey,
W. J. Gray
Abstract:
Thermal energies deposited by OB stellar clusters in starburst galaxies lead to the formation of galactic superwinds. Multi-wavelength observations of starburst-driven superwinds pointed at complex thermal and ionization structures which cannot adequately be explained by simple adiabatic assumptions. In this study, we perform hydrodynamic simulations of a fluid model coupled to radiative cooling f…
▽ More
Thermal energies deposited by OB stellar clusters in starburst galaxies lead to the formation of galactic superwinds. Multi-wavelength observations of starburst-driven superwinds pointed at complex thermal and ionization structures which cannot adequately be explained by simple adiabatic assumptions. In this study, we perform hydrodynamic simulations of a fluid model coupled to radiative cooling functions, and generate time-dependent non-equilibrium photoionization models to predict physical conditions and ionization structures of superwinds using the MAIHEM atomic and cooling package built on the program FLASH. Time-dependent ionization states and physical conditions produced by our simulations are used to calculate the emission lines of superwinds for various parameters, which allow us to explore implications of non-equilibrium ionization for starburst regions with potential radiative cooling.
△ Less
Submitted 17 February, 2023;
originally announced February 2023.
-
Catastrophic Cooling in Superwinds. III. Non-equilibrium Photoionization
Authors:
A. Danehkar,
M. S. Oey,
W. J. Gray
Abstract:
Observations of some starburst-driven galactic superwinds suggest that strong radiative cooling could play a key role in the nature of feedback and the formation of stars and molecular gas in star-forming galaxies. These catastrophically cooling superwinds are not adequately described by adiabatic fluid models, but they can be reproduced by incorporating non-equilibrium radiative cooling functions…
▽ More
Observations of some starburst-driven galactic superwinds suggest that strong radiative cooling could play a key role in the nature of feedback and the formation of stars and molecular gas in star-forming galaxies. These catastrophically cooling superwinds are not adequately described by adiabatic fluid models, but they can be reproduced by incorporating non-equilibrium radiative cooling functions into the fluid model. In this work, we have employed the atomic and cooling module MAIHEM implemented in the framework of the FLASH hydrodynamics code to simulate the formation of radiatively cooling superwinds as well as their corresponding non-equilibrium ionization (NEI) states for various outflow parameters, gas metallicities, and ambient densities. We employ the photoionization program CLOUDY to predict radiation- and density-bounded photoionization for these radiatively cooling superwinds, and we predict UV and optical line emission. Our non-equilibrium photoionization models built with the NEI states demonstrate the enhancement of C IV, especially in metal-rich, catastrophically cooling outflows, and O VI in metal-poor ones.
△ Less
Submitted 24 August, 2022;
originally announced August 2022.
-
Non-equilibrium Photoionization and Hydrodynamic Simulations of Starburst-driven Outflows
Authors:
A. Danehkar,
M. S. Oey,
W. J. Gray
Abstract:
Starburst-driven galactic outflows in star-forming galaxies have been observed to contain complex thermal structures and emission line features that are difficult to explain by adiabatic fluid models and plasmas in photoionization equilibrium (PIE) and collisional ionization equilibrium (CIE). We previously performed hydrodynamic simulations of starburst-driven outflows using the MAIHEM module for…
▽ More
Starburst-driven galactic outflows in star-forming galaxies have been observed to contain complex thermal structures and emission line features that are difficult to explain by adiabatic fluid models and plasmas in photoionization equilibrium (PIE) and collisional ionization equilibrium (CIE). We previously performed hydrodynamic simulations of starburst-driven outflows using the MAIHEM module for non-equilibrium atomic chemistry and radiative cooling functions in the hydrodynamics code FLASH, and calculated emission lines in combined CIE and PIE conditions. In the present study, we consider time-dependent non-equilibrium ionization (NEI) states produced by our MAIHEM simulations. Through extensive CLOUDY calculations made with the NEI states from our hydrodynamic simulations, we predict the UV and optical emission line profiles for starburst-driven outflows in time-evolving non-equilibrium photoionization conditions. Our hydrodynamic results demonstrate applications of non-equilibrium radiative cooling for H II regions in starburst galaxies hosting cool outflows.
△ Less
Submitted 3 November, 2021;
originally announced November 2021.
-
Catastrophic Cooling in Superwinds. II. Exploring the Parameter Space
Authors:
Ashkbiz Danehkar,
M. S. Oey,
William J. Gray
Abstract:
Superwinds and superbubbles driven by mechanical feedback from super star clusters (SSCs) are common features in many star-forming galaxies. While the adiabatic fluid model can well describe the dynamics of superwinds, several observations of starburst galaxies revealed the presence of compact regions with suppressed superwinds and strongly radiative cooling, i.e., catastrophic cooling. In the pre…
▽ More
Superwinds and superbubbles driven by mechanical feedback from super star clusters (SSCs) are common features in many star-forming galaxies. While the adiabatic fluid model can well describe the dynamics of superwinds, several observations of starburst galaxies revealed the presence of compact regions with suppressed superwinds and strongly radiative cooling, i.e., catastrophic cooling. In the present study, we employ the non-equilibrium atomic chemistry and cooling package MAIHEM, built on the FLASH hydrodynamics code, to generate a grid of models investigating the dependence of cooling modes on the metallicity, SSC outflow parameters, and ambient density. While gas metallicity plays a substantial role, catastrophic cooling is more sensitive to high mass-loading and reduced kinetic heating efficiency. Our hydrodynamic simulations indicate that the presence of a hot superbubble does not necessarily imply an adiabatic outflow, and vice versa. Using CLOUDY photoionization models, we predict UV and optical line emission for both adiabatic and catastrophic cooling outflows, for radiation-bounded and partially density-bounded models. Although the line ratios predicted by our radiation-bounded models agree well with observations of star-forming galaxies, they do not provide diagnostics that unambiguously distinguish the parameter space of catastrophically cooling flows. Comparison with observations suggests the possibility of minor density bounding, non-equilibrium ionization, and/or observational bias toward the central outflow regions.
△ Less
Submitted 4 November, 2021; v1 submitted 21 June, 2021;
originally announced June 2021.
-
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.
-
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.
-
Atomic Modeling of Photoionization Fronts in Nitrogen Gas
Authors:
William J. Gray,
P. A. Keiter,
H. Lefevre,
C. R. Patterson,
J. S. Davis,
K. G. Powell,
C. C. Kuranz,
R. P. Drake
Abstract:
Photoionization fronts play a dominant role in many astrophysical environments, but remain difficult to achieve in a laboratory experiment. Recent papers have suggested that experiments using a nitrogen medium held at ten atmospheres of pressure that is irradiated by a source with a radiation temperature of T$_{\rm R}\sim$ 100 eV can produce viable photoionization fronts. We present a suite of one…
▽ More
Photoionization fronts play a dominant role in many astrophysical environments, but remain difficult to achieve in a laboratory experiment. Recent papers have suggested that experiments using a nitrogen medium held at ten atmospheres of pressure that is irradiated by a source with a radiation temperature of T$_{\rm R}\sim$ 100 eV can produce viable photoionization fronts. We present a suite of one-dimensional numerical simulations using the \helios\ multi-material radiation hydrodynamics code that models these conditions and the formation of a photoionization front. We study the effects of varying the atomic kinetics and radiative transfer model on the hydrodynamics and ionization state of the nitrogen gas, finding that more sophisticated physics, in particular a multi-angle long characteristic radiative transfer model and a collisional-radiative atomics model, dramatically changes the atomic kinetic evolution of the gas. A photoionization front is identified by computing the ratios between the photoionization rate, the electron impact ionization rate, and the total recombination rate. We find that due to the increased electron temperatures found using more advanced physics that photoionization fronts are likely to form in our nominal model. We report results of several parameter studies. In one of these, the nitrogen pressure is fixed at ten atmospheres and varies the source radiation temperature while another fixes the temperature at 100 eV and varied the nitrogen pressure. Lower nitrogen pressures increase the likelihood of generating a photoionization front while varying the peak source temperature has little effect.
△ Less
Submitted 18 April, 2019;
originally announced April 2019.
-
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.
-
Challenges and Techniques for Simulating Line Emission
Authors:
Karen P. Olsen,
Andrea Pallottini,
Aida Wofford,
Marios Chatzikos,
Mitchell Revalski,
Francisco Guzmán,
Gergö Popping,
Enrique Vázquez-Semadeni,
Georgios E. Magdis,
Mark L. A. Richardson,
Michaela Hirschmann,
William J. Gray
Abstract:
Modeling emission lines from the millimeter to the UV and producing synthetic spectra is crucial for a good understanding of observations, yet it is an art filled with hazards. This is the proceedings of "Walking the Line", a 3-day conference held in 2018 that brought together scientists working on different aspects of emission line simulations, in order to share knowledge and discuss the methodol…
▽ More
Modeling emission lines from the millimeter to the UV and producing synthetic spectra is crucial for a good understanding of observations, yet it is an art filled with hazards. This is the proceedings of "Walking the Line", a 3-day conference held in 2018 that brought together scientists working on different aspects of emission line simulations, in order to share knowledge and discuss the methodology. Emission lines across the spectrum from the millimeter to the UV were discussed, with most of the focus on the interstellar medium, but also some topics on the circumgalactic medium. The most important quality of a useful model is a good synergy with observations and experiments. Challenges in simulating line emission are identified, some of which are already being worked upon, and others that must be addressed in the future for models to agree with observations. Recent advances in several areas aiming at achieving that synergy are summarized here, from micro-physical to galactic and circum-galactic scale.
△ Less
Submitted 24 August, 2018;
originally announced August 2018.
-
Laboratory Photoionization Fronts in Nitrogen Gas: A Numerical Feasibility and Parameter Study
Authors:
William J Gray,
P. A. Keiter,
H. Lefevre,
C. R. Patterson,
J. S. Davis,
B. van Der Holst,
K. G. Powell,
R. P. Drake
Abstract:
Photoionization fronts play a dominant role in many astrophysical situations, but remain difficult to achieve in a laboratory experiment. We present the results from a computational parameter study evaluating the feasibility of the photoionization experiment presented in the design paper by Drake, R. P., Hazak, G., Keiter, P. A., Davis, J. S., Patterson, C. R., Frank, A., Blackman, E. G., & Busque…
▽ More
Photoionization fronts play a dominant role in many astrophysical situations, but remain difficult to achieve in a laboratory experiment. We present the results from a computational parameter study evaluating the feasibility of the photoionization experiment presented in the design paper by Drake, R. P., Hazak, G., Keiter, P. A., Davis, J. S., Patterson, C. R., Frank, A., Blackman, E. G., & Busquet, M. 2016, ApJ, 833, 249 in which a photoionization front is generated in a nitrogen medium . The nitrogen gas density and the Planckian radiation temperature of the x-ray source define each simulation. Simulations modeled experiments in which the x-ray flux is generated by a laser-heated gold foil, suitable for experiments using many kJ of laser energy, and experiments in which the flux is generated by a "z-pinch" device, which implodes a cylindrical shell of conducting wires. The models are run using CRASH, our block-adaptive-mesh code for multi-material radiation hydrodynamics. The radiative transfer model uses multi-group, flux-limited diffusion with thirty radiation groups. In addition, electron heat conduction is modeled using a single-group, flux-limited diffusion. In the theory, a photoionization front can exist only when the ratios of the electron recombination rate to the photoionization rate and the electron impact ionization rate to the recombination rate lie in certain ranges. These ratios are computed for several ionization states of nitrogen. Photoionization fronts are found to exist for laser driven models with moderate nitrogen densities ($\sim$10$^{21}$ cm$^{-3}$) and radiation temperatures above 90 eV. For "z-pinch" driven models, lower nitrogen densities are preferred ($<$10$^{21}$ cm$^{-3}$). We conclude that the proposed experiments are likely to generate photoionization fronts.
△ Less
Submitted 4 April, 2018;
originally announced April 2018.
-
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.
-
Effect of Angular Momentum Alignment and Strong Magnetic Fields on the Formation of Protostellar Disks
Authors:
William J. Gray,
Christopher F. McKee,
Richard I. Klein
Abstract:
Star forming molecular clouds are observed to be both highly magnetized and turbulent. Consequently the formation of protostellar disks is largely dependent on the complex interaction between gravity, magnetic fields, and turbulence. Studies of non-turbulent protostellar disk formation with realistic magnetic fields have shown that these fields are efficient in removing angular momentum from the f…
▽ More
Star forming molecular clouds are observed to be both highly magnetized and turbulent. Consequently the formation of protostellar disks is largely dependent on the complex interaction between gravity, magnetic fields, and turbulence. Studies of non-turbulent protostellar disk formation with realistic magnetic fields have shown that these fields are efficient in removing angular momentum from the forming disks, preventing their formation. However, once turbulence is included, disks can form in even highly magnetized clouds, although the precise mechanism remains uncertain. Here we present several high resolution simulations of turbulent, realistically magnetized, high-mass molecular clouds with both aligned and random turbulence to study the role that turbulence, misalignment, and magnetic fields have on the formation of protostellar disks. We find that when the turbulence is artificially aligned so that the angular momentum is parallel to the initial uniform field, no rotationally supported disks are formed, regardless of the initial turbulent energy. We conclude that turbulence and the associated misalignment between the angular momentum and the magnetic field are crucial in the formation of protostellar disks in the presence of realistic magnetic fields.
△ Less
Submitted 21 September, 2017; v1 submitted 15 September, 2017;
originally announced September 2017.
-
Shadows of our Former Companions: How the Single-Degenerate Binary Type Ia Supernova Scenario Affects Remnants
Authors:
William J Gray,
Cody Raskin,
J. Michael Owen
Abstract:
Here we present three-dimensional high resolution simulations of Type Ia supernova in the presence of a non-degenerate companion. We find that the presence of a nearby companion leaves a long-lived hole in the supernova ejecta. In particular, we aim to study the long term evolution of this hole as the supernova ejecta interacts with the surrounding interstellar medium. Using estimates for the x-ra…
▽ More
Here we present three-dimensional high resolution simulations of Type Ia supernova in the presence of a non-degenerate companion. We find that the presence of a nearby companion leaves a long-lived hole in the supernova ejecta. In particular, we aim to study the long term evolution of this hole as the supernova ejecta interacts with the surrounding interstellar medium. Using estimates for the x-ray emission, we find that the hole generated by the companion remains for many centuries after the interaction between the ejecta and the interstellar medium. We also show that the hole is discernible over a wide range of viewing angles and companion masses.
△ Less
Submitted 20 September, 2016;
originally announced September 2016.
-
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.
-
Atomic Chemistry In Turbulent Astrophysical Media I: Effect of Atomic Cooling
Authors:
William J Gray,
Evan Scannapieco,
Daniel Kasen
Abstract:
We carry out direct numerical simulations of turbulent astrophysical media that explicitly track ionizations, recombinations, and species-by-species radiative cooling. The simulations assume solar composition and follows the evolution of hydrogen, helium, carbon, oxygen, sodium, and magnesium, but they do not include the presence of an ionizing background. In this case, the medium reaches a global…
▽ More
We carry out direct numerical simulations of turbulent astrophysical media that explicitly track ionizations, recombinations, and species-by-species radiative cooling. The simulations assume solar composition and follows the evolution of hydrogen, helium, carbon, oxygen, sodium, and magnesium, but they do not include the presence of an ionizing background. In this case, the medium reaches a global steady state that is purely a function of the one-dimensional turbulent velocity dispersion, $σ_{\rm 1D},$ and the product of the mean density and the driving scale of turbulence, $n L.$ Our simulations span a grid of models with $σ_{\rm 1D}$ ranging from 6 to 58 km s$^{-1}$ and $n L$ ranging from 10$^{16}$ to 10$^{20}$ cm$^{-2},$ which correspond to turbulent Mach numbers from $M=0.2$ to 10.6. The species abundances are well described by single-temperature estimates whenever $M$ is small, but local equilibrium models can not accurately predict the global equilibrium abundances when $M \gtrsim 1.$ To allow future studies to account for nonequilibrium effects in turbulent media, we gather our results into a series of tables, which we will extend in the future to encompass a wider range of elements, compositions, and ionizing processes.
△ Less
Submitted 4 February, 2015; v1 submitted 3 February, 2015;
originally announced February 2015.
-
Formation of Compact Clusters from High Resolution Hybrid Cosmological Simulations
Authors:
Mark L. A. Richardson,
Evan Scannapieco,
William J. Gray
Abstract:
The early Universe hosted a large population of small dark matter `minihalos' that were too small to cool and form stars on their own. These existed as static objects around larger galaxies until acted upon by some outside influence. Outflows, which have been observed around a variety of galaxies, can provide this influence in such a way as to collapse, rather than disperse the minihalo gas. Gray…
▽ More
The early Universe hosted a large population of small dark matter `minihalos' that were too small to cool and form stars on their own. These existed as static objects around larger galaxies until acted upon by some outside influence. Outflows, which have been observed around a variety of galaxies, can provide this influence in such a way as to collapse, rather than disperse the minihalo gas. Gray & Scannapieco performed an investigation in which idealized spherically-symmetric minihalos were struck by enriched outflows. Here we perform high-resolution cosmological simulations that form realistic minihalos, which we then extract to perform a large suite of simulations of outflow-minihalo interactions including non-equilibrium chemical reactions. In all models, the shocked minihalo forms molecules through non-equilibrium reactions, and then cools to form dense chemically homogenous clumps of star-forming gas. The formation of these high-redshift clusters will be observable with the next generation of telescopes, and the largest of them should survive to the present day, having properties similar to halo globular clusters.
△ Less
Submitted 9 October, 2013; v1 submitted 26 August, 2013;
originally announced August 2013.
-
Thermal and Chemical Evolution of Collapsing Filaments
Authors:
William J Gray,
Evan Scannapieco
Abstract:
Intergalactic filaments form the foundation of the cosmic web that connect galaxies together, and provide an important reservoir of gas for galaxy growth and accretion. Here we present very high resolution two-dimensional simulations of the thermal and chemical evolution of such filaments, making use of a 32 species chemistry network that tracks the evolution of key molecules formed from hydrogen,…
▽ More
Intergalactic filaments form the foundation of the cosmic web that connect galaxies together, and provide an important reservoir of gas for galaxy growth and accretion. Here we present very high resolution two-dimensional simulations of the thermal and chemical evolution of such filaments, making use of a 32 species chemistry network that tracks the evolution of key molecules formed from hydrogen, oxygen, and carbon. We study the evolution of filaments over a wide range of parameters including the initial density, initial temperature, strength of the dissociating UV background, and metallicity. In low-redshift, $Z \approx 0.1 Z_\odot $ filaments, the evolution is determined completely by the initial cooling time. If this is sufficiently short, the center of the filament always collapses to form dense, cold core containing a substantial fraction of molecules. In high-redshift, $Z=10^{-3} Z_\odot$ filaments, the collapse proceeds much more slowly. This is due mostly to the lower initial temperatures, which leads to a much more modest increase in density before the atomic cooling limit is reached, making subsequent molecular cooling much less efficient.
Finally, we study how the gravitational potential from a nearby dwarf galaxy affects the collapse of the filament and compare this to NGC 5253, a nearby starbusting dwarf galaxy thought to be fueled by the accretion of filament gas. In contrast to our fiducial case, a substantial density peak forms at the center of the potential. This peak evolves faster than the rest of the filament due to the increased rate at which chemical species form and cooling occur. We find that we achieve similar accretion rates as NGC 5253 but our two-dimensional simulations do not recover the formation of the giant molecular clouds that are seen in radio observations.
△ Less
Submitted 3 April, 2013; v1 submitted 2 April, 2013;
originally announced April 2013.
-
Identification of a Fundamental Transition in a Turbulently-Supported Interstellar Medium
Authors:
Evan Scannapieco,
William J. Gray,
Liubin Pan
Abstract:
The interstellar medium in star-forming galaxies is a multiphase gas in which turbulent support is at least as important as thermal pressure. Sustaining this configuration requires continuous radiative cooling, such that the overall average cooling rate matches the decay rate of turbulent energy into the medium. Here we carry out a set of numerical simulations of a stratified, turbulently stirred,…
▽ More
The interstellar medium in star-forming galaxies is a multiphase gas in which turbulent support is at least as important as thermal pressure. Sustaining this configuration requires continuous radiative cooling, such that the overall average cooling rate matches the decay rate of turbulent energy into the medium. Here we carry out a set of numerical simulations of a stratified, turbulently stirred, radiatively cooled medium, which uncover a fundamental transition at a critical one-dimensional turbulent velocity of ~ 35 km/s. At turbulent velocities below ~35 km/s, corresponding to temperatures below 300,000 K, the medium is stable, as the time for gas to cool is roughly constant as a function of temperature. On the other hand, at turbulent velocities above the critical value, the gas is shocked into an unstable regime in which the cooling time increases strongly with temperature, meaning that a substantial fraction of the interstellar medium is unable to cool on a turbulent dissipation timescale. This naturally leads to runaway heating and ejection of gas from any stratified medium with a one-dimensional turbulent velocity above ~35 km/s, a result that has implications for galaxy evolution at all redshifts.
△ Less
Submitted 1 December, 2011;
originally announced December 2011.
-
Formation of Compact Stellar Clusters by High-Redshift Galaxy Outflows III: Observability and Connection to Halo Globular Clusters
Authors:
William J Gray,
Evan Scannapieco
Abstract:
The early universe hosted a large population of low-mass virialized "minihalos," that were not massive enough to form stars on their own. While most minihalos were photoevaporated by ionizing photons from star-forming galaxies, these galaxies also drove large outflows, which in some cases would have reached the minihalos in advance of ionization fronts. In the previous papers in this series, we ca…
▽ More
The early universe hosted a large population of low-mass virialized "minihalos," that were not massive enough to form stars on their own. While most minihalos were photoevaporated by ionizing photons from star-forming galaxies, these galaxies also drove large outflows, which in some cases would have reached the minihalos in advance of ionization fronts. In the previous papers in this series, we carried out high-resolution, three-dimensional adaptive mesh refinement simulations of outflow-minihalo interactions that included non-equilibrium chemistry, radiative cooling, and turbulent mixing. We found that, for a fiducial set of parameters, minihalos were transformed into dense, chemically homogenous stellar clusters. Here we conduct a suite of simulations that follow these interactions over a wide range of parameters including minihalo mass, minihalo formation redshift, outflow energy, outflow redshift, distance, concentration, and spin. In almost all cases, the shocked minihalos form molecules through nonequillibrium reactions and then cool rapidly to become compact, chemically-homogenous stellar clusters. Furthermore, we show that the unique properties of these clusters make them a prime target for direct study with the next generation of telescopes, and that there are many reasons to suspect that their low-redshift counterparts are the observed population of halo globular clusters.
△ Less
Submitted 4 October, 2011; v1 submitted 7 September, 2011;
originally announced September 2011.
-
Formation of Compact Stellar Clusters by High-Redshift Galaxy Outflows I: Nonequillibrium Coolant Formation
Authors:
William J Gray,
Evan Scannapieco
Abstract:
We use high-resolution three-dimensional adaptive mesh refinement simulations to investigate the interaction of high-redshift galaxy outflows with low-mass virialized clouds of primordial composition. While atomic cooling allows star formation in objects with virial temperatures above $10^4$ K, "minihaloes" below this threshold are generally unable to form stars by themselves. However, these objec…
▽ More
We use high-resolution three-dimensional adaptive mesh refinement simulations to investigate the interaction of high-redshift galaxy outflows with low-mass virialized clouds of primordial composition. While atomic cooling allows star formation in objects with virial temperatures above $10^4$ K, "minihaloes" below this threshold are generally unable to form stars by themselves. However, these objects are highly susceptible to triggered star formation, induced by outflows from neighboring high-redshift starburst galaxies. Here we conduct a study of these interactions, focusing on cooling through non-equilibrium molecular hydrogen (H$_2$) and hydrogen deuteride (HD) formation. Tracking the non-equilibrium chemistry and cooling of 14 species and including the presence of a dissociating background, we show that shock interactions can transform minihaloes into extremely compact clusters of coeval stars. Furthermore, these clusters are all less than $\approx 10^6 M_\odot,$ and they are ejected from their parent dark matter halos: properties that are remarkably similar to those of the old population of globular clusters.
△ Less
Submitted 27 May, 2010;
originally announced May 2010.