Cosmological phase transitions in the early Universe may produce relics in the form of a network of cosmic defects. Independently of the order of a phase transition, topology of the defects, and their global or gauge nature, the defects are expected to emit gravitational waves (GWs) as the network energy-momentum tensor adapts itself to maintaining {scaling}. We show that the evolution of any defect network (and for that matter any scaling source) emits a GW background with spectrum $\Omega_{\rm GW} \propto f^3$ for $f \ll f_0$, $\Omega_{\rm GW} \propto 1/f^2$ for $f_0 \lesssim f \lesssim f_{\rm eq}$, and $\Omega_{\rm GW} \propto~const$ (i.e.~exactly scale-invariant) for $f \gg f_{\rm eq}$, where $f_0$ and $ f_{\rm eq}$ denote respectively the frequencies corresponding to the present and matter-radiation equality horizons. This background represents an irreducible emission of GWs from any scaling network of cosmic defects, with its amplitude characterized only by the symmetry breaking scale and the nature of the defects. Using classical lattice simulations we calculate the GW signal emitted by defects created after the breaking of a global symmetry $O(N) \rightarrow O(N-1)$. We obtain the GW spectrum for $N$ between 2 and 20 with two different techniques: integrating over unequal time correlators of the energy momentum tensor, updating our previous work on smaller lattices, and for the first time, comparing the result with the real time evolution of the tensor perturbations sourced by the same defects. Our results validate the equivalence of the two techniques. Using CMB upper bounds on the defects' energy scale, we discuss the difficulty of detecting this GW background in the case of global defects.

It is commonly assumed that the energy density of the Universe was dominated by radiation between reheating after inflation and the onset of matter domination 54,000 years later. While the abundance of light elements indicates that the Universe was radiation dominated during Big Bang Nucleosynthesis (BBN), there is scant evidence that the Universe was radiation dominated prior to BBN. It is therefore possible that the cosmological history was more complicated, with deviations from the standard radiation domination during the earliest epochs. Indeed, several interesting proposals regarding various topics such as the generation of dark matter, matter-antimatter asymmetry, gravitational waves, primordial black holes, or microhalos during a nonstandard expansion phase have been recently made. In this paper, we review various possible causes and consequences of deviations from radiation domination in the early Universe - taking place either before or after BBN - and the constraints on them, as they have been discussed in the literature during the recent years.

We present a comprehensive discussion on lattice techniques for the simulation of scalar and gauge field dynamics in an expanding universe. After reviewing the continuum formulation of scalar and gauge field interactions in Minkowski and FLRW backgrounds, we introduce basic tools for the discretization of field theories, including lattice gauge invariant techniques. Following, we discuss and classify numerical algorithms, ranging from methods of $O(dt^2)$ accuracy like $staggered~leapfrog$ and $Verlet$ integration, to $Runge-Kutta$ methods up to $O(dt^4)$ accuracy, and the $Yoshida$ and $Gauss-Legendre$ higher-order integrators, accurate up to $O(dt^{10})$. We adapt these methods for their use in classical lattice simulations of the non-linear dynamics of scalar and gauge fields in an expanding grid in $3+1$ dimensions, including the case of `self-consistent' expansion sourced by the volume average of the fields' energy and pressure densities. We present lattice formulations of canonical cases of: $i)$ Interacting scalar fields, $ii)$ Abelian $U(1)$ gauge theories, and $iii)$ Non-Abelian $SU(2)$ gauge theories. In all three cases we provide symplectic integrators, with accuracy ranging from $O(dt^2)$ up to $O(dt^{10})$. For each algorithm we provide the form of relevant observables, such as energy density components, field spectra and the Hubble constraint. We note that all our algorithms for gauge theories always respect the Gauss constraint to machine precision, even in the case of 'self-consistent' expansion. As a numerical example we analyze the post-inflationary dynamics of an oscillating inflaton charged under $SU(2)\times U(1)$. We note that the present manuscript is meant as part of the theoretical basis for the code $CosmoLattice$, a multi-purpose MPI-based package for simulating the non-linear evolution of field theories in an expanding universe, to be released in 2020.

We study the energy distribution and equation of state of the universe between the end of inflation and the onset of radiation domination (RD), considering observationally consistent single-field inflationary scenarios, with a potential 'flattening' at large field values, and a monomial shape $V(\phi) \propto |\phi|^p$ around the origin. We include a quadratic interaction $g^2\phi^2X^2$ between the inflaton $\phi$ and a light scalar 'daughter' field $X$, as a proxy for (p)reheating. We capture the non-perturbative and non-linear nature of the system dynamics with lattice simulations, obtaining that: $i)$ the final energy transferred to $X$ depends only on $p$, not on $g^2$; $ii)$ the final transfer of energy is always negligible for $2 \leq p < 4$, and of order $\sim 50\%$ for $p \geq 4$; $iii)$ the system goes at late times to matter-domination for $p = 2$, and always to RD for $p > 2$. In the latter case we calculate the number of e-folds until RD, significantly reducing the uncertainty in the inflationary observables $n_s$ and $r$.

The stochastic gravitational wave background (SGWB) contains a wealth of information on astrophysical and cosmological processes. A major challenge of upcoming years will be to extract the information contained in this background and to disentangle the contributions of different sources. In this paper we provide the formalism to extract, from the correlation of three signals in the Laser Interferometer Space Antenna (LISA), information about the tensor three-point function, which characterizes the non-Gaussian properties of the SGWB. This observable can be crucial to discriminate whether a SGWB has a primordial or astrophysical origin. Compared to the two-point function, the SGWB three-point function has a richer dependence on the gravitational wave momenta and chiralities. It can be used therefore as a powerful discriminator between different models. For the first time we provide the response functions of LISA to a general SGWB three-point function. As examples, we study in full detail the cases of an equilateral and squeezed SGWB bispectra, and provide the explicit form of the response functions, ready to be convoluted with any theoretical prediction of the bispectrum to obtain the observable signal. We further derive the optimal estimator to compute the signal-to-noise ratio. Our formalism covers general shapes of non-Gaussianity, and can be extended straightaway to other detector geometries. Finally, we provide a short overview of models of the early universe that can give rise to a non-Gaussian SGWB.

From a model-building perspective, the inflationary sector might very well have no direct couplings to other species, apart from inevitable gravitational interactions. Within the context of General Relativity, a thermal universe can still emerge after inflation if: $i)$ some radiation sector is excited towards the end of inflation, and $ii)$ the post-inflationary equation of state becomes sufficiently stiff $w \geq w_{\rm RD}\gtrsim 0.57$, with $w_{\rm RD}$ a threshold depending on the inflationary scale $H_*$ and the initial radiation-to-inflaton energy ratio $\Delta_*$. Furthermore, a stiff period in the expansion history enhances significantly the inflationary gravitational wave (GW) background, making this signal (potentially) observable by aLIGO, LISA and other experiments. The very same enhancement leads however to an inconsistency of the scenario: the energy of the GWs becomes too large compared to the rest of the radiation sector, violating standard BBN and CMB bounds on GW backgrounds. Except for very special scenarios where the initial radiation sector comprises hundreds of fields with couplings tuned to specific values, our result applies independently of $w$, $H_*$ and $\Delta_*$. This suggests that in order to reheat the universe, the inflationary sector should be coupled directly to other particle species. Alternatively the inflationary sector could be implemented in modified gravity theories.

Cosmic string networks offer one of the best prospects for detection of cosmological gravitational waves (GWs). The combined incoherent GW emission of a large number of string loops leads to a stochastic GW background (SGWB), which encodes the properties of the string network. In this paper we analyze the ability of the Laser Interferometer Space Antenna (LISA) to measure this background, considering leading models of the string networks. We find that LISA will be able to probe cosmic strings with tensions $G\mu \gtrsim \mathcal{O}(10^{-17})$, improving by about $6$ orders of magnitude current pulsar timing arrays (PTA) constraints, and potentially $3$ orders of magnitude with respect to expected constraints from next generation PTA observatories. We include in our analysis possible modifications of the SGWB spectrum due to different hypotheses regarding cosmic history and the underlying physics of the string network. These include possible modifications in the SGWB spectrum due to changes in the number of relativistic degrees of freedom in the early Universe, the presence of a non-standard equation of state before the onset of radiation domination, or changes to the network dynamics due to a string inter-commutation probability less than unity. In the event of a detection, LISA's frequency band is well-positioned to probe such cosmic events. Our results constitute a thorough exploration of the cosmic string science that will be accessible to LISA.

The expansion history of the Universe between the end of inflation and the onset of radiation-domination (RD) is currently unknown. If the equation of state during this period is stiffer than that of radiation, $w > 1/3$, the gravitational wave (GW) background from inflation acquires a blue-tilt ${d\log\rho_{\rm GW}\over d\log f} = {2(w-1/3)\over (w+1/3)} > 0$ at frequencies $f \gg f_{\rm RD}$ corresponding to modes re-entering the horizon during the stiff-domination (SD), where $f_{\rm RD}$ is the frequency today of the horizon scale at the SD-to-RD transition. We characterized in detail the transfer function of the GW energy density spectrum, considering both 'instant' and smooth modelings of the SD-to-RD transition. The shape of the spectrum is controlled by $w$, $f_{\rm RD}$, and $H_{\rm inf}$ (the Hubble scale of inflation). We determined the parameter space compatible with a detection of this signal by LIGO and LISA, including possible changes in the number of relativistic degrees of freedom, and the presence of a tensor tilt. Consistency with upper bounds on stochastic GW backgrounds, however, rules out a significant fraction of the observable parameter space. We find that this renders the signal unobservable by Advanced LIGO, in all cases. The GW background remains detectable by LISA, though only in a small island of parameter space, corresponding to scenarios with an equation of state in the range $0.46 \lesssim w \lesssim 0.56$ and a high inflationary scale $H_{\rm inf} \gtrsim 10^{13}~{\rm GeV}$, but low reheating temperature $1~{\rm MeV} \lesssim T_{\rm RD} \lesssim 150~{\rm MeV}$ (equivalently, $10^{-11}~{\rm Hz} \lesssim f_{\rm RD} \lesssim 3.6\cdot10^{-9}~{\rm Hz}$). Implications for early Universe scenarios resting upon an SD epoch are briefly discussed.

We present a set of tools to assess the capabilities of LISA to detect and reconstruct the spectral shape and amplitude of a stochastic gravitational wave background (SGWB). We first provide the LISA power-law sensitivity curve and binned power-law sensitivity curves, based on the latest updates on the LISA design. These curves are useful to make a qualitative assessment of the detection and reconstruction prospects of a SGWB. For a quantitative reconstruction of a SGWB with arbitrary power spectrum shape, we propose a novel data analysis technique: by means of an automatized adaptive procedure, we conveniently split the LISA sensitivity band into frequency bins, and fit the data inside each bin with a power law signal plus a model of the instrumental noise. We apply the procedure to SGWB signals with a variety of representative frequency profiles, and prove that LISA can reconstruct their spectral shape. Our procedure, implemented in the code SGWBinner, is suitable for homogeneous and isotropic SGWBs detectable at LISA, and it is also expected to work for other gravitational wave observatories.

This Astro-2020 White Paper deals with what we might learn from future gravitational wave observations about the early universe phase transitions and their energy scales, primordial black holes, Hubble parameter, dark matter and dark energy, modified theories of gravity and extra dimensions.

We present a lattice formulation of an interaction ${\phi\over\Lambda} F\tilde F$ between an axion and some $U(1)$ gauge sector with the following properties: it reproduces the continuum theory up to $\mathcal{O}(dx_\mu^2)$ corrections, it preserves exact gauge invariance and shift symmetry on the lattice, and it is suitable for self-consistent expansion of the Universe. The lattice equations of motion can no longer be solved by explicit methods, but we propose an implicit method to overcome this difficulty, which preserves the relevant system constraints down to arbitrary (tunable) precision. As a first application we study, in a comoving grid in $(3+1)$ dimensions, the last efolds of axion-inflation with quadratic potential and the preheating stage following afterwards. We fully account for the inhomogeneity and non-linearity of the system, including the gauge field contribution to the expansion rate of the Universe and its backreaction into the axion dynamics. We characterize in detail, as a function of the coupling, the energy transfer from the axion to the gauge field. Two coupling regimes are identified, sub- and super-critical, depending on whether the final energy fraction stored in the gauge field is below or above $\sim 50\%$ of the total energy. The Universe is very efficiently reheated for super-critical couplings, rapidly entering in a radiation dominated stage. Our results on preheating confirm previously published results.

We forecast the prospective of detection for a stochastic gravitational wave background sourced by cosmological first-order phase transitions. We focus on first-order phase transitions with negligible plasma effects, and consider the experimental infrastructures built by the end of the LISA mission. We make manifest the synergy among LISA, pulsar time array experiments, and ground-based interferometers. For phase transitions above the TeV scale or below the electroweak scale, LISA can detect the corresponding gravitational wave signal together with Einstein Telescope, SKA or even aLIGO-aVIRGO-KAGRA. For phase transitions at the electroweak scale, instead, LISA can be the only experiment observing the gravitational wave signal. In case of detection, by using a parameter reconstruction method that we anticipate in this work, we show that LISA on its own has the potential to determine when the phase transition occurs and, consequently, the energy scale above which the standard model of particle physics needs to be modified. The result may likely guide the collider community in the post-LHC era.

We study the post-inflationary dynamics of the Standard Model (SM) Higgs field in the presence of a non-minimal coupling $\xi|\Phi|^2R$ to gravity, both with and without the electroweak gauge fields coupled to the Higgs. We assume a minimal scenario in which inflation and reheating are caused by chaotic inflation with a quadratic potential, and no additional new physics is relevant below the Planck scale. By using classical real-time lattice simulations with a renormalisation group improved effective Higgs potential and by demanding the stability of the Higgs vacuum after inflation, we obtain upper bounds for $\xi$, taking into account the experimental uncertainty of the top-Yukawa coupling. We compare the bounds in the absence and presence of the electroweak gauge bosons, and conclude that the addition of gauge interactions has a rather minimal impact. In the unstable cases, we parametrize the time when such instability develops. For a top-quark mass $m_t \approx173.3 {\rm GeV}$, the Higgs vacuum instability is triggered for $\xi \gtrsim 4 -5$, although a slightly lower mass of $m_t \approx 172.1 {\rm GeV}$ pushes up this limit to $\xi \gtrsim 11 - 12$. This, together with the estimation $\xi \gtrsim 0.06$ for stability during inflation, provides tight constraints to the Higgs-curvature coupling within the SM.

Gravitational waves (GWs) have a great potential to probe cosmology. We review early universe sources that can lead to cosmological backgrounds of GWs. We begin by presenting definitions of GWs in flat space-time and in a cosmological setting, and discussing the reasons why GW backgrounds from the early universe are of a stochastic nature. We recap current observational constraints on stochastic backgrounds, and discuss some of the characteristics of present and future GW detectors including advanced LIGO, advanced Virgo, the Einstein Telescope, KAGRA, LISA. We then review in detail early universe GW generation mechanisms proposed in the literature, as well as the properties of the GW backgrounds they give rise to. We classify the backgrounds in five categories: GWs from quantum vacuum fluctuations during standard slow-roll inflation, GWs from processes that operate within extensions of the standard inflationary paradigm, GWs from post-inflationary preheating and related non-perturbative phenomena, GWs from first order phase transitions (related or not to the electroweak symmetry breaking), and GWs from topological defects, in particular from cosmic strings. The phenomenology of early universe processes that can generate a stochastic background of GWs is extremely rich, and some backgrounds are within the reach of near-future GW detectors. A future detection of any of these backgrounds will provide crucial information on the underlying high energy theory describing the early universe, probing energy scales well beyond the reach of particle accelerators.

Parametric resonance is among the most efficient phenomena generating gravitational waves (GWs) in the early Universe. The dynamics of parametric resonance, and hence of the GWs, depend exclusively on the resonance parameter $q$. The latter is determined by the properties of each scenario: the initial amplitude and potential curvature of the oscillating field, and its coupling to other species. Previous works have only studied the GW production for fixed value(s) of $q$. We present an analytical derivation of the GW amplitude dependence on $q$, valid for any scenario, which we confront against numerical results. By running lattice simulations in an expanding grid, we study for a wide range of $q$ values, the production of GWs in post-inflationary preheating scenarios driven by parametric resonance. We present simple fits for the final amplitude and position of the local maxima in the GW spectrum. Our parametrization allows to predict the location and amplitude of the GW background today, for an arbitrary $q$. The GW signal can be rather large, as $h^2\Omega_{\rm GW}(f_p) \lesssim 10^{-11}$, but it is always peaked at high frequencies $f_p \gtrsim 10^{7}$ Hz. We also discuss the case of spectator-field scenarios, where the oscillatory field can be e.g.~a curvaton, or the Standard Model Higgs.

We discuss the non-conservation of fermion number (or chirality breaking, depending on the fermionic charge assignment) in Abelian gauge theories at finite temperature. We study different mechanisms of fermionic charge disappearance in the high temperature plasma, with the use of both analytical estimates and real-time classical numerical simulations. We investigate the random walk of the Chern-Simons number $N_{\rm CS} \propto \int d^4x F_{\mu\nu}{\tilde F}^{\mu\nu}$, and show that it has a diffusive behaviour in the presence of an external magnetic field $B$. This indicates that the mechanism for fermionic number non-conservation for $B \neq 0$, is due to fluctuations of the gauge fields, similarly as in the case of non-Abelian gauge theories. We determine numerically the rate of chirality non-conservation associated with this diffusion, finding it larger by a factor $\sim 60$ compared to previous theoretical estimates. We also perform numerical simulations for the system which contains a chemical potential $\mu$ representing a fermionic charge density, again both with and without an external magnetic field. When $B=0$, we observe clearly the expected instability of the system for $\mu \neq 0$, as long as the chemical potential exceeds a critical value $\mu > \mu_c(L)$, which depends on the size $L$ of the system. When $B \neq 0$, the fluctuations of bosonic fields lead to the transfer of chemical potential into Chern-Simons number for arbitrary $\mu$.

If the Standard Model (SM) Higgs is weakly coupled to the inflationary sector, the Higgs is expected to be universally in the form of a condensate towards the end of inflation. The Higgs decays rapidly after inflation - via non-perturbative effects - into an out-of-equilibrium distribution of SM species, which thermalize soon afterwards. If the post-inflationary equation of state of the universe is stiff, $w \simeq +1$, the SM species eventually dominate the total energy budget. This provides a natural origin for the relativistic thermal plasma of SM species, required for the onset of the `hot Big Bang' era. The viability of this scenario requires the inflationary Hubble scale $H_*$ to be lower than the instability scale for Higgs vacuum decay, the Higgs not to generate too large curvature perturbations at cosmological scales, and the SM dominance to occur before Big Bang Nucleosynthesis. We show that successful reheating into the SM can only be obtained in the presence of a non-minimal coupling to gravity $\xi \gtrsim 1$, with a reheating temperature of $T_{\rm RH} \gtrsim \mathcal{O}(10^{10})\xi^{3/2}(H_*/10^{14}{\rm GeV})^2~{\rm GeV}$.

Real time evolution of classical gauge fields is relevant for a number of applications in particle physics and cosmology, ranging from the early Universe to dynamics of quark-gluon plasma. We present a lattice formulation of the interaction between a $shift$-symmetric field and some $U(1)$ gauge sector, $a(x)\tilde{F}_{\mu\nu}F^{\mu\nu}$, reproducing the continuum limit to order $\mathcal{O}(dx_\mu^2)$ and obeying the following properties: (i) the system is gauge invariant and (ii) shift symmetry is exact on the lattice. For this end we construct a definition of the {\it topological number density} $Q = \tilde{F}_{\mu\nu}F^{\mu\nu}$ that admits a lattice total derivative representation $Q = \Delta_\mu^+ K^\mu$, reproducing to order $\mathcal{O}(dx_\mu^2)$ the continuum expression $Q = \partial_\mu K^\mu \propto \vec E \cdot \vec B$. If we consider a homogeneous field $a(x) = a(t)$, the system can be mapped into an Abelian gauge theory with Hamiltonian containing a Chern-Simons term for the gauge fields. This allow us to study in an accompanying paper the real time dynamics of fermion number non-conservation (or chirality breaking) in Abelian gauge theories at finite temperature. When $a(x) = a(\vec x,t)$ is inhomogeneous, the set of lattice equations of motion do not admit however a simple explicit local solution (while preserving an $\mathcal{O}(dx_\mu^2)$ accuracy). We discuss an iterative scheme allowing to overcome this difficulty.

Particle production via parametric resonance in the early Universe, is a nonperturbative, non-linear and out-of-equilibrium phenomenon. Although it is a well studied topic, whenever a new scenario exhibits parametric resonance, a full re-analysis is normally required. To avoid this tedious task, many works present often only a simplified linear treatment of the problem. In order to surpass this circumstance in the future, we provide a fitting analysis of parametric resonance through all its relevant stages: initial linear growth, non-linear evolution, and relaxation towards equilibrium. Using lattice simulations in an expanding grid in 3 + 1 dimensions, we parametrize the dynamics outcome scanning over the relevant ingredients: role of the oscillatory field, particle coupling strength, initial conditions, and background expansion rate. We emphasize the inaccuracy of the linear calculation of the decay time of the oscillatory field, and propose a more appropriate definition of this scale based on the subsequent non-linear dynamics. We provide simple fits to the relevant time scales and particle energy fractions at each stage. Our fits can be applied to post-inflationary preheating scenarios, where the oscillatory field is the inflaton, or to spectator-field scenarios, where the oscillatory field can be e.g. a curvaton, or the Standard Model Higgs.

We investigate the potential for the LISA space-based interferometer to detect the stochastic gravitational wave background produced from different mechanisms during inflation. Focusing on well-motivated scenarios, we study the resulting contributions from particle production during inflation, inflationary spectator fields with varying speed of sound, effective field theories of inflation with specific patterns of symmetry breaking and models leading to the formation of primordial black holes. The projected sensitivities of LISA are used in a model-independent way for various detector designs and configurations. We demonstrate that LISA is able to probe these well-motivated inflationary scenarios beyond the irreducible vacuum tensor modes expected from any inflationary background.

During or towards the end of inflation, the Standard Model (SM) Higgs forms a condensate with a large amplitude. Following inflation, the condensate oscillates, decaying non-perturbatively into the rest of the SM species. The resulting out-of-equilibrium dynamics converts a fraction of the energy available into gravitational waves (GW). We study this process using classical lattice simulations in an expanding box, following the energetically dominant electroweak gauge bosons $W^\pm$ and $Z$. We characterize the GW spectrum as a function of the running couplings, Higgs initial amplitude, and post-inflationary expansion rate. As long as the SM is decoupled from the inflationary sector, the generation of this background is universally expected, independently of the nature of inflation. Our study demonstrates the efficiency of GW emission by gauge fields undergoing parametric resonance. The initial energy of the Higgs condensate represents however, only a tiny fraction of the inflationary energy. Consequently, the resulting background is very suppressed, with an amplitude $h^2 \Omega_{\rm GW}^{(o)} \lesssim 10^{-29}$ today. The amplitude can be boosted to $h^2 \Omega_{\rm GW}^{(o)} \lesssim 10^{-16}$, if following inflation the universe undergoes a kination-domination stage; however the background is shifted in this case to high frequencies $f_p \lesssim 10^{11} {\rm Hz}$. In all cases the signal is out of the range of current or planned GW detectors. This background will therefore remain, most likely, as a curiosity of the SM.

We study the nonperturbative dynamics of the Standard Model (SM) after inflation, in the regime where the SM is decoupled from (or weakly coupled to) the inflationary sector. We use classical lattice simulations in an expanding box in (3+1) dimensions, modeling the SM gauge interactions with both global and Abelian-Higgs analogue scenarios. We consider different post-inflationary expansion rates. During inflation, the Higgs forms a condensate, which starts oscillating soon after inflation ends. Via nonperturbative effects, the oscillations lead to a fast decay of the Higgs into the SM species, transferring most of the energy into $Z$ and $W^{\pm}$ bosons. All species are initially excited far away from equilibrium, but their interactions lead them into a stationary stage, with exact equipartition among the different energy components. From there on the system eventually reaches equilibrium. We have characterized in detail, in the different expansion histories considered, the evolution of the Higgs and of its dominant decay products, until equipartition is established. We provide a useful mapping between simulations with different parameters, from where we derive a master formula for the Higgs decay time, as a function of the coupling constants, Higgs initial amplitude and postinflationary expansion rate.

The apparent discrepancy between the value of the tensor-to-scalar ratio reported by the BICEP2 collaboration, $r = 0.20^{+0.07}_{-0.05}$ at 68% CL, and the Planck upper limit, $r < 0.11$ at 95% CL, has attracted a great deal of attention. In this short note, we show that this discrepancy is mainly due to an `apples to oranges' comparison. The result reported by BICEP2 was measured at a pivot scale $k_* = 0.05$ Mpc$^{-1}$, assuming $n_t = 0$, whereas the Planck limit was provided at $k_* = 0.002$ Mpc$^{-1}$, assuming the slow-roll consistency relation $n_t = -r/8$. One should obviously compare the BICEP2 and Planck results under the same circumstances. By imposing $n_t = 0$, the Planck constraint at $k_* = 0.05$ Mpc$^{-1}$ becomes $r < 0.135$ at $95%$ CL, which can be compared directly with the BICEP2 result. Once a plausible dust contribution to the BICEP2 signal is taken into account (DDM2 model), $r$ is reduced to $r = 0.16^{+0.06}_{-0.05}$ and the discrepancy becomes of order $1.3\sigma$ only.

We show that self-ordering scalar fields (SOSF), i.e. non-topological cosmic defects arising after a global phase transition, cannot explain the B-mode signal recently announced by BICEP2. We compute the full $C_\ell^{B}$ angular power spectrum of B-modes due to the vector and tensor perturbations of SOSF, modeled in the large-N limit of a spontaneous broken global O(N) symmetry. We conclude that the low-$\ell$ multipoles detected by BICEP2 cannot be due mainly to SOSF, since they have the wrong spectrum at low multipoles. As a byproduct we derive the first cosmological constraints on this model, showing that the BICEP2 B-mode polarization data admits at most a 2-3% contribution from SOSF in the temperature anisotropies, similar to (but somewhat tighter than) the recently studied case of cosmic strings.

We determine the full C_l spectra and correlation functions of the temperature and polarization anisotropies in the CMB, generated by a source modeled by the large N limit of spontaneously broken global O(N)-theories. We point out a problem in the standard approach of treating the radiation-matter transition by interpolating the eigenvectors of the unequal-time correlators of the source energy-momentum tensor. This affects the CMB predictions from all type of cosmic defects. We propose a method to overcome this difficulty, and find that in the large-N global model that we study, differences in the final CMB power spectra amplitudes reach up to 25%, when compared to implementations of the eigenvector interpolation technique. We discuss as well how to optimally search for the contribution in the CMB from active sources such as cosmic defects, in experiments like Planck, COrE and PRISM.

The existence of the Standard Model (SM) Higgs implies that a gravitational wave (GW) background is generated by the decay products of the Higgs, soon after the end of inflation. Theoretically, all Yukawa and $SU(2)_L$ gauge couplings of the SM are imprinted as features in the GW spectrum. In practice, the signal from the most strongly coupled species dominates, rendering inaccesible the information on the other species. If detected, this background could be used for measuring properties of high-energy particle physics, including beyond the SM. To achieve this goal, new high frequency GW detection technology is required, beyond that of currently planned detectors.

Out-of-equilibrium fermions can be created in the early Universe by non-perturbative parametric effects, both at preheating or during the thermal era. An anisotropic stress is developed in the fermion distribution, acting as a source of a stochastic background of gravitational waves (GW). We derive a general formalism to calculate the spectrum of GW produced by an ensemble of fermions, which we apply to a variety of scenarios after inflation. We discuss in detail the regularization of the source, i.e. of the unequal-time-correlator of the fermions' transverse-traceless anisotropic stress. We discuss how the GW spectrum builds up in time and present a simple parametrization of its final amplitude and peak frequency. We find that fermions may generate a GW background with a significant amplitude at very high frequencies, similarly to the case of preheating with scalar fields. A detection of this GW background would shed light about the physics of the very early Universe, but new technology at high frequencies is required, beyond the range accessible to currently planned detectors.

When a light scalar field is present during inflation, its value will vary on superhorizon scales, modulating the preheating process at the end of inflation. Consequently, the amplitude of the gravitational wave (GW) background produced during preheating will also be modulated. The observed energy density of this background will therefore appear anisotropic at different angles in the sky. We provide a master formula for the angular power spectrum C_l of the anisotropies in the GW background from preheating, valid for any scenario where the anisotropies are due to the superhorizon modulation of a light degree of freedom. Using lattice field theory simulations of massless preheating with g^2/\lambda = 2, we find a flat angular spectrum l(l+1)C_l \approx 3x10^{-4}, which represents a strong anisotropy of order 1% variations on large angular scales. For our choice of couplings, long wavelengths are amplified most strongly during parametric resonance, which is crucial for the development of the anisotropies. If future direct detection GW observatories are capable of detecting backgrounds of cosmological origin, they should be able to detect this effect. This could eventually become a powerful tool to discriminate among inflationary and preheating scenarios.

We investigate the anisotropies in the gravitational wave (GW) background produced at preheating after inflation. Using lattice field theory simulations of a massless preheating model, we show that the GW amplitude depends sensitively on the value of the decay product field chi coupled to the inflaton phi, with the only requisite that chi is light during inflation. We find a strong anisotropy in the amplitude of the GW background on large angular scales, the details of which strongly depend on the reheating dynamics. We expect similar conclusions for a wide class of inflationary models with light scalar fields. If future direct detection GW experiments are capable of detecting the GW produced by preheating, they should also be able to detect this effect. This could eventually provide a powerful way to distinguish between different inflationary and preheating scenarios.

We demonstrate that any scaling source in the radiation era produces a background of gravitational waves with an exact scale-invariant power spectrum. Cosmic defects, created after a phase transition in the early Universe, are such a scaling source. We emphasise that the result is independent of the topology of the cosmic defects, the order of phase transition, and the nature of the symmetry broken, global or gauged. As an example, using large-scale numerical simulations, we calculate the scale invariant gravitational wave power spectrum generated by the dynamics of a global O(N) scalar theory. The result approaches the large N theoretical prediction as N^(-2), albeit with a large coefficient. The signal from global cosmic strings is O(100) times larger than the large N prediction.

We investigate in detail a model where the curvaton is coupled to the Standard Model higgs. Parametric resonance might be expected to cause a fast decay of the curvaton, so that it would not have time to build up the curvature perturbation. However, we show that this is not the case, and that the resonant decay of the curvaton may be delayed even down to electroweak symmetry breaking. This delay is due to the coupling of the higgs to the thermal background, which is formed by the Standard Model degrees of freedom created from the inflaton decay. We establish the occurrence of the delay by considering the curvaton evolution and the structure of the higgs resonances. We then provide analytical expressions for the delay time, and for the subsequent resonant production of the higgs, which ultimately leads to the curvaton effective decay width. Contrary to expectations, it is possible to obtain the observed curvature perturbation for values of the curvaton-higgs coupling as large as 0.1. Our calculations also apply in the general case of curvaton decay into any non Standard Model species coupled to the thermal background.

We compute the matter bispectrum in the presence of primordial local non-Gaussianity over a wide range of scales, including the very small nonlinear ones. We use the Halo Model approach, considering non-Gaussian corrections to the halo profiles, the halo mass function and the bias functions. We compare our results in the linear and mildly nonlinear scales to a large ensemble of Gaussian and non-Gaussian numerical simulations. We consider both squeezed and equilateral configurations, at redshift z = 0 and z = 1. For z = 0, the deviations between the Halo Model and the simulations are smaller than 10% in the squeezed limit, both in the Gaussian and non-Gaussian cases. The Halo Model allows to make predictions on scales much smaller than those reached by numerical simulations. For local non-Gaussian initial conditions with a parameter fNL = 100, we find an enhancement of the bispectrum in the squeezed configuration k = k3 = k2 >> k1 \sim 0.01 h^{-1} Mpc, of \sim 15% and \sim 25% on scales k \sim 1 h^{-1} Mpc, at z = 0 and z = 1 respectively. This is mainly due to the non-Gaussian corrections in the linear bias. Finally we provide a very simple expression valid for any scenario, i.e. for any choice of the halo profile, mass and bias functions, which allow for a fast evaluation of the bispectrum on squeezed configurations.

Preheating and other particle production phenomena in the early Universe can give rise to high- energy out-of-equilibrium fermions with an anisotropic stress. We develop a formalism to calculate the spectrum of gravitational waves due to fermions, and apply it to a variety of scenarios after inflation. We pay particular attention to regularization issues. We show that fermion production sources a stochastic background of gravitational waves with a significant amplitude, but we find that typical frequencies of this new background are not within the presently accessible direct detec- tion range. However, small-coupling scenarios might still produce a signal observable by planned detectors, and thus open a new window into the physics of the very early Universe.

It has recently been pointed out that the usual procedure employed in order to obtain the transverse-traceless (TT) part of metric perturbations in lattice simulations was inconsistent with the fact that those fields live in the lattice and not in the continuum. It was claimed that this could lead to a larger amplitude and a wrong shape for the gravitational wave (GW) spectra obtained in numerical simulations of (p)reheating. In order to address this issue, we have defined a consistent prescription in the lattice for extracting the TT part of the metric perturbations. We demonstrate explicitly that the GW spectra obtained with the old continuum-based TT projection only differ marginally in amplitude and shape with respect to the new lattice-based ones. We conclude that one can therefore trust the predictions appearing in the literature on the spectra of GW produced during (p)reheating and similar scenarios simulated on a lattice.

We consider a set of scalar fields, consisting of a single flat direction and one or several non-flat directions. We take our cue from the MSSM, considering separately D-flat and F-flat directions, but our results apply to any supersymmetric scenario containing flat directions. We study the field fluctuations during pure de Sitter Inflation, following the evolution of the infrared modes by numerically solving the appropriate Langevin equations. We demonstrate that for the Standard Model U(1), SU(2) or SU(3) gauge couplings, as well as for large enough Yukawa couplings, the fluctuations along the non-flat directions effectively block the fluctuations along the flat directions. The usual expected behaviour <{\phi}^2> \propto N, with N the number of efolds, may be strongly violated, depending on the coupling strengths. As a consequence, those cosmological considerations, which are derived assuming that during inflation flat directions fluctuate freely, should be revised.

We present a new and especially powerful signature of cosmic strings and other topological or non-topological defects in the polarization of the cosmic microwave background (CMB). We show that even if defects contribute 1% or less in the CMB temperature anisotropy spectrum, their signature in the local $\tilde{B}$-polarization correlation function at angular scales of tens of arc minutes is much larger than that due to gravitational waves from inflation, even if the latter contribute with a ratio as big as $r\simeq 0.1$ to the temperature anisotropies. We show that when going from non-local to local $\tilde{B}$-polarization, the ratio of the defect signal-to-noise with respect to the inflationary value increases by about an order of magnitude. Proposed B-polarization experiments, with a good sensitivity on arcminute scales, may either detect a contribution from topological defects produced after inflation or place stringent limits on them. Even Planck should be able to improve present constraints on defect models by at least an order of magnitude, to the level of $\ep <10^{-7}$. A future full-sky experiment like CMBpol, with polarization sensitivities of the order of $1\mu$K-arcmin, will be able to constrain the defect parameter $\ep=Gv^2$ to a few $\times10^{-9}$, depending on the defect model.

Primordial gravitational waves provide a very important stochastic background that could be detected soon with interferometric gravitational wave antennas or indirectly via the induced patterns in the polarization anisotropies of the cosmic microwave background. The detection of these waves will open a new window into the early Universe, and therefore it is important to characterize in detail all possible sources of primordial gravitational waves. In this paper we develop theoretical and numerical methods to study the production of gravitational waves from out-of-equilibrium gauge fields at preheating. We then consider models of preheating after hybrid inflation, where the symmetry breaking field is charged under a local U(1) symmetry. We analyze in detail the dynamics of the system in both momentum and configuration space, and show that gauge fields leave specific imprints in the resulting gravitational wave spectra, mainly through the appearence of new peaks at characteristic frequencies that are related to the mass scales in the problem. We also show how these new features in the spectra correlate with string-like spatial configurations in both the Higgs and gauge fields that arise due to the appearance of topological winding numbers of the Higgs around Nielsen-Olesen strings. We study in detail the time evolution of the spectrum of gauge fields and gravitational waves as these strings evolve and decay before entering a turbulent regime where the gravitational wave energy density saturates.

The Universe may harbor relics of the post-inflationary epoch in the form of a network of self-ordered scalar fields. Such fossils, while consistent with current cosmological data at trace levels, may leave too weak an imprint on the cosmic microwave background and the large-scale distribution of matter to allow for direct detection. The non-Gaussian statistics of the density perturbations induced by these fields, however, permit a direct means to probe for these relics. Here we calculate the bispectrum that arises in models of self-ordered scalar fields. We find a compact analytic expression for the bispectrum, evaluate it numerically, and provide a simple approximation that may be useful for data analysis. The bispectrum is largest for triangles that are aligned (have edges $k_1\simeq 2 k_2 \simeq 2 k_3$) as opposed to the local-model bispectrum, which peaks for squeezed triangles ($k_1\simeq k_2 \gg k_3$), and the equilateral bispectrum, which peaks at $k_1\simeq k_2 \simeq k_3$. We estimate that this non-Gaussianity should be detectable by the Planck satellite if the contribution from self-ordering scalar fields to primordial perturbations is near the current upper limit.

We discuss Preheating after an inflationary stage driven by the Standard Model (SM) Higgs field non-minimally coupled to gravity. We find that Preheating is driven by a complex process in which perturbative and non-perturbative effects occur simultaneously. The Higgs field, initially an oscillating coherent condensate, produces non-perturbatively W and Z gauge fields. These decay very rapidly into fermions, thus preventing gauge bosons to accumulate and, consequently, blocking the usual parametric resonance. The energy transferred into the fermionic species is, nevertheless, not enough to reheat the Universe, and resonant effects are eventually developed. Soon after resonance becomes effective, also backreaction from the gauge bosons into the Higgs condensate becomes relevant. We have determined the time evolution of the energy distribution among the remnant Higgs condensate and the non-thermal distribution of the SM fermions and gauge fields, until the moment in which backreaction becomes important. Beyond backreaction our approximations break down and numerical simulations and theoretical considerations beyond this work are required, in order to study the evolution of the system until thermalization.

Gravitational waves were copiously produced in the early Universe whenever the processes taking place were sufficiently violent. The spectra of several of these gravitational wave backgrounds on subhorizon scales have been extensively studied in the literature. In this paper we analyze the shape and amplitude of the gravitational wave spectrum on scales which are superhorizon at the time of production. Such gravitational waves are expected from the self ordering of randomly oriented scalar fields which can be present during a thermal phase transition or during preheating after hybrid inflation. We find that, if the gravitational wave source acts only during a small fraction of the Hubble time, the gravitational wave spectrum at frequencies lower than the expansion rate at the time of production behaves as $\Omega_{\rm GW}(f) \propto f^3$ with an amplitude much too small to be observable by gravitational wave observatories like LIGO, LISA or BBO. On the other hand, if the source is active for a much longer time, until a given mode which is initially superhorizon ($k\eta_* \ll 1$), enters the horizon, for $k\eta \gtrsim 1$, we find that the gravitational wave energy density is frequency independent, i.e. scale invariant. Moreover, its amplitude for a GUT scale scenario turns out to be within the range and sensitivity of BBO and marginally detectable by LIGO and LISA. This new gravitational wave background can compete with the one generated during inflation, and distinguishing both may require extra information.

We study the details of preheating in an inflationary scenario in which the Standard Model Higgs, strongly non-minimally coupled to gravity, plays the role of the inflaton. We find that the Universe does not reheat immediately through perturbative decays, but rather initiate a complex process in which perturbative and non-perturbative effects are mixed. The Higgs condesate starts oscillating around the minimum of its potential, producing W and Z gauge bosons non-perturbatively, due to violation of the so-called adiabaticity condition. However, during each semi-oscillation, the created gauge bosons completely decay (perturbatively) into fermions. This way, the decay of the gauge bosons prevents the development of parametric resonance, since bosons cannot accummulate significantly at the beginning. However, the energy transferred to the decay products of the bosons is not enough to reheat the universe, so after about a hundred oscillations, the resonance effect will finally dominate over the perturbative decays. Around the same time (or slightly earlier), backreaction from the gauge bosons onto the Higgs condensate will also start to be significant. Soon afterwards, the Universe is filled with the remnant condensate of the Higgs and a non-thermal distribution of Standard Model particles which redshift as radiation, while the Higgs continues to oscillate as a pressureless fluid. We compute the distribution of energy among all the species present at backreaction. From there on until thermalization, the evolution of the system is highly non-linear and non-perturbative, and will require a careful study via numerical simulations.

We discuss the cosmological degeneracy between the Hubble parameter H(z), the age of the universe and cosmological parameters describing simple variations from the minimal LCDM model. We show that independent determinations of the Hubble parameter H(z) such as those recently provided by Simon et al., Phys. Rev. D 71 (2005) 123001, combined with Cosmic Microwave Background data, can provide stringent constraints on possible deviations from the LCDM model. In particular we find that this data combination constrains at the 68% (95%) confidence level the following parameters: sum of the neutrino masses < 0.5 (1.0) eV, number of relativistic neutrino species N_rel = 4.1^{+0.4}_{-0.9} (^{+1.1}_{-1.5}), dark energy equation of state parameter w = -0.95 +- 0.17 (+- 0.32), and curvature \Omega_k = 0.002 +- 0.006 (+- 0.014), in excellent agreement with dataset combinations involving Cosmic Microwave Background, Supernovae and Baryon Acoustic Oscillations. This offers a valuable consistency check for systematic errors.

The reheating of the universe after hybrid inflation proceeds through the nucleation and subsequent collision of large concentrations of energy density in the form of bubble-like structures moving at relativistic speeds. This generates a significant fraction of energy in the form of a stochastic background of gravitational waves, whose time evolution is determined by the successive stages of reheating: First, tachyonic preheating makes the amplitude of gravity waves grow exponentially fast. Second, bubble collisions add a new burst of gravitational radiation. Third, turbulent motions finally sets the end of gravitational waves production. From then on, these waves propagate unimpeded to us. We find that the fraction of energy density today in these primordial gravitational waves could be significant for GUT scale models of inflation, although well beyond the frequency range sensitivity of gravitational wave observatories like LIGO, LISA or BBO. However, low-scale models could still produce a detectable signal at frequencies accessible to BBO or DECIGO. For comparison, we have also computed the analogous background from some chaotic inflation models and obtained similar results to those of other groups. The discovery of such a background would open a new observational window into the very early universe, where the details of the process of reheating could be explored. Thus, it could also serve as a new experimental tool for testing the Inflationary Paradigm.

The process of reheating the universe after hybrid inflation is extremely violent. It proceeds through the nucleation and subsequent collision of large concentrations of energy density in bubble-like structures, which generate a significant fraction of energy in the form of gravitational waves. We study the power spectrum of the stochastic background of gravitational waves produced at reheating after hybrid inflation. We find that the amplitude could be significant for high-scale models, although the typical frequencies are well beyond what could be reached by planned gravitational wave observatories like LIGO, LISA or BBO. On the other hand, low-scale models could still produce a detectable stochastic background at frequencies accesible to those detectors. The discovery of such a background would open a new window into the very early universe.