Metis is the first solar coronagraph designed for a space mission capable of performing simultaneous imaging of the off-limb solar corona in both visible and UV light. The observations obtained with Metis aboard the Solar Orbiter ESA-NASA observatory will enable us to diagnose, with unprecedented temporal coverage and spatial resolution, the structures and dynamics of the full corona from 1.7 $R_\odot$ to about 9 $R_\odot$. Due to the uniqueness of the Solar Orbiter mission profile, Metis will be able to observe the solar corona from a close vantage point (down to 0.28 AU), achieving out-of-ecliptic views with the increase of the orbit inclination over time. Moreover, observations near perihelion, during the phase of lower rotational velocity of the solar surface relative to the spacecraft, will allow longer-term studies of the coronal features. Thanks to a novel occultation design and a combination of a UV interference coating of the mirrors and a spectral bandpass filter, Metis images the solar corona simultaneously in the visible light band, between 580 and 640 nm, and in the UV H I Lyman-{\alpha} line at 121.6 nm. The coronal images in both the UV Lyman-{\alpha} and polarised visible light are obtained at high spatial resolution with a spatial scale down to about 2000 km and 15000 km at perihelion, in the cases of the visible and UV light, respectively. A temporal resolution down to 1 second can be achieved when observing coronal fluctuations in visible light. The Metis measurements will allow for complete characterisation of the main physical parameters and dynamics of the electron and neutral hydrogen/proton plasma components of the corona in the region where the solar wind undergoes acceleration and where the onset and initial propagation of coronal mass ejections take place, thus significantly improving our understanding of the region connecting the Sun to the heliosphere.

We explore the formation and relaxation of so-called quasi-stationary states (QSS) for particle distributions in three dimensions interacting via an attractive radial pair potential $V(r \rightarrow \infty) \sim 1/r^\gamma$ with $\gamma > 0$, and either a soft-core or hard-core regularization at small $r$. In the first part of the paper we generalize, for any spatial dimension $d \geq 2$, Chandrasekhar's approach for the case of gravity to obtain analytic estimates of the rate of collisional relaxation due to two body collisions. The resultant relaxation rates indicate an essential qualitative difference depending on the integrability of the pair force at large distances: for $\gamma >d-1$ the rate diverges in the large particle number $N$ (mean field) limit, unless a sufficiently large soft core is present; for $\gamma < d-1$, on the other hand, the rate vanishes in the same limit even in the absence of any regularization. In the second part of the paper we compare our analytical predictions with the results of extensive parallel numerical simulations in $d=3$, for a range of different exponents $\gamma$ and soft cores leading to the formation of QSS. We find, just as for the previously well studied case of gravity (which we also revisit), excellent agreement between the parametric dependence of the observed relaxation times and our analytic predictions. Further, as in the case of gravity, we find that the results indicate that, when large impact factors dominate, the appropriate cut-off is the size of the system (rather than, for example, the mean inter-particle distance). Our results provide strong evidence that the existence of QSS is robust only for long-range interactions with a large distance behavior $\gamma < d-1$; for $\gamma \geq d-1$ the existence of such states will be conditioned strongly on the short range properties of the interaction.

The purpose of our study is to understand the mathematical origin in real space of modulated and damped sinusoidal peaks observed in cosmic microwave background radiation anisotropies. We use the theory of the Fourier transform to connect localized features of the two-point correlation function in real space to oscillations in the power spectrum. We also illustrate analytically and by means of Monte Carlo simulations the angular correlation function for distributions of filled disks with fixed or variable radii capable of generating oscillations in the power spectrum. While the power spectrum shows repeated information in the form of multiple peaks and oscillations, the angular correlation function offers a more compact presentation that condenses all the information of the multiple peaks into a localized real space feature. We have seen that oscillations in the power spectrum arise when there is a discontinuity in a given derivative of the angular correlation function at a given angular distance. These kinds of discontinuities do not need to be abrupt in an infinitesimal range of angular distances but may also be smooth, and can be generated by simply distributing excesses of antenna temperature in filled disks of fixed or variable radii on the sky, provided that there is a non-null minimum radius and/or the maximum radius is constrained.

"Quasi-stationary" states are approximately time-independent out of equilibrium states which have been observed in a variety of systems of particles interacting by long-range interactions. We investigate here the conditions of their occurrence for a generic pair interaction V(r \rightarrow \infty) \sim 1/r^a with a > 0, in d>1 dimensions. We generalize analytic calculations known for gravity in d=3 to determine the scaling parametric dependences of their relaxation rates due to two body collisions, and report extensive numerical simulations testing their validity. Our results lead to the conclusion that, for a < d-1, the existence of quasi-stationary states is ensured by the large distance behavior of the interaction alone, while for a > d-1 it is conditioned on the short distance properties of the interaction, requiring the presence of a sufficiently large soft-core in the interaction potential.

We formalize a classification of pair interactions based on the convergence properties of the {\it forces} acting on particles as a function of system size. We do so by considering the behavior of the probability distribution function (PDF) P(F) of the force field F in a particle distribution in the limit that the size of the system is taken to infinity at constant particle density, i.e., in the "usual" thermodynamic limit. For a pair interaction potential V(r) with V(r) \rightarrow \infty) \sim 1/r^a defining a {\it bounded} pair force, we show that P(F) converges continuously to a well-defined and rapidly decreasing PDF if and only if the {\it pair force} is absolutely integrable, i.e., for a > d-1, where d is the spatial dimension. We refer to this case as {\it dynamically short-range}, because the dominant contribution to the force on a typical particle in this limit arises from particles in a finite neighborhood around it. For the {\it dynamically long-range} case, i.e., a \leq d-1, on the other hand, the dominant contribution to the force comes from the mean field due to the bulk, which becomes undefined in this limit. We discuss also how, for a \leq d-1 (and notably, for the case of gravity, a=d-2) P(F) may, in some cases, be defined in a weaker sense. This involves a regularization of the force summation which is generalization of the procedure employed to define gravitational forces in an infinite static homogeneous universe. We explain that the relevant classification in this context is, however, that which divides pair forces with a > d-2 (or a < d-2), for which the PDF of the {\it difference in forces} is defined (or not defined) in the infinite system limit, without any regularization. In the former case dynamics can, as for the (marginal) case of gravity, be defined consistently in an infinite uniform system.

The NEMO Collaboration installed and operated an underwater detector including prototypes of the critical elements of a possible underwater km3 neutrino telescope: a four-floor tower (called Mini-Tower) and a Junction Box. The detector was developed to test some of the main systems of the km3 detector, including the data transmission, the power distribution, the timing calibration and the acoustic positioning systems as well as to verify the capabilities of a single tridimensional detection structure to reconstruct muon tracks. We present results of the analysis of the data collected with the NEMO Mini-Tower. The position of photomultiplier tubes (PMTs) is determined through the acoustic position system. Signals detected with PMTs are used to reconstruct the tracks of atmospheric muons. The angular distribution of atmospheric muons was measured and results compared with Monte Carlo simulations.

The dynamics of infinite, asymptotically uniform, distributions of self-gravitating particles in one spatial dimension provides a simple toy model for the analogous three dimensional problem. We focus here on a limitation of such models as treated so far in the literature: the force, as it has been specified, is well defined in infinite point distributions only if there is a centre of symmetry (i.e. the definition requires explicitly the breaking of statistical translational invariance). The problem arises because naive background subtraction (due to expansion, or by "Jeans' swindle" for the static case), applied as in three dimensions, leaves an unregulated contribution to the force due to surface mass fluctuations. Following a discussion by Kiessling, we show that the problem may be resolved by defining the force in infinite point distributions as the limit of an exponentially screened pair interaction. We show that this prescription gives a well defined (finite) force acting on particles in a class of perturbed infinite lattices, which are the point processes relevant to cosmological N-body simulations. For identical particles the dynamics of the simplest toy model is equivalent to that of an infinite set of points with inverted harmonic oscillator potentials which bounce elastically when they collide. We discuss previous results in the literature, and present new results for the specific case of this simplest (static) model starting from "shuffled lattice" initial conditions. These show qualitative properties (notably its "self-similarity") of the evolution very similar to those in the analogous simulations in three dimensions, which in turn resemble those in the expanding universe.

We consider the construction of point processes from tilings, with equal volume tiles, of d-dimensional Euclidean space. We show that one can generate, with simple algorithms ascribing one or more points to each tile, point processes which are "superhomogeneous'' (or "hyperuniform''), i.e., for which the structure factor S(k) vanishes when the wavenumber k tends to zero. The exponent of the leading small-k behavior depends in a simple manner on the nature of the correlation properties of the specific tiling and on the conservation of the mass moments of the tiles. Assigning one point to the center of mass of each tile gives the exponent \gamma=4 for any tiling in which the shapes and orientations of the tiles are short-range correlated. Smaller exponents, in the range 4-d<\gamma<4 (and thus always superhomogeneous for d\leq 4), may be obtained in the case that the latter quantities have long-range correlations. Assigning more than one point to each tile in an appropriate way, we show that one can obtain arbitrarily higher exponents in both cases. We illustrate our results with explicit constructions using known deterministic tilings, as well as some simple stochastic tilings for which we can calculate S(k) exactly. Our results provide, we believe, the first explicit analytical construction of point processes with \gamma > 4. Applications to condensed matter physics, and also to cosmology, are briefly discussed.

We study how the two-point density correlation properties of a point particle distribution are modified when each particle is divided, by a stochastic process, into an equal number of identical "daughter" particles. We consider generically that there may be non-trivial correlations in the displacement fields describing the positions of the different daughters of the same "mother" particle, and then treat separately the cases in which there are, or are not, correlations also between the displacements of daughters belonging to different mothers. For both cases exact formulae are derived relating the structure factor (power spectrum) of the daughter distribution to that of the mother. These results can be considered as a generalization of the analogous equations obtained in ref. [1] (cond-mat/0409594) for the case of stochastic displacement fields applied to particle distributions. An application of the present results is that they give explicit algorithms for generating, starting from regular lattice arrays, stochastic particle distributions with an arbitrarily high degree of large-scale uniformity.

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 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 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 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.

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.

The effect of a stochastic displacement field on a statistically independent point process is analyzed. Stochastic displacement fields can be divided into two large classes: spatially correlated and uncorrelated. For both cases exact transformation equations for the two-point correlation function and the power spectrum of the point process are found, and a detailed study of them with important paradigmatic examples is done. The results are general and in any dimension. A particular attention is devoted to the kind of large scale correlations that can be introduced by the displacement field, and to the realizability of arbitrary ``superhomogeneous'' point processes.

A well known argument in cosmology gives that the power spectrum (or structure function) $P(k)$ of mass density fluctuations produced from a uniform initial state by physics which is causal (i.e. moves matter and momentum only up to a finite scale) has the behaviour $P(k) \propto k^4$ at small $k$. Noting the assumption of analyticity at $k=0$ of $P(k)$ in the standard derivation of this result, we introduce a class of solvable one dimensional models which allows us to study the relation between the behaviour of $P(k)$ at small $k$ and the properties of the probability distribution $f(l)$ for the spatial extent $l$ of mass and momentum conserving fluctuations. We find that the $k^4$ behaviour is obtained in the case that the first {\it six} moments of $f(l)$ are finite. Interestingly the condition that the fluctuations be localised - taken to correspond to the convergence of the first two moments of $f(l)$ - imposes only the weaker constraint $P(k) \propto k^n$ with $n$ anywhere in the range $0< n \leq 4$. We interpret this result to suggest that the causality bound will be loosened in this way if quantum fluctuations are permitted.

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 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.

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.

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.

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.

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.

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.

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.

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 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).