Trident: a universal tool for generating synthetic absorption spectra from astrophysical simulations

Hydrodynamical simulations are increasingly able to accurately model physical systems on stellar, galactic, and cosmological scales, however, the utility of these simulations is often limited by our ability to directly compare them with the datasets produced by observers: spectra, photometry, etc. To address this problem, we have created Trident}, a Python-based, open-source tool for post-processing hydrodynamical simulations to produce synthetic absorption spectra and related data. Trident} can (i) create absorption-line spectra for any trajectory through a simulated dataset mimicking both background quasar and down-the-barrel configurations, (ii) reproduce the spectral characteristics of common instruments like the Cosmic Origins Spectrograph, (iii) operate across the ultraviolet, optical and infrared using customizable absorption line lists, (iv) trace simulated physical structures directly to spectral features, (v) approximate the presence of ion species absent from the simulation outputs, (vi) generate column density maps for any ion, and (vii) provide support for all major astrophysical hydrodynamical codes. The focus of Trident's development is for using simulated datasets to better interpret observations of the circumgalactic medium (CGM) and intergalactic medium (IGM), but it remains a general tool applicable in other contexts.Comment: 16 pages, 13 figures, published in ApJ, Code available at http://trident-project.or

Comparing Simulations and Observations of Galaxy Evolution: Methods for Constraining the Nature of Stellar Feedback

Computational hydrodynamical simulations are a very useful tool for understanding how galaxies form and evolve over cosmological timescales not easily revealed through observations. However, they are only useful if they reproduce the sorts of galaxies that we see in the real universe. One of the ways in which simulations of this sort tend to fail is in the prescription of stellar feedback, the process by which nascent stars return material and energy to their immediate environments. Careful treatment of this interaction in subgrid models, so-called because they operate on scales below the resolution of the simulation, is crucial for the development of realistic galaxy models. Equally important is developing effective methods for comparing simulation data against observations to ensure galaxy models which mimic reality and inform us about natural phenomena. This thesis examines the formation and evolution of galaxies and the observable characteristics of the resulting systems. We employ extensive use of cosmological hydrodynamical simulations in order to simulate and interpret the evolution of massive spiral galaxies like our own Milky Way. First, we create a method for producing synthetic photometric images of grid-based hydrodynamical models for use in a direct comparison against observations in a variety of filter bands. We apply this method to a simulation of a cluster of galaxies to investigate the nature of the red-sequence/blue-cloud dichotomy in the galaxy color-magnitude diagram. Second, we implement several subgrid models governing the complex behavior of gas and stars on small scales in our galaxy models. Several numerical simulations are conducted with similar initial conditions, where we systematically vary the subgrid models, afterward assessing their efficacy through comparisons of their internal kinematics with observed systems. Third, we generate an additional method to compare observations with simulations, focusing on the tenuous circumgalactic medium. Informed by our previous studies, we investigate the sensitivity of this new mode of comparison to hydrodynamical subgrid prescription. Finally, we synthesize the results of these studies and identify future avenues of research

The AGORA High-resolution Galaxy Simulations Comparison Project. II. Isolated Disk Test

Using an isolated Milky Way-mass galaxy simulation, we compare results from nine state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt–Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low-density region, and between more diffusive and less diffusive schemes in the high-density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes

On the survival of cool clouds in the circumgalactic medium

We explore the survival of cool clouds in multiphase circumgalactic media. We revisit the ‘cloud-crushing problem’ in a large survey of simulations including radiative cooling, self-shielding, self-gravity, magnetic fields, and anisotropic Braginskii conduction and viscosity (with saturation). We explore a wide range of parameters including cloud size, velocity, ambient temperature and density, and a variety of magnetic field configurations and cloud turbulence. We find that realistic magnetic fields and turbulence have weaker effects on cloud survival; the most important physics is radiative cooling and conduction. Self-gravity and self-shielding are important for clouds that are initially Jeans-unstable, but largely irrelevant otherwise. Non-self-gravitating, realistically magnetized clouds separate into four regimes: (1) at low column densities, clouds evaporate rapidly via conduction; (2) a ‘failed pressure confinement’ regime, where the ambient hot gas cools too rapidly to provide pressure confinement for the cloud; (3) an ‘infinitely long-lived’ regime, in which the cloud lifetime becomes longer than the cooling time of gas swept up in the leading bow shock, so the cloud begins to accrete and grow; and (4) a ‘classical cloud destruction’ regime, where clouds are eventually destroyed by instabilities. In the final regime, the cloud lifetime can exceed the naive cloud-crushing time owing to conduction-induced compression. However, small and/or slow-moving clouds can also evaporate more rapidly than the cloud-crushing time. We develop simple analytic models that explain the simulated cloud destruction times in this regime

Properties of the circumgalactic medium in cosmic ray-dominated galaxy haloes

We investigate the impact of cosmic rays (CRs) on the circumgalactic medium (CGM) in FIRE-2 simulations, for ultra-faint dwarf through Milky Way (MW)-mass haloes hosting star-forming (SF) galaxies. Our CR treatment includes injection by supernovae, anisotropic streaming and diffusion along magnetic field lines, and collisional and streaming losses, with constant parallel diffusivity κ∼3×10²⁹ cm² s⁻¹ chosen to match γ-ray observations. With this, CRs become more important at larger halo masses and lower redshifts, and dominate the pressure in the CGM in MW-mass haloes at z ≲ 1–2. The gas in these ‘CR-dominated’ haloes differs significantly from runs without CRs: the gas is primarily cool (a few ∼10⁴), and the cool phase is volume-filling and has a thermal pressure below that needed for virial or local thermal pressure balance. Ionization of the ‘low’ and ‘mid’ ions in this diffuse cool gas is dominated by photoionization, with O VI columns ≳10^(14.5) cm⁻² at distances ≳150kpc⁠. CR and thermal gas pressure are locally anticorrelated, maintaining total pressure balance, and the CGM gas density profile is determined by the balance of CR pressure gradients and gravity. Neglecting CRs, the same haloes are primarily warm/hot (⁠T≳10⁵) with thermal pressure balancing gravity, collisional ionization dominates, O VI columns are lower and Ne VIII higher, and the cool phase is confined to dense filaments in local thermal pressure equilibrium with the hot phase

CloudFlex: A Flexible Parametric Model for the Small-Scale Structure of the Circumgalactic Medium

We present CloudFlex, a new open-source tool for predicting the absorption-line signatures of cool gas in galaxy halos with complex small-scale structure. Motivated by analyses of cool material in hydrodynamical simulations of turbulent, multiphase media, we model individual cool gas structures as assemblies of cloudlets with a power-law distribution of cloudlet mass $\propto m_{\rm cl}^{-\alpha}$ and relative velocities drawn from a turbulent velocity field. The user may specify $\alpha$, the lower limit of the cloudlet mass distribution ($m_{\rm cl,min}$), and several other parameters that set the total mass, size, and velocity distribution of the complex. We then calculate the MgII 2796 absorption profiles induced by the cloudlets along pencil-beam lines of sight. We demonstrate that at fixed metallicity, the covering fraction of sightlines with equivalent widths $W_{2796} < 0.3$ Ang increases significantly with decreasing $m_{\rm cl,min}$, cool cloudlet number density ($n_{\rm cl}$), and cloudlet complex size. We then present a first application, using this framework to predict the projected $W_{2796}$ distribution around ${\sim}L^*$ galaxies. We show that the observed incidences of $W_{2796}>0.3$ Ang sightlines within 10 kpc < $R_{\perp}$ < 50 kpc are consistent with our model over much of parameter space. However, they are underpredicted by models with $m_{\rm cl,min}\ge100M_{\odot}$ and $n_{\rm cl}\ge0.03$ $\rm cm^{-3}$, in keeping with a picture in which the inner cool circumgalactic medium (CGM) is dominated by numerous low-mass cloudlets ($m_{\rm cl}\lesssim100M_{\odot}$) with a volume filling factor ${\lesssim}1\%$. When used to simultaneously model absorption-line datasets built from multi-sightline and/or spatially-extended background probes, CloudFlex will enable detailed constraints on the size and velocity distributions of structures comprising the photoionized CGM.Comment: 22 pages, 7 figures. Submitted to AAS Journals, with minor modifications. Comments welcome. (1) Co-first authors who made equal contributions to this wor