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.

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.

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.

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.

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.

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.

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.

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.

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 \Delta=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}(\Delta /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.

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.

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]

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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($\Omega_\Lambda = 0.7$) and CHDM($\Omega_h = 0.2$, 1 massive $\nu$), 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.

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.

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.

We study linear power spectra and formation of large scale structures in flat cosmological models with $\Lambda \ge 0$ and cold plus hot dark matter components (MLM). The hot component consists of massive neutrinos with cosmological density $\Omega_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 $\Omega_H/\Omega_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 $\Omega_H/\Omega_M\le0.2$ and $ 0.45 \le \Omega_M \le 0.75$ for $h=0.5$. If $h=0.7$ then $ 0.3 \le \Omega_M \le 0.5$. These constraints are at $1\sigma$ level and standard MDM models are clearly ruled out. We note that the range of allowed values for $\Omega_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.

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.

In this paper, we use large $\pppm$ N-body simulations to study the three-point correlation function $\zeta$ of clusters in two theoretical models. The first model (LCDM) is a low-density flat model of $\Omega_0=0.3$, $\Lambda_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 $\zeta (r, ru, ru+rv)$ to the hierarchical sum $\xi (r)\xi (ru)+ \xi (ru) \xi (ru+rv) +\xi (ru+rv)\xi (r)$ (where $\xi$ 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)=\Theta 10^{1.3v^2}$ can fit the data of $\zeta$ very nicely with $\Theta\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 $\zeta$ 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.

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 $\gamma = S_{3} \sigma^4$, with $\gamma$ 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).

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.

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.