|
Auclair, P., Blanco-Pillado, J. J., Figueroa, D. G., Jenkins, A. C., Lewicki, M., Sakellariadou, M., et al. (2020). Probing the gravitational wave background from cosmic strings with LISA. J. Cosmol. Astropart. Phys., 04(4), 034–50pp.
Abstract: Cosmic string networks offer one of the best prospects for detection of cosmological gravitational waves (GWs). The combined incoherent GW emission of a large number of string loops leads to a stochastic GW background (SGWB), which encodes the properties of the string network. In this paper we analyze the ability of the Laser Interferometer Space Antenna (LISA) to measure this background, considering leading models of the string networks. We find that LISA will be able to probe cosmic strings with tensions G μgreater than or similar to O(10(-17)), improving by about 6 orders of magnitude current pulsar timing arrays (PTA) constraints, and potentially 3 orders of magnitude with respect to expected constraints from next generation PTA observatories. We include in our analysis possible modifications of the SGWB spectrum due to different hypotheses regarding cosmic history and the underlying physics of the string network. These include possible modifications in the SGWB spectrum due to changes in the number of relativistic degrees of freedom in the early Universe, the presence of a non-standard equation of state before the onset of radiation domination, or changes to the network dynamics due to a string inter-commutation probability less than unity. In the event of a detection, LISA's frequency band is well-positioned to probe such cosmic events. Our results constitute a thorough exploration of the cosmic string science that will be accessible to LISA.
|
|
|
Figueroa, D. G., Hindmarsh, M., Lizarraga, J., & Urrestilla, J. (2020). Irreducible background of gravitational waves from a cosmic defect network: Update and comparison of numerical techniques. Phys. Rev. D, 102(10), 103516–25pp.
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.
|
|
|
Antusch, S., Figueroa, D. G., Marschall, K., & Torrenti, F. (2020). Energy distribution and equation of state of the early Universe: Matching the end of inflation and the onset of radiation domination. Phys. Lett. B, 811, 135888–7pp.
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.
|
|
|
Figueroa, D. G., Florio, A., Torrenti, F., & Valkenburg, W. (2021). The art of simulating the early universe. Part I. Integration techniques and canonical cases. J. Cosmol. Astropart. Phys., 04(4), 035–108pp.
Abstract: We present a comprehensive discussion on lattice techniques for the simulation of scalar and gauge field dynamics in an expanding universe. After reviewing the continuum formulation of scalar and gauge field interactions in Minkowski and FLRW backgrounds, we introduce the basic tools for the discretization of field theories, including lattice gauge invariant techniques. Following, we discuss and classify numerical algorithms, ranging from methods of O(delta t(2)) accuracy like staggered leapfrog and Verlet integration, to Runge-Kutta methods up to O(delta t(4)) accuracy, and the Yoshida and Gauss-Legendre higher-order integrators, accurate up to O(delta t(10)) We adapt these methods for their use in classical lattice simulations of the non-linear dynamics of scalar and gauge fields in an expanding grid in 3+1 dimensions, including the case of 'self-consistent' expansion sourced by the volume average of the fields' energy and pressure densities. We present lattice formulations of canonical cases of: i) Interacting scalar fields, ii) Abelian U(1) gauge theories, and iii) Non-Abelian SU(2) gauge theories. In all three cases we provide symplectic integrators, with accuracy ranging from O(delta t(2)) up to O(delta t(10)) For each algorithm we provide the form of relevant observables, such as energy density components, field spectra and the Hubble constraint. We note that all our algorithms for gauge theories always respect the Gauss constraint to machine precision, including when 'self-consistent' expansion is considered. As a numerical example we analyze the post-inflationary dynamics of an oscillating inflaton charged under SU(2) x U(1). We note that the present manuscript is meant to be part of the theoretical basis for the code CosmoLattice, a multi-purpose MPI-based package for simulating the non-linear evolution of field theories in an expanding universe, publicly available at http://www.cosrnolattice.net.
|
|
|
Figueroa, D. G., Raatikainen, S., Rasanen, S., & Tomberg, E. (2021). Non-Gaussian Tail of the Curvature Perturbation in Stochastic Ultraslow-Roll Inflation: Implications for Primordial Black Hole Production. Phys. Rev. Lett., 127(10), 101302–7pp.
Abstract: We consider quantum diffusion in ultraslow-roll (USR) inflation. Using the Delta 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 similar to 10(5) compared with the Gaussian approximation.
|
|
|
Aggarwal, N. et al, & Figueroa, D. G. (2021). Challenges and opportunities of gravitational-wave searches at MHz to GHz frequencies. Living Rev. Relativ., 24(1), 4–74pp.
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.
|
|
|
Antusch, S., Figueroa, D. G., Marschall, K., & Torrenti, F. (2022). Characterizing the postinflationary reheating history: Single daughter field with quadratic-quadratic interaction. Phys. Rev. D, 105(4), 043532–36pp.
Abstract: We study the evolution of the energy distribution and equation of state of the Universe from the end of inflation until the onset of either radiation domination (RD) or a transient period of matter domination (MD). We use both analytical techniques and lattice simulations. We consider two-field models where the inflaton (/) has a monomial potential after inflation V((/)) proportional to i(/) – vip (p 4, and of order similar to 50% for p 4. The system goes to MD at late times for p = 2, while it goes to RD for p > 2. In the later case, we can calculate exactly the number of e-folds until RD as a function of g2, and hence predict accurately inflationary observables like the scalar tilt ns and the tensor-to-scalar ratio r. In the scenario (ii), the energy is always transferred completely to X for p > 2, as long as its effective mass m2X = g2((/) – v)2 is not negligible. For p = 2, the final ratio between the energy densities of X and (/) depends strongly on g2. For all p > 2, the system always goes to MD at late times.
|
|
|
Figueroa, D. G., Raatikainen, S., Rasanen, S., & Tomberg, E. (2022). Implications of stochastic effects for primordial black hole production in ultra-slow-roll inflation. J. Cosmol. Astropart. Phys., 05(5), 027–48pp.
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.
|
|
|
Figueroa, D. G., Florio, A., Loayza, N., & Pieroni, M. (2022). Spectroscopy of particle couplings with gravitational waves. Phys. Rev. D, 106(6), 063522–8pp.
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.
|
|
|
LISA Cosmology Working Group(Bartolo, N. et al), & Figueroa, D. G. (2022). Probing anisotropies of the Stochastic Gravitational Wave Background with LISA. J. Cosmol. Astropart. Phys., 11, 009–65pp.
Abstract: We investigate the sensitivity of the Laser Interferometer Space Antenna (LISA) to the anisotropies of the Stochastic Gravitational Wave Background (SGWB). We first discuss the main astrophysical and cosmological sources of SGWB which are characterized by anisotropies in the GW energy density, and we build a Signal-to-Noise estimator to quantify the sensitivity of LISA to different multipoles. We then perform a Fisher matrix analysis of the prospects of detectability of anisotropic features with LISA for individual multipoles, focusing on a SGWB with a power-law frequency profile. We compute the noise angular spectrum taking into account the specific scan strategy of the LISA detector. We analyze the case of the kinematic dipole and quadrupole generated by Doppler boosting an isotropic SGWB. We find that beta Omega(GW) similar to 2 x 10(-11) is required to observe a dipolar signal with LISA. The detector response to the quadrupole has a factor similar to 10(3) beta relative to that of the dipole. The characterization of the anisotropies, both from a theoretical perspective and from a map-making point of view, allows us to extract information that can be used to understand the origin of the SGWB, and to discriminate among distinct superimposed SGWB sources.
|
|