We consider the cosmological production of fermionic dark matter during inflation and a post-inflationary radiation dominated era. This fermion \emph{only interacts gravitationally}, has a mass $m$ much smaller than the Hubble scale during inflation (but is otherwise arbitrary) and is in its Bunch-Davies vacuum state during inflation. We focus on superhorizon modes at the end of inflation, and assume instantaneous reheating. We obtain the full energy momentum tensor discussing its renormalization, and show that the contribution from particle production is of the kinetic-fluid form near matter-radiation equality. We find \emph{exactly} the distribution function of produced particles $|B(k)|^2=\frac{1}{2}\Big[1-(1-e^{-\frac{k^2}{2mT_H}})^{1/2}\Big]$ which exhibits an " emergent temperature" $T_H=H_0\sqrt{\Omega_R}\simeq 10^{-36}(\mathrm{eV})$. The abundance of the produced particles $\Omega_{pp}$ is very similar to that of a non-relativistic degree of freedom thermalized at temperature $T_H$, $\Omega_{pp} \propto m \, (m\,T_H)^{3/2}\simeq \big(m/10^8\mathrm{GeV})\big)^{5/2}$ and "cold" equation of state $w(a) \simeq (T_H/m a^2)$, both dominated by superhorizon modes at the end of inflation. We discuss subtle aspects of isocurvature perturbations.

We study the non-adiabatic cosmological production of ultra light dark matter (ULDM) under a minimal set of assumptions: a free ultra light real scalar as a spectator field in its Bunch-Davies vacuum state during inflation and instantaneous reheating into a radiation dominated era. For (ULDM) fields minimally coupled to gravity, non-adiabatic particle production yields a \emph{distribution function} peaked at \emph{low} comoving momentum $\mathcal{N}_k \propto 1/k^3$. The infrared behavior is a remnant of the infrared enhancement of light minimally coupled fields during inflation. We obtain the full energy momentum tensor, show explicity its equivalence with the fluid-kinetic one in the adiabatic regime, and extract the abundance, equation of state and free streaming length (cutoff in the matter power spectrum). Taking the upper bound on the scale of inflation from Planck, the (UDLM) saturates the dark matter abundance for $m \simeq 1.5\,\times 10^{-5}\mathrm{eV}$ with an equation of state parameter $w \simeq 10^{-14}$ and a free streaming length $\lambda_{fs} \simeq 70\,\mathrm{pc}$. Thus this cosmologically produced (ULDM) yields a \emph{cold} dark matter particle. We argue that the abundance from non-adiabatic production yields a \emph{lower bound} on generic (ULDM) and axion-like particles that must be included in any assessment of (ULDM) as a dark matter candidate.

We study the decay of a Higgs-like scalar Yukawa coupled to massless fermions in post-inflationary cosmology, combining a non-perturbative method with an adiabatic expansion. The renormalized survival probability $\mathcal{P}_\Phi(t)$ of a (quasi) particle ``born'' at time $t_b$ and decaying at rest in the comoving frame, $\mathcal{P}_\Phi(t) = \Big[\frac{t}{t_b}\Big]^{-\frac{Y^2}{8\pi^2}}~ e^{ \frac{Y^2}{4\pi^2}\,\big(t/t_b\big)^{1/4} } \,e^{-\Gamma_0\,(t-t_b)}~ \mathcal{P}_\Phi(t_b) $, with $\Gamma_0$ the decay rate at rest in Minkowski space-time. For an ultrarelativistic particle we find $\mathcal{P}_\Phi(t) = e^{-\frac{2}{3}\Gamma_0\,t_{nr}\,(t/t_{nr})^{3/2}}~ \mathcal{P}_\Phi(t_b)$ before it becomes non-relativistic at a time $t_{nr}$ as a consequence of the cosmological redshift. For $t\gg t_{nr}$ we find $\mathcal{P}_\Phi(t) = \Big[\frac{t}{t_{nr}}\Big]^{-\frac{Y^2}{8\pi^2}}~ e^{ \frac{Y^2}{4\pi^2}\,\big(t/t_{nr}\big)^{1/4} }~\Big[\frac{t}{t_{nr}}\Big]^{\Gamma_0 t_{nr}/2} \,e^{-\Gamma_0\,(t-t_{nr})}~ \mathcal{P}_\Phi(t_{nr})$. The extra power is a consequence of the memory on the past history of the decay process. We compare these results to an S-matrix inspired phenomenological Minkowski-like decay law modified by an instantaneous Lorentz factor to account for cosmological redshift. Such phenomenological description \emph{under estimates the lifetime of the particle}. For very long lived, very weakly coupled particles, we obtain an \emph{upper bound} for the survival probability as a function of redshift $z$ valid throughout the expansion history $\mathcal{P}_\Phi(z) \gtrsim e^{-\frac{\Gamma_0}{H_0}\,\Upsilon(z,z_b)}\,\mathcal{P}_\Phi(z_b)$, where $\Upsilon(z,z_b)$ only depends on cosmological parameters and $t_{nr}$.

We study scalar particle decay during the radiation and matter dominated epochs of a standard cosmological model. An adiabatic approximation is introduced that is valid for degrees of freedom with typical wavelengths much smaller than the particle horizon ($\propto$~Hubble radius) at a given time. We implement a non-perturbative method that includes the cosmological expansion and obtain a cosmological Fermi's Golden Rule that enables one to compute the decay law of a parent particle of mass $m_1$, along with the build up of the population of daughter particles of mass $m_2$. The survival probability of the decaying particle is $P(t)=e^{-\widetilde{\Gamma}_k(t)\,t}$ with $\widetilde{\Gamma}_k(t)$ being an \emph{effective momentum and time dependent decay rate}. It features a transition time scale $t_{nr}$ between the relativistic and non-relativistic regimes and for $k \neq 0$ is always smaller than the analogous rate in Minkowski spacetime, as a consequence of (local) time dilation and the cosmological redshift. For $t \ll t_{nr}$ the decay law is a "stretched exponential" $P(t) = e^{-(t/t^*)^{3/2}}$, whereas for the non-relativistic stage with $t \gg t_{nr}$, we find $P(t) = e^{-\Gamma_0 t}\,(t/t_{nr})^{\Gamma_0\,t_{nr}/2}$. The Hubble time scale $\propto 1/H(t)$ introduces an energy uncertainty $\Delta E \sim H(t)$ which relaxes the constraints of kinematic thresholds. This opens new decay channels into heavier particles for $2\pi E_k(t) H(t) \gg 4m^2_2-m^2_1$, with $E_k(t)$ the (local) comoving energy of the decaying particle. As the expansion proceeds this channel closes and the usual two particle thresholds restrict the decay kinematics.

If the inflaton couples to other degrees of freedom that populate the post-inflationary stage, such coupling modifies the dynamics of the inflaton \emph{during} inflation. We consider light fermions Yukawa coupled to the inflaton as "unobserved" degrees of freedom integrated out of the total density matrix. Tracing out these degrees of freedom yields a \emph{mixed} density matrix whose time evolution is described by an effective field theory. We show that the coupling leads to profuse fermion pair production for super-Hubble inflaton fluctuations which lead to the \emph{growth of entanglement entropy during inflation}. The power spectrum of inflaton fluctuations features scale invariance violations $\mathcal{P}(k) = \mathcal{P}_0(k)\,\,\exp\{8\,\xi_k\}$ with corrections to the \emph{index and its running directly correlated with the entanglement entropy}: $S_{vN} = - \sum_{k} \Big[ \ln(1-\xi_k) + \frac{\xi_k\,\ln(\xi_k)}{1-\xi_k} \Big]$. For super-Hubble fluctuations we find $\xi_k = -\frac{Y^2}{48\pi^2}\Big\{2\,N_T\,\ln(k/k_f) + \ln^2(k/k_f)\Big\}$ with $Y$ the Yukawa coupling, $N_T$ the total number of e-folds during inflation, and $k_f$ a "pivot" scale corresponding to the mode that crosses the Hubble radius at the end of inflation.

In the standard model extended with a seesaw mass matrix, we study the production of sterile neutrinos from the decay of vector bosons at $T\simeq M_{W,Z}$. We derive a general quantum kinetic equation for the production of sterile neutrinos and their effective mixing angles valid in a wide range of temperature, to all orders in interactions of the standard model, and to leading order mixing angle $\theta \ll 1$. Production rates and effective mixing angles depend sensitively on helicity. Positive helicity states interact more weakly with the medium and their effective mixing angle is not modified significantly whereas the mixing angle for negative helicity is strongly suppressed by the medium. If $M_s \lesssim 8.35\,\mathrm{MeV}$, there are fewer states with negative helicity produced than those with positive helicity. There is an MSW resonance in the absence of lepton asymmetry, but is screened by the damping rate, without production enhancement. Negative helicity states freeze-out at $T^-_f\simeq 5\,\mathrm{GeV}$ and positive helicity states freeze-out at $T^+_f \simeq 8\,\mathrm{GeV}$, both distributions are far from thermal. Negative helicity states feature a broader momentum distribution than that for those with positive helicity. Sterile neutrinos produced via vector boson decay do not satisfy the abundance, lifetime and cosmological constraints to be the sole dark matter component in the universe but might solve the $^{7}Li$ problem, albeit at the very edge of the possible parameter space. A heavy sterile neutrino with a mass of a few MeV could decay into light sterile neutrinos, of a few keV in mass, that contribute to warm dark matter. We argue that heavy sterile neutrinos with lifetime $\leq 1/H_0$ reach local thermodynamic equilibrium.

Sterile neutrinos are SU(2) singlets that mix with active neutrinos via a mass matrix, its diagonalization leads to mass eigenstates that couple via standard model vertices. We study the cosmological production of heavy neutrinos via standard model charged and neutral current vertices under a minimal set of assumptions: i) the mass basis contains a hierarchy of heavy neutrinos, ii) these have very small mixing angles with the active (flavor) neutrinos, iii) standard model particles, including light (active-like) neutrinos are in thermal equilibrium. The same weak interaction processes that produce active-like neutrinos also produce the heavier species. We introduce the kinetic equations that describe their production, freeze out and decay and discuss the various processes that lead to their production in a wide range of temperatures assessing their feasibility as dark matter candidates. We identify processes in which finite temperature collective excitations lead to the production of the heavy species. As a specific example, we consider the production of heavy neutrinos from pion decay shortly after the QCD crossover including finite temperature corrections to the pion form factors and mass. We consider the different decay channels that allow for the production of heavy neutrinos showing that their frozen distribution functions exhibit effects from "kinematic entanglement" and argue for their viability as mixed dark matter candidates. We discuss abundance, phase space density and stability constraints and argue that heavy neutrinos with lifetime $\tau> 1/H_0$ freeze out of local thermal equilibrium, and \emph{conjecture} that those with lifetimes $\tau \ll 1/H_0$ may undergo cascade decay into lighter DM candidates and/or inject non-LTE neutrinos into the cosmic neutrino background. We provide a comparison with non-resonant production via active-sterile mixing.

Motivated by apparent persistent large scale anomalies in the CMB we study the influence of fermionic degrees of freedom on the dynamics of inflaton fluctuations as a possible source of violations of (nearly) scale invariance on cosmological scales. We obtain the non-equilibrium effective action of an inflaton-like scalar field with Yukawa interactions ($Y_{D,M}$) to light \emph{fermionic} degrees of freedom both for Dirac and Majorana fields in de Sitter space-time. The effective action leads to Langevin equations of motion for the fluctuations of the inflaton-like field, with self-energy corrections and a stochastic gaussian noise. We solve the Langevin equation in the super-Hubble limit implementing a dynamical renormalization group resummation. For a nearly massless inflaton its power spectrum of super Hubble fluctuations is \emph{enhanced}, $\mathcal{P}(k;\eta) = (\frac{H}{2\pi})^2\,e^{\gamma_t[-k\eta] }$ with $\gamma_t[-k\eta] = \frac{1}{6\pi^2} \Big[\sum_{i=1}^{N_D}{Y^2_{i,D}}+2\sum_{j=1}^{N_M}{Y^2_{j,M}}\Big]\,\Big\{\ln^2[-k\eta]-2 \ln[-k\eta]\ln[ -k\eta_0] \Big\} $ for $N_D$ Dirac and $N_M$ Majorana fermions, and $\eta_0$ is the renormalization scale at which the inflaton mass vanishes. The full power spectrum is shown to be renormalization group invariant. These corrections to the super-Hubble power spectrum entail a violation of scale invariance as a consequence of the coupling to the fermionic fields. The effective action is argued to be \emph{exact} in a limit of large number of fermionic fields. A cancellation between the enhancement from fermionic degrees of freedom and suppression from light scalar degrees of freedom \emph{conformally coupled to gravity} suggests the possibility of a finely tuned \emph{supersymmetry} among these fields.

We obtain the non-equilibrium effective action of an inflaton like scalar field --the system-- by tracing over sub Hubble degrees of freedom of "environmental" light scalar fields. The effective action is stochastic leading to effective Langevin equations of motion for the fluctuations of the inflaton-like field, with self-energy corrections and stochastic noise correlators that obey a de Sitter space-time analog of a fluctuation dissipation relation. We solve the Langevin equation implementing a dynamical renormalization group resummation of the leading secular terms and obtain the corrections to the power spectrum of super Hubble fluctuations of the inflaton field, $\mathcal{P}(k;\eta) = \mathcal{P}_0(k)\,e^{-\gamma(k;\eta)}$ where $\mathcal{P}_0(k)$ is the nearly scale invariant power spectrum in absence of coupling. $\gamma(k;\eta)>0$ describes the suppression of the power spectrum, it features Sudakov-type double logarithms and entails violations of scale invariance. We also obtain the effective action for the case of a heavy scalar field of mass $ M \gg H$, this case yields a local "Fermi" limit with a very weak self-interaction of the inflaton-like field and dissipative terms that are suppressed by powers of $H/M$. We conjecture on the possibility that the large scale anomalies in the CMB may originate in dissipative processes from inflaton coupling to sub-Hubble degrees of freedom.

We study the power spectrum of super-Hubble fluctuations of an inflaton-like scalar field, the "system", coupled to another scalar field, the "environment" during de Sitter inflation. We obtain the reduced density matrix for the inflaton fluctuations by integrating out the environmental degrees of freedom. These are considered to be massless and conformally coupled to gravity as a \emph{proxy} to describe degrees of freedom that remain sub-Hubble all throughout inflation. The time evolution of the density matrix is described by a quantum master equation, which describes the decay of the vacuum state, the production of particles and correlated pairs and quantum entanglement between super and sub-Hubble degrees of freedom. The quantum master equation provides a non-perturbative resummation of secular terms from self-energy (loop) corrections to the inflaton fluctuations. In the case studied here these are Sudakov-type double logarithms which result in the \emph{decay} of the power spectrum of inflaton fluctuations upon horizon crossing with a concomitant violation of scale invariance. The reduced density matrix and its quantum master equation furnish a powerful non-perturbative framework to study the effective field theory of long wavelength fluctuations by tracing short wavelength degrees of freedom.

We study the production of sterile neutrinos in the early universe from $\pi \rightarrow l \nu_s$ shortly after the QCD phase transition in the absence of a lepton asymmetry while including finite temperature corrections to the $\pi$ mass and decay constant $f_{\pi}$. Sterile neutrinos with masses $\lesssim 1 MeV$ produced via this mechanism freeze-out at $T_f \simeq 10 MeV$ with a distribution function that is highly non-thermal and features a sharp enhancement at low momentum thereby making this species \emph{cold} even for very light masses. Dark matter abundance constraints from the CMB and phase space density constraints from the most dark matter dominated dwarf spheroidal galaxies provide upper and lower bounds respectively on combinations of mass and mixing angles. For $\pi \rightarrow \mu \nu_s$, the bounds lead to a narrow region of compatibility with the latest results from the $3.55 \mathrm{KeV}$ line. The non-thermal distribution function leads to free-streaming lengths (today) in the range of $\sim \mbox{few kpc}$ consistent with the observation of cores in dwarf galaxies. For sterile neutrinos with mass $\lesssim 1 eV$ that are produced by this reaction, the most recent accelerator and astrophysical bounds on $U_{ls}$ combined with the non-thermal distribution function suggests a substantial contribution from these sterile neutrinos to $N_{eff}$.

We study corrections to the power spectra of curvature and tensor perturbations and the tensor-to-scalar ratio in single field slow roll inflation due to initial conditions imprinted by a fast-roll stage prior to slow roll. For a wide range of initial inflaton kinetic energy, this stage lasts only a few e-folds and merges smoothly with slow-roll leading to non-Bunch-Davies initial conditions for modes that exit the Hubble radius during slow roll. We describe a program that yields the dynamics in the fast-roll stage while matching to the slow roll stage independent of the inflationary potentials. Corrections to the power spectra are encoded in a transfer function $\mathcal{T}_\alpha(k)$. We obtain $\mathcal{T}_\alpha(k)$ to leading order in a Born approximation valid for modes of observational relevance today. A fit yields $\mathcal{T}_\alpha(k) =1+ A_{\alpha}k^{-p}\cos[2\pi \omega k/H_{sr}+\varphi_\alpha]$, with $1.5 \lesssim p \lesssim 2$, $\omega \simeq 1$ and $H_{sr}$ the Hubble scale during slow roll inflation, where curvature and tensor perturbations feature the same $p,\omega$ for a wide range of initial conditions. These corrections lead to both a suppression of the quadrupole and oscillatory features in both $P_R(k)$ and $r(k_0)$ with a period of the order of the Hubble scale during slow roll inflation. The results are independent of the specific inflationary potentials, depending solely on the ratio of kinetic to potential energy $\kappa$ and the slow roll parameters to leading order. For a wide range of $\kappa$ and the values of $\epsilon_V;\eta_V$ corresponding to the upper bounds from Planck, we find that the low quadrupole is consistent with the results from Planck, and the oscillations in $r(k_0)$ could be observable if the modes corresponding to the quadrupole and the pivot scale crossed the Hubble radius a few e-folds after the onset of slow roll.

If the large scale anomalies in the temperature power spectrum of the cosmic microwave background are of primordial origin, they may herald modifications to the slow roll inflationary paradigm on the largest scales. We study the possibility that the origin of the large scale power suppression is a modification of initial conditions during slow roll as a result of a pre-slow roll phase during which the inflaton evolves rapidly. This stage is manifest in a potential in the equations for the Gaussian fluctuations during slow roll and modify the power spectra of scalar perturbations via an initial condition transfer function $\mathcal{T}(k)$. We provide a general analytical study of its large and small scale properties and analyze the impact of these initial conditions on the infrared aspects of typical test scalar fields. The infrared behavior of massless minimally coupled test scalar field theories leads to the dynamical generation of mass and anomalous dimensions, both depend non-analytically on $\mathcal{T}(0)$. During inflation all quanta decay into many quanta even of the same field because of the lack of kinematic thresholds. The decay leads to a quantum entangled state of sub and superhorizon quanta with correlations across the horizon. We find the modifications of the decay width and the entanglement entropy from the initial conditions. In all cases, initial conditions from a ``fast-roll'' stage that lead to a suppression in the scalar power spectrum at large scales also result in a suppression of the dynamically generated masses, anomalous dimensions and decay widths.

In inflationary cosmology all particle states decay as a consequence of the lack of kinematic thresholds. The decay of an initial single particle state yields an \emph{entangled quantum state of the product particles}. We generalize and extend a manifestly unitary field theoretical method to obtain the time evolution of the quantum state. We consider the decay of a light scalar field with mass $M\ll H$ with a cubic coupling in de Sitter space-time. Radiative corrections feature an infrared enhancement manifest as poles in $\Delta=M^2/3H^2$ and we obtain the quantum state in an expansion in $\Delta$. To leading order in $\Delta$ the pure state density matrix describing the decay of a particle with sub-horizon wavevector is dominated by the emission of superhorizon quanta, describing \emph{entanglement between superhorizon and subhorizon fluctuations and correlations across the horizon}. Tracing over the superhorizon degrees of freedom yields a mixed state density matrix from which we obtain the entanglement entropy. Asymptotically this entropy grows with the \emph{physical} volume as a consequence of more modes of the decay products crossing the Hubble radius. A generalization to localized wave packets is provided. The cascade decay of single particle states into many particle states is discussed. We conjecture on \emph{possible} impact of these results on non-gaussianity and on the ``low multipole anomalies'' of the CMB.

We argue that in an inflationary cosmology a consequence of the lack of time translational invariance is that spontaneous breaking of a continuous symmetry and Goldstone's theorem \emph{do not} imply the existence of \emph{massless} Goldstone modes. We study spontaneous symmetry breaking in an O(2) model, and implications for O(N) in de Sitter space time. The Goldstone mode acquires a radiatively generated mass as a consequence of infrared divergences, and the continuous symmetry is spontaneously broken for any finite $N$, however there is a \emph{first order phase transition} as a function of the Hawking temperature $T_H=H/2\pi$. For O(2) the symmetry is spontaneously broken for $T_H < T_c= \lambda^{1/4} v/2.419$ where $\lambda$ is the quartic coupling and $v$ is the tree level vacuum expectation value and the Goldstone mode acquires a radiatively generated mass $\mathcal{M}^2_\pi \propto \lambda^{1/4} H$. The first order nature of the transition is a consequence of the strong infrared behavior of minimally coupled scalar fields in de Sitter space time, the jump in the order parameter at $T_H=T_c$ is $\sigma_{0c} \simeq 0.61\, {H}/{\lambda^{1/4}}$. In the strict $N\rightarrow \infty$ the symmetry cannot be spontaneously broken. Furthermore, the lack of kinematic thresholds imply that the Goldstone modes \emph{decay} into Goldstone and Higgs modes by emission and absorption of superhorizon quanta.

During de Sitter inflation massless particles of minimally coupled scalar fields acquire a mass and a decay width thereby becoming \emph{quasiparticles}. For bare massless particles non-perturbative infrared radiative corrections lead to a self-consistent generation of mass, for a quartic self interaction $M \propto \lambda^{1/4} H$, and for a cubic self-interaction the mass is induced by the formation of a non-perturbative \emph{condensate} leading to $M \propto \lambda^{1/3} H^{2/3}$. These radiatively generated masses restore de Sitter invariance and result in anomalous scaling dimensions of superhorizon fluctuations. We introduce a generalization of the non-perturbative Wigner-Weisskopf method to obtain the time evolution of quantum states that include the self-consistent generation of mass and regulate the infrared behavior. The infrared divergences are manifest as poles in $\Delta=M^2/3H^2$ in the single particle self-energies, leading to a re-arrangement of the perturbative series non-analytic in the couplings. A set of simple rules that yield the leading order infrared contributions to the decay width are obtained and implemented. The lack of kinematic thresholds entail that all particle states acquire a decay width, dominated by the emission and absorption of superhorizon quanta $\propto (\lambda/H)^{4/3}\,[H/k_{ph}(\eta)]^6 ; \lambda\,[H/k_{ph}(\eta)]^6 $ for cubic and quartic couplings respectively to leading order in $M/H$. The decay of single particle quantum states hastens as their wavevectors cross the Hubble radius and their width is related to the highly squeezed limit of the bi- or tri-spectrum of scalar fluctuations respectively.

We use a field theoretic generalization of the Wigner-Weisskopf method to study the stability of the Bunch-Davies vacuum state for a massless, conformally coupled interacting test field in de Sitter space. We find that in $\lambda \phi^4$ theory the vacuum does {\em not} decay, while in non-conformally invariant models, the vacuum decays as a consequence of a vacuum wave function renormalization that depends \emph{singularly} on (conformal) time and is proportional to the spatial volume. In a particular regularization scheme the vacuum wave function renormalization is the same as in Minkowski spacetime, but in terms of the \emph{physical volume}, which leads to an interpretation of the decay. A simple example of the impact of vacuum decay upon a non-gaussian correlation is discussed. Single particle excitations also decay into two particle states, leading to particle production that hastens the exiting of modes from the de Sitter horizon resulting in the production of \emph{entangled superhorizon pairs} with a population consistent with unitary evolution. We find a non-perturbative, self-consistent "screening" mechanism that shuts off vacuum decay asymptotically, leading to a stationary vacuum state in a manner not unlike the approach to a fixed point in the space of states.

We study the scale and redshift dependence of the power spectra for density perturbations and peculiar velocities, and the evolution of a coarse grained phase space density for (WDM) particles that decoupled during the radiation dominated stage. The (WDM) corrections are obtained in a perturbative expansion valid in the range of redshifts at which N-body simulations set up initial conditions, and for a wide range of scales. The redshift dependence is determined by the kurtosis $\beta_2$ of the distribution function at decoupling. At large redshift there is an enhancement of peculiar velocities for $\beta_2 > 1$ that contributes to free streaming and leads to further suppression of the matter power spectrum and an enhancement of the peculiar velocity autocorrelation function at scales smaller than the free streaming scale. Statistical fluctuations of peculiar velocities are also suppressed on these scales by the same effect. In the linearized approximation, the coarse grained phase space density features redshift dependent (WDM) corrections from gravitational perturbations determined by the power spectrum of density perturbations and $\beta_2$. For $\beta_2 > 25/21$ it \emph{grows logarithmically} with the scale factor as a consequence of the suppression of statistical fluctuations. Two specific models for WDM are studied in detail. The (WDM) corrections relax the bounds on the mass.

We provide a semi-analytic study of the small scale aspects of the power spectra of warm dark matter (WDM) candidates that decoupled while relativistic with arbitrary distribution functions. These are characterized by two widely different scales $k_{eq} \sim 0.01\,(\mathrm{Mpc})^{-1}$ and $k_{fs}= \sqrt{3}\,k_{eq}/2\,< V^2_{eq} >^{1/2} $ with $< V^2_{eq} >^{1/2} \ll 1 $ the velocity dispersion at matter radiation equality. Density perturbations evolve through three stages: radiation domination when the particle is relativistic and non-relativistic and matter domination. An early ISW effect during the first stage leads to an enhancement of density perturbations and a plateau in the transfer function for $k \lesssim k_{fs}$. An effective fluid description emerges at small scales which includes the effects of free streaming in initial conditions and inhomogeneities. The transfer function features \emph{WDM-acoustic oscillations} at scales $k \gtrsim 2 \,k_{fs}$. We study the power spectra for two models of sterile neutrinos with $m \sim \,\mathrm{keV}$ produced non-resonantly, at the QCD and EW scales respectively. The latter case yields acoustic oscillations on mass scales $\sim 10^{8}\,M_{\odot}$. Our results reveal a \emph{quasi-degeneracy} between the mass, distribution function and decoupling temperature suggesting caveats on the constraints on the mass of a sterile neutrino from current WDM N-body simulations and Lyman-$\alpha$ forest data. A simple analytic interpolation of the power spectra between large and small scales and its numerical implementation is given.

We study the production of sterile neutrinos in the region $T\sim M_W$ in an extension beyond the standard model with the see-saw mass matrix originating in Yukawa couplings to Higgs-like scalars with masses and vev's of the order of the electroweak scale. Sterile neutrinos are produced by the decay of scalars and standard model vector bosons. We obtain the index of refraction, dispersion relations, mixing angles in the medium and production rates including those for right-handed sterile neutrinos, from the standard model and beyond the standard model self-energies. For $1 \lesssim M_W/T \lesssim 3$ we find narrow MSW resonances with $k \lesssim T$ for both left and right handed neutrinos even in absence of a lepton asymmetry in the (active) neutrino sector, as well as very low energy ($k/T \ll |\xi|$) narrow MSW resonances in the presence of a lepton asymmetry consistent with the bounds from WMAP and BBN. For small vacuum mixing angle, consistent with observational bounds, the absorptive part of the self-energies lead to a strong damping regime very near the resonances resulting in the \emph{exact} degeneracy of the propagating modes with a concomitant breakdown of adiabaticity. We argue that cosmological expansion sweeps through the resonances, \emph{resonant and non-resonant} sterile neutrino production results in a highly \emph{non-thermal} distribution function enhanced at small momentum $k < T$, with potentially important consequences for their free streaming length and transfer function at small scales.

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 clustering properties of sterile neutrinos are studied within an extension of the minimal standard model, where these are produced via the decay of a gauge singlet scalar. The distribution function after decoupling is strongly out of equilibrium. (DM) abundance and phase space density constraints from (dSphs) constrain the mass in the $\mathrm{keV}$ range consistent with a gauge singlet with mass and vacuum expectation value $\sim 100,\mathrm{GeV}$ decoupling at this temperature. The (DM) transfer function and power spectrum are obtained from the solution of the non-relativistic Boltzmann-Vlasov equation in the matter dominated era. The small momentum enhancement of the distribution function leads to long range memory of gravitational clustering and a \emph{substantial enhancement of the power spectrum at small scales compared to a thermal relic or sterile neutrino produced via non-resonant mixing with active neutrinos}. The scale of suppression of the power spectrum for such sterile neutrino with $m\sim \mathrm{keV}$ is $\lambda \sim 488 ,\mathrm{kpc}$. At large scales $T(k)\sim 1-C, k^2/k^2_{fs}(t_{eq}) +...$ with $C \sim \mathrm{O}(1)$. At small scales $65 \mathrm{kpc} \lesssim \lambda \lesssim 500 \mathrm{kpc}$ corrections to the fluid description and memory of gravitational clustering become important, and we find $T(k) \simeq 1.902 e^{-k/k_{fs}(t_{eq})}$, where $k_{fs}(t_{eq}) \sim 0.013/\mathrm{kpc}$ is the free streaming wavevector at matter-radiation equality. The enhancement of power at small scales may provide a possible relief to the tension between the constraints from X-ray and Lyman-$\alpha$ forest data.

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.

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

Free streaming in a \emph{mixture} of collisionless non-relativistic dark matter (DM) particles is studied by implementing methods from the theory of multicomponent plasmas. The mixture includes Fermionic, condensed and non condensed Bosonic particles decoupling in equilibrium while relativistic, heavy non-relativistic thermal relics (WIMPs), and sterile neutrinos that decouple \emph{out of equilibrium} when they are relativistic. The free-streaming length $\lambda_{fs}$ is obtained from the marginal zero of the gravitational polarization function, which separates short wavelength Landau-damped from long wavelength Jeans-unstable \emph{collective} modes. At redshift $z$ we find $ \frac{1}{\lambda^2_{fs}(z)}= \frac{1}{(1+z)} \big[\frac{0.071}{\textrm{kpc}} \big]^2 \sum_{a}\nu_a g^{2/3}_{d,a}({m_a}/{\mathrm{keV}})^2 I_a $,where $0\leq \nu_a \leq 1$ are the \emph{fractions} of the respective DM components of mass $m_a$ that decouple when the effective number of ultrarelativistic degrees of freedom is $g_{d,a}$, and $I_a$ only depend on the distribution functions at decoupling, given explicitly in all cases. If sterile neutrinos produced either resonantly or non-resonantly that decouple near the QCD scale are the \emph{only} DM component,we find $\lambda_{fs}(0) \simeq 7 \mathrm{kpc} (\mathrm{keV}/m)$ (non-resonant), $\lambda_{fs}(0) \simeq 1.73 \mathrm{kpc} (\mathrm{keV}/m)$ (resonant).If WIMPs with $m_{wimp} \gtrsim 100 \mathrm{GeV}$ decoupling at $T_d \gtrsim 10 \mathrm{MeV}$ are present in the mixture with $\nu_{wimp} \gg 10^{-12}$,$\lambda_{fs}(0) \lesssim 6.5 \times 10^{-3} \mathrm{pc}$ is \emph{dominated} by CDM. If a Bose Einstein condensate is a DM component its free streaming length is consistent with CDM because of the infrared enhancement of the distribution function.

We study several aspects of the kinetic approach to sterile neutrino production via active-sterile mixing. We obtain the neutrino propagator in the medium including self-energy corrections up to $\mathcal{O}(G^2_F)$, from which we extract the dispersion relations and damping rates of the propagating modes. The dispersion relations are the usual ones in terms of the index of refraction in the medium, and the damping rates are $\Gamma_1(k) = \Gamma_{aa}(k) \cos^2\theta_m(k); \Gamma_2(k) = \Gamma_{aa}(k) \sin^2\theta_m(k)$ where $\Gamma_{aa}(k)\propto G^2_F k T^4$ is the active neutrino scattering rate and $\theta_m(k)$ is the mixing angle in the medium. We provide a generalization of the transition probability in the \emph{medium from expectation values in the density matrix}: $ P_{a\to s}(t) = \frac{\sin^22\theta_m}{4}[e^{-\Gamma_1t} + e^{-\Gamma_2 t}-2e^{-{1/2}(\Gamma_1+\Gamma_2)t} \cos\big(\Delta E t\big)] $ and study the conditions for its quantum Zeno suppression directly in real time. We find the general conditions for quantum Zeno suppression, which for $m_s\sim \textrm{keV}$ sterile neutrinos with $\sin2\theta \lesssim 10^{-3}$ \emph{may only be} fulfilled near an MSW resonance. We discuss the implications for sterile neutrino production and argue that in the early Universe the wide separation of relaxation scales far away from MSW resonances suggests the breakdown of the current kinetic approach.

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

Charged lepton mixing as a consequence of neutrino mixing is studied for two generations $e,\mu$ in the temperature regime $m_\mu \ll T \ll M_W$ in the early Universe. We state the general criteria for charged lepton mixing, critically reexamine aspects of neutrino equilibration and provide arguments to suggest that neutrinos may equilibrate as mass eigenstates in the temperature regime \emph{prior} to flavor equalization. We assume this to be the case, and that neutrino mass eigenstates are in equilibrium with different chemical potentials. Charged lepton self-energies are obtained to leading order in the electromagnetic and weak interactions. The upper bounds on the neutrino asymmetry parameters from CMB and BBN without oscillations, combined with the fit to the solar and KamLAND data for the neutrino mixing angle, suggest that for the two generation case there is resonant \emph{charged lepton} mixing in the temperature range $T \sim 5 \mathrm{GeV}$. In this range the charged lepton oscillation frequency is of the same order as the electromagnetic damping rate.

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 study the space-time evolution of ``flavor'' neutrino wave-packets at finite temperature and density in the early Universe prior to BBN. We implement non-equilibrium field theory methods and linear response to study the space-time evolution directly from the effective Dirac equation in the medium. There is a rich hierarchy of time scales associated with transverse and longitudinal dispersion and coherence. A phenomenon of ``freezing of coherence'' is a result of a competition between longitudinal dispersion and the separation of wave-packets of propagating modes in the medium. Near a resonance the coherence and oscillation time scales are enhanced by a factor $1/\sin2\theta$ compared to the vacuum. Collisional relaxation via charged and neutral currents occurs on time scales much shorter than the coherence time scale and for small vacuum mixing angle, shorter than the oscillation scale. Assuming that the momentum spread of the initial wave packet is determined by the large angle scattering mean free path of charged leptons, we find that the transverse dispersion time scale is the shortest and is responsible for a large suppression in both the survival and transition probabilities on time scales much shorter than the Hubble time. For small mixing angle the oscillation time scale is \emph{longer} than the collisional relaxation scale. The method also yields the evolution of right-handed wave packets. Corrections to the oscillation frequencies emerge from wave-packet structure as well as from the energy dependence of mixing angles in the medium.

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.

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.

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.

Neutrino collective excitations are studied in the Standard Model at high temperatures below the symmetry breaking scale. Two parameters determine the properties of the collective excitations: a mass scale $m_\nu=gT/4$ which determines the \emph{chirally symmetric} gaps in the spectrum and $\Delta=M^2_W(T)/2m_\nu T$. The spectrum consists of left handed negative helicity quasiparticles, left handed positive helicity quasiholes and their respective antiparticles. For $\Delta < \Delta_c = 1.275...$ there are two gapped quasiparticle branches and one gapless and two gapped quasihole branches, all but the higher gapped quasiparticle branches terminate at end points. For $\Delta_c < \Delta < \pi/2$ the quasiparticle spectrum features a pitchfork bifurcation and for $\Delta >\pi/2$ the collective modes are gapless quasiparticles with dispersion relation below the light cone for $k\ll m_\nu$ approaching the free field limit for $k\gg m_\nu$ with a rapid crossover between the soft non-perturbative to the hard perturbative regimes for $k\sim m_\nu$.The \emph{decay} of the vector bosons leads to a \emph{width} of the collective excitations of order $g^2$ which is explicitly obtained in the limits $k =0$ and $k\gg m_\nu \Delta$. At high temperature this damping rate is shown to be competitive with or larger than the collisional damping rate of order $G^2_F$ for a wide range of neutrino energy.

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

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.

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.

In these lectures I discuss cosmological phase transitions with the goal of establishing the possibility of observational consequences. I argue that the {\em only} phase transition amenable of experimental study within the foreseeable future is that predicted by QCD and discuss some of the potential observational cosmological consequences associated with this phase transition(s). I describe the experimental effort to study the QCD phase transition(s) at RHIC and SPS and summarize some of the recent experimental results. The possibility of novel phases of QCD in the core of pulsars is discussed along with the suggested observational consequences. A brief review of standard big bang cosmology as well as the astrophysics of compact stars sets the stage for understanding the observational cosmological and astrophysical consequences of phase transitions in the standard model.

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.

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.

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.

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.

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

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.

We solve the equations of motion for a scalar field with domain wall boundary conditions in a Friedmann-Robertson-Walker (FRW) spacetime. We find (in agreement with Basu and Vilenkin) that no domain wall solutions exist in de Sitter spacetime for h = H/m >= 1/2, where H is the Hubble parameter and m is the scalar mass. In the general FRW case we develop a systematic perturbative expansion in h to arrive at an approximate solution to the field equations. We calculate the energy momentum tensor of the domain wall configuration, and show that the energy density can become negative at the core of the defect for some values of the non-minimal coupling parameter xi. We develop a translationally invariant theory for fluctuations of the wall, obtain the effective Lagrangian for these fluctuations, and quantize them using the Bunch-Davies vacuum in the de Sitter case. Unlike previous analyses, we find that the fluctuations act as zero-mass (as opposed to tachyonic) modes. This allows us to calculate the distortion and the normal-normal correlators for the surface. The normal-normal correlator decreases logarithmically with the distance between points for large times and distances, indicating that the interface becomes rougher than in Minkowski spacetime.

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