91̽»¨

Abstract:

JWST observations are providing unprecedented constraints on the history of reionization owing to the ability to detect faint galaxies at z ≫ 6. Modelling this history requires understanding both the ionizing photon production rate (ξion) and the fraction of those photons that escape into the intergalactic medium (fesc). Observational estimates of these quantities generally rely on spectroscopy for which large samples with well-defined selection functions remain limited. To overcome this challenge, we present and release a novel implicit likelihood inference pipeline, PHOTONIOn, trained on mock photometry to predict the escaped ionizing luminosity of individual galaxies (N ion) based on photometric magnitudes and redshifts. We show that PHOTONIOn is able to reliably infer N ion from photometry. This is in contrast to traditional spectral energy distribution-fitting approaches which rely on fesc prescriptions that often overpredict N ion for Lyman Continuum (LyC)-dim galaxies, even when given access to spectroscopic data. We have deployed PHOTONIOn on a sample of 4559 high-redshift galaxies from the JWST Advanced Deep Extragalactic Survey (JADES), finding gentle redshift evolutions of log10(N ion) = (0.08 ± 0.01)z + (51.60 ± 0.06) and log10(fescξion) = (0.07 ± 0.01)z + (24.12 ± 0.07). Late-time values for the ionizing photon production rate density are consistent with both theoretical models and observations. Finally, we measure the evolution of the intergalactic medium ionized fraction to find that observed populations of star-forming galaxies are capable of driving reionization in this field to completion by z ∼ 5.3 without the need for active galactic nucleus or other exotic sources, consistent with other studies of the same field. The 20 per cent of UV-brightest galaxies (MUV < −18.5) reionize roughly 35 per cent of the survey volume, demonstrating that UV faint LyC emitters are crucial for reionization.

Abstract:

<jats:p>Star formation is a key process that governs the baryon cycle within galaxies, however, the question of how it controls their growth remains elusive due to modeling uncertainties. To understand the impact of star formation models on galaxy evolution, we performed cosmological zoom-in radiation-hydrodynamic simulations of a dwarf dark matter halo, with a virial mass of <jats:italic>M</jats:italic><jats:sub>vir</jats:sub> ∼ 10<jats:sup>9</jats:sup> <jats:italic>M</jats:italic><jats:sub>⊙</jats:sub> at <jats:italic>z</jats:italic> = 6. We compared two different star formation models: a multi-freefall model combined with a local gravo-thermo-turbulent condition and a more self-consistent model based on a sink particle algorithm, where gas accretion and star formation are directly controlled by the gas kinematics. As the first study in this series, we used cosmological zoom-in simulations with different spatial resolutions and found that star formation is more bursty in the runs with the sink algorithm, generating stronger outflows than in the runs with the gravo-thermo-turbulent model. The main reason for the increased burstiness is that the gas accretion rates on the sinks are high enough to form stars on very short timescales, leading to more clustered star formation. As a result, the star-forming clumps are disrupted more quickly in the sink run due to more coherent radiation and supernova feedback. The difference in burstiness between the two star formation models becomes even more pronounced when the supernova explosion energy is artificially increased. Our results suggest that improving the modeling of star formation on small, sub-molecular cloud scales can significantly impact the global properties of simulated galaxies.</jats:p>

Abstract:

Constrained cosmological simulations play an important role in modelling the local Universe, enabling investigation of the dark matter content of local structures and their formation. We introduce an internal method for quantifying the extent to which the variance of individual halo properties is suppressed by the constraints imposed on the initial conditions. We apply it to the Constrained Simulations in BORG (CSiBORG) suite of 101 high-resolution realizations across the posterior probability distribution of initial conditions from the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm. The method is based on the overlap of the initial Lagrangian patch of a halo in one simulation with those in another, measuring the degree to which the haloes' particles are initially coincident. This addresses the extent to which the imposed large-scale structure constraints reduce the variance of individual halo properties. We find consistent reconstructions of M≳1014M⊙h-1 haloes, indicating that the constraints from the BORG algorithm are sufficient to pin down the masses, positions, and peculiar velocities of clusters to high precision, though we do not assess how well they reproduce observations of the local Universe. The effect of the constraints tapers off towards lower mass, and the halo spins and concentrations are largely unconstrained at all masses. We document the advantages of evaluating halo consistency in the initial conditions and describe how the method may be used to quantify our knowledge of the halo field given galaxy survey data analysed through the lens of probabilistic inference machines such as BORG.

Abstract:

Aims: The true nature of primordial magnetic fields (PMFs) and their role in the formation of galaxies still remains elusive. To shed light on these unknowns, we investigate their impact by varying two sets of properties: (i) accounting for the effect of PMFs on the initial matter power spectrum, and (ii) accounting for their magneto-hydrodynamical effects on the formation of galaxies. By comparing both we can determine the dominant agent in shaping galaxy evolution.

Methods: We use the magneto-hydrodynamics code RAMSES, to generate multiple new zoom-in simulations for eight different host halos of dwarf galaxies across a wide luminosity range of 103 − 106 L⊙. These halos are selected from a ΛCDM cosmological box, tracking their evolution down to redshift z = 0. We explore a variety of primordial magnetic field (comoving) strengths Bλ ranging from 0.05 to 0.50 nG.

Results: We find magnetic fields in the interstellar medium not only modify star formation in dwarf spheroidal galaxies but also completely prevent the formation of stars in less compact ultra-faints with halo mass and stellar mass below ∼ 2.5 · 109 and 3 · 106 M⊙, respectively. At high redshifts, the impact of PMFs on host halos of dwarf galaxies through the modification of the matter power spectrum is more dominant than the influence of magneto-hydrodynamics in shaping their gaseous structure. Through the amplification of small perturbations ranging in mass from 107 to 109 M⊙ in the ΛCDM+PMFs matter power spectrum, primordial fields expedite the formation of the first dark matter halos, leading to an earlier onset and a higher star formation rate at redshifts z > 12. We investigate the evolution of various energy components and demonstrate that magnetic fields with an initial strength of Bλ ≥ 0.05 nG exhibit a strong growth of magnetic energy, accompanied by a saturation phase, that starts quickly after the growth phase. These trends persist consistently, regardless of the initial conditions, whether it is the classical ΛCDM or modified by PMFs. Lastly, we investigate the impact of PMFs on the present-time observable properties of dwarf galaxies, namely, the half light radius, V-band luminosity, mean metallicity and velocity dispersion profile. We find that PMFs with moderate strengths of Bλ ≤ 0.10 nG show great agreement with the scaling relations of the observed Local group dwarfs. However, stronger fields lead to large sizes and high velocity dispersion.