In the DARKexp framework for collisionless isotropic relaxation of self--gravitating matter, the central object is the differential energy distribution $n(E)$, which takes a maximum--entropy form proportional to $\exp[-\beta(E - \Phi(0))] - 1$, $\Phi(0)$ being the depth of the potential well and $\beta$ the standard Lagrange multiplier. Then the first and quite non--trivial problem consists in the determination of an ergodic phase--space distribution which reproduces this $n(E)$. In this work we present a very extensive and accurate numerical solution of such DARKexp problem for systems with cored mass density and finite size. This solution holds throughout the energy interval $\Phi(0)\le E\le 0$ and is double--valued for a certain interval of $\beta$. The size of the system represents a unique identifier for each member of this solution family and diverges as $\beta$ approaches a specific value. In this limit, the tail of the mass density $\rho(r)$ dies off as $r^{-4}$, while at small radii it always starts off linearly in $r$, that is $\rho(r)-\rho(0)\propto r$.

The DARKexp model represents a novel and promising attempt to solve a long standing problem of statistical mechanics, that of explaining from first principles the quasi--stationary states at the end of the collisionless gravitational collapse. The model, which yields good fits to observation and simulation data on several scales, was originally conceived to provide a theoretical basis for the $1/r$ cusp of the Navarro--Frenk--White profile. In this note we show that it also allows for cored density profiles that, when viewed in three dimensions, in the $r\to0$ limit have the conical shape characteristic of the Burkert profile. It remains to be established whether both cusps and cores, or only one of the two types, are allowed beyond the asymptotic analysis of this work.

We report the results of extended high--resolution numerical integrations of the Vlasov--Poisson equation for the collapse of spherically symmetric WDM halos. For thermal relics with mass $m=1\,$keV/$c^2$, we find collapsed halos with cores of size from $0.1$ to $0.6\,$kpc. The typical core is hollow, with the mass density decreasing towards the core center by almost three orders of magnitude from its maximum near the core radius $r_{\rm c}$. The core is in equilibrium with the diffused part of the halo but far from virialization. These properties are rooted in the conservation of the squared angular momentum and in the original excess, proper of WDM initial conditions, of kinetic energy in the core region. In a sample of more than one hundred simulated collapses, the values of $r_{\rm c}$ and of the core density $\rho_{\rm c}$ are in the range typical of dwarf spheroids, while the maximal circular velocities $V_{\rm max}$ are proper of small disk galaxies. The product $\mu_{\rm c}=\rho_{\rm c}r_{\rm c}$ takes values between $116\,M_\odot/$pc$^2$ and $283\,M_\odot/$pc$^2$, while the surface density $\mu_0$, as determined from a Burkert fit, is roughly three times larger. From these data and data obtained at smaller values of $m$, we extrapolate for one particular halo $\mu_{\rm c}=263(308)\,M_\odot/$pc$^2$ and $\mu_0=754(855)\,M_\odot/$pc$^2$ at $m=2(3.3)\,$keV/$c^2$, to be compared with the observed value $140^{+83}_{-52} M_\odot$/pc$^2$. In view of the many improvements and enhancements available, we conclude that WDM is a viable solution for explaining the presence and the size of cores in low mass galaxies.

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

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.

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.

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

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)

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.

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.

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

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.

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.