Context: In our Paper I, by using statistical deconvolution methods, extended kinematics maps of Gaia-DR2 data have been produced in a range of heliocentric distances that are a factor of two to three larger than those analyzed previously by the Gaia Collaboration with the same data. It added the range of Galactocentric distances between 13 kpc and 20 kpc to the previous maps. Aims: Here, we investigate the dynamical effects produced by different mechanisms that can explain the radial and vertical components of these extended kinematic maps, including a decomposition of bending and breathing of the vertical components. This paper as a whole tries to be a compendium of different dynamical mechanisms whose predictions can be compared to the kinematic maps. Methods: Using analytical methods or simulations, we are able to predict the main dynamical factors and compare them to the predictions of the extended kinematic maps of Gaia-DR2. Results: The gravitational influence of Galactic components that are different from the disk, such as the long bar or bulge, the spiral arms, or a tidal interaction with Sagittarius dwarf galaxy, may explain some features of the velocity maps, especially in the inner parts of the disk. However, they are not sufficient in explaining the most conspicuous gradients in the outer disk. Vertical motions might be dominated by external perturbations or mergers, although a minor component may be due to a warp whose amplitude evolves with time. Here, we show with two different methods, which analyze the dispersion of velocities, that the mass distribution of the disk is flared. Despite these partial explanations, the main observed features can only be explained in terms of out-of-equilibrium models, which are either due to external perturbers or to the fact that the disk has not had time to reach equilibrium since its formation.

The evolution of self-gravitating systems, and long-range interacting systems more generally, from initial configurations far from dynamical equilibrium is often described as a simple two phase process: a first phase of violent relaxation bringing it to a quasi-stationary state in a few dynamical times, followed by a slow adiabatic evolution driven by collisional processes. In this context the complex spatial structure evident, e.g., in spiral galaxies is understood either in terms of instabilities of quasi-stationary states, or a result of dissipative non-gravitational interactions. We illustrate here, using numerical simulations, that purely self-gravitating systems evolving from quite simple initial configurations can in fact give rise easily to structures of this kind of which the lifetime can be large compared to the dynamical characteristic time, but short compared to the collisional relaxation time scale. More specifically, for a broad range of non-spherical and non-uniform rotating initial conditions, gravitational relaxation gives rise quite generically to long-lived non-stationary structures of a rich variety, characterized by spiral-like arms, bars and even ring-like structures in special cases. These structures are a feature of the intrinsically out-of-equilibrium nature of the system's collapse, associated with a part of the system's mass while the bulk is well virialized. They are characterized by predominantly radial motions in their outermost parts, but also incorporate an extended flattened region which rotates coherently about a well virialized core of triaxial shape with an approximately isotropic velocity dispersion. We discuss the possible relevance of these simple toy models to the observed structure of real galaxies emphasizing the difference between dissipative and dissipationless disc formation.

We measure the mean Galactocentric radial component of the velocity of stars ($v_R$) in the disk at 8 kpc$<R<28$ kpc in the direction of the anticenter. For this, we use the Apache Point Galactic Evolution Experiment (APOGEE). Furthermore, we compare the result with HI maps along the same line of sight. We find an increase in positive (expansion) $v_R$ at $R\approx 9-13$ kpc, reaching a maximum of $\approx 6$ km/s, and a decrease at large values of $R$ reaching a negative (contraction) value of $\approx -10$ km/s for $R>17$ kpc. Negative velocities are also observed in 21 cm HI maps, possibly dominated by local gas emission. Among the possible dynamical causes for these non-zero $v_R$, factors such as the effect of the Galactic bar, streams, or mergers do not seem appropriate to explain our observations. An explanation might be the gravitational attraction of overdensities in a spiral arm. As a matter of fact, we see a change of regime from positive to negative velocities around $R\approx 15$ kpc, in the position where we cross the Outer spiral arm in the anticenter. The mass in spiral arms necessary to produce these velocities would be about 3\% of the mass of the disk, consistent with our knowledge of the spiral arms. Another scenario that we explore is a simple class of out-of-equilibrium systems in which radial motions are generally created by the monolithic collapse of isolated self-gravitating overdensities.

Galaxy velocity maps often show the typical pattern of a rotating disk, consistent with the dynamical model where emitters rotate in circular orbits around the galactic center. The simplest template used to fit these maps consists in the rotating disk model (RDM) where the amplitude of circular velocities is fixed by the observed velocity profile along the kinematic axis. A more sophisticated template is the rotating tilted-ring model (RTRM) that takes into account the presence of warps and allows a radius-dependent orientation of the kinematic axis. In both cases, axisymmetry is assumed and residuals between the observed and the model velocity fields are interpreted as noncircular motions. We show that if a galaxy is not axisymmetric, there is an intrinsic degeneracy between a rotational and a radial velocity field. We then introduce a new galaxy template, the radial ellipse model (REM), that is not axisymmetric and has a purely radial velocity field with an amplitude that is correlated with the major axis of the ellipse. We show that best fits to the observed two-dimensional velocity fields of 28 galaxies extracted from the THINGS sample with both the REM and the RDM give residuals with similar amplitudes, where the REM residuals trace nonradial motions. Best fits obtained with the RTRM, because of its larger number of free parameters, give the smallest residuals: however, we argue that this does not necessarily imply that the RTRM gives the most accurate representation of a galaxy velocity field. Instead, we show that this method is not able to disentangle between circular and radial motions for the case of nonaxisymmetric systems. Finally, we consider the physical motivation of the REM, and discuss how the interpretation of galactic dynamics changes if one assumes that the main component of a galaxy velocity field is modeled as a RDM/RTRM or as a REM.

CONTEXT. The Gaia Collaboration has used Gaia-DR2 sources with six-dimensional (6D) phase space information to derive kinematical maps within 5 kpc of the Sun, which is a reachable range for stars with relative error in distance lower than 20%. AIMS. Here we aim to extend the range of distances by a factor of two to three, thus adding the range of Galactocentric distances between 13 kpc and 20 kpc to the previous maps, with their corresponding error and root mean square values. METHODS. We make use of the whole sample of stars of Gaia-DR2 including radial velocity measurements, which consists in more than seven million sources, and we apply a statistical deconvolution of the parallax errors based on the Lucy's inversion method of the Fredholm integral equations of the first kind, without assuming any prior. RESULTS. The new extended maps provide lots of new and corroborated information about the disk kinematics: significant departures of circularity in the mean orbits with radial Galactocentric velocities between -20 and +20 km/s and vertical velocities between -10 and +10 km/s; variations of the azimuthal velocity with position; asymmetries between the northern and the southern Galactic hemispheres, especially towards the anticenter that includes a larger azimuthal velocity in the south; and others. CONCLUSIONS. These extended kinematical maps can be used to investigate the different dynamical models of our Galaxy, and we will present our own analyses in the forthcoming second part of this paper. At present, it is evident that the Milky Way is far from a simple stationary configuration in rotational equilibrium, but is characterized by streaming motions in all velocity components with conspicuous velocity gradients.

In this paper we attempt to answer to the question: can cosmic acceleration of the Universe have a fractal solution? We give an exact solution of a Lema\^itre-Tolman-Bondi (LTB) Universe based on the assumption that such a smooth metric is able to describe, on average, a fractal distribution of matter. While the LTB model has a center, we speculate that, when the fractal dimension is not very different from the space dimension, this metric applies to any point of the fractal structure when chosen as center so that, on average, there is not any special point or direction. We examine the observed magnitude-redshift relation of type Ia supernovae (SNe Ia), showing that the apparent acceleration of the cosmic expansion can be explained as a consequence of the fractal distribution of matter when the corresponding space-time metric is modeled as a smooth LTB one and if the fractal dimension on scales of a few hundreds Mpc is $D=2.9 \pm 0.02$.

We describe how a simple class of out of equilibrium, rotating and asymmetrical mass distributions evolve under their self-gravity to produce a quasi-planar spiral structure surrounding a virialized core, qualitatively resembling a spiral galaxy. The spiral structure is transient, but can survive tens of dynamical times, and further reproduces qualitatively noted features of spiral galaxies as the predominance of trailing two-armed spirals and large pitch angles. As our models are highly idealized, a detailed comparison with observations is not appropriate, but generic features of the velocity distributions can be identified to be potential observational signatures of such a mechanism. Indeed, the mechanism leads generically to a characteristic transition from predominantly rotational motion, in a region outside the core, to radial ballistic motion in the outermost parts. Such radial motions are excluded in our Galaxy up to 15 kpc, but could be detected at larger scales in the future by GAIA. We explore the apparent motions seen by external observers of the velocity distributions of our toy galaxies, and find that it is difficult to distinguish them from those of a rotating disc with sub-dominant radial motions at levels typically inferred from observations. These simple models illustrate the possibility that the observed apparent motions of spiral galaxies might be explained by non-trivial non-stationary mass and velocity distributions without invoking a dark matter halo or modification of Newtonian gravity. In this scenario the observed phenomenological relation between the centripetal and gravitational acceleration of the visible baryonic mass could have a simple explanation.

The determination of the resolution of cosmological N-body simulations, i.e., the range of scales in which quantities measured in them represent accurately the continuum limit, is an important open question. We address it here using scale-free models, for which self-similarity provides a powerful tool to control resolution. Such models also provide a robust testing ground for the so-called stable clustering approximation, which gives simple predictions for them. Studying large N-body simulations of such models with different force smoothing, we find that these two issues are in fact very closely related: our conclusion is that the accuracy of two point statistics in the non-linear regime starts to degrade strongly around the scale at which their behaviour deviates from that predicted by the stable clustering hypothesis. Physically the association of the two scales is in fact simple to understand: stable clustering fails to be a good approximation when there are strong interactions of structures (in particular merging) and it is precisely such non-linear processes which are sensitive to fluctuations at the smaller scales affected by discretisation. Resolution may be further degraded if the short distance gravitational smoothing scale is larger than the scale to which stable clustering can propagate. We examine in detail the very different conclusions of studies by Smith et al. (2003) and Widrow et al. (2009} and find that the strong deviations from stable clustering reported by these works are the results of over-optimistic assumptions about scales resolved accurately by the measured power spectra, and the reliance on Fourier space analysis. We emphasise the much poorer resolution obtained with the power spectrum compared to the two point correlation function..

Simulations of purely self-gravitating N-body systems are often used in astrophysics and cosmology to study the collisionless limit of such systems. Their results for macroscopic quantities should then converge well for sufficiently large N. Using a study of the evolution from a simple space of spherical initial conditions - including a region characterised by so-called "radial orbit instability" - we illustrate that the values of N at which such convergence is obtained can vary enormously. In the family of initial conditions we study, good convergence can be obtained up to a few dynamical times with N $ \sim 10^3$ - just large enough to suppress two body relaxation - for certain initial conditions, while in other cases such convergence is not attained at this time even in our largest simulations with N $\sim 10^5$. The qualitative difference is due to the stability properties of fluctuations introduced by the N-body discretisation, of which the initial amplitude depends on N. We discuss briefly why the crucial role which such fluctuations can potentially play in the evolution of the N-body system could, in particular, constitute a serious problem in cosmological simulations of dark matter.

We study the collapse of an isolated, initially cold, irregular (but almost spherical) and (slightly) inhomogeneous cloud of self-gravitating particles. The cloud is driven towards a virialized quasi-equilibrium state by a fast relaxation mechanism, occurring in a typical time $\tau_c$, whose signature is a large change in the particle energy distribution. Post-collapse particles are divided into two main species: bound and free, the latter being ejected from the system. Because of the initial system's anisotropy, the time varying gravitational field breaks spherical symmetry so that the ejected mass can carry away angular momentum and the bound system can gain a non-zero angular momentum. In addition, while strongly bound particles form a compact core, weakly bound ones may form, in a time scale of the order of $\tau_c$, several satellite sub-structures. These satellites have a finite lifetime that can be longer than $\tau_c$ and generally form a flattened distribution. Their origin and their abundance are related to the amplitude and nature of initial density fluctuations and to the initial cloud deviations from spherical symmetry, which are both amplified during the collapse phase. Satellites show a time dependent virial ratio that can be different from the equilibrium value $b\approx -1$: although they are bound to the main virialized object, they are not necessarily virially relaxed.

During the violent relaxation of a self-gravitating system a significant fraction of its mass may be ejected. If the time varying gravitational field also breaks spherical symmetry this mass can potentially carry angular momentum. Thus starting initial configurations with zero angular momentum can in principle lead to a bound virialized system with non-zero angular momentum. We explore here, using numerical simulations, how much angular momentum can be generated in a virialized structure in this way, starting from configurations of cold particles which are very close to spherically symmetric. For initial configurations in which spherical symmetry is broken only by the Poissonian fluctuations associated with the finite particle number $N$, with $N$ in range $10^3$ to $10^5$, we find that the relaxed structures have standard "spin" parameters $\lambda \sim 10^{-3}$, and decreasing slowly with $N$. For slightly ellipsoidal initial conditions, in which the finite-$N$ fluctuations break the residual reflection symmetries, we observe values $\lambda \sim 10^{-2}$, of the same order of magnitude as those reported for elliptical galaxies. The net angular momentum vector is typically aligned close to normal to the major semi-axis of the triaxial relaxed structure, and also with that of the ejected mass. This simple mechanism may provide an alternative, or complement, to "tidal torque theory" for understanding the origin of angular momentum in astrophysical structures.

An isolated, initially cold and ellipsoidal cloud of self-gravitating particles represents a relatively simple system to study the effects of the deviations from spherical symmetry in the mechanism of violent relaxation. Initial deviations from spherical symmetry are shown to play a dynamical role that is equivalent to that of density fluctuations in the case of an initially spherical cloud. Indeed, these deviations control the amount of particles energy change and thus determine the properties of the final energy distribution, particularly the appearance of two species of particles: bound and free. Ejection of mass and energy from the system together with the formation of a density profile decaying as $\rho(r) \sim r^{-4}$ and a Keplerian radial velocity dispersion profile, are the prominent features similar to those observed after the violent relaxation of spherical clouds. In addition, we find that ejected particles are characterized by highly non-spherical shapes, whose features can be traced in the initial deviations from spherical symmetry that are amplified during the dynamical evolution: particles can indeed form anisotropic configurations, like bars and/or disks, even though the initial cloud was very close to spherical.

Initially cold and spherically symmetric self-gravitating systems may give rise to a virial equilibrium state which is far from spherically symmetric, and typically triaxial. We focus here on how the degree of symmetry breaking in the final state depends on the initial density profile. We note that the most asymmetric structures result when, during the collapse phase, there is a strong injection of energy preferentially into the particles which are localized initially in the outer shells. These particles are still collapsing when the others, initially located in the inner part, are already re-expanding; the motion of particles in a time varying potential allow them to gain kinetic energy --- in some cases enough to be ejected from the system. We show that this mechanism of energy gain amplifies the initial small deviations from perfect spherical symmetry due to finite $N$ fluctuations. This amplification is more efficient when the initial density profile depends on radius, because particles have a greater spread of fall times compared to a uniform density profile, for which very close to symmetric final states are obtained}. These effects lead to a distinctive correlation of the orientation of the final structure with the distribution of ejected mass, and also with the initial (very small) angular fluctuations.

One of the main problems of observational cosmology is to determine the range in which a reliable measurement of galaxy correlations is possible. This corresponds to determine the shape of the correlation function, its possible evolution with redshift and the size and amplitude of large scale structures. Different selection effects, inevitably entering in any observation, introduce important constraints in the measurement of correlations. In the context of galaxy redshift surveys selection effects can be caused by observational techniques and strategies and by implicit assumptions used in the data analysis. Generally all these effects are taken into account by using pair-counting algorithms to measure two-point correlations. We review these methods stressing that they are based on the a-priori assumption that galaxy distribution is spatially homogeneous inside a given sample. We show that, when this assumption is not satisfied by the data, results of the correlation analysis are affected by finite size effects.In order to quantify these effects, we introduce a new method based on the computation of the gradient of galaxy counts along tiny cylinders. We show, by using artificial homogeneous and inhomogeneous point distributions, that this method is to identify redshift dependent selection effects and to disentangle them from the presence of large scale density fluctuations. We then apply this new method to several redshift catalogs and we find evidences that galaxy distribution, in those samples where selection effects are small enough, is characterized by power-law correlations with exponent $\gamma=0.9$ up to $20$ Mpc/h followed by a change of slope that, in the range [20,100] Mpc/h, corresponds to a power-law exponent $\gamma=0.25$. Whether a crossover to spatial unformity occurs at $\sim 100$ Mpc/h cannot be clarified by the present data.

The evolution of an isolated over-density represents a useful toy model to test the accuracy of a cosmological N-body code in the non linear regime as it is approximately equivalent to that of a truly isolated cloud of particles, with same density profile and velocity distribution, in a non expanding background. This is the case as long as the system size is smaller than the simulation box side, so that its interaction with the infinite copies can be neglected. In such a situation, the over-density rapidly undergoes to a global collapse forming a quasi stationary state in virial equilibrium. However, by evolving the system with a cosmological code (GADGET) for a sufficiently long time, a clear deviation from such quasi-equilibrium configuration is observed. This occurs in a time t_{LI} that depends on the values of the simulation numerical parameters such as the softening length and thetime-stepping accuracy, i.e. it is a numerical artifact related to the limited spatial and temporal resolutions. The analysis of the Layzer-Irvine cosmic energy equation, confirms that this deviation corresponds to an unphysical dynamical regime. By varying the numerical parameters of the simulation and the physical parameters of the system we show that the unphysical behaviour originates from badly integrated close scatterings of high velocity particles. We find that, while the structure may remain virialized in the unphysical regime, its density and velocity profiles are modified with respect to the quasi-equilibrium configurations, converging however to well defined shapes, the former characterised by a Navarro Frenk White-type behaviour.

Beyond convergence studies and comparison of different codes, there are essentially no controls on the accuracy in the non-linear regime of cosmological N body simulations, even in the dissipationless limit. We propose and explore here a simple test which has not been previously employed: when cosmological codes are used to simulate an isolated overdensity, they should reproduce, in physical coordinates, those obtained in open boundary conditions without expansion. In particular, the desired collisionless nature of the simulations can be probed by testing for stability in physical coordinates of virialized equilibria. We investigate and illustrate the test using a suite of simulations in an Einstein de Sitter cosmology from initial conditions which rapidly settle to virial equilibrium. We find that the criterion of stable clustering allows one to determine, for given particle number N in the "halo" and force smoothing, a maximum red-shift range over which the collisionless limit may be represented with desired accuracy. We also compare our results to the so-called Layzer Irvine test, showing that it provides a weaker, but very useful, tool to constrain the choice of numerical parameters. Finally we outline in some detail how these methods could be employed to test the choice of the numerical parameters used in a cosmological simulation.

We study the collapse and virialization of an isolated spherical cloud of self-gravitating particles initially at rest and characterised by a power-law density profile, with exponent 0<= \alpha < 3, or by a Plummer, an Hernquist, a NFW, a Gaussian profile. We find that in all cases the virialized structure formed after the collapse has a density profile decaying, at large enough radii, as ~ r^{-4}, and a radial velocity dispersion profile decaying as ~r^{-1}. We show that these profiles originate from the physical mechanism responsible of the ejection of a fraction of cloud's mass and energy during the collapse and that this same mechanism washes out the dependence on the initial conditions. When a large enough initial velocity dispersion is given to the cloud particles, ejection does not occur anymore and consequently the virialized halo density and velocity profiles display features which reflect the initial conditions.

The collapse of an isolated, uniform and spherical cloud of self-gravitating particles represents a paradigmatic example of a relaxation process leading to the formation of a quasi-stationary state in virial equilibrium. We consider several N-body simulations of such a system, with the initial velocity dispersion as a free parameter. We show that there is a clear difference between structures formed when the initial virial ratio is b_0 =2K_0/W_0 < b_0^c ~ -1/2 and b_0> b_0^c. These two sets of initial conditions give rise respectively to a mild and violent relaxation occurring in about the same time scale: however in the latter case the system contracts by a large factor, while in the former it approximately maintains its original size. Correspondingly the resulting quasi equilibrium state is characterized by a density profile decaying at large enough distances as ~1/r^4 or with a sharp cut-off. The case b_0<b_0^c can be well described by the Lynden-Bell theory of collisionless relaxation considering the system confined in a box. On the other hand the relevant feature for b_0>b_0^c is the ejection of particles and energy, which is not captured by such a theoretical approach: for this case we introduce a simple physical model to explain the formation of the power-law density profile. This model shows that the behavior n(r) ~1/r^4 is the typical density profile that is obtained when the initial conditions are cold enough that mass and energy ejection occurs. In addition, we clarify the origin of the critical value of the initial virial ratio b_0^c.

We characterize galaxy correlations in the Sloan Digital Sky Survey by measuring several moments of galaxy counts in spheres. We firstly find that the average counts grows as a power-law function of the distance with an exponent D= 2.1+- 0.05 for r in [0.5,20] Mpc/h and D = 2.8+-0.05 for r in [30,150] Mpc/h. In order to estimate the systematic errors in these measurements we consider the counts variance finding that it shows systematic finite size effects which depend on the samples sizes. We clarify, by making specific tests, that these are due to galaxy long-range correlations extending up to the largest scales of the sample. The analysis of mock galaxy catalogs, generated from cosmological N-body simulations of the standard LCDM model, shows that for r<20 Mpc/h the counts exponent is D~2.0, weakly dependent on galaxy luminosity, while D=3 at larger scales. In addition, contrary to the case of the observed galaxy samples, no systematic finite size effects in the counts variance are found at large scales, a result that agrees with the absence of large scale, r~100 Mpc/h, correlations in the mock catalogs. We thus conclude that the observed galaxy distribution is characterized by correlations, fluctuations and hence structures, which are larger, both in amplitude and in spatial extension, than those predicted by the standard model LCDM of galaxy formation.

Standard models of galaxy formation predict that matter distribution is statistically homogeneous and isotropic and characterized by (i) spatial homogeneity for r<10 Mpc/h, (ii) small-amplitude structures of relatively limited size (i.e., r<100) Mpc/h and (iii) anti-correlations for r > r_c ~ 150 Mpc/h (i.e., no structures of size larger than r_c). Whether or not the observed galaxy distribution is interpreted to be compatible with these predictions depend on the a-priori assumptions encoded in the statistical methods employed to characterize the data and on the a-posteriori hypotheses made to interpret the results. We present strategies to test the most common assumptions and we find evidences that, in the available samples, galaxy distribution is spatially inhomogeneous for r<100 Mpc/h but statistically homogeneous and isotropic. We conclude that the observed inhomogeneities pose a fundamental challenge to the standard picture of cosmology but they also represent an important opportunity which may open new directions for many cosmological puzzles.

The large scale distribution of galaxies in the universe displays a complex pattern of clusters, super-clusters, filaments and voids with sizes limited only by the boundaries of the available samples. A quantitative statistical characterization of these structures shows that galaxy distribution is inhomogeneous in these samples, being characterized by large-amplitude fluctuations of large spatial extension. Over a large range of scales, both the average conditional density and its variance show a nontrivial scaling behavior: at small scales, r<20 Mpc/h, the average (conditional) density scales as 1/r. At larger scales, the density depends only weakly (logarithmically) on the system size and density fluctuations follow the Gumbel distribution of extreme value statistics. These complex behaviors are different from what is expected in a homogeneous distribution with Gaussian fluctuations. The observed density inhomogeneities pose a fundamental challenge to the standard picture of cosmology but it also represent an important opportunity which points to new directions with respect to many cosmological puzzles. Indeed, the fact that matter distribution is not uniform, in the limited range of scales sampled by observations, rises the question of understanding how inhomogeneities affect the large-scale dynamics of the universe. We discuss several attempts which try to model inhomogeneities in cosmology, considering their effects with respect to the role and abundance of dark energy and dark matter.

We study the statistical properties of the Luminous Red Galaxies sample from the Sloan Digital Sky Survey. In particular we test, by determining the probability density function (PDF) of galaxy (conditional) counts in spheres, whether statistical properties are self-averaging within the sample. We find that there are systematic differences in the shape of the PDF and in the location of its peak, signaling that there are major systematic effects in the data which make the estimation of volume average quantities unreliable within this sample. We discuss that these systematic effects are related to the fluctuating behavior of the redshift counts which can be originated by intrinsic fluctuations in the galaxy density field or by observational selection effects. The latter possibility implies that more than 20 % of the galaxies have not been observed and that such a selection should not be a smooth function of redshift.

We report results of a study of the Newtonian dynamics of N self-gravitating particles which start in a quasi-uniform spherical configuration, without initial velocities. These initial conditions would lead to a density singularity at the origin at a finite time when N \rightarrow \infty, but this singularity is regulated at any finite N (by the associated density fluctuations). While previous studies have focussed on the behaviour as a function of N of the minimal size reached during the contracting phase, we examine in particular the size and energy of the virialized halo which results. We find the unexpected result that the structure decreases in size as N increases, scaling in proportion to N^{-1/3}, a behaviour which is associated with an ejection of kinetic energy during violent relaxation which grows in proportion to N^{1/3}. This latter scaling may be qualitatively understood, and if it represents the asymptotic behaviour in N implies that this ejected energy is unbounded above. We discuss also tests we have performed which indicate that this ejection is a mean-field phenomenon (i.e. a result of collisionless dynamics).

We study the statistical properties of the gravitational field generated by galaxy distribution observed bythe Sloan Digital Sky Survey (DR7). We characterize the probability density function of gravitational force fluctuations and relate its limiting behaviors to the correlation properties of the underlying density field. In addition, we study whether the PDF converges to an asymptotic shape within sample volumes. We consider several volume-limited samples of the Sloan Digital Sky Survey and we compute the gravitational force probability density function (PDF). The gravitational force is computed in spheres of varying radius as is its PDF. We find that (i) the PDF of the force displays features that can be understood in terms of galaxy two-point correlations and (ii) density fluctuations on the largest scales probed, i.e. r~100 Mpc/h, still contribute significantly to the amplitude of the gravitational force. Our main conclusion is that fluctuations in the gravitational force field generated by galaxy structures are also relevant on scales ~ 100 Mpc/h. By assuming that the gravitational fluctuations in the galaxy distribution reflect those in the whole matter distribution, and that peculiar velocities and accelerations are simply correlated, we may conclude that large-scale fluctuations in the galaxy density field may be the source of the large-scale flows recently observed.

Cosmological density fields are assumed to be translational and rotational invariant, avoiding any special point or direction, thus satisfying the Copernican Principle. A spatially inhomogeneous matter distribution can be compatible with the Copernican Principle but not with the stronger version of it, the Cosmological Principle which requires the additional hypothesis of spatial homogeneity. We establish criteria for testing that a given density field, in a finite sample at low redshifts, is statistically and/or spatially homogeneous. The basic question to be considered is whether a distribution is, at different spatial scales, self-averaging. This can be achieved by studying the probability density function of conditional fluctuations. We find that galaxy structures in the SDSS samples, the largest currently available, are spatially inhomogeneous but statistically homogeneous and isotropic up to ~ 100 Mpc/h. Evidences for the breaking of self-averaging are found up to the largest scales probed by the SDSS data. The comparison between the results obtained in volumes of different size allows us to unambiguously conclude that the lack of elf-averaging is induced by finite-size effects due to long-range correlated fluctuations. We finally discuss the relevance of these results from the point of view of cosmological modeling.

AIMS. Several authors have claimed to detect a significant cross-correlation between microwave WMAP anisotropies and the SDSS galaxy distribution. We repeat these analyses to determine the different cross-correlation uncertainties caused by re-sampling errors and field-to-field fluctuations. The first type of error concerns overlapping sky regions, while the second type concerns non-overlapping sky regions. METHODS. To measure the re-sampling errors, we use bootstrap and jack-knife techniques. For the field-to-field fluctuations, we use three methods: 1) evaluation of the dispersion in the cross-correlation when correlating separated regions of WMAP with the original region of SDSS; 2) use of mock Monte Carlo WMAP maps; 3) a new method (developed in this article), which measures the error as a function of the integral of the product of the self-correlations for each map. RESULTS. The average cross-correlation for b>30 deg. is significantly stronger than the re-sampling errors--both the jack-knife and bootstrap techniques provide similar results--but it is of the order of the field-to-field fluctuations. This is confirmed by the cross-correlation between anisotropies and galaxies in more than the half of the sample being null within re-sampling errors. CONCLUSIONS. Re-sampling methods underestimate the errors. Field-to-field fluctuations dominate the detected signals. The ratio of signal to re-sampling errors is larger than unity in a way that strongly depends on the selected sky region. We therefore conclude that there is no evidence yet of a significant detection of the integrated Sachs-Wolfe (ISW) effect. Hence, the value of Omega_\Lambda ~ 0.8 obtained by the authors who assumed they were observing the ISW effect would appear to have originated from noise analysis

Super-homogeneity is a property that is supposed to be satisfied by matter fluctuations in all standard theoretical models of structure formation, such as LCDM and its variants. This is a global condition on the correlation properties of the matter density field, which can be understood as a consistency constraint in the framework of FRW cosmology, and it corresponds to a very fine tuned balance between negative and positive correlations of density fluctuations and to the fastest possible decay of the normalized mass variance on large scales. By considering several galaxy samples, we discuss that these are characterized by the presence of large amplitude fluctuations with spatial extension limited only the size of the current samples. There is therefore a tension between the standard prediction of super-homogeneity and the detection of large scale inhomogeneities in the matter distribution at scales of the order of 100 Mpc/h. We discuss the theoretical implications of these results with respect to models of structure formation and to future galaxy and CMBR data, emphasizing the central role of the super-homogeneity property in the current description of fluctuations in FRW models.

In order to investigate whether galaxy structures are compatible with the predictions of the standard LCDM cosmology, we focus here on the analysis of several simple and basic statistical properties of the galaxy density field. Namely, we test whether, on large enough scales (i.e., r>10 Mpc/h), this is self-averaging, uniform and characterized by a Gaussian probability density function of fluctuations. These are three different and clear predictions of the LCDM cosmology which are fulfilled in mock galaxy catalogs generated from cosmological N-body simulations representing this model. We consider some simple statistical measurements able to tests these properties in a finite sample. We discuss that the analysis of several samples of the Two Degree Field Galaxy Redshift Survey and of the Sloan Digital Sky Survey show that galaxy structures are non self-averaging and inhomogeneous on scales of ~100 Mpc/h, and are thus intrinsically different from LCDM model predictions. Correspondingly the probability density function of fluctuations shows a "fat tail" and it is thus different from the Gaussian prediction. Finally we discuss other recent observations which are odds with LCDM predictions and which are, at least theoretically, compatible with the highly inhomogeneous nature of galaxy distribution. We point out that inhomogeneous structures can be fully compatible with statistical isotropy and homogeneity, and thus with a relaxed version of the Cosmological Principle.

Statistical analyses of finite sample distributions usually assume that fluctuations are self-averaging, i.e. that they are statistically similar in different regions of the given sample volume. By using the scale-length method, we test whether this assumption is satisfied in several samples of the Sloan Digital Sky Survey Data Release Six. We find that the probability density function (PDF) of conditional fluctuations, filtered on large enough spatial scales (i.e., r>30 Mpc/h), shows relevant systematic variations in different sub-volumes of the survey. Instead for scales r<30 Mpc/h the PDF is statistically stable, and its first moment presents scaling behavior with a negative exponent around one. Thus while up to 30 Mpc/h galaxy structures have well-defined power-law correlations, on larger scales it is not possible to consider whole sample average quantities as meaningful and useful statistical descriptors. This situation is due to the fact that galaxy structures correspond to density fluctuations which are too large in amplitude and too extended in space to be self-averaging on such large scales inside the sample volumes: galaxy distribution is inhomogeneous up to the largest scales, i.e. r ~ 100 Mpc/h, probed by the SDSS samples. We show that cosmological corrections, as K-corrections and standard evolutionary corrections, do not qualitatively change the relevant behaviors. Finally we show that the large amplitude galaxy fluctuations observed in the SDSS samples are at odds with the predictions of the standard LCDM model of structure formation.(Abridged version).

We consider the conditional galaxy density around each galaxy, and study its fluctuations in the newest samples of the Sloan Digital Sky Survey Data Release 7. Over a large range of scales, both the average conditional density and its variance show a nontrivial scaling behavior, which resembles to criticality. The density depends, for 10 < r < 80 Mpc/h, only weakly (logarithmically) on the system size. Correspondingly, we find that the density fluctuations follow the Gumbel distribution of extreme value statistics. This distribution is clearly distinguishable from a Gaussian distribution, which would arise for a homogeneous spatial galaxy configuration. We also point out similarities between the galaxy distribution and critical systems of statistical physics.

One of the most striking features predicted by standard models of galaxy formation is the presence of anti-correlations in the matter distribution at large enough scales (r>r_c). Simple arguments show that the location of the length-scale r_c, marking the transition from positive to negative correlations, is the same for any class of objects as for the full matter distribution, i.e. it is invariant under biasing. This scale is predicted by models to be at about the same distance of the scale signaling the baryonic acoustic oscillation scale r_{bao}. We test these predictions in the newest SDSS galaxy samples.We find that, in several MG samples, the correlation function remains positive at scales >250 Mpc/h, while in the concordance LCDM it should be negative beyond r_c\approx 120 Mpc/h. In other samples the correlation function becomes negative at scales <50 Mpc/h. To investigate the origin of these differences we consider in detail the propagation of errors on the sample density into the estimation of the correlation function. We conclude that these are important at large enough separations, and that they are responsible for the observed differences between different estimators and for the measured sample to sample variations of the correlation function. We conclude that, in the newest SDSS samples, the large scale behavior of the galaxy correlation function is affected by intrinsic errors andv olume-dependent systematic effects which make the detection of correlations to be only an estimate of a lower limit of their amplitude, spatial extension and statistical errors. We point out that these results represent an important challenge to LCDM models as they largely differ from its predictions.(Abridged version).

When an open system of classical point particles interacting by Newtonian gravity collapses and relaxes violently, an arbitrary amount of energy may in principle be carried away by particles which escape to infinity. We investigate here, using numerical simulations, how this released energy and other related quantities (notably the binding energy and size of the virialized structure) depends on the initial conditions, for the one parameter family of starting configurations given by randomly distributing N cold particles in a spherical volume. Previous studies have established that the minimal size reached by the system scales approximately as N^{-1/3}, a behaviour which follows trivially when the growth of perturbations (which regularize the singularity of the cold collapse in the infinite N limit) are assumed to be unaffected by the boundaries. Our study shows that the energy ejected grows approximately in proportion to N^{1/3}, while the fraction of the initial mass ejected grows only very slowly with N, approximately logarithmically, in the range of N simulated. We examine in detail the mechanism of this mass and energy ejection, showing explicitly that it arises from the interplay of the growth of perturbations with the finite size of the system. A net lag of particles compared to their uniform spherical collapse trajectories develops first at the boundaries and then propagates into the volume during the collapse. Particles in the outer shells are then ejected as they scatter through the time dependent potential of an already re-expanding central core. Using modified initial configurations we explore the importance of fluctuations at different scales, and discreteness (i.e. non-Vlasov) effects in the dynamics.

The properties of the galaxy distribution at large scales are usually studied using statistics which are assumed to be self-averaging inside a given sample. We present a new analysis able to quantitatively map galaxy large scale structures while testing for the stability of average statistical quantities in different sample regions. We find that the newest samples of the Sloan Digital Sky Survey provide unambiguous evidence that galaxy structures correspond to large amplitude density fluctuations at all scales limited only by sample sizes. The two-point correlations properties are self-averaging up to approximately 30 Mpc/h and are characterized by a fractal dimension D=2.1 +- 0.1. Then at all larger scales probed density fluctuations are too large in amplitude and too extended in space to be self-averaging inside the considered volumes. These inhomogeneities are compatible with a continuation of fractal correlations but incompatible with: (i) a homogeneity scale smaller than 100 Mpc/h, (ii) predictions of standard theoretical models, (iii) mock galaxy catalogs generated from cosmological Nbody simuations.

We study statistical properties of galaxy structures in several samples extracted from the 2dF Galaxy Redshift Survey. In particular, we measured conditional fluctuations by means of the scale-length method and determined their probability distribution. In this way we find that galaxy distribution in these samples is characterized by large amplitude fluctuations with a large spatial extension, whose size is only limited by the sample's boundaries. These fluctuations are quite typical and persistent in the sample's volumes, and they are detected in two independent regions in the northern and southern galactic caps. We discuss the relation of the scale-length method to several statistical quantities, such as counts of galaxies as a function of redshift and apparent magnitude. We confirm previous results, which have determined by magnitude and redshift counts that there are fluctuations of about 30% between the southern and the northern galactic caps and we relate explicitly these counts to structures in redshift space. We show that the estimation of fluctuation amplitude normalized to the sample density is biased by systematic effects, which we discuss in detail. We consider the type of fluctuations predicted by standard cosmological models of structure formation in the linear regime and, to study nonlinear clustering, we analyze several samples of mock-galaxy catalogs generated from the distribution of dark matter in cosmological N-body simulations. In this way we conclude that the galaxy fluctuations present in these samples are too large in amplitude and too extended in space to be compatible with the predictions of the standard models of structure formation.

We apply the scale-length method to several three dimensional samples of the Two degree Field Galaxy Redshift Survey. This method allows us to map in a quantitative and powerful way large scale structures in the distribution of galaxies controlling systematic effects. By determining the probability density function of conditional fluctuations we show that large scale structures are quite typical and correspond to large fluctuations in the galaxy density field. We do not find a convergence to homogeneity up to the samples sizes, i.e. ~ 75 Mpc/h. We then measure, at scales r <~ 40 Mpc/h, a well defined and statistically stable power-law behavior of the average number of galaxies in spheres, with fractal dimension D=2.2 +- 0.2. We point out that standard models of structure formation are unable to explain the existence of the large fluctuations in the galaxy density field detected in these samples. This conclusion is reached in two ways: by considering the scale, determined by the linear perturbation analysis of a self-gravitating fluid, below which large fluctuations are expected in standard models and through the determination of statistical properties of mock galaxy catalogs generated from cosmological N-body simulations of the Millenium consortitum.

We discuss the differences and analogies of gravitational clustering in finite and infinite systems. The process of collective, or violent, relaxation leading to the formation of quasi-stationary states is one of the distinguished features in the dynamics of self-gravitating systems. This occurs, in different conditions, both in a finite than in an infinite system, the latter embedded in a static or in an expanding background. We then discuss, by considering some simple and paradigmatic examples, the problems related to the definition of a mean-field approach to gravitational clustering, focusing on role of discrete fluctuations. The effect of these fluctuations is a basic issue to be clarified to establish the range of scales and times in which a collision-less approximation may describe the evolution of a self-gravitating system and for the theoretical modeling of the non-linear phase.

In two recent articles a detailed study has been presented of the out of equilibrium dynamics of an infinite system of self-gravitating points initially located on a randomly perturbed lattice. In this article we extend the treatment of the early time phase during which strong non-linear correlations first develop, prior to the onset of ``self-similar'' scaling in the two point correlation function. We establish more directly, using appropriate modifications of the numerical integration, that the development of these correlations can be well described by an approximation of the evolution in two phases: a first perturbative phase in which particles' displacements are small compared to the lattice spacing, and a subsequent phase in which particles interact only with their nearest neighbor. For the range of initial amplitudes considered we show that the first phase can be well approximated as a transformation of the perturbed lattice configuration into a Poisson distribution at the relevant scales. This appears to explain the ``universality'' of the spatial dependence of the asymptotic non-linear clustering observed from both shuffled lattice and Poisson initial conditions.

We discuss the large scale properties of standard cold dark matter cosmological models characterizing the main features of the power-spectrum, of the two-point correlation function and of the mass variance. Both the real-space statistics have a very well defined behavior on large enough scales, where their amplitudes become smaller than unity. The correlation function, in the range 0<\xi(r)<1, is characterized by a typical length-scale r_c, at which \xi(r_c)=0, which is fixed by the physics of the early universe: beyond this scale it becomes negative, going to zero with a tail proportional to -(r^{-4}). These anti-correlations represent thus an important observational challenge to verify models in real space. The same length scale r_c characterizes the behavior of the mass variance which decays, for r>r_c, as r^{-4}, the fastest decay for any mass distribution. The length-scale r_c defines the maximum extension of (positively correlated) structures in these models. These are the features expected for the dark matter field: galaxies, which represent a biased field, however may have differences with respect to these behaviors, which we analyze. We then discuss the detectability of these real space features by considering several estimators of the two-point correlation function. By making tests on numerical simulations we emphasize the important role of finite size effects which should always be controlled for careful measurements.

We study, using numerical simulations, the dynamical evolution of self-gravitating point particles in static euclidean space, starting from a simple class of infinite ``shuffled lattice'' initial conditions. These are obtained by applying independently to each particle on an infinite perfect lattice a small random displacement, and are characterized by a power spectrum (structure factor) of density fluctuations which is quadratic in the wave number k, at small k. For a specified form of the probability distribution function of the ``shuffling'' applied to each particle, and zero initial velocities, these initial configurations are characterized by a single relevant parameter: the variance $\delta^2$ of the ``shuffling'' normalized in units of the lattice spacing $\ell$. The clustering, which develops in time starting from scales around $\ell$, is qualitatively very similar to that seen in cosmological simulations, which begin from lattices with applied correlated displacements and incorporate an expanding spatial background. From very soon after the formation of the first non-linear structures, a spatio-temporal scaling relation describes well the evolution of the two-point correlations. At larger times the dynamics of these correlations converges to what is termed ``self-similar'' evolution in cosmology, in which the time dependence in the scaling relation is specified entirely by that of the linearized fluid theory. Comparing simulations with different $\delta$, different resolution, but identical large scale fluctuations, we are able to identify and study features of the dynamics of the system in the transient phase leading to this behavior. In this phase, the discrete nature of the system explicitly plays an essential role.

We study the evolution under their self-gravity of infinite ``shuffled lattice'' particle distributions, focussing specifically on the comparison of this evolution with that of ``daughter'' particle distributions, defined by a simple coarse-graining procedure. We consider both the case that such coarse-grainings are performed (i) on the initial conditions, and (ii) at a finite time with a specific additional prescription. In numerical simulations we observe that, to a first approximation, these coarse-grainings represent well the evolution of the two-point correlation properties over a significant range of scales. We note, in particular, that the form of the two-point correlation function in the original system, when it is evolving in the asymptotic ``self-similar'' regime, may be reproduced well in a daughter coarse-grained system in which the dynamics are still dominated by two-body (nearest neighbor) interactions. Using analytical results on the early time evolution of these systems, however, we show that small observed differences between the evolved system and its coarse-grainings at the initial time will in fact diverge as the ratio of the coarse-graining scale to the original inter-particle distance increases. The second coarse-graining studied, performed at a finite time in a specified manner, circumvents this problem. It also makes more physically transparent why gravitational dynamics from these initial conditions tends toward a ``self-similar'' evolution. We finally discuss the precise definition of a limit in which a continuum (specifically Vlasov-like) description of the observed linear and non-linear evolution should be applicable.

We discuss the estimation of galaxy correlation properties in several volume limited samples, in different sky regions, obtained from the Fourth Data Release of the Sloan Digital Sky Survey. The small scale properties are characterized through the determination of the nearest neighbor probability distribution. By using a very conservative statistical analysis, in the range of scales [0.5,~30] Mpc/h we detect power-law correlations in the conditional density in redshift space, with an exponent \gamma=1.0 \pm 0.1. This behavior is stable in all different samples we considered thus it does not depend on galaxy luminosity. In the range of scales [~30,~100] Mpc/h we find evidences for systematic unaveraged fluctuations and we discuss in detail the problems induced by finite volume effects on the determination of the conditional density. We conclude that in such range of scales there is an evidence for a smaller power-law index of the conditional density. However we cannot distinguish between two possibilities: (i) that a crossover to homogeneity (corresponding to \gamma=0 in the conditional density) occurs before 100 Mpc/h, (ii) that correlations extend to scales of order 100 Mpc/h (with a smaller exponent 0 < \gamma <1). We emphasize that galaxy distributions in these samples present large fluctuations at the largest scales probed, corresponding to the presence of large scale structures extending up to the boundaries of the present survey. Finally we discuss several differences between the behavior of the conditional density in mock galaxy catalogs built from cosmological N-body simulations and real data. We discuss some theoretical implications of such a fact considering also the super-homogeneous features of primordial density fields.

We study the statistics of the force felt by a particle in the class of spatially correlated distribution of identical point-like particles, interacting via a $1/r^2$ pair force (i.e. gravitational or Coulomb), and obtained by randomly perturbing an infinite perfect lattice. In the first part we specify the conditions under which the force on a particle is a well defined stochastic quantity. We then study the small displacements approximation, giving both the limitations of its validity, and, when it is valid, an expression for the force variance. In the second part of the paper we extend to this class of particle distributions the method introduced by Chandrasekhar to study the force probability density function in the homogeneous Poisson particle distribution. In this way we can derive an approximate expression for the probability distribution of the force over the full range of perturbations of the lattice, i.e., from very small (compared to the lattice spacing) to very large where the Poisson limit is recovered. We show in particular the qualitative change in the large-force tail of the force distribution between these two limits. Excellent accuracy of our analytic results is found on detailed comparison with results from numerical simulations. These results provide basic statistical information about the fluctuations of the interactions (i) of the masses in self-gravitating systems like those encountered in the context of cosmological N-body simulations, and (ii) of the charges in the ordered phase of the One Component Plasma.

We present a perturbative treatment of the evolution under their mutual self-gravity of particles displaced off an infinite perfect lattice, both for a static space and for a homogeneously expanding space as in cosmological N-body simulations. The treatment, analogous to that of perturbations to a crystal in solid state physics, can be seen as a discrete (i.e. particle) generalization of the perturbative solution in the Lagrangian formalism of a self-gravitating fluid. Working to linear order, we show explicitly that this fluid evolution is recovered in the limit that the initial perturbations are restricted to modes of wavelength much larger than the lattice spacing. The full spectrum of eigenvalues of the simple cubic lattice contains both oscillatory modes and unstable modes which grow slightly faster than in the fluid limit. A detailed comparison of our perturbative treatment, at linear order, with full numerical simulations is presented, for two very different classes of initial perturbation spectra. We find that the range of validity is similar to that of the perturbative fluid approximation (i.e. up to close to ``shell-crossing''), but that the accuracy in tracing the evolution is superior. The formalism provides a powerful tool to systematically calculate discreteness effects at early times in cosmological N-body simulations.

We discuss the formation of the first structures in gravitational N-body simulations. The role of two-body interaction is found to be a crucial element and an analogy with the dynamics of the Coulomb lattice, well-studied in solid state physics, is discussed.

We study galaxy correlations from samples extracted from the 2dFGRS final release. Statistical properties are characterized by studying the nearest neighbor probability density, the conditional density and the reduced two-point correlation function. The result is that the conditional density has a power-law behavior in redshift space described by an exponent \gamma=0.8 \pm 0.2 in the interval from about 1 Mpc/h, the average distance between nearest galaxies, up to about 40 Mpc/h, corresponding to radius of the largest sphere contained in the samples. These results are consistent with other studies of the conditional density and are useful to clarify the subtle role of finite-size effects on the determination of the two-point correlation function in redshift and real space

We discuss some of the basic implications of recent results on galaxy correlations published by the SDSS collaboration. In particular we focus on the evidence which has been recently presented for the scale and nature of the transition to homogeneity in the galaxy distribution, and results which describe the dependence of clustering on luminosity. The two questions are in fact strictly entangled, as the stability of the measure of the amplitude of the correlation function depends on the scale at which the mean density becomes well defined. We note that the recent results which indicate the convergence to well defined homogeneity in a volume equivalent to that of a sphere of radius 70 Mpc/h, place in doubt previous detections of ``luminosity bias'' from measures of the amplitude of the correlation function. We emphasize that the way to resolve these issues is to first use, in volume limited samples corresponding to different ranges of luminosity, the unnormalized two point statistics to establish the scale (and value) at which the mean density becomes well defined. We note also that the recent SDSS results for these statistics are in good agreement with those obtained by us through analyses of many previous samples, confirming in particular that the galaxy distribution is well described by a fractal dimension D ~ 2 up to a scale of at least 20 Mpc/h. We discuss critically the agreement of this new data with current theoretical models.

One of the main problems in the studies of large scale galaxy structures concerns the relation of the correlation properties of a certain population of objects with those of a selected subsample of it, when the selection is performed by considering physical quantities like luminosity or mass. I consider the case where the sampling is defined as in the simplest thresholding selection scheme of the peaks of a Gaussian random field as well as the case of the extraction of point distributions in high density regions from gravitational N-body simulations. I show that an invariant scale under sampling is represented by the zero-crossing scale of \xi(r). By considering recent measurements in the 2dF and SDSS galaxy surveys I note that the zero-point crossing length has not yet been clearly identified, while a dependence on the finite sample size related to the integral constraint is manifest. I show that this implies that other length scales derived from \xi(r) are also affected by finite size effects. I discuss the theoretical implications of these results, when considering the comparison of structures formed in N-body simulations and observed in galaxy samples, and different tests to study this problem.

We apply a simple linearization, well known in solid state physics, to approximate the evolution at early times of cosmological N-body simulations of gravity. In the limit that the initial perturbations, applied to an infinite perfect lattice, are at wavelengths much greater than the lattice spacing $l$ the evolution is exactly that of a pressureless self-gravitating fluid treated in the analagous (Lagrangian) linearization, with the Zeldovich approximation as a sub-class of asymptotic solutions. Our less restricted approximation allows one to trace the evolution of the discrete distribution until the time when particles approach one another (i.e. ``shell crossing''). We calculate modifications of the fluid evolution, explicitly dependent on $l$ i.e. discreteness effects in the N body simulations. We note that these effects become increasingly important as the initial red-shift is increased at fixed $l$. The possible advantages of using a body centred cubic, rather than simple cubic, lattice are pointed out.

Cosmic structures at large scales represent the earliest and most extended form of matter condensation. In this lecture we review the application of the methods and concepts of modern statistical physics to these structures. This leads to a new perspective in the field which can be tested by the many new data which are appearing in the near future. In particular, galaxy structures show fractal correlation up to the present observational limits. The cosmic microwave background radiation, which should trace the initial conditions from which these structures have emerged through gravitational dynamics, is instead extremely smooth. Understanding the relation between the complex galaxy structures and the smooth microwave background radiation represents an extremely challenging problem in the field of structure formation.

In order to investigate the properties of large-scale structures of galaxies in the universe, we present an analysis of their spatial distribution at z$<$0.033. We used the LEDA extragalactic database containing over 1 million of galaxies covering the all-sky and the SDSS data included in the public release DR1, yielding to a sample of around 134,000 galaxies having a measured redshift in two survey areas representing 690 sq. degrees. The results of the study are 2D, 3D maps and magnitude number counts of galaxies, drawn from B-band samples. Movies and high resolution figures are available on http://www-obs.univ-lyon1.fr/cosmologie/helene.html

We present an analysis of different sets of gravitational N-body simulations, all describing the dynamics of discrete particles with a small initial velocity dispersion. They encompass very different initial particle configurations, different numerical algorithms for the computation of the force, with or without the space expansion of cosmological models. Despite these differences we find in all cases that the non-linear clustering which results is essentially the same, with a well-defined simple power-law behaviour in the two-point correlations in the range from a few times the lower cut-off in the gravitational force to the scale at which fluctuations are of order one. We argue, presenting quantitative evidence, that this apparently universal behaviour can be understood by the domination of the small scale contribution to the gravitational force, coming initially from nearest neighbor particles.

In the gravitational evolution of a cold infinite particle distribution, two-body interactions can be predominant at early times: we show that, by treating the simple case of a Poisson particle distribution in a static universe as an ensemble of isolated two-body systems, one may capture the origin of the first non-linear correlated structures. The developed power-law like behavior of the two-point correlation function is then simply related to the functional form of the time evolved nearest-neighbor probability distribution, whose time dependence can be computed by using Liouville theorem for the gravitational two-body problem. We then show that a similar dynamical evolution is also found in a large-scale ordered distribution, which has striking similarities to the case of a cosmological CDM simulation which we also consider.

We discuss correlation properties of a general mass density field introducing a classification of structures based on their complexity. Standard cosmological models for primordial mass fluctuations are characterized by a sort of large-scale stochastic order, that we call super-homogeneity to highlight the fact that mass fluctuations increase as a function of scale in the slowest possible way for any stochastic mass field. On the other hand the galaxy spatial distribution show complex structures with a high degree of inhomogeneity and fractal-like spatial correlations up to some relevant cosmological scale. The theoretical problem of cosmological structure formation should then explain the growth of strongly correlated and non-linear structures from the very uniform field of density fluctuations given as standard initial condition.

We reply to some comments in astro-ph/0309381 concerning the problem of setting-up initial conditions in cosmological N-body simulations

We study the statistics of the gravitational (Newtonian) force in a particular kind of weakly correlated distribution of point-like and unitary mass particles generated by the so-called Gauss-Poisson point process. In particular we extend to these distributions the analysis a' la Chandrasekhar introduced for purely Poisson processes. In this way we can find the asymptotic behavior of the probability density function of the force for large values of the field as a generalization of the Holtzmark statistics. The validity of the introduced approximations is positively tested through a direct comparison with the analysis of the statistics of the gravitational force in numerical simulations of Gauss-Poisson processes. Moreover the statistics of the force felt by a particle due to only its first nearest neighbor is analytically a numerical studied, resulting to be the dominant contribution to the total force.

Threshold biasing of a Gaussian random field gives a linear amplification of the reduced two point correlation function at large distances. We show that for standard cosmological models this does not translate into a linear amplification of the power spectrum (PS) at small k. For standard CDM type models this means that the``turn-over'' at small k of the original PS disappears in the PS of the biased field for the physically relevant range of the threshold parameter. In real space this difference is manifest in the asymptotic behaviour of the normalised mass variance in spheres of radius R, which changes from the ``super-homogeneous'' behaviour to a Poisson-like behaviour.This qualitative change results from the intrinsic stochasticity of the threshold sampling. While our quantitative results are specific to the simplest threshold biasing model, we argue that our qualitative conclusions should be valid generically for any biasing mechanism involving a scale-dependent amplification of the correlation function. One implication is that the real-space correlation function will be a better instrument to probe for the underlying Harrison Zeldovich spectrum in the distribution of visible matter, as the characteristic asymptotic negative power-law \xi (r) \sim -r^{-4} tail is undistorted by biasing.

The initial conditions describing seed fluctuations for the formation of structure in standard cosmological models, i.e.the Harrison-Zeldovich distribution, have very characteristic ``super-homogeneous'' properties: they are statistically translation invariant, isotropic, and the variance of the mass fluctuations in a region of volume V grows slower than V. We discuss the geometrical construction of distributions of points in ${\bf R}^3$ with similar properties encountered in tiling and in statistical physics, e.g. the Gibbs distribution of a one-component system of charged particles in a uniform background (OCP). Modifications of the OCP can produce equilibrium correlations of the kind assumed in the cosmological context. We then describe how such systems can be used for the generation of initial conditions in gravitational $N$-body simulations.

We analyse with simple real-space statistics the Virgo consortium's cosmological N-body simulations. Significant clustering rapidly develops well below the initial mean interparticle separation \Lambda_i, where the gravitational force on a particle is dominated by that with its nearest neighbours. A power-law behaviour in the two point correlation function emerges, which in the subsequent evolution is continuously amplified and shifted to larger scales, in a roughly self-similar manner. We conclude that the density fluctuations at the smallest scales due to the particle-like nature of the distribution being evolved are thus essential in the development of these correlations, and not solely, as usually supposed, the very small continuous (fluid-like) fluctuations at scales larger than \Lambda_i >.

In this lecture we address three different but related aspects of the initial continuous fluctuation field in standard cosmological models. Firstly we discuss the properties of the so-called Harrison-Zeldovich like spectra. This power spectrum is a fundamental feature of all current standard cosmological models. In a simple classification of all stationary stochastic processes into three categories, we highlight with the name ``super-homogeneous'' the properties of the class to which models like this, with $P(0)=0$, belong. In statistical physics language they are well described as glass-like. Secondly, the initial continuous density field with such small amplitude correlated Gaussian fluctuations must be discretised in order to set up the initial particle distribution used in gravitational N-body simulations. We discuss the main issues related to the effects of discretisation, particularly concerning the effect of particle induced fluctuations on the statistical properties of the initial conditions and on the dynamical evolution of gravitational clustering.

After reviewing some basic relevant properties of stationary stochastic processes (SSP), we discuss the properties of the so-called Harrison-Zeldovich like spectra of mass density perturbations. These correlations are a fundamental feature of all current standard cosmological models. Examining them in real space we note they imply a "sub-poissonian" normalised variance in spheres $\sigma_M^2(R) \sim R^{-4} \ln R$. In particular this latter behaviour is at the limit of the most rapid decay ($\sim R^{-4}$) of this quantity possible for any stochastic distribution (continuous or discrete). In a simple classification of all SSP into three categories, we highlight with the name ``super-homogeneous'' the properties of the class to which models like this, with $P(0)=0$, belong. In statistical physics language they are well described as lattice or glass-like. We illustrate their properties through two simple examples: (i) the ``shuffled'' lattice and the One Component Plasma at thermal equilibrium.

After reviewing the basic relevant properties of stationary stochastic processes (SSP), defining basic terms and quantities, we discuss the properties of the so-called Harrison-Zeldovich like spectra. These correlations, usually characterized exclusively in k-space (i.e. in terms of power spectra P(k)), are a fundamental feature of all current standard cosmological models. Examining them in real space we note their characteristics to be a {\it negative} power law tail \xi(r) \sim - r^{-4} and a {\it sub-poissonian} normalised variance in spheres \sigma^2(R) \sim R^{-4} \ln R. We note in particular that this latter behaviour is at the limit of the most rapid decay (\sim R^{-4}) of this quantity possible for any stochastic distribution (continuous or discrete). This very particular characteristic is usually obscured in cosmology by the use of Gaussian spheres. In a simple classification of all SSP into three categories, we highlight with the name ``super-homogeneous'' the properties of the class to which models like this, with P(0)=0, belong. In statistical physics language they are well described as glass-like. They do not have either ``scale-invariant'' features, in the sense of critical phenomena, nor fractal properties. We illustrate their properties with some simple examples, in particular that of a ``shuffled'' lattice.

The study of the properties of cosmic structures in the universe is one of the most fascinating subject of the modern cosmology research. Far from being predicted, the large scale structure of the matter distribution is a very recent discovery, which continuosly exhibits new features and issues. We have faced such topic along two directions; from one side we have studied the correlation properties of the cosmic structures, that we have found substantially different from the commonly accepted ones. From the other side, we have studied the statistical properties of the very simplified system, in the attempt to capture the essential ingredients of the formation of the observed structures.

Cosmological N-body simulations aim to calculate the non-linear gravitational growth of structures via particle dynamics. A crucial problem concerns the setting-up of the initial particle distribution, as standard theories of galaxy formation predict the properties of the initial continuous density field with small amplitude correlated Gaussian fluctuations. The discretisation of such a field is a complex issue and particle fluctuations are especially relevant at small scales where non-linear dynamics firstly takes place. In general, most of the procedures which may discretise a continuous field, gives rise to Poisson noise, which would then dominate the non-linear small-scales dynamics due to nearest-neighbours interactions. In order to avoid such a noise, and to consider the dynamics as due only to large scale (smooth) fluctuations, an ad-hoc method (lattice or glassy system plus correlated displacements) has been introduced and used in cosmological simulations. We show that such a method gives rise to a particle distribution which does not have any of the correlation properties of the theoretical continuous density field. This is because discreteness effects, different from Poisson noise but nevertheless very important, determine particle fluctuations at any scale, making it completely different from the original continuous field. We conclude that discreteness effects play a central role in the non-linear evolution of N-body simulations.

Ideas of Statistical Physics are very relevant for cosmic structures especially considering that the field is undergoing a period of exceptional development with many new data appearing on a monthly basis. In the past years we have focused mostly on galaxy distributions and their statistical properties. This led to an interesting debate which will be resolved by the next generation of data in a couple of years. In addition here we discuss the statistical properties of the fluctuations of the cosmic microwave background which are small in amplitude but complex in structure. We finally discuss the connection between these observations and the Harrison-Zeldovich spectrum and its further implications on the theories of structure formation and the and the cosmological $N$-body simulations.

In most direct estimates of the mass density (visible or dark) of the Universe, a central input parameter is the luminosity density of the Universe. Here we consider the measurement of this luminosity density from red-shift surveys, as a function of the yet undetermined characteristic scale R_H at which the spatial distribution of visible matter tends to a well defined homogeneity. Making the canonical assumption that the cluster mass to luminosity ratio M/L is the universal one, we can estimate the total mass density as a function \Omega_m(R_H,M/L). Taking the highest estimated cluster value M/L ~300h and a conservative lower limit R_H > 20 Mpc/h, we obtain the upper bound \Omega_m < 0.1 . We note that for values of the homogeneity scale R_H in the range R_H ~ (90 +/- 45) hMpc, the value of \Omega_m may be compatible with the nucleosynthesis inferred density in baryons.

Fractal properties are usually characterized by means of various statistical tools which deal with spatial average quantities. Here we focus on the determination of fluctuations around the average counts and we develop a test for the study of galaxy distribution both in redshift and magnitude space. Fluctuations in the counts of galaxies, in a fractal distribution, are of the same order of the average number at all scales as a function of redshift and magnitude. We point out that the study of these kind of fluctuations can be a powerful test to understand the nature of galaxy clustering at very large scales.

The question of the nature of galaxy clustering and the possible homogeneity of galaxy distribution is one of the fundamental problem of cosmology. It is well established that galaxy structures are characterized, up to a certain scale, by fractal properties. The possible crossover to homogeneity is instead still matter of debate. However, independently on the specific value of the homogeneity scale, the fractal nature of galaxy clustering requires new methods and theoretical concepts developed in the area of statistical physics and complexity. In this lecture we discuss a new perspective on the problem of cosmological structures formation, both from the experimental and the theoretical points of view. This new perspective leads to a very interesting and constructive interaction between the fields of cosmic structures, statistical physics and complexity with very challenging open problems which we also discuss.

In this lecture we clarify the basic difference between the correlation properties for systems characterized by small or large fluctuations. The concepts of correlation length, homogeneity scale, scale invariance and criticality are discussed as well. We relate these concepts to the interpretation of galaxy clsutering.

We present a brief introduction to the statistical properties of systems with large fluctuations. We point out that for such systems the relevant statistical quantities are scaling exponents and the nature of fluctuations is completely different from the one belonging to system with small fluctuations. We then propose two test to be performed on galaxy counts data as a function of apparent magnitude, which may clarify in a direct way the nature of the large scale galaxy clustering

Counts of galaxies as a function of apparent magnitude are among the most time-honored observations in cosmology. In this Letter, we focus on some statistical properties of these counts which are fundamental in order to characterize the large scale correlations in the galaxy spatial distribution. There are, in fact, no longer doubts since two decades of the existence of very large scale structures. The still remaining problem concerns the correct characterization of their statistical properties. We propose to study two properties of galaxy counts data, in order to discriminate between a small scale (~ 5 - 20 Mpc/h) homogeneous distribution, and a fractal structure on large scales (~ 20 - 300 Mpc/h). Firstly, the average slope of the counts which can be associated to an eventual fractal dimension in real space by simple arguments. Secondly, we propose to study fluctuations of counts around the average behavior as a function of apparent magnitude in the whole magnitude range B_j > 11. These fluctuations can discriminate between a genuine fractal distribution and a homogeneous one. In fact, they are related to the very statistical properties of the spatial distribution, independently on cosmological corrections. More specifically, in a fractal distribution one expects to find persistent scale-invariant fluctuations around the average behavior, which do not decay with apparent magnitude. On the other hand, in an homogeneous distribution, on large enough scales, the relative variance of the counts should decrease exponentially with apparent magnitude. Smooth cosmological corrections cannot change such a behavior.

The clustering of galaxies is well characterized by fractal properties, with the presence of an eventual cross-over to homogeneity still a matter of considerable debate. In this letter we discuss the cosmological implications of a fractal distribution of matter, with a possible cross-over to homogeneity at an undetermined scale R_{homo}. Contrary to what is generally assumed, we show that, even when R_{homo} -> \infty, this possibility can be treated consistently within the framework of the expanding universe solutions of Friedmann. The fractal is a perturbation to an open cosmology in which the leading homogeneous component is the cosmic background radiation (CBR). This cosmology, inspired by the observed galaxy distributions, provides a simple explanation for the recent data which indicate the absence of deceleration in the expansion (q_o \approx 0). Correspondingly the `age problem' is also resolved. Further we show that the model can be extended back from the curvature dominated arbitrarily deep into the radiation dominated era, and we discuss qualitatively the modifications to the physics of the anisotropy of the CBR, nucleosynthesis and structure formation.

In this letter we show that in a Gaussian random field the correlation length, the typical size of correlated structures, does not change with biasing. We interpret the amplification of the correlation functions of subsets identified by different thresholds being due to the increasing sparseness of peaks over threshold. This clarifies an long-standing misconception in the literature. We also argue that this effect does not explain the observed increase of the amplitude of the correlation function xi(r) when galaxies of brighter luminosity or galaxy clusters of increasing richness are considered.

Galaxy structures are certainly fractal up to a certain crossover scale \lambda_0. A clear determination of such a scale is still missing. Usually, the conceptual and practical implications of this property are neglected and the structures are only discussed in terms of their global amplitude. Here we present a compact summary of these implications. First, we discuss the problem of the identification of the crossover scale \lambda_0 and the proper characterization of the scaling. We then consider the implications of these properties with respect to various physical phenomena and to the corresponding characteristic values, i.e. r_0, \sigma_8, \Omega, etc. These implications crucially depend on the value of \lambda_0, but they are still important for a relatively small value, say \lambda_0 \approx 50 \hmp. Finally we consider the main theoretical consequences of these results.

In the debate about galaxy correlation there are different questions which can be addressed separately: Which are the statistical methods able to properly detect scale invariance and describe, in general, the properties of irregular and regular distributions? Which are the implications for cosmology of the fractal behavior of galactic structures, up to a certain scale lambda_0? Which is the homogeneity scale lambda_0, i.e. the scale beyond which galaxy distribution has an eventual crossover to homogeneity? These are three different, but related, problems, which must be considered in different steps, from the point of view of data analysis as well as from the theoretical perspective.

In a recent analysis of number counts in the ESP survey Scaramella et al. (1998) claim to find evidence for a cross-over to homogeneity at large scales, and against a fractal behaviour with dimension $D \approx 2$. In this comment we note firstly that, if such a cross-over exists as described by the authors, the scale characterizing it is ~ 100 - 300 Mpc/h. This invalidates the ``standard'' analysis of the same catalogue given elsewhere by the authors which results in a ``correlation length'' of only r_0 = 4 Mpc/h. Furthermore we show that the evidences for a cross-over to homogeneity rely on the choice of cosmological model, and most crucially on the so called K corrections. We show that the D ~ 3 behaviour seen in the K-corrected data of Scaramella et al. is in fact unstable, increasing systematically towards D=4 as a function of the absolute magnitude limit. This behaviour can be quantitatively explained as the effect of an unphysical K-correction in the relevant range of red-shift (z ~ 0.1- 0.3). A more consistent interpretation of the number counts is that D is in the range 2 - 2.5, depending on the cosmological model, consistent with the continuation of the fractal D ~ 2 behaviour observed at scales up to ~100 Mpc/h. This implies a smaller K-correction. Given, however, the uncertainty in the effect of intrinsic fluctuations on the number counts statistic, and its sensitivity on these large scales to the uncertain K corrections, we conclude that it is premature to put a definitive constraint on the galaxy distribution using the ESP data alone.

We report our analysis of the properties of galaxy clustering for a new redshift sample of galaxies, the CfA2-South catalog, using statistical methods which do not rely on the assumption of homogeneity. We find that, up to ~ 20 Mpc/h, which is the largest scale to which correlation properties can be reliably inferred, the galaxy clustering is scale-invariant and characterized by a fractal dimension D=1.9 \pm 0.1. Further there is no statistical evidence for homogeneity at any of the larger scales (up to ~150 Mpc/h) probed more weakly by the catalog. These results means that characteristic ``correlation lengths'' for the clustering of galaxies derived using standards methods of analysis are not meaningful. Further the results are very consistent with those obtained from many other catalogs using the methods adopted here, which show the D =2 fractal continuing to beyond 100 Mpc/h. The incompleteness of the relevant data conjectured by various authors to give rise to such behaviour is therefore proved to have no significant effect (up to 20 Mpc/h) on the measured correlations.

We study the (newtonian) gravitational force distribution arising from a fractal set of sources. We show that, in the case of real structures in finite samples, an important role is played by morphological properties and finite size effects. For dimensions smaller than d-1 (being $d$ the space dimension) the convergence of the net gravitational force is assured by the fast decaying of the density, while for fractal dimension D>d-1 the morphological properties of the structure determine the eventual convergence of the force as a function of distance. We clarify the role played by the cut-offs of the distribution. Some cosmological implications are discussed.

In this letter we present an idea which reconciles a homogeneous and isotropic Friedmann universe with a fractal distribution of galaxies. We use two observational facts: The flat rotation curves of galaxies and the (still debated) fractal distribution of galaxies with fractal dimension D=2. Our idea can also be interpreted as a redefinition of the notion of bias.

The fact that galaxy distribution exhibits fractal properties is well established since twenty years. Nowadays, the controversy concerns the range of the fractal regime, the value of the fractal dimension and the eventual presence of a cross-over to homogeneity. Fractal properties maybe studied with methods which do not assume homogeneity a priori as the standard statistical methods do. We show that complementary to the adoption of new methods of analysis there are important theoretical implications for the usual scenario of galaxy formation. For example, we focus on the concept of bias and we show that it needs a basic revision even if future redshift surveys will be able to identify an eventual tendency to homogenization.

The debate on the correlation properties of galaxy structures has having an increasing interest during the last year. In this lecture we discuss the claims of different authors who have criticized our approach and results. In order to have a clear cut of the situation, we focus mainly on galaxy distribution in the intermediate range of distances ~ 100 - 200 Mpc/h. In particular we discuss: (i) the validity of the actual data and the concept of "fair sample", (ii) the shift of r_0 with sample depth and luminosity bias, (iii) the value of the fractal dimension, (iv) the problem of the counts from a single point and the case of ESP, (iv) uniformity of angular catalogs. The detection of fractal behavior up to ~ 100 - 200 Mpc/h is enough to rise serious problems to the usual statistical methods used for the characterization of galaxy correlations, standard interpretation of galaxy distribution and theoretical models develoed. The clarification of the intermediate scale behavior is very instructive for the subsequent interpretation of the very large scale galaxy distribution.

There is a general agreement that galaxy structures exhibit fractal properties, at least up to some small scale. However the presence of an eventual crossover towards homogenization, as well as the exact value of the fractal dimension, are still a matter of debate. I summarize the main points of the this discussion, considering also some galaxy surveys which have recently appeared. Further I discuss the implications for the standard picture of the observed fractal behaviour in galaxy distribution. In particular I consider the co-existence of fractal structures and the linear Hubble-law within the same scales. This fact represents a challenge for the standard cosmology where the linear Hubble law is a strict consequence of homogeneity of the expanding universe. Finally I consider the comparison of CDM-like models with the data noting that the simulations are not able to reproduce the observed properties of galaxy correlations.

Two fundamental empirical laws have been established in the analysis of galaxy space distribution. First, recent analyses have revealed that the three dimensional distribution of galaxies and clusters is characterized by large scale structures and huge voids: such a distribution shows fractal correlations up to the limits of the available samples. This has confirmed the earlier de Vaucouleurs power-law density - distance relation, now corresponding to a fractal structure with dimension $D \approx 2$, at least, in the range of scales $ \sim 1 \div 200 Mpc$ ($H_0 = 55 km/sec/Mpc$). An eventual cut-off towards homogenization has not been yet identified. Second, since Hubble's discovery, the linear redshift-distance law has been well established within $200 Mpc$ and also much deeper. The co-existence of these laws within the same scales is a challenge for the standard cosmology where the linear Hubble law is a strict consequence of homogeneity of the expanding universe. This puzzle is now sufficiently strong to raise doubts for the standard cosmology.

This comment is in response to the paper by L. Guzzo recently appeared in ``New Astronomy'' related to our work. The subject of the discussion concerns the correlation properties of galaxy distribution in the available 3-d samples. There is a general agreement that galaxy structures exhibit fractal properties, at least up to some scale. However the presence of an eventual crossover towards homogenization, as well as the exact value of the fractal dimension, are still matter of debate. Here we briefly summarize our point of view by discussing three main topics. The first one is methodological, i.e. we clarify which are the correct methods to detect the real correlations properties of the 3-d galaxy distribution. Then we discuss the results of the analysis of several samples in two ranges of scales. In the first range of scale, below $100 \div 200 \hmp$, the statistical quality of the data is rather good, and we find that galaxy distribution has fractal properties with $D \approx2$. At larger distances the statistical robustness of the present data is weaker, but, we show that there is evidence for a continuation of the fractal behavior, without any tendency towards homogenization.

We present the statistical properties of large scale galaxy distribution in the LEDA redshift database. This catalog contains more than 40,000 redshifts over all the sky. We find that LEDA, although seriously affected by incompleteness, shows quite stable statistical properties. In particular, we have considered the behaviour of the two points correlation function and of the power spectrum of the density fluctuations, and we have done several tests to check whether the incompleteness of the catalog affect these statistical quantities. Our conclusion is that galaxy distribution in this catalog has fractal properties up to the distance 200 Mpc with D = 2.1 \pm 0.2, with no sign towards homogeneity: this result is statistically stable. We analyze also the properties of the angular distributions, finding a complete agreement with the results obtained in redshift space. Finally, we compare these results with those obtained in other redshift surveys, finding that this sample well reproduces the properties of galaxy distribution found in the different catalogs.

We investigate the statistical properties of the APM-Stromlo redshift survey by using the concepts and methods of modern Statistical Physics. We find that galaxy distribution in this survey exhibits scale invariant properties with fractal dimension $D = 2.1 \pm 0.1$, up to $\sim 40 \hmp$, i.e. the limit of its statistical validity, without any tendency towards homogenization. No intrinsic characteristic scales are definitely found in this galaxy sample. We present several tests to study the statistical reliability of the results. \

Some years ago we proposed a new approach to the analysis of galaxy and cluster correlations based on the concepts and methods of modern statistical Physics. This led to the surprising result that galaxy correlations are fractal and not homogeneous up to the limits of the available catalogs. The usual statistical methods, which are based on the assumption of homogeneity, are therefore inconsistent for all the length scales probed so far, and a new, more general, conceptual framework is necessary to identifythe real physical properties of these structures. In the last few years the 3-d catalogs have been significatively improved and we have extended our methods to the analysis of number counts and angular catalogs. This has led to a complete analysis of all the available data that we present in this review. The result is that galaxy structures are highly irregular and self-similar: all the available data are consistent with each other and show fractal correlations (with dimension $D \simeq 2$) up to the deepest scales probed so far ($1000 \hmp$) and even more as indicated from the new interpretation of the number counts. The evidence for scale-invariance of galaxy clustering is very strong up to $150 \hmp$ due to the statistical robustness of the data but becomes progressively weaker (statistically) at larger distances due to the limited data. In These facts lead to fascinating conceptual implications about our knowledge of the universe and to a new scenario for the theoretical challenge in this field.

We discuss various questions which arise when one considers the central projection of three dimensional fractal sets (galaxy catalogs) onto the celestial globe. The issues are related to how fractal such projections look. First we show that the lacunarity in the projection can be arbitrarily small. Further characteristics of the projected set---in particular scaling---depend sensitively on how the apparent sizes of galaxies are taken into account. Finally, we discuss the influence of opacity of galaxies. Combining these ideas, seemingly contradictory statements about lacunarity and apparent projections can be reconciled.

We introduce the basic techniques used for the analysis of three dimensional and two dimensional galaxy samples. We report the correlation analysis of various redshift surveys which shows that the available data are consistent with each other and manifest fractal correlations (with dimension $D \simeq 2$) up to the present observational limits without any tendency towards homogenization. This result points to a new interpretation of the galaxy number counts. We show that an analysis of the small scale fluctuations allows us to reconcile the correlation analysis and the number counts in a new perspective which has a number of important implications.

The properties of spatial distribution of luminous matter are investigated analysing all the available three dimensional catalogues of galaxies. In standard view, galaxies are believed to have a fractal distribution at small scale with a crossover to an homogeneous one at large scale. However up to now, the quantitative determination of this presumed homogeneity scale is still lacking. Contrary to such expectations, observational results show, in fact, a very inhomogeneous galaxy distribution. Some years ago we criticise the standard statistical approach and proposed a new one based on the concepts and methods of modern statistical analysis. The main result of new analysis is that, contrary to the conclusion of standard methods, the distribution of galaxies in the available samples, does not show any crossover to homogeneity, but has fractal correlations (with dimension $D \approx 2$) up to the limits of present three dimensional catalogs ($\approx 1000 h^{-1} Mpc$). The very first consequence of this result is that the standard approach is incorrect for all the length scale probed until now; moreover it calls for fascinating conceptual implication for the theoretical challenge in this field.

In the last years there has been a growing interest in the understanding a vast variety of scale invariant and critical phenomena occurring in nature. Experiments and observations indeed suggest that many physical systems develop spontaneously correlations with power law behaviour both in space and time. Pattern formation, aggregation phenomena, biological systems, geological systems, disordered materials, clustering of matter in the universe are just some of the fields in which scale invariance has been observed as a common and basic feature. However, the fact that certain structures exhibit fractal and complex properties does not tell us why this happens. A crucial point to understand is therefore the origin of the general scale-invariance of natural phenomena. This would correspond to the understanding of the origin of fractal structures and of the properties of Self- Organized Criticality (SOC) from the knowledge of the microscopic physical processes at the basis of these phenomena. Fractal geometry can play a crucial role in extracting the correct physical properties from experimental data. In particular, the recent availability of complete three dimensional samples of galaxies and clusters permits a direct study of their spatial properties. We present a brief review of galaxy correlations based on the methods of modern Statistical Physics. These methods are able to identify self-similar and non-analytical properties, and allow us to test the usual homogeneity assumption of luminous matter distribution. The new analysis shows that all the available data are consistent with each other and show fractal correlations (with dimension $D \simeq 2$) up to the deepest scales probed until now ($1000 \hmp$)

We report the correlation analysis of various redshift surveys which shows that the available data are consistent with each other and manifest fractal correlations (with dimension $D \simeq 2$) up to the present observational limits ($\approx 150 \hmp$) without any tendency towards homogenization. This result points to a new interpretation of the number counts that represents the main subject of this letter. We show that an analysis of the small scale fluctuations allows us to reconcile the correlation analysis and the number counts in a new perspective which has a number of important implications.

We analyze subsamples of Abell and ACO cluster catalogs, in order to study the spatial properties of the large scale matter distribution. The subsamples analyzed are estimated to be nearly complete and are the standard ones used in the correlation analysis by various authors. The statistical analysis of cluster correlations is performed without any assumptions on the nature of the distribution itself. The cluster samples show fractal correlations up to sample limits ($\approx 70 h^{-1} Mpc$) with fractal dimension $D \approx 2$, without any tendency towards homogenization. Our analysis shows that the standard correlation methods are incorrect, since they give a finite correlation length for a distribution that does not possess one. In particular, the so called correlation length $r_{0}$ is shown to be simply a fraction of the sample size. Moreover we conclude that galaxies and clusters are two different representations of the same self similar structure and that the correlations of clusters are the continuation of those of galaxies to larger scales.

In this lecture we will try to address the "frequently asked questions about fractals" in the field of large scale galaxy distribution. This paper takes its origin from a very interesting discussion we had at this meeting. A lot of points were raised, and we try to make clear the fundamental ones.

Some years ago we proposed a new approach for the analysis of galaxy and cluster correlations based on the concepts and methods of modern statistical physics. This led to the surprising result that galaxy correlations are fractal and not homogeneous up to the limits of the available catalogues. In the meantime many more redshifts have been measured and we have extended our methods also to the analysis of number counts and angular catalogues. This has led to a complete analysis of all the available data that we are going to present in detail in this lecture. In particular we will discuss the properties of the following catalogues: CfA, Perseus-Pisces, SSRS, IRAS, LEDA, Las Campanas and ESP for galaxies and Abell and ACO for clusters. The usual statistical methods, that are based on the assumption of homogeneity, are therefore inconsistent for all the length scales probed until now. A new, more general, conceptual framework is necessary to identify the real physical properties of these structures. In the range of self-similarity theories should shift from "amplitudes" to "exponents". The new analysis shows that all the available data are consistent with each other and show fractal correlations (with dimension $D \simeq 2$) up to the deepest scales probed until now ($1000 h^{-1}Mpc$) and even more as indicated from the new interpretation of the number counts. The distribution of luminous matter in the universe is therefore fractal and not homogeneous. The evidence for this being very strong up to $150 h^{-1}Mpc$ due to the statistical roboustness of the data and progressively weaker (statistically) at larger distances due to the limited data. These facts lead to fascinating conceptual implications about our knowledge of the universe and to a new scenario for the theoretical challenge in this field.

We introduce and study two new concepts which are essential for the quantitative analysis of the statistical quality of the available galaxy samples. These are the dilution effect and the small scale fluctuations. We show that the various data that are considered as pointing to a homogenous distribution are all affected by these spurious effects and their interpretation should be completely changed. In particular, we show that finite size effects strongly affect the determination of the galaxy number counts, namely the number versus magnitude relation ($N(<m)$) as computed from the origin. When one computes $N(<m)$ averaged over all the points of a redshift survey one observes an exponent $\alpha = D/5 \approx 0.4$ compatible with the fractal dimension $D \approx 2$ derived from the full correlation analysis. Instead the observation of an exponent $\alpha \approx 0.6$ at relatively small scales, where the distribution is certainly not homogeneous, is shown to be related to finite size effects. We conclude therefore that the observed counts correspond to a fractal distribution with dimension $D \approx 2$ in the entire range $12 \ltapprox m \ltapprox 28$, that is to say the largest scales ever probed for luminous matter. In addition our results permit to clarify various problems of the angular catalogs, and to show their compatibility with the fractal behavior. We consider also the distribution of Radio-galaxies, Quasars and $\gamma$ ray burst, and we show their compatibility with a fractal structure with $D \approx 1.6 \div 1.8$. Finally we have established a quantitative criterion that allows us to define and {\em predict} the statistical validity of a galaxy catalog (angular or three dimensional).

We discuss how luminosity and space distribution of galaxies are naturally linked in view of their multifractal properties. In particular we show that the mass (luminosity) function corresponding to a multifractal distribution in a given observed volume, consists of a power law followed by an exponential cut-off. This implies that the amplitude $\phi^*$ of the Schechter function scales with the sample depth, as confirmed by various observational data. This effect is analogous to the scaling of the space density due to the fractal nature of the space distribution. This has the important consequence that the luminosity function can be properly defined only in a volume limited sample. Also the so-called "luminosity segregation" and the concept of bias correspond to a natural consequence of multifractality. This implies however that they should be considered from a different perspective with respect to the usual one. Such a concept allows us to unify the space and luminosity distributions as being shaped by a single cause:multifractality which should therefore claim a central stage in theoretical investigations.

All the recent redshift surveys show highly irregular patterns of galaxies on scales of hundreds of megaparsecs such as chains, walls and cells. One of the most powerful catalog of galaxies is represented by the LEDA database that contains more than 36,000 galaxies with redshift. We study the correlation properties of such a samplefinding that galaxy distribution shows well defined fractal nature up to R$_s \sim 150 h^{-1}Mpc $ with fractal dimension $D \approx 2$. We test the consistency of these results versus the incompleteness in the sample.

We discuss the problem of galaxy correlations by considering the various methods by which this information can be obtained. We focus in particular on the volume limited three dimensional samples and discuss a new way to increase the scale of their statistical validity. From our previous results and the most recent ones on all the available catalogues we conclude that galaxy correlations show fractal properties with dimension $ D \approx 2$ up to the present observational limits ($200-800 \hmp$ depending on the catalogue) without any tendency towards homogenization. We also discuss the compatibility of this result with the galaxy counts as a function of magnitude and the angular catalogues. A new picture emerges which changes the standard ideas about the properties of the universe and requires a corresponding change in the theoretical concepts that one should use to describe it.

We study the behaviour of the power spectrum (PS) in the case of fractal structures. We show that in this case the main observational features of the PS, the large scale flattening and the scaling of the amplitude with sample depth, are related to finite size effects, due to the fractal nature of galaxy distribution in the sample. Comparing with the recent results of the PS for the CfA2 we conclude that this catalog is consistent with the fractal description.

We analyze the spatial and the luminosity properties of the Perseus-Pisces redshift survey. We find that the two point correlation function (CF) $\Gamma(r)$ is a power law up to the sample effective depth ($\sim 30 h^{-1}Mpc$), showing the fractal nature of the galaxy distribution in this catalog. The fractal dimension turns out to be $D \approx 2$. We also consider the CF $\xi(r)$ and in particular the behavior of the "{\it correlation length}" $r_0$ ($\xi(r_{0})\equiv 1$) as function of the sample size. In this respect we find, unambiguously, that the luminosity segregation effect is not supported by any experimental evidence. In addition we have studied the galaxian number-density ($n(r)$) and number-counts ($N(m)$) in the VL subsamples finding a good agreement with the properties of a fractal distribution. In particular our conclusion is that the $n(r)$ relation allows to extend the analysis of the fractal nature up to a deeper depth than that reached by the CF analysis, and, we find evidence for fractal properties up to the limiting depth of $120-140 h^{-1} Mpc$. We stress that the $N(m)$ relation must be studied in VL subsamples rather than in magnitude limited ones as usually done. Finally we have considered the correlations between galaxy positions and luminosities by means of the multifractal analysis. We find clear evidence for self-similar behavior of the whole luminosity-space distribution. These results confirm and extend those of Coleman \& Pietronero (1992).

A crucial issue in cosmology is the determination of the fluctuation power spectrum.The standard picture of the matter clustering, the Cold Dark Matter model (and its variant), assumes that,on scales smaller than a certain ``flattening scale'' $\lambda_f$, the power spectrum increases with the scale, while on much larger scales it decreases so to match the tiny fluctuations observed in the microwave background. The standard picture also assumes that, once a statistically homogeneous sample is reached, the power spectrum amplitude is fixed, and any major variation should be attributed to luminosity segregation. However, the determination of $\lambda_f$ and of the absolute amplitude, if any, is still matter of debate. In particular, there is no consensus on whether the turnaround has been detected or not, and on the actual importance of the luminosity segregation effect. We show that, due to the finiteness of the sample the standard analysis of self-similar (fractal) distributions yields a turnaround for scales close to the survey scale, and a systematic amplitude shift with the survey scale. We point out that both features, bending and scaling, are in agreement with recent determination of the power spectrum, in particular with the CfA2 spectrum. We remark that the standard power spectrum is not a well defined statistical tool to characterize the galaxy distribution when a homogeneity scale has not been reached. In order to perform an analysis that does not imply any a priori assumption one should study the PS of the density, rather than the PS of the density contrast.

A distribution of points that satisfies the property of local isotropy is not necessarily homogeneous: homogeneity is implied by the condition of local isotropy together with the assumption of analyticity or regularity. Here we show that the evidence of dipole saturation in galaxies (and clusters) catalogues, together with a monotone growth of the monopole, is an evidence of isotropy but not of homogeneity. This is fully compatible with a fractal structure which has the property of local isotropy, but it is non-analytic and non-homogeneous.

Galaxies and clusters distributions show two major properties: (i) the positions of galaxies and clusters are characterized by a power law distribution indicating properties with respect to their positions. (ii) The distribution of masses is also characterized by a power law corresponding to self similarity of different nature. These two properties are naturally unified by the concept of multifractality. This concept naturally arises if the distribution of matter, as given by both positions and masses, has self similar properties. We discuss the experimental situation in this respect and we introduce a simple stochastic model based of the aggregation process of particles The aim of this model is to understand which characteristic properties of the aggregation probability can gives rise to multifractal distribution In particular we find that a crucial element in this respect is that the aggregation probability should depend on the environment of the aggregation region.

The basic hypothesis of a post-Copernican Cosmological theory is that {\em all the points} of the Universe have to be essentially equivalent: this hypothesis is required in order to avoid any privileged {\em observer}. This assumption has been implemented by Einstein in the so-called Cosmological Principles (CP): {\em all the positions} in the Universe have to be essentially equivalent, so that the Universe is homogeneous. This situation implies also the condition of spherical symmetry about every point, so that the Universe in also Isotropic. There is a hidden assumption in the formulation of the CP in regard to the hypothesis that all the points are equivalent In fact, the condition that all the occupied points are statistically equivalent with respect to their environment correspond to the property of of Local Isotropy. It is generally believed that the Universe cannot be isotropic about every point without being also homogeneous [1]. Actually Local Isotropy does not necessarily implies homogeneity; in fact a topology theorem states that homogeneity is implied by the condition of local isotropy together with {\em the assumption of the analyticity or regularity} for the distribution of matter. Up to the seventies analyticity was an obvious implicit assumption in any physical problem. Recently however we have learned about intrinsically irregular structures for which analyticity should be considered as a property to be tested with appropriate analysis of experiment.

Two main features of the observable distribution of visible matter are the space correlations of galaxy positions and the mass function of galaxies. As discussed in Pietronero and Sylos Labini on this issue ([1], see also [2],[3]), the concept of multifractal (MF) distribution (including the masses) naturally unifies these two properties. Hence with the knowledge of the whole MF spectrum one obtains information on the correlations in space as well as on the mass function. From a theoretical point of view one would like to identify the dynamical process that leads to such a MF distribution. In order to gain some insight into this complex problem we have developed a simple stochastic model that includes some of the basic properties of aggregation process and allows us to pose a variety of interesting question concerning the possible dynamical origin of the MF distribution. The dynamics is characterised by some parameters, that have a direct physical meaning in term of cosmological processes. In this way we can relate the input parameters of the dynamics to the properties of the final configuration and produce a sort of phase diagram.