Search

The James Webb Space Telescope (JWST) has uncovered an abundance of $z>10$ galaxies bright in the ultraviolet (UV) that has challenged traditional theoretical models at high redshifts. Recently, various new models have emerged to address this discrepancy by refining their description of star formation. Here we investigate whether modifications to the stellar initial mass function (IMF) alone can reproduce the $z>10$ UV luminosity functions (UV LFs) when the star formation rate is used as a proxy for the fraction of massive stars. We incorporate an Evolving IMF into the {\sc astraeus} galaxy evolution and reionisation simulation framework, which becomes increasingly top-heavy as the gas density in a galaxy rises above a given threshold. Our implementation accounts for the IMF's effects on supernova (SN) feedback, metal enrichment, and the UV and ionising emissivities. For this Evolving IMF model, we find that (i) the maximum UV luminosity enhancement is twice as large in massive galaxies ($\Delta M_\mathrm{UV}\simeq2.6$) than those where star formation is strongly limited by SN feedback ($\Delta M_\mathrm{UV}\simeq1.3$); (ii) it successfully reproduces the observed UV LFs at $z=5-15$; (iii) galaxies with top-heavy IMFs exhibit the highest star formation rates, driven by their location in local density peaks, which facilitates higher gas accretion rates; (iv) the $1\sigma$ variances in the UV luminosity are only slightly higher compared to when assuming a Salpeter IMF, but the $2\sigma$ variances are significantly increased by a factor of $1.4-2$ boosting the abundance of UV-bright galaxies at $z>10$; (v) reionisation begins earlier with more extended large ionised regions and fewer smaller ones during its initial stages, though these differences diminish at lower redshifts, leading to a similar end of reionisation at $z\simeq5.6$.}, Comment: submitted to A&A, 12 pages, 8 figures, comments welcome

We investigate the spatial offsets between dust and ultraviolet (UV) emission in high-redshift galaxies using the Cosmic Dawn III (CoDa III) simulation, a state-of-the-art fully coupled radiation-hydrodynamics cosmological simulation. Recent observations have revealed puzzling spatial disparities between ALMA dust continuum and UV emission as seen by HST and JWST in galaxies at z=5-7, compelling us to propose a physical interpretation of such offsets. Our simulation, which incorporates a dynamical dust model, naturally reproduces these offsets in massive, UV-bright galaxies (log$_{10}$(M$_{\rm{DM}}$/M$_{\odot}$)>11.5, M$_{\rm{AB1500}}$<-20). We find that dust-UV offsets increase with halo mass and UV brightness, reaching up to $\sim 2$ pkpc for the most massive systems, in good agreement with observational data from the ALPINE and REBELS surveys. Our analysis reveals that these offsets primarily result from severe dust extinction in galactic centers rather than a misalignment between dust and stellar mass distributions. The dust remains well-aligned with the bulk stellar component, and we predict the dust continuum should therefore align well with the stellar rest-frame NIR component, less affected by dust attenuation. This study provides crucial insights into the complex interplay between star formation, dust distribution, and observed galaxy morphologies during the epoch of reionization, highlighting the importance of dust in shaping the appearance of early galaxies at UV wavelengths., Comment: 9 pages, 5 figures, submitted to A&A

We present a new method - called HINORA (HIgh-NOise RANdom SAmple Consensus) - for the identification of regular structures in 3D point distributions. Motivated by the possible existence of the so called Council of Giants, i.e. a ring of twelve massive galaxies surrounding the Local Group in the Local Sheet with a radius of 3.75 Mpc, we apply HINORA to the Local Volume Galaxy catalogue confirming its existence. When varying the lower limit of K-band luminosity of the galaxy entering the catalogue, we further report on the existence of another ring-like structure in the Local Volume that now contains the Milky Way and M31. However, this newly found structure is dominated by low-mass (satellite) galaxies. While we here simply present the novel method as well as its first application to observational data, follow-up work using numerical simulations of cosmic structure formation shall shed light into the origin of such regular patterns in the galaxy distribution. Further, the method is equally suited to identify similar (or even different) structures in various kinds of astrophysical data (e.g. locating the actual 'baryonic-acoustic oscillation spheres' in galaxy redshift surveys)., Comment: Accepted for publication in PASA

This is the second paper in a series presenting the results from a 500 $h^{-1}$Mpc large constrained hydro-dynamical simulation of the local Universe (SLOW). The initial conditions are based on peculiar velocities derived from the CosmicFlows-2 catalogue. The inclusion of galaxy formation treatment, allows to directly predict observable properties of the Intra-Cluster Medium (ICM) within galaxy clusters. Comparing the properties of observed galaxy clusters within the local Universe with the properties of their simulated counterparts, enables us to assess the effectiveness of the initial condition constraints in accurately replicating the non-linear properties of the largest, collapsed objects within the simulation. Based on the combination of several, publicly available surveys, we identified 45 local Universe galaxy clusters in SLOW, including the 13 most massive from the Planck SZ catalog and 70% of those with $M_{500} > 2\times 10^{14}$ M$_{\odot}$. We then derived the probability of the cross identification based on mass, X-ray luminosity, temperature and Compton-y by comparing it to a random selection. In relation to previous constrained simulations of the local volume, we found in SLOW a much larger amount of replicated galaxy clusters, where their simulation based mass prediction falls within the uncertainties of the observational mass estimates. Comparing the median observed and simulated masses of our cross identified sample allows to independently deduce a hydrostatic mass bias of $(1-b)\approx0.87$. The SLOW constrained simulation of the local Universe faithfully reproduces numerous fundamental characteristics of the galaxy clusters within our local neighbourhood, opening a new avenue for studying the formation and evolution of a large set of individual galaxy clusters as well as testing our understanding of physical processes governing the ICM., Comment: Submitted to Astronomy & Astrophysics

Observations with JWST have revealed an unexpected high abundance of bright z>10 galaxy candidates. We explore whether a stellar initial mass function (IMF) that becomes increasingly top-heavy towards higher redshifts and lower gas-phase metallicities results in a higher abundance of bright objects in the early universe and how it affects the evolution of galaxy properties compared to a constant IMF. We incorporate such an evolving IMF into the Astraeus framework that couples galaxy evolution and reionisation in the first billion years. Our implementation accounts for the IMF dependence of supernova feedback, metal enrichment, ionising and ultraviolet radiation emission. We conduct two simulations: one with a Salpeter IMF and one with the evolving IMF. Compared to a constant Salpeter IMF, we find that (i) the higher abundance of massive stars in the evolving IMF results in more light per unit stellar mass, a slower build-up of stellar mass and lower stellar-to-halo mass ratio; (ii) due to the self-similar growth of the underlying dark matter halos, the evolving IMF's star formation main sequence hardly deviates from that of the Salpeter IMF; (iii) the evolving IMF's stellar mass-metallicity relation shifts to higher metallicities while its halo mass-metallicity relation remains unchanged; (iv) the evolving IMF's median dust-to-metal mass ratio is lower due to its stronger SN feedback; (v) the evolving IMF requires lower values of the escape fraction of ionising photons and exhibits a flatter median relation and smaller scatter between the ionising photons emerging from galaxies and the halo mass. Yet, the topology of the ionised regions hardly changes compared to the Salpeter IMF. These results suggest that a top-heavier IMF alone is unlikely to explain the higher abundance of bright z>10 sources, since the lower mass-to-light ratio is counteracted by the stronger stellar feedback., Comment: 19 pages, 13 figures, Accepted for publication in Astronomy & Astrophysics

The reconstruction of the large scale velocity field from the grouped Cosmicflows-4 (CF4) database is presented. The lognormal bias of the inferred distances and velocities data is corrected by the Bias Gaussianization correction (BGc) scheme, and the linear density and velocity fields are reconstructed by means of the Wiener filter (WF) and constrained realizations (CRs) algorithm. These tools are tested against a suite of random and constrained Cosmicflows-3-like mock data. The CF4 data consists of 3 main subsamples - the 6dFGS and the SDSS data - and the `others'. The individual contributions of the subsamples have been studied. The quantitative analysis of the velocity field is done mostly by the mean overdensity ($\Delta_L(R)$) and the bulk velocity ($V_{\mathrm{bulk}}(R)$) profiles of the velocity field out to $300\, h^{-1}{\rm Mpc}$. The $V_{\mathrm{bulk}}(R)$ and $\Delta_{\mathrm L}(R)$ profiles of the CF4 data without its 6dFGS component are consistent with the cosmic variance to within $1\sigma$. The 6dFGS sample dominates the $V_{\mathrm{bulk}}$ ($\Delta_{\mathrm L}$) profile beyond $\sim120\, h^{-1}{\rm Mpc}$, and drives it to roughly a $3.4\sigma$ ($-1.9\sigma$) excess (deficiency) relative to the cosmic variance at $R\sim250\ (190)\ \, h^{-1}{\rm Mpc}$. The excess in the amplitude of $V_{\mathrm{bulk}}$ is dominated by its Supergalactic X component, roughly in the direction of the Shapley Concentration. The amplitude and alignment of the inferred velocity field from the CF4 data is at $\sim(2\,-\,3)\,\sigma$ discrepancy with respect to the $\Lambda$CDM model. Namely, it is somewhat atypical but yet there is no compelling tension with the model., Comment: 17 pages, 14 figures, accepted for publication to MNRAS

We study the properties of associations of dwarf galaxies and their dependence on the environment. Associations of dwarf galaxies are extended systems composed exclusively of dwarf galaxies, considering as dwarf galaxies those galaxies less massive than $M_{\star, \rm max} = 10^{9.0}$ ${\rm M}_{\odot}\,h^{-1}$. We identify these particular systems using a semi-analytical model of galaxy formation coupled to a dark matter only simulation in the $\Lambda$ Cold Dark Matter cosmological model. To classify the environment, we estimate eigenvalues from the tidal field of the dark matter particle distribution of the simulation. We find that the majority, two thirds, of associations are located in filaments ($ \sim 67$ per cent), followed by walls ($ \sim 26 $ per cent), while only a small fraction of them are in knots ($ \sim 6 $ per cent) and voids ($ \sim 1 $ per cent). Associations located in more dense environments present significantly higher velocity dispersion than those located in less dense environments, evidencing that the environment plays a fundamental role in their dynamical properties. However, this connection between velocity dispersion and the environment depends exclusively on whether the systems are gravitational bound or unbound, given that it disappears when we consider associations of dwarf galaxies that are gravitationally bound. Although less than a dozen observationally detected associations of dwarf galaxies are currently known, our results are predictions on the eve of forthcoming large surveys of galaxies, which will enable these very particular systems to be identified and studied., Comment: 13 pages, 9 figures. Accepted for publication in MNRAS

We explore how the characteristics of the cross-correlation functions between the 21cm emission from the spin-flip transition of neutral hydrogen (HI) and early Lyman-$\alpha$ (Ly$\alpha$) radiation emitting galaxies (Ly$\alpha$ emitters, LAEs) depend on the reionisation history and topology and the simulated volume. For this purpose, we develop an analytic expression for the 21cm-LAE cross-correlation function and compare it to results derived from different Astraeus and 21cmFAST reionisation simulations covering a physically plausible range of scenarios where either low-mass ($<10^{9.5}M_\odot$) or massive ($>10^{9.5}M_\odot$) galaxies drive reionisation. Our key findings are: (i) the negative small-scale ($<2$ cMpc) cross-correlation amplitude scales with the intergalactic medium's (IGM) average HI fraction ($\langle\chi_\mathrm{HI}\rangle$) and spin-temperature weighted overdensity in neutral regions ($\langle1+\delta\rangle_\mathrm{HI}$); (ii) the inversion point of the cross-correlation function traces the peak of the size distribution of ionised regions around LAEs; (iii) the cross-correlation amplitude at small scales is sensitive to the reionisation topology, with its anti-correlation or correlation decreasing the stronger the ionising emissivity of the underlying galaxy population is correlated to the cosmic web gas distribution (i.e. the more low-mass galaxies drive reionisation); (iv) the required simulation volume to not underpredict the 21cm-LAE anti-correlation amplitude when the cross-correlation is derived via the cross-power spectrum rises as the size of ionised 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 21cm signal Square Kilometre Array and LAE galaxy observations., Comment: 13 pages, 5 figures; accepted for publication in MNRAS

The aim of cosmological simulations is to reproduce the properties of the observed Universe, serving as tools to test structure and galaxy formation models. Constrained simulations of our local cosmological region up to a few hundred Mpc/h , the local Universe, are designed to reproduce the actual cosmic web of structures as observed. A question that often arises is how to judge the quality of constrained simulations against the observations of the Local Universe. Here we introduce the Local Universe model (LUM), a new methodology, whereby many constrained simulations can be judged and the ''best'' initial conditions can be identified. By characterising the Local Universe as a set of rich clusters, the model identifies haloes that serve as simulated counterparts to the observed clusters. Their merit is determined against a null hypothesis, the probability that such a counterpart could be identified in a random, unconstrained simulation. This model is applied to 100 constrained simulations using the Cosmicflows-3 data. Cluster counterparts are found for all constrained simulations, their distribution of separation from the true observed cluster position and their mass distribution are investigated. Lastly, the ''best'' constrained simulation is selected using the LUM and discussed in more detail., Comment: 10 pages, 8 figures, submitted to MNRAS

Context: Several observations of the local Universe (LU) point towards the existence of very prominent structures. The presence of massive galaxy clusters and local super clusters on the one hand, but also large local voids and under-densities on the other hand. However, it is highly non trivial to connect such different observational selected tracers to the underlying dark matter (DM) distribution. Methods (abridged): We used a 500 Mpc/h large constrained simulation of the LU with initial conditions based on peculiar velocities derived from the CosmicFlows-2 catalogue and follow galaxy formation physics directly in the hydro-dynamical simulations to base the comparison on stellar masses of galaxies or X-ray luminosity of clusters. We also used the 2668 Mpc/h large cosmological box from the Magneticum simulations to evaluate the frequency of finding such anomalies in random patches within simulations. Results: We demonstrate that haloes and galaxies in our constrained simulation trace the local DM density field very differently. Thereby, this simulation reproduces the observed 50% under-density of galaxy clusters and groups within the sphere of ~100 Mpc when applying the same mass or X-ray luminosity limit used in the observed cluster sample (CLASSIX), which is consistent with a ~1.5$\sigma$ feature. At the same time, the simulation reproduces the observed over-density of massive galaxy clusters within the same sphere, which on its own also corresponds to a ~1.5$\sigma$ feature. Interestingly, we find that only 44 out of 15635 random realizations (i.e. 0.28%) are matching both anomalies, making the LU to be a ~3$\sigma$ environment. We finally compared a mock galaxy catalogue with the observed distribution of galaxies in the LU, finding also a match to the observed factor of two over-density at ~16 Mpc as well as the observed 15% under-density at ~40 Mpc distance., Comment: 16 pages, 12 figure, accepted for publication in A&A

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

Recent studies based on numerical models of the Local Group predict the existence of field haloes and galaxies that have visited two distinct galaxies in the past, called Hermeian haloes. This work presents an analysis of the Hermeian haloes population in two high-resolution dark matter-only N-body simulations from the MultiDark suit (the ESMDPL, VSMDPL). Hermeian haloes make up from 0.5 to 2.5 per cent of the total number of field haloes depending on their mass. Furthermore, the results of our study suggest that at a sufficiently high resolution simulation, Hermeian haloes may be found around almost every halo, making them interesting for studies of matter exchange between galaxies. We find that about half (22 out of 49) of the selected Local Group analogues contain Hermeian haloes that passed through the haloes of both the Milky Way and M31 if the distance between the two main haloes is below 1 $h^{-1} \;\mathrm{Mpc}$; this fraction drops to one-fifth (24 out of 108) for distances of up to 1.5 $h^{-1} \;\mathrm{Mpc}$. We confirm earlier findings that unlike other field haloes, Hermeians are grouped along the line connecting the primary hosts of Local Group-like systems, which should facilitate their identification in observations. The vast majority of the Hermeian haloes, whose second target is a Milky Way analogue, are currently moving away from it with increased velocity compared to remaining field halo populations. Interestingly, the obtained data admits that NGC 3109 could have passed through the Andromeda Galaxy and the Milky Way earlier., Comment: 13 pages, 13 figures, submitted to MNRAS

How the cosmic web feeds halos, and fuels galaxy formation is an open question with wide implications. This study explores the mass assembly in the Local Group within the context of the local cosmography by employing simulations whose initial conditions have been constrained to reproduce the local environment. The goal of this study is to inspect whether the direction of accretion of satellites on to the Milky Way and Andromeda galaxies, is related to the cosmic web. The analysis considers the three high-resolution simulations available in the HESTIA simulation suite, as well as the derived velocity shear and tidal tensors. We notice two eras in the Local Group accretion history, delimited by an epoch around $z \approx 0.7$. We also find that satellites can travel up to $\sim 4$ Mpc, relative to their parent halo before crossing its viral radius $R_{200}$. Finally, we observe a strong alignment of the infall direction with the axis of slowest collapse $\vec{e_3}$ of both tidal and shear tensors, implying satellites of the Local Group originated from one particular region of the cosmic web and were channeled towards us via the process of accretion.This alignment is dominated by the satellites that enter during the early infall era, i.e $z>0.7$., Comment: 10 pages, 11 figures, accepted for publication in MNRAS

Patchy cosmic reionization resulted in the ionizing UV background asynchronous rise across the Universe. The latter might have left imprints visible in present day observations. Several numerical simulation-based studies show correlations between reionization time and overdensities and object masses today. To remove the mass from the study, as it may not be the sole important parameter, this paper focuses solely on the properties of paired halos within the same mass range as the Milky Way. For this purpose, it uses CoDaII, a fully-coupled radiation hydrodynamics reionization simulation of the local Universe. This simulation holds a halo pair representing the Local Group, in addition to other pairs, sharing similar mass, mass ratio, distance separation and isolation criteria but in other environments, alongside isolated halos within the same mass range. Investigations of the paired halo reionization histories reveal a wide diversity although always inside-out given our reionization model. Within this model, halos in a close pair tend to be reionized at the same time but being in a pair does not bring to an earlier time their mean reionization. The only significant trend is found between the total energy at z = 0 of the pairs and their mean reionization time: pairs with the smallest total energy (bound) are reionized up to 50 Myr earlier than others (unbound). Above all, this study reveals the variety of reionization histories undergone by halo pairs similar to the Local Group, that of the Local Group being far from an average one. In our model, its reionization time is ~625 Myr against 660+/-4 Myr (z~8.25 against 7.87+/-0.02) on average., Comment: Accepted for publication in MNRAS, 11 pages, 5 figures

Using the ASTRAEUS (semi-numerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dark matter simUlationS) framework, we explore the impact of environmental density and radiative feedback on the assembly of galaxies and their host halos during the Epoch of Reionization. The ASTRAEUS framework allows us to study the evolution of galaxies with masses ($\rm 10^{8.2}M_\odot < M_{\rm h} < 10^{13}M_\odot$) in wide variety of environment ($-0.5 < {\rm log}(1+\delta) < 1.3$ averaged over $(2~{\rm cMpc})^3$). We find that : (i) there exists a mass- and redshift- dependent "characteristic" environment (${\rm log} (1+\delta_a(M_{\rm h}, z)) = 0.021\times (M_{\rm h}/M_\odot)^{0.16} + 0.07 z -1.12$, up to $z\sim 10$) at which galaxies are most efficient at accreting dark matter, e.g at a rate of $0.2\%$ of their mass every Myr at $z=5$; (ii) the number of minor and major mergers and their contributions to the dark matter assembly increases with halo mass at all redshifts and is mostly independent of the environment; (iii) at $z=5$ minor mergers contribute slightly more (by up to $\sim 10\%$) to the dark matter assembly while for the stellar assembly, major mergers dominate the contribution from minor mergers for $M_{\rm h}\lesssim 10^{11.5}M_\odot$ galaxies; (iv) radiative feedback quenches star formation more in low-mass galaxies ($M_{\rm h} \lesssim 10^{9.5}M_\odot$) in over-dense environments (${\rm log}(1+\delta) > 0.5$); dominated by their major branch, this yields star formation histories biased towards older ages with a slower redshift evolution., Comment: 17 pages, 15 figures, submitted to MNRAS, comments welcome !

The increase in the observed volume in cosmological surveys imposes various challenges on simulation preparations. Firstly, the volume of the simulations required increases proportionally to the observations. However, large-volume simulations are quickly becoming computationally intractable. Secondly, on-going and future large-volume survey are targeting smaller objects, e.g. emission line galaxies, compared to the earlier focus, i.e. luminous red galaxies. They require the simulations to have higher mass resolutions. In this work we present a machine learning (ML) approach to calibrate the halo catalogue of a low-resolution (LR) simulation by training with a paired high-resolution (HR) simulation with the same background white noise, thus we can build the training data by matching HR haloes to LR haloes in a one-to-one fashion. After training, the calibrated LR halo catalogue reproduces the mass-clustering relation for mass down to $2.5\times 10^{11}~h^{-1}M_\odot$ within $5~{\rm per~cent}$ at scales $k<1~h\,\rm Mpc^{-1}$. We validate the performance of different statistics including halo mass function, power spectrum, two-point correlation function, and bispectrum in both real and redshift space. Our approach generates high-resolution-like halo catalogues ($>200$ particles per halo) from low-resolution catalogues ($>25$ particles per halo) containing corrected halo masses for each object. This allows to bypass the computational burden of a large-volume real high-resolution simulation without much compromise in the mass resolution of the result. The cost of our ML approach ($\sim 1$ CPU-hour) is negligible compared to the cost of a $N$-body simulation (e.g. millions of CPU-hours), The required computing time is cut a factor of 8., Comment: 15 pages, 13 figures, Accepted for publication in MNRAS

Recent determinations of the mean free path of ionising photons (mfp) in the intergalactic medium (IGM) at $\rm z=6$ are lower than many theoretical predictions. To gain insight into this issue, we investigate the evolution of the mfp in our new massive fully coupled radiation hydrodynamics cosmological simulation of reionization: Cosmic Dawn III (CoDaIII). CoDaIII's scale ($\rm 94^3 \, cMpc^3$) and resolution ($\rm 8192^3$ grid) make it particularly suitable to study the evolution of the IGM during Reionization. The simulation was performed with RAMSES-CUDATON on Summit, and used 131072 processors coupled to 24576 GPUs, making it the largest Reionization simulation, and largest RAMSES simulation ever performed. A superior agreement with global constraints on Reionization is obtained in CoDaIII over CoDaII especially for the evolution of the neutral hydrogen fraction and the cosmic photo-ionization rate, thanks to an improved calibration, later end of reionization ($\rm z=5.6$), and higher spatial resolution. Analyzing the mfp, we find that CoDaIII reproduces the most recent observations very well, from $\rm z=6$ to $\rm z=4.6$. We show that the distribution of the mfp in CoDaIII is bimodal, with short (neutral) and long (ionized) mfp modes, respectively, due to the patchiness of reionization and the co-existence of neutral versus ionized regions during Reionization. The neutral mode peaks at sub-kpc to kpc scales of mfp, while the ionized mode peak evolves from $\rm 0.1 Mpc/h$ at $\rm z=7$ to $\sim 10$ Mpc/h at $\rm z=5.2$. Computing the mfp as the average of the ionized mode provides the best match to the recent observational determinations. The distribution reduces to a single neutral (ionized) mode at $\rm z>13$ ($\rm z<5$)., Comment: Submitted to MNRAS Accepted 2022 August 08. Received 2022 August 05; in original form 2022 February 11

In this work, the sixth of a series, we use the ASTRAEUS (semi-numerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dark matter simUlationS) 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 semi-numerical reionization scheme, to include a model for black hole formation, growth, and the production of ionizing radiation from associated AGN (active galactic nuclei). 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% of the ionizing budget at z=6 and starting to play a role only below z < 5., Comment: 15+2 pages, 14+1 figures, 1 table. Accepted in MNRAS

More than 50 per cent of present-day massive disc galaxies show a rotating stellar bar. Their formation and dynamics have been widely studied both numerically and observationally. Although numerical simulations in the $\Lambda$CDM cosmological framework predict the formation of such stellar components, there seems to be a tension between theoretical and observational results. Simulated bars are typically larger in size and have slower pattern speed than observed ones. We study the formation and evolution of barred galaxies, using two $\Lambda$CDM zoom-in hydrodynamical simulations of the CLUES project that follow the evolution of a cosmological Local Group-like volume. We found that our simulated bars, at $z = 0$, are both shorter and faster rotators than previous ones found in other studies on cosmological simulations alleviating the tension mentioned above. These bars match the short tail-end of the observed bar length distribution. In agreement with previous numerical works, we find that bars form in those systems where the disc self-gravity is dominant over the dark matter halo, making them unstable against bar formation. Our bars developed in the last 3-4 Gyr until they achieve their current length and strength; as bars grow, their lengths increase while their rotation speeds decrease. Despite this slowdown, at redshift $z = 0$ their rotation speeds and size match well the observational data., Comment: 11 pages, 9 figures. Accepted for publication in MNRAS

In this work, we have implemented a detailed physical model of galaxy chemical enrichment into the ${\it Astraeus}$ (semi-numerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dark matter simUlationS) framework which couples galaxy formation and reionization in the first billion years. Simulating galaxies spanning over 2.5 orders of magnitude in halo mass with $M_h \sim 10^{8.9-11.5} M_\odot$ ($M_h \sim 10^{8.9-12.8} M_\odot$) at $z \sim 10 ~ (5)$, we find: (i) smooth-accretion of metal-poor gas from the intergalactic medium (IGM) plays a key role in diluting the interstellar medium (ISM) metallicity which is effectively restored due to self-enrichment from star formation; (ii) a redshift averaged gas-mass loading factor that depends on the stellar mass as $\eta_g \approx 1.38 ({M_*}/{10^{10} M_\odot})^{-0.43}$; (iii) the mass-metallicity relation is already in place at $z \sim 10$ and shows effectively no redshift evolution down to $z \sim 5$; (iv) for a given stellar mass, the metallicity decreases with an increase in the star formation rate (SFR); (v) the key properties of the gas-phase metallicity (in units of 12+log(O/H), stellar mass, SFR and redshift are linked through a high-redshift fundamental plane of metallicity (HFPZ) for which we provide a functional form; (vi) the mass-metallicity-SFR relations are effectively independent of the reionization radiative feedback model for $M_* \geq 10^{6.5} M_\odot$ galaxies; (vii) while low-mass galaxies ($M_h \leq 10^9 M_\odot$) are the key contributors to the metal budget of the IGM at early times, higher mass halos provide about 50% of the metal budget at lower-redshifts., Comment: Accepted to MNRAS