Galaxy Formation Physics Research Articles

Neutrino winds on the sky

We develop a first-principles formalism to compute the distortion to the relic neutrino density field caused by the peculiar motions of large-scale structures. This distortion slows halos down due to dynamical friction, causes a local anisotropy in the neutrino-CDM cross-correlation, and reduces the global cross-correlation between neutrinos and CDM. The local anisotropy in the neutrino-CDM cross-spectrum is imprinted in the three point cross-correlations of matter and galaxies, or the bispectrum in Fourier space, producing a signal peaking at squeezed triangle configurations. This bispectrum signature of neutrino masses is not limited by cosmic variance or potential inaccuracies in the modeling of complicated nonlinear and galaxy formation physics, and it is not degenerate with the optical depth to reionization. We show that future surveys have the potential to detect the distortion bispectrum.

LIMpy: A Semianalytic Approach to Simulating Multiline Intensity Maps at Millimeter Wavelengths

Mapping of multiple lines such as the fine-structure emission from [C ii] (157.7 μm), [O iii] (52 and 88.4 μm), and rotational emission lines from CO are of particular interest for upcoming line intensity mapping (LIM) experiments at millimeter wavelengths, due to their brightness features. Several upcoming experiments aim to cover a broad range of scientific goals, from detecting signatures of the epoch of reionization to the physics of star formation and its role in galaxy evolution. In this paper, we develop a semianalytic approach to modeling line strengths as functions of the star formation rate (SFR) or infrared luminosity based on observations of local and high-z galaxies. This package, LIMpy (Line Intensity Mapping in Python), estimates the intensity and power spectra of [C ii], [O iii], and CO rotational transition lines up to the J levels (1–0) to (13–12) based both on analytic formalism and on simulations. We develop a relation among halo mass, SFR, and multiline intensities that permits us to construct a generic formula for the evolution of several line strengths up to z ∼ 10. We implement a variety of star formation models and multiline luminosity relations to estimate the astrophysical uncertainties on the intensity power spectrum of these lines. As a demonstration, we predict the signal-to-noise ratio of [C ii] detection for an EoR-Spec-like instrument on the Fred Young Submillimeter Telescope. Furthermore, the ability to use any halo catalog allows the LIMpy code to be easily integrated into existing simulation pipelines, providing a flexible tool to study intensity mapping in the context of complex galaxy formation physics.

A Universal Equation to Predict Ωm from Halo and Galaxy Catalogs

We discover analytic equations that can infer the value of Ωm from the positions and velocity moduli of halo and galaxy catalogs. The equations are derived by combining a tailored graph neural network (GNN) architecture with symbolic regression. We first train the GNN on dark matter halos from Gadget N-body simulations to perform field-level likelihood-free inference, and show that our model can infer Ωm with ∼6% accuracy from halo catalogs of thousands of N-body simulations run with six different codes: Abacus, CUBEP3M, Gadget, Enzo, PKDGrav3, and Ramses. By applying symbolic regression to the different parts comprising the GNN, we derive equations that can predict Ωm from halo catalogs of simulations run with all of the above codes with accuracies similar to those of the GNN. We show that, by tuning a single free parameter, our equations can also infer the value of Ωm from galaxy catalogs of thousands of state-of-the-art hydrodynamic simulations of the CAMELS project, each with a different astrophysics model, run with five distinct codes that employ different subgrid physics: IllustrisTNG, SIMBA, Astrid, Magneticum, SWIFT-EAGLE. Furthermore, the equations also perform well when tested on galaxy catalogs from simulations covering a vast region in parameter space that samples variations in 5 cosmological and 23 astrophysical parameters. We speculate that the equations may reflect the existence of a fundamental physics relation between the phase-space distribution of generic tracers and Ωm, one that is not affected by galaxy formation physics down to scales as small as 10 h −1 kpc.

Simulating Atomic Dark Matter in Milky Way Analogs

Dark sector theories naturally lead to multicomponent scenarios for dark matter where a subcomponent can dissipate energy through self-interactions, allowing it to efficiently cool inside galaxies. We present the first cosmological hydrodynamical simulations of Milky Way analogs where the majority of dark matter is collisionless cold dark matter (CDM) but a subcomponent (6%) is strongly dissipative minimal atomic dark matter (ADM). The simulations, implemented in GIZMO and utilizing FIRE-2 galaxy formation physics to model the standard baryonic sector, demonstrate that the addition of even a small fraction of dissipative dark matter can significantly impact galactic evolution despite being consistent with current cosmological constraints. We show that ADM gas with roughly standard model–like masses and couplings can cool to form a rotating “dark disk” with angular momentum closely aligned with the visible stellar disk. The morphology of the disk depends sensitively on the parameters of the ADM model, which affect the cooling rates in the dark sector. The majority of the ADM gas gravitationally collapses into dark “clumps” (regions of black hole or mirror star formation), which form a prominent bulge and a rotating thick disk in the central galaxy. These clumps form early and quickly sink to the inner ∼kiloparsec of the galaxy, affecting the galaxy’s star formation history and present-day baryonic and CDM distributions.

Simulating the LOcal Web (SLOW)

Context. Several observations of the Local Universe point toward the existence of very prominent structures: massive galaxy clusters and local superclusters on the one hand, but also large local voids and underdensities on the other. However, it is highly nontrivial to connect such different observational selected tracers to the underlying dark matter (DM) distribution. Aims. Therefore, constructing mock catalogs of such observable tracers using cosmological hydrodynamics simulations is needed. These simulations have to follow galaxy formation physics and also have to be constrained to reproduce the Local Universe. Such constraints should be based on observables that directly probe the full underlying gravitational field, such as the observed peculiar velocity field, to provide an independent test on the robustness of these distinctive structures. Methods. We used a 500 h−1 Mpc constrained simulation of the Local Universe to investigate the anomalies in the local density field, as found in observations. Constructing the initial conditions based on peculiar velocities derived from the CosmicFlows-2 catalog makes the predictions of the simulations completely independent from the distribution of the observed tracer population, and following galaxy formation physics directly in the hydrodynamics simulations also allows the comparison to be based directly on the stellar masses of galaxies or X-ray luminosity of clusters. We also used the 2668 h−1 Mpc large cosmological box from the Magneticum simulations to evaluate the frequency of finding such anomalies in random patches within simulations. Results. We demonstrate that halos and galaxies in our constrained simulation trace the local dark matter density field very differently. Thus, this simulation reproduces the observed 50% underdensity of galaxy clusters and groups within the sphere of ≈100 Mpc when applying the same mass or X-ray luminosity limit used in the observed cluster sample (CLASSIX), which is consistent with a ≈1.5σ feature. At the same time, the simulation reproduces the observed overdensity of massive galaxy clusters within the same sphere, which on its own also corresponds to a ≈1.5σ feature. Interestingly, we find that only 44 out of 15 635 random realizations (i.e., 0.28%) match both anomalies, thus making the Local Universe a ≈3σ environment. We finally compared a mock galaxy catalog with the observed distribution of galaxies in the Local Universe, finding a match to the observed factor of 2 overdensity at ∼16 Mpc as well as the observed 15% underdensity at ∼40 Mpc. Conclusions. Constrained simulations of the Local Universe which reproduce the main features of the local density field open a new window for local field cosmology, where the imprint of the specific density field and the impact on the bias through the observational specific tracers can be investigated in detail.

Uchuu–ν2GC galaxies and AGN: cosmic variance forecasts of high-redshift AGN for JWST, Euclid, and LSST

ABSTRACT Measurements of the luminosity function of active galactic nuclei (AGN) at high redshift (z ≳ 6) are expected to suffer from field-to-field variance, including cosmic and Poisson variances. Future surveys, such as those from the Euclid telescope and JWST, will also be affected by field variance. We use the Uchuu simulation, a state-of-the-art cosmological N-body simulation with 2.1 trillion particles in a volume of 25.7 Gpc3, combined with a semi-analytic galaxy and AGN formation model, to generate the Uchuu–ν2GC catalogue, publicly available, that allows us to investigate the field-to-field variance of the luminosity function of AGN. With this Uchuu–ν2GC model, we quantify the cosmic variance as a function of survey area, AGN luminosity, and redshift. In general, cosmic variance decreases with increasing survey area and decreasing redshift. We find that at z ∼ 6 − 7, the cosmic variance depends weakly on AGN luminosity. This is because the typical mass of dark matter haloes in which AGN reside does not significantly depend on luminosity. Due to the rarity of AGN, Poisson variance dominates the total field-to-field variance, especially for bright AGN. We also examine the effect of parameters related to galaxy formation physics on the field variance. We discuss uncertainties present in the estimation of the faint-end of the AGN luminosity function from recent observations, and extend this to make predictions for the expected number of AGN and their variance for upcoming observations with Euclid, JWST, and the Legacy Survey of Space and Time (LSST).

Study of structural parameters and systemic proper motion of Sextans dwarf spheroidal galaxy with Subaru Hyper Suprime-Cam data

ABSTRACT We use the Subaru Hyper Suprime-Cam (HSC) data to study structural parameters and systemic proper motion of the Sextans dwarf spheroidal galaxy at the heliocentric distance of 86 kpc, which is one of the most important targets for studies of dark matter nature and galaxy formation physics. Thanks to the superb image quality and wide area coverage, the HSC data enable a secure selection of member star candidates based on the colour–magnitude cut, yielding about 10 000 member candidates at magnitudes down to i ∼ 24. We use a likelihood analysis of the two-dimensional distribution of stars to estimate the structural parameters of Sextans taking into account the contamination of foreground halo stars in the Milky Way, and find that the member star distribution is well fitted by an elliptical King profile with ellipticity ϵ ≃ 0.25 and the core and tidal radii of Rc = (368.4 ± 8.5) pc and Rt = (2.54 ± 0.046) kpc, respectively. Then using the two HSC data sets of 2.66 yr time baseline on average, we find the systemic proper motions of Sextans to be (μα, μδ) = (−0.448 ± 0.075, 0.058 ± 0.078) mas yr−1, which are consistent with some of the previous works using the Gaia data of relatively bright member stars in Sextans. Thus, our results give a demonstration that ground-based, large-aperture telescope data that cover a wide solid angle of the sky and have a long time baseline, such as the upcoming data from Legacy Survey of Space and Time (LSST), can be used to study systemic proper motions of dwarf galaxies.

The halo bispectrum as a sensitive probe of massive neutrinos and baryon physics

ABSTRACT The power spectrum has been a workhorse for cosmological studies of large-scale structure. However, the present-day matter distribution is highly non-Gaussian and significant cosmological information is also contained in higher order correlation functions. Meanwhile, baryon physics (particularly active galactic nucleus feedback) has previously been shown to strongly affect the two-point statistics but there has been limited exploration of its effects on higher order functions to date. Here, we use the bahamas suite of cosmological hydrodynamical simulations to explore the effects of baryon physics and massive neutrinos on the halo bispectrum. In contrast to matter clustering which is suppressed by baryon physics, we find that the halo clustering is typically enhanced. The strength of the effect and the scale over which it extends depends on how haloes are selected. On small scales ($k \gtrsim 1$ h Mpc−1, dominated by satellites of groups/clusters), we find that the bispectrum is highly sensitive to the efficiency of star formation and feedback, making it an excellent testing ground for galaxy formation models. We show that the effects of feedback and the effects of massive neutrinos are largely separable (independent of each other) and that massive neutrinos strongly suppress the halo bispectrum on virtually all scales up to the free-streaming length (apart from the smallest scales, where baryon physics dominates). The strong sensitivity of the bispectrum to neutrinos on the largest scales and galaxy formation physics on the smallest scales bodes well for upcoming precision measurements from the next generation of wide-field surveys.

Consistent and simultaneous modelling of galaxy clustering and galaxy–galaxy lensing with subhalo abundance matching

ABSTRACT The spatial distribution of galaxies and their gravitational lensing signal offer complementary tests of galaxy formation physics and cosmology. However, their synergy can only be fully exploited if both probes are modelled accurately and consistently. In this paper, we demonstrate that this can be achieved using an extension of Sub-Halo Abundance Matching extended model (SHAMe), dubbed SHAMe. Specifically, we use mock catalogues built from the TNG300 hydrodynamical simulation to show that SHAMe can simultaneously model the multipoles of the redshift–space galaxy correlation function and galaxy–galaxy lensing, without noticeable bias within the statistical sampling uncertainties of a SDSS volume and on scales r ∈ [0.6 − 30]h−1 Mpc. Modelling the baryonic processes in galaxy–galaxy lensing with a baryonification scheme allows SHAMe’s range of validity to be extended to r ∈ [0.1 − 30] h−1 Mpc. Remarkably, our model achieves this level of precision with just five free parameters beyond those describing the baryonification model. At fixed cosmology, we find that galaxy–galaxy lensing provides a general consistency test but little additional information on galaxy modelling parameters beyond that encoded in the redshift-space multipoles. It does, however, improve constraints if only the projected correlation function is available, as in surveys with only photometric redshifts. We expect SHAMe to have a higher fidelity across a wider range of scales than more traditional methods such as Halo Occupation Distribution modelling. Thus it should provide a significantly more powerful and more robust tool for analysing next-generation large-scale surveys.

Warm dark matter constraints using Milky Way satellite observations and subhalo evolution modeling

Warm dark matter (WDM) can potentially explain small-scale observations that currently challenge the cold dark matter (CDM) model, as warm particles suppress structure formation due to free-streaming effects. Observing small-scale matter distribution provides a valuable way to distinguish between CDM and WDM. In this work, we use observations from the Dark Energy Survey and PanSTARRS1, which observe 270 Milky-Way satellites after completeness corrections. We test WDM models by comparing the number of satellites in the Milky Way with predictions derived from the Semi-Analytical SubHalo Inference ModelIng (SASHIMI) code, which we develop based on the extended Press-Schechter formalism and subhalos' tidal evolution prescription. We robustly rule out WDM with masses lighter than 4.4 keV at 95% confidence level for the Milky-Way halo mass of $10^{12} M_\odot$. The limits are a weak function of the (yet uncertain) Milky-Way halo mass, and vary as $m_{\rm WDM}>3.6$-$5.1$ keV for $(0.6$-$2.0) \times 10^{12} M_\odot$. For the sterile neutrinos that form a subclass of WDM, we obtain the constraints of $m_{\nu_s}>11.6$ keV for the Milky-Way halo mass of $10^{12} M_{\odot}$. These results based on SASHIMI do not rely on any assumptions of galaxy formation physics or are not limited by numerical resolution. The models, therefore, offer a robust and fast way to constrain the WDM models. By applying a satellite forming condition, however, we can rule out the WDM mass lighter than 9.0 keV for the Milky-Way halo mass of $10^{12} M_\odot$.

Constraining galactic baryon cycle using the galaxy stellar-to-halo mass relations

ABSTRACT Galaxies display several well-behaved scaling relations between their properties, such as the star formation rate–stellar mass relation (the main sequence, MS) and the stellar mass–halo mass relation (SHMR). In principle, these scaling relations could imply different star formation histories (SFHs) of galaxies and different constraints on galaxy formation physics. In this paper, we derive the SFHs of galaxies by assuming that they always follow the SHMRs at different redshifts and use an empirical model to constrain key processes in their baryon cycle. It is found that, besides cold accretion due to halo growth, outflow of gas produced by stellar feedback has to be recycled to sustain the derived SFHs of galaxies. The recycled fraction is strongly affected by the baryon fraction in accreted low-mass haloes and the mass loading factor that quantifies the ratio between the galactic outflow rate and star formation rate. Our fiducial model predicts that around 20–60 per cent of outflow is recycled in $\sim 0.5\!-\!4\, \mathrm{Gyr}$, while simulations predict a slightly higher recycle fraction and a lower recycle time. We argue that strong constraints on the baryon cycle process can be obtained from future observation of the circum-galactic medium (CGM) of galaxies, such as the gas cooling rate of CGM. We also find that the implied SFHs from the SHMRs indicate that galaxies stay on the MS only for part of their lifetimes. Our model reproduces the evolution of the mass–metallicity relation as well.

Mesh-free hydrodynamics in pkdgrav3 for galaxy formation simulations

ABSTRACT We extend the state-of-the-art N-body code pkdgrav3 with the inclusion of mesh-free gas hydrodynamics for cosmological simulations. Two new hydrodynamic solvers have been implemented, the mesh-less finite volume and mesh-less finite mass methods. The solvers manifestly conserve mass, momentum, and energy, and have been validated with a wide range of standard test simulations, including cosmological simulations. We also describe improvements to pkdgrav3 that have been implemented for performing hydrodynamic simulations. These changes have been made with efficiency and modularity in mind, and provide a solid base for the implementation of the required modules for galaxy formation and evolution physics and future porting to GPUs. The code is released in a public repository, together with the documentation, and all the test simulations presented in this work.

Diffstar: a fully parametric physical model for galaxy assembly history

ABSTRACT We present Diffstar , a smooth parametric model for the in situ star formation history (SFH) of galaxies. The Diffstar model is distinct from traditional SFH models because it is parametrized directly in terms of basic features of galaxy formation physics. Diffstar includes ingredients for: the halo mass assembly history; the accretion of gas into the dark matter halo; the fraction of gas that is eventually transformed into stars, ϵms; the time-scale over which this transformation occurs, τcons; and the possibility that some galaxies will experience a quenching event at time tq, and may subsequently experience rejuvenated star formation. We show that our model is sufficiently flexible to describe the average stellar mass histories of galaxies in both the IllustrisTNG (TNG) and UniverseMachine (UM) simulations with an accuracy of ∼0.1 dex across most of cosmic time. We use Diffstar to compare TNG to UM in common physical terms, finding that: (i) star formation in UM is less efficient and burstier relative to TNG; (ii) UM galaxies have longer gas consumption time-scales, relative to TNG; (iii) rejuvenated star formation is ubiquitous in UM, whereas quenched TNG galaxies rarely experience sustained rejuvenation; and (iv) in both simulations, the distributions of ϵms, τcons, and tq share a common characteristic dependence upon halo mass, and present significant correlations with halo assembly history. We conclude with a discussion of how Diffstar can be used in future applications to fit the SEDs of individual observed galaxies, as well as in forward-modelling applications that populate cosmological simulations with synthetic galaxies.

Unveiling the main sequence of galaxies at z ≥ 5 with the JWST: predictions from simulations

ABSTRACT We use two independent galaxy-formation simulations, flares, a cosmological hydrodynamical simulation, and shark, a semi-analytic model, to explore how well the JWST will be able to uncover the existence and parameters of the star-forming main sequence (SFS) at z = 5 → 10, i.e. shape, scatter, normalization. Using two independent simulations allows us to isolate predictions (e.g. stellar mass, star formation rate, SFR, luminosity functions) that are robust to or highly dependent on the implementation of the physics of galaxy formation. Both simulations predict that JWST can observe ≥70–90 per cent (for shark and flares, respectively) of galaxies up to z ∼ 10 (down to stellar masses of ${\approx}10^{8.3}\rm M_{\odot }$ and SFRs of ${\approx}10^{0.5}{\rm M}_{\odot }\,{\rm yr}^{-1}$) in modest integration times and given current proposed survey areas (e.g. the Web COSMOS 0.6 deg2) to accurately constrain the parameters of the SFS. Although both simulations predict qualitatively similar distributions of stellar mass and SFR. There are important quantitative differences, such as the abundance of massive, star-forming galaxies with flares predicting a higher abundance than shark; the early onset of quenching as a result of black hole growth in flares (at z ≈ 8), not seen in shark until much lower redshifts; and the implementation of synthetic photometry with flares predicting more JWST-detected galaxies (∼90 per cent) than shark (∼70 per cent) at z = 10. JWST observations will distinguish between these models, leading to a significant improvement upon our understanding of the formation of the very first galaxies.

Spin transfer from dark matter to gas during halo formation

ABSTRACT In the protogalactic density field, diffuse gas and collision-less cold dark matter (DM) are often assumed sufficiently mixed that both components experience identical tidal torques. However, haloes in cosmological simulations consistently end up with a higher specific angular momentum (sAM) in gas, even in simulations without radiative cooling and galaxy formation physics. We refine this result by analysing the spin distributions of gas and DM in ∼50 000 well-resolved haloes in a non-radiative cosmological simulation from the SURFS suite. The sAM of the halo gas on average ends up ∼40 per cent above that of the DM. This can be pinned down to an excess AM in the inner halo (<50 per cent virial radius), paralleled by a more coherent rotation pattern in the gas. We uncover the leading driver for this AM difference through a series of control simulations of a collapsing ellipsoidal top-hat, where gas and DM are initially well mixed. These runs reveal that the pressurized inner gas shells collapse more slowly, causing the DM ellipsoid to spin ahead of the gas ellipsoid. The arising torque generally transfers AM from the DM to the gas. The amount of AM transferred via this mode depends on the initial spin, the initial axes ratios, and the collapse factor. These quantities can be combined in a single dimensionless parameter, which robustly predicts the AM transfer of the ellipsoidal collapse. This simplistic model can quantitatively explain the average AM excess of the gas found in the more complex non-radiative cosmological simulation.

Hot-mode accretion and the physics of thin-disc galaxy formation

ABSTRACT We use FIRE simulations to study disc formation in z ∼ 0, Milky Way-mass galaxies, and conclude that a key ingredient for the formation of thin stellar discs is the ability for accreting gas to develop an aligned angular momentum distribution via internal cancellation prior to joining the galaxy. Among galaxies with a high fraction ($\gt 70{{\ \rm per\ cent}}$) of their young stars in a thin disc (h/R ∼ 0.1), we find that: (i) hot, virial-temperature gas dominates the inflowing gas mass on halo scales (≳20 kpc), with radiative losses offset by compression heating; (ii) this hot accretion proceeds until angular momentum support slows inward motion, at which point the gas cools to $\lesssim 10^4\, {\rm K}$; (iii) prior to cooling, the accreting gas develops an angular momentum distribution that is aligned with the galaxy disc, and while cooling transitions from a quasi-spherical spatial configuration to a more-flattened, disc-like configuration. We show that the existence of this ‘rotating cooling flow’ accretion mode is strongly correlated with the fraction of stars forming in a thin disc, using a sample of 17 z ∼ 0 galaxies spanning a halo mass range of 1010.5 M⊙ ≲ Mh ≲ 1012 M⊙ and stellar mass range of 108 M⊙ ≲ M⋆ ≲ 1011 M⊙. Notably, galaxies with a thick disc or irregular morphology do not undergo significant angular momentum alignment of gas prior to accretion and show no correspondence between halo gas cooling and flattening. Our results suggest that rotating cooling flows (or, more generally, rotating subsonic flows) that become coherent and angular momentum-supported prior to accretion on to the galaxy are likely a necessary condition for the formation of thin, star-forming disc galaxies in a ΛCDM universe.

The local PNG bias of neutral Hydrogen, HI

We use separate universe simulations with the IllustrisTNG galaxy formation model to predict the local PNG bias parameters bΦ and bΦδ of atomic neutral hydrogen, HI. These parameters and their relation to the linear density bias parameter b 1 play a key role in observational constraints of the local PNG parameter f NL using the HI power spectrum and bispectrum. Our results show that the popular calculation based on the universality of the halo mass function overpredicts the bΦ (b 1) and bΦδ (b 1) relations measured in the simulations. In particular, our results show that at z ≲ 1 the HI power spectrum is more sensitive to f NL compared to previously thought (bΦ is more negative), but is less sensitive at other epochs (bΦ is less positive). We discuss how this can be explained by the competition of physical effects such as that large-scale gravitational potentials with local PNG (i) accelerate the conversion of hydrogen to heavy elements by star formation, (ii) enhance the effects of baryonic feedback that eject the gas to regions more exposed to ionizing radiation, and (iii) promote the formation of denser structures that shield the HI more efficiently. Our numerical results can be used to revise existing forecast studies on f NL using 21 cm line-intensity mapping data. Despite this first step towards predictions for the local PNG bias parameters of HI, we emphasize that more work is needed to assess their sensitivity on the assumed galaxy formation physics and HI modeling strategy.

RAMSES-RTZ: non-equilibrium metal chemistry and cooling coupled to on-the-fly radiation hydrodynamics

ABSTRACT Emission and absorption lines from elements heavier than helium (metals) represent one of our strongest probes of galaxy formation physics across nearly all redshifts accessible to observations. The vast majority of simulations that model these metal lines often assume either collisional or photoionization equilibrium, or a combination of the two. For the few simulations that have relaxed these assumptions, a redshift-dependent meta-galactic UV background or fixed spectrum is often used in the non-equilibrium photoionization calculation, which is unlikely to be accurate in the interstellar medium where the gas can self-shield, as well as in the high-redshift circumgalactic medium, where locally emitted radiation may dominate over the UV background. In this work, we relax this final assumption by coupling the ionization states of individual metals to the radiation hydrodynamics solver present in ramses-rt. Our chemical network follows radiative recombination, dielectronic recombination, collisional ionization, photoionization, and charge transfer, and we use the ionization states to compute non-equilibrium optically thin metal-line cooling. The fiducial model solves for the ionization states of C, N, O, Mg, Si, S, Fe, and Ne in addition to H, He, and H2, but can be easily extended for other ions. We provide interfaces to two different ODE solvers that are competitive in both speed and accuracy. The code has been benchmarked across a variety of gas conditions to reproduce results from cloudy when equilibrium is reached. We show an example isolated galaxy simulation with on-the-fly radiative transfer that demonstrates the utility of our code for translating between simulations and observations without the use of idealized photoionization models.

UNITSIM-Galaxies: data release and clustering of emission-line galaxies

ABSTRACT New surveys such as European Space Agencys (ESA’s) Euclid mission are planned to map with unprecedented precision the large-scale structure of the Universe by measuring the 3D positions of tens of millions of galaxies. It is necessary to develop theoretically modelled galaxy catalogues to estimate the expected performance and to optimize the analysis strategy of these surveys. We populate two pairs of (1 h−1 Gpc)3 volume dark matter-only simulations from the UNIT project with galaxies using the Semi-Analytic Galaxy Evolution semi-analytic model of galaxy formation, coupled to the photoionization model get_emlines to estimate their H α emission. These catalogues represent a unique suite that includes galaxy formation physics and – thanks to the fixed-pair technique used – an effective volume of $\sim\!(5\,h^{-1}\,\rm {Gpc})^3$, which is several times larger than the Euclid survey. We present the performance of these data and create five additional emission-line galaxy (ELG) catalogues by applying a dust-attenuation model as well as adjusting the flux threshold as a function of redshift in order to reproduce Euclid-forecast dN/dz values. As a first application, we study the abundance and clustering of those model H α ELGs: for scales greater than ∼5 h−1 Mpc, we find a scale-independent bias with a value of b ∼ 1 at redshift z ∼ 0.5, that can increase nearly linearly to b ∼ 4 at z ∼ 2, depending on the ELG catalogue. Model galaxy properties, including their emission-line fluxes (with and without dust extinction) are publicly available.

How cosmological merger histories shape the diversity of stellar haloes

ABSTRACT We introduce and apply a new approach to probe the response of galactic stellar haloes to the interplay between cosmological merger histories and galaxy formation physics. We perform dark matter-only, zoomed simulations of two Milky Way-mass hosts and make targeted, controlled changes to their cosmological histories using the genetic modification technique. Populating each history’s stellar halo with a semi-empirical, particle tagging approach then enables a controlled study, with all instances converging to the same large-scale structure, dynamical and stellar mass at z = 0 as their reference. These related merger scenarios alone generate an extended spread in stellar halo mass fractions (1.5 dex) comparable to the observed population, with the largest scatter achieved by growing late (z ≤ 1) major mergers that spread out existing stars to create massive, in-situ dominated stellar haloes. Increasing a last major merger at z ∼ 2 brings more accreted stars into the inner regions, resulting in smaller scatter in the outskirts which are predominantly built by subsequent minor events. Exploiting the flexibility of our semi-empirical approach, we show that the diversity of stellar halo masses across scenarios is reduced by allowing shallower slopes in the stellar mass–halo mass relation for dwarf galaxies, while it remains conserved when central stars are born with hotter kinematics across cosmic time. The merger-dependent diversity of stellar haloes thus responds distinctly to assumptions in modelling the central and dwarf galaxies respectively, opening exciting prospects to constrain star formation and feedback at different galactic mass-scales with the coming generation of deep, photometric observatories.