BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Weiguang Cui

Weiguang Cui appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-ph.HEastro-ph.IMastro-ph.SR

Catalog footprint

What is connected

49works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

49 published item(s)

preprint2026arXiv

The Advanced X-ray Imaging Satellite (AXIS) Community Science Book

The AXIS Community Science Book represents the collective effort of 592 scientists worldwide to define the transformative science enabled by the Advanced X-ray Imaging Satellite (AXIS), a next-generation X-ray mission selected by NASA's Astrophysics Probe Program for Phase A study. AXIS will advance the legacy of high-angular-resolution X-ray astronomy with ~1.5'' imaging over a wide 24' field of view and an order of magnitude greater collecting area than Chandra in the 0.3-12 keV band. Combining sharp imaging, high throughput, and rapid response capabilities, AXIS will open new windows on virtually every aspect of modern astrophysics, exploring the birth and growth of supermassive black holes, the feedback processes that shape galaxies, the life cycles of stars and exoplanet environments, and the nature of compact stellar remnants, supernova remnants, and explosive transients. This book compiles 138 community-contributed science cases developed by five Science Working Groups focused on AGN and supermassive black holes, galaxy evolution and feedback, compact objects and supernova remnants, stellar physics and exoplanets, and time-domain and multi-messenger astrophysics. Together, these studies establish the scientific foundation for next-generation X-ray exploration in the 2030s and highlight strong synergies with facilities of the 2030s, such as JWST, Roman, Rubin/LSST, SKA, ALMA, ngVLA, and next-generation gravitational-wave and neutrino networks.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Deriving accurate galaxy cluster masses using X-ray thermodynamic profiles and graph neural networks

Precise determination of galaxy cluster masses is crucial for establishing reliable mass-observable scaling relations in cluster cosmology. We employ graph neural networks (GNNs) to estimate galaxy cluster masses from radially sampled profiles of the intra-cluster medium (ICM) inferred from X-ray observations. GNNs naturally handle inputs of variable length and resolution by representing each ICM profile as a graph, enabling accurate and flexible modeling across diverse observational conditions. We trained and tested GNN model using state-of-the-art hydrodynamical simulations of galaxy clusters from The Three Hundred Project. The mass estimates using our method exhibit no systematic bias compared to the true cluster masses in the simulations. Additionally, we achieve a scatter in recovered mass versus true mass of about 6%, which is a factor of six smaller than obtained from a standard hydrostatic equilibrium approach. Our algorithm is robust to both data quality and cluster morphology and it is capable of incorporating model uncertainties alongside observational uncertainties. Finally, we apply our technique to XMM-Newton observed galaxy cluster samples and compare the GNN derived mass estimates with those obtained with $Y_{\rm SZ}$-M$_{500}$ scaling relations. Our results provide strong evidence, at 5$σ$ level, for a mass-dependent bias in SZ derived masses, with higher mass clusters exhibiting a greater degree of deviation. Furthermore, we find the median bias to be $(1-b)=0.85_{-0.14}^{+0.34}$, albeit with significant dispersion due to its mass dependence. This work takes a significant step towards establishing unbiased observable mass scaling relations by integrating X-ray, SZ and optical datasets using deep learning techniques, thereby enhancing the role of galaxy clusters in precision cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

The Physical Nature of Circumgalactic Medium Absorbers in Simba

We study the nature of the low-redshift CGM in the Simba cosmological simulations as traced by ultraviolet absorption lines around galaxies in bins of stellar mass ($M_\star>10^{10}M_\odot$) for star-forming, green valley and quenched galaxies at impact parameters $r_\perp\leq 1.25r_{200}$. We generate synthetic spectra for HI, MgII, CII, SiIII, CIV, and OVI, fit Voigt profiles to obtain line properties, and estimate the density, temperature, and metallicity of the absorbing gas. We find that CGM absorbers are most abundant around star forming galaxies with $M_\star < 10^{11}M_\odot$, while the abundance of green valley galaxies show similar behaviour to those of quenched galaxies, suggesting that the CGM "quenches" before star formation ceases. HI absorbing gas exists across a broad range of cosmic phases (condensed gas, diffuse gas, hot halo gas and Warm-Hot Intergalactic Medium), while essentially all low-ionisation metal absorption arises from condensed gas. OVI absorbers are split between hot halo gas and the WHIM. The fraction of collisionally ionised CGM absorbers is $\sim 25-55\%$ for CIV and $\sim 80-95\%$ for OVI, depending on stellar mass and impact parameter. In general, the highest column density absorption features for each ion arise from dense gas. Satellite gas, defined as that within $10r_{1/2,\star},$ contributes $\sim 3\%$ of overall HI absorption but $\sim 30\%$ of MgII absorption, with the fraction from satellites decreasing with increasing ion excitation energy.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Brightest Cluster Galaxies Trace Weak Lensing Mass Bias and Halo Triaxiality in The Three Hundred Project

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Effects of Active Galactic Nucleus Feedback on Cold Gas Depletion and Quenching of Central Galaxies

We investigate the influence of active galactic nucleus (AGN) feedback on the galaxy cold gas content and its connection to galaxy quenching in three hydrodynamical simulations of Illustris, IllustrisTNG and SIMBA. By comparing to the observed atomic and molecular neutral hydrogen measurements for central galaxies, we find that Illustris over-predicts the cold gas masses in star-forming galaxies and significantly under-predicts them for quenched galaxies. IllustrisTNG performs better in this comparison than Illustris, but quenched galaxies retain too much cold gas compared with observations. SIMBA shows good agreement with observations, by depleting the global cold gas reservoir for quenched galaxies. We find that the discrepancies in IllustrisTNG are caused by its weak kinetic AGN feedback that only redistributes the cold gas from the inner disks to the outer regions and reduces the inner cold gas densities. It agrees with observations much better when only the cold gas within the stellar disk is considered to infer the star formation rates. From dependences of cold gas reservoir on the black hole mass and Eddington ratio, we find that the cumulative energy release during the black hole growth is the dominant reason for the cold gas depletion and thus the galaxy quenching. We further measure the central stellar surface density within 1 kpc ($Σ_1$) for the high-resolution run of IllustrisTNG and find a tight correlation between $Σ_1$ and black hole mass. It suggests that the observed decreasing trend of cold gas mass with $Σ_1$ is also a reflection of the black hole growth.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

ELUCID VII: Using Constrained Hydro Simulations to Explore the Gas Component of the Cosmic Web

Using reconstructed initial conditions in the SDSS survey volume, we carry out constrained hydrodynamic simulations in three regions representing different types of the cosmic web: the Coma cluster of galaxies; the SDSS great wall; and a large low-density region at $z\sim 0.05$. These simulations, which include star formation and stellar feedback but no AGN formation and feedback, are used to investigate the properties and evolution of intergalactic and intra-cluster media. About half of the warm-hot intergalactic gas is associated with filaments in the local cosmic web. Gas in the outskirts of massive filaments and halos can be heated significantly by accretion shocks generated by mergers of filaments and halos, respectively, and there is a tight correlation between gas temperature and the strength of the local tidal field. The simulations also predict some discontinuities associated with shock fronts and contact edges, which can be tested using observations of the thermal SZ effect and X-rays. A large fraction of the sky is covered by Ly$α$ and OVI absorption systems, and most of the OVI systems and low-column density HI systems are associated with filaments in the cosmic web. The constrained simulations, which follow the formation and heating history of the observed cosmic web, provide an important avenue to interpret observational data. With full information about the origin and location of the cosmic gas to be observed, such simulations can also be used to develop observational strategies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Groups and protocluster candidates in the CLAUDS and HSC-SSP joint deep surveys

Using the extended halo-based group finder developed by Yang et al. (2021), which is able to deal with galaxies via spectroscopic and photometric redshifts simultaneously, we construct galaxy group and candidate protocluster catalogs in a wide redshift range ($0 < z < 6$) from the joint CFHT Large Area $U$-band Deep Survey (CLAUDS) and Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) deep data set. Based on a selection of 5,607,052 galaxies with $i$-band magnitude $m_{i} < 26$ and a sky coverage of $34.41\ {\rm deg}^2$, we identify a total of 2,232,134 groups, within which 402,947 groups have at least three member galaxies. We have visually checked and discussed the general properties of those richest groups at redshift $z>2.0$. By checking the galaxy number distributions within a $5-7\ h^{-1}\mathrm{Mpc}$ projected separation and a redshift difference $Δz \le 0.1$ around those richest groups at redshift $z>2$, we identified a list of 761, 343 and 43 protocluster candidates in the redshift bins $2\leq z<3$, $3\leq z<4$ and $z \geq 4$, respectively. In general, these catalogs of galaxy groups and protocluster candidates will provide useful environmental information in probing galaxy evolution along the cosmic time.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

High-ionisation oxygen absorption from the Warm-Hot Intergalactic Medium in Simba

We examine the physical conditions, environments, and statistical properties of intergalactic Ovi, Ovii and Oviii absorbers in the Simba cosmological hydrodynamic simulation suite. The goal is to understand the nature of these high ionisation absorbers, and test Simba's surprising prediction that $\sim 70\%$ of cosmic baryons at $z=0$ are in the Warm-Hot Intergalactic Medium (WHIM) driven by jet feedback from active galactic nuclei (AGN). By comparing a full-physics Simba run versus one with jets turned off, we find that jet feedback causes widespread heating that impacts the absorption morphology particularly of the higher ions. However, the distribution of the physical properties of detectable absorbers are not dramatically affected. Higher ionisation absorbers probe hotter gas as expected, but in Simba all ions arise at similar overdensities (typically $δ\sim20-30$), similar environments (predominantly filaments), and similar nearest-halo distances (typically $\sim2-3r_{200c}$). Simba matches the observed Ovi column density distribution function (CDDF) fairly well, but under-predicts the CDDF preliminarily derived from two detected intergalactic Ovii absorbers. Predicted CDDFs are very similar at $z=1$ with or without jets, but show differences by $z=0$ particularly at the high-column end. Despite some discrepancies, Simba reproduces available observations as well as or better than other comparable simulations, suggesting that Simba's widespread jet heating cannot be ruled out by these data. These results offer hope that future X-ray and ultraviolet facilities could provide significant constraints on galactic feedback models from high-ionisation IGM metal absorbers.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Three Hundred Project: dissecting the fundamental plane of galaxy clusters up to $z=1$

We perform a systematic study of the recently discovered fundamental plane of galaxy clusters (CFP) using ~250 simulated clusters from The Three Hundred project, focusing on the stability of the plane against different temperature definitions and its dependence on the dynamical relaxation state of clusters. The CFP is characterised by $T\propto M_\mathrm{s}^α\,r_\mathrm{s}^β$, defined with the gas temperature ($T$) and the characteristic halo scale radius and mass ($r_\mathrm{s}$ and $M_\mathrm{s}$) assuming a Navarro-Frenk-White halo description. We explore two definitions of weighted temperatures, namely mass-weighted and spectroscopic-like temperatures, in three radial ranges. The Three Hundred project clusters at $z=0$ lie on a thin plane whose parameters ($α,β$) and dispersion ($0.015-0.030$ dex) depend on the gas temperature definition. The CFP for mass-weighted temperatures is closer to the virial equilibrium expectation ($α=1,β=-1$) with a smaller dispersion. For gas temperatures measured within $500h^{-1}$ kpc, the resulting CFP deviates the most from the virial expectation and shifts towards the similarity solution for a secondary infall model ($α=1.5,β=-2$). Independently of the temperature definition, we find that clusters at $z=1$ and relaxed clusters form a CFP similar to the virial expectation, unlike disturbed clusters exhibiting stronger evolution. Only systems formed over the last 4 Gyr present a CFP that is closer to the similarity solution. All these findings are compatible with the CFP obtained for a CLASH subsample excluding the hottest clusters with $T_\mathrm{X}>12$ keV.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Understanding the relation between thermal Sunyaev-Zeldovich decrement and halo mass using the SIMBA and TNG simulations

The relation between the integrated thermal Sunyaev-Zeldovich (tSZ) $y$-decrement versus halo mass ($Y$--$M$) can potentially constrain galaxy formation models, if theoretical and observational systematics can be properly assessed. We investigate the $Y$--$M$ relation in the SIMBA and IllustrisTNG-100 cosmological hydrodynamic simulations, quantifying the effects of feedback, line-of-sight projection, and beam convolution. We find that SIMBA's AGN jet feedback generates strong deviations from self-similar expectations for the $Y$--$M$ relation, especially at $M_{500}<10^{13}M_{\odot}$. In SIMBA, this is driven by suppressed in-halo $y$ contributions owing to lowered halo baryon fractions. IllustrisTNG results more closely resemble SIMBA without jets. Projections of line-of-sight structures weaken these model differences slightly, but they remain significant -- particularly at group and lower halo masses. In contrast, beam smearing at $\textit{Planck}$ resolution makes the models indistinguishable, and both models appear to agree well with $\textit{Planck}$ data down to the lowest masses probed. We show that the arcminute resolution expected from forthcoming facilities would retain the differences between model predictions, and thereby provide strong constraints on AGN feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

What to expect from dynamical modelling of cluster haloes II. Investigating dynamical state indicators with Random Forest

We investigate the importances of various dynamical features in predicting the dynamical state (DS) of galaxy clusters, based on the Random Forest (RF) machine learning approach. We use a large sample of galaxy clusters from the Three Hundred Project of hydrodynamical zoomed-in simulations, and construct dynamical features from the raw data as well as from the corresponding mock maps in the optical, X-ray, and Sunyaev-Zel'dovich (SZ) channels. Instead of relying on the impurity based feature importance of the RF algorithm, we directly use the out-of-bag (OOB) scores to evaluate the importances of individual features and different feature combinations. Among all the features studied, we find the virial ratio, $η$, to be the most important single feature. The features calculated directly from the simulations and in 3-dimensions carry more information on the DS than those constructed from the mock maps. Compared with the features based on X-ray or SZ maps, features related to the centroid positions are more important. Despite the large number of investigated features, a combination of up to three features of different types can already saturate the score of the prediction. Lastly, we show that the most sensitive feature $η$ is strongly correlated with the well-known half-mass bias in dynamical modelling. Without a selection in DS, cluster halos have an asymmetric distribution in $η$, corresponding to an overall positive half-mass bias. Our work provides a quantitative reference for selecting the best features to discriminate the DS of galaxy clusters in both simulations and observations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Exploring the hydrostatic mass bias in MUSIC clusters: application to the NIKA2 mock sample

Clusters of galaxies are useful tools to constrain cosmological parameters, only if their masses can be correctly inferred from observations. In particular, X-ray and Sunyaev-Zeldovich (SZ) effect observations can be used to derive masses within the framework of the hydrostatic equilibrium. Therefore, it is crucial to have a good control of the possible mass biases that can be introduced when this hypothesis is not valid. In this work, we analyzed a set of 260 synthetic clusters from the MUSIC simulation project, at redshifts $0 \leq z \leq 0.82$. We estimate the hydrostatic mass of the MUSIC clusters from X-ray only (temperature and density) and from X-ray and SZ (density and pressure). Then, we compare them with the true 3D dynamical mass. The biases are of the order of 20%. We find that using the temperature instead of the pressure leads to a smaller bias, although the two values are compatible within 1$σ$. Non-thermal contributions to the total pressure support, arising from bulk motion and turbulence of the gas, are also computed and show that they are sufficient to account for this bias. We also present a study of the correlation between the mass bias and the dynamical state of the clusters. A clear correlation is shown between the relaxation state of the clusters and the bias factor. We applied the same analysis on a subsample of 32 objects, already selected for supporting the NIKA2 SZ Large Program.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Galaxy Velocity Bias in Cosmological Simulations: Towards Percent-level Calibration

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

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Mass Estimation of Planck Galaxy Clusters using Deep Learning

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Morphological analysis of SZ and X-ray maps of galaxy clusters with Zernike polynomials

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

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The hydrostatic mass bias in The Three Hundred clusters

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Three Hundred Project: The gas disruption of infalling objects in cluster environments

We analyse the gas content evolution of infalling haloes in cluster environments from THE THREE HUNDRED project, a collection of 324 numerically modelled galaxy clusters. The haloes in our sample were selected within $5R_{200}$ of the main cluster halo at $z=0$ and have total halo mass $M_{200}\geq10^{11} h^{-1} M_{\odot}$. We track their main progenitors and study their gas evolution since their crossing into the infall region, which we define as $1-4R_{200}$. Studying the radial trends of our populations using both the full phase space information and a line-of-sight projection, we confirm the Arthur et al. (2019) result and identify a characteristic radius around $1.7R_{200}$ in 3D and at $R_{200}$ in projection at which infalling haloes lose nearly all of the gas prior their infall. Splitting the trends by subhalo status we show that subhaloes residing in group-mass and low-mass host haloes in the infall region follow similar radial gas-loss trends as their hosts, whereas subhaloes of cluster-mass host haloes are stripped of their gas much further out. Our results show that infalling objects suffer significant gaseous disruption that correlates with time-since-infall, cluster-centric distance and host mass, and that the gaseous disruption they experience is a combination of subhalo pre-processing and object gas depletion at a radius which behaves like an accretion shock.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Three Hundred: Cluster Dynamical States and Relaxation Time Scale

Using the galaxy clusters from The Three Hundred Project, we define a new parameter: $λ_{DS}$ to describe the dynamical state of clusters, which assumes a double-Gaussian distribution in logarithm scale for our mass-complete cluster sample at $z=0$ from the dark-matter-only (DMO) run. Therefore, the threshold for distinguishing relaxed and unrelaxed clusters is naturally determined by the crossing point of the double-Gaussian fitting which has a value of $λ_{DS} = 3.424$. By applying $λ_{DS}$ with the same parameters from the DMO run to the hydro-dynamically simulated clusters (Gadget-X run and GIZMO-SIMBA run), we investigate the effect of baryons on the cluster dynamical state. We find a weak baryon-model dependence for the $λ_{DS}$ parameter. Finally, we study the evolution of $λ_{DS}$ along with clusters mass accretion history. We notice an upper limit of halo mass change $\frac{ΔM_{200}}{M_{200}} \sim 0.12$ that don't alter the cluster dynamical state, i.e. from relaxed to unrelaxed. We define relaxation period (from the most relaxed state to disturb and relaxed again) which reflects how long the dynamical state of a cluster restores its relaxation state, and propose a correlation between this relaxation period and the strength of halo mass change $\frac{ΔM_{200}}{M_{200}}$. With the proposed fitting to such correlation, we verify the relaxation period can be estimated from $\frac{ΔM_{200}}{M_{200}}$ (including multi mass change peaks) with considerably small error.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The THREEHUNDRED project: the effect of baryon processes at galaxy cluster scale

The role of baryon models played in hydrodynamic simulations is still unclear. Future surveys that use cluster statistics to precisely constrain cosmology models require a better understanding of that. With the hydro-simulated galaxy clusters from different baryon models (Gadget-MUSIC, Gadget-X and Gizmo-SIMBA) from the THREEHUNDRED project, we can look into more details of this question. We find that the galaxy cluster mass change due to different baryon models is at a few per cent level. However, the mass changes can be positive or negative, which is depending on the baryon models. Such a small mass change leaves a weak influence (slightly larger compared to the mass changes) on both the cumulative halo numbers and the differential halo mass function (HMF) above the mass completeness. Agreed to the halo mass change, the halo mass (or HMF) can be increased or decreased with respect to the dark-matter-only (DMO) run depending on the baryon models.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Velocity dispersion vs cluster mass: a new scaling law with The Three Hundred clusters

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Mapping and characterisation of cosmic filaments in galaxy cluster outskirts: strategies and forecasts for observations from simulations

Upcoming wide-field surveys are well-suited to studying the growth of galaxy clusters by tracing galaxy and gas accretion along cosmic filaments. We use hydrodynamic simulations of volumes surrounding 324 clusters from \textsc{The ThreeHundred} project to develop a framework for identifying and characterising these filamentary structures, and associating galaxies with them. We define 3-dimensional reference filament networks reaching $5R_{200}$ based on the underlying gas distribution and quantify their recovery using mock galaxy samples mimicking observations such as those of the WEAVE Wide-Field Cluster Survey. Since massive galaxies trace filaments, they are best recovered by mass-weighting galaxies or imposing a bright limit (e.g. $>L^*$) on their selection. We measure the transverse gas density profile of filaments, derive a characteristic filament radius of $\simeq0.7$--$1~h^{-1}\rm{Mpc}$, and use this to assign galaxies to filaments. For different filament extraction methods we find that at $R>R_{200}$, $\sim15$--$20%$ of galaxies with $M_*>3 \times 10^9 M_{\odot}$ are in filaments, increasing to $\sim60%$ for galaxies more massive than the Milky-Way. The fraction of galaxies in filaments is independent of cluster mass and dynamical state, and is a function of cluster-centric distance, increasing from $\sim13$% at $5R_{200}$ to $\sim21$% at $1.5R_{200}$. As a bridge to the design of observational studies, we measure the purity and completeness of different filament galaxy selection strategies. Encouragingly, the overall 3-dimensional filament networks and $\sim67$% of the galaxies associated with them are recovered from 2-dimensional galaxy positions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Protoclusters at z=5.7: A view from the MultiDark galaxies

Protoclusters, which will yield galaxy clusters at lower redshift, can provide valuable information on the formation of galaxy clusters. However, identifying progenitors of galaxy clusters in observations is not an easy task, especially at high redshift. Different priors have been used to estimate the overdense regions that are thought to mark the locations of protoclusters. In this paper, we use mimicked Ly$α$-emitting galaxies at $z=5.7$ to identify protoclusters in the MultiDark galaxies, which are populated by applying three different semi-analytic models to the 1 $Gpc h^{-1}$ MultiDark Planck2 simulation. To compare with observational results, we extend the criterion 1 (a Ly$α$ luminosity limited sample), to criterion 2 (a match to the observed mean galaxy number density). To further statistically study the finding efficiency of this method, we enlarge the identified protocluster sample (criterion 3) to about 3500 at $z=5.7$ and study their final mass distribution. The number of overdense regions and their selection probability depends on the semi-analytic models and strongly on the three selection criteria (partly by design). The protoclusters identified with criterion 1 are associated with a typical final cluster mass of $2.82\pm0.92 \times 10^{15} M_\odot$ which is in agreement with the prediction (within $\pm 1 σ$) of an observed massive protocluster at $z=5.7$. Identifying more protoclusters allows us to investigate the efficiency of this method, which is more suitable for identifying the most massive clusters: completeness ($\mathbb{C}$) drops rapidly with decreasing halo mass. We further find that it is hard to have a high purity ($\mathbb{P}$) and completeness simultaneously.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Three Hundred Project: Backsplash galaxies in simulations of clusters

In the outer regions of a galaxy cluster, galaxies may be either falling into the cluster for the first time, or have already passed through the cluster centre at some point in their past. To investigate these two distinct populations, we utilise TheThreeHundred project, a suite of 324 hydrodynamical resimulations of galaxy clusters. In particular, we study the 'backsplash population' of galaxies; those that have passed within $R_{200}$ of the cluster centre at some time in their history, but are now outside of this radius. We find that, on average, over half of all galaxies between $R_{200}$ and $2R_{200}$ from their host at $z=0$ are backsplash galaxies, but that this fraction is dependent on the dynamical state of a cluster, as dynamically relaxed clusters have a greater backsplash fraction. We also find that this population is mostly developed at recent times ($z\leq0.4$), and is dependent on the recent history of a cluster. Finally, we show that the dynamical state of a given cluster, and thus the fraction of backsplash galaxies in its outskirts, can be predicted based on observational properties of the cluster.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Three Hundred Project: quest of clusters of galaxies morphology and dynamical state through Zernike Polynomials

The knowledge of the dynamical state of galaxy clusters allows to alleviate systematics when observational data from these objects are applied in cosmological studies. Evidence of correlation between the state and the morphology of the clusters is well studied. The morphology can be inferred by images of the surface brightness in the X-ray band and of the thermal component of the Sunyaev-Zel'dovich (tSZ) effect in the millimetre range. For this purpose, we apply, for the first time, the Zernike polynomial decomposition, a common analytic approach mostly used in adaptive optics to recover aberrated radiation wavefronts at the telescopes pupil plane. With this novel way we expect to correctly infer the morphology of clusters and so possibly, their dynamical state. To verify the reliability of this new approach we use more than 300 synthetic clusters selected in THE THREE HUNDRED project at different redshifts ranging from 0 up to 1.03. Mock maps of the tSZ, quantified with the Compton parameter, $y$-maps, are modelled with Zernike polynomials inside $R_{500}$, the cluster reference radius. We verify that it is possible to discriminate the morphology of each cluster by estimating the contribution of the different polynomials to the fit of the map. The results of this new method are correlated with those of a previous analysis made on the same catalogue, using two parameters that combine either morphological or dynamical-state probes. We underline that instrumental angular resolution of the maps has an impact mainly when we extend this approach to high-redshift clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Three Hundred project: shapes and radial alignment of satellite, infalling, and backsplash galaxies

Using 324 numerically modelled galaxy clusters we investigate the radial and galaxy-halo alignment of dark matter subhaloes and satellite galaxies orbiting within and around them. We find that radial alignment depends on distance to the centre of the galaxy cluster but appears independent of the dynamical state of the central host cluster. Furthermore, we cannot find a relation between radial alignment of the halo or galaxy shape with its own mass. We report that backsplash galaxies, i.e. objects that have already passed through the cluster radius but are now located in the outskirts, show a stronger radial alignment than infalling objects. We further find that there exists a population of well radially aligned objects passing very close to the central cluster's centre which were found to be on highly radial orbit.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Three Hundred Project: the stellar and gas profiles

Using the catalogues of galaxy clusters from The Three Hundred project, modelled with both hydrodynamic simulations, (Gadget-X and Gadget-MUSIC), and semi-analytic models (SAMs), we study the scatter and self-similarity of the profiles and distributions of the baryonic components of the clusters: the stellar and gas mass, metallicity, the stellar age, gas temperature, and the (specific) star formation rate. Through comparisons with observational results, we find that the shape and the scatter of the gas density profiles matches well the observed trends including the reduced scatter at large radii which is a signature of self-similarity suggested in previous studies. One of our simulated sets, Gadget-X, reproduces well the shape of the observed temperature profile, while Gadget-MUSIC has a higher and flatter profile in the cluster centre and a lower and steeper profile at large radii. The gas metallicity profiles from both simulation sets, despite following the observed trend, have a relatively lower normalisation. The cumulative stellar density profiles from SAMs are in better agreement with the observed result than both hydrodynamic simulations which show relatively higher profiles. The scatter in these physical profiles, especially in the cluster centre region, shows a dependence on the cluster dynamical state and on the cool-core/non-cool-core dichotomy. The stellar age, metallicity and (s)SFR show very large scatter, which are then presented in 2D maps. We also do not find any clear radial dependence of these properties. However, the brightest central galaxies have distinguishable features compared to the properties of the satellite galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

nIFTY galaxy cluster simulations III: The Similarity & Diversity of Galaxies & Subhaloes

We examine subhaloes and galaxies residing in a simulated LCDM galaxy cluster ($M^{\rm crit}_{200}=1.1\times10^{15}M_\odot/h$) produced by hydrodynamical codes ranging from classic Smooth Particle Hydrodynamics (SPH), newer SPH codes, adaptive and moving mesh codes. These codes use subgrid models to capture galaxy formation physics. We compare how well these codes reproduce the same subhaloes/galaxies in gravity only, non-radiative hydrodynamics and full feedback physics runs by looking at the overall subhalo/galaxy distribution and on an individual objects basis. We find the subhalo population is reproduced to within $\lesssim10\%$ for both dark matter only and non-radiative runs, with individual objects showing code-to-code scatter of $\lesssim0.1$ dex, although the gas in non-radiative simulations shows significant scatter. Including feedback physics significantly increases the diversity. Subhalo mass and $V_{max}$ distributions vary by $\approx20\%$. The galaxy populations also show striking code-to-code variations. Although the Tully-Fisher relation is similar in almost all codes, the number of galaxies with $10^{9}M_\odot/h\lesssim M_*\lesssim 10^{12}M_\odot/h$ can differ by a factor of 4. Individual galaxies show code-to-code scatter of $\sim0.5$ dex in stellar mass. Moreover, strong systematic differences exist, with some codes producing galaxies $70\%$ smaller than others. The diversity partially arises from the inclusion/absence of AGN feedback. Our results combined with our companion papers demonstrate that subgrid physics is not just subject to fine-tuning, but the complexity of building galaxies in all environments remains a challenge. We argue even basic galaxy properties, such as the stellar mass to halo mass, should be treated with errors bars of $\sim0.2-0.4$ dex.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

nIFTy Galaxy Cluster simulations IV: Quantifying the Influence of Baryons on Halo Properties

Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern SPH approaches. For each code represented we have a dark matter only (DM) and non-radiative (NR) version of the cluster, as well as a full physics (FP) version for a subset of the codes. We compare both radial mass and kinematic profiles, as well as global measures of the cluster (e.g. concentration, spin, shape), in the NR and FP runs with that in the DM runs. Our analysis reveals good consistency (<= 20%) between global properties of the cluster predicted by different codes when integrated quantities are measured within the virial radius R200. However, we see larger differences for quantities within R2500, especially in the FP runs. The radial profiles reveal a diversity, especially in the cluster centre, between the NR runs, which can be understood straightforwardly from the division of codes into classic SPH and non-classic SPH (including the modern SPH, adaptive and moving mesh codes); and between the FP runs, which can also be understood broadly from the division of codes into those that include AGN feedback and those that do not. The variation with respect to the median is much larger in the FP runs with different baryonic physics prescriptions than in the NR runs with different hydrodynamics solvers.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

nIFTy galaxy cluster simulations V: Investigation of the Cluster Infall Region

We examine the properties of the galaxies and dark matter haloes residing in the cluster infall region surrounding the simulated $Λ$CDM galaxy cluster studied by Elahi et al. (2016) at z=0. The $1.1\times10^{15}h^{-1}\text{M}_{\odot}$ galaxy cluster has been simulated with eight different hydrodynamical codes containing a variety of hydrodynamic solvers and subgrid schemes. All models completed a dark-matter only, non-radiative and full-physics run from the same initial conditions. The simulations contain dark matter and gas with mass resolution $m_{\text{DM}}=9.01\times 10^8h^{-1}\text{M}_{\odot}$ and $m_{\text{gas}}=1.9\times 10^8h^{-1}\text{M}_{\odot}$ respectively. We find that the synthetic cluster is surrounded by clear filamentary structures that contain ~60% of haloes in the infall region with mass ~$10^{12.5} - 10^{14} h^{-1}\text{M}_{\odot}$, including 2-3 group-sized haloes ($> 10^{13}h^{-1}\text{M}_{\odot}$). However, we find that only ~10% of objects in the infall region are subhaloes residing in haloes, which may suggest that there is not much ongoing preprocessing occurring in the infall region at z=0. By examining the baryonic content contained within the haloes, we also show that the code-to-code scatter in stellar fraction across all halo masses is typically ~2 orders of magnitude between the two most extreme cases, and this is predominantly due to the differences in subgrid schemes and calibration procedures that each model uses. Models that do not include AGN feedback typically produce too high stellar fractions compared to observations by at least ~1 order of magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the dynamical state of galaxy clusters: insights from cosmological simulations II

Using a suite of cosmology simulations of a sample of $> 120$ galaxy clusters with $\log(M_{DM, vir}) \le 14.5$. We compare clusters that form in purely dark matter run and their counterparts in hydro runs and investigate 4 independent parameters, that are normally used to classify dynamical state. We find that the virial ratio $η$ in hydro-dynamical runs is $\sim 10$ per cent lower than in the DM run, and there is no clear separation between the relaxed and unrelaxed clusters for any parameter. Further, using the velocity dispersion deviation parameter $ζ$, which is defined as the ratio between cluster velocity dispersion $σ$ and the theoretical prediction $σ_t = \sqrt{G M_{total}/R}$, we find that there is a linear correlation between the virial ratio $η$ and this $ζ$ parameter. We propose to use this $ζ$ parameter, which can be easily derived from observed galaxy clusters, as a substitute of the $η$ parameter to quantify the cluster dynamical state.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

How does our choice of observable influence our estimation of the centre of a galaxy cluster? Insights from cosmological simulations

Galaxy clusters are an established and powerful test-bed for theories of both galaxy evolution and cosmology. Accurate interpretation of cluster observations often requires robust identification of the location of the centre. Using a statistical sample of clusters drawn from a suite of cosmological simulations in which we have explored a range of galaxy formation models, we investigate how the location of this centre is affected by the choice of observable - stars, hot gas, or the full mass distribution as can be probed by the gravitational potential. We explore several measures of cluster centre: the minimum of the gravitational potential, which would expect to define the centre if the cluster is in dynamical equilibrium; the peak of the density; the centre of BCG; and the peak and centroid of X-ray luminosity. We find that the centre of BCG correlates more strongly with the minimum of the gravitational potential than the X-ray defined centres, while AGN feedback acts to significantly enhance the offset between the peak X-ray luminosity and minimum gravitational potential. These results highlight the importance of centre identification when interpreting clusters observations, in particular when comparing theoretical predictions and observational data.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

nIFTy Cosmology: Comparison of Galaxy Formation Models

We present a comparison of 14 galaxy formation models: 12 different semi-analytical models and 2 halo-occupation distribution models for galaxy formation based upon the same cosmological simulation and merger tree information derived from it. The participating codes have proven to be very successful in their own right but they have all been calibrated independently using various observational data sets, stellar models, and merger trees. In this paper we apply them without recalibration and this leads to a wide variety of predictions for the stellar mass function, specific star formation rates, stellar-to- halo mass ratios, and the abundance of orphan galaxies. The scatter is much larger than seen in previous comparison studies primarily because the codes have been used outside of their native environment within which they are well tested and calibrated. The purpose of the `nIFTy comparison of galaxy formation models' is to bring together as many different galaxy formation modellers as possible and to investigate a common approach to model calibration. This paper provides a unified description for all participating models and presents the initial, uncalibrated comparison as a baseline for our future studies where we will develop a common calibration framework and address the extent to which that reduces the scatter in the model predictions seen here.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

nIFTy galaxy cluster simulations I: dark matter & non-radiative models

We have simulated the formation of a galaxy cluster in a $Λ$CDM universe using twelve different codes modeling only gravity and non-radiative hydrodynamics (\art, \arepo, \hydra\ and 9 incarnations of GADGET). This range of codes includes particle based, moving and fixed mesh codes as well as both Eulerian and Lagrangian fluid schemes. The various GADGET implementations span traditional and advanced smoothed-particle hydrodynamics (SPH) schemes. The goal of this comparison is to assess the reliability of cosmological hydrodynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. The underlying gravitational framework can be aligned very accurately for all the codes allowing a detailed investigation of the differences that develop due to the various gas physics implementations employed. As expected, the mesh-based codes ART and AREPO form extended entropy cores in the gas with rising central gas temperatures. Those codes employing traditional SPH schemes show falling entropy profiles all the way into the very centre with correspondingly rising density profiles and central temperature inversions. We show that methods with modern SPH schemes that allow entropy mixing span the range between these two extremes and the latest SPH variants produce gas entropy profiles that are essentially indistinguishable from those obtained with grid based methods.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

nIFTy galaxy cluster simulations II: radiative models

We have simulated the formation of a massive galaxy cluster (M$_{200}^{\rm crit}$ = 1.1$\times$10$^{15}h^{-1}M_{\odot}$) in a $Λ$CDM universe using 10 different codes (RAMSES, 2 incarnations of AREPO and 7 of GADGET), modeling hydrodynamics with full radiative subgrid physics. These codes include Smoothed-Particle Hydrodynamics (SPH), spanning traditional and advanced SPH schemes, adaptive mesh and moving mesh codes. Our goal is to study the consistency between simulated clusters modeled with different radiative physical implementations - such as cooling, star formation and AGN feedback. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. We find that, with respect to non-radiative simulations, dark matter is more centrally concentrated, the extent not simply depending on the presence/absence of AGN feedback. The scatter in global quantities is substantially higher than for non-radiative runs. Intriguingly, adding radiative physics seems to have washed away the marked code-based differences present in the entropy profile seen for non-radiative simulations in Sembolini et al. (2015): radiative physics + classic SPH can produce entropy cores. Furthermore, the inclusion/absence of AGN feedback is not the dividing line -as in the case of describing the stellar content- for whether a code produces an unrealistic temperature inversion and a falling central entropy profile. However, AGN feedback does strongly affect the overall stellar distribution, limiting the effect of overcooling and reducing sensibly the stellar fraction.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The source-lens clustering effect in the context of lensing tomography and its self-calibration

Cosmic shear can only be measured where there are galaxies. This source-lens clustering (SLC) effect has two sources, intrinsic source clustering and cosmic magnification (magnification/size bias). Lensing tomography can suppress the former. However, this reduction is limited by the existence of photo-z error and nonzero redshift bin width. Furthermore, SLC induced by cosmic magnification cannot be reduced by lensing tomography. Through N-body simulations, we quantify the impact of SLC on the lensing power spectrum in the context of lensing tomography. We consider both the standard estimator and the pixel-based estimator. We find that none of them can satisfactorily handle both sources of SLC. (1) For the standard estimator, SLC induced by both sources can bias the lensing power spectrum by O(1)-O(10)%. Intrinsic source clustering also increases statistical uncertainties in the measured lensing power spectrum. However, the standard estimator suppresses intrinsic source clustering in the cross-spectrum. (2) In contrast, the pixel-based estimator suppresses SLC through cosmic magnification. However, it fails to suppress SLC through intrinsic source clustering and the measured lensing power spectrum can be biased low by O(1)-O(10)%. In short, for typical photo-z errors (sigma_z/(1+z)=0.05) and photo-z bin sizes (Delta_z^P=0.2), SLC alters the lensing E-mode power spectrum by 1-10%, with ell~10^3$ and z_s~1 being of particular interest to weak lensing cosmology. Therefore the SLC is a severe systematic for cosmology in Stage-IV lensing surveys. We present useful scaling relations to self-calibrate the SLC effect.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The effect of AGN feedback on the halo mass function

[Abridged.] We investigate baryon effects on the halo mass function (HMF), with emphasis on the role played by AGN feedback. Halos are identified with both Friends-of-Friends (FoF) and Spherical Overdensity (SO) algorithms. We embed the standard SO algorithm into a memory-controlled frame program and present the {\bf P}ython spher{\bf I}c{\bf A}l {\bf O}verdensity code --- {\small PIAO}. For both FoF and SO halos, the effect of AGN feedback is that of suppressing the HMFs to a level even below that of Dark Matter simulations. The ratio between the HMFs in the AGN and in the DM simulations is $\sim 0.8$ at overdensity $Δ_c=500$, a difference that increases at higher overdensity $Δ_c=2500$, with no significant redshift and mass dependence. A decrease of the halo masses ratio with respect to the DM case induces the decrease of the HMF in the AGN simulation. The shallower inner density profiles of halos in the AGN simulation witnesses that mass reduction is induced by the sudden displacement of gas induced by thermal AGN feedback. We provide fitting functions to describe halo mass variations at different overdensities, which can recover the HMFs with a residual random scatter $\lt 5$ per cent for halo masses larger than $10^{13} ~h^{-1}{\rm M_\odot}$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Brightest cluster galaxies in cosmological simulations: achievements and limitations of AGN feedback models

We analyze the basic properties of Brightest Cluster Galaxies (BCGs) produced by state of the art cosmological zoom-in hydrodynamical simulations. These simulations have been run with different sub-grid physics included. Here we focus on the results obtained with and without the inclusion of the prescriptions for supermassive black hole (SMBH) growth and of the ensuing Active Galactic Nuclei (AGN) feedback. The latter process goes in the right direction of decreasing significantly the overall formation of stars. However, BCGs end up still containing too much stellar mass, a problem that increases with halo mass, and having an unsatisfactory structure. This is in the sense that their effective radii are too large, and that their density profiles feature a flattening on scales much larger than observed. We also find that our model of thermal AGN feedback has very little effect on the stellar velocity dispersions, which turn out to be very large. Taken together, these problems, which to some extent can be recognized also in other numerical studies typically dealing with smaller halo masses, indicate that on one hand present day sub-resolution models of AGN feedback are not effective enough in diminishing the global formation of stars in the most massive galaxies, but on the other hand they are relatively too effective in their centers. It is likely that a form of feedback generating large scale gas outflows from BCGs precursors, and a more widespread effect over the galaxy volume, can alleviate these difficulties.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Characterizing Diffused Stellar Light in simulated galaxy clusters

[Abridged] In this paper, we carry out a detailed analysis of the performance of two different methods to identify the diffuse stellar light in cosmological hydrodynamical simulations of galaxy clusters. One method is based on a dynamical analysis of the stellar component. The second method is closer to techniques commonly employed in observational studies. Both the dynamical method and the method based on the surface brightness limit criterion are applied to the same set of hydrodynamical simulations for a large sample about 80 galaxy clusters. We find significant differences between the ICL and DSC fractions computed with the two corresponding methods, which amounts to about a factor of two for the AGN simulations, and a factor of four for the CSF set. We also find that the inclusion of AGN feedback boosts the DSC and ICL fractions by a factor of 1.5-2, respectively, while leaving the BCG+ICL and BCG+DSC mass fraction almost unchanged. The sum of the BCG and DSC mass stellar mass fraction is found to decrease from ~80 per cent in galaxy groups to ~60 per cent in rich clusters, thus in excess of what found from observational analysis. We identify the average surface brightness limits that yields the ICL fraction from the SBL method close to the DSC fraction from the dynamical method. These surface brightness limits turn out to be brighter in the CSF than in the AGN simulations. This is consistent with the finding that AGN feedback makes BCGs to be less massive and with shallower density profiles than in the CSF simulations. The BCG stellar component, as identified by both methods, are slightly older and more metal-rich than the stars in the diffuse component.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nonlinearities in modified gravity cosmology. II. Impacts of modified gravity on the halo properties

The statistics of dark matter halos is an essential component of understanding the nonlinear evolution in modified gravity cosmology. Based on a series of modified gravity N-body simulations, we investigate the halo mass function, concentration and bias. We model the impact of modified gravity by a single parameter ζ, which determines the enhancement of particle acceleration with respect to GR, given the identical mass distribution (ζ=1 in GR). We select snapshot redshifts such that the linear matter power spectra of different gravity models are identical, in order to isolate the impact of gravity beyond modifying the linear growth rate. At the baseline redshift corresponding to z_S=1.2 in the standard ΛCDM, for a 10% deviation from GR(|ζ-1|=0.1), the measured halo mass function can differ by about 5-10%, the halo concentration by about 10-20%, while the halo bias differs significantly less. These results demonstrate that the halo mass function and/or the halo concentration are sensitive to the nature of gravity and may be used to make interesting constraints along this line.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gaussianizing the non-Gaussian lensing convergence field II: the applicability to noisy data

In paper I (Yu et al. [1]), we show through N-body simulation that a local monotonic Gaussian transformation can significantly reduce non-Gaussianity in a noise-free lensing convergence field. This makes the Gaussianization a promising theoretical tool to understand high-order lensing statistics. Here we present a study of its applicability in lensing data analysis, in particular when shape measurement noise is presented in lensing convergence maps. (i) We find that shape measurement noise significantly degrades the Gaussianization performance and the degradation increases for shallower surveys. (ii) The Wiener filter is efficient in reducing the impact of shape measurement noise. The Gaussianization of the Wiener-filtered lensing maps is able to suppress skewness, kurtosis, and the 5th- and 6th-order cumulants by a factor of 10 or more. It also works efficiently to reduce the bispectrum to zero.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The effects of baryons on the halo mass function

We present an analysis of the effects of baryon physics on the halo mass function. The analysis is based on simulations of a cosmological volume. Besides a Dark Matter (DM) only simulation, we also carry out two other hydrodynamical simulations. We identified halos using a spherical overdensity algorithm and their masses are computed at three different overdensities (with respect to the critical one), $Δ_c=200$, 500 and 1500. We find the fractional difference between halo masses in the hydrodynamical and in the DM simulations to be almost constant, at least for halos more massive than $\log (M_{Δ_c} / \hMsun)\geq 13.5$. In this range, mass increase in the hydrodynamical simulations is of about 4-5 per cent at $Δ_c=500$ and $\sim 1$ - 2 per cent at $Δ_c=200$. Quite interestingly, these differences are nearly the same for both radiative and non-radiative simulations. Such variations of halo masses induce corresponding variations of the halo mass function (HMF). At $z=0$, the HMFs for GH and CSF simulations are close to the DM one, with differences of $\mincir 3$ per cent at $Δ_c = 200$, and $\simeq 7$ per cent at $Δ_c=500$, with $\sim 10$ - 20 per cent differences reached at $Δ_c = 1500$. At this higher overdensity, the increase of the HMF for the radiative case is larger by about a factor 2 with respect to the non--radiative case. Assuming a constant mass shift to rescale the HMF from the hydrodynamic to the DM simulations, brings the HMF difference with respect to the DM case to be consistent with zero. Our results have interesting implications to bracket uncertainties in the mass function calibration associated to the uncertain baryon physics, in view of cosmological applications of future large surveys of galaxy clusters. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The halo mass function in interacting Dark Energy models

We present a detailed investigation of the effects that a direct interaction between Dark Energy (DE) and Cold Dark Matter (CDM) particles imprints on the Halo Mass Function (HMF) of groups and clusters of galaxies. Making use of the public halo catalogs of the {\small CoDECS} simulations, we derive the HMF for several different types of coupled DE scenarios both based on the FoF algorithm and on the SO halo identification for different values of the overdensity threshold $Δ_{c}$. We compare the computed HMFs for coupled DE cosmologies with $Λ$CDM as well as with the predictions of the standard analytic fitting functions. Our results show that the standard fitting functions still reproduce reasonably well both the FoF and the SO HMFs of interacting DE cosmologies at intermediate masses and at low redshifts, once rescaled to the characteristic amplitude of linear density perturbations of each specific model as given by $σ_{8}$. However, we also find that such apparent degeneracy with $σ_{8}$ is broken both by the high-mass tail and by the redshift evolution of our HMFs, with deviations beyond $\sim 10%$ for most of the models under investigation. Furthermore, the discrepancy with respect to the predictions of standard fitting functions rescaled with the characteristic value of $σ_{8}$ shows -- for some models -- a strong dependence on the spherical overdensity threshold $Δ_{c}$ used for the halo identification. We find that such effect is due to a significant increase of halo concentration at low redshifts in these models, that is however absent in the majority of the cosmological scenarios considered in this work. We can therefore conclude that the universality of the HMF is violated by cosmological models that feature a direct interaction between DE and CDM.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Gaussianizing the non-Gaussian lensing convergence field I: the performance of the Gaussianization

Motivated by recent works of Neyrinck et al. 2009 and Scherrer et al. 2010, we proposed a Gaussianization transform to Gaussianize the non-Gaussian lensing convergence field $κ$. It performs a local monotonic transformation $κ\rightarrow y$ pixel by pixel to make the unsmoothed one-point probability distribution function of the new variable $y$ Gaussian. We tested whether the whole $y$ field is Gaussian against N-body simulations. (1) We found that the proposed Gaussianization suppresses the non-Gaussianity by orders of magnitude, in measures of the skewness, the kurtosis, the 5th- and 6th-order cumulants of the $y$ field smoothed over various angular scales relative to that of the corresponding smoothed $κ$ field. The residual non-Gaussianities are often consistent with zero within the statistical errors. (2) The Gaussianization significantly suppresses the bispectrum. Furthermore, the residual scatters around zero, depending on the configuration in the Fourier space. (3) The Gaussianization works with even better performance for the 2D fields of the matter density projected over $\sim 300 \mpch$ distance interval centered at $z\in(0,2)$, which can be reconstructed from the weak lensing tomography. (4) We identified imperfectness and complexities of the proposed Gaussianization. We noticed weak residual non-Gaussianity in the $y$ field. We verified the widely used logarithmic transformation as a good approximation to the Gaussianization transformation. However, we also found noticeable deviations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Properties of fossil groups in cosmological simulations and galaxy formation models

It has been a long-standing question whether fossil groups are just sampling the tail of the distribution of ordinary groups, or whether they are a physically distinct class of objects, characterized by an unusual and special formation history. To study this question, we here investigate fossil groups identified in the hydrodynamical simulations of the GIMIC project, which consists of resimulations of five regions in the Millennium Simulation (MS) that are characterized by different large-scale densities, ranging from a deep void to a proto-cluster region. For comparison, we also consider semi-analytic models built on top of the MS, as well as a conditional luminosity function approach. We identify galaxies in the GIMIC simulations as groups of stars and use a spectral synthesis code to derive their optical properties. The X-ray luminosity of the groups is estimated in terms of the thermal bremsstrahlung emission of the gas in the host halos, neglecting metallicity effects. We focus on comparing the properties of fossil groups in the theoretical models and observational results, highlighting the differences between them, and trying to identify possible dependencies on environment for which our approach is particularly well set-up. We find that the optical fossil fraction in all of our theoretical models declines with increasing halo mass, and there is no clear environmental dependence. Combining the optical and X-ray selection criteria for fossil groups, the halo mass dependence of the fossil groups seen in optical vanishes. Over the GIMIC halo mass range we resolve best, 9.0\times1012 \sim 4.0\times1013 h-1 M, the central galaxies in the fossil groups show similar properties as those in ordinary groups, in terms of age, metallicity, color, concentration, and mass-to-light ratio. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Nonlinearities in modified gravity cosmology I: signatures of modified gravity in the nonlinear matter power spectrum

A large fraction of cosmological information on dark energy and gravity is encoded in the nonlinear regime. Precision cosmology thus requires precision modeling of nonlinearities in general dark energy and modified gravity models. We modify the Gadget-2 code and run a series of N-body simulations on modified gravity cosmology to study the nonlinearities. The modified gravity model that we investigate in the present paper is characterized by a single parameter ζ, which determines the enhancement of particle acceleration with respect to general relativity (GR), given the identical mass distribution (ζ= 1 in GR). The first nonlinear statistics we investigate is the nonlinear matter power spectrum at k < 3h/Mpc, which is the relevant range for robust weak lensing power spectrum modeling at l < 2000. In this study, we focus on the relative difference in the nonlinear power spectra at corresponding redshifts where different gravity models have the same linear power spectra. This particular statistics highlights the imprint of modified gravity in the nonlinear regime and the importance to include the nonlinear regime in testing GR. By design, it is less susceptible to the sample variance and numerical artifacts. We adopt a mass assignment method based on wavelet to improve the power spectrum measurement. We run a series of tests to determine the suitable simulation specifications (particle number, box size and initial redshift). We find that, the nonlinear power spectra can differ by ~30% for 10% deviation from GR (|ζ-1| = 0.1) where the rms density fluctuations reach 10. This large difference, on one hand, shows the richness of information on gravity in the corresponding scales, and on the other hand, invalidates simple extrapolations of some existing fitting formulae to modified gravity cosmology.

0 citations0 reviewsNo code linkedOpen access