Abstract: The first direct measurement of gravitational waves by the LIGO and Virgo collaborations has opened up new avenues to explore our Universe. This white paper outlines the challenges and gains expected in gravitational-wave searches at frequencies above the LIGO/Virgo band, with a particular focus on Ultra High-Frequency Gravitational Waves (UHF-GWs), covering the MHz to GHz range. The absence of known astrophysical sources in this frequency range provides a unique opportunity to discover physics beyond the Standard Model operating both in the early and late Universe, and we highlight some of the most promising gravitational sources. We review several detector concepts that have been proposed to take up this challenge, and compare their expected sensitivity with the signal strength predicted in various models. This report is the summary of the workshop “Challenges and opportunities of high-frequency gravitational wave detection” held at ICTP Trieste, Italy in October 2019, that set up the stage for the recently launched Ultra-High-Frequency Gravitational Wave (UHF-GW) initiative.

Abstract: We study the energy distribution and equation of state of the universe between the end of inflation and the onset of radiation domination (RD), considering observationally consistent single-field inflationary scenarios, with a potential 'flattening' at large field values, and a monomial shape V(phi) proportional to vertical bar phi vertical bar(p) around the origin. As a proxy for (p)reheating, we include a quadratic interaction g(2)phi X-2(2) between the inflaton phi and a light scalar 'daughter' field X, with g(2) > 0. We capture the non-perturbative and non-linear nature of the system dynamics with lattice simulations, obtaining that: i) the final energy transferred to X depends only on p, not on g(2); ii) the final transfer of energy is always negligible for 2 <= p < 4, and of order similar to 50% for p >= 4; iii) the system goes at late times to matter-domination for p = 2, and always to RD for p > 2. In the latter case we calculate the number of e-folds until RD, significantly reducing the uncertainty in the inflationary observables Tl-s and r.

Abstract: We use the latest constraints on the population of stellar origin binary black holes (SOBBH) from LIGO/Virgo/KAGRA (LVK) observations, to estimate the stochastic gravitational wave background (SGWB) they generate in the frequency band of LISA. In order to account for the faint and distant binaries, which contribute the most to the SGWB, we extend the merger rate at high redshift assuming that it tracks the star formation rate. We adopt different methods to compute the SGWB signal: we perform an analytical evaluation, we use Monte Carlo sums over the SOBBH population realisations, and we account for the role of the detector by simulating LISA data and iteratively removing the resolvable signals until only the confusion noise is left. The last method allows the extraction of both the expected SGWB and the number of resolvable SOBBHs. Since the latter are few for signal-to-noise ratio thresholds larger than five, we confirm that the spectral shape of the SGWB in the LISA band agrees with the analytical prediction of a single power law. We infer the probability distribution of the SGWB amplitude from the LVK GWTC-3 posterior of the binary population model: at the reference frequency of 0.003 Hz it has an interquartile range of h2ΩGW(f = 3 × 10-3 Hz) ∈ [5.65, 11.5] × 10-13, in agreement with most previous estimates. We then perform a MC analysis to assess LISA's capability to detect and characterise this signal. Accounting for both the instrumental noise and the galactic binaries foreground, with four years of data, LISA will be able to detect the SOBBH SGWB with percent accuracy, narrowing down the uncertainty on the amplitude by one order of magnitude with respect to the range of possible amplitudes inferred from the population model. A measurement of this signal by LISA will help to break the degeneracy among some of the population parameters, and provide interesting constraints, in particular on the redshift evolution of the SOBBH merger rate.

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.

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.