Showing 1–50 of 274 results for author: Yepes, G

The Three Hundred project: Estimating the dependence of gas filaments on the mass of galaxy clusters

Authors: Sara Santoni, Marco De Petris, Gustavo Yepes, Antonio Ferragamo, Matteo Bianconi, Meghan E. Gray, Ulrike Kuchner, Frazer R. Pearce, Weiguang Cui, Stefano Ettori

Abstract: 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… ▽ More 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. △ Less

Submitted 27 May, 2024; originally announced May 2024.

Comments: 12 pages, 11 figures

Generating Galaxy Clusters Mass Density Maps from Mock Multiview Images via Deep Learning

Authors: Daniel de Andres, Weiguang Cui, Gustavo Yepes, Marco De Petris, Gianmarco Aversano, Antonio Ferragamo, Federico De Luca, A. Jiménez Muñoz

Abstract: 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 l… ▽ More 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. △ Less

Submitted 9 April, 2024; v1 submitted 8 April, 2024; originally announced April 2024.

Comments: 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

3D scaling laws and projection effects in The300-NIKA2 Sunyaev-Zeldovich Large Program Twin Samples

Authors: A. Paliwal, W. Cui, D. de Andrés, M. De Petris, A. Ferragamo, C. Hanser, J. -F. Macías-Pérez, F. Mayet, A. Moyer-Anin, M. Muñoz-Echeverría, L. Perotto, E. Rasia, G. Yepes

Abstract: The abundance of galaxy clusters with mass and redshift is a well-known cosmological probe. The cluster mass is a key parameter for studies that aim to constrain cosmological parameters using galaxy clusters, making it critical to understand and properly account for the errors in its estimates. Subsequently, it becomes important to correctly calibrate scaling relations between observables like the… ▽ More The abundance of galaxy clusters with mass and redshift is a well-known cosmological probe. The cluster mass is a key parameter for studies that aim to constrain cosmological parameters using galaxy clusters, making it critical to understand and properly account for the errors in its estimates. Subsequently, it becomes important to correctly calibrate scaling relations between observables like the integrated Compton parameter and the mass of the cluster. The NIKA2 Sunyaev-Zeldovich Large program (LPSZ) enables one to map the intracluster medium profiles in the mm-wavelength band with great details (resolution of $11 \ \mathrm{\&}\ 17^{\prime \prime}$ at $1.2 \ \mathrm{\&}\ 2 $ mm, respectively) and hence, to estimate the cluster hydrostatic mass more precisely than previous SZ observations. However, there are certain systematic effects which can only be accounted for with the use of simulations. For this purpose, we employ THE THREE HUNDRED simulations which have been modelled with a range of physics modules to simulate galaxy clusters. The so-called twin samples are constructed by picking synthetic clusters of galaxies with properties close to the observational targets of the LPSZ. In particular, we use the Compton parameter maps and projected total mass maps of these twin samples along 29 different lines of sight. We investigate the scatter that projection induces on the total masses. Eventually, we consider the statistical values along different lines of sight to construct a kind of 3D scaling law between the integrated Compton parameter, total mass, and overdensity of the galaxy clusters to determine the overdensity that is least impacted by the projection effect. △ Less

Submitted 4 April, 2024; originally announced April 2024.

Comments: 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

An improved Halo Occupation Distribution prescription from UNITsim H_alpha Emission Line Galaxies: conformity and modified radial profile

Authors: Guillermo Reyes-Peraza, Santiago Avila, Violeta Gonzalez-Perez, Daniel Lopez-Cano, Alexander Knebe, Sujatha Ramakrishnan, Gustavo Yepes

Abstract: 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 th… ▽ More 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. △ Less

Submitted 13 March, 2024; v1 submitted 20 December, 2023; originally announced December 2023.

Comments: 17 pages, 12 figures and 5 tables

PNG-UNITsims: Halo clustering response to primordial non-Gaussianities as a function of mass

Authors: Adrian G. Adame, Santiago Avila, Violeta Gonzalez-Perez, Gustavo Yepes, Marcos Pellejero, Mike S. Wang, Chia-Hsun Chuang, Yu Feng, Juan Garcia-Bellido, Alexander Knebe

Abstract: 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.… ▽ More 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σ$. Combining all the mass bins, we find $p$ consistent with a value of $0.955\pm0.013$, which is $3σ$ 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_φ$ 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. △ Less

Submitted 30 May, 2024; v1 submitted 19 December, 2023; originally announced December 2023.

Comments: 19 pages, 11 figures. Accepted for publication in A&A

The impact of an evolving stellar initial mass function on early galaxies and reionisation

Authors: Elie Rasmussen Cueto, Anne Hutter, Pratika Dayal, Stefan Gottlöber, Kasper E. Heintz, Charlotte Mason, Maxime Trebitsch, Gustavo Yepes

Abstract: 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 I… ▽ More 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. △ Less

Submitted 25 March, 2024; v1 submitted 19 December, 2023; originally announced December 2023.

Comments: 19 pages, 13 figures, Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 686, A138 (2024)

The hydrostatic-to-lensing mass bias from resolved X-ray and optical-IR data

Authors: M. Muñoz-Echeverría, J. F. Macías-Pérez, G. W. Pratt, E. Pointecouteau, I. Bartalucci, M. De Petris, A. Ferragamo, C. Hanser, F. Kéruzoré, F. Mayet, A. Moyer-Anin, A. Paliwal, L. Perotto, G. Yepes

Abstract: An accurate reconstruction of galaxy cluster masses is key to use this population of objects as a cosmological probe. In this work we present a study on the hydrostatic-to-lensing mass scaling relation for a sample of 53 clusters whose masses were reconstructed homogeneously in a redshift range between $z= 0.05$ and $1.07$. The $M_{500}$ mass for each cluster was indeed inferred from the mass prof… ▽ More An accurate reconstruction of galaxy cluster masses is key to use this population of objects as a cosmological probe. In this work we present a study on the hydrostatic-to-lensing mass scaling relation for a sample of 53 clusters whose masses were reconstructed homogeneously in a redshift range between $z= 0.05$ and $1.07$. The $M_{500}$ mass for each cluster was indeed inferred from the mass profiles extracted from the X-ray and lensing data, without using a priori observable-mass scaling relations. We assessed the systematic dispersion of the masses estimated with our reference analyses with respect to other published mass estimates. Accounting for this systematic scatter does not change our main results, but enables the propagation of the uncertainties related to the mass reconstruction method or used dataset. Our analysis gives a hydrostatic-to-lensing mass bias of $(1-b) =0.739^{+0.075}_{-0.070}$ and no evidence of evolution with redshift. These results are robust against possible subsample differences. △ Less

Submitted 2 December, 2023; originally announced December 2023.

Journal ref: Astronomy & Astrophysics, Volume 682, id.A147, 27 pp., February 2024

The Three Hundred Project: Mapping The Matter Distribution in Galaxy Clusters Via Deep Learning from Multiview Simulated Observations

Authors: Daniel de Andres, Weiguang Cui, Gustavo Yepes, Marco De Petris, Antonio Ferragamo, Federico De Luca, Gianmarco Aversano, Douglas Rennehan

Abstract: 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 inf… ▽ More 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. △ Less

Submitted 16 January, 2024; v1 submitted 4 November, 2023; originally announced November 2023.

Comments: 15 pages, 13 figures, published in MNRAS

Improving and extending non-Poissonian distributions for satellite galaxies sampling in HOD: applications to eBOSS ELGs

Authors: Bernhard Vos-Ginés, Santiago Avila, Violeta Gonzalez-Perez, Gustavo Yepes

Abstract: 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 Dis… ▽ More 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. △ Less

Submitted 27 October, 2023; originally announced October 2023.

Comments: Submitted to MNRAS, 19 pages

The three hundred project: thermodynamical properties, shocks and gas dynamics in simulated galaxy cluster filaments and their surroundings

Authors: Agustín M. Rost, Sebastián E. Nuza, Federico Stasyszyn, Ulrike Kuchner, Matthias Hoeft, Charlotte Welker, Frazer Pearce, Meghan Gray, Alexander Knebe, Weiguang Cui, Gustavo Yepes

Abstract: 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… ▽ More 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. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: 16 pages, 13 figures. Accepted for publication in MNRAS

Galaxy clusters morphology with Zernike polynomials: the first application on $\textit{Planck}$ Compton parameter maps

Authors: Valentina Capalbo, Marco De Petris, Weiguang Cui, Antonio Ferragamo, Florian Ruppin, Gustavo Yepes

Abstract: 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 mo… ▽ More 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. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Comments: 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

Impact of filaments on galaxy cluster properties in The Three Hundred simulation

Authors: Sara Santoni, Marco De Petris, Antonio Ferragamo, Gustavo Yepes, Weiguang Cui

Abstract: 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… ▽ More 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 $χ$. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Comments: 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

The NIKA2 Sunyaev-Zeldovich Large Program: Sample and upcoming product public release

Authors: L. Perotto, R. Adam, P. Ade, H. Ajeddig, P. André, E. Artis, H. Aussel, R. Barrena, I. Bartalucci, A. Beelen, A. Benoît, S. Berta, L. Bing, O. Bourrion, M. Calvo, A. Catalano, M. De Petris, F. -X. Désert, S. Doyle, E. F. C. Driessen, G. Ejlali, A. Ferragamo, A. Gomez, J. Goupy, C. Hanser , et al. (30 additional authors not shown)

Abstract: The NIKA2 camera operating at the IRAM 30 m telescope excels in high-angular resolution mapping of the thermal Sunyaev-Zeldovich effect towards galaxy clusters at intermediate and high-redshift. As part of the NIKA2 guaranteed time, the SZ Large Program (LPSZ) aims at tSZ-mapping a representative sample of SZ-selected galaxy clusters in the catalogues of the Planck satellite and of the Atacama Cos… ▽ More The NIKA2 camera operating at the IRAM 30 m telescope excels in high-angular resolution mapping of the thermal Sunyaev-Zeldovich effect towards galaxy clusters at intermediate and high-redshift. As part of the NIKA2 guaranteed time, the SZ Large Program (LPSZ) aims at tSZ-mapping a representative sample of SZ-selected galaxy clusters in the catalogues of the Planck satellite and of the Atacama Cosmology Telescope, and also observed in X-ray with XMM Newton or Chandra. Having completed observations in January 2023, we present tSZ maps of 38 clusters spanning the targeted mass ($3 < M_{500}/10^{14} M_{\odot} < 10$) and redshift ($0.5 < z < 0.9$) ranges. The first in depth studies of individual clusters highlight the potential of combining tSZ and X-ray observations at similar angular resolution for accurate mass measurements. These were milestones for the development of a standard data analysis pipeline to go from NIKA2 raw data to the thermodynamic properties of galaxy clusters for the upcoming LPSZ data release. Final products will include unprecedented measurements of the mean pressure profile and mass observable scaling relation using a distinctive SZ-selected sample, which will be key for ultimately improving the accuracy of cluster based cosmology. △ Less

Submitted 6 October, 2023; originally announced October 2023.

Comments: 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

Journal ref: EPJ Web of Conferences 293 (2024) 00040

The Three Hundred: $M_{sub}-V_{circ}$ relation

Authors: Atulit Srivastava, Weiguang Cui, Massimo Meneghetti, Romeel Dave, Alexander Knebe, Antonio Ragagnin, Carlo Giocoli, Francesco Calura, Giulia Despali, Lauro Moscardini, Gustavo Yepes

Abstract: 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 $Λ$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 o… ▽ More 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 $Λ$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 $Λ$CDM cosmology. △ Less

Submitted 12 September, 2023; originally announced September 2023.

Comments: see the companion paper, Meneghetti et al. (2023), on today's arxiv list

A persistent excess of galaxy-galaxy strong lensing observed in galaxy clusters

Authors: Massimo Meneghetti, Weiguang Cui, Elena Rasia, Gustavo Yepes, Ana Acebron, Giuseppe Angora, Pietro Bergamini, Stefano Borgani, Francesco Calura, Giulia Despali, Carlo Giocoli, Giovanni Granata, Claudio Grillo, Alexander Knebe, Andrea Macciò, Amata Mercurio, Lauro Moscardini, Priyamvada Natarajan, Antonio Ragagnin, Piero Rosati, Eros Vanzella

Abstract: Previous studies have revealed that the estimated probability of galaxy-galaxy strong lensing in observed galaxy clusters exceeds the expectations from the $Λ$ 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 ado… ▽ More Previous studies have revealed that the estimated probability of galaxy-galaxy strong lensing in observed galaxy clusters exceeds the expectations from the $Λ$ 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. △ Less

Submitted 13 September, 2023; v1 submitted 11 September, 2023; originally announced September 2023.

Comments: 6 pages, 2 figures, accepted for publication on A&A Letters. See companion paper Srivastava et al. (2023)

Journal ref: A&A 678, L2 (2023)

The Three Hundred : contrasting clusters galaxy density in hydrodynamical and dark matter simulations

Authors: A. Jiménez Muñoz, J. F. Macías-Pérez, G. Yepes, M. De Petris, A. Ferragamo, W. Cui, J. S. Gómez

Abstract: Cluster number counts will be a key cosmological probe in the next decade thanks to the Euclid satellite mission. For this purpose, cluster detection algorithm performance, which are sensitive to the spatial distribution of the cluster galaxy members and their luminosity function, need to be accurately characterized. Using The Three Hundred hydrodynamical and dark matter only simulations we study… ▽ More Cluster number counts will be a key cosmological probe in the next decade thanks to the Euclid satellite mission. For this purpose, cluster detection algorithm performance, which are sensitive to the spatial distribution of the cluster galaxy members and their luminosity function, need to be accurately characterized. Using The Three Hundred hydrodynamical and dark matter only simulations we study a complete sample of massive clusters beyond 7 (5) $\times$ 10$^{14}$ M$_{\odot}$ at redshift 0 (1) on a $(1.48 \ \mathrm{Gpc})^3$ volume. We find that the mass resolution of the current hydrodynamical simulations (1.5 $\times$ 10$^9$ M$_{\odot}$) is not enough to characterize the luminosity function of the sample in the perspective of Euclid data. Nevertheless, these simulations are still useful to characterize the spatial distribution of the cluster substructures assuming a common relative mass threshold for the different flavours and resolutions. By comparing with the dark matter only version of these simulations, we demonstrate that baryonic physics preserves significantly low mass subhalos (galaxies) as have also been observed in previous studies with less statistics. Furthermore, by comparing the hydro simulations with higher resolution dark matter only simulations of the same objects and taking the same limit in subhalo mass we find significantly more cuspy galaxy density profiles towards the center of the clusters, where the low mass substructures would tend to concentrate. We conclude that using dark matter only simulation may lead to some biases on the spatial distribution and density of galaxy cluster members. Based on the preliminary analysis of few high resolution hydro simulations we conclude that a mass resolution of 1.8 $\times$ 10$^8$ h$^{-1}$ M$_{\odot}$ will be needed for The Three Hundred simulations to approach the expected magnitude limits for the Euclid survey. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Journal ref: A&A 686, A257 (2024)

Environmental effects on associations of dwarf galaxies

Authors: C. Yamila Yaryura, Mario G. Abadi, Stefan Gottlöber, Noam I. Libeskind, Sofía A. Cora, Andrés N. Ruiz, Cristian A. Vega-Martínez, Gustavo Yepes

Abstract: 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… ▽ More 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 $Λ$ 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. △ Less

Submitted 26 July, 2023; originally announced July 2023.

Comments: 13 pages, 9 figures. Accepted for publication in MNRAS

On the general nature of 21cm-Lyman-$α$ emitters cross-correlations during reionisation

Authors: Anne Hutter, Caroline Heneka, Pratika Dayal, Stefan Gottlöber, Andrei Mesinger, Maxime Trebitsch, Gustavo Yepes

Abstract: 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-$α$ (Ly$α$) radiation emitting galaxies (Ly$α$ 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… ▽ More 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-$α$ (Ly$α$) radiation emitting galaxies (Ly$α$ 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χ_\mathrm{HI}\rangle$) and spin-temperature weighted overdensity in neutral regions ($\langle1+δ\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. △ Less

Submitted 8 August, 2023; v1 submitted 5 June, 2023; originally announced June 2023.

Comments: 13 pages, 5 figures; accepted for publication in MNRAS

Galaxy cluster mass bias from projected mass maps: The Three Hundred-NIKA2 LPSZ twin samples

Authors: M. Muñoz-Echeverría, J. F. Macías-Pérez, E. Artis, W. Cui, D. de Andres, F. De Luca, M. De Petris, A. Ferragamo, C. Giocoli, C. Hanser, F. Mayet, M. Meneghetti, A. Moyer-Anin, A. Paliwal, L. Perotto, E. Rasia, G. Yepes

Abstract: The determination of the mass of galaxy clusters from observations is subject to systematic uncertainties. Beyond the errors due to instrumental and observational systematic effects, in this work we investigate the bias introduced by modelling assumptions. In particular, we consider the reconstruction of the mass of galaxy clusters from convergence maps employing spherical mass density models. We… ▽ More The determination of the mass of galaxy clusters from observations is subject to systematic uncertainties. Beyond the errors due to instrumental and observational systematic effects, in this work we investigate the bias introduced by modelling assumptions. In particular, we consider the reconstruction of the mass of galaxy clusters from convergence maps employing spherical mass density models. We made use of The Three Hundred simulations, selecting clusters in the same redshift and mass range as the NIKA2 Sunyaev-Zel'dovich Large Programme sample: $3 \leq M_{500}/ 10^{14} \mathrm{M}_{\odot} \leq 10$ and $0.5 \leq z \leq 0.9$. We studied different modelling and intrinsic uncertainties that should be accounted for when using the single cluster mass estimates for scaling relations. We confirm that the orientation of clusters and the radial ranges considered for the fit have an important impact on the mass bias. The effect of the projection adds uncertainties to the order of $10\%$ to $16\%$ to the mass estimates. We also find that the scatter from cluster to cluster in the mass bias when using spherical mass models is less than $9\%$ of the true mass of the clusters. △ Less

Submitted 2 December, 2023; v1 submitted 24 May, 2023; originally announced May 2023.

Journal ref: Astronomy & Astrophysics, Volume 682, id.A124, 21 pp., February 2024

Galaxy pairs in The Three Hundred simulations II: studying bound ones and identifying them via machine learning

Authors: Ana Contreras-Santos, Alexander Knebe, Weiguang Cui, Roan Haggar, Frazer Pearce, Meghan Gray, Marco De Petris, Gustavo Yepes

Abstract: 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… ▽ More 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. △ Less

Submitted 18 April, 2023; originally announced April 2023.

Comments: 20 pages, 8 figures. Published in MNRAS

Cosmic Variance and the Inhomogeneous UV Luminosity Function of Galaxies During Reionization

Authors: Taha Dawoodbhoy, Paul R. Shapiro, Pierre Ocvirk, Joseph S. W. Lewis, Dominique Aubert, Jenny G. Sorce, Kyungjin Ahn, Ilian T. Iliev, Hyunbae Park, Romain Teyssier, Gustavo Yepes

Abstract: 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- a… ▽ More 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. △ Less

Submitted 8 August, 2023; v1 submitted 16 February, 2023; originally announced February 2023.

Comments: 16 pages, 16 figures, accepted by MNRAS 07/20/23, comments welcome

Euclid preparation. XXXII. Evaluating the weak lensing cluster mass biases using the Three Hundred Project hydrodynamical simulations

Authors: Euclid Collaboration, C. Giocoli, M. Meneghetti, E. Rasia, S. Borgani, G. Despali, G. F. Lesci, F. Marulli, L. Moscardini, M. Sereno, W. Cui, A. Knebe, G. Yepes, T. Castro, P. -S. Corasaniti, S. Pires, G. Castignani, L. Ingoglia, T. Schrabback, G. W. Pratt, A. M. C. Le Brun, N. Aghanim, L. Amendola, N. Auricchio, M. Baldi , et al. (191 additional authors not shown)

Abstract: The photometric catalogue of galaxy clusters extracted from ESA Euclid data is expected to be very competitive for cosmological studies. Using state-of-the-art hydrodynamical simulations, we present systematic analyses simulating the expected weak lensing profiles from clusters in a variety of dynamic states and at wide range of redshifts. In order to derive cluster masses, we use a model consiste… ▽ More The photometric catalogue of galaxy clusters extracted from ESA Euclid data is expected to be very competitive for cosmological studies. Using state-of-the-art hydrodynamical simulations, we present systematic analyses simulating the expected weak lensing profiles from clusters in a variety of dynamic states and at wide range of redshifts. In order to derive cluster masses, we use a model consistent with the implementation within the Euclid Consortium of the dedicated processing function and find that, when jointly modelling mass and the concentration parameter of the Navarro-Frenk-White halo profile, the weak lensing masses tend to be, on average, biased low by 5-10% with respect to the true mass, up to z=0.5. Using a fixed value for the concentration $c_{200} = 3$, the mass bias is diminished below 5%, up to z=0.7, along with its relative uncertainty. Simulating the weak lensing signal by projecting along the directions of the axes of the moment of inertia tensor ellipsoid, we find that orientation matters: when clusters are oriented along the major axis, the lensing signal is boosted, and the recovered weak lensing mass is correspondingly overestimated. Typically, the weak lensing mass bias of individual clusters is modulated by the weak lensing signal-to-noise ratio, related to the redshift evolution of the number of galaxies used for weak lensing measurements: the negative mass bias tends to be larger toward higher redshifts. However, when we use a fixed value of the concentration parameter, the redshift evolution trend is reduced. These results provide a solid basis for the weak-lensing mass calibration required by the cosmological application of future cluster surveys from Euclid and Rubin. △ Less

Submitted 18 October, 2023; v1 submitted 1 February, 2023; originally announced February 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 681, A67 (2024)

A study of the hydrostatic mass bias dependence and evolution within The Three Hundred clusters

Authors: Giulia Gianfagna, Elena Rasia, Weiguang Cui, Marco De Petris, Gustavo Yepes, Ana Contreras-Santos, Alexander Knebe

Abstract: 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… ▽ More 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. △ Less

Submitted 16 November, 2022; v1 submitted 15 November, 2022; originally announced November 2022.

Comments: 11 pages, 8 figures. Accepted for publication in MNRAS

The Three Hundred project: connection between star formation quenching and dynamical evolution in and around simulated galaxy clusters

Authors: Tomás Hough, Sofía A. Cora, Roan Haggar, Cristian Vega-Martínez, Ulrike Kuchner, Frazer Pearce, Meghan Gray, Alexander Knebe, Gustavo Yepes

Abstract: 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 o… ▽ More 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. △ Less

Submitted 8 November, 2022; originally announced November 2022.

Comments: 23 pages, 13 figures + Supplementary material. Accepted for publication in MNRAS

Astraeus VIII: A new framework for Lyman-$α$ emitters applied to different reionisation scenarios

Authors: Anne Hutter, Maxime Trebitsch, Pratika Dayal, Stefan Gottlöber, Gustavo Yepes, Laurent Legrand

Abstract: We use the {\sc astraeus} framework to investigate how the visibility and spatial distribution of Lyman-$α$ (Ly$α$) 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… ▽ More We use the {\sc astraeus} framework to investigate how the visibility and spatial distribution of Lyman-$α$ (Ly$α$) 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$α$ luminosities of galaxies for three different analytic Ly$α$ line profiles and associated Ly$α$ escape fraction ($f_\mathrm{esc}^\mathrm{Lyα}$) models: importantly, we introduce two novel analytic Ly$α$ line profile models that describe the surrounding interstellar medium (ISM) as dusty gas clumps. They are based on parameterising results from radiative transfer simulations, with one of them relating $f_\mathrm{esc}^\mathrm{Lyα}$ 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$α$ 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$α$ luminosity functions reflect the Ly$α$ emission from more massive galaxies, there is a degeneracy between the $f_\mathrm{esc}$-dependent intrinsic Ly$α$ luminosity and the Ly$α$ attenuation by dust in the ISM if $f_\mathrm{esc}$ does not exceed $\sim50\%$. △ Less

Submitted 8 August, 2023; v1 submitted 29 September, 2022; originally announced September 2022.

Comments: 25 pages, 9 figures; accepted for publication in MNRAS

The Three Hundred project: Galaxy groups do not survive cluster infall

Authors: Roan Haggar, Ulrike Kuchner, Meghan E. Gray, Frazer R. Pearce, Alexander Knebe, Gustavo Yepes, Weiguang Cui

Abstract: 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… ▽ More 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. △ Less

Submitted 13 October, 2022; v1 submitted 27 September, 2022; originally announced September 2022.

Comments: 21 pages, 11 figures, published in MNRAS. v2: fixed inconsistent notation

A Deep Learning Approach to Infer Galaxy Cluster Masses from Planck Compton$-y$ parameter maps

Authors: Daniel de Andres, Weiguang Cui, Florian Ruppin, Marco De Petris, Gustavo Yepes, Giulia Gianfagna, Ichraf Lahouli, Gianmarco Aversano, Romain Dupuis, Mahmoud Jarraya, Jesús Vega-Ferrero

Abstract: 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… ▽ More 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. △ Less

Submitted 18 October, 2022; v1 submitted 21 September, 2022; originally announced September 2022.

Comments: 17 pages, 3 Figures and Supplementary material (+11 figures). Published in Nature Astronomy

Multi-probe analysis of the galaxy cluster CL J1226.9+3332: Hydrostatic mass and hydrostatic-to-lensing bias

Authors: M. Muñoz-Echeverría, J. F. Macías-Pérez, G. W. Pratt, R. Adam, P. Ade, H. Ajeddig, P. André, M. Arnaud, E. Artis, H. Aussel, I. Bartalucci, A. Beelen, A. Benoît, S. Berta, L. Bing, O. Bourrion, M. Calvo, A. Catalano, M. De Petris, F. -X. Désert, S. Doyle, E. F. C. Driessen, A. Ferragamo, A. Gomez, J. Goupy , et al. (28 additional authors not shown)

Abstract: The precise estimation of the mass of galaxy clusters is a major issue for cosmology. Large galaxy cluster surveys rely on scaling laws that relate cluster observables to their masses. From the high resolution observations of ~ 45 galaxy clusters with NIKA2 and XMM-Newton instruments, the NIKA2 SZ Large Program should provide an accurate scaling relation between the thermal Sunyaev-Zel'dovich effe… ▽ More The precise estimation of the mass of galaxy clusters is a major issue for cosmology. Large galaxy cluster surveys rely on scaling laws that relate cluster observables to their masses. From the high resolution observations of ~ 45 galaxy clusters with NIKA2 and XMM-Newton instruments, the NIKA2 SZ Large Program should provide an accurate scaling relation between the thermal Sunyaev-Zel'dovich effect and the hydrostatic mass. In this paper, we present an exhaustive analysis of the hydrostatic mass of the well known galaxy cluster CL J1226.9+3332, the highest-redshift cluster in the NIKA2 SZ Large Program at z = 0.89. We combine the NIKA2 observations with thermal Sunyaev-Zel'dovich data from NIKA, Bolocam and MUSTANG instruments and XMM-Newton X-ray observations and test the impact of the systematic effects on the mass reconstruction. We conclude that slight differences in the shape of the mass profile can be crucial when defining the integrated mass at R500, which demonstrates the importance of the modeling in the mass determination. We prove the robustness of our hydrostatic mass estimates by showing the agreement with all the results found in the literature. Another key information for cosmology is the bias of the masses estimated assuming hydrostatic equilibrium hypothesis. Based on the lensing convergence maps from the Cluster Lensing And Supernova survey with Hubble (CLASH) data, we obtain the lensing mass estimate for CL J1226.9+3332. From this we are able to measure the hydrostatic-to-lensing mass bias for this cluster, that spans from 1 - bHSE/lens ~ 0.7 to 1, presenting the impact of data-sets and mass reconstruction models on the bias. △ Less

Submitted 15 September, 2022; originally announced September 2022.

Journal ref: A&A 671, A28 (2023)

Hermeian haloes in cosmological volumes

Authors: Anastasiia Osipova, Sergey Pilipenko, Stefan Gottlöber, Noam I. Libeskind, Oliver Newton, Jenny G. Sorce, Gustavo Yepes

Abstract: 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 pe… ▽ More 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. △ Less

Submitted 15 September, 2022; originally announced September 2022.

Comments: 13 pages, 13 figures, submitted to MNRAS

Journal ref: Physics of the Dark Universe 42 (2023) 101328

Galaxy pairs in The Three Hundred simulations: a study on the performance of observational pair-finding techniques

Authors: Ana Contreras-Santos, Alexander Knebe, Weiguang Cui, Roan Haggar, Frazer Pearce, Meghan Gray, Marco De Petris, Gustavo Yepes

Abstract: 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 a… ▽ More 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. △ Less

Submitted 31 August, 2022; v1 submitted 27 July, 2022; originally announced July 2022.

Comments: 15 pages, 10 figures. Published in MNRAS; typos corrected

Reionization time of the Local Group and Local-Group-like halo pairs

Authors: Jenny G. Sorce, Pierre Ocvirk, Dominique Aubert, Stefan Gottloeber, Paul R. Shapiro, Taha Dawoodbhoy, Gustavo Yepes, Kyungjin Ahn, Ilian T. Iliev, Joseph S. W. Lewis

Abstract: 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… ▽ More 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. △ Less

Submitted 26 July, 2022; originally announced July 2022.

Comments: Accepted for publication in MNRAS, 11 pages, 5 figures

The Three Hundred project: A Machine Learning method to infer clusters of galaxies mass radial profiles from mock Sunyaev-Zel'dovich maps

Authors: A. Ferragamo, D. de Andres, A. Sbriglio, W. Cui, M. De Petris, G. Yepes, R. Dupuis, M. Jarraya, I. Lahouli, F. De Luca, G. Gianfagna, E. Rasia

Abstract: We develop a machine learning algorithm to infer the 3D cumulative radial profiles of total and gas mass in galaxy clusters from thermal Sunyaev-Zel'dovich effect maps. We generate around 73,000 mock images along various lines of sight using 2,522 simulated clusters from the \thethreehundred{} project at redshift $z< 0.12$ and train a model that combines an autoencoder and a random forest. Without… ▽ More We develop a machine learning algorithm to infer the 3D cumulative radial profiles of total and gas mass in galaxy clusters from thermal Sunyaev-Zel'dovich effect maps. We generate around 73,000 mock images along various lines of sight using 2,522 simulated clusters from the \thethreehundred{} project at redshift $z< 0.12$ and train a model that combines an autoencoder and a random forest. Without making any prior assumptions about the hydrostatic equilibrium of the clusters, the model is capable of reconstructing the total mass profile as well as the gas mass profile, which is responsible for the SZ effect. We show that the recovered profiles are unbiased with a scatter of about $10\%$, slightly increasing towards the core and the outskirts of the cluster. We selected clusters in the mass range of $10^{13.5} \leq M_{200} /(\hMsun) \leq 10^{15.5}$, spanning different dynamical states, from relaxed to disturbed halos. We verify that both the accuracy and precision of this method show a slight dependence on the dynamical state, but not on the cluster mass. To further verify the consistency of our model, we fit the inferred total mass profiles with an NFW model and contrast the concentration values with those of the true profiles. We note that the inferred profiles are unbiased for higher concentration values, reproducing a trustworthy mass-concentration relation. The comparison with a widely used mass estimation technique, such as hydrostatic equilibrium, demonstrates that our method recovers the total mass that is not biased by non-thermal motions of the gas. △ Less

Submitted 1 February, 2023; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: Accepted for publication on Monthly Notices of the Royal Astronomical Society, 8 pages, 10 figures

Astraeus VII: The environmental-dependent assembly of galaxies in the Epoch of Reionization

Authors: Laurent Legrand, Pratika Dayal, Anne Hutter, Stefan Gottlöber, Gustavo Yepes, Maxime Trebitsch

Abstract: 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 (… ▽ More 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+δ) < 1.3$ averaged over $(2~{\rm cMpc})^3$). We find that : (i) there exists a mass- and redshift- dependent "characteristic" environment (${\rm log} (1+δ_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+δ) > 0.5$); dominated by their major branch, this yields star formation histories biased towards older ages with a slower redshift evolution. △ Less

Submitted 15 December, 2022; v1 submitted 14 July, 2022; originally announced July 2022.

Comments: 17 pages, 15 figures, submitted to MNRAS, comments welcome !

A multi-simulation study of relativistic SZ temperature scalings in galaxy clusters and groups

Authors: Elizabeth Lee, Dhayaa Anbajagane, Priyanka Singh, Jens Chluba, Daisuke Nagai, Scott T. Kay, Weiguang Cui, Klaus Dolag, Gustavo Yepes

Abstract: 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 mas… ▽ More 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. △ Less

Submitted 26 September, 2022; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: Accepted by MNRAS

A stochastic model to reproduce the star-formation history of individual galaxies in hydrodynamic simulations

Authors: Yang Wang, Nicola R. Napolitano, Weiguang Cui, Xiao-Dong Li, Alexander Knebe, Chris Power, Frazer Pearce, Lin Tang, Gustavo Yepes, Xi Kang

Abstract: 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 galax… ▽ More 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. △ Less

Submitted 27 July, 2022; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: 22 pages, 18 figures, 2 tables

Massive merging cluster PSZ2G091 as seen by the NIKA2 camera

Authors: E. Artis, R. Adam, P. Ade, H. Ajeddig, P. André, M. Arnaud, H. Aussel, I. Bartalucci, A. Beelen, A. Benoît, S. Berta, L. Bing, O. Bourrion, M. Calvo, A. Catalano, M. De Petris, F. -X. Désert, S. Doyle, E. F. C. Driessen, A. Ferragamo, A. Gomez, J. Goupy, C. Hanser, F. Kéruzoré, C. Kramer , et al. (27 additional authors not shown)

Abstract: PSZ2 G091.83+26.11 is a galaxy cluster with M500 = 7.43 x 10^14 Msun at z = 0.822 1. This object exhibits a complex morphology with a clear bimodality observed in X-rays. However, it was detected and analysed in the Planck sample as a single, spherical cluster following a universal profile 2. This model can lead to miscalculations of thermodynamical quantities, like the pressure profile. As future… ▽ More PSZ2 G091.83+26.11 is a galaxy cluster with M500 = 7.43 x 10^14 Msun at z = 0.822 1. This object exhibits a complex morphology with a clear bimodality observed in X-rays. However, it was detected and analysed in the Planck sample as a single, spherical cluster following a universal profile 2. This model can lead to miscalculations of thermodynamical quantities, like the pressure profile. As future multiwavelength cluster experiments will detect more and more objects at high redshifts, it is crucial to quantify this systematic effect. In this work, we use high-resolution observations of the NIKA2 camera3,4,5,6 to integrate the morphological characteristics of the cluster in our modelling. This is achieved by fitting a two-halo model to the SZ image and then by reconstruction of the resulting projected pressure profile. We then compare these results with the spherical assumption. △ Less

Submitted 29 April, 2022; originally announced April 2022.

Comments: 4 pages, 3 figures Contribution to the 2022 Cosmology session of the 56th Rencontres de Moriond. arXiv admin note: substantial text overlap with arXiv:2111.05082

Machine Learning methods to estimate observational properties of galaxy clusters in large volume cosmological N-body simulations

Authors: Daniel de Andres, Gustavo Yepes, Federico Sembolini, Gonzalo Martínez-Muñoz, Weiguang Cui, Francisco Robledo, Chia-Hsun Chuang, Elena Rasia

Abstract: 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… ▽ More 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. △ Less

Submitted 10 November, 2022; v1 submitted 22 April, 2022; originally announced April 2022.

Comments: 20 pages, 10 figures, published in MNRAS

A Machine Learning Approach to Correct for Mass Resolution Effects in Simulated Halo Clustering Statistics

Authors: Daniel Forero-Sánchez, Chia-Hsun Chuang, Sergio Rodríguez-Torres, Gustavo Yepes, Stefan Gottlöber, Cheng Zhao

Abstract: 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 galaxie… ▽ More 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. △ Less

Submitted 29 April, 2022; v1 submitted 23 March, 2022; originally announced March 2022.

Comments: 15 pages, 13 figures, Accepted for publication in MNRAS

\textsc{The Three Hundred} project: The \textsc{Gizmo-Simba} run

Authors: Weiguang Cui, Romeel Dave, Alexander Knebe, Elena Rasia, Meghan Gray, Frazer Pearce, Chris Power, Gustavo Yepes, Dhayaa Anbajagane, Daniel Ceverino, Ana Contreras-Santos, Daniel de Andres, Marco De Petris, Stefano Ettori, Roan Haggar, Qingyang Li, Yang Wang, Xiaohu Yang, Stefano Borgani, Klaus Dolag, Ying Zu, Ulrike Kuchner, Rodrigo Cañas, Antonio Ferragamo, Giulia Gianfagna

Abstract: 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 ga… ▽ More 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. △ Less

Submitted 31 May, 2022; v1 submitted 28 February, 2022; originally announced February 2022.

Comments: 20 pages, 12 figures, matched to published version in MNRAS

The short ionizing photon mean free path at z=6 in Cosmic Dawn III, a new fully-coupled radiation-hydrodynamical simulation of the Epoch of Reionization

Authors: Joseph S. W. Lewis, Pierre Ocvirk, Jenny G. Sorce, Yohan Dubois, Dominique Aubert, Luke Conaboy, Paul R. Shapiro, Taha Dawoodbhoy, Romain Teyssier, Gustavo Yepes, Stefan Gottlöber, Yann Rasera, Kyungjin Ahn, Ilian T. Iliev, Hyunbae Park, Émilie Thélie

Abstract: 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$) a… ▽ More 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$). △ Less

Submitted 22 August, 2022; v1 submitted 11 February, 2022; originally announced February 2022.

Comments: Submitted to MNRAS Accepted 2022 August 08. Received 2022 August 05; in original form 2022 February 11

Astraeus VI: Hierarchical assembly of AGN and their large-scale effect during the Epoch of Reionization

Authors: Maxime Trebitsch, Anne Hutter, Pratika Dayal, Stefan Gottlöber, Laurent Legrand, Gustavo Yepes

Abstract: 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… ▽ More 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. △ Less

Submitted 26 July, 2022; v1 submitted 4 February, 2022; originally announced February 2022.

Comments: 15+2 pages, 14+1 figures, 1 table. Accepted in MNRAS

The Three Hundred project: galaxy cluster mergers and their impact on the stellar component of brightest cluster galaxies

Authors: Ana Contreras-Santos, Alexander Knebe, Frazer Pearce, Roan Haggar, Meghan Gray, Weiguang Cui, Gustavo Yepes, Marco De Petris, Federico De Luca, Chris Power, Robert Mostoghiu, Sebastián E. Nuza, Matthias Hoeft

Abstract: 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 in… ▽ More 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. △ Less

Submitted 28 January, 2022; originally announced January 2022.

Comments: 18 pages, 12 figures. Accepted for publication in MNRAS

Numerical simulations of bar formation in the Local Group

Authors: Ornela F. Marioni, Mario G. Abadi, Stefan Gottlöber, Gustavo Yepes

Abstract: 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 $Λ$CDM cosmological framework predict the formation of such stellar components, there seems to be a tension between theoretical and observational results. Simulated bars are typicall… ▽ More 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 $Λ$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 $Λ$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. △ Less

Submitted 13 January, 2022; originally announced January 2022.

Comments: 11 pages, 9 figures. Accepted for publication in MNRAS

Covariance matrices for variance-suppressed simulations

Authors: Tony Zhang, Chia-Hsun Chuang, Risa H. Wechsler, Shadab Alam, Joseph DeRose, Yu Feng, Francisco-Shu Kitaura, Marcos Pellejero-Ibanez, Sergio Rodríguez-Torres, Chun-Hao To, Gustavo Yepes, Cheng Zhao

Abstract: 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 propose… ▽ More 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. △ Less

Submitted 17 March, 2023; v1 submitted 20 December, 2021; originally announced December 2021.

Comments: 9 pages, 7 figures; in v2 we incorporate minor modifications from peer review and copyediting

Journal ref: Monthly Notices of the Royal Astronomical Society, 518, 3737 (2023)

Astraeus V: The emergence and evolution of metallicity scaling relations during the Epoch of Reionization

Authors: Graziano Ucci, Pratika Dayal, Anne Hutter, Chiaki Kobayashi, Stefan Gottloeber, Gustavo Yepes, Leslie Hunt, Laurent Legrand, Crescenzo Tortora

Abstract: 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… ▽ More 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 $η_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. △ Less

Submitted 14 September, 2022; v1 submitted 3 December, 2021; originally announced December 2021.

Comments: Accepted to MNRAS

An inventory of galaxies in cosmic filaments feeding galaxy clusters: galaxy groups, backsplash galaxies, and pristine galaxies

Authors: Ulrike Kuchner, Roan Haggar, Alfonso Aragón-Salamanca, Frazer R. Pearce, Meghan E. Gray, Agustín Rost, Weiguang Cui, Alexander Knebe, Gustavo Yepes

Abstract: 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… ▽ More 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. △ Less

Submitted 22 November, 2021; originally announced November 2021.

Comments: 13 pages, 9 figures, accepted 2021 November 21

PSZ2G091:A massive double cluster at z=0.822 observed by the NIKA2 camera

Authors: E. Artis, R. Adam, P. Ade, H. Ajeddig, P. André, M. Arnaud, H. Aussel, I. Bartalucci, A. Beelen, A. Benoît, S. Berta, L. Bing, O. Bourrion, M. Calvo, A. Catalano, M. De Petris, F. -X. Désert, S. Doyle, E. F. C. Driessen, A. Ferragamo, A. Gomez, J. Goupy, F. Kéruzoré, C. Kramer, B. Ladjelate , et al. (26 additional authors not shown)

Abstract: PSZ2 G091.83+26.11 is a massive galaxy cluster with M500 = 7.43 x 10^14 Msun at z = 0.822. This object exhibits a complex morphology with a clear bimodality observed in X-rays. However, it was detected and analysed in the Planck sample as a single, spherical cluster following a universal profile [1]. This model can lead to miscalculations of thermodynamical quantities, like the pressure profile. A… ▽ More PSZ2 G091.83+26.11 is a massive galaxy cluster with M500 = 7.43 x 10^14 Msun at z = 0.822. This object exhibits a complex morphology with a clear bimodality observed in X-rays. However, it was detected and analysed in the Planck sample as a single, spherical cluster following a universal profile [1]. This model can lead to miscalculations of thermodynamical quantities, like the pressure profile. As future multiwavelength cluster experiments will detect more and more objects at higher redshifts (where we expect the fraction of merging objects to be higher), it is crucial to quantify this systematic effect. In this work, we use high-resolution observations of PSZ2 G091.83+26.11 by the NIKA2 camera to integrate the morphological characteristics of the cluster in our modelling. This is achieved by fitting a two-halo model to the SZ image and then by reconstruction of the resulting projected pressure profile. We then compare these results with the spherical assumption. △ Less

Submitted 9 November, 2021; originally announced November 2021.

Comments: "To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences."

Mass Estimation of Planck Galaxy Clusters using Deep Learning

Authors: Daniel de Andres, Weiguang Cui, Florian Ruppin, Marco De Petris, Gustavo Yepes, Ichraf Lahouli, Gianmarco Aversano, Romain Dupuis, Mahmoud Jarraya

Abstract: 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… ▽ More 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. △ Less

Submitted 3 December, 2021; v1 submitted 2 November, 2021; originally announced November 2021.

Comments: To appear in the Proceedings of the International Conference entitled "mm Universe @NIKA2", Rome(Italy), June 2021, EPJ Web of conferences

The Three Hundred-NIKA2 Sunyaev-Zeldovich Large Program twin samples: Synthetic clusters to support real observations

Authors: A. Paliwal, E. Artis, W. Cui, M. De Petris, F. -X. Désert, A. Ferragamo, G. Gianfagna, F. Kéruzoré, J. -F. Macías-Pérez, F. Mayet, M. Muñoz-Echeverría, L. Perotto, E. Rasia, F. Ruppin, G. Yepes

Abstract: The simulation database of THE THREE HUNDRED Project has been used to pick synthetic clusters of galaxies with properties close to the observational targets of the NIKA2 camera Sunyaev-Zeldovich (SZ) Large Program. Cross-matching of cluster parameters such as mass and redshift of the cluster in the two databases has been implemented to generate the so-called twin samples for the Large Program. Thi… ▽ More The simulation database of THE THREE HUNDRED Project has been used to pick synthetic clusters of galaxies with properties close to the observational targets of the NIKA2 camera Sunyaev-Zeldovich (SZ) Large Program. Cross-matching of cluster parameters such as mass and redshift of the cluster in the two databases has been implemented to generate the so-called twin samples for the Large Program. This SZ Large Program is observing a selection of galaxy clusters at intermediate and high redshift $\left( 0.5 < z < 0.9 \right)$, covering one order of magnitude in mass. These are SZ-selected clusters from the Planck and Atacama Cosmology Telescope catalogs, wherein the selection is based on their integrated Compton parameter values, $Y_{500}$: the value of the parameter within the characteristics radius $R_{500}$. THE THREE HUNDRED hydrodynamical simulations provide us with hundreds of clusters satisfying these redshift, mass, and $Y_{500}$ requirements. In addition to the standard post-processing analysis of the simulation, mock observational maps are available mimicking X-ray, optical, gravitational lensing, radio, and SZ observations of galaxy clusters. The primary goal of employing the twin samples is to compare different cluster mass proxies from synthetic X-ray, SZ effect and optical maps (via the velocity dispersion of member galaxies and lensing $κ$-maps) of the clusters. Eventually, scaling laws between different mass proxies and the cluster mass will be cross-correlated to reduce the scatter on the inferred mass and the mass bias will be related to various physical parameters. △ Less

Submitted 1 December, 2021; v1 submitted 2 November, 2021; originally announced November 2021.

Comments: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences

The hydrostatic mass bias in The Three Hundred clusters

Authors: Giulia Gianfagna, Elena Rasia, Weiguang Cui, Marco De Petris, Gustavo Yepes

Abstract: 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… ▽ More 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. △ Less

Submitted 20 December, 2021; v1 submitted 2 November, 2021; originally announced November 2021.

Comments: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences