Present and future of CosmoLattice.

We discuss the present state and planned updates of CosmoLattice, a cutting-edge code for lattice simulations of non-linear dynamics of scalar-gauge field theories in an expanding background. We first review current capabilities of the code, including the simulation of interacting singlet scalars and of Abelian and non-Abelian scalar-gauge theories. We also comment on new features recently implemented, such as the simulation of gravitational waves from scalar and gauge fields. Secondly, we discuss new extensions of CosmoLattice that we plan to release publicly. On the one hand, we comment on new physics modules, which include axion-gauge interactions φFF̃, non-minimal gravitational couplings φ^2R, creation and evolution of cosmic defect networks, and magneto-hydro-dynamics (MHD). On the other hand, we discuss new technical features, including evolvers for non-canonical interactions, arbitrary initial conditions, simulations in 2+1 dimensions, and higher accuracy spatial derivatives. .

Cosmological Background Interpretation of Pulsar Timing Array Data.

We discuss the interpretation of the detected signal by pulsar timing array (PTA) observations as a gravitational wave background of cosmological origin. We combine NANOGrav 15-years and EPTA-DR2new datasets and confront them against backgrounds from supermassive black hole binaries (SMBHBs), and cosmological signals from inflation, cosmic (super)strings, first-order phase transitions, Gaussian and non-Gaussian large scalar fluctuations, and audible axions. We find that scalar-induced, and to a lesser extent audible axion and cosmic superstring signals, provide a better fit than SMBHBs. These results depend, however, on modeling assumptions, so further data and analysis are needed to reach robust conclusions. Independently of the signal origin, the data strongly constrain the parameter space of cosmological signals, for example, setting an upper bound on primordial non-Gaussianity at PTA scales as |f_{nl}|≲2.34 at 95% C.L.

Strong Backreaction Regime in Axion Inflation.

We study the nonlinear dynamics of axion inflation, capturing for the first time the inhomogeneity and full dynamical range during strong backreaction, till the end of inflation. Accounting for inhomogeneous effects leads to a number of new relevant results, compared to spatially homogeneous studies: (i) the number of extra efoldings beyond slow-roll inflation increases very rapidly with the coupling, (ii) oscillations of the inflaton velocity are attenuated, (iii) the tachyonic gauge field helicity spectrum is smoothed out (i.e., the spectral oscillatory features disappear), broadened, and shifted to smaller scales, and (iv) the nontachyonic helicity is excited, reducing the chiral asymmetry, now scale dependent. Our results are expected to impact strongly on the phenomenology and observability of axion inflation, including gravitational wave generation and primordial black hole production.

Non-Gaussian Tail of the Curvature Perturbation in Stochastic Ultraslow-Roll Inflation: Implications for Primordial Black Hole Production.

We consider quantum diffusion in ultraslow-roll (USR) inflation. Using the ΔN formalism, we present the first stochastic calculation of the probability distribution P(R) of the curvature perturbation during USR. We capture the nonlinearity of the system, solving the coupled evolution of the coarse-grained background with random kicks from the short wavelength modes, simultaneously with the mode evolution around the stochastic background. This leads to a non-Markovian process from which we determine the highly non-Gaussian tail of P(R). Studying the production of primordial black holes in a viable model, we find that stochastic effects during USR increase their abundance by a factor of ∼10^{5} compared with the Gaussian approximation.

Anisotropies in the gravitational wave background from preheating.

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 χ coupled to the inflaton φ, with the only requisite that χ 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.

Exact scale-invariant background of gravitational waves from cosmic defects.

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

Stochastic background of gravitational waves from hybrid preheating.

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 bubblelike 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. On the other hand, low-scale models could still produce a detectable stochastic background at frequencies accessible to those detectors. The discovery of such a background would open a new window into the very early Universe.