-
An N-body/hydrodynamical simulation study of the merging cluster El Gordo: A compelling case for self-interacting dark matter?
Authors:
R. Valdarnini
Abstract:
We use a large set N-body/hydrodynamical simulations to study the physical properties of the merging cluster El Gordo. We find that the observed X-ray structures, along with other data, can be matched fairly well by simulations with collision velocities 2,000 kms <= V <= 2,500 kms and impact parameters 600 kpc <= P <= 800 kpc. The mass of the primary is constrained to be between 10^{15} M_sun and…
▽ More
We use a large set N-body/hydrodynamical simulations to study the physical properties of the merging cluster El Gordo. We find that the observed X-ray structures, along with other data, can be matched fairly well by simulations with collision velocities 2,000 kms <= V <= 2,500 kms and impact parameters 600 kpc <= P <= 800 kpc. The mass of the primary is constrained to be between 10^{15} M_sun and ~ 1.6 10^{15} M_sun, in accordance with recent lensing-based mass measurements. Moreover, a returning, post-apocenter, scenario is not supported by our head-on simulations. We considered merger models that incorporate dark matter self-interactions. The simulation results show that the observed spatial offsets between the different mass components are well reproduced in self-interacting dark matter models with an elastic cross-section in the range σ_DM/m_X ~ 4 -5 cm^2/gr. In addition, the mean relative line-of-sight radial velocity between the two brightest cluster galaxies is found to be on the order of several hundred km/s. We argue that these findings provide an unambiguous signature of a dark matter behavior that exhibits collisional properties in a very energetic high-redshift cluster collision. The range of allowed values we find for sigma_DM/m_X is, however, inconsistent with present upper limits. To resolve this tension we suggest the possibility that the self-interacting dark matter model used here be considered as only a low order approximation, and that the underlying physical processes that describe the interaction of dark matter in major cluster mergers are more complex than can be adequately represented by the commonly assumed approach based on the scattering of dark matter particles.
△ Less
Submitted 9 February, 2024; v1 submitted 19 September, 2023;
originally announced September 2023.
-
A study of cool core resiliency and entropy mixing in simulations of galaxy cluster mergers
Authors:
R. Valdarnini,
C. L. Sarazin
Abstract:
We present results from a suite of binary merging cluster simulations. The hydrodynamical cluster simulations are performed employing a smoothed particle hydrodynamics (SPH) formulation in which gradient errors are strongly reduced by means of an integral approach. We consider adiabatic as well as radiative simulations, in which we include gas cooling, star formation and energy feedback from super…
▽ More
We present results from a suite of binary merging cluster simulations. The hydrodynamical cluster simulations are performed employing a smoothed particle hydrodynamics (SPH) formulation in which gradient errors are strongly reduced by means of an integral approach. We consider adiabatic as well as radiative simulations, in which we include gas cooling, star formation and energy feedback from supernovae. We explore the effects of merging on the thermodynamic structure of the intracluster gas of the final merger remnant. In particular, we study how core entropy is generated during the merging and the stability properties of the initial cool-core profile against disruption. To this end, we consider a range of initial mass ratio and impact parameters. Final entropy profiles of our adiabatic merging simulations are in good accord with previous findings (ZuHone 2011), with cool-cores being disrupted for all of the initial merging setups. For equal-mass off-axis mergers, we find that a significant contribution to the final primary core entropy is due to hydrodynamic instabilities generated by rotational motions, which are induced by tidal torques during the first pericenter passage. In radiative simulations, cool-cores are more resilient against heating processes; nonetheless, they are able to maintain their integrity only in the case of off-axis mergers with very unequal masses. We suggest that these results are robust against changes in the gas physical modeling, in particular to the inclusion of AGN thermal feedback.
Our findings support the view that the observed core cluster morphology emerges naturally in a merging cluster context, and conclude that the merging angular momentum is a key parameter in shaping the thermodynamical properties of the final merger remnant.
△ Less
Submitted 11 May, 2021; v1 submitted 16 April, 2021;
originally announced April 2021.
-
Einstein, Planck and Vera Rubin: relevant encounters between the Cosmological and the Quantum Worlds
Authors:
Paolo Salucci,
Giampiero Esposito,
Gaetano Lambiase,
Emmanuele Battista,
Micol Benetti,
Donato Bini,
Lumen Boco,
Gauri Sharma,
Valerio Bozza,
Luca Buoninfante,
Antonio Capolupo,
Salvatore Capozziello,
Giovanni Covone,
Rocco D'Agostino,
Mariafelicia DeLaurentis,
Ivan De Martino,
Giulia De Somma,
Elisabetta Di Grezia,
Chiara Di Paolo,
Lorenzo Fatibene,
Viviana Gammaldi,
Andrea Geralico,
Lorenzo Ingoglia,
Andrea Lapi,
Giuseppe G. Luciano
, et al. (16 additional authors not shown)
Abstract:
In Cosmology and in Fundamental Physics there is a crucial question like: where the elusive substance that we call Dark Matter is hidden in the Universe and what is it made of?, that, even after 40 years from the Vera Rubin seminal discovery does not have a proper answer. Actually, the more we have investigated, the more this issue has become strongly entangled with aspects that go beyond the esta…
▽ More
In Cosmology and in Fundamental Physics there is a crucial question like: where the elusive substance that we call Dark Matter is hidden in the Universe and what is it made of?, that, even after 40 years from the Vera Rubin seminal discovery does not have a proper answer. Actually, the more we have investigated, the more this issue has become strongly entangled with aspects that go beyond the established Quantum Physics, the Standard Model of Elementary particles and the General Relativity and related to processes like the Inflation, the accelerated expansion of the Universe and High Energy Phenomena around compact objects. Even Quantum Gravity and very exotic DM particle candidates may play a role in framing the Dark Matter mystery that seems to be accomplice of new unknown Physics. Observations and experiments have clearly indicated that the above phenomenon cannot be considered as already theoretically framed, as hoped for decades. The Special Topic to which this review belongs wants to penetrate this newly realized mystery from different angles, including that of a contamination of different fields of Physics apparently unrelated. We show with the works of this ST that this contamination is able to guide us into the required new Physics. This review wants to provide a good number of these "paths or contamination" beyond/among the three worlds above; in most of the cases, the results presented here open a direct link with the multi-scale dark matter phenomenon, enlightening some of its important aspects. Also in the remaining cases, possible interesting contacts emerges.
△ Less
Submitted 16 November, 2020;
originally announced November 2020.
-
A multifiltering study of turbulence in a large sample of simulated galaxy clusters
Authors:
R. Valdarnini
Abstract:
We present results from a large set of N-body/SPH hydrodynamical cluster simulations aimed at studying the statistical properties of turbulence in the ICM. The numerical hydrodynamical scheme employs a SPH formulation in which gradient errors are strongly reduced by using an integral approach. We consider both adiabatic and radiative simulations. We construct clusters subsamples according to the c…
▽ More
We present results from a large set of N-body/SPH hydrodynamical cluster simulations aimed at studying the statistical properties of turbulence in the ICM. The numerical hydrodynamical scheme employs a SPH formulation in which gradient errors are strongly reduced by using an integral approach. We consider both adiabatic and radiative simulations. We construct clusters subsamples according to the cluster dynamical status or gas physical modeling, from which we extract small-scale turbulent velocities obtained by applying to cluster velocities different multiscale filtering methods. The velocity power spectra of non-radiative relaxed clusters are mostly solenoidal and steeper than Kolgomorov. Cooling runs are distinguished by much shallower spectra, a feature which we interpret as the injection of turbulence at small scales due to the interaction of compact cool gas cores with the ICM. Turbulence in galaxy clusters is then characterized by multiple injection scales, with the small scale driving source acting in addition to the large scale injection mechanisms. Cooling runs of relaxed clusters exhibit enstrophy profiles with a power-law behavior over more than two decades in radius, and a turbulent-to-thermal energy ratio ~1 %. In accord with Hitomi observations, in the core of a highly relaxed cluster we find low level of gas motions. In addition, the estimated cluster radial profile of the sloshing oscillation period exhibits an associated Froude number satisfying Fr ~ 0.1 within r / r_{200} <~ 0.1. Our findings suggest that in cluster cores ICM turbulence approaches a stratified anisotropic regime, with weak stirring motions dominated by gravity buoyancy forces and strongly suppressed along the radial direction. We conclude that turbulent heating cannot be considered the main heating source in cluster cores.
△ Less
Submitted 19 February, 2019;
originally announced February 2019.
-
Improved Performances in Subsonic Flows of an SPH Scheme with Gradients Estimated using an Integral Approach
Authors:
R. Valdarnini
Abstract:
In this paper we present results from a series of hydrodynamical tests aimed at validating the performance of a smoothed particle hydrodynamics (SPH) formulation in which gradients are derived from an integral approach. We specifically investigate the code behavior with subsonic flows, where it is well known that zeroth-order inconsistencies present in standard SPH make it particularly problematic…
▽ More
In this paper we present results from a series of hydrodynamical tests aimed at validating the performance of a smoothed particle hydrodynamics (SPH) formulation in which gradients are derived from an integral approach. We specifically investigate the code behavior with subsonic flows, where it is well known that zeroth-order inconsistencies present in standard SPH make it particularly problematic to correctly model the fluid dynamics. In particular we consider the Gresho-Chan vortex problem, the growth of Kelvin-Helmholtz instabilities, the statistics of driven subsonic turbulence and the cold Keplerian disc problem. We compare simulation results for the different tests with those obtained, for the same initial conditions, using standard SPH. We also compare the results with the corresponding ones obtained previously with other numerical methods, such as codes based on a moving-mesh scheme or Godunov-type Lagrangian meshless methods. We quantify code performances by introducing error norms and spectral properties of the particle distribution, in a way similar to what was done in other works. We find that the new SPH formulation exhibits strongly reduced gradient errors and outperforms standard SPH in all of the tests considered. In fact, in terms of accuracy we find good agreement between the simulation results of the new scheme and those produced using other recently proposed numerical schemes. These findings suggest that the proposed method can be successfully applied for many astrophysical problems in which the presence of subsonic flows previously limited the use of SPH, with the new scheme now being competitive in these regimes with other numerical methods.
△ Less
Submitted 31 August, 2016; v1 submitted 30 August, 2016;
originally announced August 2016.
-
Overlapping inflows as catalysts of AGN activity - II: Relative importance of turbulence and inflow-disc interaction
Authors:
Juan Manuel Carmona-Loaiza,
Monica Colpi,
Massimo Dotti,
Riccardo Valdarnini
Abstract:
The main challenge for understanding the fuelling of supermassive black holes in active galactic nuclei is not to account for the source of fuel, but rather to explain its delivery from the boundaries of the black hole sphere of influence (10-100 pc) down to sub-parsec scales. In this work, we report on a series of numerical experiments aimed at exploring in further depth our model of "overlapping…
▽ More
The main challenge for understanding the fuelling of supermassive black holes in active galactic nuclei is not to account for the source of fuel, but rather to explain its delivery from the boundaries of the black hole sphere of influence (10-100 pc) down to sub-parsec scales. In this work, we report on a series of numerical experiments aimed at exploring in further depth our model of "overlapping inflow events" as catalysts for rapid accretion, seeding a turbulent field in the infalling gas. We initially set a gaseous shell in non-equilibrium rotation around a supermassive black hole. After infall, the shell stalls in a disc-like structure. A second shell is then set in either co-rotation or counter-rotation with respect to the first and is let to impinge on the previously-formed disc. We find that combined turbulence and overlap significantly enhance accretion in counter-rotating inflows, while turbulence dominates for co-rotating inflows. The leftovers of overlapping inflows are warped nuclear discs, whose morphology depends on the relative orientation and angular momentum of the disc and the shell. Overlapping inflows leave observational signatures in the gas rotation curves.
△ Less
Submitted 31 July, 2015;
originally announced August 2015.
-
The role of the artificial conductivity in SPH simulations of galaxy clusters: effects on the ICM properties
Authors:
Veronica Biffi,
Riccardo Valdarnini
Abstract:
We study the thermal structure of the intra-cluster medium (ICM) in a set of cosmological hydrodynamical cluster simulations performed with an SPH numerical scheme employing an artificial conductivity (AC) term. We explore the effects of this term on the ICM temperature and entropy profiles, thermal distribution, velocity field and expected X-ray emission. We find that in adiabatic runs the artifi…
▽ More
We study the thermal structure of the intra-cluster medium (ICM) in a set of cosmological hydrodynamical cluster simulations performed with an SPH numerical scheme employing an artificial conductivity (AC) term. We explore the effects of this term on the ICM temperature and entropy profiles, thermal distribution, velocity field and expected X-ray emission. We find that in adiabatic runs the artificial conductivity favours (i) the formation of an entropy core, raising and flattening the central entropy profiles, in better agreement with findings from Eulerian codes; and (ii) a systematic reduction of the cold gas component. In fact, the cluster large-scale structure and dynamical state are preserved across different runs, but the improved gas mixing enabled by the AC term strongly increases the stripping rate of gas from the cold clumps moving through the ICM. This in turn reduces the production of turbulence generated by the instabilities which develop because of the interaction between clumps and ambient ICM. We then find that turbulent motions, enhanced by the time-dependent artificial viscosity scheme we use, are rather damped by the AC term. The ICM synthetic X-ray emission substantially mirrors the changes in its thermo-dynamical structure, stressing the robustness of the AC impact. All these effects are softened by the introduction of radiative cooling but still present, especially a partial suppression of cold gas. Therefore, not only the physics accounted for, but also the numerical approach itself can have an impact in shaping the ICM thermo-dynamical structure and ultimately in the use of SPH cluster simulations for cosmological studies.
△ Less
Submitted 30 October, 2014;
originally announced October 2014.
-
The MUSIC of Galaxy Clusters II: X-ray global properties and scaling relations
Authors:
Veronica Biffi,
Federico Sembolini,
Marco De Petris,
Riccardo Valdarnini,
Gustavo Yepes,
Stefan Gottlöber
Abstract:
We present the X-ray properties and scaling relations of a large sample of clusters extracted from the Marenostrum MUltidark SImulations of galaxy Clusters (MUSIC) dataset. We focus on a sub-sample of 179 clusters at redshift z~0.11, with 3.2e14M_sun/h<M_vir<2e15Msun/h, complete in mass. We employed the X-ray photon simulator PHOX to obtain synthetic Chandra Observations and derive observable-like…
▽ More
We present the X-ray properties and scaling relations of a large sample of clusters extracted from the Marenostrum MUltidark SImulations of galaxy Clusters (MUSIC) dataset. We focus on a sub-sample of 179 clusters at redshift z~0.11, with 3.2e14M_sun/h<M_vir<2e15Msun/h, complete in mass. We employed the X-ray photon simulator PHOX to obtain synthetic Chandra Observations and derive observable-like global properties of the intracluster medium (ICM), as X-ray temperature (T_X) and luminosity (L_X). T_X is found to slightly under-estimate the true mass-weighted temperature, although tracing fairly well the cluster total mass. We also study the effects of T_X on scaling relations with cluster intrinsic properties: total (M_500) and gas (M_g500) mass; integrated Compton parameter (Y_SZ) of the Sunyaev-Zel'dovich (SZ) thermal effect; Y_X=M_g500 T_X. We confirm that Y_X is a very good mass proxy, with a scatter on M_500-Y_X and Y_SZ-Y_X lower than 5%. The study of scaling relations among X-ray, intrinsic and SZ properties indicates that MUSIC clusters reasonably resemble the self-similar prediction, especially for correlations involving T_X. The observational approach also allows for a more direct comparison with real clusters, from which we find deviations mainly due to the physical description of the ICM, affecting T_X and, particularly, L_X.
△ Less
Submitted 13 January, 2014;
originally announced January 2014.
-
Overlapping Inflow Events as Catalysts for Supermassive Black Hole Growth
Authors:
Juan Manuel Carmona-Loaiza,
Monica Colpi,
Massimo Dotti,
Riccardo Valdarnini
Abstract:
One of the greatest issues in modelling black hole fuelling is our lack of understanding of the processes by which gas loses angular momentum and falls from galactic scales down to the nuclear region where an accretion disc forms, subsequently guiding the inflow of gas down to the black hole horizon. It is feared that gas at larger scales might still retain enough angular momentum and settle into…
▽ More
One of the greatest issues in modelling black hole fuelling is our lack of understanding of the processes by which gas loses angular momentum and falls from galactic scales down to the nuclear region where an accretion disc forms, subsequently guiding the inflow of gas down to the black hole horizon. It is feared that gas at larger scales might still retain enough angular momentum and settle into a larger scale disc with very low or no inflow to form or replenish the inner accretion disc (on ~0.01 pc scales). In this paper we report on hydrodynamical simulations of rotating infalling gas shells impacting at different angles onto a pre-existing, primitive large scale (~10 pc) disc around a super-massive black hole. The aim is to explore how the interaction between the shell and the disc redistributes the angular momentum on scales close to the black hole's sphere of influence. Angular momentum redistribution via hydrodynamical shocks leads to inflows of gas across the inner boundary, enhancing the inflow rate by more than 2-3 orders of magnitude. In all cases, the gas inflow rate across the inner parsec is higher than in the absence of the interaction, and the orientation of the angular momentum of the flow in the region changes with time due to gas mixing. Warped discs or nested misaligned rings form depending on the angular momentum content of the infalling shell relative to the disc. In the cases in which the shell falls in near counter-rotation, part of the resulting flows settle into an inner dense disc which becomes more susceptible to mass transfer.
△ Less
Submitted 22 December, 2013; v1 submitted 9 December, 2013;
originally announced December 2013.
-
Hydrodynamic capabilities of an SPH code incorporating an artificial conductivity term with a gravity-based signal velocity
Authors:
R. Valdarnini
Abstract:
This paper investigates the hydrodynamic performances of an SPH code incorporating an artificial heat conductivity term in which the adopted signal velocity is applicable when gravity is present. In accordance with previous findings it is shown that the performances of SPH to describe the development of Kelvin-Helmholtz instabilities depend strongly on the consistency of the initial condition set-…
▽ More
This paper investigates the hydrodynamic performances of an SPH code incorporating an artificial heat conductivity term in which the adopted signal velocity is applicable when gravity is present. In accordance with previous findings it is shown that the performances of SPH to describe the development of Kelvin-Helmholtz instabilities depend strongly on the consistency of the initial condition set-up and on the leading error in the momentum equation due to incomplete kernel sampling. An error and stability analysis shows that the quartic B-spline kernel (M_5) possesses very good stability properties and we propose its use with a large neighbor number, between ~50 (2D) to ~ 100 (3D), to improve convergence in simulation results without being affected by the so-called clumping instability. SPH simulations of the blob test show that in the regime of strong supersonic flows an appropriate limiting condition, which depends on the Prandtl number, must be imposed on the artificial conductivity SPH coefficients in order to avoid an unphysical amount of heat diffusion. Results from hydrodynamic simulations that include self-gravity show profiles of hydrodynamic variables that are in much better agreement with those produced using mesh-based codes. In particular, the final levels of core entropies in cosmological simulations of galaxy clusters are consistent with those found using AMR codes. Finally, results of the Rayleigh-Taylor instability test demonstrate that in the regime of very subsonic flows the code has still several difficulties in the treatment of hydrodynamic instabilities. These problems being intrinsically due to the way in which in standard SPH gradients are calculated and not to the implementation of the artificial conductivity term.
△ Less
Submitted 30 July, 2012;
originally announced July 2012.
-
The impact of numerical viscosity in SPH simulations of galaxy clusters
Authors:
R. Valdarnini
Abstract:
A SPH code employing a time-dependent artificial viscosity scheme is used to construct a large set of N-body/SPH cluster simulations for studying the impact of artificial viscosity on the thermodynamics of the ICM and its velocity field statistical properties. Spectral properties of the gas velocity field are investigated by measuring for the simulated clusters the velocity power spectrum E(k). Th…
▽ More
A SPH code employing a time-dependent artificial viscosity scheme is used to construct a large set of N-body/SPH cluster simulations for studying the impact of artificial viscosity on the thermodynamics of the ICM and its velocity field statistical properties. Spectral properties of the gas velocity field are investigated by measuring for the simulated clusters the velocity power spectrum E(k). The longitudinal component E_c(k) exhibits over a limited range a Kolgomorov-like scaling k^{-5/3}, whilst the solenoidal power spectrum component E_s(k) is strongly influenced by numerical resolution effects. The dependence of the spectra E(k) on dissipative effects is found to be significant at length scales 100-300Kpc, with viscous damping of the velocities being less pronounced in those runs with the lowest artificial viscosity. The turbulent energy density radial profile E_{turb}(r) is strongly affected by the numerical viscosity scheme adopted in the simulations, with the turbulent-to-total energy density ratios being higher in the runs with the lowest artificial viscosity settings and lying in the range between a few percent and ~10%. These values are in accord with the corresponding ratios extracted from previous cluster simulations realized using mesh-based codes. At large cluster radii, the mass correction terms to the hydrostatic equilibrium equation are little affected by the numerical viscosity of the simulations, showing that the X-ray mass bias is already estimated well in standard SPH simulations. Finally, simulations in which the gas can cool radiatively are characterized by the presence in the cluster inner regions of high levels of turbulence, generated by the interaction of the compact cool gas core with the ambient medium.
△ Less
Submitted 16 October, 2010;
originally announced October 2010.
-
LoCuSS: A Comparison of Cluster Mass Measurements from XMM-Newton and Subaru - Testing Deviation from Hydrostatic Equilibrium and Non-Thermal Pressure Support
Authors:
Yu-Ying Zhang,
Nobuhiro Okabe,
Alexis Finoguenov,
Graham P. Smith,
Rocco Piffaretti,
Riccardo Valdarnini,
Arif Babul,
August E. Evrard,
Pasquale Mazzotta,
Alastair J. R. Sanderson,
Daniel P. Marrone
Abstract:
We compare X-ray hydrostatic and weak-lensing mass estimates for a sample of 12 clusters that have been observed with both XMM-Newton and Subaru. At an over-density of Δ=500, we obtain 1-M^{X}/M^{WL}=0.01+/-0.07 for the whole sample. We also divided the sample into undisturbed and disturbed sub-samples based on quantitative X-ray morphologies using asymmetry and fluctuation parameters, obtaining…
▽ More
We compare X-ray hydrostatic and weak-lensing mass estimates for a sample of 12 clusters that have been observed with both XMM-Newton and Subaru. At an over-density of Δ=500, we obtain 1-M^{X}/M^{WL}=0.01+/-0.07 for the whole sample. We also divided the sample into undisturbed and disturbed sub-samples based on quantitative X-ray morphologies using asymmetry and fluctuation parameters, obtaining 1-M^{X}/M^{WL}=0.09+/-0.06 and -0.06+/-0.12 for the undisturbed and disturbed clusters, respectively. In addition to non-thermal pressure support, there may be a competing effect associated with adiabatic compression and/or shock heating which leads to overestimate of X-ray hydrostatic masses for disturbed clusters, for example, in the famous merging cluster A1914. Despite the modest statistical significance of the mass discrepancy, on average, in the undisturbed clusters, we detect a clear trend of improving agreement between M^{X} and M^{WL} as a function of increasing over-density, M^{X}/M^{WL}=(0.908+/-0.004)+(0.187+/-0.010) \log_{10}(Δ/500). We also examine the gas mass fractions, f_{gas}=M^{gas}/M^{WL}, finding that they are an increasing function of cluster radius, with no dependence on dynamical state, in agreement with predictions from numerical simulations. Overall, our results demonstrate that XMM-Newton and Subaru are a powerful combination for calibrating systematic uncertainties in cluster mass measurements.
△ Less
Submitted 23 February, 2010; v1 submitted 5 January, 2010;
originally announced January 2010.
-
Measurement of the dark matter velocity anisotropy in galaxy clusters
Authors:
Ole Host,
Steen H. Hansen,
Rocco Piffaretti,
Andrea Morandi,
Stefano Ettori,
Scott T. Kay,
Riccardo Valdarnini
Abstract:
The internal dynamics of a dark matter structure may have the remarkable property that the local temperature in the structure depends on direction. This is parametrized by the velocity anisotropy beta which must be zero for relaxed collisional structures, but has been shown to be non-zero in numerical simulations of dark matter structures. Here we present a method to infer the radial profile of…
▽ More
The internal dynamics of a dark matter structure may have the remarkable property that the local temperature in the structure depends on direction. This is parametrized by the velocity anisotropy beta which must be zero for relaxed collisional structures, but has been shown to be non-zero in numerical simulations of dark matter structures. Here we present a method to infer the radial profile of the velocity anisotropy of the dark matter halo in a galaxy cluster from X-ray observables of the intracluster gas. This non-parametric method is based on a universal relation between the dark matter temperature and the gas temperature which is confirmed through numerical simulations. We apply this method to observational data and we find that beta is significantly different from zero at intermediate radii. Thus we find a strong indication that dark matter is effectively collisionless on the dynamical time-scale of clusters, which implies an upper limit on the self-interaction cross-section per unit mass sigma/m < 1 cm2/g. Our results may provide an independent way to determine the stellar mass density in the central regions of a relaxed cluster, as well as a test of whether a cluster is in fact relaxed.
△ Less
Submitted 14 August, 2008;
originally announced August 2008.
-
Total mass biases in X-ray galaxy clusters
Authors:
Rocco Piffaretti,
Riccardo Valdarnini
Abstract:
The exploitation of clusters of galaxies as cosmological probes relies on accurate measurements of their total gravitating mass. X-ray observations provide a powerful means of probing the total mass distribution in galaxy clusters, but might be affected by observational biases and rely on simplistic assumptions originating from our limited understanding of the intracluster medium physics. This p…
▽ More
The exploitation of clusters of galaxies as cosmological probes relies on accurate measurements of their total gravitating mass. X-ray observations provide a powerful means of probing the total mass distribution in galaxy clusters, but might be affected by observational biases and rely on simplistic assumptions originating from our limited understanding of the intracluster medium physics. This paper is aimed at elucidating the reliability of X-ray total mass estimates in clusters of galaxies by properly disentangling various biases of both observational and physical origin. We use N-body/SPH simulation of a large sample of ~100 galaxy clusters and investigate total mass biases by comparing the mass reconstructed adopting an observational-like approach with the true mass in the simulations. X-ray surface brightness and temperature profiles extracted from the simulations are fitted with different models and adopting different radial fitting ranges in order to investigate modeling and extrapolation biases. Different theoretical definitions of gas temperature are used to investigate the effect of spectroscopic temperatures and a power ratio analysis of the surface brightness maps allows us to assess the dependence of the mass bias on cluster dynamical state. Moreover, we perform a study on the reliability of hydrostatic and hydrodynamical equilibrium mass estimates using the full three-dimensional information in the simulation.[abridged]
△ Less
Submitted 7 August, 2008;
originally announced August 2008.
-
X-ray temperature spectroscopy of simulated cooling clusters
Authors:
R. Valdarnini
Abstract:
Results from a large sample of hydrodynamical/N-body simulations of galaxy clusters in a LCDM cosmology are used to simulate cluster X-ray observations as expected from Chandra observations. The physical modeling of the gas includes radiative cooling, star formation, energy feedback and metal enrichment. The biasing of spectral temperatures with respect to mass-weighted temperatures is found to…
▽ More
Results from a large sample of hydrodynamical/N-body simulations of galaxy clusters in a LCDM cosmology are used to simulate cluster X-ray observations as expected from Chandra observations. The physical modeling of the gas includes radiative cooling, star formation, energy feedback and metal enrichment. The biasing of spectral temperatures with respect to mass-weighted temperatures is found to be influenced by two independent processes. The first scale dependency is absent in adiabatic runs and is due to cooling, whose efficiency to transform cold gas into stars is higher for cool clusters and this in turn implies a strong dependency of the spectral versus mass-weighted temperature relation on the cluster mass. The second dependency is due to photon emission because of cool gas which is accreted during merging events and biases the spectral fits. These events have been quantified according to the power ratio method and a robust correlation is found to exist between the spectral bias and the amount of cluster substructure. The shape of the simulated temperature profiles is not universal and it is steeper at the cluster center for cool clusters than for the massive ones. The profiles are in good agreement with data in the radial range between $\sim 0.1 r_{vir}$ and $\sim 0.4 r_{vir}$; at small radii ($r< 0.1 r_{vir}$) the cooling runs fail to reproduce the shape of the observed profiles. The fit is improved if one considers a hierarchical merging scenario in which cluster cores can accrete cooler gas through merging with cluster subclumps, though the shape of the temperature profiles is modified in a significant way only in the regime where the mass of the substructure is a large fraction of the cluster mass.
△ Less
Submitted 17 August, 2006; v1 submitted 20 June, 2006;
originally announced June 2006.
-
X-ray Cluster Properties in SPH Simulations of Galaxy Clusters
Authors:
R. Valdarnini
Abstract:
Results from a large set of hydrodynamical SPH simulations of galaxy clusters in a flat LCDM cosmology are used to investigate cluster X-ray properties. The physical modeling of the gas includes radiative cooling, star formation, energy feedback and metal enrichment that follows from the explosions of SNe type II and Ia. The metallicity dependence of the cooling function is also taken into accou…
▽ More
Results from a large set of hydrodynamical SPH simulations of galaxy clusters in a flat LCDM cosmology are used to investigate cluster X-ray properties. The physical modeling of the gas includes radiative cooling, star formation, energy feedback and metal enrichment that follows from the explosions of SNe type II and Ia. The metallicity dependence of the cooling function is also taken into account. It is found that the luminosity-temperature relation of simulated clusters is in good agreement with the data, and the X-ray properties of cool clusters are unaffected by the amount of feedback energy that has heated the intracluster medium (ICM). The fraction of hot gas $f_g$ at the virial radius increases with $T_X$ and the distribution obtained from the simulated cluster sample is consistent with the observational ranges.
△ Less
Submitted 8 October, 2003;
originally announced October 2003.
-
Parallelization of a treecode
Authors:
R. Valdarnini
Abstract:
I describe here the performance of a parallel treecode with individual particle timesteps. The code is based on the Barnes-Hut algorithm and runs cosmological N-body simulations on parallel machines with a distributed memory architecture using the MPI message-passing library. For a configuration with a constant number of particles per processor the scalability of the code was tested up to P=128…
▽ More
I describe here the performance of a parallel treecode with individual particle timesteps. The code is based on the Barnes-Hut algorithm and runs cosmological N-body simulations on parallel machines with a distributed memory architecture using the MPI message-passing library. For a configuration with a constant number of particles per processor the scalability of the code was tested up to P=128 processors on an IBM SP4 machine. In the large $P$ limit the average CPU time per processor necessary for solving the gravitational interactions is $\sim 10 %$ higher than that expected from the ideal scaling relation. The processor domains are determined every large timestep according to a recursive orthogonal bisection, using a weighting scheme which takes into account the total particle computational load within the timestep. The results of the numerical tests show that the load balancing efficiency $L$ of the code is high ($>=90%$) up to P=32, and decreases to $L\sim 80%$ when P=128. In the latter case it is found that some aspects of the code performance are affected by machine hardware, while the proposed weighting scheme can achieve a load balance as high as $L\sim 90%$ even in the large $P$ limit.
△ Less
Submitted 18 March, 2003;
originally announced March 2003.
-
Performance characteristics of a parallel treecode
Authors:
R. Valdarnini
Abstract:
I describe here the performances of a parallel treecode with individual particle timesteps. The code is based on the Barnes-Hut algorithm and runs cosmological N-body simulations on parallel machines with a distributed memory architecture using the MPI message passing library. For a configuration with a constant number of particles per processor the scalability of the code has been tested up to…
▽ More
I describe here the performances of a parallel treecode with individual particle timesteps. The code is based on the Barnes-Hut algorithm and runs cosmological N-body simulations on parallel machines with a distributed memory architecture using the MPI message passing library. For a configuration with a constant number of particles per processor the scalability of the code has been tested up to P=32 processors. The average CPU time per processor necessary for solving the gravitational interactions is within $\sim 10 %$ of that expected from the ideal scaling relation. The load balancing efficiency is high ($\simgt90%$) if the processor domains are determined every large timestep according to a weighting scheme which takes into account the total particle computational load within the timestep.
△ Less
Submitted 11 December, 2002;
originally announced December 2002.
-
Iron abundances in hydrodynamical simulations of galaxy clusters
Authors:
R. Valdarnini
Abstract:
Hydrodynamical SPH simulations of galaxy clusters are used to investigate the metal enrichment of the intracluster medium. The final metallicity abundances of the simulated clusters are determined according to the numerical resolution and a number of model parameters. For a fiducial set of model prescriptions the results of the simulations indicate iron abundances in broad agreement with data. F…
▽ More
Hydrodynamical SPH simulations of galaxy clusters are used to investigate the metal enrichment of the intracluster medium. The final metallicity abundances of the simulated clusters are determined according to the numerical resolution and a number of model parameters. For a fiducial set of model prescriptions the results of the simulations indicate iron abundances in broad agreement with data. Final X-ray properties are not sensitive to the heating of the ICM. This supports a scenario where the ICM evolution of cool clusters is driven by radiative cooling.
△ Less
Submitted 6 November, 2002;
originally announced November 2002.
-
Iron abundances and heating of the ICM in hydrodynamical simulations of galaxy clusters
Authors:
R. Valdarnini
Abstract:
Results from a large set of hydrodynamical SPH simulations of galaxy clusters in a flat LCDM cosmology are used to investigate the metal enrichment and heating of the ICM. The physical modeling of the gas includes radiative cooling, star formation, energy feedback and metal enrichment that follow from the explosions of SNe of type II and Ia. The metallicity dependence of the cooling function is…
▽ More
Results from a large set of hydrodynamical SPH simulations of galaxy clusters in a flat LCDM cosmology are used to investigate the metal enrichment and heating of the ICM. The physical modeling of the gas includes radiative cooling, star formation, energy feedback and metal enrichment that follow from the explosions of SNe of type II and Ia. The metallicity dependence of the cooling function is also taken into account. For a fiducial set of model prescriptions the results indicate radial iron profiles in broad agreement with observations; global iron abundances are also consistent with data. It is found that the iron distribution in the ICM is critically dependent on the shape of the metal deposition profile. For low temperatu re clusters simulations yield iron abundances below the allowed observational range, unless it is introduced a minimum diffusion length of metals in the ICM. The simulated emission-weighted radial temperature profiles are in good agreement with data for cooling flow clusters, but at very small distances from the cluster centres ($\sim 2%$ of the virial radii) the temperatures are a factor $\sim$ two higher than the measured spectral values. The luminosity-temperature relation is in excellent agreement with the data, cool clusters ($T_X\sim 1keV$) have a core excess entropy of $\sim 200 keVcm^2$ and their X-ray properties are unaffected by the amount of feedback energy that has heated the ICM. The fraction of hot gas $f_g$ at the virial radius increases with $T_X$ and the distribution obtained from the simulated cluster sample is consistent with the observational ranges.
△ Less
Submitted 29 October, 2002; v1 submitted 11 October, 2002;
originally announced October 2002.
-
Radiative cooling in SPH hydrodynamical simulations of X-ray clusters
Authors:
R. Valdarnini
Abstract:
The results from hydrodynamical TREESPH simulations of galaxy clusters are used to investigate the dependence of the final cluster X-ray properties upon the numerical resolution and the assumed star formation models for the cooled gas. When cold gas particles are allowed to convert into stars the final gas profiles show a well defined core radius and the temperature profiles are nearly flat. Fin…
▽ More
The results from hydrodynamical TREESPH simulations of galaxy clusters are used to investigate the dependence of the final cluster X-ray properties upon the numerical resolution and the assumed star formation models for the cooled gas. When cold gas particles are allowed to convert into stars the final gas profiles show a well defined core radius and the temperature profiles are nearly flat. Final X-ray luminosities are found to be numerically stable, with uncertainties of a factor 2.
△ Less
Submitted 7 November, 2001;
originally announced November 2001.
-
Detection of non-random patterns in large-scale structure
Authors:
R. Valdarnini
Abstract:
A new method for analyzing the morphological features of point patterns is presented. The method is taken from the study of molecular liquids, where it has been introduced for making a statistical description of anisotropic distributions. The statistical approach is based on the spherical harmonic expansion of angular correlations.
A new method for analyzing the morphological features of point patterns is presented. The method is taken from the study of molecular liquids, where it has been introduced for making a statistical description of anisotropic distributions. The statistical approach is based on the spherical harmonic expansion of angular correlations.
△ Less
Submitted 6 November, 2001;
originally announced November 2001.
-
Numerical Convergence of Hydrodynamical SPH Simulations of Cooling Clusters
Authors:
R. Valdarnini
Abstract:
The results from hydrodynamical TREESPH simulations of galaxy clusters are used to investigate the dependence of the final cluster X-ray properties upon the numerical resolution and the assumed star formation models for the cooled gas. A comparison between runs with different star formation methods shows that the results of simulations, based on star formation methods in which gas conversion int…
▽ More
The results from hydrodynamical TREESPH simulations of galaxy clusters are used to investigate the dependence of the final cluster X-ray properties upon the numerical resolution and the assumed star formation models for the cooled gas. A comparison between runs with different star formation methods shows that the results of simulations, based on star formation methods in which gas conversion into stars is controlled by an efficiency parameter c_{star}, are sensitive to the simulation numerical resolution. In this respect star formation methods based instead on a local density threshold, are shown to give more stable results. Final X-ray luminosities are found to be numerically stable, with uncertainties of a factor 2.
△ Less
Submitted 5 November, 2001;
originally announced November 2001.
-
Numerical convergence of physical variables in hydrodynamical simulations of cooling clusters
Authors:
R. Valdarnini
Abstract:
Results from hydrodynamical SPH simulations of galaxy clusters are used to investigate the dependence of the final cluster X-ray properties on the numerical resolution and the assumed models for the physical gas processes. Two different spatially flat cosmological models have been considered: a low-density cold dark matter universe with a vacuum energy density Omega_L=0.7 (LCDM) and a cold+hot d…
▽ More
Results from hydrodynamical SPH simulations of galaxy clusters are used to investigate the dependence of the final cluster X-ray properties on the numerical resolution and the assumed models for the physical gas processes. Two different spatially flat cosmological models have been considered: a low-density cold dark matter universe with a vacuum energy density Omega_L=0.7 (LCDM) and a cold+hot dark matter model (CHDM). These simulations first include radiative cooling and then also conversion of cold gas particles into stars. When cold gas particles are allowed to convert into stars the final gas profiles show a well defined core radius and the temperature profiles are nearly flat. For the most massive test cluster in the LCDM model, these simulations show a prominent cooling flow in the cluster core. This cluster was analyzed in detail, running simulations with different star formation methods and increasing numerical resolution. A comparison between different runs shows that the results of simulations, based on star formation methods in which gas conversion into stars is controlled by an efficiency parameter c_star, are sensitive to the numerical resolution of the simulation. In this respect star formation methods based instead on a local density threshold, as in Navarro and White (1993), are shown to give more stable results. Final X-ray luminosities are found to be numerically stable, with uncertainties of a factor two.
△ Less
Submitted 25 October, 2001;
originally announced October 2001.
-
The morphological and dynamical evolution of simulated galaxy clusters
Authors:
C. Beisbart,
R. Valdarnini,
T. Buchert
Abstract:
We explore the morphological and dynamical evolution of galaxy clusters in simulations using scalar and vector-valued Minkowski valuations and the concept of fundamental plane relations. In this context, three questions are of fundamental interest: 1. How does the average cluster morphology depend on the cosmological background model? 2. Is it possible to discriminate between different cosmologi…
▽ More
We explore the morphological and dynamical evolution of galaxy clusters in simulations using scalar and vector-valued Minkowski valuations and the concept of fundamental plane relations. In this context, three questions are of fundamental interest: 1. How does the average cluster morphology depend on the cosmological background model? 2. Is it possible to discriminate between different cosmological models using cluster substructure in a statistically significant way? 3. How is the dynamical state of a cluster, especially its distance from a virial equilibrium, correlated to its visual substructure? To answer these questions, we quantify cluster substructure using a set of morphological order parameters constructed on the basis of the Minkowski valuations (MVs). The dynamical state of a cluster is described using global cluster parameters: in certain spaces of such parameters fundamental band-like structures are forming indicating the emergence of a virial equilibrium. We find that the average distances from these fundamental structures are correlated to the average amount of cluster substructure for our cluster samples during the time evolution. Furthermore, significant differences show up between the high- and the low-Omega models. We pay special attention to the redshift evolution of morphological characteristics and find large differences between the cosmological models even for higher redshifts.
△ Less
Submitted 25 September, 2001;
originally announced September 2001.
-
Detection of nonrandom patterns in cosmological gravitational clustering
Authors:
R. Valdarnini
Abstract:
A new method for analyzing point patterns produced by the evolution of gravitational clustering is presented. The method is taken from the study of molecular liquids, where it has been introduced for making a statistical description of anisotropic distributions. The statistical approach is based on the spherical harmonic expansion of angular correlations. A general introduction to the method is…
▽ More
A new method for analyzing point patterns produced by the evolution of gravitational clustering is presented. The method is taken from the study of molecular liquids, where it has been introduced for making a statistical description of anisotropic distributions. The statistical approach is based on the spherical harmonic expansion of angular correlations. A general introduction to the method is given; a theoretical analysis shows that it is partially connected with previous harmonic analyses applied to galaxy catalogs. The effectiveness of the statistical analysis in quantifying clustering morphology is illustrated by applying the statistical estimators to point distributions produced by an ensemble of cosmological $N-$ body simulations with a CDM spectrum. The results demonstrate that the statistical method is able to detect anisotropies in non-random patterns, with different scales being probed according to the expansion coefficients.
△ Less
Submitted 16 November, 2000;
originally announced November 2000.
-
Morphological evolution of X-ray clusters using hydrodynamical simulations
Authors:
R. Valdarnini,
S. Ghizzardi,
S. Bonometto
Abstract:
A large set of TREESPH simulations is used to test the global morphology of galaxy clusters and its evolution against X-ray data. A powerful method to investigate substructures in galaxy clusters are the power ratios introduced by Buote & Tsai. We consider three flat cosmological models: CDM, LCDM($Ω_Λ= 0.7$) and CHDM($Ω_h = 0.2$, 1 massive $ν$), all normalized so to fit the observed number of c…
▽ More
A large set of TREESPH simulations is used to test the global morphology of galaxy clusters and its evolution against X-ray data. A powerful method to investigate substructures in galaxy clusters are the power ratios introduced by Buote & Tsai. We consider three flat cosmological models: CDM, LCDM($Ω_Λ= 0.7$) and CHDM($Ω_h = 0.2$, 1 massive $ν$), all normalized so to fit the observed number of clusters. For each model we built 40 clusters, using a TREESPH code, and performed a statistical comparison with a data sample including nearby clusters observed with ROSAT PSPC instrument. The comparison disfavors the LCDM model, as clusters appear too relaxed, while CDM and CHDM clusters, in which a higher degree of complexity occurs, seem to be closer to observations. A better fit of data can be expected for some different DM mix. If DM distributions are used instead of baryons, we find substructures more pronounced than in gas and models have a different score. Using hydrodynamical simulations is therefore essential to our aims.
△ Less
Submitted 4 January, 2000;
originally announced January 2000.
-
Global cluster morphology and its evolution: X-ray data vs CDM, LCDM and mixed models
Authors:
R. Valdarnini,
S. Ghizzardi,
S. A. Bonometto
Abstract:
The global structure of galaxy clusters and its evolution are tested within a large set of TREESPH simulations, so to allow a fair statistical comparison with available X-ray data. Structure tests are based on the "power ratios", introduced by Buote & Tsai. Cosmological models considered are CDM, LCDM (Omega_L=0.7) and CHDM (1 mass.neu., Omega_h = 0.2). All models are normalized to provide a fai…
▽ More
The global structure of galaxy clusters and its evolution are tested within a large set of TREESPH simulations, so to allow a fair statistical comparison with available X-ray data. Structure tests are based on the "power ratios", introduced by Buote & Tsai. Cosmological models considered are CDM, LCDM (Omega_L=0.7) and CHDM (1 mass.neu., Omega_h = 0.2). All models are normalized to provide a fair number density of clusters. For each model we run a P3M simulation in a large box, where we select the most massive 40 clusters. Going back to the initial redshift we run a hydro-TREESPH simulation for each of them. In this way we perform a statistical comparison of the global morphology of clusters, for each cosmological model, with ROSAT data, using Student t-test, F-test and K-S test. The last test and its generalization to 2--D distributions are also used to compare the joint distributions of 2 or 3 power ratios. We find that, using DM distribution, instead of gas, as done by some authors, leads to biased results, as baryons are distributed in a less structured way than DM. We also find that the cosmological models considered have different behaviours in these tests: LCDM has the worst performance. CDM and our CHDM have similar scores. The general trend of power ratio distributions is already fit by these models, but a further improvement is expected either from a different DM mix or a non-flat CDM model.
△ Less
Submitted 12 January, 1999;
originally announced January 1999.
-
Large scale structure formation in mixed dark matter models with a cosmological constant
Authors:
R. Valdarnini,
T. Kahniashvili,
B. Novosyadlyj
Abstract:
We study linear power spectra and formation of large scale structures in flat cosmological models with $Λ\ge 0$ and cold plus hot dark matter components (MLM). The hot component consists of massive neutrinos with cosmological density $Ω_H$. The linearized Einstein-Boltzmann equations for the evolution of the metric and density perturbations are integrated for a set of values of the cosmological…
▽ More
We study linear power spectra and formation of large scale structures in flat cosmological models with $Λ\ge 0$ and cold plus hot dark matter components (MLM). The hot component consists of massive neutrinos with cosmological density $Ω_H$. The linearized Einstein-Boltzmann equations for the evolution of the metric and density perturbations are integrated for a set of values of the cosmological parameters. For all the considered models we assume a scale-invariant primeval spectrum. The density weighted final linear power spectra are normalized to the four year COBE data and have been used to constrain the parameter space by a comparison of linear predictions with the current observational data on large scales. The consistency of MLM predictions with the observable data set is best obtained for models with one specie of massive neutrinos and $Ω_H/Ω_M \le0.2$. For this class of MLM models we obtain constraints from linear data on the present matter density. Consistency with the estimated cluster abundance can be achieved for COBE normalized MLM models with $Ω_H/Ω_M\le0.2$ and $ 0.45 \le Ω_M \le 0.75$ for $h=0.5$. If $h=0.7$ then $ 0.3 \le Ω_M \le 0.5$. These constraints are at $1σ$ level and standard MDM models are clearly ruled out. We note that the range of allowed values for $Ω_M$, that we obtain for MLM models from linear analysis, is also approximately the same range that is needed in order to consistently satisfy a variety of independent observational constraints.
△ Less
Submitted 6 April, 1998;
originally announced April 1998.
-
Evolution of Cluster Morphology: X-ray data vs. CDM models
Authors:
R. Valdarnini,
S. Ghizzardi,
S. Bonometto
Abstract:
We compare the evolution of morphology in X-ray cluster data and in clusters obtained in a simulation of flat CDM, normalized to the observed cluster abundance. We find that the evolution rate in model clusters is significantly higher than in data clusters. The results is stricking, as our cluster data are all just for $z \simlt 0.2$, but is not contrary to expectations, as is known that the CDM…
▽ More
We compare the evolution of morphology in X-ray cluster data and in clusters obtained in a simulation of flat CDM, normalized to the observed cluster abundance. We find that the evolution rate in model clusters is significantly higher than in data clusters. The results is stricking, as our cluster data are all just for $z \simlt 0.2$, but is not contrary to expectations, as is known that the CDM spectrum is too flat, around the cluster scale, and therefore induces the presence of too many substructures. The cluster simulations were run using a TREESPH code and therefore include hydrodynamical effects. The test we performed, which turns out to be very sensitive, starts from the so--called power ratios introduced by Buote and Tsai, and makes use of the 2--dimensional generalization of the Kolmogorov--Smirnov test. The discrepancy between data and model is however directly visible through linear regression techniques and using simpler statistical tests.
△ Less
Submitted 24 February, 1998;
originally announced February 1998.
-
A global descriptor of spatial pattern interaction in the galaxy distribution
Authors:
Martin Kerscher,
Maria Jesus Pons-Borderia,
Jens Schmalzing,
Roberto Trasarti-Battistoni,
Thomas Buchert,
Vicent J. Martinez,
Riccardo Valdarnini
Abstract:
We present the function J as a morphological descriptor for point patterns formed by the distribution of galaxies in the Universe. This function was recently introduced in the field of spatial statistics, and is based on the nearest neighbor distribution and the void probability function. The J descriptor allows to distinguish clustered (i.e. correlated) from ``regular'' (i.e. anti-correlated) p…
▽ More
We present the function J as a morphological descriptor for point patterns formed by the distribution of galaxies in the Universe. This function was recently introduced in the field of spatial statistics, and is based on the nearest neighbor distribution and the void probability function. The J descriptor allows to distinguish clustered (i.e. correlated) from ``regular'' (i.e. anti-correlated) point distributions. We outline the theoretical foundations of the method, perform tests with a Matern cluster process as an idealised model of galaxy clustering, and apply the descriptor to galaxies and loose groups in the Perseus-Pisces Survey. A comparison with mock-samples extracted from a mixed dark matter simulation shows that the J descriptor can be profitably used to constrain (in this case reject) viable models of cosmic structure formation.
△ Less
Submitted 28 September, 1998; v1 submitted 6 December, 1997;
originally announced December 1997.
-
Three-point correlation function of galaxy clusters in cosmological models: a strong dependence on triangle shape
Authors:
Y. P. Jing,
G. Boerner,
R. Valdarnini
Abstract:
In this paper, we use large $\pppm$ N-body simulations to study the three-point correlation function $ζ$ of clusters in two theoretical models. The first model (LCDM) is a low-density flat model of $Ω_0=0.3$, $Λ_0=0.7$ and $h=0.75$, and the second model (PIM) is an Einstein-de-Sitter model with its linear power spectrum obtained from observations. We found that the scaled function $Q(r,u,v)$, wh…
▽ More
In this paper, we use large $\pppm$ N-body simulations to study the three-point correlation function $ζ$ of clusters in two theoretical models. The first model (LCDM) is a low-density flat model of $Ω_0=0.3$, $Λ_0=0.7$ and $h=0.75$, and the second model (PIM) is an Einstein-de-Sitter model with its linear power spectrum obtained from observations. We found that the scaled function $Q(r,u,v)$, which is defined as the ratio of $ζ(r, ru, ru+rv)$ to the hierarchical sum $ξ(r)ξ(ru)+ ξ(ru) ξ(ru+rv) +ξ(ru+rv)ξ(r)$ (where $ξ$ is the two-point correlation function of clusters), depends weakly on $r$ and $u$, but very strongly on $v$. $Q(r,u,v)$ is about 0.2 at $v=0.1$ and 1.8 at $v=0.9$. A model of $Q(r,u,v)=Θ10^{1.3v^2}$ can fit the data of $ζ$ very nicely with $Θ\approx 0.14$. This model is found to be universal for the LCDM clusters and for the PIM clusters. Furthermore, $Q(r,u,v)$ is found to be insensitive to the cluster richness. We have compared our N-body results with simple analytical theories of cluster formation, like the peak theories or the local maxima theories. We found that these theories do not provide an adequate description for the three-point function of clusters. We have also examined the observational data of $ζ$ presently available, and have not found any contradiction between the observations and our model predictions. The $v$-dependence of $q$ in a projected catalogue of clusters is shown to be much weaker than the $v$-dependence of $Q$ found in the three-dimensional case. It would be important to search for the $v$-dependence of $Q$ in redshift samples of rich clusters.
△ Less
Submitted 28 June, 1995;
originally announced June 1995.
-
The 1-Point Cluster Distribution Function and its Moments
Authors:
M. Plionis,
R. Valdarnini
Abstract:
We derive the 1-point probability density function of the smoothed 3-D Abell-ACO cluster density field and we compare it with that of artificial cluster samples, generated as high peaks of a Gaussian field in such a way that they reproduce the low-order (2- and 3-point) correlation functions and the observed cluster selection functions. We find that both real and simulated {\em pdf}'s are well a…
▽ More
We derive the 1-point probability density function of the smoothed 3-D Abell-ACO cluster density field and we compare it with that of artificial cluster samples, generated as high peaks of a Gaussian field in such a way that they reproduce the low-order (2- and 3-point) correlation functions and the observed cluster selection functions. We find that both real and simulated {\em pdf}'s are well approximated by a log-normal distribution even when the Gaussian smoothing radius is as large as 40 $h^{-1}$ Mpc. Furthermore the low-order moments of the {\em pdf} are found to obey a relation $γ= S_{3} σ^4$, with $γ$ being the skewness and $S_{3}\approx 1.8$. These results are consistent with clusters being high-peaks of an underlying initial Gaussian density field. A by-product of our analysis is that when we rescale the {\em pdf} cluster moments to those of the QDOT-IRAS galaxies, using linear biasing with $b_{cI}\sim 4.5$ and for the common smoothing radius of 20 $h^{-1}$ Mpc, we find them to be significantly smaller than those directly estimated from the QDOT data by Saunders et al. (1991).
△ Less
Submitted 2 September, 1994;
originally announced September 1994.
-
Is there any Scaling in the Cluster Distribution ?
Authors:
S. Borgani,
V. J. Martinez,
M. A. Perez,
R. Valdarnini
Abstract:
We apply fractal analysis methods to investigate the scaling properties in the Abell and ACO catalogs of rich galaxy clusters. We also discuss different technical aspects of the method when applied to data sets with small number of points as the cluster catalogs. Results are compared with simulations based on the Zel'dovich approximation. We limit our analysis to scales less than 100 $\hm$. The…
▽ More
We apply fractal analysis methods to investigate the scaling properties in the Abell and ACO catalogs of rich galaxy clusters. We also discuss different technical aspects of the method when applied to data sets with small number of points as the cluster catalogs. Results are compared with simulations based on the Zel'dovich approximation. We limit our analysis to scales less than 100 $\hm$. The cluster distribution show a scale invariant multifractal behavior in a limited scale range. For the Abell catalog this range is 15--60$\hm$, while for the ACO sample it extends to smaller scales. Despite this difference in the extension of the scale--range where scale--invariant clustering takes place, both samples are characterized by remarkably similar multifractal spectra in the corresponding scaling regime. In particular, the correlation dimension turns out to be $D_2\simeq 2.2$ for both Abell and ACO clusters. Although it is difficult to point out the scale at which homogeneity is reached with the present size of these redshift surveys, our results indicate that the cluster distribution shows a tendency to homogeneity at large scales, disproving the picture of a pure scale invariant fractal structure extending to arbitrarily large distances.
△ Less
Submitted 7 March, 1994;
originally announced March 1994.
-
Modelling the Power Spectrum of Density Fluctuations: A Phenomenological Approach
Authors:
Enzo Branchini,
Luigi Guzzo,
Riccardo Valdarnini
Abstract:
We show how, based on considerations on the observed form of the galaxy 2-point spatial correlation function xi(r), a very simplified -- yet surprisingly effective -- model for the linear density fluctuations power spectrum can be constructed. We first relate the observed large-scale shape of xi(r) to a power-law form for the power spectrum, P(k)\propto k^{-2.2}. For a plausible value of the bia…
▽ More
We show how, based on considerations on the observed form of the galaxy 2-point spatial correlation function xi(r), a very simplified -- yet surprisingly effective -- model for the linear density fluctuations power spectrum can be constructed. We first relate the observed large-scale shape of xi(r) to a power-law form for the power spectrum, P(k)\propto k^{-2.2}. For a plausible value of the bias parameter b = 1/sigma_8 ~ 1.8, one has (delta_rho / rho)_{rms} ~ 1 r ~ 3.5/h Mpc, suggesting that the change of slope observed in xi(r) around this scale marks the transition between the linear and nonlinear gravitational regimes. Under this working hypothesis, we use a simple analytical form to fit the large-scale correlations constraints together with the COBE CMB anisotropy measurement, thus constructing a simple phenomenological model for the linear power spectrum. Despite its simplicity, the model fits remarkably well directly estimated power spectra from different optical galaxy samples, and when evolved through an N-body simulation it provides a good match to the observed galaxy correlations. One of the most interesting features of the model is the small-scale one-dimensional velocity dispersion produced: sigma_{1d} = 450 Km s^{-1} at 0.5/h Mpc and sigma_{1d} = 350 Km s^{-1} for separations larger than ~ 2/h Mpc.
△ Less
Submitted 12 January, 1994;
originally announced January 1994.