BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Alexander J. Richings

Alexander J. Richings appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.CO

Catalog footprint

What is connected

9works
2topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2026arXiv

Modelling the evolution and influence of dust in cosmological simulations that include the cold phase of the interstellar medium

While marginal in mass terms, dust grains play an outsized role in both the physics and observation of the interstellar medium (ISM). However, explicit modelling of this ISM constituent remains uncommon in large cosmological simulations. In this work, we present a model for the life-cycle of dust in the ISM that couples to the forthcoming COLIBRE galaxy formation model, which explicitly simulates the cold ISM. We follow 6 distinct grain types: 3 chemical species, including carbon and two silicate grains, with 2 size bins each. Our dust model accounts for seeding of grains from stellar ejecta, self-consistent element-by-element metal yields and growth by accretion, grain size transfer (shattering and coagulation) and destruction of dust by thermal sputtering in the ISM. We detail the calibration of this model, particularly the use of a clumping factor, to account for unresolved gas clouds in which dust readily evolves. We present a fiducial run in a 25$^3$~cMpc$^3$ cosmological volume that displays good agreement with observations of the cosmic evolution of dust density, as well as the $z=0$ galaxy dust mass function and dust scaling relations. We highlight known tensions between observational datasets of the dust-to-gas ratio as a function of metallicity depending on which metallicity calibrator is used; our model favours higher-normalisation metallicity calibrators, which agree with the observations within 0.1~dex for stellar masses $>10^9 \; {\rm M_\odot}$. We compare the grain size distribution to observations of local galaxies, and find that our simulation suggests a higher concentration of small grains, associated with more diffuse ISM and the warm-neutral medium (WNM), which both play a key role in boosting H$_2$ content. Putting these results and modelling approaches in context, we set the stage for upcoming insights into the dusty ISM of galaxies using the COLIBRE simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Rapid disc settling and the transition from bursty to steady star formation in Milky Way-mass galaxies

Recent observations and simulations indicate substantial evolution in the properties of galaxies with time, wherein rotationally-supported and steady thin discs (like those frequently observed in the local universe) emerge from galaxies that are clumpy, irregular, and have bursty star formation rates (SFRs). To better understand the progenitors of local disc galaxies we carry out an analysis of three FIRE-2 simulated galaxies with a mass similar to the Milky Way at redshift z=0. We show that all three galaxies transition from bursty to steady SFRs at a redshift between z=0.5 and z=0.8, and that this transition coincides with a rapid (< ~1 Gyr) emergence of a rotationally-supported interstellar medium (ISM).In the late phase with steady SFR, the rotational energy comprises > ~90% of the total kinetic + thermal energy in the ISM, and is roughly half the gravitational energy. By contrast, during the early phase with bursty star formation, the ISM has a quasi-spheroidal morphology and its energy budget is dominated by quasi-isotropic flows including turbulence and coherent inflows/outflows. This result, that rotational support is subdominant at early times, challenges the common application of equilibrium disc models to the high-redshift progenitors of Milky Way-like galaxies. We further find that the formation of a rotation-supported ISM coincides with the formation of a thermal energy-supported inner circumgalactic medium (CGM). Before this transition, the inner CGM is also supported by turbulence and coherent flows, indicating that at early times there is no clear boundary between the ISM and inner CGM.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The effects of local stellar radiation and dust depletion on non-equilibrium interstellar chemistry

Interstellar chemistry is important for galaxy formation, as it determines the rate at which gas can cool, and enables us to make predictions for observable spectroscopic lines from ions and molecules. We explore two central aspects of modelling the chemistry of the interstellar medium (ISM): (1) the effects of local stellar radiation, which ionises and heats the gas, and (2) the depletion of metals onto dust grains, which reduces the abundance of metals in the gas phase. We run high-resolution (400 M$_\odot$ per baryonic particle) simulations of isolated disc galaxies, from dwarfs to Milky Way-mass, using the FIRE galaxy formation models together with the CHIMES non-equilibrium chemistry and cooling module. In our fiducial model, we couple the chemistry to the stellar fluxes calculated from star particles using an approximate radiative transfer scheme, and we implement an empirical density-dependent prescription for metal depletion. For comparison, we also run simulations with a spatially uniform radiation field, and without metal depletion. Our fiducial model broadly reproduces observed trends in HI and H2 mass with stellar mass, and in line luminosity versus star formation rate for [CII] 158$μ$m, [OI] 63$μ$m, [OIII] 88$μ$m, [NII] 122$μ$m and H$α$ 6563A. Our simulations with a uniform radiation field predict fainter luminosities, by up to an order of magnitude for [OIII] 88$μ$m and H$α$ 6563A, while ignoring metal depletion increases the luminosity of carbon and oxygen lines by a factor $\approx$2. However, the overall evolution of the galaxy is not strongly affected by local stellar fluxes or metal depletion, except in dwarf galaxies where the inclusion of local fluxes leads to weaker outflows and hence higher gas fractions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Unravelling the physics of multiphase AGN winds through emission line tracers

Observations of emission lines in Active Galactic Nuclei (AGN) often find fast (~1000 km s^-1) outflows extending to kiloparsec scales, seen in ionised, neutral atomic and molecular gas. In this work we present radiative transfer calculations of emission lines in hydrodynamic simulations of AGN outflows driven by a hot wind bubble, including non-equilibrium chemistry, to explore how these lines trace the physical properties of the multiphase outflow. We find that the hot bubble compresses the line-emitting gas, resulting in higher pressures than in the ambient ISM or that would be produced by the AGN radiation pressure. This implies that observed emission line ratios such as [OIV] 25 $μ$m / [NeII] 12 $μ$m , [NeV] 14 $μ$m / [NeII] 12 $μ$m and [NIII] 57 $μ$m / [NII] 122 $μ$m constrain the presence of the bubble and hence the outflow driving mechanism. However, the line-emitting gas is under-pressurised compared to the hot bubble itself, and much of the line emission arises from gas that is out of pressure, thermal and/or chemical equilibrium. Our results thus suggest that assuming equilibrium conditions, as commonly done in AGN line emission models, is not justified if a hot wind bubble is present. We also find that >50 per cent of the mass outflow rate, momentum flux and kinetic energy flux of the outflow are traced by lines such as [NII] 122 $μ$m and [NeIII] 15 $μ$m (produced in the 10^4 K phase) and [CII] 158 $μ$m (produced in the transition from 10^4 K to 100 K).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Bimodality of low-redshift circumgalactic O VI in non-equilibrium EAGLE zoom simulations

We introduce a series of 20 cosmological hydrodynamical simulations of Lstar (M_200 =10^11.7 - 10^12.3 Msol) and group-sized (M_200 = 10^12.7 - 10^13.3 Msol) haloes run with the model used for the EAGLE project, which additionally includes a non-equilibrium ionization and cooling module that follows 136 ions. The simulations reproduce the observed correlation, revealed by COS-Halos at z~0.2, between O VI column density at impact parameters b < 150 kpc and the specific star formation rate (sSFR=SFR/Mstar) of the central galaxy at z~0.2. We find that the column density of circumgalactic O VI is maximal in the haloes associated with Lstar galaxies, because their virial temperatures are close to the temperature at which the ionization fraction of O VI peaks (T~10^5.5 K). The higher virial temperature of group haloes (> 10^6 K) promotes oxygen to higher ionization states, suppressing the O VI column density. The observed NO VI-sSFR correlation therefore does not imply a causal link, but reflects the changing characteristic ionization state of oxygen as halo mass is increased. In spite of the mass-dependence of the oxygen ionization state, the most abundant circumgalactic oxygen ion in both Lstar and group haloes is O VII; O VI accounts for only 0.1% of the oxygen in group haloes and 0.9-1.3% with Lstar haloes. Nonetheless, the metals traced by O VI absorbers represent a fossil record of the feedback history of galaxies over a Hubble time; their characteristic epoch of ejection corresponds to z > 1 and much of the ejected metal mass resides beyond the virial radius of galaxies. For both Lstar and group galaxies, more of the oxygen produced and released by stars resides in the circumgalactic medium (within twice the virial radius) than in the stars and ISM of the galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Chemical evolution of giant molecular clouds in simulations of galaxies

We present an analysis of Giant Molecular Clouds (GMCs) within hydrodynamic simulations of isolated, low-mass (M* ~ 10^9 M_sol) disc galaxies. We study the evolution of molecular abundances and the implications for CO emission and the X_CO conversion factor in individual clouds. We define clouds either as regions above a density threshold n_H,min = 10 cm^-3, or using an observationally motivated CO intensity threshold of 0.25 K km s^-1. Our simulations include a non-equilibrium chemical model with 157 species, including 20 molecules. We also investigate the effects of resolution and pressure floors (i.e. Jeans limiters). We find cloud lifetimes up to ~40 Myr, with a median of 13 Myr, in agreement with observations. At one tenth solar metallicity, young clouds (<10-15 Myr) are underabundant in H2 and CO compared to chemical equilibrium, by factors of ~3 and 1-2 orders of magnitude, respectively. At solar metallicity, GMCs reach chemical equilibrium faster (within ~1 Myr). We also compute CO emission from individual clouds. The mean CO intensity, I_CO, is strongly suppressed at low dust extinction, A_v, and possibly saturates towards high A_v, in agreement with observations. The I_CO - A_v relation shifts towards higher A_v for higher metallicities and, to a lesser extent, for stronger UV radiation. At one tenth solar metallicity, CO emission is weaker in young clouds (<10-15 Myr), consistent with the underabundance of CO. Consequently, X_CO decreases by an order of magnitude from 0 to 15 Myr, albeit with a large scatter.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The effects of metallicity, UV radiation and non-equilibrium chemistry in high-resolution simulations of galaxies

We present a series of hydrodynamic simulations of isolated galaxies with stellar mass of $10^{9} \, \rm{M}_{\odot}$. The models use a resolution of $750 \, \rm{M}_{\odot}$ per particle and include a treatment for the full non-equilibrium chemical evolution of ions and molecules (157 species in total), along with gas cooling rates computed self-consistently using the non-equilibrium abundances. We compare these to simulations evolved using cooling rates calculated assuming chemical (including ionisation) equilibrium, and we consider a wide range of metallicities and UV radiation fields, including a local prescription for self-shielding by gas and dust. We find higher star formation rates and stronger outflows at higher metallicity and for weaker radiation fields, as gas can more easily cool to a cold (few hundred Kelvin) star forming phase under such conditions. Contrary to variations in the metallicity and the radiation field, non-equilibrium chemistry generally has no strong effect on the total star formation rates or outflow properties. However, it is important for modelling molecular outflows. For example, the mass of H$_{2}$ outflowing with velocities $> 50 \, \rm{km} \, \rm{s}^{-1}$ is enhanced by a factor $\sim 20$ in non-equilibrium. We also compute the observable line emission from CII and CO. Both are stronger at higher metallicity, while CII and CO emission are higher for stronger and weaker radiation fields respectively. We find that CII is generally unaffected by non-equilibrium chemistry. However, emission from CO varies by a factor of $\sim 2 - 4$. This has implications for the mean $X_{\rm{CO}}$ conversion factor between CO emission and H$_{2}$ column density, which we find is lowered by up to a factor $\sim 2.3$ in non-equilibrium, and for the fraction of CO-dark molecular gas.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Non-equilibrium chemistry and cooling in the diffuse interstellar medium - I. Optically thin regime

An accurate treatment of the multiphase interstellar medium (ISM) in hydrodynamic galaxy simulations requires that we follow not only the thermal evolution of the gas, but also the evolution of its chemical state, including its molecular chemistry, without assuming chemical (including ionisation) equilibrium. We present a reaction network that can be used to solve for this thermo-chemical evolution. Our model follows the evolution of all ionisation states of the 11 elements that dominate the cooling rate, along with important molecules such as H2 and CO, and the intermediate molecular species that are involved in their formation (20 molecules in total). We include chemical reactions on dust grains, thermal processes involving dust, cosmic ray ionisation and heating and photochemical reactions. We focus on conditions typical for the diffuse ISM, with densities of 10^-2 cm^-3 < nH < 10^4 cm^-3 and temperatures of 10^2 K < T < 10^4 K, and we consider a range of radiation fields, including no UV radiation. In this paper we consider only gas that is optically thin, while paper II considers gas that becomes shielded from the radiation field. We verify the accuracy of our model by comparing chemical abundances and cooling functions in chemical equilibrium with the photoionisation code Cloudy. We identify the major coolants in diffuse interstellar gas to be CII, SiII and FeII, along with OI and H2 at densities nH > 10^2 cm^-3. Finally, we investigate the impact of non-equilibrium chemistry on the cooling functions of isochorically or isobarically cooling gas. We find that, at T < 10^4 K, recombination lags increase the electron abundance above its equilibrium value at a given temperature, which can enhance the cooling rate by up to two orders of magnitude. The cooling gas also shows lower H2 abundances than in equilibrium, by up to an order of magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Non-equilibrium chemistry and cooling in the diffuse interstellar medium - II. Shielded gas

We extend the non-equilibrium model for the chemical and thermal evolution of diffuse interstellar gas presented in Richings et al. (2014) to account for shielding from the UV radiation field. We attenuate the photochemical rates by dust and by gas, including absorption by HI, H2, HeI, HeII and CO where appropriate. We then use this model to investigate the dominant cooling and heating processes in interstellar gas as it becomes shielded from the UV radiation. We consider a one-dimensional plane-parallel slab of gas irradiated by the interstellar radiation field, either at constant density and temperature or in thermal and pressure equilibrium. The dominant thermal processes tend to form three distinct regions in the clouds. At low column densities cooling is dominated by ionised metals such as SiII, FeII, FeIII and CII, which are balanced by photoheating, primarily from HI. Once the hydrogen-ionising radiation becomes attenuated by neutral hydrogen, photoelectric dust heating dominates, while CII becomes dominant for cooling. Finally, dust shielding triggers the formation of CO and suppresses photoelectric heating. The dominant coolants in this fully shielded region are H2 and CO. The column density of the HI-H2 transition predicted by our model is lower at higher density (or at higher pressure for gas clouds in pressure equilibrium) and at higher metallicity, in agreement with previous PDR models. We also compare the HI-H2 transition in our model to two prescriptions for molecular hydrogen formation that have been implemented in hydrodynamic simulations.

0 citations0 reviewsNo code linkedOpen access