We investigate the importance of local gravity during preheating, the non-linear dynamics that may be responsible for starting the process of reheating the universe after inflation. We introduce three numerical methods that study a simple preheating scenario while relaxing gravitational assumptions, culminating in studying the process in full numerical relativity. We confirm that perturbation theory is no longer valid when one considers modes whose wavelengths are comparable to the size of the horizon at the end of inflation; however, this breakdown does not necessarily lead to a breakdown of the preheating process in non-linear gravity. For the specific model we test we find no evidence for the creation of primordial black holes from the instabilities in this model. Finally, we remark on the opportunity for future numerical study of non-linear gravitational dynamics in the early universe.

We demonstrate that gravitational waves generated by efficient gauge preheating after axion inflation generically contribute significantly to the effective number of relativistic degrees of freedom $N_\mathrm{eff}$. We show that, with existing Planck limits, gravitational waves from preheating already place the strongest constraints on the inflaton's possible axial coupling to gauge fields. We demonstrate that gauge preheating can completely reheat the Universe regardless of the inflationary potential. Further, we quantify the variation of the efficiency of gravitational wave production from model to model and show that it is correlated with the tensor-to-scalar ratio. In particular, when combined with constraints on models whose tensor-to-scalar ratios would be detected by next-generation cosmic microwave background experiments, $r\gtrsim 10^{-3}$, constraints from $N_\mathrm{eff}$ will probe or rule out the entire coupling regime for which gauge preheating is efficient.

We study gravitational wave production from gauge preheating in a variety of inflationary models, detailing its dependence on both the energy scale and the shape of the potential. We show that gauge preheating generically leads to a large gravitational wave background that contributes significantly to the effective number of relativistic degrees of freedom in the early Universe, $N_\mathrm{eff}$. We demonstrate that the efficiency of gravitational wave production is correlated with the tensor-to-scalar ratio, $r$. In particular, we show that efficient gauge preheating in models whose tensor-to-scalar ratio would be detected by next-generation cosmic microwave background experiments ($r \gtrsim 10^{-3}$) will either be detected through its contribution to $N_\mathrm{eff}$ or ruled out. Furthermore, we show that bounds on $N_\mathrm{eff}$ provide the most sensitive probe of the possible axial coupling of the inflaton to gauge fields regardless of the potential.

We examine how cosmological expansion arises in a universe containing a lattice of spinning black holes. We study averaged expansion properties as a function of fundamental properties of the black holes, including the bare mass of the black holes and black hole spin. We then explore how closely the expansion properties correspond to properties of a corresponding matter-dominated FLRW universe. As residual radiation present in the initial data decays, we find good agreement with a matter-dominated FLRW solution, and the effective density in the volume is well-described by the horizon mass of the black hole.

The next generation of instruments designed to measure the polarization of the cosmic microwave background (CMB) will provide a historic opportunity to open the gravitational wave window to the primordial Universe. Through high sensitivity searches for primordial gravitational waves, and tighter limits on the energy released in processes like phase transitions, the CMB polarization data of the next decade has the potential to transform our understanding of the laws of physics underlying the formation of the Universe.

Standard cosmological models rely on an approximate treatment of gravity, utilizing solutions of the linearized Einstein equations as well as physical approximations. In an era of precision cosmology, we should ask: are these approximate predictions sufficiently accurate for comparison to observations, and can we draw meaningful conclusions about properties of our Universe from them? In this work we examine the accuracy of linearized gravity in the presence of collisionless matter and a cosmological constant utilizing fully general relativistic simulations. We observe the gauge-dependence of corrections to linear theory, and note the amplitude of these corrections. For perturbations whose amplitudes are in line with expectations from the standard $\Lambda$CDM model, we find that the full, general relativistic metric is well-described by linear theory in Newtonian and harmonic gauges, while the metric in comoving-synchronous gauge is not. For the largest observed structures in our Universe, our results suggest that corrections to linear gravitational theory can reach or surpass the percent-level.

We develop a full four-dimensional numerical code to study scalar gravitational radiation emitted from binary systems and probe the Vainshtein mechanism in situations that break the static and spherical symmetry, relevant for binary pulsars as well as black holes and neutron stars binaries. The present study focuses on the cubic Galileon which arises as the decoupling limit of massive theories of gravity. Limitations associated with the numerical methods prevent us from reaching a physically realistic hierarchy of scales; nevertheless, within this context we observe the same power law scaling of the radiated power as previous analytic estimates, and confirm a strong suppression of the power emitted in the monopole and dipole as compared with quadrupole radiation. Following the trend to more physically realistic parameters, we confirm the suppression of the power emitted in scalar gravitational radiation and the recovery of General Relativity with good accuracy. This paves the way for future numerical work, probing more generic, physically relevant situations and sets of interactions that may exhibit the Vainshtein mechanism.

We study gravitational wave production during Abelian gauge-field preheating following inflation. We consider both scalar and pseudoscalar inflaton models coupled directly to Abelian gauge fields via either a dilatonic coupling to the gauge-field kinetic term or an axial coupling to a Chern-Simons term. In both cases gravitational waves are produced efficiently during the preheating phase, with a signature louder than most cosmological signals. These gravitational waves can contribute to the radiation energy budget of Universe at a level which will be probed by upcoming cosmic microwave background experiments through $N_{\rm eff}$. For axially coupled fields the resulting gravitational wave spectrum is helically polarized---a unique feature that can be used to differentiate it from other stochastic gravitational wave backgrounds. We compute the gravitational topological charge and demonstrate that gauge preheating following axion inflation may be responsible for the matter-antimatter asymmetry of the Universe via gravitational leptogenesis.

It has been recently suggested that oscillons produced in the early universe from certain asymmetric potentials continue to emit gravitational waves for a number of $e$-folds of expansion after their formation, leading to potentially detectable gravitational wave signals. We revisit this claim by conducting a convergence study using graphics processing unit (GPU)-accelerated lattice simulations and show that numerical errors accumulated with time are significant in low-resolution scenarios, or in scenarios where the run-time causes the resolution to drop below the relevant scales in the problem. Our study determines that the dominant, growing high frequency peak of the gravitational wave signals in the fiducial "hill-top model" in [arXiv:1607.01314] is a numerical artifact. This finding prompts the need for a more careful analysis of the numerical validity of other similar results related to gravitational waves from oscillon dynamics.

We study preheating in models where a scalar inflaton is directly coupled to a non-Abelian $SU(2)$ gauge field. In particular, we examine $m^2 \phi^2$ inflation with a conformal, dilaton-like coupling to the non-Abelian sector. We describe a numerical scheme that combines lattice gauge theory with standard finite difference methods applied to the scalar field. We show that a significant tachyonic instability allows for efficient preheating, which is parametrically suppressed by increasing the non-Abelian self-coupling. Additionally, we comment on the technical implementation of the evolution scheme and setting initial conditions.

We compute the weak lensing convergence power spectrum, $C^{\kappa\kappa}(\theta)$, in a dust-filled universe using fully non-linear general relativistic simulations. The spectrum is then compared to more standard, approximate calculations by computing the Bardeen (Newtonian) potentials in linearized gravity and utilizing the Born approximation. We find corrections to the angular power spectrum amplitude of order ten percent at very large angular scales, $\ell ~ 2-3$, and percent-level corrections at intermediate angular scales of $\ell ~ 20-30$.

We investigate whether Effective Field Theory (EFT) approaches, which have been useful in examining inflation and dark energy, can also be used to establish a systematic approach to inflationary reheating. We consider two methods. First, we extend Weinberg's background EFT to the end of inflation and reheating. We establish when parametric resonance and decay of the inflaton occurs, but also find intrinsic theoretical limitations, which make it difficult to capture some reheating models. This motivates us to next consider Cheung, et. al.'s EFT approach, which instead focuses on perturbations and the symmetry breaking induced by the cosmological background. Adapting the latter approach to reheating implies some new and important differences compared to the EFT of Inflation. In particular, there are new hierarchical scales, and we must account for inflaton oscillations during reheating, which lead to discrete symmetry breaking. Guided by the fundamental symmetries, we construct the EFT of reheating, and as an example of its usefulness we establish a new class of reheating models and the corresponding predictions for gravity wave observations. In this paper we primarily focus on the first stages of preheating. We conclude by discussing challenges for the approach and future directions. This paper builds on ideas first proposed in the note arXiv:1507.06651.

In this work we compute the production of magnetic fields in models of axion inflation coupled to the hypercharge sector of the Standard Model through a Chern-Simons interaction term. We make the simplest choice of a quadratic inflationary potential and use lattice simulations to calculate the magnetic field strength, helicity and correlation length at the end of inflation. For small values of the axion-gauge field coupling strength the results agree with no-backreaction calculations and estimates found in the literature. For larger couplings the helicity of the magnetic field differs from the no-backreaction estimate and depends strongly on the comoving wavenumber. We estimate the post-inflationary evolution of the magnetic field based on known results for the evolution of helical and non-helical magnetic fields. The magnetic fields produced by axion inflation with large couplings to $U(1)_Y$ can reach $B_{\rm eff} \gtrsim 10^{-16}\, {\rm G}$, exhibiting a field strength $B_{\rm phys} \approx 10^{-13}\, {\rm G}$ and a correlation length $\lambda_{\rm phys}\approx10\, {\rm pc}$. This result is insensitive to the exact value of the coupling, as long as the coupling is large enough to allow for instantaneous preheating. Depending on the assumptions for the physical processes that determine blazar properties, these fields can be found consistent with blazar observations based on the value of $B_{\rm eff}$. Finally, the intensity of the magnetic field for large coupling can be enough to satisfy the requirements for a recently proposed baryogenesis mechanism, which utilizes the chiral anomaly of the Standard Model.

We study the onset of the reheating epoch at the end of axion-driven inflation where the axion is coupled to an Abelian, $U(1)$, gauge field via a Chern-Simons interaction term. We focus primarily on $m^2\phi^2$ inflation and explore the possibility that preheating can occur for a range of coupling values consistent with recent observations and bounds on the overproduction of primordial black holes. We find that for a wide range of parameters preheating is efficient. In certain cases the inflaton is seen to transfer all its energy to the gauge fields within a few oscillations. In most cases, we find that the gauge fields on sub-horizon scales end preheating in an unpolarized state due to the existence of strong rescattering between the inflaton and gauge-field modes. We also present a preliminary study of an axion monodromy model coupled to $U(1)$ gauge fields, seeing a similarly efficient preheating behavior as well as indications that the coupling strength has an effect on the creation of oscillons.

How does inhomogeneity affect our interpretation of cosmological observations? It has long been wondered to what extent the observable properties of an inhomogeneous universe differ from those of a corresponding Friedman-Lemaitre-Robertson-Walker (FLRW) model, and how the inhomogeneities affect that correspondence. Here, we use numerical relativity to study the behavior of light beams traversing an inhomogeneous universe and construct the resulting Hubble diagrams. The universe that emerges exhibits an average FLRW behavior, but inhomogeneous structures contribute to deviations in observables across the observer's sky. We also investigate the relationship between angular diameter distance and the angular extent of a source, finding deviations that grow with source redshift. These departures from FLRW are important path-dependent effects with implications for using real observables in an inhomogeneous universe such as our own.

We present cosmological-scale numerical simulations of an evolving universe in full general relativity (GR) and introduce a new numerical tool, {\sc CosmoGRaPH}, which employs the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formalism on a 3-dimensional grid. Using {\sc CosmoGRaPH}, we calculate the effect of an inhomogeneous matter distribution on the evolution of a spacetime. We also present the results of a set of standard stability tests to demonstrate the robustness of our simulations.

While the use of numerical general relativity for modeling astrophysical phenomena and compact objects is commonplace, the application to cosmological scenarios is only just beginning. Here, we examine the expansion of a spacetime using the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formalism of numerical relativity in synchronous gauge. This work represents the first numerical cosmological study that is fully relativistic, non-linear and without symmetry. The universe that emerges exhibits an average Friedmann-Lema\"itre-Robertston-Walker (FLRW) behavior, however this universe also exhibits locally inhomogeneous expansion beyond that expected in linear perturbation theory around a FLRW background.

Quantum fluctuations of the gravitational field in the early Universe, amplified by inflation, produce a primordial gravitational-wave background across a broad frequency band. We derive constraints on the spectrum of this gravitational radiation, and hence on theories of the early Universe, by combining experiments that cover 29 orders of magnitude in frequency. These include Planck observations of cosmic microwave background temperature and polarization power spectra and lensing, together with baryon acoustic oscillations and big bang nucleosynthesis measurements, as well as new pulsar timing array and ground-based interferometer limits. While individual experiments constrain the gravitational-wave energy density in specific frequency bands, the combination of experiments allows us to constrain cosmological parameters, including the inflationary spectral index, $n_t$, and the tensor-to-scalar ratio, $r$. Results from individual experiments include the most stringent nanohertz limit of the primordial background to date from the Parkes Pulsar Timing Array, $\Omega_{\rm gw}(f)<2.3\times10^{-10}$. Observations of the cosmic microwave background alone limit the gravitational-wave spectral index at 95\% confidence to $n_t\lesssim5$ for a tensor-to-scalar ratio of $r = 0.11$. However, the combination of all the above experiments limits $n_t<0.36$. Future Advanced LIGO observations are expected to further constrain $n_t<0.34$ by 2020. When cosmic microwave background experiments detect a non-zero $r$, our results will imply even more stringent constraints on $n_t$ and hence theories of the early Universe.

The recent claim by BICEP2 of evidence for primordial gravitational waves from inflation has focused interest on the potential for early-Universe cosmology using observations of gravitational waves. In addition to cosmic microwave background detectors, efforts are underway to carry out gravitational-wave astronomy over a wide range of frequencies including pulsar timing arrays (nHz), space-based detectors (mHz), and terrestrial detectors ($\sim$10-2000 Hz). This multiband effort will probe a wide range of times in the early Universe (each corresponding to a different energy scale), during which gravitational-wave backgrounds may have been produced through processes such as phase transitions or preheating. In this letter, we derive a rule of thumb (not quite so strong as an upper limit) governing the maximum energy density of cosmological backgrounds. For most cosmological scenarios, we expect the energy density spectrum to peak at values of $\Omega_\text{gw}(f)\lesssim10^{-12\pm2}$. We discuss the applicability of this rule of thumb and the implications for gravitational-wave astronomy.

First-order phase transitions are a source of stochastic gravitational radiation. Precision calculations of the gravitational waves emitted during these processes, sourced by both the degrees of freedom undergoing the transition and anisotropic stress of the coupled, ambient constituents, have reached an age of maturity. Here we present high-resolution numerical simulations of a scalar field coupled to a fluid and parameterize the final gravitational wave spectrum as a function of the ratio of the energies of the two sectors and the coupling between the two sectors for a set of models that represent different types of first-order phase transitions. In most cases, the field sector is the dominant source of gravitational radiation, but it is possible in certain scenarios for the fluid to have the most important contribution.

We present the first 3+1-dimensional numerical simulations of scalar fields with non-minimal kinetic terms. As an example, we examine the existence and stability of preheating in the presence of a Dirac-Born-Infeld (DBI) inflaton coupled to a canonical matter field. The simulations represent the full non-linear theory in the presence of an expanding Universe. We show that parametric resonance in the matter field, along with self-resonance in the inflaton, repopulate the Universe with matter particles as efficiently as in traditional preheating.

Here we consider the possibility of preheating the Universe via the parametric amplification of a massless, U(1) abelian gauge field. We assume that the gauge field is coupled to the inflaton via a conformal factor with one free parameter. We present the results of high-resolution three-dimensional simulations of this model and show this mechanism efficiently preheats the Universe to a radiation-dominated final state.

We suggest that stable states of positronium might exist in the jets of active galactic nuclei (AGN). Electrons and positrons are created near the accretion disks of supermassive black holes at the centers of AGN and are accelerated along magnetic field lines while within the {\sl Alfv\'en radius}. The conditions in this region are ideal for the creation of bound states of positronium which are stable against annihilation. Traveling at relativistic speeds along the jet, the helical magnetic field enables the atoms to survive for great distances.

It is believed that first-order phase transitions at or around the GUT scale will produce high-frequency gravitational radiation. This radiation is a consequence of the collisions and coalescence of multiple bubbles during the transition. We employ high-resolution lattice simulations to numerically evolve a system of bubbles using only scalar fields, track the anisotropic stress during the process and evolve the metric perturbations associated with gravitational radiation. Although the radiation produced during the bubble collisions has previously been estimated, we find that the coalescence phase enhances this radiation even in the absence of a coupled fluid or turbulence. We comment on how these simulations scale and propose that the same enhancement should be found at the Electroweak scale; this modification should make direct detection of a first-order electroweak phase transition easier.

We discuss preheating after multi-field inflation in the presence of several preheat matter fields that become light in the vicinity of (but not at) the inflatons' VEV, at distinct extra-species-points (ESP); this setup is motivated by inflationary models that include particle production during inflation, e.g.trapped inflation, grazing ESP encounters or modulated trapping, among others. While de-phasing of inflatons tends to suppress parametric resonance, we find two new effects leading to efficient preheating: particle production during the first in-fall (efficient if many preheat matter fields are present) and a subsequent (narrow) resonance phase (efficient if an ESP happens to be at one of several distinct distances from the inflatons' VEV). Particles produced during the first in-fall are comprised of many species with low occupation number, while the latter are made up of a few species with high occupation number. We provide analytic descriptions of both phases in the absence of back-reaction, which we test numerically. We further perform lattice simulations to investigate the effects of back-reaction. We find resonances to be robust and the most likely cause of inflaton decay in multi-field trapped inflation if ESP distributions are dense.

We present the simplest model for classical transitions in flux vacua. A complex field with a spontaneously broken U(1) symmetry is embedded in $M_2\times S_1$. We numerically construct different winding number vacua, the vortices interpolating between them, and simulate the collisions of these vortices. We show that classical transitions are generic at large boosts, independent of whether or not vortices miss each other in the compact $S_1$.

Global second-order phase transitions are expected to produce scale-invariant gravitational wave spectra. In this manuscript we explore the dynamics of a symmetry-breaking phase transition using lattice simulations. We explicitly calculate the stochastic gravitational wave background produced during the transition and subsequent self-ordering phase. We comment on this signal as it compares to the scale-invariant spectrum produced during inflation.

We numerically compute features in the power-spectrum that originate from the decay of fields during inflation. Using a simple, phenomenological, multi-field setup, we increase the number of fields from a few to thousands. Whenever a field decays, its associated potential energy is transferred into radiation, causing a jump in the equation of state parameter and mode mixing at the perturbed level. We observe discrete steps in the power-spectrum if the number of fields is low, in agreement with analytic arguments in the literature. These features become increasingly smeared out once many fields decay within a given Hubble time. In this regime we confirm the validity of the analytic approach to staggered inflation, which is based on a coarse-graining procedure. Our numerical approach bridges the aforementioned analytic treatments, and can be used in more complicated scenarios.

Parametric resonances provide a mechanism by which particles can be created just after inflation. Thus far, attention has focused on a single or many inflaton fields coupled to a single scalar field. However, generically we expect the inflaton to couple to many other relativistic degrees of freedom present in the early universe. Using simulations in an expanding Friedmann-Lema\^itre-Robertson-Walker spacetime, in this paper we show how preheating is affected by the addition of multiple fields coupled to the inflaton. We focus our attention on gravitational wave production--an important potential observational signature of the preheating stage. We find that preheating and its gravitational wave signature is robust to the coupling of the inflaton to more matter fields.

It has recently been shown in high resolution numerical simulations that relativistic collisions of bubbles in the context of a multi-vacua potential may lead to the creation of bubbles in a new vacuum. In this paper, we show that scalar fields with only potential interactions behave like free fields during high-speed collisions; the kick received by them in a collision can be deduced simply by a linear superposition of the bubble wall profiles. This process is equivalent to the scattering of solitons in 1+1 dimensions. We deduce an expression for the field excursion (shortly after a collision), which is related simply to the field difference between the parent and bubble vacua, i.e. contrary to expectations, the excursion cannot be made arbitrarily large by raising the collision energy. There is however a minimum energy threshold for this excursion to be realized. We verify these predictions using a number of 3+1 and 1+1 numerical simulations. A rich phenomenology follows from these collision induced excursions - they provide a new mechanism for scanning the landscape, they might end/begin inflation, and they might constitute our very own big bang, leaving behind a potentially observable anisotropy.

A population of very light primordial black holes which evaporate before nucleosynthesis begins is unconstrained unless the decaying black holes leave stable relics. We show that gravitons Hawking radiated from these black holes would source a substantial stochastic background of high frequency gravititational waves ($10^{12}$ Hz or more) in the present universe. These black holes may lead to a transient period of matter dominated expansion. In this case the primordial universe could be temporarily dominated by large clusters of "Hawking stars" and the resulting gravitational wave spectrum is independent of the initial number density of primordial black holes.

Given a scalar field with metastable minima, bubbles nucleate quantum mechanically. When bubbles collide, energy stored in the bubble walls is converted into kinetic energy of the field. This kinetic energy can facilitate the classical nucleation of new bubbles in minima that lie below those of the "parent" bubbles. This process is efficient and classical, and changes the dynamics and statistics of bubble formation in models with multiple vacua, relative to that derived from quantum tunneling.

We investigate the feasibility of explosive particle production via parametric resonance or tachyonic preheating in multi-field inflationary models by means of lattice simulations. We observe a strong suppression of resonances in the presence of four-leg interactions between the inflaton fields and a scalar matter field, leading to insufficient preheating when more than two inflatons couple to the same matter field. This suppression is caused by a dephasing of the inflatons that increases the effective mass of the matter field. Including three-leg interactions leads to tachyonic preheating, which is not suppressed by an increase in the number of fields. If four-leg interactions are sub-dominant, we observe a slight enhancement of tachyonic preheating. Thus, in order for preheating after multi-field inflation to be efficient, one needs to ensure that three-leg interactions are present. If no tachyonic contributions exist, we expect the old theory of reheating to be applicable.

We consider the impact of thermal inflation -- a short, secondary period of inflation that can arise in supersymmetric scenarios -- on the stochastic gravitational wave background. We show that while the primordial inflationary gravitational wave background is essentially unchanged at CMB scales, it is massively diluted at solar system scales and would be unobservable by a BBO style experiment. Conversely, bubble collisions at the end of thermal inflation can generate a new stochastic background. We calculate the likely properties of the bubbles created during this phase transition, and show that the expected amplitude and frequency of this signal would fall within the BBO range.

Parametric resonance or preheating is a plausible mechanism for bringing about the transition between the inflationary phase and a hot, radiation dominated universe. This epoch results in the rapid production of heavy particles far from thermal equilibrium and could source a significant stochastic background of gravitational radiation. Here, we present a numerical algorithm for computing the contemporary power spectrum of gravity waves generated in this post-inflationary phase transition for a large class of scalar-field driven inflationary models. We explicitly calculate this spectrum for both quartic and quadratic models of chaotic inflation, and low-scale hybrid models. In particular, we consider hybrid models with an ``inverted'' potential. These models have a very short and intense period of resonance which is qualitatively different from previous examples studied in this context, but we find that they lead to a similar spectrum of gravitational radiation.

We consider gravitational wave production due to parametric resonance at the end of inflation, or ``preheating''. This leads to large inhomogeneities which source a stochastic background of gravitational waves at scales inside the comoving Hubble horizon at the end of inflation. We confirm that the present amplitude of these gravitational waves need not depend on the inflationary energy scale. We analyze an explicit model where the inflationary energy scale is ~10^9 GeV, yielding a signal close to the sensitivity of Advanced LIGO and BBO. This signal highlights the possibility of a new observational ``window'' into inflationary physics, and provides significant motivation for searches for stochastic backgrounds of gravitational waves in the Hz to GHz range, with an amplitude on the order of \Omega_{gw}(k)h^2 ~ 10^-11. Finally, the strategy used in our numerical computations is applicable to the gravitational waves generated by many inhomogeneous processes in the early universe.

We examine the dynamics of inflation driven by multiple, interacting scalar fields and derive a multi field version of the Hubble slow roll expansion. We show that the properties of this expansion naturally generalize those of the single field case. We present an analogous hierarchy of slow roll parameters, and derive the system of ``flow'' equations that describes their evolution, and show that when this system is truncated at finite order, it can always be solved exactly. Lastly, we express the scalar and tensor perturbation spectra in terms of the slow roll parameters.