We generalize the Thomas-Fermi approach to galaxy structure to include self-consistently and non-linearly central supermassive black holes. This approach naturally incorporates the quantum pressure of the warm dark matter (WDM) particles and shows its full powerful and clearness in the presence of supermassive black holes (SPMHs). We find the main galaxy and central black hole magnitudes: halo radius r_h , halo mass M_h, black hole mass M_BH, velocity dispersion, phase space density, with their realistic astrophysical values, masses and sizes over a wide galaxy range. The SMBH masses arise naturally in this framework. Our extensive numerical calculations and detailed analytic resolution show that with SMBH's, both WDM regimes: classical (Boltzmann dilute) and quantum (compact) do necessarily co-exist in any galaxy: from the smaller and compact galaxies to the largest ones. The transition from the quantum to the classical region occurs precisely at the same point r_A where the chemical potential vanishes. A novel halo structure with three regions shows up: A small quantum compact core of radius r_A around the SMBH, followed by a less compact region till the BH influence radius r_i, and then for r> r_i the known halo galaxy shows up with its astrophysical size. Three representative families of galaxy plus central SMBH solutions are found and analyzed:small, medium and large galaxies having SMBH masses of 10^5, 10^7 and 10^9 M_sun respectively. A minimum galaxy size and mass ~ 10^7 M_sun larger than the one without SMBH is found. Small galaxies in the range 10^4 M_sun < M_h < 10^7 M_sun cannot harbor central SMBHs. We find novel scaling M_BH - r_h - M_h relations. The galaxy equation of state is derived: The pressure P(r) takes huge values in the SMBH vecinity and then sharply decreases entering the classical region following a local perfect gas behaviour.(Abridged)

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

We find the distribution function f(E) for dark matter (DM) halos in galaxies and the corresponding equation of state from the (empirical) DM density profiles derived from observations. We solve for DM in galaxies the analogous of the Eddington equation originally used for the gas of stars in globular clusters. The observed density profiles are a good realistic starting point and the distribution functions derived from them are realistic. We do not make any assumption about the DM nature, the methods developed here apply to any DM kind, though all results are consistent with Warm DM. With these methods we find: (i) Cored density profiles behaving quadratically for small distances rho(r) r -> 0 = rho(0) - K r^2 produce distribution functions which are finite and positive at the halo center while cusped density profiles always produce divergent distribution functions at the center. (ii) Cored density profiles produce approximate thermal Boltzmann distribution functions for r < 3 r_h where r_h is the halo radius. (iii) Analytic expressions for the dispersion velocity and the pressure are derived yielding an ideal DM gas equation of state with local temperature T(r) = m v^2(r)/3. T(r) turns to be constant in the same region where the distribution function is thermal and exhibits the same temperature within the percent. The self-gravitating DM gas can thermalize despite being collisionless because it is an ergodic system. (iv) The DM halo can be consistently considered at local thermal equilibrium with: (a) a constant temperature T(r) = T_0 for r < 3 \; r_h, (b) a space dependent temperature T(r) for 3 r_h < r < R_{virial}, which slowly decreases with r. That is, the DM halo is realistically a collisionless self-gravitating thermal gas for r < R_{virial}. (v) T(r) outside the halo radius nicely follows the decrease of the circular velocity squared.

Light scalars (as the axion) with mass m ~ 10^{-22} eV forming a Bose-Einstein condensate (BEC) exhibit a Jeans length in the kpc scale and were therefore proposed as dark matter (DM) candidates. Our treatment here is generic, independent of the particle physics model and applies to all DM BEC, in or out of equilibrium. Two observed quantities crucially constrain DM in an inescapable way: the average DM density rho_{DM} and the phase-space density Q. The observed values of rho_{DM} and Q in galaxies today constrain both the possibility to form a BEC and the DM mass m. These two constraints robustly exclude axion DM that decouples just after the QCD phase transition. Moreover, the value m ~ 10^{-22} eV can only be obtained with a number of ultrarelativistic degrees of freedom at decoupling in the trillions which is impossible for decoupling in the radiation dominated era. In addition, we find for the axion vacuum misalignment scenario that axions are produced strongly out of thermal equilibrium and that the axion mass in such scenario turns to be 17 orders of magnitude too large to reproduce the observed galactic structures. Moreover, we also consider inhomogenous gravitationally bounded BEC's supported by the bosonic quantum pressure independently of any particular particle physics scenario. For a typical size R ~ kpc and compact object masses M ~ 10^7 Msun they remarkably lead to the same particle mass m ~ 10^{-22} eV as the BEC free-streaming length. However, the phase-space density for the gravitationally bounded BEC's turns to be more than sixty orders of magnitude smaller than the galaxy observed values. We conclude that the BEC's and the axion cannot be the DM particle. However, an axion in the mili-eV scale may be a relevant source of dark energy through the zero point cosmological quantum fluctuations.

Warm dark matter (WDM) means DM particles with mass m in the keV scale. For large scales, (structures beyond ~ 100 kpc) WDM and CDM yield identical results which agree with observations. For intermediate scales, WDM gives the correct abundance of substructures. Inside galaxy cores, below ~ 100 pc, N-body WDM classical physics simulations are incorrect because at such scales quantum WDM effects are important. WDM quantum calculations (Thomas-Fermi approach) provide galaxy cores, galaxy masses, velocity dispersions and density profiles in agreement with the observations. For a dark matter particle decoupling at thermal equilibrium (thermal relic), all evidences point out to a 2 keV particle. Remarkably enough, sterile neutrinos decouple out of thermal equilibrium with a primordial power spectrum similar to a 2 keV thermal relic when the sterile neutrino mass is about 7 keV. Therefore, WDM can be formed by 7 keV sterile neutrinos. Excitingly enough, Bulbul et al. (2014) announced the detection of a cluster X-ray emission line that could correspond to the decay of a 7.1 keV sterile neutrino and to a neutrino decay mixing angle of \sin^2 2 \theta ~ 7 10^{-11} . This is a further argument in favour of sterile neutrino WDM. Baryons, represent 10 % of DM or less in galaxies and are expected to give a correction to pure WDM results. The detection of the DM particle depends upon the particle physics model. Sterile neutrinos with keV scale mass (the main WDM candidate) can be detected in beta decay for Tritium and Renium and in the electron capture in Holmiun. The sterile neutrino decay into X rays can be detected observing DM dominated galaxies and through the distortion of the black-body CMB spectrum. So far, not a single valid objection arose against WDM.

The Thomas-Fermi approach to galaxy structure determines selfconsistently the fermionic warm dark matter (WDM) gravitational potential given the distribution function f(E). This framework is appropriate for macroscopic quantum systems: neutron stars, white dwarfs and WDM galaxies. Compact dwarf galaxies follow from the quantum degenerate regime, while dilute and large galaxies from the classical Boltzmann regime. We find analytic scaling relations for the main galaxy magnitudes as halo radius r_h, mass M_h and phase space density. The observational data for a large variety of galaxies are all well reproduced by these theoretical scaling relations. For the compact dwarfs, our results show small deviations from the scaling due to quantum macroscopic effects. We contrast the theoretical curves for the circular velocities and density profiles with the observational ones. All these results are independent of any WDM particle physics model, they only follow from the gravity interaction of the WDM particles and their fermionic nature. The theory rotation and density curves reproduce very well for r < r_h the observations of 10 different and independent sets of data for galaxy masses from 5x10^9 Msun till 5x10^{11} Msun. Our normalized circular velocity curves turn to be universal functions of r/r_h for all galaxies and reproduce very well the observational curves for r < r_h. Conclusion: the Thomas-Fermi approach correctly describes the galaxy structures (Abridged).

The Thomas-Fermi approach to galaxy structure determines selfconsistently and nonlinearly the gravitational potential of the fermionic WDM particles given their quantum distribution function f(E). Galaxy magnitudes as the halo radius r_h, mass M_h, velocity dispersion and phase space density are obtained. We derive the general equation of state for galaxies (relation between the pressure and the density), and provide an analytic expression. This clearly exhibits two regimes: (i) Large diluted galaxies for M_h > 2.3 10^6 Msun corresponding to temperatures T_0 > 0.017 K, described by the classical self gravitating WDM Boltzman regime and (ii) Compact dwarf galaxies for 1.6 10^6 Msun > M_h>M_{h,min}=30000 (2keV/m)^{16/5} Msun, T_0<0.011 K described by the quantum fermionic WDM regime. The T_0=0 degenerate quantum limit predicts the most compact and smallest galaxy (minimal radius and mass M_{h,min}). All magnitudes in the diluted regime exhibit square root of M_h scaling laws and are universal functions of r/r_h when normalized to their values at the origin or at r_h. We find that universality in galaxies (for M_h > 10^6 Msun) reflects the WDM perfect gas behaviour. These theoretical results contrasted to robust and independent sets of galaxy data remarkably reproduce the observations. For the small galaxies, 10^6>M_h>M_{h,min} corresponding to effective temperatures T_0 < 0.017 K, the equation of state is galaxy dependent and the profiles are no more universal. These non-universal properties in small galaxies account to the quantum physics of the WDM fermions in the compact regime. Our results are independent of any WDM particle physics model, they only follow from the gravitational interaction of the WDM particles and their fermionic quantum nature.

Analytic formulas reproducing the warm dark matter (WDM) primordial spectra are obtained for WDM particles decoupling in and out of thermal equilibrium which provide the initial data for WDM non-linear structure formation. We compute and analyze the corresponding WDM overdensities and compare them to the CDM case. We consider the ratio of the WDM to CDM primordial spectrum and the WDM to CDM overdensities: they turn to be self-similar functions of k/k_{1/2} and R/R_{1/2} respectively, k_{1/2} and R_{1/2} being the wavenumber and length where the WDM spectrum and overdensity are 1/2 of the respective CDM magnitudes. Both k_{1/2} and R_{1/2} show scaling as powers of the WDM particle mass m while the self-similar functions are independent of m. The WDM primordial spectrum sharply decreases around k_{1/2} with respect to the CDM spectrum, while the WDM overdensity slowly decreases around R_{1/2}. The nonlinear regions where WDM structure formation takes place are shown and compared to those in CDM: the WDM non-linear structures start to form later than in CDM, and as a general trend, decreasing the DM particle mass delays the onset of the non-linear regime. The non-linear regime starts earlier for smaller objects than for larger ones; smaller objects can form earlier both in WDM and CDM. We compute and analyze the differential mass function dN/dM for WDM at redshift z in the Press-Schechter approach. The WDM suppression effect of small scale structure increases with the redshift z. Our results for dN/dM are useful to be contrasted with observations, in particular for 4 < z < 12. We perfom all these studies for the most popular WDM particle physics models. Contrasting them to observations should point out the precise value of the WDM particle mass in the keV scale, and help to single out the best WDM particle physics model (Abridged).

Warm Dark Matter(WDM), considerably clarifies and simplifies galaxies and galaxy formation in agreement with observations. WDM essentially works, naturally reproducing the astronomical observations over all scales, small (galactic) as well as large and cosmological scales.Evidence that CDM, CDM+baryons and proposed tailored cures do not work in galaxies is staggering, The Chalonge Meudon Workshop 2012 approached DM in a fourfold coherent way: astronomical observations (galaxy and cluster properties, haloes, rotation curves, density profiles, surface density), LambdaWDM N-body simulations, WDM theory (Boltzmann-Vlasov evolution, halo mass functions, halo models, improved perturbative approachs), quantum WDM fermions forming the observed cores, WDM particle and nuclear physics (sterile neutrinos) and its experimental search. N Amorisco, P Biermann, S Das, H J de Vega, A Kamada, E Ferri(MARE),I D Karanchetsev, W Liao, M Lovell, M Papastergis, N G Sanchez, P Valageas,C Watson, J Zavala,He Zhang present their Highlights.Inside galaxy cores, N-body classical physics simulations are incorrect for WDM because of important quantum effects at such scales.Quantum WDM calculations (Thomas-Fermi) provide galaxy cores, galaxy masses, velocity dispersions and density profiles in agreement with observations.Baryons (16% of DM) are expected to give a correction to pure WDM results. The summary and conclusions by H. J. de Vega and N. G. Sanchez stress that all evidences point to a DM particle mass around 2 keV. Peter Biermann in his live minutes concludes that a few keV sterile neutrino is the most serious DM candidate. MARE -and hopefully KATRIN- could provide a sterile neutrino signal.There is a formidable WDM work to perform ahead of us, these highlights point research directions worthwhile to pursue.Photos of the Workshop are included (Abridged).

Quantum mechanics is necessary to compute galaxy structures at kpc scales and below. This is so because near the galaxy center, at scales below 10 - 100 pc, warm dark matter (WDM) quantum effects are important: observations show that the interparticle distance is of the order of, or smaller than the de Broglie wavelength for WDM. This explains why all classical (non-quantum) WDM N-body simulations fail to explain galactic cores and their sizes. We describe fermionic WDM galaxies in an analytic semiclassical framework based on the Thomas-Fermi approach, we resolve it numerically and find the main physical galaxy magnitudes: mass, halo radius, phase-space density, velocity dispersion, fully consistent with observations, including compact dwarf galaxies. Namely, fermionic WDM treated quantum mechanically, as it must be, reproduces the observed galaxy DM cores and their sizes. [In addition, as is known, WDM simulations produce the right DM structures in agreement with observations for scales > kpc]. We show that compact dwarf galaxies are natural quantum macroscopic objects supported against gravity by the fermionic WDM quantum pressure (quantum degenerate fermions) with a minimal galaxy mass and minimal velocity dispersion. Interestingly enough, the minimal galaxy mass implies a minimal mass m_{min} for the WDM particle. The lightest known dwarf galaxy (Willman I) implies m > m_{min} = 1.91 keV. These results and the observed halo radius and mass of the compact galaxies provide further indication that the WDM particle mass m is approximately around 2 keV.

We derive the main physical galaxy properties: mass, halo radius, phase space density and velocity dispersion from a semiclassical gravitational approach in which fermionic WDM is treated quantum mechanically. They turn out to be compatible with observations. Pauli Principle implies for the fermionic DM phase-space density Q(r) = rho(r)/sigma^3(r) the quantum bound Q(r) < K m^4/hbar^3, where m is the DM particle mass, sigma(r) is the DM velocity dispersion and K is a pure number of order one which we estimate. N-body galaxy simulations produce a divergent Q(r) at r = 0 violating this quantum bound. Combining this bound with the behaviour of Q(r) from simulations, the virial and galaxy data on Q implies lower bounds on the halo radius and a minimal distance r_{min} at which classical dynamics for DM fermions breaks down. This quantum bound rules out the presence of galaxy cusps for fermionic WDM, in agreement with astronomical observations, which show that the DM halos are cored. We show that compact dwarf galaxies are quantum objects supported against gravity by the fermionic WDM quantum pressure. Quantum mechanical calculations become necessary to compute galaxy structures at kpc scales and below. Classical N-body simulations are not valid at scales below r_{min}. We apply the Thomas-Fermi semiclassical approach to fermionic WDM galaxies and find the physical magnitudes: mass, halo radius, phase-space density, velocity dispersion, consistent with observations especially for compact dwarf galaxies. Namely, fermionic WDM treated quantum mechanically, as it must be, reproduces the sizes of the observed cores. The lightest known galaxy Willman I implies a lower bound for the WDM particle mass m > 0.96 keV. These results and the observed galaxies with halo radius > 30 pc and halo mass > 4 10^5 M_sun provide further indication that the WDM particle mass m is in the range 1-2 keV.

Sterile neutrinos with mass in the range of one to a few keV are important as extensions of the Standard Model of particle physics and are serious dark matter (DM) candidates. This DM mass scale (warm DM) is in agreement with both cosmological and galactic observations. We study the role of a keV sterile neutrino through its mixing with a light active neutrino in Rhenium 187 and Tritium beta decays. We pinpoint the energy spectrum of the beta particle, 0 < T_e < (Q_{beta} - m_s), as the region where a sterile neutrino could be detected and where its mass m_s could be measured. This energy region is at least 1 keV away rom the region suitable to measure the mass of the light active neutrino, located near the endpoint Q_{beta} . The emission of a keV sterile neutrino in a beta decay could show up as a small kink in the spectrum of the emitted beta particle. With this in view, we perform a careful calculation of the Rhenium and Tritium beta spectra and estimate the size of this perturbation by means of the dimensionless ratio R of the sterile neutrino to the active neutrino contributions. We comment on the possibility of searching for sterile neutrino signatures in two experiments which are currently running at present, MARE and KATRIN, focused on the Rhenium 187 and Tritium beta decays respectively.

We combine observed properties of galaxies as the core density and radius with the theoretical linear evolution of density fluctuations computed from first principles since the end of inflation till today. The halo radius r_0 is computed in terms of cosmological parameters. The theoretical density profiles rho(r)/rho(0) have an universal shape as a function of r/r_0 which reproduces the observations. We show that the linear approximation to the Boltzmann-Vlasov equation is valid for very large galaxies and correctly provides universal quantities which are common to all galaxies, as the surface density and density profile. By matching the theoretically computed surface density to its observed value we obtain (i) the decreasing of the phase-space density during the MD era (ii) the mass of the dark matter particle which turns to be between 1 and 2 keV and the decoupling temperature T_d which turns to be above 100 GeV (iii) the core vs. cusp discrimination: keV dark matter particles produce cored density profiles while wimps (m \sim 100 GeV, T_d \sim 5 GeV) produce cusped profiles at scales about 0.003 pc. These results are independent of the particle model and vary very little with the statistics of the dark matter particle. Non-universal galaxy quantities (which need to include non-linear effects as mergers and baryons) are reproduced in the linear approximation up to a factor of order one for the halo radius r_0, galaxy mass M_{gal}, halo central density rho_{0} and velocity dispersion sqrt{{\bar {v^2}}_{halo}} in the limiting case of large galaxies (both r_0 and M_{gal} large). This shows the power of the linear approximation scheme: although it cannot capture the whole content of the structure formation, it correctly provides universal quantities which as well as the main non-universal galaxy properties.

The Chalonge 15th Paris Cosmology Colloquium 2011 was held on 20-22 July in the historic Paris Observatory's Perrault building, in the Chalonge School spirit combining real cosmological/astrophysical data and hard theory predictive approach connected to them in the Warm Dark Matter Standard Model of the Universe: News and reviews from Herschel, QUIET, Atacama Cosmology Telescope (ACT), South Pole Telescole (SPT), Planck, PIXIE, the JWST, UFFO, KATRIN and MARE experiments; astrophysics, particle and nuclear physics warm dark matter (DM) searches and galactic observations, related theory and simulations, with the aim of synthesis, progress and clarification. Philippe Andre, Peter Biermann, Pasquale Blasi, Daniel Boyanovsky, Carlo Burigana, Hector de Vega, Joanna Dunkley, Gerry Gilmore, Alexander Kashlinsky, Alan Kogut, Anthony Lasenby, John Mather, Norma Sanchez, Alexei Smirnov, Sylvaine Turck-Chieze present here their highlights of the Colloquium. Ayuki Kamada and Sinziana Paduroiu present here their poster highlights. LambdaWDM (Warm Dark Matter) is progressing impressively over LambdaCDM whose galactic scale crisis and decline are staggering. The International School Daniel Chalonge issued an statement of strong support to the James Webb Space Telescope (JSWT). The Daniel Chalonge Medal 2011 was awarded to John C. Mather, Science PI of the JWST. Summary and conclusions are presented by H. J. de Vega, M. C. Falvella and N. G. Sanchez. Overall, LambdaWDM and keV scale DM particles deserve dedicated astronomical and laboratory experimental searches, theoretical work and simulations. KATRIN experiment in the future could perhaps adapt its set-up to look to keV scale sterile neutrinos. It will be a a fantastic discovery to detect dark matter in a beta decay. Photos of the Colloquium are included. (Abridged)

We solve the cosmological evolution of warm dark matter (WDM) density fluctuations with the Volterra integral equations of paper I. In the absence of neutrinos, the anisotropic stress vanishes and the Volterra equations reduce to a single integral equation. We solve numerically this equation both for DM fermions decoupling at equilibrium and DM sterile neutrinos decoupling out of equilibrium. We give the exact analytic solution for the density fluctuations and gravitational potential at zero wavenumber. We compute the density contrast as a function of the scale factor a for a wide range of wavenumbers k. At fixed a, the density contrast grows with k for k < k_c while it decreases for k > k_c, (k_c ~ 1.6/Mpc). The density contrast depends on k and a mainly through the product k a exhibiting a self-similar behavior. Our numerical density contrast for small k gently approaches our analytic solution for k = 0. For fixed k < 1/(60 kpc), the density contrast generically grows with a while for k > 1/(60 kpc) it exhibits oscillations since the RD era which become stronger as k grows. We compute the transfer function of the density contrast for thermal fermions and for sterile neutrinos in: a) the Dodelson-Widrow (DW) model and b) in a model with sterile neutrinos produced by a scalar particle decay. The transfer function grows with k for small k and then decreases after reaching a maximum at k = k_c reflecting the time evolution of the density contrast. The integral kernels in the Volterra equations are nonlocal in time and their falloff determine the memory of the past evolution since decoupling. This falloff is faster when DM decouples at equilibrium than when it decouples out of equilibrium. Although neutrinos and photons can be neglected in the MD era, they contribute in the MD era through their memory from the RD era.

We study the complete cosmological evolution of dark matter (DM) density fluctuations for DM particles that decoupled being ultrarelativistic during the radiation dominated era which is the case of keV scale warm DM (WDM). The new framework presented here can be applied to other types of DM and in particular we extend it to cold DM (CDM). The collisionless and linearized Boltzmann-Vlasov equations (B-V) for WDM and neutrinos in the presence of photons and coupled to the linearized Einstein equations are studied in detail in the presence of anisotropic stress with the Newtonian potential generically different from the spatial curvature perturbations. We recast this full system of B-V equations for DM and neutrinos into a system of coupled Volterra integral equations. These Volterra-type equations are valid both in the radiation dominated (RD) and matter dominated (MD) eras during which the WDM particles are ultrarelativistic and then nonrelativistic. This generalizes the so-called Gilbert integral equation only valid for nonrelativistic particles in the MD era. We succeed to reduce the system of four Volterra integral equations for the density and anisotropic stress fluctuations of DM and neutrinos into a system of only two coupled Volterra equations. The kernels and inhomogeneities in these equations are explicitly given functions. Combining the Boltzmann-Vlasov equations and the linearized Einstein equations constrain the initial conditions on the distribution functions and gravitational potentials. In the absence of neutrinos the anisotropic stress vanishes and the Volterra-type equations reduce to a single integral equation. These Volterra integral equations provide a useful and precise framework to compute the primordial WDM fluctuations over a wide range of scales including small scales up to k ~ 1/5 kpc.

Warm Dark Matter (WDM) research is progressing fast, the subject is new and WDM essentially works, naturally reproducing the astronomical observations over all scales: small (galactic) and large (cosmological) scales (LambdaWDM). Evidence that Cold Dark Matter (LambdaCDM) and its proposed tailored cures do not work at small scales is staggering. Fedor Bezrukov, Pier-Stefano Corasaniti, Hector J. de Vega, Stefano Ettori, Frederic Hessmann, Ayuki Kamada, Marco Lombardi, Alexander Merle, Christian Moni Bidin, Angelo Nucciotti on behalf of the MARE collaboration, Sinziana Paduroiu, Henri Plana, Norma Sanchez, Patrick Valageas, Shun Zhou present here their highlights of the Workshop. LambdaWDM simulations with keV particles remarkably reproduce the observations, small and large structures and velocity functions. Cored DM halos and WDM are clearly determined from theory and astronomical observations, they naturally produce the observed structures at all scales. keV sterile neutrinos are the leading candidates, they naturally appear extensions of the standard model of particle physics. Astrophysical constraints including Lyman alpha bounds put its mass in the range 1< m <13 keV. Predictions for EUCLID and PLANCK have been presented. MARE and an adapted KATRIN experiment could detect a keV sterile neutrino. It will be a fantastic discovery to detect dark matter in a beta decay. A formidable WDM work to perform is ahead of us, these highlights point out some relevant research directions to put the effort. Photos of the Workshop are included (Abridged).

The MCMC analysis of the CMB+LSS data in the context of the Ginsburg-Landau approach to inflation indicated that the fourth degree double--well inflaton potential best fits the present CMB and LSS data. This provided a lower bound for the ratio r of the tensor to scalar fluctuations and as most probable value r = 0.05, within reach of the forthcoming CMB observations. We systematically analyze here the effects of arbitrary higher order terms in the inflaton potential on the CMB observables: spectral index ns and ratio r. Furthermore, we compute in close form the inflaton potential dynamically generated when the inflaton field is a fermion condensate in the inflationary universe. This inflaton potential turns to belong to the Ginsburg-Landau class too. The theoretical values in the (ns,r) plane for all double well inflaton potentials in the Ginsburg-Landau approach (including the potential generated by fermions) turn to be inside a universal banana-shaped region B. The upper border of the banana-shaped region B is given by the fourth order double--well potential and provides an upper bound for the ratio r.The lower border of B is defined by the quadratic plus an infinite barrier inflaton potential and provides a lower bound for the ratio r. For example, the current best value of the spectral index ns = 0.964, implies r is in the interval: 0.021 < r < 0.053. Interestingly enough, this range is within reach of forthcoming CMB observations.

Cumulative observational evidence confirm that the surface gravity of dark matter (DM) halo mu_{0 D} = r_0 rho_0 where r_0 and rho_0 are the halo core radius and central density, respectively, is nearly constant and independent of galaxy luminosity for a high number of galactic systems (spirals, dwarf irregular and spheroidals, elliptics) spanning over 14 magnitudes in luminosity and of different Hubble types. Remarkably, its numerical value mu_{0 D} = 140 M_{sun}/pc^2 = (18.6 Mev)^3 is approximately the same (up to a factor of two) in all these systems. First, we present the physical consequences of the independence of mu_{0 D} on r_0: the energy scales as the volume sim r_0^3 while the mass and the entropy scale as the surface ~ r_0^2 and the surface times log r_0, respectively. Namely, the entropy scales similarly to the black-hole entropy but with a much smaller coefficient. Second, we compute the surface gravity and the density profile for small scales from first principles and the evolution of primordial density fluctuations since the end of inflation till today using the linearized Boltzmann-Vlasov equation. The density profile rho_{lin}(r) obtained in this way decreases as r^{-1-n_s/2} for intermediate scales where n_s = 0.964 is the primordial spectral index. This scaling is in remarkable agreement with the empirical behaviour found observationally and in N-body simulations: r^{-1.6\pm 0.4}. The observed value of mu_{0 D} indicates that the DM particle mass m is in the keV scale. The theoretically derived density profiles rho_{lin}(r) turn to be cored for m in the keV scale and they look as cusped for m in the GeV scale or beyond. We consider both fermions and bosons as DM particles decoupling either ultrarelativistically or non-relativistically. Our results do not use any particle physics model and vary slightly with the statistics of the DM particle.

The Chalonge 14th Paris Cosmology Colloquium was held on 22-24 July 2010 in Paris Observatory on the Standard Model of the Universe: News from WMAP7, BICEP, QUAD, SPT, AMI, ACT, Planck, QUIJOTE and Herschel; dark matter (DM) searches and galactic observations; related theory and simulations. %aiming synthesis, progress and clarification. P Biermann, D Boyanovsky, A Cooray, C Destri, H de Vega, G Gilmore, S Gottlober, E Komatsu, S McGaugh, A Lasenby, R Rebolo, P Salucci, N Sanchez and A Tikhonov present here their highlights of the Colloquium. Inflection points emerged: LambdaWDM (Warm DM) emerges impressively over LambdaCDM whose galactic scale problems are ever-increasing. Summary and conclusions by H. J. de Vega, M. C. Falvella and N. G. Sanchez stress among other points: (i) Primordial CMB gaussianity is confirmed. Inflation effective theory predicts a tensor to scalar ratio 0.05-0.04 at reach/border line of next CMB observations, early fast-roll inflation provides lowest multipoles depression. SZ amplitudes are smaller than expected: CMB and X-ray data agree but intracluster models need revision and relaxed/non-relaxed clusters distinction. (ii) cosmic ray positron excess is explained naturally by astrophysical processes, annihilating/decaying dark matter needs growing tailoring. (iii) Cored (non cusped) DM halos and warm (keV scale mass) DM are increasingly favored from theory and observations, naturally producing observed small scale structures, wimps turn strongly disfavoured. LambdaWDM 1 keV simulations well reproduce observations. Evidence that LambdaCDM does not work at small scales is staggering. P Biermann presents his live minutes of the Colloquium and concludes that a keV sterile neutrino is the most interesting DM candidate. Photos of the Colloquium are included.

The Planck satellite is right now measuring with unprecedented accuracy the primary Background CMB anisotropies. The Standard Model of the Universe (including inflation) provides the context to analyze the CMB and other data. The Planck performance for r, the tensor to scalar ratio related to primordial B mode polarization, will depend on the quality of the data analysis. The Ginsburg Landau approach to inflation allows to take high benefit of the CMB data. The fourth degree double well inflaton potential gives an excellent fit to the current CMB+LSS data. We evaluate the Planck precision to the recovery of cosmological parameters within a reasonable toy model for residuals of systematic effects of instrumental and astrophysical origin based on publicly available information.We use and test two relevant models: the LambdaCDMr model, i.e. the standard LambdaCDM model augmented by r, and the LambdaCDMrT model, where the scalar spectral index, n_s, and r are related through the theoretical `banana-shaped' curve r = r(n_s) coming from the double-well inflaton potential. In the latter case, r = r(n_s) is imposed as a hard constraint in the MCMC data analysis. We take into account the white noise sensitivity of Planck in the 70, 100 and 143 GHz channels as well as the residuals from systematics errors and foregrounds. Foreground residuals turn to affect only the cosmological parameters sensitive to the B modes. The best value for r in the presence of residuals turns to be about r simeq 0.04 for both the LambdaCDMr and the LambdaCDMrT models. We compute the B mode detection probability by the most sensitive HFI-143 channel. At the level of foreground residual equal to 30% of our toy model only a 68% CL detection of r is very likely. For a 95% CL detection the level of foreground residual should be reduced to 10% or lower of the adopted toy model (ABRIDGED).

The 13th Chalonge Paris Cosmology Colloquium was held at the historic Perrault building of Observatoire de Paris on 23-25 July 2009. The Colloquium was held in the spirit of the Chalonge School, putting together real cosmological and astrophysical data and hard theory predictive approach in the framework of the Standard Model of the Universe. Peter Biermann, James Bullock, Claudio Destri, Hector de Vega, Massimo Giovannini, Sasha Kashlinsky, Eiichiro Komatsu, Anthony Lasenby and Norma G. Sanchez present here their highlights of the Colloquium. The summary and conclusions by H. J. de Vega, M. C. Falvella and N. G. Sanchez stress among other points: (i) the primordial CMB fluctuations are almost gaussian, large primordial non-gaussianity and large running index are strongly disfavored. The primordial graviton ratio r is predicted in the effective theory of inflation to be between 0.021 and 0.053, at reach of the next CMB observations. (ii) Dark energy observations are consistent with the cosmological constant. (iii) The cosmic ray positron excess recently observed is explained by natural astrophysical processes. (iv) The heavy dark matter particle candidates ( > 1 GeV, wimps) became strongly disfavored, while cored (non cusped) dark matter halos and light (keV scale mass) dark matter are being increasingly favored from theory and astrophysical observations. The Daniel Chalonge Medal 2009 was awarded to Peter Biermann during the Colloquium. Photos of the Colloquium are included.

The CIAS Chalonge Workshop `Dark Matter in the Universe and Universal Properties of Galaxies: Theory and Observations', was held at the Meudon Ch^ateau of Observatoire de Paris on 8-11 June 2010. The Workshop approached DM in a fourfold way: astronomical observations of DM structures (galaxy properties, halos, rotation curves and density profiles), DM numerical simulations (with and without baryons), theoretical astrophysics and cosmology (kinetic theory, Boltzmann-Vlasov evolution), astroparticle physics. Peter Biermann, Alfonso Cavaliere, Hector J. de Vega, Gianfranco Gentile, Chandra Jog, Andrea Lapi, Paolo Salucci, Norma G. Sanchez, Pasquale Serpico, Rainer Stiele, Janine van Eymeren and Markus Weber present here their highlights of the Workshop. The summary and conclusions by H. J. de Vega and N. G. Sanchez stress among other points the growing evidence that DM particles have a mass in the keV scale and that those keV scale particles naturally produce the small scale structures observed in galaxies. Wimps (DM particles heavier than 1 GeV) are strongly disfavoured combining theory with galaxy astronomical observations. Peter Biermann presents his live minutes of the Workshop and concludes that a right-handed sterile neutrino of mass of a few keV is the most interesting DM candidate. Photos of the Workshop are included.

We present here the effective theory of inflation `a la Ginsburg-Landau in which the inflaton potential is a polynomial. The slow-roll expansion becomes a systematic 1/N expansion where N ~ 60. The spectral index and the ratio of tensor/scalar fluctuations are n_s - 1 = O(1/N), r = O(1/N) while the running turns to be d n_s/d \ln k = O(1/N^2) and can be neglected. The energy scale of inflation M ~ 0.7 10^{16} GeV is completely determined by the amplitude of the scalar adiabatic fluctuations. A complete analytic study plus the Monte Carlo Markov Chains (MCMC) analysis of the available CMB+LSS data showed: (a) the spontaneous breaking of the phi -> - phi symmetry of the inflaton potential. (b) a lower bound for r: r > 0.023 (95% CL) and r > 0.046 (68% CL). (c) The preferred inflation potential is a double well, even function of the field with a moderate quartic coupling yielding as most probable values: n_s = 0.964, r = 0.051. This value for r is within reach of forthcoming CMB observations. We investigate the DM properties using cosmological theory and the galaxy observations. Our DM analysis is independent of the particle physics model for DM and it is based on the DM phase-space density rho_{DM}/sigma^3_{DM}. We derive explicit formulas for the DM particle mass m and for the number of ultrarelativistic degrees of freedom g_d (hence the temperature) at decoupling. We find that m turns to be at the keV scale. The keV scale DM is non-relativistic during structure formation, reproduces the small and large scale structure but it cannot be responsible of the e^+ and pbar excess in cosmic rays which can be explained by astrophysical mechanisms (Abridged).

We study the entire coupled evolution of the inflaton and the scale factor for general initial conditions at a given initial time. The generic early universe evolution has three stages: decelerated fast-roll followed by inflationary fast roll and then inflationary slow-roll. This evolution is valid for all regular inflaton potentials. In addition, we find a special (extreme) slow-roll solution starting at t = -infty in which the fast-roll stages are absent. At some time t = t_*, the generic evolution backwards in time reaches a mathematical singu- larity where a(t) vanishes and Hubble becomes singular. We find the general behaviour near the singularity. The classical inflaton description is valid for t-t_* > 10 t_{Planck} well before the beginning of inflation, quantum loop effects are negligible there. The singularity is never reached in the validity region of the classical treatment and therefore it is not a real physical phenomenon here. The whole evolution of the fluctuations is computed. The Bunch-Davies initial conditions (BDic) are generalized for the present case. The power spectrum gets dynamically modified by the effect of the fast-roll eras and the BDic choice at a finite time through the transfer function D(k) of initial conditions. D(0) = 0. D(k) presents a first peak for k ~ 2/eta_0 (eta_0 being the conformal initial time), then oscillates with decreasing amplitude and vanishes asymptotically for k -> infty. The transfer function D(k) affects the low CMB multipoles C_l: the change Delta C_l/C_l for l=1-5 is computed as a function of the starting instant of the fluctuations t_0. CMB quadrupole observations give large suppressions which are well reproduced here(Abridged)

We present a model independent analysis of dark matter (DM) both decoupling ultra relativistic (UR) and non-relativistic (NR) based in the phase-space density D = rho_{DM}/sigma^3_{DM}. We derive explicit formulas for the DM particle mass m and for the number of ultra relativistic degrees of freedom g_d at decoupling. We find that for DM particles decoupling UR both at local thermal equilibrium (LTE) and out of LTE, m turns to be at the keV scale. For example, for DM Majorana fermions decoupling at LTE the mass results m ~ 0.85 keV. For DM particles decoupling NR, sqrt{m T_d} results in the keV scale (T_d is the decoupling temperature) and the m value is consistent with the keV scale. In all cases, DM turns to be cold DM (CDM). Also, lower and upper bounds on the DM annihilation cross-section for NR decoupling are derived. We evaluate the free-streaming (Jeans') length and Jeans' mass: they result independent of the type of DM except for the DM self-gravity dynamics. The free-streaming length today results in the kpc range. These results are based on our theoretical analysis, astronomical observations of dwarf spheroidal satellite galaxies in the Milky Way and N-body numerical simulations. We analyze and discuss the results on D from analytic approximate formulas both for linear fluctuations and the (non-linear) spherical model and from N-body simulations results. We obtain in this way upper bounds for the DM particle mass which all result below the 100 keV range.

Inflation is part of the Standard Model of the Universe supported by CMB and large scale structure LSS datasets. This review presents new developments of inflation in three main chapters. (I): The effective theory of inflation a la Ginsburg-Landau (GL): the inflaton potential is a polynomial with universal form making explicit the inflation energy scale M, the Planck mass and the inflation e-folds number N ~ 60. The slow-roll expansion becomes a systematic 1/N expansion and the inflaton couplings are naturally small as powers of (M/M_{Pl})^2. The spectral index (n_s - 1) and the ratio of tensor/scalar fluctuations r are O(1/N), the running index is O(1/N^2). M ~ 0.7 10^{16} GeV is completely determined by the scalar adiabatic fluctuations amplitude. (II): A Monte Carlo Markov Chains (MCMC) analysis of the CMB+LSS data (including WMAP5) with our analytic theoretical results yields: a lower bound for r (new inflation): r > 0.023 (95%CL), r > 0.046 (68%CL); the preferred inflation potential is a double well, even function of the field yielding as most probable values n_s ~ 0.964, r ~ 0.051. This value for r is within reach of forthcoming CMB observations. Slow-roll inflation is generically preceded by a short fast-roll stage which leads to a suppression of the CMB quadrupoles. MCMC analysis of the WMAP+SDSS data shows that fast-roll fits the TT, TE and EE modes well reproducing the quadrupole suppression and fixes the total number of efolds of inflation to be N_{total} ~ 64. (III) Quantum loop corrections are very small and controlled by powers of (H /M_{Pl})^2 ~ 10^{-9} which validates the effective theory of inflation. We show how powerful is the GL theory of inflation in predicting observables.

The transfer function $T(k)$ of dark matter (DM) perturbations during matter domination is obtained by solving the collisionless Boltzmann-Vlasov equation. We find an \emph{exact} expression for $T(k)$ for \emph{arbitrary} distribution functions of decoupled particles and initial conditions}. We find a remarkably accurate and simple approximation valid on all scales of cosmological relevance for structure formation in the linear regime. The natural scale of suppression is the free streaming wavevector at matter-radiation equality, $ k_{fs}(t_{eq}) = [{4\pi\rho_{0M}}/{[< \vec{V}^2> (1+z_{eq})]} ]^\frac12 $. An important ingredient is a non-local kernel determined by the distribution functions of the decoupled particles which describes the \emph{memory of the initial conditions and gravitational clustering} and yields a correction to the fluid description. Distribution functions that favor the small momentum region lead to an \emph{enhancement of power at small scales} $ k > k_{fs}(t_{eq}) $. For DM thermal relics that decoupled while ultrarelativistic we find $ k_{fs}(t_{eq}) \simeq 0.003 (g_d/2)^\frac13 (m/\mathrm{keV}) [\mathrm{kpc}]^{-1} $, where $ g_d $ is the number of degrees of freedom at decoupling. For WIMPS we obtain $ k_{fs}(t_{eq}) = 5.88 (g_d/2)^\frac13 (m/100 \mathrm{GeV})^\frac12 (T_d/10 \mathrm{MeV})^\frac12 [\mathrm{pc}]^{-1} $. For $k\ll k_{fs}(t_{eq})$, $T(k) \sim 1-\mathrm{C}[k/k_{fs}(t_{eq})]^2 $ where $C =\mathrm{O}(1)$ for all cases considered and simple and accurate fits for \emph{small} scales.

Quantum fast-roll initial conditions for the inflaton which are different from the classical fast-roll conditions and from the quantum slow-roll conditions can lead to inflation that last long enough. These quantum fast-roll initial conditions for the inflaton allow for kinetic energies of the same order of the potential energies and nonperturbative inflaton modes with nonzero wavenumbers. Their evolution starts with a transitory epoch where the redshift due to the expansion succeeds to assemble the quantum excited modes of the inflaton in a homogeneous (zero mode) condensate, and the large value of the Hubble parameter succeeds to overdamp the fast-roll of the redshifted inflaton modes. After this transitory stage the effective classical slow-roll epoch is reached. Most of the efolds are produced during the slow-roll epoch and we recover the classical slow-roll results for the scalar and tensor metric perturbations plus corrections. These corrections are important, both for scalar and for tensor perturbations, if scales which are horizon-size today exited the horizon by the end of the transitory stage and as a consequence the lower CMB multipoles get suppressed (fast-roll) or enhanced (precondensate). These two types of corrections can compete and combine in a scale dependent manner. They arise as natural consequences of the quantum nonperturbative inflaton dynamics, and provide a consistent and contrastable model for the origin of the suppression of the quadrupole and for other departures of the low CMB multipoles from the slow-roll inflation-LambdaCMB model which are to be contrasted to the TE and EE multipoles and to the forthcoming and future CMB data.

Generically,the classical evolution of the inflaton has a brief fast roll stage that precedes the slow roll regime. The fast roll stage leads to a purely attractive potential for the curvature and tensor perturbations (this potential is purely repulsive in slowroll). This attractive potential depresses the CMB quadrupole moment for the curvature and B-mode spectra. A single new cosmological parameter emerges: the comoving wavenumber k_1 characteristic scale of this attractive potential.k_1 happens to exit the horizon precisely at the transition from fast to slowroll.The fastroll stage dynamically modifies the initial power spectrum by a transfer function D(k). We compute D(k) by solving the inflaton evolution equations. D(k) suppresses the primordial power for k < k_1 and enjoys the scaling property D(k) = Psi(k/k_1) where Psi(x) is an universal function. We perform a MCMC analysis of the WMAP/SDSS data including the fast-roll stage and find k_1 = 0.266/Gpc. The quadrupole mode k_Q = 0.242/Gpc exits the horizon 1/10 of an efold before k_1. We compare the fastroll fit with a fit with a sharp lower cutoff on the primordial power. Fastroll provides a slightly better fit than a sharp cutoff for the TT,TE and EE modes.Moreover, our fits provide non-zero lower bounds for r, while for the other cosmological parameters we essentially get those of the LCDM model. The fact that k_Q exits the horizon before the slowroll stage implies an upper bound in the total number of efolds N_{tot} during inflation.Combining this with estimates during the radiation dominated era we obtain N_{tot} ~ 66, with the bounds 62 < N_{tot} < 82.We repeated the analysis with the WMAP-5/ACBAR-2007/SDSS data confirming the overall picture.

Mass bounds on dark matter (DM) candidates are obtained for particles decoupling in or out of equilibrium with {\bf arbitrary} isotropic and homogeneous distribution functions. A coarse grained Liouville invariant primordial phase space density $ \mathcal D $ is introduced. Combining its value with recent photometric and kinematic data on dwarf spheroidal satellite galaxies in the Milky Way (dShps), the DM density today and $N$-body simulations, yields upper and lower bounds on the mass, primordial phase space densities and velocity dispersion of the DM candidates. The mass of the DM particles is bound in the few keV range. If chemical freeze out occurs before thermal decoupling, light bosonic particles can Bose-condense. Such Bose-Einstein {\it condensate} is studied as a dark matter candidate. Depending on the relation between the critical($T_c$)and decoupling($T_d$)temperatures, a BEC light relic could act as CDM but the decoupling scale must be {\it higher} than the electroweak scale. The condensate tightens the upper bound on the particle's mass. Non-equilibrium scenarios that describe particle production and partial thermalization, sterile neutrinos produced out of equilibrium and other DM models are analyzed in detail obtaining bounds on their mass, primordial phase space density and velocity dispersion. Light thermal relics with $ m \sim \mathrm{few} \mathrm{keV}$ and sterile neutrinos lead to a primordial phase space density compatible with {\bf cored} dShps and disfavor cusped satellites. Light Bose condensed DM candidates yield phase space densities consistent with {\bf cores} and if $ T_c\gg T_d $ also with cusps. Phase space density bounds from N-body simulations suggest a potential tension for WIMPS with $ m \sim 100 \mathrm{GeV},T_d \sim 10 \mathrm{MeV} $.

We perform a MCMC (Monte Carlo Markov Chains) analysis of the available CMB and LSS data (including the three years WMAP data) with single field slow-roll new inflation and chaotic inflation models. We do this within our approach to inflation as an effective field theory in the GinsburgLandau spirit with fourth degree trinomial potentials in the inflaton field phi.We derive explicit formu- lae and study in detail the spectral index ns of the adiabatic fluctuations the ratio r of tensor to scalar fluctuations and the running index dn_s/dln k. We use these analytic formulas as hard constraints on n_s and r in the MCMC analysis.Our analysis differs in this crucial aspect from previous MCMC studies in the literature involving the WMAP3 data. Our results are as follow: (i) The data strongly indicate the breaking (whether spontaneous or explicit) of the phi -> -phi symmetry of the inflaton potentials both for new and for chaotic inflation.(ii)Trinomial new inflation naturally satisfies this requirement and provides an excellent fit to the data.(iii)Trinomial chaotic inflation produces the best fit in a very narrow corner of the parameter space.(iv) The chaotic symmetric trinomial potential is almost certainly ruled out(at 95% CL).In tri- nomial chaotic inflation the MCMC runs go towards a potential in the boundary of the parameter space and which ressembles a spontaneously symmetry broken potential of new inflation. (v) The above results and further physical analysis here lead us to conclude that new inflation gives the best description of the data.(vi) We find a lower bound for r within trinomial new inflation potentials r > 0.016 (95% CL) and r > 0.049 (68% CL). (vii) The preferred new inflation trinomial potential is a double well, even function of the field yielding as most probable values: n_s ~ 0.958, r ~ 0.055. Such r could be observed soon.

The CMB power spectra are studied for different families of single field new and chaotic inflation models in the effective field theory approach to inflation. We implement a systematic expansion in 1/N_e where N_e ~ 50 is the number of e-folds before the end of inflation. We study the dependence of the observables (n_s, r and dn_s/d\ln k) on the degree of the potential (2 n) and confront them to the WMAP3 and large scale structure data: this shows in general that fourth degree potentials (n=2) provide the best fit to the data; the window of consistency with the WMAP3 and LSS data narrows for growing n. New inflation models yield a good fit to the r and n_s data in a wide range of field and parameter space. Small field inflation yields r<0.16 while large field inflation yields r>0.16 (for N_e=50). All members of the new inflation family predict a small but negative running -4(n+1)10^{-4}< dn_s/d\ln k <-2 10^{-4}. A reconstruction program is carried out suggesting quite generallythat for n_s consistent with the WMAP3 and LSS data and r<0.1 the symmetry breaking scale for new inflation is |\phi_0| ~ 10 M_{Pl} while the field scale at Hubble crossing is |\phi_{50}| ~ M_{Pl}.The family of chaotic models feature r>0.16 (for N_e=50) and only a restricted subset of chaotic models are consistent with the combined WMAP3 bounds on r, n_s,dn_s/d\ln k with a narrow window in field amplitude around |\phi_{50}| ~ 15 M_{Pl}.A measurement of r<0.16 (for N_e =50) will distinctly rule out a large class of chaotic scenarios and will favor new inflationary models. As a general consequence, new inflation emerges more favoured than chaotic inflation.

We uncover the general mechanism producing the dark energy(DE).This is only based on well known quantum physics and cosmology. We show that the observed DE originates from the cosmological quantum vacuum of light particles which provides a continuous energy distribution able to reproduce the data. Bosons give positive contributions to the DE while fermions yield negative contri- butions. As usual in field theory, ultraviolet divergences are subtracted from the physical quantities. The subtractions respect the symmetries of the theory and we normalize the physical quantities to be zero for the Minkowski vacuum. The resulting finite contributions to the energy density and the pressure from the quantum vacuum grow as log a(t) where a(t) is the scale factor, while the particle contributions dilute as 1/a^3(t), as it must be for massive particles. The DE equation of state P=w(z)H turns to be w(z)<-1 with w(z) asymptotically reaching the value -1 from below.A scalar particle can produce the observed DE through its quantum cosmological vacuum provided:(i)its mass is of the order of 10^{-3}eV=1 meV,(ii) it is very weakly coupled and (iii) it is stable on the time scale of the age of the universe.The axion vacuum thus appears as a natural candidate. The neutrino vacuum (especially the lightest mass eigenstate) can give negative contributions to the DE. We find that w(z=0) is slightly below -1 by an amount ranging from -1.5 10^{-3}to -8 10^{-3} while the axion mass results between 4 and 5 meV.We find that the universe will expand in the future faster than the de Sitter universe, as an exponential in the square of the cosmic time. DE arises from the quantum vacua of light particles in FRW cosmological space time in an analogous way to the Casimir effect in Minkowski spacetime with non trivial boundaries.

We investigate the issue of initial conditions of curvature and tensor perturba- tions at the beginning of slow roll inflation and their effect on the power spectra. Renormalizability and small back reaction constrain the high k behavior of the Bogoliubov coefficients that define these initial conditions.We introduce a transfer function D(k) which encodes the effect of generic initial conditions on the power spectra. The constraint from renormalizability and small back reaction entails that D(k) < mu^2/k^2 for large k, implying that observable effects from initial conditions are more prominent in the low multipoles. This behavior affects the CMB quadrupole by the observed amount \~10-20% when mu is of the order of the energy scale of inflation. The effects on high l-multipoles are suppressed by a factor ~1/l^2 due to the fall off of D(k) for large wavevectors k. We show that the determination of generic initial conditions for the fluc- tuations is equivalent to the scattering problem by a potential V(eta) localized just prior to the slow roll stage. Such potential leads to a transfer function D(k) which automatically obeys the renormalizability and small backreaction constraints. We find that an attractive potential V(eta) yields a suppression of the lower CMB multipoles.Both for curvature and tensor modes, the quadrupole suppression depends only on the energy scale of V(eta) and on the time interval where V(eta) is nonzero. A suppression of the quadrupole for curvature pertur- bations consistent with the data is obtained when the scale of the potential is of the order of k^2_Q where k_Q is the wavevector whose physical wavelength is the Hubble radius today.

Within the effective field theory of inflation, an initialization of the classical dynamics of the inflaton with approximate equipartition between the kinetic and potential energy of the inflaton leads to a brief fast roll stage that precedes the slow roll regime. The fast roll stage leads to an attractive potential in the wave equations for the mode functions of curvature and tensor perturbations. The evolution of the inflationary perturbations is equivalent to the scattering by this potential and a useful dictionary between the scattering data and observables is established.Implementing methods from scattering theory we prove that this attractive potential leads to a suppression of the quadrupole moment for CMB and B-mode angular power spectra. The scale of the potential is determined by the Hubble parameter during slow roll. Within the effective field theory of inflation at the grand unification (GUT) energy scale we find that if inflation lasts a total number of efolds N_{tot} ~ 59, there is a 10-20% suppression of the CMB quadrupole and about 2-4% suppression of the tensor quadrupole. The suppression of higher multipoles is smaller, falling off as 1/l^2. The suppression is much smaller for N_{tot} > 59, therefore if the observable suppression originates in the fast roll stage, there is the upper bound N_{tot} ~ 59.

The self-gravitating gas in thermal equilibrium is studied using a Newtonian potential regularized at short distances. This short distance cutoff permits us to obtain a complete description of the gas including its collapsed phase. We give a field theory description of the N-body regularized self-gravitating gas in the canonical ensemble.The corresponding functional integral is dominated in the N -> infty limit by saddle points which provide a mean field description. The well-known dilute solutions (isothermal spheres) are recovered. We find new solutions which are regular in the regularized theory but become singular in the zero cutoff limit. They describe collapsed configurations where the particles are densely concentrated in a region of the size of the cutoff. These collapsed solutions provide the absolute minimum of the free energy. We find further new solutions which interpolate between the collapsed and the dilute configurations and describe tunneling processes where the gas collapses. The transition proba- bility for such collapse processes turns out to be extremely small for large N. That is, the dilute solutions are in practice stable in the regime where they are locally stable.

We study the single field slow-roll inflation models that better fit the available CMB and LSS data including the three years WMAP data: new inflation and hybrid inflation. We study them as effective field theories in the Ginsburg-Landau context: a trinomial potential turns out to be a simple and well motivated model. The compute the spectral index n_s of the adiabatic fluctuations, the ratio r of tensor to scalar fluctuations and the running index d n_s/dln k, derive explicit formulae and provide relevant plots. In new inflation, and for the three years WMAP central value n_s = 0.95, we predict 0.03<r<0.04 and -0.00070<d n_s/d ln k<-0.00055. In hybrid inflation, and for n_s = 0.95, we predict r = 0.2 and dn_s/dln k=-0.001.We find in new inflation n_s bounded from above by n_{smax} = 0.9615 and that r is a two valued function of n_s in the interval 0.96<n_s<0.9615. In the first branch we find r<r_{max} = 0.1148.In hybrid inflation we find a critical value mu_{0 crit}^2 for the mass parameter mu_0^2 of the field sigma coupled to the inflaton.For mu_0^2< Lambda_0 M_{Pl}^2/192, hybrid inflation is ruled out by the WMAP three years data since it yields n_s>1. Hybrid inflation for mu_0^2>Lambda_0 M_{Pl}^2/192 can fullfill all the present CMB+LSS data. Even if chaotic inflation predicts n_s values compatible with the data, chaotic inflation is disfavoured since it predicts a too high value for the ratio r=0.27. The model which best fits the current data and which best prepares the way to the expected data r < 0.1, is the trinomial potential with negative mass term: new inflation.

The evolution of the Universe is the ultimate laboratory to study fundamental physics across energy scales that span about 25 orders of magnitude: from the grand unification scale through particle and nuclear physics scales down to the scale of atomic physics. The standard models of cosmology and particle physics provide the basic understanding of the early and present Universe and predict a series of phase transitions that occurred in succession during the expansion and cooling history of the Universe. We survey these phase transitions, highlighting the equilibrium and non-equilibrium effects as well as their observational and cosmological consequences. We discuss the current theoretical and experimental programs to study phase transitions in QCD and nuclear matter in accelerators along with the new results on novel states of matter as well as on multi- fragmentation in nuclear matter. A critical assessment of similarities and differences between the conditions in the early universe and those in ultra- relativistic heavy ion collisions is presented. Cosmological observations and accelerator experiments are converging towards an unprecedented understanding of the early and present Universe.

We show that the self-gravitating gas at thermal equilibrium has an infinite volume limit in the three ensembles (GCE, CE, MCE) when (N, V) -> infty, keeping N/V^{1/3} fixed, that is, with eta = G m^2 N/[ V^{1/3} T] fixed. We develop MonteCarlo simulations, analytic mean field methods (MF) and low density expansions. We compute the equation of state and find it to be locally p(r) = T rho_V(r), that is a local ideal gas equation of state. The system is in a gaseous phase for eta < eta_T = 1.51024...and collapses into a very dense object for eta > eta_T in the CE with the pressure becoming large and negative. The isothermal compressibility diverges at eta = eta_T. We compute the fluctuations around mean field for the three ensembles. We show that the particle distribution can be described by a Haussdorf dimension 1 < D < 3.

Slow-roll inflation is studied as an effective field theory.We find as consistent form of the inflaton potential V(phi)=N M^4 w(phi/[sqrt{N}M_P]) where phi is the inflaton field, M the inflation energy scale, M_P the Planck mass, and N~50 the number of efolds since the relevant modes exited the horizon till the end of inflation. The dimensionless function w(chi) and field chi are O(1). The WMAP value for the amplitude of scalar adiabatic fluctuations |\Delta_{k ad}^(S)| fixes the inflation scale M ~ 0.77 10^16 GeV precisely at the GUT scale. This general form of the potential makes manifest that the slow roll expansion is an expansion in 1/N. Powers of 1/N count the orders in the slow roll expansion.This form of the inflaton potential suggests that the super symmetry breaking scale is at the inflation and GUT scales.A Ginzburg-Landau realization of this inflaton potential reveals that Hubble, inflaton mass and non-linear couplings are of the see-saw form in terms of the small ratio M/M_P. For example, the quartic coupling lambda ~ 1/N (M/M_P)^4.The smallness of the non-linear couplings is not a result of fine tuning but a natural consequence of the validity of the effective field theory. We clarify the Lyth bound which relates the tensor/scalar ratio and the value of phi/M_P.Effective field theory is valid for V(phi)<<M_P^4 allowing inflaton fields phi well beyond M_P.Thus, no bounds based on the value of phi/M_P appear on r.In case one wishes to restrict oneself to chi < 1 in the above potential, values for r larger than the previous bounds are allowed thanks to the new factor sqrt{N}~7.Finally,it is conjectured that inflation is hovering near a trivial gaussian infrared fixed point during slow roll. Slow-roll inflation would then enjoy universality properties as those of critical phenomena.

We provide a complete picture to the self-gravitating 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) -> infty, keeping N/V^{1/3} fixed. We compute the equation of state (we do not assume it as is customary in hydrodynamics), as well as the energy, free energy, entropy, chemical potential, specific heats, compressibilities 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 -> infty and V -> infty limit. The system is in a gaseous phase for eta < eta_T and collapses into a dense object for eta > eta_T in the CE with the pressure becoming large and negative. At eta \simeq eta_T the isothermal compressibility diverges and the gas collapses. Our Monte Carlo simulations yield eta_T \simeq 1.515. We find that PV/[NT] = f(eta). 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 gaseous phase. Both collapse phase transitions (in the CE and in the MCE) are of zeroth order since the Gibbs free energy has a jump at the transitions.

Slow-roll inflation can be studied as an effective field theory. The form of the inflaton potential consistent with the data is V(phi) = N M^4 w(phi/[sqrt{N} M_{Pl}]) where phi is the inflaton field, M is the inflation energy scale, and N ~ 50 the number of efolds. The dimensionless function w(chi) and field chi are O(1). This form of the potential encodes the slow-roll expansion as an expansion in 1/N.A The Hubble parameter, inflaton mass and non-linear couplings are of the see-saw form in terms of M/M_{Pl}. The quartic coupling is lambda~1/N (M/M_{Pl})^4. The smallness of the non-linear couplings is not a result of fine tuning but a natural consequence of the validity of the effective field theory and slow roll approximation. Quantum corrections to slow roll inflation are computed and turn to be an expansion in powers (H/M_{Pl})^2. The corrections to the inflaton effective potential and its equation of motion are computed, as well as the quantum corrections to the observable power spectra. The near scale invariance of the fluctuations introduces a strong infrared behavior naturally regularized by Delta=(n_s -1)/2+r/8. We consider scalar curvature and tensor perturbations as well as light scalars and Dirac fermions coupled to the inflaton.The subhorizon part is completely specified by the trace anomaly of the fields with different spins and is solely determined by the space-time geometry. This inflationary effective potential is strikingly different from the usual Minkowski space-time result.Quantum corrections to the power spectra are expressed in terms of the CMB observables. Trace anomalies (especially the graviton part) dominate these quantum corrections in a definite direction: they enhance the scalar curvature fluctuations and reduce the tensor fluctuations.

We investigate the dynamics of thermalization and the approach to equilibrium in the classical phi^4 theory in 3+1 spacetime dimensions. The non-equilibrium dynamics is studied by numerically solving the equations of motion in a light- cone-like discretization of the model for a broad range of initial conditions and energy densities.A smooth cascade of energy towards the ultraviolet is found to be the basic mechanism of thermalization.After an initial transient stage,at a time scale of several hundreds inverse masses,the squared of the field gradient becomes larger than the nonlinear term and a stage of universal cascade emerges. As the cascade progresses, the modes with higher wavenumbers exhibit weaker and weaker nonlinearities well described by the Hartree approximation while the infrared modes retain strong selfinteractions. Two timescales for equilibration appears.For k^2><phi^2>(t) we observe an effective thermalization with a time scale in the thousands of inverse masses and the Hartree approximation holds. By effective thermalization we mean that the observable acquires the equilibrium functional form with an effective time dependent temperature Teff, which slowly decreases with time. Infrared modes with k^2 < <phi^2>(t) equilibrate only by time scales in the millions of inverse masses. Infrared modes with k^2 < <phi^2>(t) equilibrate only by time scales in the millions.Virialization and the equation of state start to set much earlier than effective thermalization.The applicability of these results in quantum field theory for large occupation numbers and small coupling is analyzed.

Neutrino oscillations in the early Universe prior to the epoch of primordial nucleosynthesis is studied by implementing real time non-equilibrium field theory methods. We focus on two flavors of Dirac neutrinos, however, the formulation is general. We obtain the equations of motion for neutrino wavepackets of either chirality and helicity in the plasma allowing for CP asymmetry. Contributions non-local in space-time to the self-energy dominate over the asymmetry for $T \gtrsim 3-5 \mathrm{MeV}$ if the lepton and neutrino asymmetries are of the same order as the baryon asymmetry. We find a new contribution which cannot be interpreted as the usual effective potential. The mixing angles and dispersion relations in the medium depend on \emph{helicity}. We find that resonant transitions are possible in the temperature range $ 10 \lesssim T \ll 100 \mathrm{MeV} $. Near a resonance in the mixing angle, the oscillation time scale in the medium as compared to the vacuum is \emph{slowed-down} substantially for small vacuum mixing angle. The time scale of oscillations \emph{speeds-up} for off resonance high energy neutrinos for which the mixing angle becomes vanishingly small. The equations of motion reduce to the familiar oscillation formulae for negative helicity ultrarelativistic neutrinos, but include consistently both the \emph{mixing angle and the oscillation frequencies in the medium}. These equations of motion also allow to study the dynamics of right handed as well as positive helicity neutrinos.

Precise cosmological data from WMAP and forthcoming CMB experiments motivate the study of the quantum corrections to the slowroll inflationary parameters.We find the quantum (loop) corrections to the equations of motion of the classical inflaton, its quantum fluctuations and the Friedmann equation in general single field slow roll inflation.We implement a renormalized effective field theory EFT approach based on an expansion in (H/M_{Pl})^2 and slow roll parameters epsilon_V,eta_V,sigma_V, xi_V.We find that the leading order quantum corrections to the inflaton effective potential and its equation of motion are determined by the power spectrum of scalar fluctuations. Its near scale invariance introduces a strong infrared behavior naturally regularized by the slow roll parameter Delta = eta_V-epsilon_V=(n_s-1)/2+r/8.To leading order in the EFT and slow roll expansions we find V_{eff}(Phi_0)=V_R(Phi_0)[1+(Delta^2_T/32)(n_s-1+3r/8) /(n_s-1+r/4)+higher orders]where n_s and r=Delta^2_T/Delta^2_R are the CMB observables that depend implicitly on Phi_0, and V_R(Phi_0) is the renormalized classical inflaton potential.This effective potential during slow roll inflation is strikingly different from the Minkowski space-time result.Superhorizon scalar field fluctuations grow for late times eta -> 0^- as |\eta|^{-1+Delta-d_} where d_ is a novel quantum correction to the scaling exponent related to the self decay of superhorizon inflaton fluctuations eta is the conformal time. We generalize this to the case of the inflaton interacting with a light scalar field. These quantum corrections arising from interactions will compete with higher order slow-roll corrections and must be taken into account for the precision determination of inflationary parameters.

We obtain the effective inflaton potential during slow roll inflation by including the one loop quantum corrections to the energy momentum tensor from scalar curvature and tensor perturbations as well as quantum fluctuations from light scalars and light Dirac fermions generically coupled to the inflaton. During slow roll inflation there is a clean and unambiguous separation between superhorizon and subhorizon contributions to the energy momentum tensor. The superhorizon part is determined by the curvature perturbations and scalar field fluctuations: both feature infrared enhancements as the inverse of a combination of slow roll parameters which measure the departure from scale invariance in each case.Fermions and gravitons do not exhibit infrared divergences. The subhorizon part is completely specified by the trace anomaly of the fields with different spins and is solely determined by the space-time geometry. The one-loop quantum corrections to the amplitude of curvature and tensor perturbations are obtained to leading order in slow-roll and in the (H/M_PL)^2 expansion. This study provides a complete assessment of the backreaction problem up to one loop including bosonic and fermionic degrees of freedom. The result validates the effective field theory description of inflation and confirms the robustness of the inflationary paradigm to quantum fluctuations. Quantum corrections to the power spectra are expressed in terms of the CMB observables:n_s, r and dn_s/dln k. Trace anomalies (especially the graviton part) dominate these quantum corrections in a definite direction: they enhance the scalar curvature fluctuations and reduce the tensor fluctuations.

We clarify inflaton models by considering them as effective field theories in the Ginzburg-Landau spirit.In this new approach, the precise form of the inflationary potential is constructed from the present WMAP data, and a useful scheme is prepared to confront with the forthcoming data. In this approach, the WMAP statement excluding the pure phi^4 potential implies the presence of an inflaton mass term at the scale m sim 10^{13}GeV. Chaotic, new and hybrid inflation is studied in an unified way. In all cases the inflaton potential takes the form V(phi) = m^2 M_{Pl}^2 v(phi/M_{Pl}), where all coefficients in the polynomial v(x) are of order one. If such potential corresponds to super symmetry breaking, the susy breaking scale is sqrt{m M_{Pl}} \sim 10^{16}GeV which turns to coincide with the GUT scale. The inflaton mass is therefore given by a see-saw formula m sim M_{GUT}^2/M_{Pl}. The observables turn to be two valued functions: one branch corresponds to new inflation and the other to chaotic inflation,the branch point being the pure quadratic potential.For red tilt spectrum, the potential which fits the best the present data and which best prepares the way for the forthcoming data is a trinomial polynomial withnegative quadratic term (new inflation).For blue tilt spectrum, hybrid inflation turns to be the best choice. In both cases we find a formula relating the inflaton mass with the ratio r of tensor/scalar perturbations and the spectral index ns of scalar perturbations: 10^6 m/M_{Pl}= 127 sqrt{r|1-n_s|} ;(the coefficient 127 follows from the WMAP amplitude.Implications for string theory are discussed.

The self-gravitating gas in the presence of a positive cosmological constant Lambda is studied in thermal equilibrium by Monte Carlo simulations and by the mean field approach. We find excellent agreement between both approaches already for N = 1000 particles on a volume V [The mean field is exact in the infinite N limit]. The domain of stability of the gas is found to increase when the cosmological constant increases. The particle density is shown to be an increasing (decreasing) function of the distance when the dark energy dominates over self-gravity (and vice-versa).We confirm the validity of the thermodynamic limit: N, V -> infty with N/V^{1/3} and Lambda V^{2/3} fixed. In such dilute limit extensive thermodynamic quantities like energy, free energy, entropy turn to be proportional to N. We find that the gas is stable till the isothermal compressibility diverges. Beyond this point the gas becomes a extremely dense object whose properties are studied by Monte Carlo.

We clarify classical inflaton models by considering them as effective field theories `a la Ginzburg-Landau. In this approach, the WMAP statement excluding the pure phi^4 potential implies the presence of an inflaton mass term at the scale m ~ 10^{13}GeV.Chaotic, new and hybrid inflation models are studied in an unified manner.In all cases the inflaton potential takes the form V(phi)=m^2 M_{Pl}^2 v(phi/M_{Pl}), where all coefficients in v(varphi) are of order one.If such potential corresponds to supersymmetry breaking, the susy scale is sqrt{m M_{Pl}}~10^{16}GeV which coincides with the GUT scale. For red tilted spectrum new inflation turns to be favoured while for blue tilted spectrum, hybrid inflation turns to be the best choice.In both cases, we find an analytic formula relating the inflaton mass with r (tensor/scalar ratio) and the scalar spectral index ns: 10^6 m/M_{Pl} = 127 sqrt{r|1-n_s|} where the numerical coefficient follows from the WMAP amplitude. We review quantum phenomena during inflation which contribute to relevant observables in the CMB anisotropies and polarization and we focus on inflaton decay.The deviation from the scale invariant power spectrum is measured by a small parameter Delta.Delta regulates the infrared divergences too.In slow roll inflation, Delta is a simple function of the slow roll parameters. The quantum fluctuations can selfdecay in the inflationary expansion through processes forbidden in Minkowski. We compute the self-decay of the inflaton quantum fluctuations during slow roll:superhorizon fluctuations decay as eta^{Gamma} in conformal time where the new scaling dimension Gamma is expressed in terms of the amplitude, n_s and r. we discuss the implications for scalar and tensor perturbations as well as for non gaussianities.

The universe cools down monotonically following its expansion.This generates a sequence of phase transitions. If a second order phase transition happens during the radiation dominated era with a charged order parameter, spinodal unstabilities generate large numbers of charged particles. These particles hence produce magnetic fields.We use out of equilibrium field theory methods to study the dynamics in a mean field or large N setup.The dynamics after the transition features two distinct stages: a spinodal regime dominated by linear long wave length instabilities, and a scaling stage in which the non-linearities and backreaction of the scalar fields are dominant. This second stage describes the growth of horizon sized domains. We implement a formulation based on the non equilibrium Schwinger-Dyson equations to obtain the spectrum of magnetic fields that includes the dissipative effects of the plasma. We find that large scale magnetogenesis is efficient during the scaling regime. Charged scalar field fluctuations with wavelengths of the order of the Hubble radius induce large scale magnetogenesis via loop effects.The leading processes are:pair production, pair annihilation and low energy bremsstrahlung, these processes while forbidden in equilibrium are allowed strongly out of equilibrium. The ratio between the energy density on scales larger than L and that in the background radiation r(L,T)= rho_B(L,T)/rho_{cmb}(T) is r(L,T) ~ 10^{-34} at the Electroweak scale and r(L,T) ~ 10^{-14} at the QCD scale for L sim 1 Mpc. The resulting spectrum is insensitive to the magnetic diffusion length and equipartition between electric and magnetic fields does not hold. We conjecture that a similar mechanism could be operative after the QCD chiral phase transition.

Particle decay during inflation is studied by implementing a dynamical renormalization group resummation combined with a small Delta expansion. Delta measures the deviation from the scale invariant power spectrum and regulates the infrared. In slow roll inflation, Delta is a simple function of the slow roll parameters epsilon_V, eta_V.We find that quantum fluctuations can self-decay as a consequence of the inflationary expansion through processes which are forbidden in Minkowski space-time. We compute the self-decay of the inflaton quantum fluctuations during slow roll inflation.For wavelengths deep inside the Hubble radius the decay is enhanced by the emission of ultrasoft collinear quanta, i.e. bremsstrahlung radiation of superhorizon quanta which becomes the leading decay channel for physical wavelengths H<<k_{ph}(eta)<<H/(eta_V-eps_V). The decay of short wavelength fluctuations hastens as the physical wave vector approaches the horizon. Superhorizon fluctuations decay with a power law eta^Gamma in conformal time where in terms of the amplitude of curvature perturbations Delta^2_R, the scalar spectral index n_s, the tensor to scalar ratio r and slow roll parameters: Gamma \simeq [32 xi^2_V Delta^2_R]/ /(n_s-1+r/4)^2.The behavior of the growing mode eta^{eta_V-epsilon_V+Gamma}/eta features an anomalous scaling dimension Gamma. We discuss the implications of these results for scalar and tensor perturbations and for non-gaussianities in the power spectrum. The recent WMAP data suggests Gamma >3.6 10^{-9}.

We compute the primordial scalar, vector and tensor metric perturbations arising from quantum field inflation. Quantum field inflation takes into account the nonperturbative quantum dynamics of the inflaton consistently coupled to the dynamics of the (classical) cosmological metric. For chaotic inflation, the quantum treatment avoids the unnatural requirements of an initial state with all the energy in the zero mode. For new inflation it allows a consistent treatment of the explosive particle production due to spinodal instabilities. Quantum field inflation (under conditions that are the quantum analog of slow roll) leads, upon evolution, to the formation of a condensate starting a regime of effective classical inflation. We compute the primordial perturbations taking the dominant quantum effects into account. The results for the scalar, vector and tensor primordial perturbations are expressed in terms of the classical inflation results. For a N-component field in a O(N) symmetric model, adiabatic fluctuations dominate while isocurvature or entropy fluctuations are negligible. The results agree with the current WMAP observations and predict corrections to the power spectrum in classical inflation.Such corrections are estimated to be of the order of m^2/[N H^2] where m is the inflaton mass and H the Hubble constant at horizon crossing. This turns to be about 4% for the cosmologically relevant scales. This quantum field treatment of inflation provides the foundations to the classical inflation and permits to compute quantum corrections to it.

We investigate the relaxation and decay of a particle during inflation by implementing the dynamical renormalization group. This investigation allows us to give a meaningful definition for the decay rate in an expanding universe. As a prelude to a more general scenario, the method is applied here to study the decay of a particle in de Sitter inflation via a trilinear coupling to massless conformally coupled particles, both for wavelengths much larger and much smaller than the Hubble radius. For superhorizon modes we find that the decay is of the form eta^{Gamma1} with eta being conformal time and we give an explicit expression for Gamma1 to leading order in the coupling which has a noteworthy interpretation in terms of the Hawking temperature of de Sitter space-time. We show that if the mass M of the decaying field is << H then the decay rate during inflation is enhanced over the Minkowski spacetime result by a factor 2H/[pi M]. For wavelengths much smaller than the Hubble radius we find that the decay law is e^{-alpha/[k H C(eta)} with C(eta) the scale factor and alpha determined by the strength of the trilinear coupling. This result suggests a suppression of power for long wavelength modes upon horizon crossing. In all cases we find a substantial enhancement in the decay law as compared to Minkowski space-time. These results suggest potential implications for the spectrum of scalar density fluctuations as well as non-gaussianities.

We develop the cluster expansion and the Mayer expansion for the self-gravitating thermal gas and prove the existence and stability of the thermodynamic limit N, V to infty with N/V^{1/3} fixed. The essential (dimensionless) variable is here eta = [G m^2 N]/[V^{1/3} T] (which is kept fixed in the thermodynamic limit). We succeed in this way to obtain the expansion of the grand canonical partition function in powers of the fugacity. The corresponding cluster coefficients behave in the thermodynamic limit as [eta/N]^{j-1} c_j where c_j are pure numbers. They are expressed as integrals associated to tree cluster diagrams. A bilinear recurrence relation for the coefficients c_j is obtained from the mean field equations in the Abel form. In this way the large j behaviour of the c_j is calculated. This large j behaviour provides the position of the nearest singularity which corresponds to the critical point (collapse) of the self-gravitating gas in the grand canonical ensemble. Finally, we discuss why other attempts to define a thermodynamic limit for the self-gravitating gas fail.

The non-relativistic self-gravitating gas in thermal equilibrium in the presence of a positive cosmological constant Lambda (dark energy) is investigated. The dark energy introduces a force pushing outward all particles with strength proportional to their distance to the center of mass.We consider the statistical mechanics of the self-gravitating gas of N particles in a volume V at thermal equilibrium in the presence of Lambda. It is shown that the thermodynamic limit exists and is described by the mean field equations provided N, V -> inf with N/V^{1/3} fixed and Lambda V^{2/3} fixed. That is, Lambda -> 0 for N, V -> inf. The case of Lambda fixed and N, V -> inf is solved too. We solve numerically the mean field equation for spherical symmetry obtaining an isothermal sphere for Lambda>0. The particle distribution turns to flatten compared with the Lambda = 0 case.Moreover, the particle density increases with the distance when the cosmological constant dominates. There is a bordering case with uniform density.The density contrast between the center and the boundary is significatively reduced. In addition, the critical point associated to the collapse (Jeans') phase transition is pushed towards higher values of N/[T V^{1/3}] by the presence of Lambda > 0.The nature and the behaviour near the critical points is not affected by the presence of Lambda>0.

We study the dynamics of thermalization and the approach to equilibrium in the classical phi^4 theory in 1+1 spacetime dimensions. At thermal equilibrium we exploit the equivalence between the classical canonical averages and transfer matrix quantum traces of the anharmonic oscillator to obtain exact results for the temperature dependence of several observables, which provide a set of criteria for thermalization. We find that the Hartree approximation is remarkably accurate in equilibrium. The non-equilibrium dynamics is studied by numerically solving the equations of motion in light-cone coordinates for a broad range of initial conditions and energy densities.The time evolution is described by several stages with a cascade of energy towards the ultraviolet. After a transient stage, the spatio-temporal gradient terms become larger than the nonlinear term and a stage of universal cascade emerges.This cascade starts at a time scale t_0 independent of the initial conditions (except for very low energy density). Here the power spectra feature universal scaling behavior and the front of the cascade k(t) grows as a power law k(t) sim t^alpha with alpha lesssim 0.25. The wake behind the cascade is described as a state of Local Thermodynamic Equilibrium (LTE) with all correlations being determined by the equilibrium functional form with an effective time dependent temperatureTeff(t) which slowly decreases as sim t^{-alpha}.Two well separated time scales emerge while Teff(t) varies slowly, the wavectors in the wake with k < k(t) attain LTE on much shorter time scales.This universal scaling stage ends when the front of the cascade reaches the cutoff at a time t_1 sim a^{-1/alpha}. Virialization starts to set much earlier than LTE. We find that strict thermalization is achieved only for an infinite time scale.

This is a pedagogical introduction to early cosmology and the host of fundamental physics involved in it (particle physics, grand unification and general relativity). Inflation and the inflaton field are the central theme of this review. The quantum field treatment of the inflaton is presented including its out of equilibrium evolution and the use of nonperturbative methods. The observational predictions for the CMB anisotropies are briefly discussed. Finally, open problems and future perspectives in connection with dark energy and string theory are overviewed.

We study the generation of large scale primordial magnetic fields by a cosmological phase transition during the radiation dominated era. The setting is a theory of N charged scalar fields coupled to an abelian gauge field, that undergoes a phase transition at a critical temperature much larger than the electroweak scale. The dynamics after the transition features two distinct stages: a spinodal regime dominated by linear long-wavelength instabilities, and a scaling stage in which the non-linearities and backreaction of the scalar fields are dominant.This second stage describes the growth of horizon sized domains.We implement a formulation based on the non-equilibrium Schwinger-Dyson equations to obtain the spectrum of magnetic fields that includes the dissipative effects of the plasma. We find that large scale magnetogenesis is efficient during the scaling regime.Charged scalar field fluctuations with wavelengths of the order of the Hubble radius induce large scale magnetogenesis via loop effects. The leading processes are: pair production, pair annihilation and low energy bremsstrahlung, these processes while forbidden in equilibrium are allowed strongly out of equilibrium.The ratio between the energy density on scales larger than L and that in the background radiation r(L,T)= rho_B(L,T)/ rho_{cmb}(T) is r(L,T) sim 10^{-34} at the Electroweak scale and r(L,T) sim 10^{-14} at the QCD scale for L sim 1 Mpc.The resulting spectrum is insensitive to the magnetic diffusion length and equipartition between electric and magnetic fields does not hold. We conjecture that a similar mechanism could be operative after the QCD chiral phase transition.

We study the generation of large scale primordial magnetic fields by a cosmological phase transition during the radiation dominated era. The setting is a theory of N charged scalar fields coupled to an abelian gauge field, that undergoes a phase transition at a critical temperature much larger than the electroweak scale. The dynamics after the transition features two distinct stages: a spinodal regime dominated by linear long-wavelength instabilities, and a scaling stage in which the non-linearities and backreaction of the scalar fields are dominant. This second stage describes the growth of horizon sized domains. We implement a recently introduced formulation to obtain the spectrum of magnetic fields that includes the dissipative effects of the plasma. We find that large scale magnetogenesis is very efficient during the scaling regime. The ratio between the energy density on scales larger than L and that in the background radiation r(L,T) = rho_B(L,T)/rho_{cmb}(T) is r(L,T) \sim 10^{-34} at the Electroweak scale and r(L,T) \sim 10^{-14} at the QCD scale for L \sim 1 Mpc. The resulting spectrum is insensitive to the magnetic diffusion length. We conjecture that a similar mechanism could be operative after the QCD chiral phase transition.

We present an overview on extreme energy cosmic rays (EECR) and the fundamental physics connected with them. The top-down and bottom-up scenarii are contrasted. We summarize the essential features underlying the top-down scenarii for EECR, namely, the lifetime and the mass {\bf imposed} to the heavy relics whatever they be: topological and non-topological solitons, X-particles, cosmic defects, microscopic black-holes, fundamental strings. An unified formula for the quantum decay rate of all these objects was provided in hep-ph/0202249. The key point in the top-down scenarii is the necessity to {\bf adjust} the lifetime of the heavy object to the age of the universe. The natural lifetimes of such heavy objects are, however, microscopic times associated to the GUT energy scale (sim 10^{-28} sec. or shorter); such heavy objects could have been abundantly formed by the end of inflation and it seems natural they decayed shortly after being formed. The arguments produced to {\bf fine tune} the relics lifetime to the age of the universe are critically analyzed. The annihilation scenario (`Wimpzillas') is analyzed too. Top-down scenarii based on networks of topological defects are strongly disfavored at the light of the recent CMB anisotropy observations. We discuss the acceleration mechanisms of cosmic rays,their possible astrophysical sources and the main open physical problems and difficulties in the context of bottom-up scenarii, and we conclude by outlining the expectations from future observatories like EUSO and where the theoretical effort should be placed.

We study the generation of magnetic fields during the stage of particle production resulting from spinodal instabilities during phase transitions out of equilibrium. The main premise is that long-wavelength instabilities that drive the phase transition lead to strong non-equilibrium charge and current fluctuations which generate electromagnetic fields. We present a formulation based on the non-equilibrium Schwinger-Dyson equations that leads to an exact expression for the spectrum of electromagnetic fields valid for general theories and cosmological backgrounds and whose main ingredient is the transverse photon polarization out of equilibrium. This formulation includes the dissipative effects of the conductivity in the medium. As a prelude to cosmology we study magnetogenesis in Minkowski space-time in a theory of N charged scalar fields to lowest order in the gauge coupling and to leading order in the large N within two scenarios of cosmological relevance. The long-wavelength power spectrum for electric and magnetic fields at the end of the phase transition is obtained explicitly. It follows that equipartition between electric and magnetic fields does not hold out of equilibrium. In the case of a transition from a high temperature phase, the conductivity of the medium severely hinders the generation of magnetic fields, however the magnetic fields generated are correlated on scales of the order of the domain size, which is much larger than the magnetic diffusion length. Implications of the results to cosmological phase transitions driven by spinodal unstabilities are discussed.

The physics during the inflationary stage of the universe is of quantum nature involving extremely high energy densities. Moreover, it is out of equilibrium on a fastly expanding dynamical geometry. We complement here the 1999 Chalonge Lectures on out of equilibrium fields in self-consistent inflationary dynamics [astro-ph/0006446] investigating inflation driven by the evolution of highly excited quantum states. These states are characterized by a non-perturbatively large number of quanta in a band of momenta and with zero or nonzero expectation value of the inflaton scalar field. They represent the situation in which initially a non-perturbatively large energy density is localized in a band of high energy quantum modes and are coined tsunami-waves. The self- consistent evolution of this quantum state and the scale factor is studied analytically and numerically. It is shown that the time evolution of these quantum states lead to two consecutive stages of inflation under conditions that are the quantum analogue of slow-roll. The evolution of the scale factor during the first stage has new features that are characteristic of the quantum state. During this initial stage the quantum fluctuations in the highly excited band build up an effective homogeneous condensate with a non-perturbatively large amplitude as a consequence of the large number of quanta.The second stage of inflation is similar to the usual classical chaotic scenario but driven by this effective condensate.The excited quantum modes are superhorizon in the first stage and do not affect the power spectrum of scalar perturbations.Thus, this tsunami quantum state provides a field theoretical justification for chaotic scenarios driven by a classical homogeneous scalar field of large amplitude.

The essential features underlying the top-down scenarii for UHECR are discussed, namely, the stability (or lifetime) imposed to the heavy objects (particles) whatever they be: topological and non-topological solitons, X-particles, cosmic defects, microscopic black-holes, fundamental strings. We provide an unified formula for the quantum decay rate of all these objects as well as the particle decays in the standard model. The key point in the top-down scenarii is the necessity to adjust the lifetime of the heavy object to the age of the universe. This ad-hoc requirement needs a very high dimensional operator to govern its decay and/or an extremely small coupling constant. The natural lifetimes of such heavy objects are, however, microscopic times associated to the GUT energy scale (sim 10^{-28} sec. or shorter). It is at this energy scale (by the end of inflation) where they could have been abundantly formed in the early universe and it seems natural that they decayed shortly after being formed.

We study the statistical mechanics of the self-gravitating gas at thermal equilibrium with two kinds of particles. We start from the partition function in the canonical ensemble which we express as a functional integral over the densities of the two kinds of particles for a large number of particles. The system is shown to possess an infinite volume limit when (N_1,N_2,V)->infty, keeping N_1/V^{1/3} and N_2/V^{1/3} fixed. The saddle point approximation becomes here exact for (N_1,N_2,V)->infty.It provides a nonlinear differential equation on the particle densities. For the spherically symmetric case, we compute the densities as functions of two dimensionless physical parameters: eta_1=G m_1^2 N_1/[V^{1/3} T] and eta_2=G m_2^2 N_2/[V^{1/3} T] (where G is Newton's constant, m_1 and m_2 the masses of the two kinds of particles and T the temperature). According to the values of eta_1 and eta_2 the system can be either in a gaseous phase or in a highly condensed phase.The gaseous phase is stable for eta_1 and eta_2 between the origin and their collapse values. The gas is inhomogeneous and the mass M(R) inside a sphere of radius R scales with R as M(R) propto R^d suggesting a fractal structure. The value of d depends in general on eta_1 and eta_2 except on the critical line for the canonical ensem- ble where it takes the universal value d simeq 1.6 for all values of N_1/N_2. The equation of state is computed.It is found to be locally a perfect gas equation of state. Thermodynamic functions are computed as functions of eta_1 and eta_2. They exhibit a square root Riemann sheet with the branch points on the critical canonical line. This treatment is further generalized to the self-gravitating gas with n-types of particles.

We investigate inflation driven by the evolution of highly excited quantum states within the framework of out of equilibrium field dynamics. These states are characterized by a non-perturbatively large number of quanta in a band of momenta but with vanishing expectation value of the scalar field.They represent the situation in which initially a non-perturbatively large energy density is localized in a band of high energy quantum modes and are coined tsunami-waves. The self-consistent evolution of this quantum state and the scale factor is studied analytically and numerically. It is shown that the time evolution of these quantum states lead to two consecutive stages of inflation under conditions that are the quantum analogue of slow-roll. The evolution of the scale factor during the first stage has new features that are characteristic of the quantum state. During this initial stage the quantum fluctuations in the highly excited band build up an effective homogeneous condensate with a non- perturbatively large amplitude as a consequence of the large number of quanta. The second stage of inflation is similar to the usual classical chaotic scenario but driven by this effective condensate.The excited quantum modes are already superhorizon in the first stage and do not affect the power spectrum of scalar perturbations. Thus, this tsunami quantum state provides a field theoretical justification for chaotic scenarios driven by a classical homogeneous scalar field of large amplitude.

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 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 investigate the out of equilibrium dynamics of global chiral supersymmetry at finite energy density. We concentrate on two specific models. The first is the massive Wess-Zumino model which we study in a selfconsistent one-loop approximation. We find that for energy densities above a certain threshold, the fields are driven dynamically to a point in field space at which the fermionic component of the superfield is massless. The state, however is found to be unstable, indicating a breakdown of the one-loop approximation. To investigate further, we consider an O(N) massive chiral model which is solved exactly in the large $N$ limit. For sufficiently high energy densities, we find that for late times the fields reach a nonperturbative minimum of the effective potential degenerate with the perturbative minimum. This minimum is a true attractor for O(N) invariant states at high energy densities, and this provides a mechanism for determining which of the otherwise degenerate vacua is chosen by the dynamics. The final state for large energy density is a cloud of massless particles (both bosons and fermions) around this new nonperturbative supersymmetric minimum. By introducing boson masses which softly break the supersymmetry, we demonstrate a see-saw mechanism for generating small fermion masses. We discuss some of the cosmological implications of our results.

The non-equilibrium dynamics of a Yukawa theory with N fermions coupled to a scalar field is studied in the large N limit with the goal of comparing the dynamics predicted from the renormalization group improved effective potential to that obtained including the fermionic backreaction. The effective potential is of the Coleman-Weinberg type. Its renormalization group improvement is unbounded from below and features a Landau pole. When viewed self-consistently, the initial time singularity does not arise. The different regimes of the dynamics of the fully renormalized theory are studied both analytically and numerically. Despite the existence of a Landau pole in the model, the dynamics of the mean field is smooth as it passes the location of the pole. This is a consequence of a remarkable cancellation between the effective potential and the dynamical chiral condensate. The asymptotic evolution is effectively described by a quartic upright effective potential. In all regimes, profuse particle production results in the formation of a dense fermionic plasma with occupation numbers nearly saturated up to a scale of the order of the mean field. This can be interpreted as a chemical potential. We discuss the implications of these results for cosmological preheating.

In normal degenerate quark matter, the exchange of dynamically screened transverse gluons introduces infrared divergences in the quark self-energies that lead to the breakdown of the Fermi liquid description. If the core of neutron stars are composed of quark matter with a normal component, cooling by direct quark Urca processes may be modified by non-Fermi liquid corrections. We find that while the quasiparticle density of states is finite and non-zero at the Fermi surface, its frequency derivative diverges and results in non-Fermi liquid corrections to the specific heat of the normal, degenerate component of quark matter. We study these non-perturbative non-Fermi liquid corrections to the specific heat and the temperature dependence of the chemical potential and show that these lead to a reduction of the specific heat.

This work studies the hydrodynamics of self-gravitating compressible isothermal fluids. We show that the hydrodynamic evolution equations in absence of viscosity are scale covariant. We study the evolution of the time dependent fluctuations around singular and regular isothermal spheres. We linearize the fluid equations around such stationary solutions and apply Laplace transform to solve them. We find that the system is stable below a critical size (X ~ 9.0 in dimensionless variables) and unstable above; this size is the same critical size found in the study of the thermodynamical stability in the canonical ensemble and associated to a center-to-border density ratio of 32.1 . We prove that the value of this critical size is independent of the Reynolds number of the system. Furthermore, we give a detailed description of the series of successive dynamical instabilities that appear at higher and higher sizes following the geometric progression X_n ~ 10.7^n. We turn then to study exact solutions of the hydrodynamic equations without viscosity and we provide analytic and numerical axisymmetric soliton-type solutions. The stability of exact solutions corresponding to a collapsing filament is studied by computing linear fluctuations. Radial fluctuations growing faster than the background are found for all sizes of the system. However, a critical size (X ~ 4.5) appears, separating a weakly from a strongly unstable regime.

Infrared divergences from the exchange of dynamically screened magnetic gluons (photons) lead to the breakdown of the Fermi liquid description of the {\em normal} state of cold and dense QCD and QED. We implement a resummation of these divergences via the renormalization group to obtain the spectral density, dispersion relation, widths and wave function renormalization of single quasiparticles near the Fermi surface. We find that all feature scaling with anomalous dimensions: $\omega_p({k}) \propto |k-k_F|^{\frac{1}{1-2\lambda}} ~ ; ~ \Gamma(k) \propto |k-k_F|^{\frac{1}{1-2\lambda}} ~;~ Z_p({k}) \propto |k-k_F|^{\frac{2\lambda}{1-2\lambda}}$ with $\lambda = \frac{\alpha}{6\pi} ~ {for QED} \vspace{0.5 ex} ~,~ \frac{\alpha_s}{6\pi} \frac{N^2_c-1}{2N_c} \~~{for QCD with}$. The discontinuity of the distribution function for quasiparticles near the Fermi surface vanishes. The dynamical renormalization group is implemented to study the relaxation of quasiparticles in real time. Quasiparticles with Fermi momentum have vanishing group velocity and relax with a power law with a coupling dependent anomalous dimension.

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 physics during the inflationary stage of the universe is of quantum nature involving extremely high energy densities. Moreover, it is out of equilibrium on a fastly expanding dynamical geometry.We present in these lectures non-perturbative out of equilibrium field theoretical methods in cosmological universes. We then study the non-linear dynamics of quantum fields in matter and radiation dominated FRW and de Sitter universes. We investigate the explosive particle production due to spinodal instabilities and parametric amplification in FRW and deSitter universes with and without symmetry breaking. We show how the particle production is sensitive to the expansion of the universe.We present a complete renormalization scheme for the equation of motion and the energy momentum tensor in flat cosmologies. We then consider an O(N) inflaton model coupled self-consistently to gravity in the semiclassical approximation, with `new inflation' type initial conditions. We study the dynamics self-consistently and non-perturbatively with non-equilibrium field theory methods in the large N limit. We find that spinodal instabilities drive the growth of non-perturbatively large quantum fluctuations which shut off the inflationary growth of the scale factor. A very specific combination of these large quantum fluctuations plus the inflaton zero mode assemble into a new effective field. This new field behaves classically and it is the object which actually rolls down. The metric perturbations during inflation are computed using this effective field and the Bardeen variable for superhorizon modes during inflation. We compute the amplitude and index for the spectrum of scalar density and tensorperturbations and find for these models that the spinodal instabilities are responsible for a `red' primordial spectrum.

Critical slowing down of the relaxation of the order parameter is relevant both in early the universe and in ultrarelativistic heavy ion collisions. We study the relaxation rate of the order parameter in an O(N) scalar theory near the critical point to model the non-equilibrium dynamics of critical fluctuations near the chiral phase transition.A lowest order perturbative calculation (two loops in the coupling lambda) reveals the breakdown of perturbation theory for long-wavelength fluctuations in the critical region and the emergence of a hierarchy of scales with hard q>T, semisoft T >> q >> lambda T and soft lambda T>>q loop momenta which are widely separated for weak coupling. The non-perturbative resummation in the large N limit reveals the renormalization of the interaction and the crossover to an effective 3D-theory for soft momenta.The effective 3D coupling goes to the Wilson-Fisher 3D fixed point in the soft limit.The relaxation rate of the order parameter for wave vectors lambda T >>k>> k_{us} or near the critical temperature lambda T>>m_T>>k_{us} with the ultra soft scale k_{us} = [(lambda T)/(4pi)] exp[-(4pi/lambda)] is dominated by classical semisoft loop momentum leading to Gamma(k,T) = lambda T/(2 pi N). For wavectors k<< k_{us} the damping rate is dominated by hard loop momenta and given by Gamma(k,T)=4 pi T/[3N ln(T/k)]. Analogously, for homogeneous fluctuations in the ultracritical region m_T<<k_{us} the damping rate is given by Gamma_0(m_T,T)=4 pi T/[3N ln(T/m_T)]. Thus critical slowing down emerges for ultrasoft fluctuations where the rate is lambda-independent. The strong coupling regime and the shortcomings of the quasiparticle interpretation are discussed.

The dynamics of a symmetry breaking phase transition is studied in a radiation and matter dominated spatially flat FRW cosmology in the large N limit of a scalar field theory.The quantum density matrix is evolved from an initial state of quasiparticles in thermal equilibrium at a temperature higher than the critical. The cosmological expansion decreases the temperature and triggers the phase transition. We identify three different time scales: an early regime dominated by linear instabilities and the exponential growth of long-wavelength fluctuations,an intermediate scale when the field fluctuations probe the broken symmetry states and an asymptotic scale wherein a scaling regime emerges for modes of wavelength comparable to or larger than the horizon.The scaling regime is characterized by a dynamical physical correlation length xi_{phys} = d_H(t) with d_H(t) the size of the causal horizon, thus there is one correlated region per causal horizon. Inside these correlated regions the field fluctuations sample the broken symmetry states. The amplitude of the long-wavelength fluctuations becomes non-perturbatively large due to the early times instabilities and a semiclassical but stochastic description emerges in the asymptotic regime. In the scaling regime, the power spectrum is peaked at zero momentum revealing the onset of a Bose-Einstein condensate.The scaling solution results in that the equation of state of the scalar fields is the same as that of the background fluid. This implies a Harrison-Zeldovich spectrum of scalar density perturbations for long-wavelengths. We discuss the corrections to scaling as well as the universality of the scaling solution and the differences and similarities with the classical non-linear sigma model.

The dynamics of symmetry breaking during out of equilibrium phase transitions is a topic of great importance in many disciplines, from condensed matter to particle physics and early Universe cosmology with definite experimental impact. In these notes we provide a summary of the relevant aspects of the dynamics of symmetry breaking in many different fields with emphasis on the experimental realizations. In condensed matter we address the dynamics of phase ordering, the emergence of condensates, coarsening and dynamical scaling. In QCD the possibility of disoriented chiral condensates of pions emerging during a strongly out of equilibrium phase transition is discussed. We elaborate on the dynamics of phase ordering in phase transitions in the Early Universe, in particular the emergence of condensates and scaling in FRW cosmologies. We mention some experimental efforts in different fields that study this wide ranging phenomena and offer a quantitative theoretical description both at the phenomenological level in condensed matter, introducing the scaling hypothesis as well as at a microscopic level in quantum field theories. The emergence of semiclassical condensates and a dynamical length scale is shown in detail, in quantum field theory this length scale is constrained by causality. The large N limit provides a natural bridge to compare the solutions in different settings and to establish similarities and differences.

Non-singular two and three dimensional string cosmologies are constructed using the exact conformal field theories corresponding to SO(2,1)/SO(1,1) and SO(2,2)/SO(2,1). All semi-classical curvature singularities are canceled in the exact theories for both of these cosets, but some new quantum curvature singularities emerge. However, considering different patches of the global manifolds, allows the construction of non-singular spacetimes with cosmological interpretation. In both two and three dimensions, we construct non-singular oscillating cosmologies, non-singular expanding and inflationary cosmologies including a de Sitter (exponential) stage with positive scalar curvature as well as non-singular contracting and deflationary cosmologies. We analyse these cosmologies in detail with respect to the behaviour of the scale factors, the scalar curvature and the string-coupling. The sign of the scalar curvature is changed by the quantum corrections in oscillating cosmologies and evolves with time in the non-oscillating cases. Similarities between the two and three dimensional cases suggest a general picture for higher dimensional coset cosmologies: Anisotropy seems to be a generic unavoidable feature, cosmological singularities are generically avoided and it is possible to construct non-singular cosmologies where some spatial dimensions are experiencing inflation while the others experience deflation.

These lectures address the dynamics of phase ordering out of equilibrium in condensed matter and in quantum field theory in cosmological settings, emphasizing their similarities and differences. In condensed matter we describe the phenomenological approach based on the Time Dependent Ginzburg-Landau (TDGL) description. After a general discussion of the main experimental and theoretical features of phase ordering kinetics and the description of linear (spinodal) instabilities we introduce the scaling hypothesis and show how a dynamical correlation length emerges in the large N limit in condensed matter systems. The large N approximation is a powerful tool in quantum field theory that allows the study of non-perturbative phenomena in a consistent manner. We study the exact solution to the dynamics after a quench in this limit in Minkowski space time and in radiation dominated Friedman-Robertson-Walker Cosmology. There are some remarkable similarities between these very different settings such as the emergence of a scaling regime and of a dynamical correlation length at late times that describe the formation and growth of ordered regions. In quantum field theory and cosmology this length scale is constrained by causality and its growth in time is also associated with coarsening and the onset of a condensate. We provide a density matrix interpretation of the formation of defects and the classicalization of quantum fluctuations.

We analyze the non-equilibrium dynamics of the O(N) Phi^4 model in the large N limit and for states of large energy density. The dynamics is dramatically different when the energy density is above the top of the tree level potential V_0 than when it is below it.When the energy density is below V_0, we find that non-perturbative particle production through spinodal instabilities provides a dynamical mechanism for the Maxwell construction. The asymptotic values of the order parameter only depend on the initial energy density and all values between the minima of the tree level potential are available, the asymptotic dynamical `effective potential' is flat between the minima. When the energy density is larger than V_0, the evolution samples ergodically the broken symmetry states, as a consequence of non-perturbative particle production via parametric amplification. Furthermore, we examine the quantum dynamics of phase ordering into the broken symmetry phase and find novel scaling behavior of the correlation function. There is a crossover in the dynamical correlation length at a time scale t_s \sim \ln(1/lambda). For t < t_s the dynamical correlation length \xi(t) \propto \sqrt{t} and the evolution is dominated by spinodal instabilities, whereas for t>t_s the evolution is non-linear and dominated by the onset of non-equilibrium Bose-Einstein condensation of long-wavelength Goldstone bosons.In this regime a true scaling solution emerges with a non- perturbative anomalous scaling length dimension z=1/2 and a dynamical correlation length \xi(t) \propto (t-t_s). The equal time correlation function in this scaling regime vanishes for r>2(t-t_s) by causality. For t > t_s the equal time correlation function falls of as 1/r. A semiclassical but stochastic description emerges for time scales t > t_s.

We study the generation of gravitational radiation by sources moving in the de Sitter background. Exploiting the maximal symmetry and the conformal flatness of de Sitter space-time we prove that the derivation of this gravitational radiation can be done along the same lines as in Minkowski space-time. A gauge is chosen in which all the physical and unphysical modes of the graviton are those of a minimally coupled massless scalar field in de Sitter space-time and a massless field in Minkowski space-time. The graviton retarded Green's function and the Schwinger commutator function are computed in this gauge using Quantum Field Theory techniques. We obtain closed formulae for the spectral decomposition in frequencies of the linaearized gravitational field produced by the source, in terms of a suitable spectral decomposition of the source energy-momentum tensor. This spectral decomposition is dictated by the free (sourceless) gravitational wave modes in the de Sitter background

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) and large structure of the universe. 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 scalar field phi(x) with exponential self-interaction. We analyze this field theory perturbatively and non-perturbatively through the renormalization group(RG).We show scaling behaviour (critical) for a continuous range of the physical parameters as the temperature. We derive in this framework the scaling relation M(R) \sim R^{d_H} for the mass on a region of size R, and Delta v \sim R^\frac12(d_H -1) for the velocity dispersion. For the density-density correlations we find a power-law behaviour for large distances \sim |r_1 - r_2|^{2D - 6}.The fractal dimension D turns to be related with the critical exponent \nu by D = 1/ \nu. Mean field theory yields \nu = 1/2, D = 2. Both the Ising and the mean field values are compatible with the present ISM observational data:1.4\leq D \leq 2. We develop a field theoretical approach to the galaxy distribution considering a gas of self-gravitating masses on the FRW background, in quasi-thermal equi- librium. We show that it exhibits scaling behaviour by RG methods. The galaxy correlations are computed without assuming homogeneity. We find <\rho({\vec r_0})\rho({\vec r_0} + {\vec r}) > \sim r^{D-3} $. The theory allows to compute the three and higher density correlators without any assumption.We find that the connected N-points density scales as r_1^{N(D-3)}, when $ r_1 >> r_i, 2\leq i \leq N $. There are no free parameters in this theory.

The energy and time scales during the inflationary stage of the universe calls for an out of equilibrium quantum field treatment. Moreover, the high energy densities involved make necessary the use of non-perturbative approaches as large N and Hartree methods. We start these lectures by introducing such non-perturbative out of equilibrium methods in cosmological universes. We discuss the renormalization procedure and the choice of initial conditions. We then study the nonlinear dynamics of quantum fields in matter and radiation dominated FRW and de Sitter universes. For a variety of initial conditions, we compute the evolution of the inflaton,its quantum fluctuations and the equation of state. We investigate the explosive particle production due to spinodal unstabilities and parametric amplification in FRW and de Sitter universes with and without symmetry breaking.We find that the particle production is sensitive to the expansion of the universe.For symmetry breaking scenarios, we determine generic late time fields behavior for FRW and deSitter cosmologies.We find that quantum fluctuations damp in FRW as the square of the scale factor while the order parameter approaches a minimum of the potential at the same rate.We con- sider an O(N) inflaton model coupled self-consistently to gravity in the semi- classical approximation for a `new inflation' scenario. We find that spinodal instabilities drive the growth of non-perturbatively large quantum fluctuations which shut off inflation. These fluctuations assemble with the inflaton zero mode yielding a new effective field that actually rolls down behaving classi- cally. We compute the amplitude and index for scalar density and tensor perturbations. In all models of this type the spinodal instabilities produce a `red' spectrum of primordial scalar density perturbations.

While the universe becomes more and more homogeneous at large scales, statistical analysis of galaxy catalogs have revealed a fractal structure at small-scales (\lambda < 100 h^{-1} Mpc), with a fractal dimension D=1.5-2 (Sylos Labini et al 1996). We study the thermodynamics of a self-gravitating system with the theory of critical phenomena and finite-size scaling and show that gravity provides a dynamical mechanism to produce this fractal structure. We develop a field theoretical approach to compute the galaxy distribution, assuming them to be in quasi-isothermal equilibrium. Only a limited, (although large), range of scales is involved, between a short-distance cut-off below which other physics intervene, and a large-distance cut-off, where the thermo- dynamic equilibrium is not satisfied. The galaxy ensemble can be considered at critical conditions, with large density fluctuations developping at any scale. From the theory of critical phenomena, we derive the two independent critical exponents nu and eta and predict the fractal dimension D = 1/nu to be either 1.585 or 2, depending on whether the long-range behaviour is governed by the Ising or the mean field fixed points, respectively. Both set of values are compatible with present observations. In addition, we predict the scaling behaviour of the gravitational potential to be r^{-(1 + eta)/2}. That is, r^{-0.5} for mean field or r^{- 0.519} for the Ising fixed point. The theory allows to compute the three and higher density correlators without any assumption or Ansatz. We find that the N-points density scales as r_1^{(N-1)(D-3)}, when r_1 >> r_i, 2 leq i leq N . There are no free parameters in this theory.

We analyze the dynamics of dissipation and relaxation in the unbroken and broken symmetry phases of scalar theory in the nonlinear regime for large initial energy densities, and after linear unstabilities (parametric or spinodal) are shut-off by the quantum backreaction. A new time scale emerges that separates the linear from the non-linear regimes. This scale is non- perturbative in the coupling and initial amplitude. The non-perturbative evolution is studied in the large the N limit for the O(N) vector model. A combination of numerical analysis and the multitime scale analysis reveals the presence of unstable bands in the nonlinear regime. These are associated with power law growth of quantum fluctuations, that result in power law relaxation and dissipation with non-universal and non-perturbative dynamical anomalous exponents.We find substantial particle production during this non-linear evolu- tion which is of the same order as that in the linear regime. The expectation value of the scalar field vanishes asymptotically transferring all of the initial energy into produced particles via the non-linear resonances in the unbroken symmetry phase.The effective mass for the quantum modes tends asympto- tically to a constant plus oscillating O(1/t) terms. This slow decay causes the power behaviour in the modes which become free for t=infty. We derive a simple expression for the equation of state for the fluid of produced particles that interpolates between radiation-type and dust-type equations according to the initial value of the order parameter for unbroken symmetry. For broken symmetry the produced particles are massless Goldstone bosons with a radiation-type equation of state. We find the onset of a novel form of dynamical Bose condensation in the collisionless regime in the absence of thermalization.

We consider an O(N) model coupled self-consistently to gravity in the semiclassical approximation, where the field is subject to `new inflation' type initial conditions. We study the dynamics self-consistently and non-perturbatively with non-equilibrium field theory methods in the large N limit. We find that spinodal instabilities drive the growth of non-perturbatively large quantum fluctuations which shut off the inflationary growth of the scale factor. We find that a very specific combination of these large fluctuations plus the inflaton zero mode assemble into a new effective field. This new field behaves classically and it is the object which actually rolls down. We show how this reinterpretation saves the standard picture of how metric perturbations are generated during inflation and that the spinodal growth of fluctuations dominates the time dependence of the Bardeen variable for superhorizon modes during inflation. We compute the amplitude and index for the spectrum of scalar density and tensor perturbations and argue that in all models of this type the spinodal instabilities are responsible for a `red' spectrum of primordial scalar density perturbations. A criterion for the validity of these models is provided and contact with the reconstruction program is established validating some of the results within a non-perturbative framework. The decoherence aspects and the quantum to classical transition through inflation are studied in detail by following the full evolution of the density matrix and relating the classicality of cosmological perturbations to that of long-wavelength matter fluctuations.

The physics of the inflationary universe requires the study of the out of equilibrium evolution of quantum fields in curved spacetime. We present the evolution for both the geometry and the matter (described by the quantum inflaton field) by means of the non-perturbative large N limit combined with semi-classical gravitational dynamics including the back-reaction of quantum fluctuations self-consistently for a new inflation scenario. We provide a criterion for the validity of the classical approximation and a full analysis of the case in which spinodal quantum fluctuations drive the evolution of the scale factor. Under carefully determined conditions, we show that the full field equations may be well approximated by those of a single composite field which obeys the classical equation of motion in all cases. The de Sitter stage is found to be followed by a matter dominated phase. We compute the spectrum of scalar density perturbations and argue that the spinodal instabilities are responsible for a `red' spectrum with more power at longer wavelengths. A criterion for the validity of these models is provided and contact with the reconstruction program is established.

We study the non-linear dynamics of quantum fields in matter and radiation dominated universes, using the non-equilibrium field theory approach combined with the non-perturbative Hartree and the large N approximations. We examine the phenomenon of explosive particle production due to spinodal instabilities and parametric amplification in expanding universes with and without symmetry breaking. For a variety of initial conditions, we compute the evolution of the inflaton, its quantum fluctuations, and the equation of state. We find explosive growth of quantum fluctuations, although particle production is somewhat sensitive to the expansion of the universe. In the large N limit for symmetry breaking scenarios, we determine generic late time solutions for any flat Friedman-Robertson-Walker cosmology. We also present a complete and numerically implementable renormalization scheme for the equation of motion and the energy momentum tensor in flat FRW cosmologies. In this scheme the renormalization constants are independent of time and of the initial conditions.

At the end of the inflationary stage of the early universe, profuse particle production leads to the reheating of the universe. Such explosive particle production is due to parametric amplification of quantum fluctuations for the unbroken symmetry case (appropriate for chaotic inflation), or spinodal instabilities in the broken symmetry phase(which is the case in new inflation). The self-consistent methods presented in these lectures are the only approaches, so far, that lead to reliable quantitative results on the reheating mechanism in the inflationary universe. They fully use the field theoretical Schwinger-Keldysh out of equilibrium method in a non-perturbative way. These approaches take into account the non-linear interaction between the quantum modes and exactly conserve energy (covariantly). Simplified analysis that do not include the full backreaction and do not conserve energy, result in unbound particle production and lead to quantitatively erroneous results. For spontaneously broken theories the issue of whether the symmetry may be restored or not by the quantum fluctuations is analyzed. The precise criterion for symmetry restoration is presented. The field dynamics is symmetric when the energy density in the initial state is larger than the top of the tree level potential. When the initial energy density is below the top of the tree level potential, the symmetry is broken.

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.

We analyze the phenomenon of preheating,i.e. explosive particle production due to parametric amplification of quantum fluctuations in the unbroken case, or spinodal instabilities in the broken phase, using the Minkowski space $O(N)$ vector model in the large $N$ limit to study the non-perturbative issues involved. We give analytic results for weak couplings and times short compared to the time at which the fluctuations become of the same order as the tree level,as well as numerical results including the full backreaction.In the case where the symmetry is unbroken, the analytic results agree spectacularly well with the numerical ones in their common domain of validity. In the broken symmetry case, slow roll initial conditions from the unstable minimum at the origin, give rise to a new and unexpected phenomenon: the dynamical relaxation of the vacuum energy.That is, particles are abundantly produced at the expense of the quantum vacuum energy while the zero mode comes back to almost its initial value.In both cases we obtain analytically and numerically the equation of state which turns to be written in terms of an effective polytropic index that interpolates between vacuum and radiation-like domination. We find that simplified analysis based on harmonic behavior of the zero mode, giving rise to a Mathieu equation forthe non-zero modes miss important physics. Furthermore, analysis that do not include the full backreaction do not conserve energy, resulting in unbound particle production. Our results do not support the recent claim of symmetry restoration by non-equilibrium fluctuations.Finally estimates of the reheating temperature are given,as well as a discussion of the inconsistency of a kinetic approach to thermalization when a non-perturbatively large number of particles is created.

The nonlinear time evolution of the quantum fields is studied in the O(N) model for large N in a radiation dominated FRW universe, with a view towards the phenomenon of explosive particle production due to either spinodal instabilities or parametric amplification, i.e. preheating. Quantum backreaction effects due to the produced particles are included consistently within the large N approximation. We find that preheating persists when the expansion is included, although the amount of particle production is reduced compared to the values found in Minkowski space. We also see that the behavior of the evolving zero mode is very different from that in Minkowski space, though the late time behavior in all cases is determined by a sum rule that implies the existence of Goldstone bosons in the final state.

Recent progress in the preheating phenomena for inflationary cosmology is reviewed. We first discuss estimates of the preheating time scale and particle production at the early stages of parametric amplification within the Mathieu and Lam'e approximations and we analyze their precision and limitations. The necessity of self-consistent calculations including the non-linearity of the field theory equations in an energy conserving scheme is stressed. The large N calculations including the field back-reaction are reviewed. For spontaneously broken theories the issue of symmetry restoration is analyzed. A discussion of the possibility and criterion for symmetry restoration is presented.

The gas clouds of the interstellar medium have a fractal structure, the origin of which has generally been thought to lie in turbulence. The energy of turbulence could come from galactic rotation at large-scale, then cascade down to be dissipated on small-scales by viscosity; it has been suggested that such turbulence helps to prevent massive molecular clouds from collapsing in response to their own gravity. Here we show that, on the contrary, self-gravity itself may be the dominant factor in making clouds fractal. We develop a field-theory approach to the structure of clouds, assuming them to be isothermal, and with only gravitational interactions; we find that the observed fractal dimension of the clouds arise naturally from this approach. Although this result does not imply that turbulence is not important, it does demonstrate that the fractal structure can be understood without it.

Using the O(4) linear $\sigma$ model, we address the topic of non-equilibrium relaxation of an inhomogeneous initial configuration due to quantum and thermal fluctuations. The space-time evolution of an inhomogeneous fluctuation of the condensate in the isoscalar channel decaying via the emission of pions in the medium is studied within the context of disoriented chiral condensates. We use out of equilibrium closed time path methods in field theory combined with the amplitude expansion. We give explicit expressions for the asymptotic space-time evolution of an initial inhomogeneous configuration including the contribution of thresholds at zero and non-zero temperature. At non-zero temperature we find new relaxational processes due to thermal cuts that have no counterpart in the homogeneous case. Within the one-loop approximation, we find that the space time evolution of such inhomogeneous configuration out of equilibrium is effectively described in terms of a rapidity dependent temperature $T(\vartheta)=T/\cosh[\vartheta]$ as well as a rapidity dependent decay rate $\Gamma(\vartheta, T(\vartheta))$. This rate is to be interpreted as the production minus absorption rate of pions in the medium and approaches the zero temperature value at large rapidities. An initial configuration localized on a bounded region spreads and decays in spherical waves with slower relaxational dynamics at large rapidity.

Recent progress on string theory in curved spacetimes is reviewed. The string dynamics in cosmological and black hole spacetimes is investigated.The methods to solve the string equations of motion in curved spacetimes are described.That is, the perturbation approach, the null approach, the $\tau$-expansion, and the construction of global solutions.The behaviour of strings in FRW and inflatio- nary spacetimes is obtained from the various types of string solutions. Three different types of behaviour appear:{\bf unstable, dual} to unstable and {\bf stable}.For the unstable strings, the energy and size grow proportional to $R$ for large scale factors $R \to \infty$. For the dual to unstable strings, the energy and size blow up for R\to 0 as 1/R. For stable strings, the energy and size are bounded. (In Minkowski, all solutions are of the stable type). The self-consistent solution to the Einstein equations for string dominated universes is reviewed. The energy-momentum tensor for a gas of strings is taken as source and from the above behaviours the string equation of state is {\bf derived}. The self-consistent string solution exhibits realistic matter dominated behaviour for large times and radiation dominated behaviour for early times. We report on the {\bf exact integrability} of the string dynamics in de Sitter spacetime that allows to systematically find {\bf exact} string solutions by soliton methods. {\bf Multistring solutions} are a new feature in curved spacetimes. That is, a single world-sheet simultaneously describes many different and independent strings. This phenomenon has no analogue in flat spacetime and follows from the coupling of the strings with the geometry. Finally, the string dynamics next and inside a Schwarzschild black hole is analyzed and their physical properties discussed.

We consider the non-equilibrium evolution of the inflaton field coupled to both lighter scalars and fermions. The dissipational dynamics of this field is studied and found to be quite different than that believed in inflationary models. In particular, the damping time scale for the expectation value of the zero momentum mode of the inflaton can be much shorter than that given by the single particle decay rate when the inflaton amplitudes are large, as in chaotic inflation scenarios. We find that the reheating temperature may depart considerably from the usual estimates.

The exact general solution of circular strings in $2+1$ dimensional de Sitter spacetime is described completely in terms of elliptic functions. The novel feature here is that one single world-sheet generically describes {\it infinitely many} (different and independent) strings. This has no analogue in flat spacetime. The circular strings are either oscillating ("stable") or indefinitely expanding ("unstable"). We then compute the {\it exact} equation of state of circular strings in the $2+1$ dimensional de Sitter (dS) and anti de Sitter (AdS) spacetimes, and analyze its properties for the different (oscillating, contracting and expanding) strings. We finally perform a semi-classical quantization of the oscillating circular strings. We find the mass formula $\alpha'm^2_{\mbox{dS}}\approx 4n-5H^2\alpha'n^2,\;(n\in N_0),$ and a {\it finite} number of states $N_{\mbox{dS}}\approx 0.34/(H^2\alpha')$ in de Sitter spacetime; $m^2_{\mbox{AdS}}\approx H^2n^2$ (large $n\in N$) and $N_{\mbox{AdS}}=\infty$ in anti de Sitter spacetime. The level spacing grows with $n$ in AdS spacetime, while is approximately constant (although smaller than in Minkowski spacetime and slightly decreasing) in dS spacetime.

We find exact solutions of the string equations of motion and constraints describing the {\em classical}\ splitting of a string into two. We show that for the same Cauchy data, the strings that split have {\bf smaller} action than the string without splitting. This phenomenon is already present in flat space-time. The mass, energy and momentum carried out by the strings are computed. We show that the splitting solution describes a natural decay process of one string of mass $M$ into two strings with a smaller total mass and some kinetic energy. The standard non-splitting solution is contained as a particular case. We also describe the splitting of a closed string in the background of a singular gravitational plane wave, and show how the presence of the strong gravitational field increases (and amplifies by an overall factor) the negative difference between the action of the splitting and non-splitting solutions.

We compute the {\it exact} equation of state of circular strings in the (2+1) dimensional de Sitter (dS) and anti de Sitter (AdS) spacetimes, and analyze its properties for the different (oscillating, contracting and expanding) strings. The string equation of state has the perfect fluid form $P=(\gamma-1)E,$ with the pressure and energy expressed closely and completely in terms of elliptic functions, the instantaneous coefficient $\gamma$ depending on the elliptic modulus. We semi-classically quantize the oscillating circular strings. The string mass is $m=\sqrt{C}/(\pi H\alpha'),\;C$ being the Casimir operator, $C=-L_{\mu\nu}L^{\mu\nu},$ of the $O(3,1)$-dS [$O(2,2)$-AdS] group, and $H$ is the Hubble constant. We find $\alpha'm^2_{\mbox{dS}}\approx 5.9n,\;(n\in N_0),$ and a {\it finite} number of states $N_{\mbox{dS}}\approx 0.17/(H^2\alpha')$ in de Sitter spacetime; $m^2_{\mbox{AdS}}\approx 4H^2n^2$ (large $n\in N_0$) and $N_{\mbox{AdS}}=\infty$ in anti de Sitter spacetime. The level spacing grows with $n$ in AdS spacetime, while is approximately constant (although larger than in Minkowski spacetime) in dS spacetime. The massive states in dS spacetime decay through tunnel effect and the semi-classical decay probability is computed. The semi-classical quantization of {\it exact} (circular) strings and the canonical quantization of generic string perturbations around the string center of mass strongly agree.

The non-equilibrium dynamics of the first stage of the reheating process, that is dissipation via particle production is studied in scalar field theories in the unbroken and in the broken symmetry phase. We begin with a perturbative study to one loop and show explicitly that the mechanism of dissipation via particle production cannot be explained with a simple derivative term in the equation of motion. The dissipative contribution is non-local and there does not exist a local (Markovian) limit at zero temperature. Furthermore, we show that both an amplitude as well as a one-loop calculation present instabilities, requiring a non-perturbative resummation. Within the same approximations, we study an O(2) linear sigma model that allows to study dissipation by Goldstone bosons. We find infrared divergences that require non-perturbative resummation in order to understand the long-time dynamics. We obtain a perturbative Langevin equation that exhibits a generalized fluctuation-dissipation relation, with non-Markovian kernels and colored noise. We then study a Hartree approximation and clearly exhibit dissipative effects related to the thresholds to particle production. The asymptotic dynamics depends on the coupling and initial conditions but does not seem to lead to exponential relaxation. The effect of dissipation by Goldstone bosons is studied non-perturbatively in the large N limit in an O(N) theory. Dissipation produced by Goldstone bosons dramatically changes the picture of the phase transition. We find the remarkable result that for ``slow-roll'' initial conditions (with the expectation value of the field initially near the origin) the final value of the expectation value of the scalar field is very close to its initial value. We find that the minima of the effective action depend on the initial conditions. We provide extensive numerical analysis of the dynamics.

We compute the string energy-momentum tensor and {\bf derive} the string equation of state from exact string dynamics in cosmological spacetimes. $1+1,~2+1$ and $D$-dimensional universes are treated for any expansion factor $R$. Strings obey the perfect fluid relation $ p = (\gamma -1) \rho $ with three different behaviours: (i) {\it Unstable} for $ R \to \infty $ with growing energy density $ \rho \sim R^{2-D} $, {\bf negative} pressure, and $ \gamma =(D-2)/(D-1) $; (ii){\it Dual} for $ R \to 0 $, with $ \rho \sim R^{-D} $, {\bf positive} pressure and $\gamma = D/(D-1) $ (as radiation); (iii) {\it Stable} for $ R \to \infty $ with $ \rho \sim R^{1-D} $, {\bf vanishing} pressure and $\gamma = 1 $ (as cold matter). We find the back reaction effect of these strings on the spacetime and we take into account the quantum string decay through string splitting. This is achieved by considering {\bf self-consistently} the strings as matter sources for the Einstein equations, as well as for the complete effective string equations. String splitting exponentially suppress the density of unstable strings for large $R$. The self-consistent solution to the Einstein equations for string dominated universes exhibits the realistic matter dominated behaviour $ R \sim (X^0)^{2/(D-1)}\; $ for large times and the radiation dominated behaviour $ R \sim (X^0)^{2/D}\; $ for early times. De Sitter universe does not emerge as solution of the effective string equations. The effective string action (whatever be the dilaton, its potential and the central charge term) is not the appropriate framework in which to address the question of string driven inflation.

We derive the effective equations for the out of equilibrium time evolution of the order parameter and the fluctuations of a scalar field theory in spatially flat FRW cosmologies. After setting the problem in general we propose a non-perturbative, self-consistent Hartree approximation. The method consists of evolving an initial functional thermal density matrix in time and is suitable for studying phase transitions out of equilibrium. The renormalization aspects are studied in detail and we find that the counterterms depend on the initial state. We investigate the high temperature expansion and show that it breaks down at long times. The infinite time limit is computed for de Sitter spacetime. We obtain the time evolution of the initial Boltzmann distribution functions, and argue that in the Hartree approximation, the time evolved state is a ``squeezed'' state. We illustrate the departure from thermal equilibrium by studying the case of a free field in de Sitter and radiation dominated cosmologies. It is found that a suitably defined non-equilibrium entropy per mode increases linearly with comoving time in a de Sitter cosmology, whereas it is {\it not} a monotonically increasing function in the radiation dominated case. This work has been done in collaboration with D. Boyanovsky and R. Holman.

We derive the effective equations for the out of equilibrium time evolution of the order parameter and the fluctuations of a scalar field theory in spatially flat FRW cosmologies.The calculation is performed both to one-loop and in a non-perturbative, self-consistent Hartree approximation.The method consists of evolving an initial functional thermal density matrix in time and is suitable for studying phase transitions out of equilibrium. The renormalization aspects are studied in detail and we find that the counterterms depend on the initial state. We investigate the high temperature expansion and show that it breaks down at long times. We also obtain the time evolution of the initial Boltzmann distribution functions, and argue that to one-loop order or in the Hartree approximation, the time evolved state is a ``squeezed'' state. We illustrate the departure from thermal equilibrium by numerically studying the case of a free massive scalar field in de Sitter and radiation dominated cosmologies. It is found that a suitably defined non-equilibrium entropy per mode increases linearly with comoving time in a de Sitter cosmology, whereas it is {\it not} a monotonically increasing function in the radiation dominated case.

{\bf Exact} solutions of the string equations of motion and constraints are {\bf systematically} constructed in de Sitter spacetime using the dressing method of soliton theory. The string dynamics in de Sitter spacetime is integrable due to the associated linear system. We start from an exact string solution $q_{(0)}$ and the associated solution of the linear system $\Psi^{(0)} (\lambda)$, and we construct a new solution $\Psi(\lambda)$ differing from $\Psi^{(0)}(\lambda)$ by a rational matrix in $\lambda$ with at least four poles $\lambda_{0},1/\lambda_{0},\lambda_{0}^*,1/\lambda_{0}^*$. The periodi- city condition for closed strings restrict $\lambda _{0}$ to discrete values expressed in terms of Pythagorean numbers. Here we explicitly construct solu- tions depending on $(2+1)$-spacetime coordinates, two arbitrary complex numbers (the 'polarization vector') and two integers $(n,m)$ which determine the string windings in the space. The solutions are depicted in the hyperboloid coor- dinates $q$ and in comoving coordinates with the cosmic time $T$. Despite of the fact that we have a single world sheet, our solutions describe {\bf multi- ple}(here five) different and independent strings; the world sheet time $\tau$ turns to be a multivalued function of $T$.(This has no analogue in flat space- time).One string is stable (its proper size tends to a constant for $T\to\infty $, and its comoving size contracts); the other strings are unstable (their proper sizes blow up for $T\to\infty$, while their comoving sizes tend to cons- tants). These solutions (even the stable strings) do not oscillate in time. The interpretation of these solutions and their dynamics in terms of the sinh- Gordon model is particularly enlighting.

The string equations of motion and constraints are solved for a ring shaped Ansatz in cosmological and black hole spacetimes. In FRW universes with arbitrary power behavior [$R(X^0) = a\;|X^0|^{\a}\, $], the asymptotic form of the solution is found for both $X^0 \to 0$ and $X^0 \to \infty$ and we plot the numerical solution for all times. Right after the big bang ($X^0 = 0$), the string energy decreasess as $ R(X^0)^{-1} $ and the string size grows as $ R(X^0) $ for $ 0 < \a < 1 $ and as $ X^0 $ for $ \a < 0 $ and $ \a > 1 $. Very soon [ $ X^0 \sim 1 $] , the ring reaches its oscillatory regime with frequency equal to the winding and constant size and energy. This picture holds for all values of $ \a $ including string vacua (for which, asymptotically, $ \a = 1$). In addition, an exact non-oscillatory ring solution is found. For black hole spacetimes (Schwarzschild, Reissner-Nordstr\oo m and stringy), we solve for ring strings moving towards the center. Depending on their initial conditions (essentially the oscillation phase), they are are absorbed or not by Schwarzschild black holes. The phenomenon of particle transmutation is explicitly observed (for rings not swallowed by the hole). An effective horizon is noticed for the rings. Exact and explicit ring solutions inside the horizon(s) are found. They may be interpreted as strings propagating between the different universes described by the full black hole manifold.

Exact and explicit string solutions in de Sitter spacetime are found. (Here, the string equations reduce to a sinh-Gordon model). A new feature without flat spacetime analogy appears: starting with a single world-sheet, several (here two) strings emerge. One string is stable and the other (unstable) grows as the universe grows. Their invariant size and energy either grow as the expansion factor or tend to constant. Moreover, strings can expand (contract) for large (small) universe radius with a different rate than the universe.

Progress on the physics of strings in curved spacetime are comprehensively reviewed.We start by showing through renormalization group arguments that a meaningful quantum theory of gravity must be finite and must include all particle physics.Then, we review classical and quantum string propagation in curved spacetimes.We start by the general expansion method proposed by de Vega and S\'anchez in 1987. The particle transmutation phenomena in asymptotically flat spacetimes are detailed including fermion-boson transitions in supergravity backgrounds.The next chapters review the exactly solvable cases of string propagation: shock waves, singular plane waves, conical spacetimes and de Sitter cosmological spacetime. The calculation of various physical quantities like the string mass and the energy-momentum tensor shows that classical and quantum string propagation in shock-waves and singular plane waves is physically meaningful and full of interesting new phenomena.The important phenomenom of {\bf string stretching} that takes place when strings fall into spacetime singularities and in expanding universes is analyzed.We conclude by reporting on strings in de Sitter spacetime, where the string equations are integrable and reduce to the sinh-Gordon equation and to integrable generalizations of it. Lectures delivered at the ERICE SCHOOL ``STRING QUANTUM GRAVITY AND PHYSICS AT THE PLANCK ENERGY SCALE'', 21-28 June 1992 , to appear in the Proceedings edited by N. S\'anchez, World Scientific.