The $\Lambda$CDM cosmological scenario predicts that our Galaxy should contain hundreds of stellar streams at the solar vicinity, fossil relics of the merging history of the Milky Way and more generally of the hierarchical growth of galaxies. Because of the mixing time scales in the inner Galaxy, it has been claimed that these streams should be difficult to detect in configuration space but can still be identifiable in kinematic-related spaces like the energy/angular momenta spaces, E-Lz and Lperp-Lz, or spaces of orbital/velocity parameters. By means of high-resolution, dissipationless N-body simulations, containing between 25$\times10^6$ and 35$\times10^6$ particles, we model the accretion of a series of up to four 1:10 mass ratio satellites then up to eight 1:100 satellites and we search systematically for the signature of these accretions in these spaces. In all spaces considered (1) each satellite gives origin to several independent overdensities; (2) overdensities of multiple satellites overlap; (3) satellites of different masses can produce similar substructures; (4) the overlap between the in-situ and the accreted population is considerable everywhere; (5) in-situ stars also form substructures in response to the satellite(s) accretion. These points are valid even if the search is restricted to kinematically-selected halo stars only. As we are now entering the 'Gaia era', our results warn that an extreme caution must be employed before interpreting overdensities in any of those spaces as evidence of relics of accreted satellites. Reconstructing the accretion history of our Galaxy will require a substantial amount of accurate spectroscopic data, that, complemented by the kinematic information, will possibly allow us to (chemically) identify accreted streams and measure their orbital properties. (abridged)

We present a class of models where both the primordial inflation and the late times de Sitter phase are driven by simple phenomenological agegraphic potentials. In this context, a possible new scenario for a smooth exit from inflation to the radiation era is discussed by resorting the kination (stiff) era but without the inefficient radiation production mechanism of these models. This is done by considering rapidly decreasing expressions for $V(t)$ soon after inflation. We show that the parameters of our models can reproduce the scalar spectral parameter $n_s$ predicted by Planck data in particular for models with concave potentials. Finally, according to the recent BICEP2 data, all our models allow a huge amount of primordial gravitational waves.

The distribution of matter in the universe shows a complex pattern, formed by cluster of galaxies, voids and filaments denoted as cosmic web. Different approaches have been proposed to model such structure in the framework of the general relativity. Recently, one of us has proposed a generalization ($\Lambda$FB model) of the Fractal Bubble model, proposed by Wiltshire, which accounts for such large scale structure. The $\Lambda$FB model is an evolution of FB model and includes in a consistent way a description of inhomogeneous matter distribution and a $\Lambda$ term. Here we analyze the $\Lambda$FB model focusing on the relation between cosmological parameters. The main result is the consistency of $\Lambda$CDM model values for $\Omega_{\Lambda0}$ ($\approx 0.7$) and $\Omega_{k0}$ ($|\Omega_{k0}|<\approx 0.01$) with a large fraction of voids. This allows to quantify to which extent the inhomogeneous structure could account for $\Lambda$ constant consistently with standard values of the other cosmological parameters.

In this paper we report a study (see Mastrobuono-Battisti et al., 2012) about the formation and characteristics of the tidal tails around Palomar 5 along its orbit in the Milky Way potential, by means of direct N-body simulations and simplified numerical models. Unlike previous findings, we are able to reproduce the substructures observed in the stellar streams of this cluster, without including any lumpiness in the dark matter halo. We show that overdensities similar to those observed in Palomar 5 can be reproduced by the epicyclic motion of stars along its tails, i.e. a simple local accumulation of orbits of stars that escaped from the cluster with very similar positions and velocities. This process is able to form stellar clumps at distances of several kiloparsecs from the cluster, so it is not a phenomenon confined to the inner part of Palomar 5's tails, as previously suggested.

By means of direct N-body simulations and simplified numerical models, we study the formation and characteristics of the tidal tails around Palomar 5, along its orbit in the Milky Way potential. Unlike previous findings, we are able to reproduce the substructures observed in the stellar streams of this cluster, without including any lumpiness in the dark matter halo. We show that overdensities similar to those observed in Palomar 5 can be reproduced by the epicyclic motion of stars along its tails, i.e. a simple local accumulation of orbits of stars that escaped from the cluster with very similar positions and velocities. This process is able to form stellar clumps at distances of several kiloparsecs from the cluster, so it is not a phenomenon confined to the inner part of Palomar 5's tails, as previously suggested. Our models can reproduce the density contrast between the clumps and the surrounding tails found in the observed streams, without including any lumpiness in the dark halo, suggesting new upper limits on its granularity.

We present the results of a detailed analysis of the projected velocity dispersion of the globular cluster Palomar 14 performed using recent high-resolution spectroscopic data and extensive Monte Carlo simulations. The comparison between the data and a set of dynamical models (differing in fraction of binaries, degree of anisotropy, mass-to-light ratio M/L, cluster orbit and theory of gravity) shows that the observed velocity dispersion of this stellar system is well reproduced by Newtonian models with a fraction of binaries f_b<30% and a M/L compatible with the predictions of stellar evolution models. Instead, models computed with a large fraction of binaries systematically overestimate the cluster velocity dispersion. We also show that, across the parameter space sampled by our simulations, models based on the Modified Newtonian Dynamics theory can be reconciled with observations only assuming values of M/L lower than those predicted by stellar evolution models under standard assumptions.

It has been known since many decades that galaxy interactions can induce star formation (hereafter SF) enhancements and that one of the driving mechanisms of this enhancement is related to gas inflows into the central galaxy regions, induced by asymmetries in the stellar component, like bars. In the last years many evidences have been accumulating, showing that interacting pairs have central gas-phase metallicities lower than those of field galaxies, by {\sim} 0.2-0.3 dex on average. These diluted ISM metallicities have been explained as the result of inflows of metal-poor gas from the outer disk to the galaxy central regions. A number of questions arises: What's the timing and the duration of this dilution? How and when does the SF induced by the gas inflow enrich the circumnuclear gas with re-processed material? Is there any correlation between the timing and strength of the dilution and the timing and intensity of the SF? By means of Tree-SPH simulations of galaxy major interactions, we have studied the effect that gas inflows have on the ISM dilution, and the effect that the induced SF has, subsequently, in re-enriching the nuclear gas. In this contribution, we present the main results of this study.

Strong inflows of gas from the outer disk to the inner kiloparsecs are induced during the interaction of disk galaxies. This inflow of relatively low-metallicity gas dilutes the metallicity of the circumnuclear gas. We have investigated several aspects of the process as the timing and duration of the dilution and its correlation with the induced star formation. We analysed major (1:1) gas-rich interactions and mergers, spanning a range of initial orbital characteristics. Star formation and metal enrichment from SNe are included in our model. Our results show that the strongest trend is between the star formation rate and the dilution of the metals in the nuclear region; i.e., the more intense the central burst of star formation, the more the gas is diluted. This trend comes from strong inflows of relatively metal-poor gas from the outer regions of both disks, which fuels the intense star formation and lowers the overall metallicity for a time. The strong inflows happen on timescales of about 10^8 years or less, and the most intense star formation and lowest gas phase metallicities are seen generally after the first pericentre passage. As the star formation proceeds and the merger advances, the dilution reduces and enrichment becomes dominant - ultimately increasing the metallicity of the circumnuclear gas to a level higher than the initial metallicities of the merging galaxies. The "fly-bys" - pairs that interact but do not merge - also cause some dilution. We even see some dilution early in the merger or in the "fly-bys" and thus do not observe a strong trend between the nuclear metallicities and separation in our simulations until the merger is well advanced. We also analyse the O and Fe enrichment of the ISM, and show that the evolution of the alpha/Fe ratios, as well as the dilution of the central gas metallicity, can be used as a clock for "dating" the interaction.

We present the results of detailed N-body simulations of clusters moving in a realistic Milky Way (MW) potential. The strong interaction with the bulge and the disk of the Galaxy leads to the formation of tidal tails, emanating from opposite sides of the cluster. Their orientation and morphology may be interpreted easily in terms of a comoving frame of coordinates.

We present the results of detailed N-body simulations of clusters moving in a realistic Milky Way (MW) potential. The strong interaction with the bulge and the disk of the Galaxy leads to the formation of tidal tails, emanating from opposite sides of the cluster. Some characteristic features in the morphology and orientation of these streams are recognized and intepreted. The tails have a complex morphology, in particular when the cluster approaches its apogalacticon, showing multiple ``arms'' in remarkable similarity to the structures observed around NGC 288 and Willman 1. Actually, the tails are generally good tracers of the cluster path quite far from the cluster center (>7--8 tidal radii), while on the smaller scale they are mainly pointing in the direction of the Galaxy center. In particular, the orientation of the inner part of the tails is highly correlated to the cluster orbital phase and to the local orbital angular acceleration. This implies that, in general, the orbital path cannot be estimated directly from the orientation of the tails, unless a sufficient large field around the cluster is available.

In the last decade, observational studies have shown the existence of tidal streams in the outer part of many galactic globular clusters. The most striking examples of clusters with well defined tidal tails are represented by Palomar 5 and NGC 5466 (both observed in the framework of the Sloan Digital Sky Survey), which show structures elongated for 4 kpc and 1 kpc in length, respectively. Unfortunately, most of the observational studies about globular clusters (GCs) do not cover such a large field of the sky as the SDSS does. In this framework, by mean of a parallel, adaptive tree-code, we performed detailed N-body simulations of GCs moving in a realistic three-components (bulge, disk and halo) Milky Way potential, in order to clarify whether and to what extent tails in the clusters outer regions (few tidal radii) are tracers of the local orbits and, also, if some kind of correlation exists among the cluster orbital phase and the orientation of such streams.

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.

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 present a cosmological model in which the metric allows for an inhomogeneous Universe with no intrinsic symmetries (Stephani models), providing the ideal features to describe a fractal distribution of matter. Constraints on the metric functions are derived using the expansion and redshift relations and allowing for scaling number counts, as expected in a fractal set. The main characteristics of such a cosmological model are discussed.

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.

Erosion by flowing water is one of the major forces shaping the surface of Earth. Studies in the last decade have shown, in particular, that the drainage region of rivers, where water is collected, exhibits scale invariant features characterized by exponents that are the same for rivers around the world. Here we show that from the data obtained by the MOLA altimeter of the Mars Global Surveyor one can perform the same analysis for mountain sides on Mars. We then show that in some regions fluid erosion might have played a role in the present martian landscape.

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

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 study the angular correlations of various galaxy catalogs (CfA1, SSRS1, Perseus-Pisces, APM Bright Galaxies and Zwicky). We find that the angular correlation exponent is $\gamma_a= 0.1 \pm 0.1$ rather than $\gamma_a=0.7$ as usually found by the standard correlation function $\omega(\theta)$. We identify the problem in the artificial decay of $\omega(\theta)$. Moreover we find that no characteristic angular scale is present in any of the analyzed catalogs. Finally we show that all the available data are consistent with each other and the angular distribution of galaxies is quite naturally compatible with a fractal structure with $D \approx 2$.

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

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