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

S. Borgani

S. Borgani 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-phastro-ph.HEastro-ph.IMgr-qcastro-ph.SR

Catalog footprint

What is connected

79works
7topics
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

79 published item(s)

preprint2026arXiv

Euclid preparation. Calibrated intrinsic galaxy alignments in the Euclid Flagship simulation

Intrinsic alignments of galaxies are potentially a major contaminant of cosmological analyses of weak gravitational lensing. We construct a semi-analytic model of galaxy ellipticities and alignments in the \Euclid Flagship simulation to predict this contamination in Euclid's weak lensing observations. Galaxy shapes and orientations are determined by the corresponding properties of the host haloes in the underlying $N$-body simulation, as well as the relative positions of galaxies within their halo. Alignment strengths are moderated via stochastic misalignments, separately for central and satellite galaxies and conditional on the galaxy's redshift, luminosity, and rest-frame colour. The resulting model is calibrated against galaxy ellipticity statistics from the COSMOS Survey, selected alignment measurements based on Sloan Digital Sky Survey samples, and galaxy orientations extracted from the Horizon-AGN hydrodynamic simulation at redshift $z=1$. The best-fit model has a total of 12 alignment parameters and generally reproduces the calibration data sets well within the $1σ$ statistical uncertainties of the observations and the \flagship simulation, with notable exceptions for the most luminous sub-samples on small physical scales. The statistical power of the calibration data and the volume of the single \flagship realisation are still too small to provide informative prior ranges for intrinsic alignment amplitudes in relevant galaxy samples. As a first application, we predict that \Euclid end-of-mission tomographic weak gravitational lensing two-point statistics are modified by up to order $10\,\%$ due to intrinsic alignments.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Euclid preparation. Galaxy 2-point correlation function modelling in redshift space

The Euclid satellite will measure spectroscopic redshifts for tens of millions of emission-line galaxies. In the context of Stage-IV surveys, the 3-dimensional clustering of galaxies plays a key role in providing cosmological constraints. In this paper, we conduct a model comparison for the multipole moments of the galaxy 2-point correlation function (2PCF) in redshift space. We test state-of-the-art models, in particular the effective field theory of large-scale structure (EFT), one based on the velocity difference generating function (VDG$_{\infty}$), and different variants of Lagrangian perturbation theory (LPT) models, such as convolutional LPT (CLPT) and its effective-field-theory extension (CLEFT). We analyse the first three even multipoles of the 2PCF in the Flagship 1 simulation, which consists of four snapshots at $z\in\{0.9,1.2,1.5,1.8\}$. We study both template-fitting and full-shape approaches and find that with the template-fitting approach, only the VDG$_{\infty}$ model is able to reach a minimum fitting scale of $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ at $z=0.9$ without biasing the recovered parameters. Indeed, the EFT model becomes inaccurate already at $s_{\rm min}=30\,h^{-1}\,{\rm Mpc}$. Conversely, in the full-shape analysis, the CLEFT and VDG$_{\infty}$ models perform similarly well, but only the CLEFT model can reach $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ while the VDG$_{\infty}$ model is unbiased down to $s_{\rm min}=25\,h^{-1}\,{\rm Mpc}$ at the lowest redshift. Overall, in order to achieve the accuracy required by Euclid, non-perturbative modelling such as in the VDG$_{\infty}$ or CLEFT models should be considered. At $z=1.8$, the CLPT model is sufficient to describe the data with high figure of merit. This comparison selects baseline models that perform best in ideal conditions and sets the stage for an optimal analysis of Euclid data in configuration space.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Euclid preparation. Testing analytic models of galaxy intrinsic alignments in the Euclid Flagship simulation

We model intrinsic alignments (IA) in Euclid's Flagship simulation to investigate its impact on Euclid's weak lensing signal. Our IA implementation in the Flagship simulation takes into account photometric properties of galaxies as well as their dark matter host halos. We compare simulations against theory predictions, determining the parameters of two of the most widely used IA models: the Non Linear Alignment (NLA) and the Tidal Alignment and Tidal Torquing (TATT) models. We measure the amplitude of the simulated IA signal as a function of galaxy magnitude and colour in the redshift range $0.1<z<2.1$. We find that both NLA and TATT can accurately describe the IA signal in the simulation down to scales of $6$-$7 \,h^{-1}\,$Mpc. We measure alignment amplitudes for red galaxies comparable to those of the observations, with samples not used in the calibration procedure. For blue galaxies, our constraints are consistent with zero alignments in our first redshift bin $0.1 < z < 0.3$, but we detect a non-negligible signal at higher redshift, which is, however, consistent with the upper limits set by observational constraints. Additionally, several hydrodynamical simulations predict alignment for spiral galaxies, in agreement with our findings. Finally, the evolution of alignment with redshift is realistic and comparable to that determined in the observations. However, we find that the commonly adopted redshift power-law for IA fails to reproduce the simulation alignments above $z=1.1$. A significantly better agreement is obtained when a luminosity dependence is included, capturing the intrinsic luminosity evolution with redshift in magnitude-limited surveys. We conclude that the Flagship IA simulation is a useful tool for translating current IA constraints into predictions for IA contamination of Euclid-like samples.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Euclid: Galaxy SED reconstruction in the PHZ processing function: impact on the PSF and the role of medium-band filters

Weak lensing surveys require accurate correction for the point spread function (PSF) when measuring galaxy shapes. For a diffraction-limited PSF, as arises in space-based missions, this correction depends on each galaxy SED. In the Euclid mission, galaxy SED reconstruction, a tasks of the photometric-redshift processing function (PHZ PF), relies on broad- and medium-band ancillary photometry. The limited wavelength sampling of the Euclid VIS passband and signal-to-noise ratio may affect the reconstruction accuracy and translate into biases in the weak lensing measurements. In this study, we present the methodology, which is employed in the Euclid PHZ PF, for reconstructing galaxy SEDs at 55 wavelengths, sampling the VIS passband every 10 nm, and we assess whether it fulfils the accuracy requirements imposed on the Euclid PSF model. We employ both physics- and data-driven methods, focusing on a new approach of template-based flux correction and Gaussian processes, and we introduce an SED metric whose bias propagates into PSF quadrupole moment errors. Our findings demonstrate that Gaussian processes and template fitting meet the requirements only in specific, but complementary, redshift intervals. We therefore propose a hybrid approach, which leverages both methods. This solution proves to be effective in meeting the Euclid accuracy requirements for most of the redshift range of the survey. Finally, we investigate the impact on the SED reconstruction of a new set of 16 evenly-spaced medium-band filters for the Subaru telescope, providing quasi-spectroscopic coverage of the VIS passband. This study shows promising results, ensuring accurate SED reconstruction and meeting the mission PSF requirements. This work thus provides not only the methodological foundation of galaxy SED reconstruction in the Euclid PHZ PF, but also a roadmap for future improvements using a new medium-band survey.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Euclid preparation. Simulating thousands of Euclid spectroscopic skies

We present two extensive sets of 3500+1000 simulations of dark matter haloes on the past light cone, and two corresponding sets of simulated (`mock') galaxy catalogues that represent the Euclid spectroscopic sample. The simulations were produced with the latest version of the PINOCCHIO code, and provide the largest, public set of simulated skies. Mock galaxy catalogues were obtained by populating haloes with galaxies using an halo occupation distribution (HOD) model extracted from the Flagship galaxy catalogue provided by Euclid Collaboration. The Geppetto set of 3500 simulated skies was obtained by tiling a 1.2 Gpc/h box to cover a light-cone whose sky footprint is a circle of 30 deg radius, for an area of 2763 deg$^2$ and a minimum halo mass of $1.5\times10^{11}$ Msun/h. The relatively small box size makes this set unfit for measuring very large scales. The EuclidLargeBox set consists of 1000 simulations of 3.38 Gpc/h, with the same mass resolution and a footprint that covers half of the sky, excluding the Milky Way zone of avoidance. From this we produced a set of 1000 EuclidLargeMocks on the 30 deg radius footprint, whose comoving volume is fully contained in the simulation box. We validated the two sets of catalogues by analysing number densities, power spectra, and 2-point correlation functions, showing that the Flagship spectroscopic catalogue is consistent with being one of the realisations of the simulated sets, although we noticed small deviations limited to the quadrupole at k>0.2 h/Mpc. We show cosmological parameter inference from these catalogues and demonstrate that using one realisation of EuclidLargeMocks in place of the Flagship mock produces the same posteriors, to within the expected shift given by sample variance. These simulated skies will be used for the galaxy clustering analysis of Euclid's Data Release 1 (DR1).

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Euclid preparation: XIII. Forecasts for galaxy morphology with the Euclid Survey using Deep Generative Models

We present a machine learning framework to simulate realistic galaxies for the Euclid Survey. The proposed method combines a control on galaxy shape parameters offered by analytic models with realistic surface brightness distributions learned from real Hubble Space Telescope observations by deep generative models. We simulate a galaxy field of $0.4\,\rm{deg}^2$ as it will be seen by the Euclid visible imager VIS and show that galaxy structural parameters are recovered with similar accuracy as for pure analytic Sérsic profiles. Based on these simulations, we estimate that the Euclid Wide Survey will be able to resolve the internal morphological structure of galaxies down to a surface brightness of $22.5\,\rm{mag}\,\rm{arcsec}^{-2}$, and $24.9\,\rm{mag}\,\rm{arcsec}^{-2}$ for the Euclid Deep Survey. This corresponds to approximately $250$ million galaxies at the end of the mission and a $50\,\%$ complete sample for stellar masses above $10^{10.6}\,\rm{M}_\odot$ (resp. $10^{9.6}\,\rm{M}_\odot$) at a redshift $z\sim0.5$ for the wide (resp. deep) survey. The approach presented in this work can contribute to improving the preparation of future high-precision cosmological imaging surveys by allowing simulations to incorporate more realistic galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Euclid preparation: XIX. Impact of magnification on photometric galaxy clustering

We investigate the importance of lensing magnification for estimates of galaxy clustering and its cross-correlation with shear for the photometric sample of Euclid. Using updated specifications, we study the impact of lensing magnification on the constraints and the shift in the estimation of the best fitting cosmological parameters that we expect if this effect is neglected. We follow the prescriptions of the official Euclid Fisher matrix forecast for the photometric galaxy clustering analysis and the combination of photometric clustering and cosmic shear. The slope of the luminosity function (local count slope), which regulates the amplitude of the lensing magnification, and the galaxy bias have been estimated from the Euclid Flagship simulation.We find that magnification significantly affects both the best-fit estimation of cosmological parameters and the constraints in the galaxy clustering analysis of the photometric sample. In particular, including magnification in the analysis reduces the 1$σ$ errors on $Ω_{\text{m},0}, w_{0}, w_a$ at the level of 20-35%, depending on how well we will be able to independently measure the local count slope. In addition, we find that neglecting magnification in the clustering analysis leads to shifts of up to 1.6$σ$ in the best-fit parameters. In the joint analysis of galaxy clustering, cosmic shear, and galaxy-galaxy lensing, magnification does not improve precision, but it leads to an up to 6$σ$ bias if neglected. Therefore, for all models considered in this work, magnification has to be included in the analysis of galaxy clustering and its cross-correlation with the shear signal ($3\times2$pt analysis) for an accurate parameter estimation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Euclid preparation. XVIII. The NISP photometric system

Euclid will be the first space mission to survey most of the extragalactic sky in the 0.95-2.02 $μ$m range, to a 5$σ$ point-source median depth of 24.4 AB mag. This unique photometric data set will find wide use beyond Euclid's core science. In this paper, we present accurate computations of the Euclid Y_E, J_E and H_E passbands used by the Near-Infrared Spectrometer and Photometer (NISP), and the associated photometric system. We pay particular attention to passband variations in the field of view, accounting among others for spatially variable filter transmission, and variations of the angle of incidence on the filter substrate using optical ray tracing. The response curves' cut-on and cut-off wavelengths - and their variation in the field of view - are determined with 0.8 nm accuracy, essential for the photometric redshift accuracy required by Euclid. After computing the photometric zeropoints in the AB mag system, we present linear transformations from and to common ground-based near-infrared photometric systems, for normal stars, red and brown dwarfs, and galaxies separately. A Python tool to compute accurate magnitudes for arbitrary passbands and spectral energy distributions is provided. We discuss various factors from space weathering to material outgassing that may slowly alter Euclid's spectral response. At the absolute flux scale, the Euclid in-flight calibration program connects the NISP photometric system to Hubble Space Telescope spectrophotometric white dwarf standards; at the relative flux scale, the chromatic evolution of the response is tracked at the milli-mag level. In this way, we establish an accurate photometric system that is fully controlled throughout Euclid's lifetime.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Machine Learning to identify ICL and BCG in simulated galaxy clusters

Nowadays, Machine Learning techniques offer fast and efficient solutions for classification problems that would require intensive computational resources via traditional methods. We examine the use of a supervised Random Forest to classify stars in simulated galaxy clusters after subtracting the member galaxies. These dynamically different components are interpreted as the individual properties of the stars in the Brightest Cluster Galaxy (BCG) and IntraCluster Light (ICL). We employ matched stellar catalogues (built from the different dynamical properties of BCG and ICL) of 29 simulated clusters from the DIANOGA set to train and test the classifier. The input features are cluster mass, normalized particle cluster-centric distance, and rest-frame velocity. The model is found to correctly identify most of the stars, while the larger errors are exhibited at the BCG outskirt, where the differences between the physical properties of the two components are less obvious. We investigate the robustness of the classifier to numerical resolution, redshift dependence (up to $z=1$), and included astrophysical models. We claim that our classifier provides consistent results in simulations for $z<1$, at different resolution levels and with significantly different subgrid models. The phase-space structure is examined to assess whether the general properties of the stellar components are recovered: (i) the transition radius between BCG-dominated and ICL-dominated region is identified at $0.04$ \r200; (ii) the BCG outskirt ($> 0.1$ \r200) is significantly affected by uncertainties in the classification process. In conclusion, this work suggests the importance of employing Machine Learning to speed up a computationally expensive classification in simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The 700 ks Chandra Spiderweb Field I: evidence for widespread nuclear activity in the Protocluster

(Abridged) We present an analysis of the 700 ks Chandra ACIS-S observation of the field around the Spiderweb Galaxy at z=2.156, focusing on the nuclear activity in the associated large-scale environment. We identify unresolved X-ray sources down to flux limits of 1.3X10^{-16} and 3.9X10^{-16} erg/s/cm^2 in the soft and hard band, respectively. We search for counterparts in the optical, NIR and submm bands to identify X-ray sources belonging to the protocluster. We detect 107 X-ray unresolved sources within 5 arcmin (corresponding to 2.5 Mpc) of J1140-2629, among which 13 have optical counterparts with spectroscopic redshift 2.11<z<2.20, and 1 source with photometric redshift consistent with this range. Our X-ray spectral analysis shows that their intrinsic spectral slope is consistent with an average <Γ>~1.84+-0.04. The best-fit intrinsic absorption for 5 protocluster X-ray members is N_H>10^{23} cm^{-2}, while other 6 have upper limits of the order of fewX10^{22} cm^{-2}. Two sources can only be fitted with very flat Γ<=1, and are therefore considered Compton-thick candidates. Their 0.5-10 keV rest frame luminosities are larger than 2X10^{43} erg/s, significantly greater than X-ray luminosities expected from star formation activity. The X-ray luminosity function of AGN in the volume associated to the Spiderweb protocluster in the range 10^{43}<L_X<10^{44.5} erg/s, is at least 10 times higher than that in the field at the same redshift and significantly flatter. The X-ray AGN fraction is measured to be (25.5+-4.5)% in the stellar mass range log(M*/M_sun)>10.5, corresponding to an enhancement of 6.0^{+9.0}_{-3.0} with respect to the COSMOS field at comparable redshifts and stellar mass range. We conclude that the galaxy population in the Spiderweb Protocluster is characterized by enhanced X-ray nuclear activity triggered by environmental effects on Mpc scales.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The BEHOMO project: $Λ$LTB $N$-body simulations

Our Universe may feature large-scale inhomogeneities and anisotropies which cannot be explained by the standard model of cosmology, that is, the homogeneous and isotropic FLRW metric, on which the $Λ$CDM model is built, may not describe accurately observations. Currently, there is not a satisfactory understanding of the evolution of the large-scale structure on an inhomogeneous background. We start the cosmology beyond homogeneity and isotropy (BEHOMO) project and study the inhomogeneous $Λ$LTB model with the methods of numerical cosmology. Understanding the evolution of the large-scale structure is a necessary step to constrain inhomogeneous models with present and future observables and place the standard model on more solid grounds. We perform Newtonian $N$-body simulations, whose accuracy in describing the background evolution is checked against the general relativistic solution. The large-scale structure of the corresponding $Λ$CDM simulation is also validated. We obtain the first set of simulations of the $Λ$LTB model ever produced. The data products consist of 11 snapshots between redshift 0 and 3.7 for each of the 68 simulations that have been performed, together with halo catalogs and lens planes relative to 21 snapshots, between redshift 0 and 4.2, for a total of approximately 180 TB of data. We plan to study the growth of perturbations at the linear and nonlinear level, gravitational lensing, cluster abundances and proprieties. Data can be obtained upon request. Further information is available at valerio-marra.github.io/BEHOMO-project .

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

X-ray Emission and Radio Emission from the Jets and Lobes of the Spiderweb Radio Galaxy

Deep Chandra and VLA imaging reveals a clear correlation between X-ray and radio emission on scales $\sim 100$~kpc in the Spiderweb radio galaxy at z=2.16. The X-ray emission associated with the extended radio source is likely dominated by inverse Compton up-scattering of cosmic microwave background photons by the radio emitting relativistic electrons. For regions dominated by high surface brightness emission, such as hot spots and jet knots, the implied magnetic fields are $\sim 50~μ$G to $70~μ$G. The non-thermal pressure is these brighter regions is then $\sim 9\times 10^{-10}$ dyne cm$^{-2}$, or three times larger than the non-thermal pressure derived assuming minimum energy conditions, and an order of magnitude larger than the thermal pressure in the ambient cluster medium. Assuming ram pressure confinement implies an average advance speed for the radio source of $\sim 2400$ km s$^{-1}$, and a source age of $\sim 3\times 10^7$ years. Considering the lower surface brightness, diffuse radio emitting regions, we identify an evacuated cavity in the Ly$α$ emission coincident with the tail of the eastern radio lobe. Making reasonable assumptions for the radio spectrum, we find that the relativistic electrons and fields in the lobe are plausibly in pressure equilibrium with the thermal gas, and close to a minimum energy configuration. The radio morphology suggests that the Spiderweb is a high-$z$ example of the rare class of hybrid morphology radio sources (or HyMoRS), which we attribute to interaction with the asymmetric gaseous environment indicated by the Ly$α$ emission.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

$Euclid$ preparation: XV. Forecasting cosmological constraints for the $Euclid$ and CMB joint analysis

The combination and cross-correlation of the upcoming $Euclid$ data with cosmic microwave background (CMB) measurements is a source of great expectation since it will provide the largest lever arm of epochs, ranging from recombination to structure formation across the entire past light cone. In this work, we present forecasts for the joint analysis of $Euclid$ and CMB data on the cosmological parameters of the standard cosmological model and some of its extensions. This work expands and complements the recently published forecasts based on $Euclid$-specific probes, namely galaxy clustering, weak lensing, and their cross-correlation. With some assumptions on the specifications of current and future CMB experiments, the predicted constraints are obtained from both a standard Fisher formalism and a posterior-fitting approach based on actual CMB data. Compared to a $Euclid$-only analysis, the addition of CMB data leads to a substantial impact on constraints for all cosmological parameters of the standard $Λ$-cold-dark-matter model, with improvements reaching up to a factor of ten. For the parameters of extended models, which include a redshift-dependent dark energy equation of state, non-zero curvature, and a phenomenological modification of gravity, improvements can be of the order of two to three, reaching higher than ten in some cases. The results highlight the crucial importance for cosmological constraints of the combination and cross-correlation of $Euclid$ probes with CMB data.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Euclid preparation: I. The Euclid Wide Survey

Euclid is an ESA mission designed to constrain the properties of dark energy and gravity via weak gravitational lensing and galaxy clustering. It will carry out a wide area imaging and spectroscopy survey (EWS) in visible and near-infrared, covering roughly 15,000 square degrees of extragalactic sky on six years. The wide-field telescope and instruments are optimized for pristine PSF and reduced straylight, producing very crisp images. This paper presents the building of the Euclid reference survey: the sequence of pointings of EWS, Deep fields, Auxiliary fields for calibrations, and spacecraft movements followed by Euclid as it operates in a step-and-stare mode from its orbit around the Lagrange point L2. Each EWS pointing has four dithered frames; we simulate the dither pattern at pixel level to analyse the effective coverage. We use up-to-date models for the sky background to define the Euclid region-of-interest (RoI). The building of the reference survey is highly constrained from calibration cadences, spacecraft constraints and background levels; synergies with ground-based coverage are also considered. Via purposely-built software optimized to prioritize best sky areas, produce a compact coverage, and ensure thermal stability, we generate a schedule for the Auxiliary and Deep fields observations and schedule the RoI with EWS transit observations. The resulting reference survey RSD_2021A fulfills all constraints and is a good proxy for the final solution. Its wide survey covers 14,500 square degrees. The limiting AB magnitudes ($5σ$ point-like source) achieved in its footprint are estimated to be 26.2 (visible) and 24.5 (near-infrared); for spectroscopy, the H$_α$ line flux limit is $2\times 10^{-16}$ erg cm$^{-2}$ s$^{-1}$ at 1600 nm; and for diffuse emission the surface brightness limits are 29.8 (visible) and 28.4 (near-infrared) mag arcsec$^{-2}$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Calibration of systematics in constraining modified gravity models with galaxy cluster mass profiles

Joint lensing and dynamical mass profile determinations of galaxy clusters are an excellent tool to constrain modification of gravity at cosmological scales. However, search for tiny departures from General Relativity calls for an accurate control of the systematics affecting the method. In this analysis we concentrate on the systematics in the reconstruction of mass profiles from the dynamics of cluster member galaxies, while assuming that lensing provides unbiased mass profile reconstructions. In particular, in the case study of linear $f(R)$ gravity, we aim at veryfying whether in realistic simulations of cluster formation a spurious detection of departure from GR can be detected due to violation of the main assumptions (e.g. dynamical equilibrium and spherical symmetry) on which the method is based. We aim at identifying and calibrating the impact of those systematics by analyzing a set of Dark Matter halos taken from $Λ$CDM N-body cosmological simulations performed with the GADGET-3 code. [...] If no selection criteria are applied, $\sim 60\%$ of clusters in a $Λ$CDM Universe (where GR is assumed) produce a spurious detection of modified gravity. We find that the probability of finding cluster in agreement with GR predictions $P_{GR}$ mainly depends on the properties of the halo's projected phase-space and on shape orientation of the cluster along the line-of-sight projection. We define two observational criteria which correlate with the probability to find clusters in agreement with GR predictions and which can be used to select [...] those objects that are more suitable for the application of the proposed method. In particular, we find that according to these criteria the percentage of spurious detection can be lowered down to $\sim 20\%$ in the best case. Our results are relevant in view of data that will be available with the next generation surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Three Hundred Project: correcting for the hydrostatic-equilibrium mass bias in X-ray and SZ surveys

Accurate and precise measurements of masses of galaxy clusters are key to derive robust constraints on cosmological parameters. Rising evidence from observations, however, confirms that X-ray masses, obtained under the assumption of hydrostatic equilibrium, might be underestimated, as previously predicted by cosmological simulations. We analyse more than 300 simulated massive clusters, from `The Three Hundred Project', and investigate the connection between mass bias and several diagnostics extracted from synthetic X-ray images of these simulated clusters. We find that the azimuthal scatter measured in 12 sectors of the X-ray flux maps is a statistically significant indication of the presence of an intrinsic (i.e. 3D) clumpy gas distribution. We verify that a robust correction to the hydrostatic mass bias can be inferred when estimates of the gas inhomogeneity from X-ray maps (such as the azimuthal scatter or the gas ellipticity) are combined with the asymptotic external slope of the gas density or pressure profiles, which can be respectively derived from X-ray and millimetric (Sunyaev-Zeldovich effect) observations. We also obtain that mass measurements based on either gas density and temperature or gas density and pressure result in similar distributions of the mass bias. In both cases, we provide corrections that help reduce both the dispersion and skewness of the mass bias distribution. These are effective even when irregular clusters are included leading to interesting implications for the modelling and correction of hydrostatic mass bias in cosmological analyses of current and future X-ray and SZ cluster surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

CLASH-VLT: Testing the Nature of Gravity with Galaxy Cluster Mass Profiles

We use high-precision kinematic and lensing measurements of the total mass profile of the dynamically relaxed galaxy cluster MACS J1206.2-0847 at $z=0.44$ to estimate the value of the ratio $η=Ψ/Φ$ between the two scalar potentials in the linear perturbed Friedmann-Lemaitre-Robertson-Walker metric.[...] Complementary kinematic and lensing mass profiles were derived from exhaustive analyses using the data from the Cluster Lensing And Supernova survey with Hubble (CLASH) and the spectroscopic follow-up with the Very Large Telescope (CLASH-VLT). Whereas the kinematic mass profile tracks only the time-time part of the perturbed metric (i.e. only $Φ$), the lensing mass profile reflects the contribution of both time-time and space-space components (i.e. the sum $Φ+Ψ$). We thus express $η$ as a function of the mass profiles and perform our analysis over the radial range $0.5\,Mpc\le r\le r_{200}=1.96\,Mpc$. Using a spherical Navarro-Frenk-White mass profile, which well fits the data, we obtain $η(r_{200})=1.01\,_{-0.28}^{+0.31}$ at the 68\% C.L. We discuss the effect of assuming different functional forms for mass profiles and of the orbit anisotropy in the kinematic reconstruction. Interpreting this result within the well-studied $f(R)$ modified gravity model, the constraint on $η$ translates into an upper bound to the interaction length (inverse of the scalaron mass) smaller than 2 Mpc. This tight constraint on the $f(R)$ interaction range is however substantially relaxed when systematic uncertainties in the analysis are considered. Our analysis highlights the potential of this method to detect deviations from general relativity, while calling for the need of further high-quality data on the total mass distribution of clusters and improved control on systematic effects.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical simulations challenged on the prediction of massive subhalo abundance in galaxy clusters: the case of Abell 2142

In this Letter we compare the abundance of member galaxies of a rich, nearby ($z=0.09$) galaxy cluster, Abell 2142, with that of halos of comparable virial mass extracted from sets of state-of-the-art numerical simulations, both collisionless at different resolutions and with the inclusion of baryonic physics in the form of cooling, star formation, and feedback by active galactic nuclei. We also use two semi-analytical models to account for the presence of orphan galaxies. The photometric and spectroscopic information, taken from the Sloan Digital Sky Survey Data Release 12 (SDSS DR12) database, allows us to estimate the stellar velocity dispersion of member galaxies of Abell 2142. This quantity is used as proxy for the total mass of secure cluster members and is properly compared with that of subhalos in simulations. We find that simulated halos have a statistically significant ($\gtrsim 7$ sigma confidence level) smaller amount of massive (circular velocity above $200\,{\rm km\, s^{-1}}$) subhalos, even before accounting for the possible incompleteness of observations. These results corroborate the findings from a recent strong lensing study of the Hubble Frontier Fields galaxy cluster MACS J0416 \citep{grillo2015} and suggest that the observed difference is already present at the level of dark matter (DM) subhalos and is not solved by introducing baryonic physics. A deeper understanding of this discrepancy between observations and simulations will provide valuable insights into the impact of the physical properties of DM particles and the effect of baryons on the formation and evolution of cosmological structures.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the nature of hydrostatic equilibrium in galaxy clusters

In this paper we investigate the level of hydrostatic equilibrium (HE) in the intra-cluster medium of simulated galaxy clusters, extracted from state-of-the-art cosmological hydrodynamical simulations performed with the Smoothed-Particle-Hydrodynamic code GADGET-3. These simulations include several physical processes, among which stellar and AGN feedback, and have been performed with an improved version of the code that allows for a better description of hydrodynamical instabilities and gas mixing processes. Evaluating the radial balance between the gravitational and hydrodynamical forces, via the gas accelerations generated, we effectively examine the level of HE in every object of the sample, its dependence on the radial distance from the center and on the classification of the cluster in terms of either cool-coreness or dynamical state. We find an average deviation of 10-20% out to the virial radius, with no evident distinction between cool-core and non-cool-core clusters. Instead, we observe a clear separation between regular and disturbed systems, with a more significant deviation from HE for the disturbed objects. The investigation of the bias between the hydrostatic estimate and the total gravitating mass indicates that, on average, this traces very well the deviation from HE, even though individual cases show a more complex picture. Typically, in the radial ranges where mass bias and deviation from HE are substantially different, the gas is characterized by a significant amount of random motions (>~30 per cent), relative to thermal ones. As a general result, the HE-deviation and mass bias, at given interesting distance from the cluster center, are not very sensitive to the temperature inhomogeneities in the gas.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Radiative feedback and cosmic molecular gas: the role of different radiative sources

We present results from multifrequency radiative hydrodynamical chemistry simulations addressing primordial star formation and related stellar feedback from various populations of stars, stellar energy distributions (SEDs) and initial mass functions. Spectra for massive stars, intermediate-mass stars and regular solar-like stars are adopted over a grid of 150 frequency bins and consistently coupled with hydrodynamics, heavy-element pollution and non-equilibrium species calculations. Powerful massive population III stars are found to be able to largely ionize H and, subsequently, He and He$^+$, causing an inversion of the equation of state and a boost of the Jeans masses in the early intergalactic medium. Radiative effects on star formation rates are between a factor of a few and 1 dex, depending on the SED. Radiative processes are responsible for gas heating and photoevaporation, although emission from soft SEDs has minor impacts. These findings have implications for cosmic gas preheating, primordial direct-collapse black holes, the build-up of "cosmic fossils" such as low-mass dwarf galaxies, the role of AGNi during reionization, the early formation of extended disks and angular-momentum catastrophe.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Velocity segregation effects in galaxy clusters at 0.4<~z<~1.5

Our study is meant to extend our knowledge of the galaxy color and luminosity segregation in velocity space (VCS and VLS, resp.), to clusters at intermediate and high redshift. Our sample is a collection of 41 clusters in the 0.4<~z<~1.5 redshift range, for a total of 4172 galaxies, 1674 member galaxies within 2R200 with photometric or spectroscopic information, as taken from the literature. We pay attention to perform homogeneous procedures to select cluster members, compute global cluster properties, in particular the LOS velocity dispersion sigmaV, and separate blue from red galaxies. We find evidence of VCS in clusters out to z~0.8 (at the 97%-99.99% c.l., depending on the test), in the sense that the blue galaxy population has a 10-20% larger sigmaV than the red galaxy population. Poor or no VCS is found in the High-z sample at z>=0.8. For the first time, we detect VLS in non-local clusters and confirm that VLS only affects the very luminous galaxies, with brighter galaxies having lower velocities. The threshold magnitude of VLS is ~m3+0.5, where m3 is the magnitude of the third brightest cluster galaxy, and current data suggest that the threshold value moves to fainter magnitudes at higher redshift. We also detect (marginal) evidence of VLS for blue galaxies. We conclude that the segregation effects, when their study is extended to distant clusters, can be important tracers of the galaxy evolution and cluster assembly and discuss the poor/no evidence of VCS at high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An improved SPH scheme for cosmological simulations

We present an implementation of smoothed particle hydrodynamics (SPH) with improved accuracy for simulations of galaxies and the large-scale structure. In particular, we combine, implement, modify and test a vast majority of SPH improvement techniques in the latest instalment of the GADGET code. We use the Wendland kernel functions, a particle wake-up time-step limiting mechanism and a time-dependent scheme for artificial viscosity, which includes a high-order gradient computation and shear flow limiter. Additionally, we include a novel prescription for time-dependent artificial conduction, which corrects for gravitationally induced pressure gradients and largely improves the SPH performance in capturing the development of gas-dynamical instabilities. We extensively test our new implementation in a wide range of hydrodynamical standard tests including weak and strong shocks as well as shear flows, turbulent spectra, gas mixing, hydrostatic equilibria and self-gravitating gas clouds. We jointly employ all modifications; however, when necessary we study the performance of individual code modules. We approximate hydrodynamical states more accurately and with significantly less noise than standard SPH. Furthermore, the new implementation promotes the mixing of entropy between different fluid phases, also within cosmological simulations. Finally, we study the performance of the hydrodynamical solver in the context of radiative galaxy formation and non-radiative galaxy cluster formation. We find galactic disks to be colder, thinner and more extended and our results on galaxy clusters show entropy cores instead of steadily declining entropy profiles. In summary, we demonstrate that our improved SPH implementation overcomes most of the undesirable limitations of standard SPH, thus becoming the core of an efficient code for large cosmological simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

CLASH-VLT: A Highly Precise Strong Lensing Model of the Galaxy Cluster RXC J2248.7-4431 (Abell S1063) and Prospects for Cosmography

We perform a comprehensive study of the total mass distribution of the galaxy cluster RXCJ2248 ($z=0.348$) with a set of high-precision strong lensing models, which take advantage of extensive spectroscopic information on many multiply lensed systems. In the effort to understand and quantify inherent systematics in parametric strong lensing modelling, we explore a collection of 22 models where we use different samples of multiple image families, parametrizations of the mass distribution and cosmological parameters. As input information for the strong lensing models, we use the CLASH HST imaging data and spectroscopic follow-up observations, carried out with the VIMOS and MUSE spectrographs, to identify bona-fide multiple images. A total of 16 background sources, over the redshift range $1.0-6.1$, are multiply lensed into 47 images, 24 of which are spectroscopically confirmed and belong to 10 individual sources. The cluster total mass distribution and underlying cosmology in the models are optimized by matching the observed positions of the multiple images on the lens plane. We show that with a careful selection of a sample of spectroscopically confirmed multiple images, the best-fit model reproduces their observed positions with a rms of $0.3$ in a fixed flat $Λ$CDM cosmology, whereas the lack of spectroscopic information lead to biases in the values of the model parameters. Allowing cosmological parameters to vary together with the cluster parameters, we find (at $68\%$ confidence level) $Ω_m=0.25^{+0.13}_{-0.16}$ and $w=-1.07^{+0.16}_{-0.42}$ for a flat $Λ$CDM model, and $Ω_m=0.31^{+0.12}_{-0.13}$ and $Ω_Λ=0.38^{+0.38}_{-0.27}$ for a universe with $w=-1$ and free curvature. Using toy models mimicking the overall configuration of RXCJ2248, we estimate the impact of the line of sight mass structure on the positional rms to be $0.3\pm 0.1$.(ABRIDGED)

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cool Core Clusters from Cosmological Simulations

We present results obtained from a set of cosmological hydrodynamic simulations of galaxy clusters, aimed at comparing predictions with observational data on the diversity between cool-core (CC) and non-cool-core (NCC) clusters. Our simulations include the effects of stellar and AGN feedback and are based on an improved version of the smoothed particle hydrodynamics code GADGET-3, which ameliorates gas mixing and better captures gas-dynamical instabilities by including a suitable artificial thermal diffusion. In this Letter, we focus our analysis on the entropy profiles, the primary diagnostic we used to classify the degree of cool-coreness of clusters, and on the iron profiles. In keeping with observations, our simulated clusters display a variety of behaviors in entropy profiles: they range from steadily decreasing profiles at small radii, characteristic of cool-core systems, to nearly flat core isentropic profiles, characteristic of non-cool-core systems. Using observational criteria to distinguish between the two classes of objects, we find that they occur in similar proportions in both simulations and in observations. Furthermore, we also find that simulated cool-core clusters have profiles of iron abundance that are steeper than those of NCC clusters, which is also in agreement with observational results. We show that the capability of our simulations to generate a realistic cool-core structure in the cluster population is due to AGN feedback and artificial thermal diffusion: their combined action allows us to naturally distribute the energy extracted from super-massive black holes and to compensate for the radiative losses of low-entropy gas with short cooling time residing in the cluster core.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Next Generation Cosmology: Constraints from the Euclid Galaxy Cluster Survey

We study the characteristics of the galaxy cluster samples expected from the European Space Agency's Euclid satellite and forecast constraints on cosmological parameters describing a variety of cosmological models. The method used in this paper, based on the Fisher Matrix approach, is the same one used to provide the constraints presented in the Euclid Red Book (Laureijs et al.2011). We describe the analytical approach to compute the selection function of the photometric and spectroscopic cluster surveys. Based on the photometric selection function, we forecast the constraints on a number of cosmological parameter sets corresponding to different extensions of the standard LambdaCDM model. The dynamical evolution of dark energy will be constrained to Delta w_0=0.03 and Delta w_a=0.2 with free curvature Omega_k, resulting in a (w_0,w_a) Figure of Merit (FoM) of 291. Including the Planck CMB covariance matrix improves the constraints to Delta w_0=0.02, Delta w_a=0.07 and a FoM=802. The amplitude of primordial non-Gaussianity, parametrised by f_NL, will be constrained to Δf_NL ~ 6.6 for the local shape scenario, from Euclid clusters alone. Using only Euclid clusters, the growth factor parameter γ, which signals deviations from GR, will be constrained to Delta γ=0.02, and the neutrino density parameter to Delta Omega_ν=0.0013 (or Delta \sum m_ν=0.01). We emphasise that knowledge of the observable--mass scaling relation will be crucial to constrain cosmological parameters from a cluster catalogue. The Euclid mission will have a clear advantage in this respect, thanks to its imaging and spectroscopic capabilities that will enable internal mass calibration from weak lensing and the dynamics of cluster galaxies. This information will be further complemented by wide-area multi-wavelength external cluster surveys that will already be available when Euclid flies. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Semi-Analytic Model Predictions of the Galaxy Population in Proto-clusters

We investigate the galaxy population in simulated proto-cluster regions using a semi-analytic model of galaxy formation, coupled to merger trees extracted from N-body simulations. We select the most massive clusters at redshift $z=0$ from our set of simulations, and follow their main progenitors back in time. The analysis shows that proto-cluster regions are dominated by central galaxies and their number decreases with time as many become satellites, clustering around the central object. In agreement with observations, we find an increasing velocity dispersion with cosmic time, the increase being faster for satellites. The analysis shows that proto-clusters are very extended regions, $\gtrsim 20 \, Mpc$ at $z \gtrsim 1$. The fraction of galaxies in proto-cluster regions that are not progenitor of cluster galaxies varies with redshift, stellar mass and area considered. It is about 20-30 per cent for galaxies with stellar mass $\sim 10^9\,{\rm M}_{\sun}$, while negligible for the most massive galaxies considered. Nevertheless, these objects have properties similar to those of progenitors. We investigate the building-up of the passive-sequence in clusters, and find that their progenitors are on average always active at any redshift of interest of proto-clusters. The main mechanism which quenches their star formation is the removal of the hot gas reservoir at the time of accretion. The later galaxies are accreted (become satellite), and the more the cold gas available, the longer the time spent as active. Central galaxies are active over all redshift range considered, although a non-negligible fraction of them become passive at redshift $z<1$, due to strong feedback from Active Galactic Nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Chandra deep observation of XDCP J0044.0-2033, a massive galaxy cluster at z>1.5

We report the analysis of the Chandra observation of XDCP J0044.0-2033, a massive, distant (z=1.579) galaxy cluster discovered in the XDCP survey. The total exposure time of 380 ks with Chandra ACIS-S provides the deepest X-ray observation currently achieved on a massive, high redshift cluster. Extended emission from the Intra Cluster Medium (ICM) is detected at a very high significance level (S/N~20) on a circular region with a 44" radius, corresponding to $R_{ext}=375$ kpc at the cluster redshift. We perform an X-ray spectral fit of the ICM emission modeling the spectrum with a single-temperature thermal mekal model. Our analysis provides a global temperature $kT=6.7^{+1.3}_{-0.9}$ keV, and a iron abundance $Z_{Fe} = 0.41_{-0.26}^{+0.29}Z_{Fe_\odot}$ (error bars correspond to 1 $σ$). We fit the background-subtracted surface brightness profile with a single $β$-model out to 44", finding a rather flat profile with no hints of a cool core. We derive the deprojected electron density profile and compute the ICM mass within the extraction radius $R_{ext}=375$ kpc to be $M_{ICM}(r<R_{ext}) = (1.48 \pm 0.20) \times 10^{13} M_\odot$. Under the assumption of hydrostatic equilibrium and assuming isothermality within $R_{ext}$, the total mass is $M_{2500}= 1.23_{-0.27}^{+0.46} \times 10 ^{14} M_\odot$ for $R_{2500} = 240_{-20}^{+30}$ kpc. Extrapolating the profile at radii larger than the extraction radius $R_{ext}$ we find $M_{500} = 3.2_{-0.6}^{+0.9} \times 10 ^{14}M_\odot$ for $R_{500} = 562_{-37}^{+50}$ kpc. This analysis establishes the existence of virialized, massive galaxy clusters at redshift $z\sim 1.6$, paving the way to the investigation of the progenitors of the most massive clusters today. Given its mass and the XDCP survey volume, XDCP J0044.0-2033 does not create significant tension with the WMAP-7 $Λ$CDM cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

CLASH-VLT: Insights on the mass substructures in the Frontier Fields Cluster MACS J0416.1-2403 through accurate strong lens modeling

We present a detailed mass reconstruction and a novel study on the substructure properties in the core of the CLASH and Frontier Fields galaxy cluster MACS J0416.1-2403. We show and employ our extensive spectroscopic data set taken with the VIMOS instrument as part of our CLASH-VLT program, to confirm spectroscopically 10 strong lensing systems and to select a sample of 175 plausible cluster members to a limiting stellar mass of log(M_*/M_Sun) ~ 8.6. We reproduce the measured positions of 30 multiple images with a remarkable median offset of only 0.3" by means of a comprehensive strong lensing model comprised of 2 cluster dark-matter halos, represented by cored elliptical pseudo-isothermal mass distributions, and the cluster member components. The latter have total mass-to-light ratios increasing with the galaxy HST/WFC3 near-IR (F160W) luminosities. The measurement of the total enclosed mass within the Einstein radius is accurate to ~5%, including systematic uncertainties. We emphasize that the use of multiple-image systems with spectroscopic redshifts and knowledge of cluster membership based on extensive spectroscopic information is key to constructing robust high-resolution mass maps. We also produce magnification maps over the central area that is covered with HST observations. We investigate the galaxy contribution, both in terms of total and stellar mass, to the total mass budget of the cluster. When compared with the outcomes of cosmological $N$-body simulations, our results point to a lack of massive subhalos in the inner regions of simulated clusters with total masses similar to that of MACS J0416.1-2403. Our findings of the location and shape of the cluster dark-matter halo density profiles and on the cluster substructures provide intriguing tests of the assumed collisionless, cold nature of dark matter and of the role played by baryons in the process of structure formation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Neutrino constraints: what large-scale structure and CMB data are telling us?

(Abridged) We discuss the reliability of neutrino mass constraints, either active or sterile, from the combination of WMAP 9-year or Planck CMB data with BAO measurements from BOSS DR11, galaxy shear measurements from CFHTLenS, SDSS Ly-$α$ forest constraints and galaxy cluster mass function from Chandra observations. To avoid model dependence of the constraints we perform a full likelihood analysis for all the datasets employed. As for the cluster data analysis we rely on to the most recent calibration of massive neutrino effects in the halo mass function and we explore the impact of the uncertainty in the mass bias and re-calibration of the halo mass function due to baryonic feedback processes on cosmological parameters. We find that none of the low redshift probes alone provide evidence for massive neutrinos in combination with CMB measurements, while a larger than $2σ$ detection of non zero neutrino mass, either active or sterile, is achieved combining cluster or shear data with CMB and BAO measurements. The preference for massive neutrino is larger in the sterile neutrino scenario, and for the combination of Planck, BAO, shear and cluster datasets we find that the vanilla $Λ$CDM model is rejected at more than $3σ$ and a sterile neutrino mass as motivated by accelerator anomaly is within the $2σ$ errors. Finally, results from the full data combination reflect the tension between the $σ_8$ constraints obtained from cluster and shear data and that inferred from Ly-$α$ forest measurements; in the active neutrino scenario for both CMB datasets employed, the full data combination yields only an upper limits on $\sum m_ν$, while assuming an extra sterile neutrino we still get preference for non-vanishing mass, $m_s^{\rm eff}=0.26^{+0.22}_{-0.24}$ eV, and dark contribution to the radiation content, $ΔN_{\rm eff}=0.82\pm0.55$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Swift X-ray Telescope Cluster Survey II. X-ray spectral analysis

(Abridged) We present a spectral analysis of a new, flux-limited sample of 72 X-ray selected clusters of galaxies identified with the X-ray Telescope (XRT) on board the Swift satellite down to a flux limit of ~10-14 erg/s/cm2 (SWXCS, Tundo et al. 2012). We carry out a detailed X-ray spectral analysis with the twofold aim of measuring redshifts and characterizing the properties of the Intra-Cluster Medium (ICM). Optical counterparts and spectroscopic or photometric redshifts are obtained with a cross-correlation with NED. Additional photometric redshifts are computed with a dedicated follow-up program with the TNG and a cross-correlation with the SDSS. We also detect the iron emission lines in 35% of the sample, and hence obtain a robust measure of the X-ray redshift zX. We use zX whenever the optical redshift is not available. Finally, for all the sources with measured redshift, background-subtracted spectra are fitted with a mekal model. We perform extensive spectral simulations to derive an empirical formula to account for fitting bias. The bias-corrected values are then used to investigate the scaling properties of the X-ray observables. Overall, we are able to characterize the ICM of 46 sources. The sample is mostly constituted by clusters with temperatures between 3 and 10 keV, plus 14 low-mass clusters and groups with temperatures below 3 keV. The redshift distribution peaks around z~0.25 and extends up to z~1, with 60% of the sample at 0.1<z<0.4. We derive the Luminosity-Temperature relation for these 46 sources, finding good agreement with previous studies. The quality of the SWXCS sample is comparable to other samples available in the literature and obtained with much larger X-ray telescopes. Our results have interesting implications for the design of future X-ray survey telescopes, characterised by good-quality PSF over the entire field of view and low background.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An accurate tool for the fast generation of dark matter halo catalogs

We present a new parallel implementation of the PINpointing Orbit Crossing-Collapsed HIerarchical Objects (PINOCCHIO) algorithm, a quick tool, based on Lagrangian Perturbation Theory, for the hierarchical build-up of Dark Matter halos in cosmological volumes. To assess its ability to predict halo correlations on large scales, we compare its results with those of an N-body simulation of a 3 Gpc/h box sampled with 2048^3 particles taken from the MICE suite, matching the same seeds for the initial conditions. Thanks to the FFTW libraries and to the relatively simple design, the code shows very good scaling properties. The CPU time required by PINOCCHIO is a tiny fraction (~1/2000) of that required by the MICE simulation. Varying some of PINOCCHIO numerical parameters allows one to produce a universal mass function that lies in the range allowed by published fits, although it underestimates the MICE mass function of FoF halos in the high mass tail. We compare the matter-halo and the halo-halo power spectra with those of the MICE simulation and find that these 2-point statistics are well recovered on large scales. In particular, when catalogs are matched in number density, agreement within ten per cent is achieved for the halo power spectrum. At scales k>0.1 h/Mpc, the inaccuracy of the Zel'dovich approximation in locating halo positions causes an underestimate of the power spectrum that can be modeled as a Gaussian factor with a damping scale of d=3 Mpc/h at z=0, decreasing at higher redshift. Finally, a remarkable match is obtained for the reduced halo bispectrum, showing a good description of nonlinear halo bias. Our results demonstrate the potential of PINOCCHIO as an accurate and flexible tool for generating large ensembles of mock galaxy surveys, with interesting applications for the analysis of large galaxy redshift surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Chandra View of the Warm-Hot IGM toward 1ES 1553+113: Absorption Line Detections and Identifications (Paper I)

We present the first results from our pilot 500 ks Chandra-LETG Large Program observation of the soft X-ray brightest source in the z>=0.4 sky, the blazar 1ES 1553+113, aimed to secure the first uncontroversial detections of the missing baryons in the X-rays. We identify a total of 11 possible absorption lines, with single-line statistical significances between 2.2-4.1σ. Six of these lines are detected at high single-line statistical significance (3.6 <= sigma <= 4.1), while the remaining five are regarded as marginal detections in association with either other X-ray lines detected at higher significance and/or Far-Ultraviolet (FUV) signposts. In particular, five of these possible intervening absorption lines, are identified as CV and CVI Kαabsorbers belonging to three WHIM systems at z_X = 0.312, z_X = 0.237 and <z_X> = 0.133, which also produce broad HI (and OVI for the z_X = 0.312 system) absorption in the FUV. For two of these systems (z_X = 0.312 and 0.237), the Chandra X-ray data led the a-posteriori discovery of physically consistent broad HI associations in the FUV, so confirming the power of the X-ray-FUV synergy for WHIM studies. The true statistical significances of these three X-ray absorption systems, after properly accounting for the number of redshift trials, are 5.8 sigma (z_X = 0.312; 6.3 sigma if the low-significance OV and CV K-beta associations are considered), 3.9 sigma (z_X = 0.237), and 3.8 sigma (\langle z_X \rangle = 0.133), respectively.

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

Constraining neutrino properties with a Euclid-like galaxy cluster survey

We perform a forecast analysis on how well a Euclid-like photometric galaxy cluster survey will constrain the total neutrino mass and effective number of neutrino species. We base our analysis on the Monte Carlo Markov Chains technique by combining information from cluster number counts and cluster power spectrum. We find that combining cluster data with CMB measurements from Planck improves by more than an order of magnitude the constraint on neutrino masses compared to each probe used independently. For the LCDM+m_nu model the 2 sigma upper limit on total neutrino mass shifts from M_nu < 0.35 eV using cluster data alone to M_nu < 0.031 eV when combined with CMB data. When a non-standard model with N_eff number of neutrino species is considered, we estimate N_eff<3.14 (95% CL), while the bounds on neutrino mass are relaxed to M_nu < 0.040 eV. This accuracy would be sufficient for a 2 sigma detection of neutrino mass even in the minimal normal hierarchy scenario. We also consider scenarios with a constant dark energy equation of state and a non-vanishing curvature. When these models are considered the error on M_nu is only slightly affected, while there is a larger impact of the order of ~ 15 % and ~ 20% respectively on the 2 sigma error bar of N_eff with respect to the standard case. We also treat the LCDM+m_nu+N_eff case with free nuisance parameters, which parameterize the uncertainties on the cluster mass determination. In this case, the upper bounds on M_nu are relaxed by a factor larger than two, M_nu < 0.083 eV (95% CL), hence compromising the possibility of detecting the total neutrino mass with good significance. We thus confirm the potential that a large optical/near-IR cluster survey, like that to be carried out by Euclid, could have in constraining neutrino properties [abridged].

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Large-scale inhomogeneities of the intracluster medium: improving mass estimates using the observed azimuthal scatter

Using a set of hydrodynamical simulations of 62 galaxy clusters and groups we study the ICM of inhomogeneities, focusing on the ones on the large scale that, unlike clumps, are the most difficult to identify. To this purpose we introduce the concept of residual clumpiness, C_R, that quantifies the large-scale inhomogeneity of the ICM. After showing that this quantity can be robustly defined for relaxed systems, we characterize how it varies with radius, mass and dynamical state of the halo. Most importantly, we observe that it introduces an overestimate in the determination of the density profile from the X-ray emission, which translates into a systematic overestimate of 6 (12)% in the measurement of M_gas at R_200 for our relaxed (perturbed) cluster sample. At the same time, the increase of C_R with radius introduces also a ~2% systematic underestimate in the measurement of the hydrostatic-equilibrium mass (M_he), which adds to the previous one generating a systematic ~8.5% overestimate in f_gas in our relaxed sample. Since the residual clumpiness of the ICM is not directly observable, we study its correlation with the azimuthal scatter in the X-ray surface brightness of the halo and in the y-parameter profiles. We find that their correlation is highly significant (r_S = 0.6-0.7), allowing to define the azimuthal scatter measured in the X-ray surface brightness profile and in the y-parameter as robust proxies of C_R. After providing a function that connects the two quantities, we obtain that correcting the observed gas density profiles using the azimuthal scatter eliminates the bias in the measurement of M_gas for relaxed objects, which becomes (0+/-2)% up to 2R_200, and reduces it by a factor of 3 for perturbed ones. This method allows also to eliminate the systematics on the measurements of M_he and f_gas, although a significant halo to halo scatter remains. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the formation and physical properties of the Intra-Cluster Light in hierarchical galaxy formation models

We study the formation of the Intra-Cluster Light (ICL) using a semi-analytic model of galaxy formation, coupled to merger trees extracted from N-body simulations of groups and clusters. We assume that the ICL forms by (1) stellar stripping of satellite galaxies and (2) relaxation processes that take place during galaxy mergers. The fraction of ICL in groups and clusters predicted by our models ranges between 10 and 40 per cent, with a large halo-to-halo scatter and no halo mass dependence. We note, however, that our predicted ICL fractions depend on the resolution: for a set of simulations with particle mass one order of magnitude larger than that adopted in the high resolution runs used in our study, we find that the predicted ICL fractions are ~30-40 per cent larger than those found in the high resolution runs. On cluster scale, large part of the scatter is due to a range of dynamical histories, while on smaller scale it is driven by individual accretion events and stripping of very massive satellites, $M_{*} \gtrsim 10^{10.5} M_{\odot}$, that we find to be the major contributors to the ICL. The ICL in our models forms very late (below $z\sim 1$), and a fraction varying between 5 and 25 per cent of it has been accreted during the hierarchical growth of haloes. In agreement with recent observational measurements, we find the ICL to be made of stars covering a relatively large range of metallicity, with the bulk of them being sub-solar.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the role of AGN feedback on the thermal and chemodynamical properties of the hot intra-cluster medium

We present an analysis of the properties of the ICM in an extended set of cosmological hydrodynamical simulations of galaxy clusters and groups performed with the TreePM+SPH GADGET-3 code. Besides a set of non-radiative simulations, we carried out two sets of simulations including radiative cooling, star formation, metal enrichment and feedback from supernovae, one of which also accounts for the effect of feedback from AGN resulting from gas accretion onto super-massive black holes. These simulations are analysed with the aim of studying the relative role played by SN and AGN feedback on the general properties of the diffuse hot baryons in galaxy clusters and groups: scaling relations, temperature, entropy and pressure radial profiles, and ICM chemical enrichment. We find that simulations including AGN feedback produce scaling relations that are in good agreement with X-ray observations at all mass scales. However, our simulations are not able to account for the observed diversity between CC and NCC clusters: unlike for observations, we find that temperature and entropy profiles of relaxed and unrelaxed clusters are quite similar and resemble more the observed behaviour of NCC clusters. As for the pattern of metal enrichment, we find that an enhanced level of iron abundance is produced by AGN feedback with respect to the case of purely SN feedback. As a result, while simulations including AGN produce values of iron abundance in groups in agreement with observations, they over-enrich the ICM in massive clusters. The efficiency of AGN feedback in displacing enriched gas from halos into the inter-galactic medium at high redshift also creates a widespread enrichment in the outskirts of clusters and produces profiles of iron abundance whose slope is in better agreement with observations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Searching for large-scale structures around high-redshift radio galaxies with Herschel

This paper presents the first results of a far-infrared search for protocluster-associated galaxy overdensities using the SPIRE instrument on-board the {\it Herschel} Space Observatory. Large ($\sim$400 arcmin$^{2}$) fields surrounding 26 powerful high-redshift radio galaxies ($2.0 < z < 4.1$; $L_{\rm 500 MHz} > 10^{28.5}$ WHz$^{-1}$) are mapped at 250, 350 and 500\mic to give a unique wide-field sample. On average the fields have a higher than expected, compared to blank fields, surface density of 500\mic sources within 6 comoving Mpc of the radio galaxy. The analysis is then restricted to potential protocluster members only, which are identified using a far-infrared colour selection; this reveals significant overdensities of galaxies in 2 fields, neither of which are previously known protoclusters. The probability of finding 2 overdensities of this size by chance, given the number of fields observed is $5 \times 10^{-4}$. Overdensities here exist around radio galaxies with $L_{\rm 500 MHz} \gtrsim 10^{29}$ WHz$^{-1}$ and $z < 3$. The radial extent of the average far-infrared overdensity is found to be $\sim$6 comoving Mpc. Comparison with predictions from numerical simulations shows that the overdensities are consistent with having masses $> 10^{14}$Msolar. However, the large uncertainty in the redshift estimation means that it is possible that these far-infrared overdensities consist of several structures across the redshift range searched.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Hot and Energetic Universe: The astrophysics of galaxy groups and clusters

As the nodes of the cosmic web, clusters of galaxies trace the large-scale distribution of matter in the Universe. They are thus privileged sites in which to investigate the complex physics of structure formation. However, the complete story of how these structures grow, and how they dissipate the gravitational and non-thermal components of their energy budget over cosmic time, is still beyond our grasp. Fundamental questions such as How do hot diffuse baryons accrete and dynamically evolve in dark matter potentials? How and when was the energy that we observe in the ICM generated and distributed? Where and when are heavy elements produced and how are they circulated? are still unanswered. Most of the cluster baryons exists in the form of a diffuse, hot, metal-enriched plasma that radiates primarily in the X-ray band (the intracluster medium, ICM), allowing the X-ray observations of the evolving cluster population to provide a unique opportunity to address these topics. Athena+ with its large collecting area and unprecedented combination of high spectral and angular resolution offers the only way to make major advances in answering these questions. Athena+ will show how the baryonic gas evolves in the dark matter potential wells by studying the motions and turbulence in the ICM. Athena+ will be able to resolve the accreting region both spatially and spectroscopically, probing the true nature and physical state of the X-ray emitting plasma. Athena+ has the capabilities to permit a definitive understanding of the formation and evolution of large-scale cosmic structure through the study of the cluster population.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Hot and Energetic Universe: The evolution of galaxy groups and clusters

Major astrophysical questions related to the formation and evolution of structures, and more specifically of galaxy groups and clusters, will still be open in the coming decade and beyond: what is the interplay of galaxy, supermassive black hole, and intergalactic gas evolution in the most massive objects in the Universe - galaxy groups and clusters? What are the processes driving the evolution of chemical enrichment of the hot diffuse gas in large-scale structures? How and when did the first galaxy groups in the Universe, massive enough to bind more than 10^7 K gas, form? Focussing on the period when groups and clusters assembled (0.5<z<2.5), we show that, due to the continuum and line emission of this hot intergalactic gas at X-ray wavelengths, Athena+, combining high sensitivity with excellent spectral and spatial resolution, will deliver breakthrough observations in view of the aforementioned issues. Indeed, the physical and chemical properties of the hot intra-cluster gas, and their evolution across time, are a key to understand the co-evolution of galaxy and supermassive black hole within their environments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Lensing and X-ray mass estimates of clusters (SIMULATION)

[Abridged] We present a comparison between weak-lensing (WL) and X-ray mass estimates of a sample of numerically simulated clusters. The sample consists on the 20 most massive objects at redshift z=0.25 and Mvir > 5 x 10^{14} Msun h^{-1}. They were found in a cosmological simulation of volume 1 h^{-3} Gpc^3, evolved in the framework of a WMAP-7 normalized cosmology. Each cluster has been resimulated at higher resolution and with more complex gas physics. We processed it thought Skylens and X-MAS to generate optical and X-ray mock observations along three orthogonal projections. The optical simulations include lensing effects on background sources. Standard observational tools and methods of analysis are used to recover the mass profiles of each cluster projection from the mock catalogues. Given the size of our sample, we could also investigate the dependence of the results on cluster morphology, environment, temperature inhomogeneity, and mass. We confirm previous results showing that WL masses obtained from the fit of the cluster tangential shear profiles with NFW functionals are biased low by ~ 5-10% with a large scatter (~10-25%). We show that scatter could be reduced by optimally selecting clusters either having regular morphology or living in substructure-poor environment. The X-ray masses are biased low by a large amount (~25-35%), evidencing the presence of both non-thermal sources of pressure in the ICM and temperature inhomogeneity, but they show a significantly lower scatter than weak-lensing-derived masses. The X-ray mass bias grows from the inner to the outer regions of the clusters. We find that both biases are weakly correlated with the third-order power ratio, while a stronger correlation exists with the centroid shift. Finally, the X-ray bias is strongly connected with temperature inhomogeneities.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Planck intermediate results. III. The relation between galaxy cluster mass and Sunyaev-Zeldovich signal

We examine the relation between the galaxy cluster mass M and Sunyaev-Zeldovich (SZ) effect signal D_A^2 Y for a sample of 19 objects for which weak lensing (WL) mass measurements obtained from Subaru Telescope data are available in the literature. Hydrostatic X-ray masses are derived from XMM-Newton archive data and the SZ effect signal is measured from Planck all-sky survey data. We find an M_WL-D_A^2 Y relation that is consistent in slope and normalisation with previous determinations using weak lensing masses; however, there is a normalisation offset with respect to previous measures based on hydrostatic X-ray mass-proxy relations. We verify that our SZ effect measurements are in excellent agreement with previous determinations from Planck data. For the present sample, the hydrostatic X-ray masses at R_500 are on average ~ 20 per cent larger than the corresponding weak lensing masses, at odds with expectations. We show that the mass discrepancy is driven by a difference in mass concentration as measured by the two methods, and, for the present sample, the mass discrepancy and difference in mass concentration is especially large for disturbed systems. The mass discrepancy is also linked to the offset in centres used by the X-ray and weak lensing analyses, which again is most important in disturbed systems. We outline several approaches that are needed to help achieve convergence in cluster mass measurement with X-ray and weak lensing observations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Planck Intermediate Results. IV. The XMM-Newton validation programme for new Planck galaxy clusters

We present the final results from the XMM-Newton validation follow-up of new Planck galaxy cluster candidates. We observed 15 new candidates, detected with signal-to-noise ratios between 4.0 and 6.1 in the 15.5-month nominal Planck survey. The candidates were selected using ancillary data flags derived from the ROSAT All Sky Survey (RASS) and Digitized Sky Survey all-sky maps, with the aim of pushing into the low SZ flux, high-z regime and testing RASS flags as indicators of candidate reliability. 14 new clusters were detected by XMM, including 2 double systems. Redshifts lie in the range 0.2 to 0.9, with 6 clusters at z>0.5. Estimated M500 range from 2.5 10^14 to 8 10^14 Msun. We discuss our results in the context of the full XMM validation programme, in which 51 new clusters have been detected. This includes 4 double and 2 triple systems, some of which are chance projections on the sky of clusters at different z. We find that association with a RASS-BSC source is a robust indicator of the reliability of a candidate, whereas association with a FSC source does not guarantee that the SZ candidate is a bona fide cluster. Nevertheless, most Planck clusters appear in RASS maps, with a significance greater than 2 sigma being a good indication that the candidate is a real cluster. The full sample gives a Planck sensitivity threshold of Y500 ~ 4 10^-4 arcmin^2, with indication for Malmquist bias in the YX-Y500 relation below this level. The corresponding mass threshold depends on z. Systems with M500 > 5 10^14 Msun at z > 0.5 are easily detectable with Planck. The newly-detected clusters follow the YX-Y500 relation derived from X-ray selected samples. Compared to X-ray selected clusters, the new SZ clusters have a lower X-ray luminosity on average for their mass. There is no indication of departure from standard self-similar evolution in the X-ray versus SZ scaling properties. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Planck Intermediate Results. V. Pressure profiles of galaxy clusters from the Sunyaev-Zeldovich effect

Taking advantage of the all-sky coverage and broad frequency range of the Planck satellite, we study the Sunyaev-Zeldovich (SZ) and pressure profiles of 62 nearby massive clusters detected at high significance in the 14-month nominal survey. Careful reconstruction of the SZ signal indicates that most clusters are individually detected at least out to R500. By stacking the radial profiles, we have statistically detected the radial SZ signal out to 3 x R500, i.e., at a density contrast of about 50-100, though the dispersion about the mean profile dominates the statistical errors across the whole radial range. Our measurement is fully consistent with previous Planck results on integrated SZ fluxes, further strengthening the agreement between SZ and X-ray measurements inside R500. Correcting for the effects of the Planck beam, we have calculated the corresponding pressure profiles. This new constraint from SZ measurements is consistent with the X-ray constraints from XMM-Newton in the region in which the profiles overlap (i.e., [0.1-1]R500), and is in fairly good agreement with theoretical predictions within the expected dispersion. At larger radii the average pressure profile is slightly flatter than most predictions from numerical simulations. Combining the SZ and X-ray observed profiles into a joint fit to a generalised pressure profile gives best-fit parameters [P0, c500, gamma, alpha, beta] = [6.41, 1.81, 0.31, 1.33, 4.13]. Using a reasonable hypothesis for the gas temperature in the cluster outskirts we reconstruct from our stacked pressure profile the gas mass fraction profile out to 3 x R500. Within the temperature driven uncertainties, our Planck constraints are compatible with the cosmic baryon fraction and expected gas fraction in halos.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Probing dark energy with the next generation X-ray surveys of galaxy clusters

We present forecasts on the capability of future wide-area high-sensitivity X-ray surveys of galaxy clusters to yield constraints on the parameters defining the Dark Energy (DE) equation of state (EoS). Our analysis is carried out for future X-ray surveys which have enough sensitivity to provide accurate measurements of X-ray mass proxies and Fe-line based redshifts for about 2x10^4 clusters. We base our analysis on the Fisher Matrix formalism, by combining information on the cluster number counts and power spectrum, also including, for the first time in the analysis of the large scale cluster distribution, the effect of linear redshift-space distortions (RSDs). This study is performed with the main purpose of dissecting the cosmological information provided by geometrical and growth tests, which are both included in the analysis of number counts and clustering of galaxy clusters. We compare cosmological constraints obtained by assuming different levels of prior knowledge of the parameters which define the observable-mass X-ray relation. This comparison further demonstrates the fundamental importance of having a well calibrated observable-mass relation and, most importantly, its redshift evolution. Such a calibration can be achieved only by having at least $\sim 10^3$ net photon counts for each cluster included in the survey. We show that RSDs in the power spectrum analysis carry important cosmological information also when traced with galaxy clusters and the DE FoM increases by a factor of 8. Besides confirming the potential that large cluster surveys have in constraining the nature of DE, our analysis emphasizes that a full exploitation of the cosmological information carried by such surveys requires not only a large statistic but also a robust measurement of the mass proxies and redshifts for a significant fraction of the cluster sample, derived from the same X-ray survey data.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The effect of feedback on the emission properties of the Warm-Hot Intergalactic Medium

At present, 30-40 per cent of the baryons in the local Universe is still undetected. According to theoretical predictions, this gas should reside in filaments filling the large-scale structure (LSS) in the form of a Warm-Hot Intergalactic Medium (WHIM), at temperatures of 10^5 - 10^7 K, thus emitting in the soft X-ray energies via free-free interaction and line emission from heavy elements. In this work we characterize the properties of the X-ray emission of the WHIM, and the LSS in general, focusing on the influence of different physical mechanisms, namely galactic winds (GWs), black-hole feedback and star-formation, and providing estimates of possible observational constraints. To this purpose we use a set of cosmological hydrodynamical simulations that include a self-consistent treatment of star-formation and chemical enrichment of the intergalactic medium, that allows us to follow the evolution of different metal species. We construct a set of simulated light-cones to make predictions of the emission in the 0.3-10 keV energy range. We obtain that GWs increase by a factor of 2 the emission of both galaxy clusters and WHIM. The amount of oxygen at average temperature and, consequently, the amount of expected bright Ovii and Oviii lines is increased by a factor of 3 due to GWs and by 20 per cent when assuming a top-heavy IMF. We compare our results with current observational constraints and find that the emission from faint groups and WHIM should account from half to all of the unresolved X-ray background in the 1-2 keV band.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A non-ideal MHD Gadget: Simulating massive galaxy clusters

Magnetic fields in the intra-cluster medium of galaxy clusters have been studied in the past years through different methods. In the next years the up-coming generation of radio telescopes is going to provide new data that have the potential of setting constraints on the properties of magnetic fields in galaxy clusters. Here we present zoomed-in simulations for a set of massive galaxy clusters (M_v > 10^15 M_sun/h). This is an ideal sample to study the evolution of magnetic field during the process of structure formation in detail. Turbulent motions of the gas within the ICM will manifest themselves in a macroscopic magnetic resistivity eta_m, which has to be taken explicitly into account, especially at scales below the resolution limit. We have adapted the MHD GADGET code by Dolag & Stasyszyn (2009) to include the treatment of the magnetic resistivity and for the first time we have included non-ideal MHD equations to better follow the evolution of the magnetic field within galaxy clusters. We investigate which value of the magnetic resistivity eta_m is required to match the magnetic field profile derived from radio observations. We find that a value of eta_m ~ 6 10^27 cm^2/s is necessary to recover the shape of the magnetic field profile inferred from radio observations of the Coma cluster. This value agrees well with the expected level of turbulent motions within the ICM at our resolution limit. The magnetic field profiles can be fitted by a beta-model like profile (Cavaliere & Fusco-Femiano 1976), with small dispersion of the parameters. We find also that that the temperature, density and entropy profiles of the clusters depend on the magnetic resistivity constant, having flatter profiles in the inner regions when the magnetic resistivity increases.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Euclid Definition Study Report

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Fossil Groups Origins: I. RX J105453.3+552102 a very massive and relaxed system at z~0.5

The most accepted scenario for the origin of fossil groups (FGs) is that they are galaxy associations in which the merging rate was fast and efficient. These systems have assembled half of their mass at early epoch of the Universe, subsequently growing by minor mergers. They could contain a fossil record of the galaxy structure formation. We have started a project in order to characterize a large sample of FGs. In this paper we present the analysis of the fossil system RX J105453.3+552102. Optical deep images were used for studying the properties of the brightest group galaxy and for computing the photometric luminosity function of the group. We have also performed a detail dynamical analysis of the system based on redshift data for 116 galaxies. This galaxy system is located at z=0.47, and shows a quite large line-of-sight velocity dispersion σ_{v}~1000 km/s. Assuming the dynamical equilibrium, we estimated a virial mass of M ~ 10^{15} h_{70} M_{\odot}. No evidence of substructure was found within 1.4 Mpc radius. We found a statistically significant departure from Gaussianity of the group members velocities in the most external regions of the group. This could indicate the presence of galaxies in radial orbits in the external region of the group. We also found that the photometrical luminosity function is bimodal, showing a lack of M_{r} ~ -19.5 galaxies. The brightest group galaxy shows low Sersic parameter (n~2) and a small peculiar velocity. Indeed, our accurate photometry shows that the difference between the brightest and the second brightest galaxies is 1.9 mag in the r-band, while the classical definition of FGs is based on a magnitude gap of 2. We conclude that this fossil system does not follow the empirical definition of FGs. Nevertheless, it is a massive, old and undisturbed galaxy system with little infall of L^{*} galaxies since its initial collapse.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

ORIGIN: Metal Creation and Evolution from the Cosmic Dawn

ORIGIN is a proposal for the M3 mission call of ESA aimed at the study of metal creation from the epoch of cosmic dawn. Using high-spectral resolution in the soft X-ray band, ORIGIN will be able to identify the physical conditions of all abundant elements between C and Ni to red-shifts of z=10, and beyond. The mission will answer questions such as: When were the first metals created? How does the cosmic metal content evolve? Where do most of the metals reside in the Universe? What is the role of metals in structure formation and evolution? To reach out to the early Universe ORIGIN will use Gamma-Ray Bursts (GRBs) to study their local environments in their host galaxies. This requires the capability to slew the satellite in less than a minute to the GRB location. By studying the chemical composition and properties of clusters of galaxies we can extend the range of exploration to lower redshifts (z ~ 0.2). For this task we need a high-resolution spectral imaging instrument with a large field of view. Using the same instrument, we can also study the so far only partially detected baryons in the Warm-Hot Intergalactic Medium (WHIM). The less dense part of the WHIM will be studied using absorption lines at low redshift in the spectra for GRBs.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Planck Early Results XXVI: Detection with Planck and confirmation by XMM-Newton of PLCK G266.6-27.3, an exceptionally X-ray luminous and massive galaxy cluster at z~1

We present first results on PLCK G266.6-27.3, a galaxy cluster candidate detected at a signal-to-noise ratio of 5 in the Planck All Sky survey. An XMM-Newton validation observation has allowed us to confirm that the candidate is a bona fide galaxy cluster. With these X-ray data we measure an accurate redshift, z = 0.94 +/- 0.02, and estimate the cluster mass to be M_500 = (7.8 +/- 0.8)e+14 solar masses. PLCK G266.6-27.3 is an exceptional system: its luminosity of L_X(0.5-2.0 keV)=(1.4 +/- 0.05)e+45 erg/s, equals that of the two most luminous known clusters in the z > 0.5 universe, and it is one of the most massive clusters at z~1. Moreover, unlike the majority of high-redshift clusters, PLCK G266.6-27.3 appears to be highly relaxed. This observation confirms Planck's capability of detecting high-redshift, high-mass clusters, and opens the way to the systematic study of population evolution in the exponential tail of the mass function.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Pointing to the minimum scatter: the generalized scaling relations for galaxy clusters

We introduce a generalized scaling law, M_tot = 10^K A^a B^b, to look for the minimum scatter in reconstructing the total mass of hydrodynamically simulated X-ray galaxy clusters, given gas mass M_gas, luminosity L and temperature T. We find a locus in the plane of the logarithmic slopes $a$ and $b$ of the scaling relations where the scatter in mass is minimized. This locus corresponds to b_M = -3/2 a_M +3/2 and b_L = -2 a_L +3/2 for A=M_gas and L, respectively, and B=T. Along these axes, all the known scaling relations can be identified (at different levels of scatter), plus a new one defined as M_tot ~ (LT)^(1/2). Simple formula to evaluate the expected evolution with redshift in the self-similar scenario are provided. In this scenario, no evolution of the scaling relations is predicted for the cases (b_M=0, a_M=1) and (b_L=7/2, a_L=-1), respectively. Once the single quantities are normalized to the average values of the sample under considerations, the normalizations K corresponding to the region with minimum scatter are very close to zero. The combination of these relations allows to reduce the number of free parameters of the fitting function that relates X-ray observables to the total mass and includes the self-similar redshift evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Statistics of Substructures in Dark Matter Haloes

We study the amount and distribution of dark matter substructures within dark matter haloes, using a large set of high-resolution simulations ranging from group size to cluster size haloes, and carried our within a cosmological model consistent with WMAP 7-year data. In particular, we study how the measured properties of subhaloes vary as a function of the parent halo mass, the physical properties of the parent halo, and redshift. The fraction of halo mass in substructures increases with increasing mass. There is, however, a very large halo-to-halo scatter that can be explained only in part by a range of halo physical properties, e.g. concentration. At given halo mass, less concentrated haloes contain significantly larger fractions of mass in substructures because of the reduced strength of tidal disruption. Most of the substructure mass is located at the outskirts of the parent haloes, in relatively few massive subhaloes. This mass segregation appears to become stronger at increasing redshift, and should reflect into a more significant mass segregation of the galaxy population at different cosmic epochs. When haloes are accreted onto larger structures, their mass is significantly reduced by tidal stripping. Haloes that are more massive at the time of accretion (these should host more luminous galaxies) are brought closer to the centre on shorter time-scales by dynamical friction, and therefore suffer of a more significant stripping. The halo merger rate depends strongly on the environment with substructure in more massive haloes suffering more important mergers than their counterparts residing in less massive systems. This should translate into a different morphological mix for haloes of different mass.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

X-ray mass proxies from hydrodynamic simulations of galaxy clusters (paper I)

We present a detailed study of scaling relations between total cluster mass and three mass proxies based on X-ray observables: temperature of the intra-cluster medium, gas mass and the product of the two, Y_X. Our analysis is based on two sets of high-resolution hydrodynamical simulations performed with the TreePM-SPH GADGET code. The first set includes about 140 clusters with masses above 5x10^13 M_sun/h (30 having mass above 10^15 M_sun/h), that have been simulated with (i) non-radiative physics and including (ii) cooling, star formation, chemical enrichment and the effect of supernova feedback triggering galactic ejecta. This large statistics is used to quantify the robustness of the scaling relations, to determine their redshift evolution and to calibrate their intrinsic scatter and its distribution. We use a smaller set of clusters including 18 halos with masses above 5x10^13 M_sun/h to test the robustness of mass proxies against changing the physical processes included in simulations (thermal conduction, artificial viscosity, cooling and star formation, galactic winds and AGN feedback). We find the M-Y_X scaling relation to be the least sensitive one to variations of the ICM physics, with its slope and redshift evolution close to the self-similar model predictions. The distribution of the scatter around the best-fitting relations is close to a log-normal one. M_gas has the smallest scatter in mass, with values of sigma_lnM = 0.04-0.06, depending on the physics included in the simulation, and with a mild dependence on redshift. The M-T relation is the one with the largest scatter, with sigma_lnM > 0.1 at z=0, increasing to > 0.15 at z=1. The intrinsic scatter in the M-Y_X relation is slightly larger than in the M-M_gas relation. These results confirm that both Y_X and M_gas mass proxies are well suited for cosmological applications of future large X-ray surveys. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Astrophysics and cosmology with galaxy clusters: the WFXT perspective

We discuss the central role played by the X-ray study of hot baryons within galaxy clusters to reconstruct the assembly of cosmic structures and to trace the past history of star formation and accretion onto supermassive Black Holes (BHs). We shortly review the progress in this field contributed by the current generation of X-ray telescopes. Then, we focus on the outstanding scientific questions that have been opened by observations carried out in the last years and that represent the legacy of Chandra and XMM: (a) When and how is entropy injected into the inter-galactic medium (IGM)? (b) What is the history of metal enrichment of the IGM? (c) What physical mechanisms determine the presence of cool cores in galaxy clusters? (d) How is the appearance of proto-clusters at z~2 related to the peak of star formation activity and BH accretion? (e) What do galaxy clusters tell us about the nature of primordial density perturbations and on the history of their growth? We show that the most efficient observational strategy to address these questions is to carry out a large-area X-ray survey, reaching a sensitivity comparable to that of deep Chandra and XMM pointings, but extending over several thousands of square degrees. A similar survey can only be carried out with a Wide-Field X-ray Telescope (WFXT), which combines a high survey speed with a sharp PSF across the entire FoV. We emphasize the important synergies that WFXT will have with a number of future ground-based and space telescopes, covering from the radio to the X-ray bands. Finally, we discuss the immense legacy value that such a mission will have for extragalactic astronomy at large.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dynamical difference between the cD galaxy and the stellar diffuse component in simulated galaxy clusters

Member galaxies within galaxy clusters nowadays can be routinely identified in cosmological, hydrodynamical simulations using methods based on identifying self bound, locally over dense substructures. However, distinguishing the central galaxy from the stellar diffuse component within clusters is notoriously difficult, and in the center it is not even clear if two distinct stellar populations exist. Here, after subtracting all member galaxies, we use the velocity distribution of the remaining stars and detect two dynamically, well-distinct stellar components within simulated galaxy clusters. These differences in the dynamics can be used to apply an un-binding procedure which leads to a spatial separation of the two components into a cD and a diffuse stellar component (DSC). Applying our new algorithm to a cosmological, hydrodynamical simulation we find that -- in line with previous studies -- these two components have clearly distinguished spatial and velocity distributions as well as different star formation histories. We show that the DSC fraction -- which can broadly be associated with the observed intra cluster light -- does not depend on the virial mass of the galaxy cluster and is much more sensitive to the formation history of the cluster. We conclude that the separation of the cD and the DSC in simulations, based on our dynamical criteria, is more physically motivated than current methods which depend on implicit assumptions on a length scale associated with the cD galaxy and therefore represent a step forward in understanding the different stellar components within galaxy clusters. Our results also show the importance of analyzing the dynamics of the DSC to characterize its properties and understand its origin.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Gas cooling in semi-analytic models and SPH simulations: are results consistent?

We present a detailed comparison between the galaxy populations within a massive cluster, as predicted by hydrodynamical SPH simulations and by a semi-analytic model (SAM) of galaxy formation. Both models include gas cooling and a simple prescription of star formation, which consists in transforming instantaneously any cold gas available into stars, while neglecting any source of energy feedback. We find that, in general, galaxy populations from SAMs and SPH have similar statistical properties, in agreement with previous studies. However, when comparing galaxies on an object-by-object basis, we find a number of interesting differences: a) the star formation histories of the brightest cluster galaxies (BCGs) from SAM and SPH models differ significantly, with the SPH BCG exhibiting a lower level of star formation activity at low redshift, and a more intense and shorter initial burst of star formation with respect to its SAM counterpart; b) while all stars associated with the BCG were formed in its progenitors in the semi-analytic model used here, this holds true only for half of the final BCG stellar mass in the SPH simulation, the remaining half being contributed by tidal stripping of stars from the diffuse stellar component associated with galaxies accreted on the cluster halo; c) SPH satellites can loose up to 90 per cent of their stellar mass at the time of accretion, due to tidal stripping, a process not included in the semi-analytic model used in this study; d) in the SPH simulation, significant cooling occurs on the most massive satellite galaxies and this lasts for up to 1 Gyr after accretion. This physical process is not included in the semi-analytic model used in our study, as well as in most of the models discussed in the recent literature.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Simulating the WFXT sky

We investigate the scientific impact of the Wide Field X-ray Telescope mission. We present simulated images and spectra of X-ray sources as observed from the three surveys planned for the nominal 5-year WFXT lifetime. The goal of these simulations is to provide WFXT images of the extragalactic sky in different energy bands based on accurate description of AGN populations, normal and star forming galaxies, groups and clusters of galaxies. The images are realized using a detailed PSF model, instrumental and physical backgrounds/foregrounds, accurate model of the effective area and the related vignetting effect. Thanks to this comprehensive modelization of the WFXT properties, the simulated images can be used to evaluate the flux limits for detection of point and extended sources, the effect of source confusion at very faint fluxes, and in general the efficiency of detection algorithms. We also simulate the spectra of the detected sources, in order to address specific science topics which are unique to WFXT. Among them, we focus on the characterization of the Intra Cluster Medium (ICM) of high-z clusters, and in particular on the measurement of the redshift from the ICM spectrum in order to build a cosmological sample of galaxy clusters. The end-to-end simulation procedure presented here, is a valuable tool in optimizing the mission design. Therefore, these simulations can be used to reliably characterize the WFXT discovery space and to verify the connection between mission requirements and scientific goals. Thanks to this effort, we can conclude on firm basis that an X-ray mission optimized for surveys like WFXT is necessary to bring X-ray astronomy at the level of the optical, IR, submm and radio wavebands as foreseen in the coming decade.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The potential of X-ray cluster surveys to constrain primordial non-Gaussianity

We present forecasts for constraints on deviations from Gaussian distribution of primordial density perturbations from future X-ray surveys of galaxy clusters. Our analysis is based on computing the Fisher-Matrix for number counts and large-scale power spectrum of clusters. We consider a survey with high-sensitivity and wide-area to detect about 2.5 x 10^5 extended sources. Based on the self-calibration approach, and including Planck priors in our analysis, we constrain 9 cosmological parameters and 4 nuisance parameters, which define the relation between cluster mass and X-ray flux. Because of the scale dependence of large-scale bias induced by local-shape non-Gaussianity, we find that the power spectrum provides strong constraints on the non-Gaussianity f_nl parameter, which complement the stringent constraints on the power spectrum normalization, σ_8, from the number counts. To quantify the joint constraints on theese two parameters, that specify the timing of structure formation for a fixed background expansion, we define the figure-of-merit FoM_SFT = (det[Cov(sigma_8,f_nl)])^{-1/2}. We find that our surveys constrain deviations from Gaussianity with a precision of Δf_nl~10 at 1 σconfidence level, with FoM_SFT~39. We point out that constraints on f_nl are weakly sensitive to the uncertainties in the knowledge of nuisance parameters. As an application of non-Gaussian constraints from available data, we analyse the impact of positive skewness on the occurrence of XMMU-J2235, a massive distant cluster discovered at z~1.4. We confirm that in a WMAP-7 Gaussian Lambda CDM cosmology, within the survey volume, ~5x10^{-3} objects like this are expected to be found: to increase this probability by a factor of at least 10, one needs to evade either the available constraints on f_nl or sigma_8.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Wide Field X-ray Telescope: Mission Overview

The Wide Field X-Ray Telescope (WFXT) is a medium-class mission designed to be 2-orders-of-magnitude more sensitive than any previous or planned X-ray mission for large area surveys and to match in sensitivity the next generation of wide-area optical, IR and radio surveys. Using an innovative wide-field X-ray optics design, WFXT provides a field of view of 1 square degree (10 times Chandra) with an angular resolution of 5" (Half Energy Width, HEW) nearly constant over the entire field of view, and a large collecting area (up to 1 m^2 at 1 keV, > 10x Chandra) over the 0.1-7 keV band. WFXTs low-Earth orbit also minimizes the particle background. In five years of operation, WFXT will carry out three extragalactic surveys at unprecedented depth and address outstanding questions in astrophysics, cosmology and fundamental physics. In this article, we illustrate the mission concept and the connection between science requirements and mission parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Gas entropy in a representative sample of nearby X-ray galaxy clusters (REXCESS): relationship to gas mass fraction

(Abridged) We examine the radial entropy distribution and its scaling using 31 nearby galaxy clusters from the Representative XMM-Newton Cluster Structure Survey (REXCESS). The entropy profiles are robustly measured at least out to R_1000 in all systems and out to R_500 in 13 systems. Compared to theoretical expectations, the observed distributions show a radial and mass-dependent excess entropy that is greater and extends to larger radii in lower mass systems. At R_500, the mass dependence and entropy excess are both negligible within the uncertainties. Mirroring this behaviour, the scaling of gas entropy is shallower than self-similar in the inner regions, but steepens with radius, becoming consistent with self-similar at R_500. The dispersion in scaled entropy in the inner regions is linked to the presence of cool cores and dynamical activity; at larger radii the dispersion decreases by a factor of two and the dichotomy between subsamples disappears. Parameterising the profiles with a power law plus constant model, there are two peaks in central entropy K_0; however, we cannot distinguish between a bimodal or a left-skewed distribution. The outer slopes are correlated with system temperature; their distribution is unimodal with a median value of 0.98. Renormalising the dimensionless entropy profiles by the gas mass fraction profile f_gas(< R), leads to a remarkable reduction in the scatter, implying that gas mass fraction variations with radius and mass are the cause of the observed entropy properties. We discuss a tentative scenario to explain the behaviour of the entropy and gas mass fraction in the REXCESS sample, in which extra heating and merger mixing maintains an elevated central entropy level in the majority of the population, and a smaller fraction of systems develops a cool core.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Lyman Alpha Emitter Evolution in the Reionization Epoch

Combining cosmological SPH simulations with a previously developed Lyman Alpha production/transmission model and the Early Reionization Model (ERM, reionization ends at redshift z~7), we obtain Lyman Alpha and UV Luminosity Functions (LFs) for Lyman Alpha Emitters (LAEs) for redshifts between 5.7 and 7.6. Matching model results to observations at z~5.7 requires escape fractions of Lyman Alpha, f_alpha=0.3, and UV (non-ionizing) continuum photons, f_c=0.22, corresponding to a color excess, E(B-V)=0.15. We find that (i) f_c increases towards higher redshifts, due the decreasing mean dust content of galaxies, (ii) the evolution of f_alpha/f_c hints at the dust content of the ISM becoming progressively inhomogeneous/clumped with decreasing redshift. The clustering photoionization boost is important during the initial reionization phases but has little effect on the Lyman Alpha LF for a highly ionized IGM. Halo (stellar) masses are in the range 10.0 < \log M_h < 11.8 (8.1 < \log M_* < 10.4) with M_h \propto M_*^{0.64}. The star formation rates are between 3-120 solar masses per year, mass-weighted mean ages are greater than 20 Myr at all redshifts, while the mean stellar metallicity increases from Z=0.12 to 0.22 solar metallicity from z~7.6 to z~5.7; both age and metallicity positively correlate with stellar mass. The brightest LAEs are all characterized by large star formation rates and intermediate ages (~200 Myr), while objects in the faint end of the Lyman Alpha LF show large age and star formation rate spreads. With no more free parameters, the Spectral Energy Distributions of three LAE at z~5.7 observed by Lai et al. (2007) are well reproduced by an intermediate age (182-220 Myr) stellar population and the above E(B-V) value.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Simulating the effect of AGN feedback on the metal enrichment of galaxy clusters

We present a study of the effect of AGN feedback on metal enrichment and thermal properties of the intracluster medium (ICM) in hydrodynamical simulations. The cosmological simulations are performed for a set of clusters using a version of the TreePM-SPH Gadget code that follows chemo-dynamical evolution by accounting for metal enrichment by different stellar populations. Besides runs not including any efficient form of energy feedback, we carry out simulations including: (i) kinetic feedback in the form of galactic winds triggered by supernova explosions; (ii) AGN feedback from gas accretion onto super-massive black holes (BHs); (iii) AGN feedback in which a 'radio mode' is included. We find that AGN feedback is able to quench star formation in the brightest cluster galaxies at z<4 and provides correct temperature profiles in the central regions of galaxy groups. However, its effect is not sufficient to create cool cores in massive clusters. AGN feedback creates a widespread enrichment in the outskirts of clusters, thanks to its efficiency in displacing enriched gas from galactic halos to the inter-galactic medium at relatively high redshift. Iron abundance profiles are in better agreement with observations, with a more pristine enrichment of the ICM around and beyond the cluster virial regions. From the pattern of the relative abundances of Silicon and Iron, we conclude that a significant fraction of ICM enrichment in simulations is contributed by a diffuse population of intra-cluster stars. Our simulations also predict that profiles of Z_Si/Z_Fe abundance ratio do not increase at least out to 0.5 R_vir. Our results clearly show that different sources of energy feedback leave distinct imprints in the enrichment pattern of the ICM, that are more evident when looking at cluster external regions.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Substructure of the galaxy clusters in the REXCESS sample: observed statistics and comparison to numerical simulations

We study the substructure statistics of a representative sample of galaxy clusters by means of two currently popular substructure characterisation methods, power ratios and centroid shifts. We use the 31 clusters from the REXCESS sample, compiled from the southern ROSAT All-Sky cluster survey REFLEX with a morphologically unbiased selection in X-ray luminosity and redshift, all of which have been reobserved with XMM-Newton. We investigate the uncertainties of the substructure parameters and examine the dependence of the results on projection effects, finding that the uncertainties of the parameters can be quite substantial. Thus while the quantification of the dynamical state of individual clusters with these parameters should be treated with extreme caution, these substructure measures provide powerful statistical tools to characterise trends of properties in large cluster samples. The centre shift parameter, w, is found to be more sensitive in general. For the REXCESS sample neither the occurence of substructure nor the presence of cool cores depends on cluster mass. There is a significant anti-correlation between the existence of substantial substructure and cool cores. The simulated clusters show on average larger substructure parameters than the observed clusters, a trend that is traced to the fact that cool regions are more pronounced in the simulated clusters, leading to stronger substructure measures in merging clusters and clusters with offset cores. Moreover, the frequency of cool regions is higher in the simulations than in the observations, implying that the description of the physical processes shaping cluster formation in the simulations requires further improvement.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Substructures in hydrodynamical cluster simulations

The abundance and structure of dark matter subhalos has been analyzed extensively in recent studies of dark matter-only simulations, but comparatively little is known about the impact of baryonic physics on halo substructures. We here extend the SUBFIND algorithm for substructure identification such that it can be reliably applied to dissipative hydrodynamical simulations that include star formation. This allows, in particular, the identification of galaxies as substructures in simulations of clusters of galaxies, and a determination of their content of gravitationally bound stars, dark matter, and hot and cold gas. Using a large set of cosmological cluster simulations, we present a detailed analysis of halo substructures in hydrodynamical simulations of galaxy clusters, focusing in particular on the influence both of radiative and non-radiative gas physics, and of non-standard physics such as thermal conduction and feedback by galactic outflows. We also examine the impact of numerical nuisance parameters such as artificial viscosity parameterizations. We find that diffuse hot gas is efficiently stripped from subhalos when they enter the highly pressurized cluster atmosphere. This has the effect of decreasing the subhalo mass function relative to a corresponding dark matter-only simulation. These effects are mitigated in radiative runs, where baryons condense in the central subhalo regions and form compact stellar cores. However, in all cases, only a very small fraction, of the order of one percent, of subhalos within the cluster virial radii preserve a gravitationally bound hot gaseous atmosphere. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The cluster gas mass fraction as a cosmological probe: a revised study

(Abriged) We present the analysis of the baryonic content of 52 X-ray luminous galaxy clusters observed with Chandra in the redshift range 0.3-1.273. We use the deprojected X-ray surface brightness profiles and the measured values of the gas temperature to recover the gas and total mass profiles. By assuming that galaxy clusters are representative of the cosmic baryon budget, the distribution of the cluster baryon fraction in the hottest (T> 4 keV) systems as a function of redshift is used to constrain the cosmological parameters. We discuss how our constraints are affected by several systematics, namely the isothermality, the assumed baryon fraction in stars, the depletion parameter and the sample selection. By using only the cluster baryon fraction as a proxy for the cosmological parameters, we obtain that Omega is very well constrained at the value of 0.35 with a relative statistical uncertainty of 11% (1 sigma level; w=-1) and a further systematic error of about (-6,+7)%. On the other hand, constraints on Lambda (without the prior of flat geometry) and w (using the prior of flat geometry) are definitely weaker due to the presence of larger statistical and systematic uncertainties (of the order of 40 per cent on Lambda and larger than 50 per cent on w). If the WMAP 5-year best-fit results are assumed to fix the cosmological parameters, we limit the contributions expected from non-thermal pressure support and ICM clumpiness to be lower than about 10 per cent, leaving also room to accommodate baryons not accounted for either in the X-ray emitting plasma or in stars of the order of 18 per cent of the total cluster baryon budget.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Tracing the evolution in the iron content of the ICM

We present a Chandra analysis of the X-ray spectra of 56 clusters of galaxies at z>0.3, which cover a temperature range of 3>kT>15 keV. Our analysis is aimed at measuring the iron abundance in the ICM out to the highest redshift probed to date. We find that the emission-weighted iron abundance measured within (0.15-0.3)R_vir in clusters below 5 keV is, on average, a factor of ~2 higher than in hotter clusters, following Z(T)~0.88T^-(0.47)Z_o, which confirms the trend seen in local samples. We made use of combined spectral analysis performed over five redshift bins at 0.3>z>1.3 to estimate the average emission weighted iron abundance. We find a constant average iron abundance Z_Fe~0.25Z_o as a function of redshift, but only for clusters at z>0.5. The emission-weighted iron abundance is significantly higher (Z_Fe~0.4Z_o) in the redshift range z~0.3-0.5, approaching the value measured locally in the inner 0.15R_vir radii for a mix of cool-core and non cool-core clusters in the redshift range 0.1<z<0.3. The decrease in Z_Fe with redshift can be parametrized by a power law of the form ~(1+z)^(-1.25). The observed evolution implies that the average iron content of the ICM at the present epoch is a factor of ~2 larger than at z=1.2. We confirm that the ICM is already significantly enriched (Z_Fe~0.25Z_o) at a look-back time of 9 Gyr. Our data provide significant constraints on the time scales and physical processes that drive the chemical enrichment of the ICM.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

The Intracluster Medium in z > 1 Galaxy Clusters

The Chandra X-ray Observatory was used to obtain a 190 ks image of three high redshift galaxy clusters in one observation. The results of our analysis of these data are reported for the two z > 1 clusters in this Lynx field, including the most distant known X-ray selected cluster. Spatially-extended X-ray emission was detected from both these clusters, indicating the presence of hot gas in their intracluster media. A fit to the X-ray spectrum of RX J0849+4452, at z=1.26, yields a temperature of kT = 5.8^{+2.8}_{-1.7} keV. Using this temperature and the assumption of an isothermal sphere, the total mass of RX J0849+4452 is found to be 4.0^{+2.4}_{-1.9} X 10^{14} h_{65}^{-1} M_{\sun} within r = 1 h_{65}^{-1} Mpc. The T_x for RX J0849+4452 approximately agrees with the expectation based on its L_{bol} = 3.3^{+0.9}_{-0.5} X 10^{44}$ erg s$^{-1} according to the low redshift L_x - T_x relation. The very different distributions of X-ray emitting gas and of the red member galaxies in the two z > 1 clusters, in contrast to the similarity of the optical/IR colors of those galaxies, suggests that the early-type galaxies mostly formed before their host clusters.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Correlation Analysis of SFI Peculiar Velocities

We present results of a statistical analysis of the SFI catalog of peculiar velocities, a recently completed survey of spiral field galaxies with I-band Tully-Fisher distances (Haynes et al. 1999). The velocity field statistic utilized is the velocity correlation function, $ψ_1(r)$ (Gorski et al. 1989). The analysis is performed in redshift space, so as to circumvent potential ambiguities connected with inhomogeneous Malmquist bias corrections. The results from the SFI sample are compared with linear-theory predictions. We generate a large set of mock samples, extracted from N-body simulations, which are used to assess the reliability of our analysis and to estimate the associated uncertainties. Defining $η_8=σ_8Ω_0^{0.6}$, we find that the measured $ψ_1(r)$ implies a degenerate constraint in the $η_8-Γ$ plane, with $η_8=(0.3 +/- 0.1) (Γ/0.2)^{0.5}$, at the $2σ$ level, for the inverse Tully-Fisher (ITF) calibration presented in this paper. Model constraints are quite sensitive to the ITF calibration. The other ITF calibrations by Giovanelli et al. (1997) and da Costa et al. (1998) both give, for $Γ=0.2$, $η_8\simeq 0.6$ as the best-fitting value.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

The Observational Mass Function of Nearby Galaxy Clusters

We present a new determination of the mass function of galaxy clusters, based on optical virial mass estimates for a large sample of 152 nearby (z\le0.15) Abell-ACO clusters, as provided by Girardi et al. (1998). This sample includes both data from the literature and the new ENACS data. The resulting mass function is reliably estimated for masses larger than M_{lim}\simeq 4\times 10^{14} \msun, while it is affected by sample incompleteness at smaller masses. We find N(>M_{lim})=(6.3\pm 1.2) 10^{-6} (h^{-1} Mpc)^{-3} for cluster masses estimated within a 1.5 h^{-1} radius. Our mass function is intermediate between the two previous estimates by Bahcall & Cen (1993) and by Biviano et al. (1993). Based on the Press-Schechter approach, we use this mass function to constrain the amplitude of the fluctuation power spectrum at the cluster scale. After suitably convolving the PS predictions with observational errors on cluster masses and COBE-normalizing the fluctuation power spectrum, we find sigma_8=(0.60\pm 0.04) Omega_0^{-0.46+0.09 Omega_0} for flat low-density models and sigma_8=(0.60\pm 0.04) Omega_0^{-0.48+0.17 Omega_0} for open models (at the 90% c.l.).

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Minkowski Functionals of Abell/ACO Clusters

We determine the Minkowski functionals for a sample of Abell/ACO clusters, 401 with measured and 16 with estimated redshifts. The four Minkowski functionals (including the void probability function and the mean genus) deliver a global description of the spatial distribution of clusters on scales from $10$ to $60\hMpc$ with a clear geometric interpretation. Comparisons with mock catalogues of N--body simulations using different variants of the CDM model demonstrate the discriminative power of the description. The standard CDM model and the model with tilted perturbation spectrum cannot generate the Minkowski functionals of the cluster data, while a model with a cosmological constant and a model with breaking of the scale invariance of perturbations (BSI) yield compatible results.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

MASS SEGREGATION IN DARK MATTER MODELS.

We use the moments of counts of neighbors as given by the Generalized Correlation Integrals, to study the clustering properties of Dark Matter Halos (DH) in Cold Dark Matter (CDM) and Cold+Hot Dark Matter (CHDM) models. We compare the results with those found in the CfA and SSRS galaxy catalogs. We show that if we apply the analysis in redshift space, both models reproduce equally well the observed clustering of galaxies. Mass segregation is also found in the models: more massive DHs are more clustered compared with less massive ones. In redshift space, this mass segregation is reduced by a factor 2-3 due to the peculiar velocities. Observational catalogs give an indication of luminosity and size segregation, which is consistent with the predictions of the models. Because the mass segregation is smaller in redshift space, it is suggestive that the real luminosity or size segregation of galaxies could be significantly larger than what it is found in redshift catalogs.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Sampling Effects on Cosmological Dipoles

We use numerical simulations to investigate the behaviour of the dipole moment of the spatial distribution of different kinds of mass tracers. We select density peaks of the simulated matter distribution with mean separations of 38 and 30 $h^{-1}$ Mpc to represent two samples of rich clusters of galaxies. We extract, from the same simulations, samples selected to mimic the full 3D galaxy distribution of IR galaxies, and the flux--limited IRAS and QDOT galaxy samples. We compare the dipole moments of these ``galaxy'' and ``cluster'' samples in order to assess the effects of sampling uncertainties and shot--noise on the relationship between the ``true'' underlying galaxy dipole and the dipoles obtained for clusters and for the flux--limited galaxy samples. The results of this analysis demonstrate that the dipoles of both the IRAS and QDOT--like catalogues should trace the full 3D dipole shape fairly accurately, with the loss however of about 15--20$\%$ of the total 3D dipole amplitude. Furthermore, using a simple argument based on linear perturbation theory, on the linear biasing assumption and on the amplitude of the cluster dipole relative to that of galaxies, we can estimate their relative biasing factors quite accurately and in agreement with results obtained by other methods.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

The Cluster Distribution as a Test for Dark Matter Models. I: Clustering Properties

We present simulations of the cluster distribution in several dark matter models, using an optimized version of the truncated Zel'dovich approximation (TZA). We compare them with N-body cluster simulations and find that the TZA provides a very accurate description of the cluster distribution as long as fluctuations on the cluster mass scale are in the mildly non-linear regime. The simulated dark matter models are: Standard CDM (SCDM), Tilted CDM (TCDM) with n=0.7, Cold+Hot DM (CHDM) with 30% of hot component, low Hubble constant (h=0.3) CDM (LOWH) and a spatially flat low-density CDM model with Ω_0=0.2. We compare the simulations with a redshift sample of Abell/ACO clusters, using the integral of the 2-point correlation function and the probability density function. We find that the best models at reproducing the data are CHDM and LCDM. All the other models are ruled out. The reduced skewness S_3 is fairly constant with S_3=1.9, independent of the DM model and consistent with observational data. The abundances of clusters predicted using the Press--Schechter theory provide strong constraints: only the CHDM, LOWH and LCDM models appear to produce the correct number-density of clusters.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

The Cluster Distribution as a Test for Dark Matter Models. II: The Dipole Structure

We study the dipole for a redshift sample of Abell/ACO clusters. To elucidate the constraints it places on dark matter models, we use numerical simulations based on the Zel'dovich approximation. We run 20 realizations of each of six different dark matter models: four of these have a density parameter $Ω_0=1$, while the other two have $Ω_0=0.2$, one with and one without a cosmological constant term. For the Abell/ACO sample, we evaluate the parameter $β=Ω_0^{0.6}/b_{ cl}$, which reaches its asymptotic value at $R=160 Mpc/h$. We find $β=0.15\pm 0.04$ and $0.25\pm 0.06$ in LG and CMB frames respectively. After identifying in the simulations LG-like observers, we construct mock cluster samples. We find that an alignment between the cluster dipole and observer velocity, such as that observed (about 20 degrees) should not be expected necessarily in the simulations. This, together with the large observer-to-observer variance estimates of $β$, makes it difficult to place any firm constraints on cosmological models. Our results demonstrate that the large amplitude and convergence depth of the cluster dipole cannot be taken as strong evidence either for or against a low--density Universe.

0 citations0 reviewsNo code linkedOpen access
preprint1993arXiv

Angular Distribution of Clustersin Skewed CDM Models

We perform a detailed investigation of the statistical properties of the projected distribution of galaxy clusters obtained in Cold Dark Matter (CDM) models with both Gaussian and skewed primordial density fluctuations. We use N-body simulations to construct a set artificial Lick maps. An objective cluster--finding algorithm is used to identify clusters of different richness. For Gaussian models, the overall number of clusters is too small in the standard CDM case, but a model with higher normalisation fares much better; non--Gaussian models with negative skewness also fit faily well. We apply several statistical tests to compare real and simulated cluster samples, such as the 2-point correlation function, the minimal spanning tree construction, the multifractal analysis and the skewness of cell counts. The emerging picture is that Gaussian models, even with a higher normalization, are in trouble. Skew-positive models are also ruled out, while skew-negative models can reproduce the observed clustering of galaxy clusters in the CDM framework.

0 citations0 reviewsNo code linkedOpen access
preprint1993arXiv

Cluster Correlations in the Zel'dovich Approximation

We show how to simulate the clustering of rich clusters of galaxies using a technique based on the Zel'dovich approximation. This method well reproduces the spatial distribution of clusters obtainable from full N-body simulations at a fraction of the computational cost. We use an ensemble of large--scale simulations to assess the level and statistical significance of cluster clustering in open, tilted and flat versions of the Cold Dark Matter (CDM) model, as well as a model comprising a mixture of Cold and Hot Dark Matter (CHDM). We find the open and flat CDM models are excluded by the data. The tilted CDM model with a slight tilt is in marginal agreement, while larger tilt produces the right amount of clustering; CHDM is the best of all our models at reproducing the observations of cluster clustering. We find that {\em all} our models display a systematically weaker relationship between clustering length and mean cluster separation than seems to be implied by observations. We also note that the cluster bias factor, is not constant in any of the models, showing that one needs to be very careful when relating cluster clustering statistics to primordial density fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint1993arXiv

Non-Linear Clustering in the Cold+Hot Dark Matter Model

We use high resolution PM N-body simulations to follow the development of non-linear clustering in a flat Universe, dominated by Cold + Hot Dark Matter (CHDM) with 60% of CDM, 30% of HDM and 10% of baryons; a simulation box of 100 Mpc a side ($h=0.5$) is used. We analyze two CHDM simulations with $b =1.5$ (COBE normalization). We also compare them with CDM simulations with $b=1.5$ and $b=1$. We evaluate high-order correlation functions and the void-probability-function (VPF). Correlation functions are obtained both from counts in cells and counts of neighbors. The analysis is made for DM particles and for galaxies, identified as massive halos in the evolved density field. We also check the effects of dynamical evolution and redshift space distortions. We find that clustering of DM particles exhibits deviations from the hierarchical scaling, which decrease somewhat in redshift space. Galaxies follow hierarchical scaling far more closely, with coefficients S_3=2.5 and S_4=7.5, in general agreement with observational results, just marginally affected by redshift distortions and on the choices for the initial spectra. The hierarchical scaling of galaxy clustering is confirmed by the VPF analysis and in all the cases considered a good agreement with observational results is obtained.

0 citations0 reviewsNo code linkedOpen access