51ÁÔÆæÈë¿Ú

Abstract:

We study the stellar mass-iron metallicity relation of dwarf galaxies in the new high-resolution megatron cosmological radiation-hydrodynamics simulations. These simulations model galaxy formation up to z ≈ 8 in a region that will collapse into a Milky-Way-like galaxy at z = 0, while self-consistently tracking Population III and II (Pop. III, Pop. II) star formation, feedback and chemical enrichment. megatron dwarf galaxies are in excellent agreement with the observed stellar mass-metallicity relation at z = 0, including an over-abundance of dwarfs along a flat plateau in metallicity (⟨[Fe/H]⟩ ≈ −2.5) at low stellar masses (M⋆ ≤ 10⁵ M⊙). We tie this feature to the chemical enrichment of dwarf galaxies by Pop. III pair-instability supernova (PISN) explosions. The strong Lyman-Werner background (LW) from the protogalaxy ensures that PISNe occur in haloes massive enough (≈ 10⁷ M⊙) to retain their ejecta. We also predict a tail of ≈ 20% of iron-deficient (⟨[Fe/H]⟩ ≤ −3) dwarf galaxies. We show that both plateau and tail (i) are robust to large variations in Pop. II feedback assumptions, and (ii) survive in bound satellites surrounding the central galaxy at z = 0.

Abstract:

We present the MEGATRON suite of cosmological radiation hydrodynamics simulations following the formation of Milky Way-mass galaxies from the earliest cosmic epochs when Population III stars form to Cosmic Noon. The suite represents the first set of cosmological simulations that couples a vast nonequilibrium thermochemistry network of primordial species, metals, and molecules to multifrequency, on-the-fly radiation transport, allowing us to directly predict the spectral properties of early galaxies. By initializing the simulations at zero metallicity, resolving haloes well below the atomic cooling threshold, reaching parsec-scale resolution, and modeling a Milky Way-mass environment, we aim to address four key science themes: 1) Star formation at cosmic dawn, 2) Galaxy formation and the interstellar medium in the epoch of reionization, 3) The circumgalactic medium towards cosmic noon, and 4) Reionization in a local volume environment and near-field cosmology. In this introductory work, we present an overview of the physical characteristics of high-redshift MEGATRON galaxies and their environment at z > 8. We present a library of > 175, 000 simulated galaxy spectra and demonstrate how much of the diversity of galaxy spectra seen by JWST is naturally reproduced in the context of a ΛCDM cosmology. Caveats are discussed, such as the lack of AGN in our simulations and the limitations of our adopted stellar population and chemical yield models. This project represents a step towards making more direct comparisons between simulations and observations and is particularly applicable for optimizing methods to infer galaxy properties from existing high-redshift JWST spectra and imaging data.

Abstract:

Abstract Despite the vast amount of energy released by active galactic nuclei (AGN), their role in early galaxy formation and in regulating the growth of supermassive black holes (SMBHs) remains poorly understood. Through new high-resolution zoom-in cosmological simulations, we follow the co-evolution of 105 M⊙ black hole seeds with their host dwarf galaxy. We model ionizing feedback from a Pop III.1 progenitor, applicable to a wide range of internally or externally irradiated SMBH formation scenarios. The simulated suite progressively spans physics ranging from no AGN feedback to more complex setups including thermal, kinetic and radiative feedback – explored for both low and enhanced AGN power. Across all our models, we find that black hole seeds efficiently reach masses of ∼107 M⊙ within a ∼1010 M⊙ halo by z = 8. Although they exhibit notably different mass growth histories, these latter seem unimpeded by the presence of AGN feedback. The simulation including radiative feedback is the most distinct, with super-Eddington episodes driving fast and mass-loaded gas outflows (exceeding 2500 km s−1) up to ∼50 kpc, along with minor stellar mass suppression in the host galaxy. Our measurements are in broad agreement with moderate luminosity quasars recently observed by JWST, producing overmassive black holes (SMBH-to-galaxy mass ratios 0.01 − 1), dynamical masses of ∼109.5 M⊙, stellar masses of ∼108.5 M⊙, and high, though short-lived, Eddington fraction accretion rates. These results advocate for a scenario where AGN feedback allows for rapid SMBH growth during the reionisation era, while driving winds that extend deep into the intergalactic medium – shaping host galaxies as well as more distant surroundings.

Abstract:

Abstract The peculiar velocity field of the local Universe provides direct insights into its matter distribution and the underlying theory of gravity, and is essential in cosmological analyses for modelling deviations from the Hubble flow. Numerous methods have been developed to reconstruct the density and velocity fields at z ≲ 0.05, typically constrained by redshift-space galaxy positions or by direct distance tracers such as the Tully–Fisher relation, the fundamental plane, or Type Ia supernovae. We introduce a validation framework to evaluate the accuracy of these reconstructions against catalogues of direct distance tracers. Our framework assesses the goodness-of-fit of each reconstruction using Bayesian evidence, residual redshift discrepancies, velocity scaling, and the need for external bulk flows. Applying this framework to a suite of reconstructions—including those derived from the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm and from linear theory—we find that the non-linear BORG reconstruction consistently outperforms others. We highlight the utility of such a comparative approach for supernova or gravitational wave cosmological studies, where selecting an optimal peculiar velocity model is essential. Additionally, we present calibrated bulk flow curves predicted by the reconstructions and perform a density–velocity cross-correlation using a linear theory reconstruction to constrain the growth factor, yielding S8 = 0.793 ± 0.035. The result is in good agreement with both weak lensing and Planck, but is in strong disagreement with some peculiar velocity studies.