Search

The Voids-within-Voids-within-Voids project used dark-matter-only (DMO) simulations to study the abundance and structure of dark matter (DM) haloes over the full mass range populated in the standard Lambda cold dark matter cosmology. Here, we explore how baryonic effects modify these results for |$z=0$| halo masses in the range |$10^4$| – |$10^7~\mathrm{M_\odot }$|⁠ , below the threshold for galaxy formation. Our main study focuses on three simulations from identical initial conditions at |$z=127$|⁠ , one following DMO, one including non-radiative gas, and one additionally including the baryonic physics relevant in this halo mass range (cooling and photoheating). In the non-radiative simulation, above |$10^{5.5}~\mathrm{M_\odot }$|⁠ , halo abundance and internal structure are very similar to the DMO simulation, and the baryon to DM ratio is everywhere close to the cosmic value. At lower mass, this ratio drops and haloes are less concentrated and less massive in the non-radiative case. Test simulations at higher resolution show this to be mainly a resolution effect; the expected drop in baryon content due to residual pressure effects only becomes substantial for |$z=0$| haloes below |${\sim}10^{2.7}~\mathrm{M_\odot }$|⁠. However, gas is heated by reionization at |$z=6$| in our 'full physics' run, and this results in almost complete expulsion of gas from all haloes in our simulated mass range. This suppresses the halo mass function by |${\sim}30{{\ \rm per\ cent}}$|⁠ , lowers halo concentration, and consequently weakens the DM annihilation signal by |${\sim}40{-}60{{\ \rm per\ cent}}$|⁠. [ABSTRACT FROM AUTHOR]

We present an extended suite of the Auriga cosmological gravo-magnetohydrodynamical 'zoom-in' simulations of 40 Milky Way-mass haloes and 26 dwarf galaxy–mass haloes run with the moving-mesh code arepo. Auriga adopts the Lambda cold dark matter cosmogony and includes a comprehensive galaxy formation physics model following the coupled cosmic evolution of dark matter, gas, stars, and supermassive black holes which has been shown to produce numerically well-converged galaxy properties for Milky Way-mass systems. We describe the first public data release of this augmented suite of Auriga simulations, which includes raw snapshots, group catalogues, merger trees, initial conditions, and supplementary data, as well as public analysis tools with worked examples of how to use the data. To demonstrate the value and robustness of the simulation predictions, we analyse a series of low-redshift global properties that compare well with many observed scaling relations, such as the Tully–Fisher relation, the star-forming (SF) main sequence, and H  i gas fraction/disc thickness. Finally, we show that SF gas discs appear to build rotation and velocity dispersion rapidly for |$z\gtrsim 3$| before they 'settle' into ever-increasing rotation-dispersion ratios (⁠|$V/\sigma$|⁠). This evolution appears to be in rough agreement with some kinematic measurements from H |$\alpha$| observations, and demonstrates an application of how to utilize the released data. [ABSTRACT FROM AUTHOR]

We present a machine-learning-based model for the total density profiles of subhaloes with masses |$M \gtrsim 7\times 10^8\, h^{-1}{\rm M}_\odot$| in the IllustrisTNG100 simulation. The model is based on an interpretable variational encoder (IVE) which returns the independent factors of variation in the density profiles within a low-dimensional representation, as well as the predictions for the density profiles themselves. The IVE returns accurate and unbiased predictions on all radial ranges, including the outer region profile where the subhaloes experience tidal stripping; here its fit accuracy exceeds that of the commonly used Einasto profile. The IVE discovers three independent degrees of freedom in the profiles, which can be interpreted in terms of the formation history of the subhaloes. In addition to the two parameters controlling the normalization and inner shape of the profile, the IVE discovers a third parameter that accounts for the impact of tidal stripping on to the subhalo outer profile; this parameter is sensitive to the mass loss experienced by the subhalo after its infall on to its parent halo. Baryonic physics in the IllustrisTNG galaxy formation model does not impact the number of degrees of freedom identified in the profile compared to the pure dark matter expectations, nor their physical interpretation. Our newly proposed profile fit can be used in strong lensing analyses or other observational studies which aim to constrain cosmology from small-scale structures. [ABSTRACT FROM AUTHOR]

An important characteristic of cosmic hydrogen reionization is the growth of ionized gas bubbles surrounding early luminous objects. Ionized bubble sizes are beginning to be probed using Lyman α emission from high-redshift galaxies, and will also be probed by upcoming 21 cm maps. We present results from a study of bubble sizes using the state-of-the-art thesan radiation-hydrodynamics simulation suite, which self-consistently models radiation transport and realistic galaxy formation. We employ the mean free path method and track the evolution of the effective ionized bubble size at each point (R eff) throughout the Epoch of Reionization. We show that there is a slow growth period for regions ionized early, but a rapid 'flash ionization' process for regions ionized later as they immediately enter a large, pre-existing bubble. We also find that bright sources are preferentially in larger bubbles, and find consistency with recent observational constraints at z ≳ 9, but tension with idealized Lyman α damping-wing models at z ≈ 7. We find that high-overdensity regions have larger characteristic bubble sizes, but the correlation decreases as reionization progresses, likely due to runaway formation of large percolated bubbles. Finally, we compare the redshift at which a region transitions from neutral to ionized (z reion) with the time it takes to reach a given bubble size and conclude that z reion is a reasonable local probe of small-scale bubble size statistics (⁠|$R_\text{eff} \lesssim 1\, \rm {cMpc}$|⁠). However, for larger bubbles, the correspondence between z reion and size statistics weakens due to the time delay between the onset of reionization and the expansion of large bubbles, particularly at high redshifts. [ABSTRACT FROM AUTHOR]

Cosmological simulations serve as invaluable tools for understanding the Universe. However, the technical complexity and substantial computational resources required to generate such simulations often limit their accessibility within the broader research community. Notable exceptions exist, but most are not suited for simultaneously studying the physics of galaxy formation and cosmic reionization during the first billion years of cosmic history. This is especially relevant now that a fleet of advanced observatories (e.g. James Webb Space Telescope, Nancy Grace Roman Space Telescope , SPHEREx, ELT, SKA) will soon provide an holistic picture of this defining epoch. To bridge this gap, we publicly release all simulation outputs and post-processing products generated within the thesan simulation project at www.thesan-project.com. This project focuses on the z ≥ 5.5 Universe, combining a radiation-hydrodynamics solver (arepo-rt), a well-tested galaxy formation model (IllustrisTNG) and cosmic dust physics to provide a comprehensive view of the Epoch of Reionization. The thesan suite includes 16 distinct simulations, each varying in volume, resolution, and underlying physical models. This paper outlines the unique features of these new simulations, the production and detailed format of the wide range of derived data products, and the process for data retrieval. Finally, as a case study, we compare our simulation data with a number of recent observations from the James Webb Space Telescope , affirming the accuracy and applicability of thesan. The examples also serve as prototypes for how to utilize the released data set to perform comparisons between predictions and observations. [ABSTRACT FROM AUTHOR]

The properties of galaxies in low-density regions of the universe suggest an interplay between galaxy formation and environment. However, the specific reason why this particular large-scale environment influences the evolution of galaxies remains unclear. This paper examines the properties and evolutionary paths of galaxies within cosmic voids using the Illustris TNG300 simulation. The population of void galaxies at z = 0 has a higher star formation rate, a smaller stellar-to-halo-mass ratio, higher gas metallicity, and lower stellar metallicity in comparison with non-void galaxies at fixed stellar mass. Our analysis shows that these differences are mainly due to the characteristics of galaxies classified as satellites, for which the largest differences between void and non-void samples are found. Although the mean number of mergers is similar between void and non-void samples at a fixed stellar mass, void galaxies tend to experience mergers at later times, resulting in a more recent accumulation of accreted stellar mass. While the mean net accreted mass is comparable for high-mass galaxies, low-mass void galaxies tend to exhibit higher fractions of accreted stars than non-void galaxies. This finding challenges the common notion that void galaxies predominantly experience growth with infrequent mergers or interactions. With this paper, we also publicly release our void catalogue as part of the IllustrisTNG project. [ABSTRACT FROM AUTHOR]

We investigate the growth of stellar discs in Milky Way-mass galaxies using the magnetohydrodynamical simulations of the Auriga Project in a full cosmological context. We focus on the gas accretion process along the discs, calculating the net, infall and outflow rates as a function of galactocentric distance, and investigate the relation between them and the star formation activity. The stellar distributions of around 70 per cent of the simulated galaxies exhibit an 'inside–out' pattern, with older (younger) stellar populations preferentially located in the inner (outer) disc regions. In all cases, we find a very tight correlation between the infall, outflow, and net accretion rates, as well as between these three quantities and the star formation rate. This is because the amount of gas which is ultimately available for star formation in each radial ring depends not only on the infall rates, but also on the amount of gas leaving the disc in outflows, which directly relates to the local star formation level. Therefore, any of these rates can be used to identify galaxies with inside–out growth. For these galaxies, the correlation between the dominant times of accretion/star formation and disc radius is well fitted by a linear function. We also find that, when averaged over galaxies with formation histories similar to the Milky Way, the simulated accretion rates show a similar evolution (both temporally and radially integrated) to the usual accretion prescriptions used in chemical evolution models, although some major differences arise at early times and in the inner disc regions. [ABSTRACT FROM AUTHOR]

We take advantage of the unprecedented dynamical range provided by the 'Cosmic-Zoom' project to study the mass accretion history (MAH) of present-day dark matter haloes over the entire mass range present in the Lambda cold dark matter paradigm when the dark matter is made of weakly interacting massive particles of mass 100 GeV. In particular, we complement previous studies by exploring the MAHs of haloes with mass from |$10^8\ h^{-1}\,\mathrm{{\rm M}_{\odot }}$| down to Earth mass, |$10^{-6}\ h^{-1}\,\mathrm{{\rm M}_{\odot }}$|⁠. The formation redshift of low-mass haloes anticorrelates weakly with mass, peaking at z = 3 for haloes of mass |$10^{-4}\ h^{-1}\,\mathrm{{\rm M}_{\odot }}$|⁠. Even lower masses are affected by the free-streaming cut-off in the primordial spectrum of density fluctuations and form at lower redshift. We compare MAHs in our simulations with predictions from two analytical models based on the extended Press–Schechter theory (EPS), and three empirical models derived by fitting and extrapolating either results from cosmological N -body simulations or Monte Carlo realizations of halo growth. All models fit our simulations reasonably well over the mass range for which they were calibrated. While the empirical models match better for more massive haloes, |$M\gt 10^{10}\ h^{-1}\,\mathrm{{\rm M}_{\odot }}$|⁠ , the analytical models do better when extrapolated down to Earth mass. At the higher masses, we explore the correlation between local environment density and MAH, finding that biases are relatively weak, with typical MAHs for haloes in extremely low-density and in typical regions differing by less than 20 per cent at high redshift. If this result can be extrapolated to lower halo masses, we conclude that EPS theory predicts the hierarchical build up of dark matter haloes quite well over the entire halo mass range. [ABSTRACT FROM AUTHOR]

We analyse the internal structure and dynamics of cosmic-web filaments connecting massive high- z haloes. Our analysis is based on a high-resolution arepo cosmological simulation zooming-in on three Mpc-scale filaments feeding three massive haloes of |$\sim 10^{12}\, \text{M}_\odot$| at z ∼ 4, embedded in a large-scale sheet. Each filament is surrounded by a cylindrical accretion shock of radius |$r_{\rm shock} \sim 50 \, {\rm kpc}$|⁠. The post-shock gas is in virial equilibrium within the potential well set by an isothermal dark-matter filament. The filament line-mass is |$\sim 9\times 10^8\, \text{M}_\odot \, {\rm kpc}^{-1}$|⁠ , the gas fraction within r shock is the universal baryon fraction, and the virial temperature is ∼7 × 105 K. These all match expectations from analytical models for filament properties as a function of halo mass and redshift. The filament cross-section has three radial zones. In the outer 'thermal' (T) zone, |$r \ge 0.65 \, r_{\rm shock}$|⁠ , inward gravity, and ram-pressure forces are overbalanced by outward thermal pressure forces, decelerating the inflowing gas and expanding the shock outwards. In the intermediate 'vortex' (V) zone, 0.25 ≤ r / r shock ≤ 0.65, the velocity field is dominated by a quadrupolar vortex structure due to offset inflow along the sheet through the post-shock gas. The outward force is dominated by centrifugal forces associated with these vortices, with additional contributions from global rotation and thermal pressure. Shear and turbulent forces associated with the vortices act inwards. The inner 'stream' (S) zone, |$r \lt 0.25 \, r_{\rm shock}$|⁠ , is a dense isothermal core, |$T\sim 3 \times 10^4 \, {\rm K}$| and |$n_{\rm H}\sim 0.01 \, {\rm cm^{-3}}$|⁠ , defining the cold streams that feed galaxies. The core is formed by an isobaric cooling flow and is associated with a decrease in outward forces, though exhibiting both inflows and outflows. [ABSTRACT FROM AUTHOR]

Cosmic rays (CRs) are dynamically important in the evolution of galaxies by regulating star formation and powering galactic outflows. However, to what extent CRs regulate galaxy formation depends on the coupling strength of CRs with the ambient plasma and the effective CR transport speed. Moreover, both properties sensitively depend on the CR momentum, which is largely unexplored in three-dimensional hydrodynamical simulations. We perform magnetohydrodynamical simulations of entire galaxies with masses ranging from 1010 to |$10^{12}\, \mathrm{M}_\odot$| and compare dynamically coupled CRs in the grey approximation with a spectrally resolved model that includes CR momenta from |$0.1\, \mathrm{GeV} c^{-1}$| to |$100\, \mathrm{TeV} c^{-1}$|⁠. We find that hadronic cooling of CRs dominates over Alfvén cooling, with the latter emulating CR losses as a result of streaming of CRs down their pressure gradient. While star formation rates and galaxy morphologies are only mildly affected by the spectral CR modelling, mass loading factors of galactic outflows can differ by up to a factor of 4 in dwarf galaxies. All simulated low-mass haloes (M  = 1010, 1011, and |$3\times 10^{11}\, \mathrm{M}_\odot$|⁠) drive strong outflows, where CR transport is temporally dominated by advection. In contrast, the Milky Way-mass galaxy with |$M=10^{12}\, \mathrm{M}_\odot$| does not drive sustained outflows, so that CR transport is entirely dominated by diffusion. The effective energy weighted diffusion coefficients vary by two orders of magnitude from the canonical energy-weighted values of |$\langle {D}\rangle _{e_\mathrm{cr}}\sim 10^{28}\, \mathrm{cm^2\, s^{-1}}$| in the disc up to |$3\times 10^{29}\, \mathrm{cm^2\, s^{-1}}$| in the circumgalactic medium, where we observe substantial temperature and CR pressure differences between our grey and spectral CR models. [ABSTRACT FROM AUTHOR]