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.

$SU(2)$ gauge fields and axions can have a stable, isotropic and homogeneous configuration during inflation. However, couplings to other matter species lead to particle production, which in turn induces backreaction on and destabilization of the non-abelian and axion background. In this paper, we first study the particle production by a $SU(2)$ gauge field coupled to a massive Dirac doublet. To carry out this calculation we have made two technical improvements compared to what has been done in the literature. First, we apply the anti-symmetrization of the operators to treat particles and anti-particles on equal footing, second, to deal with the UV divergences, we apply instantaneous subtraction. We find that, the backreaction of produced fermions on the $SU(2)$ background is negligible for model parameters of observational interest. Next, we consider production of fermions due to coupling to the axion. The tree-level backreaction on the gauge fields, as well as on the axion, is vanishingly small. We also provide an estimate for the loop effects.

In this paper, we examine the viability of inflation models with a spectator axion field coupled to both gravitational and SU(2) gauge fields via Chern-Simons couplings. Requiring phenomenological success of the axion-SU(2) sector constrains the coupling strength of the gravitational Chern-Simons term. We find that the impact of this term on the production and propagation of gravitational waves can be as large as fifty percent enhancement for the helicity that is not sourced by the gauge field, if the cut-off scale is as low as $\Lambda$ = 20H. The effect becomes smaller for a larger value of $\Lambda$, while the impact on the helicity sourced by the gauge field is negligible regardless of $\Lambda$.

We investigate the potential of the galaxy power spectrum to constrain compensated isocurvature perturbations (CIPs), primordial fluctuations in the baryon density that are compensated by fluctuations in CDM density to ensure an unperturbed total matter density. We show that CIPs contribute to the galaxy overdensity at linear order, and if they are close to scale-invariant, their effects are nearly perfectly degenerate with the local PNG parameter $f_{\rm nl}$ if they correlate with the adiabatic perturbations. This degeneracy can however be broken by analyzing multiple galaxy samples with different bias parameters, or by taking CMB priors on $f_{\rm nl}$ into account. Parametrizing the amplitude of the CIP power spectrum as $P_{\sigma\sigma} = A_{\rm CIP}^2P_{\mathcal{R}\mathcal{R}}$ (where $P_{\mathcal{R}\mathcal{R}}$ is the adiabatic power spectrum) we find, for a number of fiducial galaxy samples in a simplified forecast setup, that constraints on $A_{\rm CIP}$, relative to those on $f_{\rm nl}$, of order $\sigma_{A_{\rm CIP}}/\sigma_{f_{\rm nl}} \approx 1-2$ are achievable for CIPs correlated with adiabatic perturbations, and $\sigma_{A_{\rm CIP}}/\sigma_{f_{\rm nl}} \approx 5$ for the uncorrelated case. These values are independent of survey volume, and suggest that current galaxy data are already able to improve significantly on the tightest existing constraints on CIPs from the CMB. Future galaxy surveys that aim to achieve $\sigma_{f_{\rm nl}} \sim 1$ have the potential to place even stronger bounds on CIPs.

We present a new numerical algorithm and code, ${\sf GFiRe}$, for solving the non-linear evolution of Abelian gauge fields coupled to complex scalar fields in homogeneous and isotropic spacetimes. We adopt a hybrid approach to solving the system: the spatial derivatives are discretized using standard Lattice Gauge Field Theory techniques, whereas the time evolution of the fields and scale factor is implemented with explicit, composite, symplectic integrators. An important property of our compound algorithm is that the discretized Gauss constraint is respected exactly, regardless of the order of the symplectic integrator. This remains true even when the background expansion is computed "self-consistently"; that is, when the expansion history is computed using spatial averaged components of the energy momentum tensor in the simulation volume. Hence, our code can also be used in cases where the fields dominate the energy density of the universe, for example, during reheating after inflation. We test the algorithm in scenarios of reheating where the inflaton is a complex scalar field with a potential $\propto(2|\varphi|^2-v^2)^2$ and is coupled to an Abelian gauge field. Tracing the evolution of the system through complex dynamics (including resonant excitation of fields, backreaction, formation of solitons, and changes in the equation of state) in a self-consistently expanding universe, we find the energy conservation violation ($<10^{-4}$) to be very stable and the Gauss constraint violation ($<10^{-6}$) to be dominated by differencing noise.

This lecture note provides a generic introduction to reheating after inflation, covering theoretical, phenomenological and observational aspects of the process.

We study the scalar and tensor perturbations generated by the fragmentation of the inflaton condensate into oscillons or transients after inflation, using nonlinear classical lattice simulations. Without including the backreaction of metric perturbations, we find that the magnitude of scalar metric perturbations never exceeds a few $\times 10^{-3}$, whereas the maximal strength of the gravitational wave signal today is $\mathcal{O}(10^{-9})$ for standard post-inflationary expansion histories. We provide parameter scalings for the $\alpha$-attractor models of inflation, which can be easily applied to other models. We also discuss the likelihood of primordial black hole formation, as well as conditions under which the gravitational wave signal can be at observationally interesting frequencies and amplitudes. Finally, we provide an upper bound on the frequency of the peak of the gravitational wave signal, which applies to all preheating scenarios.

The Higgs potential appears to be fine-tuned, hence very sensitive to values of other scalar fields that couple to the Higgs. We show that this feature can lead to a new epoch in the early universe featuring violent dynamics coupling the Higgs to a scalar modulus. The oscillating modulus drives tachyonic Higgs particle production. We find a simple parametric understanding of when this process can lead to rapid modulus fragmentation, resulting in gravitational wave production. A nontrivial equation-of-state arising from the nonlinear dynamics also affects the time elapsed from inflation to the CMB, influencing fits of inflationary models. Supersymmetric theories automatically contain useful ingredients for this picture.

Non-Abelian gauge fields may exist during inflation. We study the Schwinger effect by an $SU(2)$ gauge field coupled to a charged scalar doublet in a (quasi) de Sitter background and the possible backreaction of the generated charged particles on the homogeneous dynamics. Contrary to the Abelian $U(1)$ case, we find that both the Schwinger pair production and the induced current decrease as the interaction strength increases. The reason for this suppression is the isotropic vacuum expectation value of the $SU(2)$ field which generates a (three times) greater effective mass for the scalar field than the $U(1)$. In the weak interaction limit, the above effect is negligible and both the $SU(2)$ and $U(1)$ cases exhibit a linear increase of the current and a constant conductivity with the interaction strength. We conclude that the Schwinger effect does not pose a threat to the dynamics of inflationary models involving an $SU(2)$ gauge field.

Homogeneous oscillations of the inflaton after inflation can be unstable to small spatial perturbations even without coupling to other fields. We show that for inflaton potentials $\propto |\phi|^{2n}$ near $|\phi|=0$ and flatter beyond some $|\phi|=M$, the inflaton condensate oscillations can lead to self-resonance, followed by its complete fragmentation. We find that for non-quadratic minima ($n>1$), shortly after backreaction, the equation of state parameter, $w\rightarrow1/3$. If $M\ll m_{pl}$, radiation domination is established within less than an e-fold of expansion after the end of inflation. In this case self-resonance is efficient and the condensate fragments into transient, localised spherical objects which are unstable and decay, leaving behind them a virialized field with mean kinetic and gradient energies much greater than the potential energy. This end-state yields $w=1/3$. When $M\ll m_{pl}$ we observe slow and steady, self-resonace that can last many {\it e}-folds before backreaction eventually shuts it off, followed by fragmentation and $w\rightarrow 1/3$. We provide analytical estimates for the duration to $w\rightarrow 1/3$ after inflation, which can be used as an upper bound (under certain assumptions) on the duration of the transition between the inflationary and the radiation dominated states of expansion. This upper bound can reduce uncertainties in CMB observables such as the spectral tilt $n_{\rm{s}}$, and the tensor-to-scalar ratio $r$. For quadratic minima ($n=1$), $w\rightarrow0$ regardless of the value of $M$.

We calculate the equation of state after inflation and provide an upper bound on the duration before radiation domination by taking the nonlinear dynamics of the fragmented inflaton field into account. A broad class of single-field inflationary models with observationally consistent flattening of the potential at a scale $M$ away from the origin, $V(\phi)\propto |\phi|^{2n}$ near the origin, and where the couplings to other fields are ignored are included in our analysis. We find that the equation of state parameter $w\rightarrow 0$ for $n=1$ and $w\rightarrow 1/3$ (after sufficient time) for $n\gtrsim 1$. We calculate how the number of $e$-folds to radiation domination depends on both $n$ and $M$ when $M\sim m_{\rm pl}$, whereas when $M\ll m_{\rm pl}$, we find that the duration to radiation domination is negligible. Our results are explained in terms of a linear instability analysis in an expanding universe, scaling arguments, and are supported by detailed 3+1 dimensional lattice simulations. We show how our work significantly reduces the uncertainty in inflationary observables, even after including couplings to additional light fields.

We calculate particle production during inflation and in the early stages of reheating after inflation in models with a charged scalar field coupled to Abelian and non-Abelian gauge fields. A detailed analysis of the power spectra of primordial electric fields, magnetic fields and charge fluctuations at the end of inflation and preheating is provided. We carefully account for the Gauss constraints during inflation and preheating, and clarify the role of the longitudinal components of the electric field. We calculate the timescale for the back-reaction of the produced gauge fields on the inflaton condensate, marking the onset of non-linear evolution of the fields. We provide a prescription for initial conditions for lattice simulations necessary to capture the subsequent nonlinear dynamics. On the observational side, we find that the primordial magnetic fields generated are too small to explain the origin of magnetic fields on galactic scales and the charge fluctuations are well within observational bounds for the models considered in this paper.

The dynamics at the end of inflation can generate an asymmetry between particles and anti-particles of the inflaton field. This asymmetry can be transferred to baryons via decays, generating a baryon asymmetry in our Universe. We explore this idea in detail for a complex inflaton governed by an observationally consistent -"flatter than quadratic"- potential with a weakly broken global U(1) symmetry. We find that most of the inflaton asymmetry is locked in non-topological soliton like configurations (oscillons) produced copiously at the end of inflation. These solitons eventually decay into baryons and generate the observed matter-antimatter asymmetry for a range of model parameters. Through a combination of three dimensional lattice simulations and a detailed linearized analysis, we show how the inflaton asymmetry depends on the fragmentation, the magnitude of the symmetry breaking term and initial conditions at the end of inflation. We discuss the final decay into baryons, but leave a detailed analysis of the inhomogeneous annihilation, reheating and thermalization to future work. As part of our work, we pay particular attention to generating multifield initial conditions for the field fluctuations (including metric perturbations) at the end of inflation for lattice simulations.

Results are presented from a 500ks long XMM-Newton observation of the Narrow-Line Seyfert 1 galaxy IRAS13224-3809. The source is rapidly variable on timescales down to a few 100s. The spectrum shows strong broad Fe-K and L emission features which are interpreted as arising from reflection from the inner parts of an accretion disc around a rapidly spinning black hole. Assuming a power-law emissivity for the reflected flux and that the innermost radius corresponds to the innermost stable circular orbit, the black hole spin is measured to be 0.988 with a statistical precision better than one per cent. Systematic uncertainties are discussed. A soft X-ray lag of 100s confirms this scenario. The bulk of the power-law continuum source is located at a radius of 2-3 gravitational radii.