Search for: All records

We investigate galaxy sizes at redshift $z\gtrsim 6$ with the cosmological radiation-magnetohydrodynamic simulation suite thesan(-hr). These simulations simultaneously capture reionization of the large-scale intergalactic medium and resolved galaxy properties. The intrinsic sizes ($r^{\ast }_{1/2}$) of simulated galaxies increase moderately with stellar mass at $M_{\ast } \lesssim 10^{8}{\, \rm M_\odot}$ and decrease fast at larger masses, resulting in a hump feature at $M_{\ast }\sim 10^{8}{\, \rm M_\odot}$ that is insensitive to redshift. Low-mass galaxies are in the initial phase of size growth and are better described by a spherical shell model with feedback-driven outflows competing with the cold inflowing gas streams. In contrast, massive galaxies fit better with the disc formation model. They generally experience a phase of rapid compaction and gas depletion, likely driven by internal disc instability rather than external processes. We identify four compact quenched galaxies in the $(95.5\, {\rm cMpc})^{3}$ volume of thesan-1 at $z\simeq 6$ and their quenching follows reaching a characteristic stellar surface density akin to the massive compact galaxies at cosmic noon. Compared to observations, we find that the median ultraviolet effective radius ($R^{\rm UV}_{\rm eff}$) of simulated galaxies is at least three times larger than the observed ones at $M_{\ast }\lesssim 10^{9}{\, \rm M_\odot}$ or $M_{\rm UV}\gtrsim -20$ at $6 \lesssim z \lesssim 10$. The population of compact galaxies ($R^{\rm UV}_{\rm eff}\lesssim 300\, {\rm pc}$) galaxies at $M_{\ast }\sim 10^{8}{\, \rm M_\odot}$ is missing in our simulations. This inconsistency persists across many other cosmological simulations with different galaxy formation models and demonstrates the potential of using galaxy morphology to constrain physics of galaxy formation at high redshifts.

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 (Reff) 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 (zreion) with the time it takes to reach a given bubble size and conclude that zreion 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 zreion and size statistics weakens due to the time delay between the onset of reionization and the expansion of large bubbles, particularly at high redshifts.

The distribution of offsets between the brightest cluster galaxies of galaxy clusters and the centroid of their dark matter distributions is a promising probe of the underlying dark matter physics. In particular, since this distribution is sensitive to the shape of the potential in galaxy cluster cores, it constitutes a test of dark matter self-interaction on the largest mass scales in the universe. We examine these offsets in three suites of modern cosmological simulations; IllustrisTNG, MillenniumTNG and BAHAMAS. For clusters above , we examine the dependence of the offset distribution on gravitational softening length, the method used to identify centroids, redshift, mass, baryonic physics, and establish the stability of our results with respect to various nuisance parameter choices. We find that offsets are overwhelmingly measured to be smaller than the minimum converged length scale in each simulation, with a median offset of in the highest resolution simulation considered, TNG300-1, which uses a gravitational softening length of . We also find that centroids identified via source extraction on smoothed dark matter and stellar particle data are consistent with the potential minimum, but that observationally relevant methods sensitive to cluster strong gravitational lensing scales, or those using the the “light traces mass” approach, in this context meaning gas is used as a tracer for the potential, can overestimate offsets by factors of$\sim 10$and$\sim 30$, respectively. This has the potential to reduce tensions with existing offset measurements which have served as evidence for a nonzero dark matter self-interaction cross section.

It has been claimed that traditional models struggle to explain the tentative detection of the 21 cm absorption trough centered atz∼ 17 measured by the EDGES collaboration. On the other hand, it has been shown that the EDGES results are consistent with an extrapolation of a declining UV luminosity density, following a simple power law of deep Hubble Space Telescope observations of 4 <z< 9 galaxies. We here explore the conditions by which the EDGES detection is consistent with current reionization and post-reionization observations, including the neutral hydrogen fraction atz∼ 6–8, Thomson-scattering optical depth, and ionizing emissivity atz∼ 5. By coupling a physically motivated source model derived from radiative transfer hydrodynamic simulations of reionization to a Markov Chain Monte Carlo sampler, we find that it is entirely possible to reconcile existing high-redshift (cosmic dawn) and low-redshift (reionization) constraints. In particular, we find that high contributions from low-mass halos along with high photon escape fractions are required to simultaneously reproduce cosmic dawn and reionization constraints. Our analysis further confirms that low-mass galaxies produce a flatter emissivity evolution, which leads to an earlier onset of reionization with a gradual and longer duration, resulting in a higher optical depth. While the models dominated by faint galaxies successfully reproduce the measured globally averaged quantities over the first one billion years, they underestimate the late redshift-instantaneous measurements in efficiently star-forming and massive systems. We show that our (simple) physically motivated semianalytical prescription produces results that are consistent with the (sophisticated) state-of-the-artTHESANradiation-magnetohydrodynamic simulation of the reionization.

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.

JWST observations have revealed a population of galaxies bright enough that potentially challenge standard galaxy formation models in the Λ cold dark matter (ΛCDM) cosmology. Using a minimal empirical framework, we investigate the influence of variability on the rest-frame ultra-violet (UV) luminosity function of galaxies at z ≥ 9. Our study differentiates between the median UV radiation yield and the variability of UV luminosities of galaxies at a fixed dark matter halo mass. We primarily focus on the latter effect, which depends on halo assembly and galaxy formation processes and can significantly increase the abundance of UV-bright galaxies due to the upscatter of galaxies in lower-mass haloes. We find that a relatively low level of variability, σUV ≈ 0.75 mag, matches the observational constraints at z ≈ 9. However, increasingly larger σUV is necessary when moving to higher redshifts, reaching $\sigma _{\rm UV} \approx 2.0\, (2.5)\, {\rm mag}$ at z ≈ 12 (16). This implied variability is consistent with expectations of physical processes in high-redshift galaxies such as bursty star formation and dust clearance during strong feedback cycles. Photometric constraints from JWST at z ≳ 9 therefore can be reconciled with a standard ΛCDM-based galaxy formation model calibrated at lower redshifts without the need for adjustments to the median UV radiation yield.

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.

Using high-resolution cosmological radiation-hydrodynamic (RHD) simulations (thesan-hr), we explore the impact of alternative dark matter (altDM) models on galaxies during the Epoch of Reionization. The simulations adopt the IllustrisTNG galaxy formation model. We focus on altDM models that exhibit small-scale suppression of the matter power spectrum, namely warm dark matter (WDM), fuzzy dark matter (FDM), and interacting dark matter (IDM) with strong dark acoustic oscillations (sDAO). In altDM scenarios, both the halo mass functions and the ultraviolet luminosity functions at z ≳ 6 are suppressed at the low-mass/faint end, leading to delayed global star formation and reionization histories. However, strong non-linear effects enable altDM models to ‘catch up’ with cold dark matter (CDM) in terms of star formation and reionization. The specific star formation rates are enhanced in halos below the half-power mass in altDM models. This enhancement coincides with increased gas abundance, reduced gas depletion times, more compact galaxy sizes, and steeper metallicity gradients at the outskirts of the galaxies. These changes in galaxy properties can help disentangle altDM signatures from a range of astrophysical uncertainties. Meanwhile, it is the first time that altDM models have been studied in RHD simulations of galaxy formation. We uncover significant systematic uncertainties in reionization assumptions on the faint-end luminosity function. This underscores the necessity of accurately modeling the small-scale morphology of reionization in making predictions for the low-mass galaxy population. Upcoming James Webb Space Telescope imaging surveys of deep lensed fields hold potential for uncovering the faint low-mass galaxy population, which could provide constraints on altDM models.

The feedback loop between the galaxies producing the background radiation field for reionization and their growth is crucial, particularly for low-mass haloes. Despite this, the vast majority of galaxy formation studies employ a spatially uniform, time-varying reionizing background, with the majority of reionization studies employing galaxy formation models only required to work at high redshift. This paper uses the well-studied TNG galaxy formation model, calibrated at low redshift, coupled to the arepo-rt code, to self-consistently solve the coupled problems of galaxy evolution and reionization, evaluating the impact of patchy (and slow) reionization on early galaxies. thesan-hr is an extension of the thesan project to higher resolution (a factor of 50 increase, with a baryonic mass of mb ≈ 104 M⊙), to additionally enable the study of ‘mini-haloes’ with virial temperatures Tvir < 104 K. Comparing the self-consistent model to a uniform UV background, we show that galaxies in thesan-hr are predicted to be larger in physical extent (by a factor ∼2), less metal enriched (by ∼0.2 dex), and less abundant (by a factor ∼10 at M1500 =   − 10) by z = 5. We show that differences in star formation and enrichment patterns lead to significantly different predictions for star formation in low mass haloes, low-metallicity star formation, and even the occupation fraction of haloes. We posit that cosmological galaxy formation simulations aiming to study early galaxy formation (z ≳ 3) must employ a spatially inhomogeneous UV background to accurately reproduce galaxy properties.

Upcoming large galaxy surveys will subject the standard cosmological model, Lambda Cold Dark Matter, to new precision tests. These can be tightened considerably if theoretical models of galaxy formation are available that can predict galaxy clustering and galaxy–galaxy lensing on the full range of measurable scales, throughout volumes as large as those of the surveys, and with sufficient flexibility that uncertain aspects of the underlying astrophysics can be marginalized over. This, in particular, requires mock galaxy catalogues in large cosmological volumes that can be directly compared to observation, and can be optimized empirically by Monte Carlo Markov Chains or other similar schemes, thus eliminating or estimating parameters related to galaxy formation when constraining cosmology. Semi-analytic galaxy formation methods implemented on top of cosmological dark matter simulations offer a computationally efficient approach to construct physically based and flexibly parametrized galaxy formation models, and as such they are more potent than still faster, but purely empirical models. Here, we introduce an updated methodology for the semi-analytic L-Galaxies code, allowing it to be applied to simulations of the new MillenniumTNG project, producing galaxies directly on fully continuous past lightcones, potentially over the full sky, out to high redshift, and for all galaxies more massive than $\sim 10^8\, {\rm M}_\odot$. We investigate the numerical convergence of the resulting predictions, and study the projected galaxy clustering signals of different samples. The new methodology can be viewed as an important step towards more faithful forward-modelling of observational data, helping to reduce systematic distortions in the comparison of theory to observations.