Abstract: The Laser Interferometer Space Antenna (LISA) has two scientific objectives of cosmological focus: to probe the expansion rate of the universe, and to understand stochastic gravitational-wave backgrounds and their implications for early universe and particle physics, from the MeV to the Planck scale. However, the range of potential cosmological applications of gravitational-wave observations extends well beyond these two objectives. This publication presents a summary of the state of the art in LISA cosmology, theory and methods, and identifies new opportunities to use gravitational-wave observations by LISA to probe the universe.

Abstract: We study the impact of stochastic noise on the generation of primordial black hole (PBH) seeds in ultra-slow-roll (USR) inflation with numerical simulations. We consider the non-linearity of the system by consistently taking into account the noise dependence on the inflaton perturbations, while evolving the perturbations on the coarse-grained background affected by the noise. We capture in this way the non-Markovian nature of the dynamics, and demonstrate that non-Markovian effects are subleading. Using the Delta N formalism, we find the probability distribution P(R) of the comoving curvature perturbation R. We consider inflationary potentials that fit the CMB and lead to PBH dark matter with i) asteroid, ii) solar, or iii) Planck mass, as well as iv) PBHs that form the seeds of supermassive black holes. We find that stochastic effects enhance the PBH abundance by a factor of O(10)-O(10(8)), depending on the PBH mass. We also show that the usual approximation, where stochastic kicks depend only on the Hubble rate, either underestimates or overestimates the abundance by orders of magnitude, depending on the potential. We evaluate the gauge dependence of the results, discuss the quantum-to-classical transition, and highlight open issues of the application of the stochastic formalism to USR inflation.

Abstract: 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(GW) proportional to f(3) for f << f(0), Omega(GW) proportional to 1/f(2) for f(0) less than or similar to f less than or similar to feq, and Omega(GW) proportional to const (i.e., exactly scale invariant) for f >> f(eq), where f(0) and f(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) -> 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 cosmic microwave background upper bounds on the defects' energy scale, we discuss the difficulty of detecting this GW background in the case of global defects.

Abstract: This paper describes CosmoGattice, a modern package for lattice simulations of the dynamics of interacting scalar and gauge fields in an expanding universe. CosmoGattice incorporates a series of features that makes it very versatile and powerful: i) it is written in C++ fully exploiting the object oriented programming paradigm, with a modular structure and a clear separation between the physics and the technical details, ii) it is MPI-based and uses a discrete Fourier transform parallelized in multiple spatial dimensions, which makes it specially appropriate for probing scenarios with well -separated scales, running very high resolution simulations, or simply very long ones, iii) it introduces its own symbolic language, defining field variables and operations over them, so that one can introduce differential equations and operators in a manner as close as possible to the continuum, iv) it includes a library of numerical algorithms, ranging from O(delta t(2)) to O(delta t(10)) methods, suitable for simulating global and gauge theories in an expanding grid, including the case of 'self-consistent' expansion sourced by the fields themselves. Relevant observables are provided for each algorithm (e.g. energy densities, field spectra, lattice snapshots) and we note that, remarkably, all our algorithms for gauge theories (Abelian or non-Abelian) always respect the Gauss constraint to machine precision. Program summary Program Title:: CosmoGattice CPC Library link to program files: https://doi .org /10 .17632 /44vr5xssc6 .1 Developer's repository link: http://github .com /cosmolattice /cosmolattice Licensing provisions: MIT Programming language: C++, MPI Nature of problem: The phenomenology of high energy physics in the early universe is typically characterized by non-linear dynamics, which cannot be captured accurately with analytical techniques. In order to fully understand the non-linearities developed in a given scenario, one needs to carry out lattice simulations. A number of public packages for lattice simulations have appeared over the years, but most of them are only capable of simulating scalar fields. However, realistic models of particle physics do contain other kind of field species, such as (Abelian or non-Abelian) gauge fields, whose non-linear dynamics can also play a relevant role in the early universe. Tensor modes representing gravitational waves are also naturally expected in many scenarios. Solution method: CosmoGattice represents a modern code for lattice simulations of scalar-gauge field theories in an expanding universe. It allows for the simulation of the evolution of interacting (singlet) scalar fields, charged scalar fields under U(1) and/or SU(2) gauge groups, and the corresponding associated Abelian and/or non-Abelian gauge fields. From version 1.1 onward, CosmoGattice also allows to simulate the production of gravitational waves. Simulations can be done either in a flat space-time background, or in a homogeneous and isotropic (spatially flat) expanding FLRW background. CosmoGattice provides symplectic integrators, with accuracy ranging from O (delta t(2)) up to O(delta t(10)), to simuate the non-linear dynamics of the appropriate fields in comoving three-dimensional lattices. The code is parallelized with MPI, and uses a discrete Fourier Transform parallelized in multiple spatial dimensions, which makes it a very powerful code for probing physical problems with well-separated scales. Moreover, the code has been designed as a `platform' to implement any system of dynamical equations suitable for discretization on a lattice.

Abstract: We discuss the possibility to measure particle couplings with stochastic gravitational wave backgrounds (SGWBs). Under certain circumstances a sequence of peaks of different amplitude and frequency-a stairway-emerges in a SGWB spectrum, with each peak probing a different coupling. The detection of such signature opens the possibility to reconstruct couplings (spectroscopy) of particle species involved in high energy phenomena generating SGWBs. Stairwaylike signatures may arise in causally produced backgrounds in the early Universe, e.g., from preheating or first order phase transitions. As a proof of principle we study a preheating scenario with an inflaton 0 coupled to multiple daughter fields f chi jg with different coupling strengths. As a clear stairway signature is imprinted in the SGWB spectrum, we reconstruct the relevant couplings with various detectors.