We provide a complete picture to the selfgravitating non-relativistic gas at thermal equilibrium using Monte Carlo simulations, analytic mean field methods (MF) and low density expansions. The system is shown to possess an infinite volume limit in the grand canonical (GCE), canonical (CE) and microcanonical (MCE) ensembles when(N, V) --> infinity, keeping N/ V^{1/3} fixed. We compute the equation of state (we do not assume it as is customary), as well as the energy, free energy, entropy, chemical potential, specific heats, compressibi- lities and speed of sound;we analyze their properties, signs and singularities. All physical quantities turn out to depend on a single variable eta = G m^2 N/ [V^{1/3} T] that is kept fixed in the N--> infinity and V --> infinity limit. The system is in a gaseous phase for eta < eta_T and collapses into a dense objet for eta > \eta_T in the CE with the pressure becoming large and negative. At eta simeq eta_T the isothermal compressibility diverges. Our Monte Carlo simulations yield eta_T simeq 1.515. PV/[NT] = f(eta) and all physical magni- tudes exhibit a square root branch point at eta = eta_C > eta_T. The MF for spherical symmetry yields eta_C = 1.561764.. while Monte Carlo on a cube yields eta_C simeq 1.540.The function f(eta) has a second Riemann sheet which is only physically realized in the MCE.In the MCE, the collapse phase transition takes place in this second sheet near eta_MC = 1.26 and the pressure and temperature are larger in the collapsed phase than in the gas phase.Both collapse phase transitions (CE and MCE) are of zeroth order since the Gibbs free energy jumps at the transitions. f(eta), obeys in MF a first order non-linear differential equation of first kind Abel's type.The MF gives an extremely accurate picture in agreement with Monte Carlo both in the CE and MCE.

We complete our study of the selfgravitating gas by computing the fluctuations around the saddle point solution for the three statistical ensembles. Although the saddle point is the same for the three ensembles, the fluctuations change from one ensemble to another. The zeroes of the fluctuations determinant determine the position of the critical points for each ensemble.This yields the domains of validity of the mean field approach. Only the S wave determinant exhibits critical points.Closed formulae for the S and P wave determinants are derived.The local properties of the selfgravitating gas in thermodynamic equilibrium are studied in detail. The pressure, energy density, particle density and speed of sound are computed and analyzed as functions of the position. The equation of state turns out to be locally p(r) = T rho(r) as for the ideal gas. Starting from the partition function of the selfgravitating gas, we prove in this microscopic calculation that the hydrostatic description yielding locally the ideal gas equation of state is exact in the N = infinity limit. The dilute nature of the thermodynamic limit(N, L -> infinity with N/L fixed) together with the long range nature of the gravitational forces play a crucial role in the obtention of such ideal gas equation. The self-gravitating gas being inhomogeneous, we have PV/[NT] = f(eta) leq 1 for any finite volume V. The inhomogeneous particle distribution in the ground state suggests a fractal distribution with Haussdorf dimension D, D is slowly decreasing with increasing density, 1 lesssim D < 3. The average distance between particles is computed in Monte Carlo simulations andanalytically in the mean field approach. A dramatic drop at the phase transition is exhibited, clearly illustrating the properties of the collapse.

We provide a complete picture of the self-gravitating non-relativistic gas at thermal equilibrium using Monte Carlo simulations (MC), analytic mean field methods (MF) and low density expansions. The system is shown to possess an infinite volume limit, both in the canonical (CE) and in the microcanonical ensemble (MCE) when N, V \to \infty, keeping N/ V^{1/3} fixed. We {\bf compute} the equation of state (we do not assume it as is customary), the entropy, the free energy, the chemical potential, the specific heats, the compressibilities, the speed of sound and analyze their properties, signs and singularities. The MF equation of state obeys a {\bf first order} non-linear differential equation of Abel type. The MF gives an accurate picture in agreement with the MC simulations both in the CE and MCE. The inhomogeneous particle distribution in the ground state suggest a fractal distribution with Haussdorf dimension D with D slowly decreasing with increasing density, 1 \lesssim D < 3.

The density of mass levels \rho(m) and the critical temperature for strings in de Sitter space-time are found. QFT and string theory in de Sitter space are compared. A `Dual'-transform is introduced which relates classical to quantum string lengths, and more generally, QFT and string domains. Interestingly, the string temperature in De Sitter space turns out to be the Dual transform of the QFT-Hawking-Gibbons temperature. The back reaction problem for strings in de Sitter space is addressed selfconsistently in the framework of the `string analogue' model (or thermodynamical approach), which is well suited to combine QFT and string study.We find de Sitter space-time is a self-consistent solution of the semiclassical Einstein equations in this framework. Two branches for the scalar curvature R(\pm) show up: a classical, low curvature solution (-), and a quantum high curvature solution (+), enterely sustained by the strings. There is a maximal value for the curvature R_{\max} due to the string back reaction. Interestingly, our Dual relation manifests itself in the back reaction solutions: the (-) branch is a classical phase for the geometry with intrinsic temperature given by the QFT-Hawking-Gibbons temperature.The (+) is a stringy phase for the geometry with temperature given by the intrinsic string de Sitter temperature. 2 + 1 dimensions are considered, but conclusions hold generically in D dimensions.

The self-gravitating thermal gas (non-relativistic particles of mass m at temperature T) is exactly equivalent to a field theory with a single scalar field phi(x) and exponential self-interaction. We build up perturbation theory around a space dependent stationary point phi_0(r) in a finite size domain delta \leq r \leq R ,(delta << R), which is relevant for astrophysical applica- tions (interstellar medium,galaxy distributions).We compute the correlations of the gravitational potential (phi) and of the density and find that they scale; the latter scales as 1/r^2. A rich structure emerges in the two-point correl- tors from the phi fluctuations around phi_0(r). The n-point correlators are explicitly computed to the one-loop level.The relevant effective coupling turns out to be lambda=4 pi G m^2 / (T R). The renormalization group equations (RGE) for the n-point correlator are derived and the RG flow for the effective coupling lambda(tau) [tau = ln(R/delta), explicitly obtained.A novel dependence on tau emerges here.lambda(tau) vanishes each time tau approaches discrete values tau=tau_n = 2 pi n/sqrt7-0, n=0,1,2, ...Such RG infrared stable behavior [lambda(tau) decreasing with increasing tau] is here connected with low density self-similar fractal structures fitting one into another.For scales smaller than the points tau_n, ultraviolet unstable behaviour appears which we connect to Jeans' unstable behaviour, growing density and fragmentation. Remarkably, we get a hierarchy of scales and Jeans lengths following the geometric progression R_n=R_0 e^{2 pi n /sqrt7} = R_0 [10.749087...]^n . A hierarchy of this type is expected for non-spherical geometries,with a rate different from e^{2 n/sqrt7}.

Fractal structures are observed in the universe in two very different ways. Firstly, in the gas forming the cold interstellar medium in scales from 10^{-4} pc till 100 pc. Secondly, the galaxy distribution has been observed to be fractal in scales up to hundreds of Mpc. We give here a short review of the statistical mechanical (and field theoretical) approach developed by us. We consider a non-relativistic self-gravitating gas in thermal equilibrium at temperature T inside a volume V. The statistical mechanics of such system has special features and, as is known, the thermodynamical limit does not exist in its customary form. Moreover, the treatments through microcanonical, canonical and grand canonical ensembles yield different results.We present here for the first time the equation of state for the self-gravitating gas in the canonical ensemble. We find that it has the form p = [N T/ V] f(eta), where p is the pressure, N is the number of particles and \eta \equiv {G m^2 N \over V^{1/3} T}. The N \to\infty and V \to\infty limit exists keeping \eta fixed. We compute the function f(\eta) using Monte Carlo simulations and for small eta analytically. We compute the thermodynamic quantities of the system as free energy, entropy, chemical potential, specific heat, compressibility and speed of sound. We reproduce the well-known gravitational phase transition associated to the Jeans' instability. Namely, a gaseous phase for eta < eta_c and a condensed phase for eta > eta_c. Moreover, we derive the precise behaviour of the physical quantities near the transition. In particular, the pressure vanishes as p \sim(eta_c-eta)^B with B \sim 0.2 and eta_c \sim 1.6 and the energy fluctuations diverge as \sim(eta_c-eta)^{B-1}. The speed of sound decreases monotonically and approaches the value sqrt{T/6} at the transition.

We develop a field theoretical approach to the cold interstellar medium (ISM). We show that a non-relativistic self-gravitating gas in thermal equilibrium with variable number of atoms or fragments is exactly equivalent to a field theory of a single scalar field \phi({\vec x}) with exponential self interaction.We analyze this field theory perturbatively and non-perturbatively through the renormalization group approach.We show scaling behaviour(critical) for a continuous range of the temperature and of the other physical parameters. We derive in this framework the scaling relation Delta M(R) \sim R^{d_H} for the mass on a region of size R, and Delta v \sim R^q for the velocity dispersion where q = (d_H -1)/2. For the density-density correlations we find a power-law behaviour for large distances \sim|{\vec r_1} -{\vec r_2}|^{2 d_H -6}. The fractal dimension d_H turns to be related with the critical exponent nu of the correlation length by d_H = 1/nu. The renormalization group approach for a single component scalar field in three dimensions states that the long- distance critical behaviour is governed by the (non-perturbative) Ising fixed point. The corresponding values of the scaling exponents are nu = 0.631..., d_H = 1.585... and q = 0.293.... Mean field theory yields for the scaling exponents nu =1/2, d_H = 2 and q = 1/2. Both the Ising and the mean field values are compatible with the present ISM observational data: 1.4 \leq d_H \leq 2, 0.3 \leq q \leq 0.6.As typical in critical phenomena, the scaling behaviour and critical exponents of the ISM can be obtained without dwelling into the dynamical (time-dependent) behaviour. The relevant r\^ole of self- gravity is stressed by the authors in a Letter to Nature, September 5, 1996.