Search

Galaxy clusters are located in the densest areas of the universe and are intricately connected to larger structures through the filamentary network of the Cosmic Web. In this scenario, matter flows from areas of lower density to higher density. As a result, the properties of galaxy clusters are deeply influenced by the filaments that are attached to them, which are quantified by a parameter known as connectivity. We explore the dependence of gas-traced filaments connected to galaxy clusters on the mass and dynamical state of the cluster. Moreover, we evaluate the effectiveness of the cosmic web extraction procedure from the gas density maps of simulated cluster regions. Using the DisPerSE cosmic web finder, we identify filamentary structures from 3D gas particle distribution in 324 simulated regions of $30 \, h^{-1}$ Mpc side from The Three Hundred hydrodynamical simulation at redshifts z=0, 1, and 2. We estimate the connectivity at various apertures for $\sim3000$ groups and clusters spanning a mass range from $10^{13} \, h^{-1} \, M_{\odot}$ to $10^{15} \, h^{-1} \, M_{\odot}$. Relationships between connectivity and cluster properties like radius, mass, dynamical state and hydrostatic mass bias are explored. We show that the connectivity is strongly correlated with the mass of galaxy clusters, with more massive clusters being on average more connected. This finding aligns with previous studies in literature, both from observational and simulated data sets. Additionally, we observe a dependence of the connectivity on the aperture at which it is estimated. We find that connectivity decreases with cosmic time, while no dependencies on the dynamical state and hydrostatic mass bias of the cluster are found. Lastly, we observe a significant agreement between the connectivity measured from gas-traced and mock-galaxies-traced filaments in the simulation., Comment: 12 pages, 11 figures

Galaxy clusters are composed of dark matter, gas and stars. Their dark matter component, which amounts to around 80\% of the total mass, cannot be directly observed but traced by the distribution of diffused gas and galaxy members. In this work, we aim to infer the cluster's projected total mass distribution from mock observational data, i.e. stars, Sunyaev-Zeldovich, and X-ray, by training deep learning models. To this end, we have created a multiview images dataset from {\sc{The Three Hundred}} simulation that is optimal for training Machine Learning models. We further study deep learning architectures based on the U-Net to account for single-input and multi-input models. We show that the predicted mass distribution agrees well with the true one., Comment: To appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

Emission line galaxies (ELGs) are targeted by the new generation of spectroscopic surveys to make unprecedented measurements in cosmology from their distribution. Accurately interpreting this data requires understanding the imprints imposed by the physics of galaxy formation and evolution on galaxy clustering. In this work we utilize a semi-analytical model of galaxy formation (SAGE) to explore the necessary components for accurately reproducing the clustering of ELGs. We focus on developing a Halo Occupation Distribution (HOD) prescription able to reproduce the clustering of SAGE galaxies. Typically, HOD models assume that satellite and central galaxies of a given type are independent events. We investigate the need for conformity, i.e. whether the average satellite occupation depends on the existence of a central galaxy of a given type. Incorporating conformity into HOD models is crucial for reproducing the clustering in the reference galaxy sample. Another aspect we investigate is the radial distribution of satellite galaxies within haloes. The traditional density profile models, NFW and Einasto profiles, fail to accurately replicate the small-scale clustering measured for SAGE satellite galaxies. To overcome this limitation, we propose a generalization of the NFW profile, thereby enhancing our understanding of galaxy clustering., Comment: 17 pages, 12 figures and 5 tables

We present the largest full N-body simulation to date with local primordial non-Gaussianities (L-PNG), the \textsc{PNG-UNITsim}. It tracks the evolution of $4096^3$ particles within a periodic box with $L_{\rm box} = 1 \; h^{-1}\,{\rm Gpc}$, leading to a mass resolution of $m_{p} = 1.24\times 10^{9}\; h^{-1}\,M_\odot$. This is enough to resolve galaxies targeted by stage-IV spectroscopic surveys. The \textsc{PNG-UNIT} has \textit{Fixed} initial conditions whose phases are also \textit{Matched} to the pre-existing \textsc{UNIT} simulation. These two features in the simulations reduce our uncertainty significantly so we use 100 \textsc{FastPM} mocks to estimate this reduction. The amplitude of the non-Gaussianities used to set the initial conditions of this new simulation is $f_{\rm NL}^{\rm local} = 100$. In this first study, we use mass selected dark matter haloes from the \textsc{PNG-UNIT} simulation to constrain the local PNG parameters. PNG induce a scale dependent bias, parameterised through \bp or $p$, which might depend on the type of cosmological tracer. Those cases when $p=1$ are referred to as the {\it universality relation}. We measure $p$ as a function of the halo mass. Haloes with masses between $1\times 10^{12}$ and $2\times 10^{13} \, h^{-1} M_\odot$ are well described by the {\it universality relation}. For haloes with masses between $2\times 10^{10}$ and $1\times 10^{12} \, h^{-1} M_\odot$ we find that $p<1$ at $3\sigma$. Combining all the mass bins, we find $p$ consistent with a value of $0.955\pm0.013$, which is $3\sigma$ away from \textit{universality}, as low mass haloes are more numerous. We also study the effect of using priors on $p$ when constraining $f_{\rm NL}$. Using the values we obtain for $b_\phi$ as priors, we forecast that a DESI-like (stage-IV) survey will be able to constrain $f_{\rm NL}$ better than if the universality relation is assumed., Comment: 19 pages, 11 figures. Accepted for publication in A&A

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

A galaxy cluster as the most massive gravitationally-bound object in the Universe, is dominated by Dark Matter, which unfortunately can only be investigated through its interaction with the luminous baryons with some simplified assumptions that introduce an un-preferred bias. In this work, we, {\it for the first time}, propose a deep learning method based on the U-Net architecture, to directly infer the projected total mass density map from idealised observations of simulated galaxy clusters at multi-wavelengths. The model is trained with a large dataset of simulated images from clusters of {\sc The Three Hundred Project}. Although Machine Learning (ML) models do not depend on the assumptions of the dynamics of the intra-cluster medium, our whole method relies on the choice of the physics implemented in the hydrodynamic simulations, which is a limitation of the method. Through different metrics to assess the fidelity of the inferred density map, we show that the predicted total mass distribution is in very good agreement with the true simulated cluster. Therefore, it is not surprising to see the integrated halo mass is almost unbiased, around 1 per cent for the best result from multiview, and the scatter is also very small, basically within 3 per cent. This result suggests that this ML method provides an alternative and more accessible approach to reconstructing the overall matter distribution in galaxy clusters, which can complement the lensing method., Comment: 15 pages, 13 figures, published in MNRAS

Halo Occupation Distribution (HOD) models help us to connect observations and theory, by assigning galaxies to dark matter haloes. In this work we study one of the components of HOD models: the probability distribution function (PDF), which is used to assign a discrete number of galaxies to a halo, given a mean number of galaxies. For satellite galaxies, the most commonly used PDF is a Poisson Distribution. PDFs with super-Poisson variances have also been studied, allowing for continuous values of variances. This has not been the case for sub-Poisson variances, for which only the Nearest Integer distribution, with a single variance, has been used in the past. In this work we propose a distribution based on the binomial one, which provides continuous sub-Poisson variances. We have generated mock galaxy catalogues from two dark-matter only simulations, UNIT and OUTERIM, with HOD models assuming different PDFs. We show that the variance of the PDF for satellite galaxies affects the one-halo term of the projected correlation function, and the Count-In-Cells (CIC) one point statistics. We fit the clustering of eBOSS Emission Line Galaxies, finding a preference for a sub-poissonian PDF, when we only vary the parameter controlling the PDF variance and the fraction of satellites. Using a mock catalogue as a reference, we have also included both the clustering and CIC to constrain the parameters of the HOD model. CIC can provide strong constraints to the PDF variance of satellite galaxies., Comment: Submitted to MNRAS, 19 pages

Using cosmological simulations of galaxy cluster regions from The Three Hundred project we study the nature of gas in filaments feeding massive clusters. By stacking the diffuse material of filaments throughout the cluster sample, we measure average gas properties such as density, temperature, pressure, entropy and Mach number and construct one-dimensional profiles for a sample of larger, radially-oriented filaments to determine their characteristic features as cosmological objects. Despite the similarity in velocity space between the gas and dark matter accretion patterns onto filaments and their central clusters, we confirm some differences, especially concerning the more ordered radial velocity dispersion of dark matter around the cluster and the larger accretion velocity of gas relative to dark matter in filaments. We also study the distribution of shocked gas around filaments and galaxy clusters, showing that the surrounding shocks allow an efficient internal transport of material, suggesting a laminar infall. The stacked temperature profile of filaments is typically colder towards the spine, in line with the cosmological rarefaction of matter. Therefore, filaments are able to isolate their inner regions, maintaining lower gas temperatures and entropy. Finally, we study the evolution of the gas density-temperature phase diagram of our stacked filament, showing that filamentary gas does not behave fully adiabatically through time but it is subject to shocks during its evolution, establishing a characteristic z = 0, entropy-enhanced distribution at intermediate distances from the spine of about 1 - 2 $h^{-1}$ Mpc for a typical galaxy cluster in our sample., Comment: 16 pages, 13 figures. Accepted for publication in MNRAS

The study of the morphology of 2D projected maps of galaxy clusters is a suitable approach to infer, from real data, the dynamical state of those systems. We recently developed a new method to recover the morphological features in galaxy cluster maps which consists of an analytical modelling through the Zernike polynomials. After the first validation of this approach on a set of high-resolution mock maps of the Compton parameter, $y$, from hydrodynamically simulated galaxy clusters in THE THREE HUNDRED project, we apply the Zernike modelling on $y$-maps of local ($z < 0.1$) galaxy clusters observed by the $Planck$ satellite. With a single parameter collecting the main information of the Zernike modelling, we classify their morphology. A set of mock $Planck$-like $y$-maps, generated from THE THREE HUNDRED clusters, is also used to validate our indicator with a proper dynamical state classification. This approach allows us to test the efficiency of the Zernike morphological modelling in evaluating the dynamical population in the real $Planck$ sample., Comment: To appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

Galaxy clusters and their filamentary outskirts reveal useful laboratories to test cosmological models and investigate Universe composition and evolution. Their environment, in particular the filaments of the Cosmic Web to which they are connected, plays an important role in shaping the properties of galaxy clusters. In this project, we analyse the gas filamentary structures present in 324 regions of The Three Hundred hydrodynamical simulation extracted with the DisPerSE filament finder. We estimate the number of gas filaments globally connected to several galaxy clusters, i.e. the connectivity k, with a mass range of $10^{13} \leq M_{200} \, h^{-1} \, M_{\odot} \leq 10^{15} $ at redshift $z=0$. We study the positive correlation between the connectivity and mass of galaxy clusters. Moreover, we explore the impact of filaments on the dynamical state of clusters, quantified by the degree of relaxation parameter $\chi$., Comment: to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

In this study, we investigate a recent finding based on strong lensing observations, which suggests that the sub-halos observed in clusters exhibit greater compactness compared to those predicted by $\Lambda$CDM simulations. To address this discrepancy, we performed a comparative analysis by comparing the cumulative mass function of sub-halos and the $M_{\text{sub}}$-$V_{\text{circ}}$ relation between observed clusters and 324 simulated clusters from The Three Hundred project, focusing on re-simulations using GADGET-X and GIZMO-SIMBA baryonic models. The sub-halos' cumulative mass function of the GIZMO-SIMBA simulated clusters agrees with observations, while the GADGET-X simulations exhibit discrepancies in the lower sub-halo mass range possibly due to its strong SuperNova feedback. Both GADGET-X and GIZMO-SIMBA simulations demonstrate a redshift evolution of the sub-halo mass function and the $V_{max}$ function, with slightly fewer sub-halos observed at lower redshifts. Neither the GADGET-X nor GIZMO-SIMBA(albeit a little closer) simulated clusters' predictions for the $M_{\text{sub}}$-$V_{\text{circ}}$ relation align with the observational result. Further investigations on the correlation between sub-halo/halo properties and the discrepancy in the $M_{\text{sub}}$-$V_{\text{circ}}$ relation reveals that the sub-halo's half mass radius and galaxy stellar age, the baryon fraction and sub-halo distance from the cluster's centre, as well as the halo relaxation state play important roles on this relation. Nevertheless, we think it is still challenging in accurately reproducing the observed $M_{\text{sub}}$-$V_{\text{circ}}$ relation in our current hydrodynamic cluster simulation under the standard $\Lambda$CDM cosmology., Comment: see the companion paper, Meneghetti et al. (2023), on today's arxiv list

Previous studies have revealed that the estimated probability of galaxy-galaxy strong lensing in observed galaxy clusters exceeds the expectations from the $\Lambda$ Cold Dark Matter cosmological model by one order of magnitude. We aim to understand the origin of this excess by analyzing a larger set of simulated galaxy clusters and investigating how the theoretical expectations vary under different adopted prescriptions and numerical implementations of star formation and feedback in simulations. We perform a ray-tracing analysis of 324 galaxy clusters from the Three Hundred project, comparing the Gadget-X and Gizmo-Simba runs. These simulations, which start from the same initial conditions, are performed with different implementations of hydrodynamics and galaxy formation models tailored to match different observational properties of the Intra-Cluster-Medium and cluster galaxies. We find that galaxies in the Gizmo-Simba simulations develop denser stellar cores than their Gadget-X counterparts. Consequently, their probability for galaxy-galaxy strong lensing is higher by a factor of $\sim 3$. This increment is still insufficient to fill the gap with observations, as a discrepancy by a factor $\sim 4$ still persists. In addition, we find that several simulated galaxies have Einstein radii that are too large compared to observations. We conclude that a persistent excess of galaxy-galaxy strong lensing exists in observed galaxy clusters. The origin of this discrepancy with theoretical predictions is still unexplained in the framework of the cosmological hydrodynamical simulations. This might signal a hitherto unknown issue with either the simulation methods or our assumptions regarding the standard cosmological model., Comment: 6 pages, 2 figures, accepted for publication on A&A Letters. See companion paper Srivastava et al. (2023)

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

Using the data set of The Three Hundred project, i.e. 324 hydrodynamical resimulations of cluster-sized haloes and the regions of radius 15 $h^{-1}$Mpc around them, we study galaxy pairs in high-density environments. By projecting the galaxies' 3D coordinates onto a 2D plane, we apply observational techniques to find galaxy pairs. Based on a previous theoretical study on galaxy groups in the same simulations, we are able to classify the observed pairs into "true" or "false", depending on whether they are gravitationally bound or not. We find that the fraction of true pairs (purity) crucially depends on the specific thresholds used to find the pairs, ranging from around 30 to more than 80 per cent in the most restrictive case. Nevertheless, in these very restrictive cases, we see that the completeness of the sample is low, failing to find a significant number of true pairs. Therefore, we train a machine learning algorithm to help us to identify these true pairs based on the properties of the galaxies that constitute them. With the aid of the machine learning model trained with a set of properties of all the objects, we show that purity and completeness can be boosted significantly using the default observational thresholds. Furthermore, this machine learning model also reveals the properties that are most important when distinguishing true pairs, mainly the size and mass of the galaxies, their spin parameter, gas content and shape of their stellar components., Comment: 20 pages, 8 figures. Published in MNRAS

When the first galaxies formed and starlight escaped into the intergalactic medium to reionize it, galaxy formation and reionization were both highly inhomogeneous in time and space, and fully-coupled by mutual feedback. To show how this imprinted the UV luminosity function (UVLF) of reionization-era galaxies, we use our large-scale, radiation-hydrodynamics simulation CoDa II to derive the time- and space-varying halo mass function and UVLF, from $z\simeq6$-15. That UVLF correlates strongly with local reionization redshift: earlier-reionizing regions have UVLFs that are higher, more extended to brighter magnitudes, and flatter at the faint end than later-reionizing regions observed at the same $z$. In general, as a region reionizes, the faint-end slope of its local UVLF flattens, and, by $z=6$ (when reionization ended), the global UVLF, too, exhibits a flattened faint-end slope, `rolling-over' at $M_\text{UV}\gtrsim-17$. CoDa II's UVLF is broadly consistent with cluster-lensed galaxy observations of the Hubble Frontier Fields at $z=6$-8, including the faint end, except for the faintest data point at $z=6$, based on one galaxy at $M_\text{UV}=-12.5$. According to CoDa II, the probability of observing the latter is $\sim5\%$. However, the effective volume searched at this magnitude is very small, and is thus subject to significant cosmic variance. We find that previous methods adopted to calculate the uncertainty due to cosmic variance underestimated it on such small scales by a factor of 2-4, primarily by underestimating the variance in halo abundance when the sample volume is small., Comment: 16 pages, 16 figures, accepted by MNRAS 07/20/23, comments welcome

We use a set of about 300 simulated clusters from The Three Hundred Project to calculate their hydrostatic masses and evaluate the associated bias by comparing them with the true cluster mass. Over a redshift range from 0.07 to 1.3, we study the dependence of the hydrostatic bias on redshift, concentration, mass growth, dynamical state, mass, and halo shapes. We find almost no correlation between the bias and any of these parameters. However, there is a clear evidence that the scatter of the mass-bias distribution is larger for low-concentrated objects, high mass growth, and more generically for disturbed systems. Moreover, we carefully study the evolution of the bias of twelve clusters throughout a major-merger event. We find that the hydrostatic-mass bias follows a particular evolution track along the merger process: to an initial significant increase of the bias recorded at the begin of merger, a constant plateaus follows until the end of merge, when there is a dramatic decrease in the bias before the cluster finally become relaxed again. This large variation of the bias is in agreement with the large scatter of the hydrostatic bias for dynamical disturbed clusters. These objects should be avoided in cosmological studies because their exact relaxation phase is difficult to predict, hence their mass bias cannot be trivially accounted for., Comment: 11 pages, 8 figures. Accepted for publication in MNRAS

In this work, we combine the semi-analytic model of galaxy formation and evolution SAG with the $102$ relaxed simulated galaxy clusters from The Three Hundred project, and we study the link between the quenching of star formation (SF) and the physical processes that galaxies experience through their dynamical history in and around clusters. We classify galaxies in four populations based on their orbital history: recent and ancient infallers, and backsplash and neighbouring galaxies. We find that $\sim 85$ per cent of the current population of quenched galaxies located inside the clusters are ancient infallers with low or null content of hot and cold gas. The fraction of quenched ancient infallers increases strongly between the first and second pericentric passage, due to the removal of hot gas by the action of ram-pressure stripping (RPS). The majority of them quenches after the first pericentric passage, but a non-negligible fraction needs a second passage, specially galaxies with $M_\star \leq 10^{10.5} \, {\rm M_\odot}$. Recent infallers represent $\sim 15$ per cent of the quenched galaxies located inside the cluster and, on average, they contain a high proportion of hot and cold gas; moreover, pre-processing effects are the responsible for quenching the recent infallers prior to infall onto the main cluster progenitor. The $\sim 65$ per cent of quenched galaxies located around clusters are backsplash galaxies, for which the combination of RPS acting during a pre-processing stage and inside the cluster is necessary for the suppression of SF in this population., Comment: 23 pages, 13 figures + Supplementary material. Accepted for publication in MNRAS

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

Galaxy clusters grow by accreting galaxies as individual objects, or as members of a galaxy group. These groups can strongly impact galaxy evolution, stripping the gas from galaxies, and enhancing the rate of galaxy mergers. However, it is not clear how the dynamics and structure of groups are affected when they interact with a large cluster, or whether all group members necessarily experience the same evolutionary processes. Using data from TheThreeHundred project, a suite of 324 hydrodynamical resimulations of large galaxy clusters, we study the properties of 1340 groups passing through a cluster. We find that half of group galaxies become gravitationally unbound from the group by the first pericentre, typically just 0.5-1 Gyr after cluster entry. Most groups quickly mix with the cluster satellite population; only 8% of infalling group haloes later leave the cluster, although for nearly half of these, all of their galaxies have become unbound, tidally disrupted or merged into the central by this stage. The position of galaxies in group-centric phase space is also important -- only galaxies near the centre of a group ($r\lesssim0.7R_{200}$) remain bound once a group is inside a cluster, and slow-moving galaxies in the group centre are likely to be tidally disrupted, or merge with another galaxy. This work will help future observational studies to constrain the environmental histories of group galaxies. For instance, groups observed inside or nearby to clusters have likely approached very recently, meaning that their galaxies will not have experienced a cluster environment before., Comment: 21 pages, 11 figures, published in MNRAS. v2: fixed inconsistent notation

Galaxy clusters are useful laboratories to investigate the evolution of the Universe, and accurately measuring their total masses allows us to constrain important cosmological parameters. However, estimating mass from observations that use different methods and spectral bands introduces various systematic errors. This paper evaluates the use of a Convolutional Neural Network (CNN) to reliably and accurately infer the masses of galaxy clusters from the Compton-y parameter maps provided by the Planck satellite. The CNN is trained with mock images generated from hydrodynamic simulations of galaxy clusters, with Planck's observational limitations taken into account. We observe that the CNN approach is not subject to the usual observational assumptions, and so is not affected by the same biases. By applying the trained CNNs to the real Planck maps, we find cluster masses compatible with Planck measurements within a 15% bias. Finally, we show that this mass bias can be explained by the well known hydrostatic equilibrium assumption in Planck masses, and the different parameters in the Y500-M500 scaling laws. This work highlights that CNNs, supported by hydrodynamic simulations, are a promising and independent tool for estimating cluster masses with high accuracy, which can be extended to other surveys as well as to observations in other bands., Comment: 17 pages, 3 Figures and Supplementary material (+11 figures). Published in Nature Astronomy

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

Close pairs of galaxies have been broadly studied in the literature as a way to understand galaxy interactions and mergers. In observations they are usually defined by setting a maximum separation in the sky and in velocity along the line of sight, and finding galaxies within these ranges. However, when observing the sky, projection effects can affect the results, by creating spurious pairs that are not close in physical distance. In this work we mimic these observational techniques to find pairs in The Three Hundred simulations of clusters of galaxies. The galaxies' 3D coordinates are projected into 2D, with Hubble flow included for their line-of-sight velocities. The pairs found are classified into "good" or "bad" depending on whether their 3D separations are within the 2D spatial limit or not. We find that the fraction of good pairs can be between 30 and 60 per cent depending on the thresholds used in observations. Studying the ratios of observable properties between the pair member galaxies, we find that the likelihood of a pair being "good" can be increased by around 40, 20 and 30 per cent if the given pair has, respectively, a mass ratio below 0.2, metallicity ratio above 0.8, or colour ratio below 0.8. Moreover, shape and stellar-to-halo mass ratios respectively below 0.4 and 0.2 can increase the likelihood by 50 to 100 per cent. These results suggest that these properties can be used to increase the chance of finding good pairs in observations of galaxy clusters and their environment., Comment: 15 pages, 10 figures. Published in MNRAS; typos corrected

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 star formation history (SFH) of galaxies is critical for understanding galaxy evolution. Hydrodynamical simulations enable us to precisely reconstruct the SFH of galaxies and establish a link to the underlying physical processes. In this work, we present a model to describe individual galaxies' SFHs from three simulations: TheThreeHundred, Illustris-1 and TNG100-1. This model divides the galaxy SFH into two distinct components: the "main sequence" and the "variation". The "main sequence" part is generated by tracing the history of the $SFR-M_*$ main sequence of galaxies across time. The "variation" part consists of the scatter around the main sequence, which is reproduced by fractional Brownian motions. We find that: 1) The evolution of the main sequence varies between simulations; 2) fractional Brownian motions can reproduce many features of SFHs, however, discrepancies still exist; 3) The variations and mass-loss rate are crucial for reconstructing the SFHs of the simulations. This model provides a fair description of the SFHs in simulations. On the other hand, by correlating the fractional Brownian motion model to simulation data, we provide a 'standard' against which to compare simulations., Comment: 22 pages, 18 figures, 2 tables

The Sunyaev-Zeldovich (SZ) effect is a powerful tool in modern cosmology. With future observations promising ever improving SZ measurements, the relativistic corrections to the SZ signals from galaxy groups and clusters are increasingly relevant. As such, it is important to understand the differences between three temperature measures: (a) the average relativistic SZ (rSZ) temperature, (b) the mass-weighted temperature relevant for the thermal SZ (tSZ) effect, and (c) the X-ray spectroscopic temperature. In this work, we compare these cluster temperatures, as predicted by the {\sc Bahamas} \& {\sc Macsis}, {\sc Illustris-TNG}, {\sc Magneticum}, and {\sc The Three Hundred Project} simulations. Despite the wide range of simulation parameters, we find the SZ temperatures are consistent across the simulations. We estimate a $\simeq 10\%$ level correction from rSZ to clusters with $Y\simeq10^{-4}$~Mpc$^{-2}$. Our analysis confirms a systematic offset between the three temperature measures; with the rSZ temperature $\simeq 20\%$ larger than the other measures, and diverging further at higher redshifts. We demonstrate that these measures depart from simple self-similar evolution and explore how they vary with the defined radius of haloes. We investigate how different feedback prescriptions and resolution affect the observed temperatures, and discover the SZ temperatures are rather insensitive to these details. The agreement between simulations indicates an exciting avenue for observational and theoretical exploration, determining the extent of relativistic SZ corrections. We provide multiple simulation-based fits to the scaling relations for use in future SZ modelling., Comment: Accepted by MNRAS

In this paper we study the applicability of a set of supervised machine learning (ML) models specifically trained to infer observed related properties of the baryonic component (stars and gas) from a set of features of dark matter only cluster-size halos. The training set is built from THE THREE HUNDRED project which consists of a series of zoomed hydrodynamical simulations of cluster-size regions extracted from the 1 Gpc volume MultiDark dark-matter only simulation (MDPL2). We use as target variables a set of baryonic properties for the intra cluster gas and stars derived from the hydrodynamical simulations and correlate them with the properties of the dark matter halos from the MDPL2 N-body simulation. The different ML models are trained from this database and subsequently used to infer the same baryonic properties for the whole range of cluster-size halos identified in the MDPL2. We also test the robustness of the predictions of the models against mass resolution of the dark matter halos and conclude that their inferred baryonic properties are rather insensitive to their DM properties which are resolved with almost an order of magnitude smaller number of particles. We conclude that the ML models presented in this paper can be used as an accurate and computationally efficient tool for populating cluster-size halos with observational related baryonic properties in large volume N-body simulations making them more valuable for comparison with full sky galaxy cluster surveys at different wavelengths. We make the best ML trained model publicly available., Comment: 20 pages, 10 figures, published in MNRAS

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

We introduce \textsc{Gizmo-Simba}, a new suite of galaxy cluster simulations within \textsc{The Three Hundred} project. \textsc{The Three Hundred} consists of zoom re-simulations of 324 clusters with $M_{200}\gtrsim 10^{14.8}M_\odot$ drawn from the MultiDark-Planck $N$-body simulation, run using several hydrodynamic and semi-analytic codes. The \textsc{Gizmo-Simba} suite adds a state-of-the-art galaxy formation model based on the highly successful {\sc Simba} simulation, mildly re-calibrated to match $z=0$ cluster stellar properties. Comparing to \textsc{The Three Hundred} zooms run with \textsc{Gadget-X}, we find intrinsic differences in the evolution of the stellar and gas mass fractions, BCG ages, and galaxy colour-magnitude diagrams, with \textsc{Gizmo-Simba} generally providing a good match to available data at $z \approx 0$. \textsc{Gizmo-Simba}'s unique black hole growth and feedback model yields agreement with the observed BH scaling relations at the intermediate-mass range and predicts a slightly different slope at high masses where few observations currently lie. \textsc{Gizmo-Simba} provides a new and novel platform to elucidate the co-evolution of galaxies, gas, and black holes within the densest cosmic environments., Comment: 20 pages, 12 figures, matched to published version 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

Using the data set of The Three Hundred project, i.e. a suite of 324 hydrodynamical resimulations of cluster-sized haloes, we study galaxy cluster mergers and their effect on colour and luminosity changes of their brightest cluster galaxies (BCG). We track the main progenitor of each halo at z=0 and search for merger situations based on its mass accretion history, defining mergers as very rapid increases in the halo mass. Based upon the evolution of the dynamical state of the cluster we define a pre- and post-merger phase. We create a list of all these events and statistically study their mass ratio and timescales, with the former verifying that all instances are in fact major mergers. By comparing to a control sample of clusters without mergers, we study the effect mergers have on the stellar component of the BCG. Analysing the mass, age and metallicity of the BCG stellar particles, we find that the stellar content of BCGs grows significantly during mergers and, even though the main growth mechanism is the accretion of older stars, there is even a burst in star formation induced by the merger. In our simulations, BCGs in mergers form in median around 70 per cent more stars than those normally growing, although this depends on the radius considered for defining the BCG. Regarding observable properties, we see an increase in SDSS-u luminosity of 20 per cent during mergers, accompanied by a slightly slower increase of the galaxy g-r colour as compared to the control sample., Comment: 18 pages, 12 figures. Accepted for publication 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

Cosmological $N$-body simulations provide numerical predictions of the structure of the Universe against which to compare data from ongoing and future surveys, but the growing volume of the Universe mapped by surveys requires correspondingly lower statistical uncertainties in simulations, usually achieved by increasing simulation sizes at the expense of computational power. It was recently proposed to reduce simulation variance without incurring additional computational costs by adopting fixed-amplitude initial conditions. This method has been demonstrated not to introduce bias in various statistics, including the two-point statistics of galaxy samples typically used for extracting cosmological parameters from galaxy redshift survey data, but requires us to revisit current methods for estimating covariance matrices of clustering statistics for simulations. In this work, we find that it is not trivial to construct covariance matrices analytically for fixed-amplitude simulations, but we demonstrate that EZmock (Effective Zel'dovich approximation mock catalogue), the most efficient method for constructing mock catalogues with accurate two- and three-point statistics, provides reasonable covariance matrix estimates for such simulations. We further examine how the variance suppression obtained by amplitude-fixing depends on three-point clustering, small-scale clustering, and galaxy bias, and propose intuitive explanations for the effects we observe based on the EZmock bias model., Comment: 9 pages, 7 figures; in v2 we incorporate minor modifications from peer review and copyediting

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

Galaxy clusters grow by accreting galaxies from the field and along filaments of the cosmic web. As galaxies are accreted they are affected by their local environment before they enter (pre-processing), and traverse the cluster potential. Observations that aim to constrain pre-processing are challenging to interpret because filaments comprise a heterogeneous range of environments including groups of galaxies embedded within them and backsplash galaxies that contain a record of their previous passage through the cluster. This motivates using modern cosmological simulations to dissect the population of galaxies found in filaments that are feeding clusters, to better understand their history, and aid the interpretation of observations. We use zoom-in simulations from The ThreeHundred project to track halos through time and identify their environment. We establish a benchmark for galaxies in cluster infall regions that supports the reconstruction of the different modes of pre-processing. We find that up to 45% of all galaxies fall into clusters via filaments (closer than 1Mpc/h from the filament spine). 12% of these filament galaxies are long-established members of groups and between 30 and 60% of filament galaxies at R200 are backsplash galaxies. This number depends on the cluster's dynamical state and sharply drops with distance. Backsplash galaxies return to clusters after deflecting widely from their entry trajectory, especially in relaxed clusters. They do not have a preferential location with respect to filaments and cannot collapse to form filaments. The remaining pristine galaxies (30 - 60%) are environmentally effected by cosmic filaments alone., Comment: 13 pages, 9 figures, accepted 2021 November 21

Several methods are used to evaluate, from observational data, the dynamical state of galaxy clusters. Among them, the morphological analysis of cluster images is well suited for this purpose. We report a new approach to the morphology, which consists in analytically modelling the images with a set of orthogonal functions, the Zernike polynomials (ZPs). We validated the method on mock high-resolution Compton parameter maps of synthetic galaxy clusters from THE THREE HUNDRED project. To classify the maps for their morphology we defined a single parameter, $\mathcal{C}$, by combining the contribution of some ZPs in the modelling. We verify that $\mathcal{C}$ is linearly correlated with a combination of common morphological parameters and also with a proper 3D dynamical-state indicator available for the synthetic clusters we used. We also show the early results of the Zernike modelling applied on Compton parameter maps of local clusters ($z < 0.1$) observed by the $\textit{Planck}$ satellite. At last, we report the preliminary results of this kind of morphological analysis on mock X-ray maps of THE THREE HUNDRED clusters., Comment: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences

Clusters of galaxies mass can be inferred by indirect observations, see X-ray band, Sunyaev-Zeldovich (SZ) effect signal or optical. Unfortunately, all of them are affected by some bias. Alternatively, we provide an independent estimation of the cluster masses from the Planck PLSZ2 catalog of galaxy clusters using a machine-learning method. We train a Convolutional Neural Network (CNN) model with the mock SZ observations from The Three Hundred(the300) hydrodynamic simulations to infer the cluster masses from the real maps of the Planck clusters. The advantage of the CNN is that no assumption on a priory symmetry in the cluster's gas distribution or no additional hypothesis about the cluster physical state are made. We compare the cluster masses from the CNN model with those derived by Planck and conclude that the presence of a mass bias is compatible with the simulation results., Comment: To appear in the Proceedings of the International Conference entitled "mm Universe @NIKA2", Rome(Italy), June 2021, EPJ Web of conferences

Hydrostatic equilibrium (HE) is often used in observations to estimate galaxy clusters masses. We use a set of almost 300 simulated clusters from The Three Hundred Project, to estimate the cluster HE mass and the bias deriving from it. We study the dependence of the bias on several dynamical state indicators across a redshift range from 0.07 to 1.3, finding no dependence between them. Moreover, we focus our attention on the evolution of the HE bias during the merger phase, where the bias even reaches negative values due to an overestimation of the mass with HE., Comment: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences

The Planck Collaboration has shown that the number of clusters as a function of their mass and redshift is an extremely powerful tool for cosmological analyses. However, the true cluster mass is not directly measurable. Among the possible approaches, clusters mass could be related to different observables via self similar scaling law. These observables are related to the baryonic components of which a cluster is composed. However, the theoretical relations that allow the use of these proxies often are affected by observational and physical biases, which impacts on the determination of the cluster mass. Fortunately, cosmological simulations are an extremely powerful tool to assess these problems. We present our calibration of the scaling relation between mass and velocity dispersion of galaxy members from the study of the simulated clusters of \THP{} project with mass above $10^{13} M_\odot$. In order to investigate the presence of a redshift dependence, we analyzed 16 different redshifts between $z = 0$ and $z = 2$. Finally, we investigated the impact of different AGN feedback models., Comment: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences

We use the \textsc{astraeus} framework, that couples an N-body simulation with a semi-analytic model for galaxy formation and a semi-numerical model for reionization, to quantify the star formation histories (SFHs) of galaxies in the first billion years. Exploring four models of radiative feedback, we fit the SFH of each galaxy at $z>5$ as $\mathrm{log}(\mathrm{SFR}(z))=-\alpha(1 + z)+\beta$; star formation is deemed stochastic if it deviates from this fit by more than $\Delta_\mathrm{SFR}=0.6\,$dex. Our key findings are: (i) The fraction of stellar mass formed and time spent in the stochastic phase decrease with increasing stellar mass and redshift $z$. While galaxies with stellar masses of $M_\star\sim10^7M_\odot$ at $z\sim5~(10)$ form $\sim70\%~(20\%)$ of their stellar mass in the stochastic phase, this reduces to $<10\%$ at all redshifts for galaxies with $M_\star > 10^{10}M_\odot$; (ii) the fractional mass assembled and lifetime spent in the stochastic phase do not significantly change with the radiative feedback model used; (iii) at all redshifts, $\alpha$ increases (decreases for the strongest radiative feedback model) with stellar mass for galaxies with $M_\star\lesssim 10^{8.5}M_\odot$ and converges to $\sim0.18$ for more massive galaxies; $\beta$ always increases with stellar mass. Our proposed fits can reliably recover the stellar masses and mass-to-light ratios for galaxies with $M_\star\sim10^{8-10.5}M_\odot$ and $M_{UV}\sim-17~{\rm to}~-23$ at $z\sim 5-9$. This physical model can therefore be used to derive the SFHs for galaxies observed by a number of forthcoming instruments., Comment: 19 pages, 14 figures, accepted for publication in MNRAS

Galaxy cluster masses, rich with cosmological information, can be estimated from internal dark matter (DM) velocity dispersions, which in turn can be observationally inferred from satellite galaxy velocities. However, galaxies are biased tracers of the DM, and the bias can vary over host halo and galaxy properties as well as time. We precisely calibrate the velocity bias, b_v -- defined as the ratio of galaxy and DM velocity dispersions -- as a function of redshift, host halo mass, and galaxy stellar mass threshold (Mstarsat), for massive halos (M200c > 1e13.5 msun) from five cosmological simulations: IllustrisTNG, Magneticum, Bahamas + Macsis, The Three Hundred Project, and MultiDark Planck-2. We first compare scaling relations for galaxy and DM velocity dispersion across simulations; the former is estimated using a new ensemble velocity likelihood method that is unbiased for low galaxy counts per halo, while the latter uses a local linear regression. The simulations show consistent trends of b_v increasing with M200c and decreasing with redshift and Mstarsat. The ensemble-estimated theoretical uncertainty in b_v is 2-3% but becomes percent-level when considering only the three highest resolution simulations. We update the mass-richness normalization previously estimated by Farahi et al. (2016) for an SDSS redMaPPer cluster sample. The improved accuracy of our b_v estimates reduces the mass normalization uncertainty from 22% to 8%, demonstrating that dynamical estimation techniques can be competitive with weak lensing in calibrating population mean masses. We discuss necessary steps for further improving this precision. Our estimates for b_v(M200c, Mstarsat, z) are made publicly available., Comment: [v1]: 8 figures, 14 pages in main text + 5 figures, 3 pages in Appendix. [v2]: Minor text edits. Accepted in MNRAS

Recent observational and theoretical studies of the Local Group (LG) dwarf galaxies have highlighted their unique star-formation history, stellar metallicity, gas content, and kinematics. We investigate the commonality of these features by comparing constrained LG and field central dwarf halo simulations in the Numerical Investigation of a Hundred Astrophysical Objects (NIHAO) project. Our simulations, performed with NIHAO-like hydrodynamics which track the evolution of the Milky Way (MW) and M31 along with ~100 dwarfs in the LG, reveal the total gas mass and stellar properties (velocity dispersion, evolution history, etc.) of present-day LG dwarfs to be similar to field systems. However, relative to field galaxies, LG dwarfs have more cold gas in their central parts and more metal-rich gas in the halo stemming from interactions with other dwarfs living in a high-density environment like the LG. Interestingly, the direct impact of massive MW/M31 analogues on the metallicity evolution of LG dwarfs is minimal; LG dwarfs accrete high-metallicity gas mostly from other dwarfs at late times. We have also tested for the impact of metal diffusion on the chemical evolution of LG dwarfs, and found that it does not affect the stellar or gaseous content of LG dwarfs. Our simulations suggest that the stellar components of LG dwarfs offer a unique and unbiased local laboratory for galaxy-formation tests and comparisons, especially against the overall dwarf population in the Universe., Comment: 16 pages, 12 figures, Accepted for publication in MNRAS

Galaxy clusters have a triaxial matter distribution. The weak-lensing signal, an important part in cosmological studies, measures the projected mass of all matter along the line-of-sight, and therefore changes with the orientation of the cluster. Studies suggest that the shape of the brightest cluster galaxy (BCG) in the centre of the cluster traces the underlying halo shape, enabling a method to account for projection effects. We use 324 simulated clusters at four redshifts between 0.1 and 0.6 from `The Three Hundred Project' to quantify correlations between the orientation and shape of the BCG and the halo. We find that haloes and their embedded BCGs are aligned, with an average $\sim$20 degree angle between their major axes. The bias in weak lensing cluster mass estimates correlates with the orientation of both the halo and the BCG. Mimicking observations, we compute the projected shape of the BCG, as a measure of the BCG orientation, and find that it is most strongly correlated to the weak-lensing mass for relaxed clusters. We also test a 2-dimensional cluster relaxation proxy measured from BCG mass isocontours. The concentration of stellar mass in the projected BCG core compared to the total stellar mass provides an alternative proxy for the BCG orientation. We find that the concentration does not correlate to the weak-lensing mass bias, but does correlate with the true halo mass. These results indicate that the BCG shape and orientation for large samples of relaxed clusters can provide information to improve weak-lensing mass estimates., Comment: 15 pages, 9 figures, Figure 7 is key plot, Updated to match version accepted by MNRAS

We study the clustering of HI intensity maps produced from simulations with a focus on baryonic acoustic oscillations (BAO) and the effects induced by telescope beam smoothing and foreground cleaning. We start by creating a HI catalogue at $z=1.321$ based on the Semi-Analytic Galaxy Evolution (SAGE) model applied to the UNIT simulations. With this catalogue we investigate the relation between model HI and the dark matter haloes and we also study the abundance of HI, $\Omega_{\rm HI}$, predicted by this model. We then create synthetic HI intensity maps with a Nearest-Grid-Point approach. In order to simulate the telescope beam effect, a Gaussian smoothing is applied on the plane perpendicular to the line of sight. The effect of foreground removal methods is simulated by exponentially damping the largest wavelength Fourier modes on the radial direction. We study the anisotropic 2-point correlation function (2PCF) $\xi(r_\perp,r_\parallel)$ and how it is affected by the aforementioned observational effects. In order to better isolate the BAO signal, we study several 2PCF $\mu$-wedges (with a restricted range of orientations $\mu$) tailored to address the systematics effects and we compare them with different definitions of radial 2PCFs. Finally, we discuss our findings in the context of an SKA-like survey, finding a clear BAO signal in most of the estimators here proposed., Comment: published in MNRAS; 18 pages, 10 figures

Using 324 numerically modelled galaxy clusters as provided by THE THREE HUNDRED project, we study the evolution of the kinematic properties of the stellar component of haloes on first infall. We select objects with M$_{\textrm{star}}>5\times10^{10} h^{-1}M_{\odot}$ within $3R_{200}$ of the main cluster halo at $z=0$ and follow their progenitors. We find that although haloes are stripped of their dark matter and gas after entering the main cluster halo, there is practically no change in their stellar kinematics. For the vast majority of our `galaxies' -- defined as the central stellar component found within the haloes that form our sample -- their kinematic properties, as described by the fraction of ordered rotation, and their position in the specific stellar angular momentum$-$stellar mass plane $j_{\rm star}$ -- M$_{\rm star}$, are mostly unchanged by the influence of the central host cluster. However, for a small number of infalling galaxies, stellar mergers and encounters with remnant stellar cores close to the centre of the main cluster, particularly during pericentre passage, are able to spin-up their stellar component by $z=0$., Comment: 8 pages, 5 figures; accepted for publication in A&A

At ~16-17Mpc from us, the Virgo cluster is a formidable source of information to study cluster formation and galaxy evolution in rich environments. Several observationally-driven formation scenarios arose within the past decade to explain the properties of galaxies that entered the cluster recently and the nature of the last significant merger that the cluster underwent. Confirming these scenarios requires extremely faithful numerical counterparts of the cluster. This paper presents the first CLONE, Constrained LOcal and Nesting Environment, simulation of the Virgo cluster within a ~15Mpc radius sphere. This cosmological hydrodynamical simulation, with feedback from supernovae and active galactic nuclei, with a ~3x10^7Msun dark matter particle mass and a minimum cell size of 350pc in the zoom region, reproduces Virgo within its large scale environment unlike a random cluster simulation. Overall the distribution of the simulated galaxy population matches the observed one including M87. The simulated cluster formation reveals exquisite agreements with observationally-driven scenarios: within the last Gigayear, about 300 small galaxies (M*>10^7Msun) entered the cluster, most of them within the last 500Myr. The last significant merger event occurred about 2 Gigayears ago: a group with a tenth of the mass of today's cluster entered from the far side as viewed from the Milky Way. This excellent numerical replica of Virgo will permit studying different galaxy type evolution (jellyfish, backsplash, etc.) as well as feedback phenomena in the cluster core via unbiased comparisons between simulated and observed galaxies and hot gas phase profiles to understand this great physics laboratory., Comment: Accepted for publication in MNRAS, 15 pages, 14 figures, 4 tables

New surveys such as ESA's Euclid mission are planned to map with unprecedented precision the large-scale structure of the Universe by measuring the 3D positions of tens of millions of galaxies. It is necessary to develop theoretically modelled galaxy catalogues to estimate the expected performance and to optimise the analysis strategy of these surveys. We populate two pairs of (1 Gpc/h)^3 volume dark-matter-only simulations from the UNIT project with galaxies using the SAGE semi-analytic model of galaxy formation, coupled to the photoionisation model get_emlines to estimate their H_alpha emission. These catalogues represent a unique suite that includes galaxy formation physics and - thanks to the fixed-pair technique used - an effective volume of ~(5 Gpc/h)^3, which is several times larger than the Euclid survey. We present the performance of these data and create five additional emission-line galaxy (ELG) catalogues by applying a dust attenuation model as well as adjusting the flux threshold as a function of redshift in order to reproduce Euclid-forecast dN/dz values. As a first application, we study the abundance and clustering of those model H_alpha ELGs: for scales greater than ~5 Mpc/h, we find a scale-independent bias with a value of $b\sim 1$ at redshift z ~ 0.5, that can increase nearly linearly to $b\sim 4$ at $z\sim 2$, depending on the ELG catalogue. Model galaxy properties, including their emission-line fluxes (with and without dust extinction) are publicly available., Comment: 16 pages, revised version as accepted by MNRAS