Search

Dust in galaxies is an important tracer of galaxy properties and their evolution over time. The physical origin of the grain size distribution, the dust chemical composition, and, hence, the associated ultraviolet-to-optical extinctions in diverse galaxies remains elusive. To address this issue, we introduce a model for dust evolution in the RAMSES code for simulations of galaxies with a resolved multiphase interstellar medium. Dust is modelled as a fluid transported with the gas component, and is decomposed into two sizes, 5 nm and 0.1 μm, and two chemical compositions for carbonaceous and silicate grains. This dust model includes the growth of dust by accretion of elements from the gas phase and by the release of dust in stellar ejecta, the destruction by thermal sputtering, supernovae, and astration, and the exchange of dust mass between the two main populations of grain sizes by coagulation and shattering. Using a suite of isolated disc simulations with different masses and metallicities, the simulations can explore the role of these processes in shaping the key properties of dust in galaxies. The simulated Milky Way analogue reproduces the dust-to-metal mass ratio, depletion factors, size distribution and extinction curves of the Milky Way. Galaxies with lower metallicities reproduce the observed decrease in the dust-to-metal mass ratio with metallicity at around a few 0.1 Z⊙. This break in the dust-to-metal ratio corresponds to a galactic gas metallicity threshold that marks the transition from an ejecta-dominated to an accretion-dominated grain growth, and that is different for silicate and carbonaceous grains, with ≃0.1 Z⊙ and ≃0.5 Z⊙ respectively. This leads to more Magellanic Cloud-like extinction curves, i.e. with steeper slopes in the ultraviolet and a weaker bump feature at 2175 Å, in galaxies with lower masses and lower metallicities. Steeper slopes in these galaxies are caused by the combination of the higher efficiency of gas accretion by silicate relative to carbonaceous grains and by the low rates of coagulation that preserves the amount of small silicate grains. Weak bumps are due to the overall inefficient accretion growth of carbonaceous dust at low metallicity, whose growth is mostly supported by the release of large grains in SN ejecta. We also show that the formation of CO molecules is a key component to limit the ability of carbonaceous dust to grow, in particular in low-metallicity gas-rich galaxies. [ABSTRACT FROM AUTHOR]

We explore how the characteristics of the cross-correlation functions between the 21-cm emission from the spin-flip transition of neutral hydrogen (H  |$\scriptstyle \rm I $|⁠) and early Lyman α (Ly α) radiation emitting galaxies (Ly α emitters, LAEs) depend on the reionization history and topology and the simulated volume. For this purpose, we develop an analytic expression for the 21-cm-LAE cross-correlation function and compare it to results derived from different astraeus and 21cmfast reionization simulations covering a physically plausible range of scenarios where either low-mass (≲ 109.5 M⊙) or massive (≳ 109.5 M⊙) galaxies drive reionization. Our key findings are: (i) the negative small-scale (≲  2 cMpc) cross-correlation amplitude scales with the intergalactic medium's (IGM) average H  |$\scriptstyle \rm I $| fraction (〈χ H i 〉) and spin-temperature weighted overdensity in neutral regions (〈1 + δ〉 H i ); (ii) the inversion point of the cross-correlation function traces the peak of the size distribution of ionized regions around LAEs; (iii) the cross-correlation amplitude at small scales is sensitive to the reionization topology, with its anticorrelation or correlation decreasing the stronger the ionizing emissivity of the underlying galaxy population is correlated to the cosmic web gas distribution (i.e. the more low-mass galaxies drive reionization); (iv) the required simulation volume to not underpredict the 21-cm-LAE anticorrelation amplitude when the cross-correlation is derived via the cross-power spectrum rises as the size of ionized regions and their variance increases. Our analytic expression can serve two purposes: to test whether simulation volumes are sufficiently large, and to act as a fitting function when cross-correlating future 21-cm signal Square Kilometre Array and LAE galaxy observations. [ABSTRACT FROM AUTHOR]

We use the astraeus framework to investigate how the visibility and spatial distribution of Lyman-α (Lyα) emitters (LAEs) during reionization is sensitive to a halo mass-dependent fraction of ionizing radiation escaping from the galactic environment (f esc) and the ionization topology. To this end, we consider the two physically plausible bracketing scenarios of f esc increasing and decreasing with rising halo mass. We derive the corresponding observed Lyα luminosities of galaxies for three different analytic Lyα line profiles and associated Lyα escape fraction (⁠|$f_\mathrm{esc}^\mathrm{Ly\alpha }$|⁠) models: importantly, we introduce two novel analytic Lyα line profile models that describe the surrounding interstellar medium (ISM) as dusty gas clumps. They are based on parameterizing results from radiative transfer simulations, with one of them relating |$f_\mathrm{esc}^\mathrm{Ly\alpha }$| to f esc by assuming the ISM of being interspersed with low-density tunnels. Our key findings are: (i) for dusty gas clumps, the Lyα line profile develops from a central to double peak dominated profile as a galaxy's halo mass increases; (ii) LAEs are galaxies with |$M_h\gtrsim 10^{10}\, {\rm \rm M_\odot }$| located in overdense and highly ionized regions; (iii) for this reason, the spatial distribution of LAEs is primarily sensitive to the global ionization fraction and only weakly in second-order to the ionization topology or a halo mass-dependent f esc; (iv) furthermore, as the observed Lyα luminosity functions reflect the Lyα emission from more massive galaxies, there is a degeneracy between the f esc-dependent intrinsic Lyα luminosity and the Lyα attenuation by dust in the ISM if f esc does not exceed |$\sim 50~{{\ \rm per \, cent}}$|⁠. [ABSTRACT FROM AUTHOR]

The Ly  α line is a powerful probe of distant galaxies, which contains information about inflowing/outflowing gas through which Ly  α photons scatter. To develop our understanding of this probe, we post-process a zoom-in radiation-hydrodynamics simulation of a low-mass (M star ∼ 109 M⊙) galaxy to construct 22 500 mock spectra in 300 directions from z  = 3 to 4. Remarkably, we show that one galaxy can reproduce the variety of a large sample of spectroscopically observed Ly  α line profiles. While most mock spectra exhibit double-peak profiles with a dominant red peak, their shapes cover a large parameter space in terms of peak velocities, peak separation, and flux ratio. This diversity originates from radiative transfer effects at interstellar medium and circum-galactic medium (CGM) scales, and depends on galaxy inclination and evolutionary phase. Red-dominated lines preferentially arise in face-on directions during post-starburst outflows and are bright. Conversely, accretion phases usually yield symmetric double peaks in the edge-on direction and are fainter. While resonant scattering effects at <0.2 × R vir are responsible for the broadening and velocity shift of the red peak, the extended CGM acts as a screen and impacts the observed peak separation. The ability of simulations to reproduce observed Ly  α profiles and link their properties with galaxy physical parameters offers new perspectives to use Ly  α to constrain the mechanisms that regulate galaxy formation and evolution. Notably, our study implies that deeper Ly  α surveys may unveil a new population of blue-dominated lines tracing inflowing gas. [ABSTRACT FROM AUTHOR]

In this work, the sixth of a series, we use the seminumerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N -body dark-matter simUlationS (astraeus) framework to investigate the nature of the sources that reionized the Universe. We extend astraeus , which already couples a galaxy formation semi-analytical model with a detailed seminumerical reionization scheme, to include a model for black-hole formation, growth, and the production of ionizing radiation from associated active galactic nuclei (AGNs). We calibrate our fiducial AGN model to reproduce the bolometric luminosity function at z ≃ 5, and explore the role of the resulting AGN population in reionizing the Universe. We find that in all the models yielding a reasonable AGN luminosity function, galaxies dominate overwhelmingly the ionizing budget during the Epoch of Reionization, with AGN accounting for 1–10 per cent of the ionizing budget at z  = 6 and starting to play a role only below z ≲ 5. [ABSTRACT FROM AUTHOR]