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

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.

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

Precise analyses of the statistical and scaling properties of galaxy distribution are essential to elucidate the large-scale structure of the universe. Given the ongoing debate on its statistical features, the development of statistical tools permitting to discriminate accurately different spatial patterns is highly desiderable. This is specially the case when non-fractal distributions have power law two-point correlation functions, which are usually signatures of fractal properties. Here we review some possible methods used in the literature and introduce a new variable called "scaling gradient". This tool and the conditional variance are shown to be effective in providing an unambiguous way for such a distinction. Their application is expected to be of outmost importance in the analysis of upcoming galaxy catalogues.

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.

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.

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.

We study gravitational clustering of mass points in three dimensions with random initial positions and periodic boundary conditions (no expansion) by numerical simulations. Correlation properties are well defined in the system and a sort of thermodynamic limit can be defined for the transient regime of cluste ring. Structure formation proceeds along two paths: (i) fluid-like evolution of density perturbations at large scales and (ii) shift of the granular (non fluid) properties from small to large scales. The latter mechanism finally dominates at all scales and it is responsible for the self-similar characteristics of the clustering.

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.

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.

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.

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.

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.

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.

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.

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

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.