skip to main content

Title: The AGORA High-resolution Galaxy Simulations Comparison Project. VI. Similarities and Differences in the Circumgalactic Medium

We analyze the circumgalactic medium (CGM) for eight commonly-used cosmological codes in the AGORA collaboration. The codes are calibrated to use identical initial conditions, cosmology, heating and cooling, and star formation thresholds, but each evolves with its own unique code architecture and stellar feedback implementation. Here, we analyze the results of these simulations in terms of the structure, composition, and phase dynamics of the CGM. We show properties such as metal distribution, ionization levels, and kinematics are effective tracers of the effects of the different code feedback and implementation methods, and as such they can be highly divergent between simulations. This is merely a fiducial set of models, against which we will in the future compare multiple feedback recipes for each code. Nevertheless, we find that the large parameter space these simulations establish can help disentangle the different variables that affect observable quantities in the CGM, e.g., showing that abundances for ions with higher ionization energy are more strongly determined by the simulation’s metallicity, while abundances for ions with lower ionization energy are more strongly determined by the gas density and temperature.

more » « less
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; « less
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Medium: X Size: Article No. 29
["Article No. 29"]
Sponsoring Org:
National Science Foundation
More Like this

    The coupling state between ions and neutrals in the interstellar medium plays a key role in the dynamics of magnetohydrodynamic (MHD) turbulence, but is challenging to study numerically. In this work, we investigate the damping of MHD turbulence in a partially ionized medium using 3D two-fluid (ions + neutrals) simulations generated with the athenak code. Specifically, we examine the velocity, density, and magnetic field statistics of the two-fluid MHD turbulence in different regimes of neutral-ion coupling. Our results demonstrate that when ions and neutrals are strongly coupled, the velocity statistics resemble those of single-fluid MHD turbulence. Both the velocity structures and kinetic energy spectra of ions and neutrals are similar, while their density structures can be significantly different. With an excess of small-scale sharp density fluctuations in ions, the density spectrum in ions is shallower than that of neutrals. When ions and neutrals are weakly coupled, the turbulence in ions is more severely damped due to the ion-neutral collisional friction than that in neutrals, resulting in a steep kinetic energy spectrum and density spectrum in ions compared to the Kolmogorov spectrum. We also find that the magnetic energy spectrum basically follows the shape of the kinetic energy spectrum of ions, irrespective of the coupling regime. In addition, we find large density fluctuations in ions and neutrals and thus spatially inhomogeneous ionization fractions. As a result, the neutral-ion decoupling and damping of MHD turbulence take place over a range of length-scales.

    more » « less

    Modelling galaxy formation in hydrodynamic simulations has increasingly adopted various radiative transfer methods to account for photoionization feedback from young massive stars. However, the evolution of H ii regions around stars begins in dense star-forming clouds and spans large dynamical ranges in both space and time, posing severe challenges for numerical simulations in terms of both spatial and temporal resolution that depends strongly on gas density (∝n−1). In this work, we perform a series of idealized H ii region simulations using the moving-mesh radiation-hydrodynamic code arepo-rt to study the effects of numerical resolution. The simulated results match the analytical solutions and the ionization feedback converges only if the Strömgren sphere is resolved by at least 10–100 resolution elements and the size of each time integration step is smaller than 0.1 times the recombination time-scale. Insufficient spatial resolution leads to reduced ionization fraction but enhanced ionized gas mass and momentum feedback from the H ii regions, as well as degrading the multiphase interstellar medium into a diffuse, partially ionized, warm (∼8000 K) gas. On the other hand, insufficient temporal resolution strongly suppresses the effects of ionizing feedback. This is because longer time-steps are not able to resolve the rapid variation of the thermochemistry properties of the gas cells around massive stars, especially when the photon injection and thermochemistry are performed with different cadences. Finally, we provide novel numerical implementations to overcome the above issues when strict resolution requirements are not achievable in practice.

    more » « less

    A key assumption in quasar absorption-line studies of the circumgalactic medium (CGM) is that each absorption component maps to a spatially isolated ‘cloud’ structure that has single valued properties (e.g. density, temperature, metallicity). We aim to assess and quantify the degree of accuracy underlying this assumption. We used adaptive mesh refinement hydrodynamic cosmological simulations of two z = 1 dwarf galaxies and generated synthetic quasar absorption-line spectra of their CGM. For the Si ii λ1260 transition, and the C iv λλ1548, 1550 and O vi λλ1031, 1037 fine-structure doublets, we objectively determined which gas cells along a line of sight (LOS) contribute to detected absorption. We implemented a fast, efficient, and objective method to define individual absorption components in each absorption profile. For each absorption component, we quantified the spatial distribution of the absorbing gas. We studied a total of 1302 absorption systems containing a total of 7755 absorption components. 48  per cent of Si ii, 68  per cent of C iv, and 72  per cent of O vi absorption components arise from two or more spatially isolated ‘cloud’ structures along the LOS. Spatially isolated ‘cloud’ structures were most likely to have cloud–cloud LOS separations of 0.03Rvir (1.3 kpc), 0.11Rvir (4.8 kpc), and 0.13Rvir (5.6 kpc) for Si ii, C iv, and O vi, respectively. There can be very little overlap between multiphase gas structures giving rise to absorption components. If our results reflect the underlying reality of how absorption lines record CGM gas, they place tension on current observational analysis methods as they suggest that component-by-component absorption-line formation is more complex than is assumed and applied for chemical-ionization modelling.

    more » « less

    Arkenstone is a new model for multiphase, stellar feedback-driven galactic winds designed for inclusion in coarse resolution cosmological simulations. In this first paper of a series, we describe the features that allow Arkenstone to properly treat high specific energy wind components and demonstrate them using idealized non-cosmological simulations of a galaxy with a realistic circumgalactic medium (CGM), using the arepo code. Hot, fast gas phases with low mass loadings are predicted to dominate the energy content of multiphase outflows. In order to treat the huge dynamic range of spatial scales involved in cosmological galaxy formation at feasible computational expense, cosmological volume simulations typically employ a Lagrangian code or else use adaptive mesh refinement with a quasi-Lagrangian refinement strategy. However, it is difficult to inject a high specific energy wind in a Lagrangian scheme without incurring artificial burstiness. Additionally, the low densities inherent to this type of flow result in poor spatial resolution. Arkenstone addresses these issues with a novel scheme for coupling energy into the transition region between the interstellar medium (ISM) and the CGM, while also providing refinement at the base of the wind. Without our improvements, we show that poor spatial resolution near the sonic point of a hot, fast outflow leads to an underestimation of gas acceleration as the wind propagates. We explore the different mechanisms by which low and high specific energy winds can regulate the star formation rate of galaxies. In future work, we will demonstrate other aspects of the Arkenstone model.

    more » « less

    Hydrogen emission lines can provide extensive information about star-forming galaxies in both the local and high-redshift Universe. We present a detailed Lyman continuum (LyC), Lyman-α (Lyα), and Balmer line (Hα and Hβ) radiative transfer study of a high-resolution isolated Milky Way simulation using the state-of-the-art Arepo-RT radiation hydrodynamics code with the SMUGGLE galaxy formation model. The realistic framework includes stellar feedback, non-equilibrium thermochemistry accounting for molecular hydrogen, and dust grain evolution in the interstellar medium (ISM). We extend our publicly available Cosmic Lyα Transfer (COLT) code with photoionization equilibrium Monte Carlo radiative transfer and various methodology improvements for self-consistent end-to-end (non-)resonant line predictions. Accurate LyC reprocessing to recombination emission requires modelling pre-absorption by dust ($f_\text{abs} \approx 27.5\,\rm{per\,\,cent}$), helium ionization ($f_\text{He} \approx 8.7\,\rm{per\,\,cent}$), and anisotropic escape fractions ($f_\text{esc} \approx 7.9\,\rm{per\,\,cent}$), as these reduce the available budget for hydrogen line emission ($f_\text{H} \approx 55.9\,\rm{per\,\,cent}$). We investigate the role of the multiphase dusty ISM, disc geometry, gas kinematics, and star formation activity in governing the physics of emission and escape, focusing on the time variability, gas-phase structure, and spatial spectral, and viewing angle dependence of the emergent photons. Isolated disc simulations are well-suited for comprehensive observational comparisons with local Hα surveys, but would require a proper cosmological circumgalactic medium (CGM) environment as well as less dust absorption and rotational broadening to serve as analogs for high-redshift Lyα emitting galaxies. Future applications of our framework to next-generation cosmological simulations of galaxy formation including radiation-hydrodynamics that resolve ≲10 pc multiphase ISM and ≲1 kpc CGM structures will provide crucial insights and predictions for current and upcoming Lyα observations.

    more » « less