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.

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 {\lambda}1549 and O VI {\lambda}1037 emission. For density bounded photoionization, He II {\lambda}1640, {\lambda}4686, C III] {\lambda}1908, Si IV {\lambda}1206, and Si III {\lambda}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.

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.

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.

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 $\sigma_{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 $\sigma_{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.

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.

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.

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, $\sigma_{\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 $\sigma_{\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, $\sigma_{\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.

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.

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.

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.

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, $\sigma_{\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 $\sigma_{\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.

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.

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.

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.

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.

In the primordial universe, low mass structures with virial temperatures less than 10$^{4}$ K were unable to cool by atomic line transitions, leading to a strong suppression of star formation. On the other hand, these "minihalos" were highly prone to triggered star formation by interactions from nearby galaxy outflows. In Gray & Scannapieco (2010), we explored the impact of nonequilibrium chemistry on these interactions. Here we turn our attention to the role of metals, carrying out a series of high-resolution three-dimensional adaptive mesh refinement simulations that include both metal cooling and a subgrid turbulent mixing model. Despite the presence of an additional coolant, we again we find that outflow-minihalo interactions produce a distribution of dense, massive stellar clusters. We also find that these clusters are evenly enriched with metals to a final abundance of Z $\approx$ 10$^{-2}$ Z$_{\odot}$. As in our previous simulations, all of these properties suggest that these interactions may have given rise to present-day halo globular clusters.

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.