skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 5:00 PM ET until 11:00 PM ET on Friday, June 21 due to maintenance. We apologize for the inconvenience.

Title: The ASTRID simulation: galaxy formation and reionization

We introduce the Astrid  simulation, a large-scale cosmological hydrodynamic simulation in a $250 \, h^{-1}\mathrm{Mpc}$ box with 2 × 55003 particles. Astrid contains a large number of high redshift galaxies, which can be compared to future survey data, and resolves galaxies in haloes more massive than $2\times 10^9 \, \mathrm{M}_{\odot }$. Astrid  has been run from z = 99 to 3. As a particular focus is modelling the high redshift Universe, it contains models for inhomogeneous hydrogen and helium reionization, baryon relative velocities and massive neutrinos, as well as supernova and AGN feedback. The black hole model includes mergers driven by dynamical friction rather than repositioning. We briefly summarize the implemented models, and the technical choices we took when developing the simulation code. We validate the model, showing good agreement with observed ultraviolet luminosity functions, galaxy stellar mass functions and specific star formation rates (SFRs). We show that the redshift at which a given galaxy underwent hydrogen reionization has a large effect on the halo gas fraction. Finally, at z = 6, haloes with $M \sim 2\times 10^9 \, \mathrm{M}_{\odot }$ which have been reionized have an SFR 1.5 times greater than those which have not yet been reionized.

more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Medium: X Size: p. 3703-3716
["p. 3703-3716"]
Sponsoring Org:
National Science Foundation
More Like this

    We examine the dual [both black hole (BH) active] and offset (one BH active and in distinct galaxies) active galactic nucleus (AGN) population (comprising ∼ 2000 pairs at $0.5\, \text{kpc}\lesssim \Delta r\lt 30\, \text{kpc}$) at z = 2 ∼ 3 in the ASTRID simulation covering (360 cMpc)3. The dual (offset) AGN make up 3.0(0.5) per cent of all AGN at z = 2. The dual fraction is roughly constant while the offset fraction increases by a factor of 10 from z = 4 ∼ 2. Compared with the full AGN population, duals are characterized by low MBH/M* ratios, high specific star formation rates (sSFR) of $\sim 1\, \text{Gyr}^{-1}$, and high Eddington ratios (∼0.05, double that of single AGN). Dual AGNs are formed in major galaxy mergers (typically involving $M_\text{halo}\lt 10^{13}\, M_\odot$), with simular-mass BHs. At small separations (when host galaxies are in the late phase of the merger), duals become 2 ∼ 8 times brighter (albeit more obscured) than at larger separations. 80  per cent of the bright, close duals would merge within $\sim 500\, \text{Myr}$. Notably, the initially less-massive BHs in duals frequently become the brighter AGN during galaxy mergers. In offset AGN, the active BH is typically ≳ 10 times more massive than its non-active counterpart and than most BHs in duals. Offsets are predominantly formed in minor galaxy mergers with the active BH residing in the centre of massive haloes ($M_\text{ halo}\sim 10^{13-14}\, \mathrm{M}_\odot$). In these deep potentials, gas stripping is common and the secondary quickly deactivates. The stripping also leads to inefficient orbital decay amongst offsets, which stall at $\Delta r\sim 5\, \text{kpc}$ for a few hundred Myrs.

    more » « less

    Our understanding of reionization has advanced considerably over the past decade, with several results now demonstrating that the intergalactic medium transitioned from substantially neutral at z = 7 to largely reionized at z = 6. However, little remains known about the sizes of ionized bubbles at z ≳ 7 as well as the galaxy overdensities which drive their growth. Fortunately, rest-ultraviolet (UV) spectroscopic observations offer a pathway towards characterizing these ionized bubbles thanks to the resonant nature of Lyman-alpha photons. In a previous work, we presented Ly α detections from three closely separated Lyman-break galaxies at z ≃ 6.8, suggesting the presence of a large (R > 1 physical Mpc) ionized bubble in the 1.5 deg2 COSMOS field. Here, we present new deep Ly α spectra of 10 UV-bright ($\mathrm{\mathit{ M}}_{\mathrm{UV}}^{} \le -20.4$) z ≃ 6.6–6.9 galaxies in the surrounding area, enabling us to better characterize this potential ionized bubble. We confidently detect (S/N > 7) Ly α emission at z = 6.701–6.882 in nine of ten observed galaxies, revealing that the large-scale volume spanned by these sources (characteristic radius R = 3.2 physical Mpc) traces a strong galaxy overdensity (N/〈N〉 ≳ 3). Our data additionally confirm that the Ly α emission of UV-bright galaxies in this volume is significantly enhanced, with 40 per cent (4/10) showing strong Ly α emission (equivalent width >25 Å) compared to the 8–9 per cent found on average at z ∼ 7. The median Ly α equivalent width of our observed galaxies is also ≈2 times that typical at z ∼ 7, consistent with expectations if a very large (R ∼ 3 physical Mpc) ionized bubble is allowing the Ly α photons to cosmologically redshift far into the damping wing before encountering H i.

    more » « less
  3. Abstract A fundamental requirement for reionizing the Universe is that a sufficient fraction of the ionizing photons emitted by galaxies successfully escapes into the intergalactic medium. However, due to the scarcity of high-redshift observational data, the sources driving reionization remain uncertain. In this work, we calculate the ionizing escape fractions (fesc) of reionization-era galaxies from the state-of-the-art thesan simulations, which combine an accurate radiation-hydrodynamic solver (arepo-rt) with the well-tested IllustrisTNG galaxy formation model to self-consistently simulate both small-scale galaxy physics and large-scale reionization throughout a large patch of the universe ($L_\text{box} = 95.5\, \text{cMpc}$). This allows the formation of numerous massive haloes ($M_\text{halo} \gtrsim 10^{10}\, {\text{M}_{\odot }}$), which are often statistically underrepresented in previous studies but are believed to be important to achieving rapid reionization. We find that low-mass galaxies ($M_\text{stars} \lesssim 10^7\, {\text{M}_{\odot }}$) are the main drivers of reionization above z ≳ 7, while high-mass galaxies ($M_\text{stars} \gtrsim 10^8\, {\text{M}_{\odot }}$) dominate the escaped ionizing photon budget at lower redshifts. We find a strong dependence of fesc on the effective star formation rate (SFR) surface density defined as the SFR per gas mass per escape area, i.e. $\bar{\Sigma }_\text{SFR} = \text{SFR}/M_\text{gas}/R_{200}^2$. The variation in halo escape fractions decreases for higher mass haloes, which can be understood from the more settled galactic structure, SFR stability, and fraction of sightlines within each halo significantly contributing to the escaped flux. Dust is capable of reducing the escape fractions of massive galaxies, but the impact on the global fesc depends on the dust model. Finally, active galactic nuclei are unimportant for reionization in thesan and their escape fractions are lower than stellar ones due to being located near the centres of galaxy gravitational potential wells. 
    more » « less

    We present band 6 ALMA observations of a heavily obscured radio-loud (L1.4 GHz = 1025.4 W Hz−1) active galactic nucleus (AGN) candidate at zphot = 6.83 ± 0.06 found in the 1.5 deg2 COSMOS field. The ALMA data reveal detections of exceptionally strong [C ii]158 $\mu$m (z[C ii] = 6.8532) and underlying dust continuum emission from this object (COS-87259), where the [C ii] line luminosity, line width, and 158 $\mu$m continuum luminosity are comparable to those seen from z ∼ 7 sub-mm galaxies and quasar hosts. The 158 $\mu$m continuum detection suggests a total infrared luminosity of $9\times 10^{12}\, \mathrm{ L}_\odot$ with corresponding very large obscured star formation rate (1300 M⊙ yr−1) and dust mass ($2\times 10^9\, \mathrm{ M}_\odot$). The strong break seen between the VIRCam and IRAC photometry perhaps suggests that COS-87259 is an extremely massive reionization-era galaxy with $M_\ast \approx 1.7\times 10^{11}\, \mathrm{ M}_\odot$. Moreover, the MIPS, PACS, and SPIRE detections imply that this object harbours an AGN that is heavily obscured ($\tau _{_{\mathrm{9.7\,\mu m}}}=2.3$) with a bolometric luminosity of approximately $5\times 10^{13}\, \mathrm{ L}_\odot$. Such a very high AGN luminosity suggests that this object is powered by an ≈1.6 × 10$^9\, \mathrm{ M}_\odot$ black hole if accreting near the Eddington limit, and is effectively a highly obscured version of an extremely ultraviolet (UV)-luminous (M1450 ≈ −27.3) z ∼ 7 quasar. Notably, these z ∼ 7 quasars are an exceedingly rare population (∼0.001 deg−2), while COS-87259 was identified over a relatively small field. Future very wide area surveys with e.g. Roman and Euclid have the potential to identify many more extremely red yet UV-bright z ≳ 7 objects similar to COS-87259, providing richer insight into the occurrence of intense obscured star formation and supermassive black hole growth among this population.

    more » « less

    We introduce the thesan project, a suite of large volume ($L_\mathrm{box} = 95.5 \, \mathrm{cMpc}$) radiation-magnetohydrodynamic simulations that simultaneously model the large-scale statistical properties of the intergalactic medium during reionization and the resolved characteristics of the galaxies responsible for it. The flagship simulation has dark matter and baryonic mass resolutions of $3.1 \times 10^6\, {\rm M_\odot }$ and $5.8 \times 10^5\, {\rm M_\odot }$, respectively. The gravitational forces are softened on scales of 2.2 ckpc with the smallest cell sizes reaching 10 pc at z = 5.5, enabling predictions down to the atomic cooling limit. The simulations use an efficient radiation hydrodynamics solver (arepo-rt) that precisely captures the interaction between ionizing photons and gas, coupled to well-tested galaxy formation (IllustrisTNG) and dust models to accurately predict the properties of galaxies. Through a complementary set of medium resolution simulations we investigate the changes to reionization introduced by different assumptions for ionizing escape fractions, varying dark matter models, and numerical convergence. The fiducial simulation and model variations are calibrated to produce realistic reionization histories that match the observed evolution of the global neutral hydrogen fraction and electron scattering optical depth to reionization. They also match a wealth of high-redshift observationally inferred data, including the stellar-to-halo-mass relation, galaxy stellar mass function, star formation rate density, and the mass–metallicity relation, despite the galaxy formation model being mainly calibrated at z = 0. We demonstrate that different reionization models give rise to varied bubble size distributions that imprint unique signatures on the 21 cm emission, especially on the slope of the power spectrum at large spatial scales, enabling current and upcoming 21 cm experiments to accurately characterize the sources that dominate the ionizing photon budget.

    more » « less