We demonstrate that light spectator fields can source sizeable CMB anisotropies through modulated reheating even in the absence of direct couplings to the inflaton. The effect arises when the phase space of the inflaton decay is modulated by the spectator which generates masses for the decay products. We call the mechanism \textit{indirect modulation} and show that it can source perturbations even four orders of magnitude larger than the observed. Importantly, the indirect mechanism is present in the Standard Model extended with right-handed neutrinos. For a minimally coupled Higgs boson this leads to a novel lower bound on the quartic coupling and constrains the neutrino Yukawas below unity.

We use the spectral representation of the stochastic Starobinsky-Yokoyama approach to compute correlation functions in de Sitter space for a scalar field with a symmetric or asymmetric double-well potential. The terms in the spectral expansion are determined by the eigenvalues and eigenfunctions of the time-independent Fokker-Planck differential operator, and we solve them numerically. The long-distance asymptotic behaviour is given by the lowest state in the spectrum, but we demonstrate that the magnitude of the coeffients of different terms can be very different, and the correlator can be dominated by different terms at different distances. This can give rise to potentially observable cosmological signatures. In many cases the dominant states in the expansion do not correspond to small fluctuations around a minimum of the potential and are therefore not visible in perturbation theory. We discuss the physical interpretation these states, which can be present even when the potential has only one minimum.

We consider light scalar fields during inflation and show how the stochastic spectral expansion method can be used to calculate two-point correlation functions of an arbitrary local function of the field in de Sitter space. In particular, we use this approach for a massive scalar field with quartic self-interactions to calculate the fluctuation spectrum of the density contrast and compare it to other approximations. We find that neither Gaussian nor linear approximations accurately reproduce the power spectrum, and in fact always overestimate it. For example, for a scalar field with only a quartic term in the potential, $V=\lambda\phi^4/4$, we find a blue spectrum with spectral index $n-1=0.579\sqrt{\lambda}$.

We study the contribution to the primordial curvature perturbation on observational scales generated by the reheating field in massless preheating. To do so we use lattice simulations and a recent extension to the $\delta N$ formalism. The work demonstrates the functionality of these techniques for calculating the observational signatures of models in which non-perturbative reheating involves a light scalar field.

The current central experimental values of the parameters of the Standard Model give rise to a striking conclusion: metastability of the electroweak vacuum is favoured over absolute stability. A metastable vacuum for the Higgs boson implies that it is possible, and in fact inevitable, that a vacuum decay takes place with catastrophic consequences for the Universe. The metastability of the Higgs vacuum is especially significant for cosmology, because there are many mechanisms that could have triggered the decay of the electroweak vacuum in the early Universe. We present a comprehensive review of the implications from Higgs vacuum metastability for cosmology along with a pedagogical discussion of the related theoretical topics, including renormalization group improvement, quantum field theory in curved spacetime and vacuum decay in field theory.

The observed dark matter abundance in the Universe can be fully accounted for by a minimally coupled spectator scalar field that was light during inflation and has sufficiently strong self-coupling. In this scenario, dark matter was produced during inflation by amplification of quantum fluctuations of the spectator field. The self-interaction of the field suppresses its fluctuations on large scales, and therefore avoids isocurvature constraints. The scenario does not require any fine-tuning of parameters. In the simplest case of a single real scalar field, the mass of the dark matter particle would be in the range $1~{\rm GeV}\lesssim m\lesssim 10^8~{\rm GeV}$, depending on the scale of inflation, and the lower bound for the quartic self-coupling is $\lambda\gtrsim 0.45$.

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.

We revisit the question of how to calculate correlations of the curvature perturbation, $\zeta$, using the $\delta N$ formalism when one cannot employ a truncated Taylor expansion of $N$. This problem arises when one uses lattice simulations to probe the effects of isocurvature modes on models of reheating. Working in real space, we use an expansion in the cross-correlation between fields at different positions, and present simple expressions for observables such as the power spectrum and the reduced bispectrum, $f_{\rm NL}$. These take the same form as those of the usual $\delta N$ expressions, but with the derivatives of $N$ replaced by non-perturbative $\delta N$ coefficients. We test the validity of this expansion and argue that our expressions are particularly well suited for use with simulations.

We show that the dynamics of the Higgs field during inflation is not affected by metric fluctuations if the Higgs is an energetically subdominant light spectator. For Standard Model parameters we find that couplings between Higgs and metric fluctuations are suppressed by $\mathcal{O}(10^{-7})$. They are negligible compared to both pure Higgs terms in the effective potential and the unavoidable non-minimal Higgs coupling to background scalar curvature. The question of the electroweak vacuum instability during high energy scale inflation can therefore be studied consistently using the Jordan frame action in a Friedmann--Lema\^itre--Robertson--Walker metric, where the Higgs-curvature coupling enters as an effective mass contribution. Similar results apply for other light spectator scalar fields during inflation.

Magnetic monopoles, if they exist, would be produced amply in strong magnetic fields and high temperatures via the thermal Schwinger process. Such circumstances arise in heavy ion collisions and in neutron stars, both of which imply lower bounds on the mass of possible magnetic monopoles. In showing this, we construct the cross section for pair production of magnetic monopoles in heavy ion collisions, which indicates that they are particularly promising for experimental searches such as MoEDAL.

The behaviour of a massive, non-interacting and non-minimally coupled quantised scalar field in an expanding de Sitter background is investigated by solving the field evolution for an arbitrary initial state. In this approach there is no need to choose a vacuum in order to provide a definition for particle states. We conclude that the expanding de Sitter space is a stable equilibrium configuration under small perturbations of the initial conditions. Depending on the initial state, the energy density can approach its asymptotic value from above or below, the latter of which implies a violation of the weak energy condition. The backreaction of the quantum corrections can therefore lead to a phase of super-acceleration also in the non-interacting massive case.

We investigate the dynamics of the Higgs field at the end of inflation in the minimal scenario consisting of an inflaton field coupled to the Standard Model only through the non-minimal gravitational coupling $\xi$ of the Higgs field. Such a coupling is required by renormalisation of the Standard Model in curved space, and in the current scenario also by vacuum stability during high-scale inflation. We find that for $\xi\gtrsim 1$, rapidly changing spacetime curvature at the end of inflation leads to significant production of Higgs particles, potentially triggering a transition to a negative-energy Planck scale vacuum state and causing an immediate collapse of the Universe.

It has been claimed that the electroweak vacuum may be unstable during inflation due to large fluctuations of order $H$ in case of a high inflationary scale as suggested by BICEP2. We compute the Standard Model Higgs effective potential including UV-induced curvature corrections at one-loop level. We find that for a high inflationary scale a large curvature mass is generated due to RG running of non-minimal coupling $\xi$, which either stabilizes the potential against fluctuations for $\xi_{\rm EW} \gtrsim 6\cdot 10^{-2}$, or destabilizes it for $\xi_{\rm EW} \lesssim 2 \cdot 10^{-2}$ when the generated curvature mass is negative. Only in the narrow intermediate region the effect of the curvature mass may be significantly smaller.

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.

Light scalar fields present during inflation can lead to interesting observable signatures, especially in models with non-equilibrium reheating dynamics. We study a supersymmetric hybrid inflation model with a third scalar whose lightness is protected by symmetry, using analytical and numerical techniques, and demonstrate that the amplitude, the spectral index and the non-Gaussianity parameter fNL of the primordial curvature perturbation are within the Planck observational bounds for suitable parameter values.

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 consider the possibility that the primordial curvature perturbation was generated through the curvaton mechanism from a scalar field with an electric charge, or precisely the Standard Model U(1) weak hypercharge. This links the dynamics of the very early universe concretely to the Standard Model of particle physics, and because the coupling strength is known, it reduces the number of free parameters in the curvaton model. The gauge coupling also introduces several new physical effects. Charge fluctuations are generated during inflation, but they are screened by electron-positron pairs therefore do not violate observational constraints. After inflation, the curvaton interacts with thermal radiation which destroys the curvaton condensate and prevents the generation of curvature perturbations, unless the inflaton dynamics satisfy strong constraints. The curvaton also experiences a period of parametric resonance with the U(1) gauge field. Using the standard perturbative approach, we find that the model can generate the observed density perturbation for Hubble rate H_* > 10^8 GeV and curvaton mass m > 0.01 H_*, but with a level of non-Gaussianity (f_NL > 130) that violates observational constraints. However, previous studies have shown that the parametric resonance changes the predicted perturbations significantly, and therefore fully non-linear numerical field theory simulations are required.

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 study the hybrid inflationary potential in a regime where the defect field is light, and more than 60 e-folds of accelerated expansion occur after the symmetry breaking transition. Using analytic and numerical techniques, we then identify parameter values within this regime for which the statistics of the primordial curvature perturbation are significantly non-Gaussian. Focusing on this range of parameters, we provide a specific example which leads to an observationally consistent power spectrum, and a level of non-Gaussianity within current WMAP bounds and in reach of the Planck satellite. An interesting feature of this example is that the initial conditions at horizon crossing appear quite natural.

We calculate curvature perturbations in the scenario in which the curvaton field decays into another scalar field via parametric resonance. As a result of a nonlinear stage at the end of the resonance, standard perturbative calculation techniques fail in this case. Instead, we use lattice field theory simulations and the separate universe approximation to calculate the curvature perturbation as a nonlinear function of the curvaton field. For the parameters tested, the generated perturbations are highly non-Gaussian and not well approximated by the usual fNL parameterisation. Resonant decay plays an important role in the curvaton scenario and can have a substantial effect on the resulting perturbations.

When gauged cosmic strings form in a symmetry-breaking phase transition, the gauge field configuration at the time becomes imprinted in the spatial string distribution by the flux trapping mechanism. It is shown that quantum and thermal gauge field fluctuations give rise to long-range superhorizon correlations in the string network. Numerical simulations in the Abelian Higgs model confirm this finding. In contrast, the Kibble-Zurek mechanism, which most cosmic string studies are based on, only gives rise to short-distance subhorizon correlations. The potential observable effects of the correlations are discussed.

If light scalar fields are present at the end of inflation, their non-equilibrium dynamics such as parametric resonance or a phase transition can produce non-Gaussian density perturbations. We show how these perturbations can be calculated using non-linear lattice field theory simulations and the separate universe approximation. In the massless preheating model, we find that some parameter values are excluded while others lead to acceptable but observable levels of non-Gaussianity. This shows that preheating can be an important factor in assessing the viability of inflationary models.

Preheating can convert superhorizon fluctuations of light scalar fields present at the end of inflation into observable density perturbations. We show in detail how lattice field theory simulations and the separate universes approximation can be used to calculate these perturbations and make predictions for the nonlinearity parameter fNL . We also present a simple approximation scheme that can reproduce these results analytically. Applying these methods to the massless preheating model, we determine the parameter values that are ruled out by too high levels of non-Gaussianity.

We investigate the behaviour of (p,q) string networks, focusing on two aspects: (1) modelling more realistic (p,q) string networks than the Z_N networks used so far and (2) investigating the effect of long-range interactions on the evolution of the network. We model the network with no long-range interactions using two sets of fields, complex scalars coupled to gauge fields, with a potential chosen such that the two types of strings will form bound states. This way we can model junctions of 3 strings with different tension; in Z_N models used so far in simulations all the strings have identical tensions. In order to introduce long-range interactions we also study a network in which one of the scalars forms global strings. We observe that in the absence of long-range interactions the formation of bound states has a significant influence on the evolution of the network. When long-range interactions are turned on the bound states are short-lived and have a minimal effect on the network evolution.

We argue that cosmic strings with high winding numbers generally form in first-order gauge symmetry breaking phase transitions, and we demonstrate this using computer simulations. These strings are heavier than single-winding strings and therefore more easily observable. Their cosmological evolution may also be very different.

We investigate the end of the inflationary period in the recently proposed scenario of locked inflation, and consider various constraints arising from density perturbations, loop corrections, parametric resonance and defect formation. We show that in a scenario where there is one long period of locked inflation, it is not possible to satisfy all of these constraints without having a period of saddle inflation afterwards, which would wipe out all observable signatures. On the other hand, if one does not insist on satisfying the loop correction constraint, saddle inflation can be avoided, but then inflation must have ended through parametric resonance.

We carry out numerical simulations to investigate spontaneous vortex formation during a temperature quench of a superconductor film from the normal to the superconducting phase in the absence of an external magnetic field. Our results agree roughly with quantitative predictions of the flux trapping scenario: In fast quenches the agreement is almost perfect, but there appears to be some discrepancy in slower ones. In particular, our simulations demonstrate the crucial role the electromagnetic field plays in this phenomenon, making it very different from vortex formation in superfluids. Besides superconductor experiments, our findings also shed more light on the possible formation of cosmic defects in the early universe.

Thermal fluctuations of the gauge field lead to monopole formation at the grand unified phase transition in the early Universe, even if the transition is merely a smooth crossover. The dependence of the produced monopole density on various parameters is qualitatively different from theories with global symmetries, and the monopoles have a positive correlation at short distances. The number density of monopoles may be suppressed if the grand unified symmetry is only restored for a short time by, for instance, nonthermal symmetry restoration after preheating.

Topological defects are common in many everyday systems. In general, they appear if a symmetry is broken at a rapid phase transition. In this article, I explain why it is believed that they should have also been formed in the early universe and how that would have happened. If topological defects are found, this will provide a way to study observationally the first fractions of a second after the Big Bang, but their apparent absence can also tell us many things about the early universe.

We study the three-dimensional Georgi-Glashow model to demonstrate how magnetic monopoles can be studied fully non-perturbatively in lattice Monte Carlo simulations, without any assumptions about the smoothness of the field configurations. We examine the apparent contradiction between the conjectured analytic connection of the `broken' and `symmetric' phases, and the interpretation of the mass (i.e., the free energy) of the fully quantised 't Hooft-Polyakov monopole as an order parameter to distinguish the phases. We use Monte Carlo simulations to measure the monopole free energy and its first derivative with respect to the scalar mass. On small volumes we compare this to semi-classical predictions for the monopole. On large volumes we show that the free energy is screened to zero, signalling the formation of a confining monopole condensate. This screening does not allow the monopole mass to be interpreted as an order parameter, resolving the paradox.

We study the instability of a scalar field at the end of hybrid inflation, using both analytical techniques and numerical simulations. We improve previous studies by taking the inflaton field fully into account, and show that the range of unstable modes depends sensitively on the velocity of the inflaton field, and thereby on the Hubble rate, at the end of inflation. If topological defects are formed, their number density is determined by the shortest unstable wavelength. Finally, we show that the oscillations of the inflaton field amplify the inhomogeneities in the energy density, leading to local symmetry restoration and faster thermalization. We believe this explains why tachyonic preheating is so effective in transferring energy away from the inflaton zero mode.

The baryon asymmetry of the universe may originate in the phase transition at the end of hybrid inflation, provided that the reheat temperature is low enough. I show that if the field that triggers the end of inflation is the electroweak Higgs field and CP is violated, the transition leads to baryon asymmetry even if no preheating or non-thermal symmetry restoration takes place. I estimate the strength of this effect and the constraints it imposes on the inflationary model.

We consider the construction of inverted hybrid inflation models in which the vacuum energy during inflation is at the TeV scale, and the inflaton couples to the Higgs field. Such models are of interest in the context of some recently proposed models of electroweak baryogenesis. We demonstrate how constraints on these models arise from quantum corrections, and how self-consistent examples may be constructed, albeit at the expense of fine-tuning. We discuss two possible ways in which the baryon asymmetry of the universe may be produced in these models. One of them is based on preheating and a consequent non-thermal electroweak symmetry restoration and the other on the formation of Higgs winding configurations by the Kibble mechanism at the end of inflation.

We show by studying the Abelian Higgs model with numerical lattice simulations that non-thermal phase transitions arising out of preheating after inflation are possible in gauge-Higgs models under rather general circumstances. This may lead to the formation of gauged topological defects and, if the scale at which inflation ends is low enough, to electroweak baryogenesis after preheating.

Scalar electrodynamics can be used to investigate the formation of cosmic strings in the early universe. We present the results of lattice Monte Carlo simulations of an effective three-dimensional U(1)+Higgs theory that describes the equilibrium properties of finite-temperature scalar electrodynamics near the transition. A gauge-invariant criterion for the existence of a vortex is used in measuring the properties of the vortex network in the equilibrium state both in the Coulomb and in the Higgs phase of the system. The naive definition of the vortex density becomes meaningless in the continuum limit and special care is needed in extracting physical quantities. Numerical evidence for a physical discontinuity in the vortex density is given.