Search

We study the generation of gravitational waves (GWs) during a first-order cosmological phase transition (PT) using the recently introduced Higgsless approach to numerically evaluate the fluid motion induced by the PT. We present for the first time spectra from strong first-order PTs ($\alpha = 0.5$), alongside weak ($\alpha = 0.0046$) and intermediate ($\alpha = 0.05$) transitions previously considered in the literature. We test the regime of applicability of the stationary source assumption, characteristic of the sound-shell model, and show that it agrees with our numerical results when the kinetic energy, sourcing GWs, does not decay with time. However, we find in general that for intermediate and strong PTs, the kinetic energy in our simulations decays following a power law in time, and provide a theoretical framework that extends the stationary assumption to one that allows to include the time evolution of the source. This decay of the kinetic energy, potentially determined by non-linear dynamics and hence, related to the production of vorticity, modifies the usually assumed linear growth with the source duration to an integral over time of the kinetic energy fraction, effectively reducing the growth rate. We validate the novel theoretical model with the results of our simulations covering a broad range of wall velocities. We provide templates for the GW amplitude and spectral shape for a broad range of PT parameters., Comment: 55 pages, 15 figures

The Lunar Gravitational-wave Antenna (LGWA) is a proposed array of next-generation inertial sensors to monitor the response of the Moon to gravitational waves (GWs). Given the size of the Moon and the expected noise produced by the lunar seismic background, the LGWA would be able to observe GWs from about 1 mHz to 1 Hz. This would make the LGWA the missing link between space-borne detectors like LISA with peak sensitivities around a few millihertz and proposed future terrestrial detectors like Einstein Telescope or Cosmic Explorer. In this article, we provide a first comprehensive analysis of the LGWA science case including its multi-messenger aspects and lunar science with LGWA data. We also describe the scientific analyses of the Moon required to plan the LGWA mission.

We develop a tool for the analysis of stochastic gravitational wave backgrounds from cosmological first-order phase transitions with LISA: we initiate a template databank for these signals, prototype their searches, and forecast their reconstruction. The templates encompass the gravitational wave signals sourced by bubble collisions, sound waves and turbulence. Accounting for Galactic and extra-Galactic foregrounds, we forecast the region of the parameter space that LISA will reconstruct with better than $\sim 10\,\%$ accuracy, if certain experimental and theoretical uncertainties are solved by the time LISA flies. We illustrate the accuracy with which LISA can reconstruct the parameters on a few benchmark signals, both in terms of the template parameters and the phase transition ones. To show the impact of the forecasts on physics beyond the Standard Model, we map the reconstructed benchmark measurements into the parameter spaces of the singlet extension of the Standard Model and of the classically conformal invariant $U(1)_{B-L}$ model., Comment: 54 pages, 18 figures; v2 matches version published in JCAP

We compute the gravitational wave (GW) spectrum sourced by sound waves produced during a first-order phase transition in the radiation-dominated epoch. The correlator of the velocity field is evaluated in accordance with the sound shell model. In our derivation we include the effects of the expansion of the Universe, which are relevant in particular for sourcing processes whose time duration is comparable with the Hubble time. Our results show a causal growth at small frequencies, $\Omega_{\rm GW} \sim k^3$, possibly followed by a linear regime $\Omega_{\rm GW} \sim k$ at intermediate $k$, depending on the phase transition parameters. Around the peak, we find a steep growth that approaches the $k^9$ scaling found within the sound shell model. The resulting bump around the peak of the GW spectrum may represent a distinctive feature of GWs produced from acoustic motion. Nothing similar has been observed for vortical (magneto)hydrodynamic turbulence. Nevertheless, we find that the $k^9$ scaling is less extended than expected in the literature, and it does not necessarily appear. The dependence on the duration of the source, $\delta \tau_{\rm fin}$, is quadratic at small frequencies $k$, and proportional to $\ln^2 (1 + \delta \tau_{\rm fin} H_*)$ for an expanding Universe. At frequencies around the peak, the growth is suppressed by a factor $\Upsilon = 1 - 1/(1 + \delta \tau_{\rm fin} {H}_*)$ that becomes linear when the GW source is short. We discuss in which cases the dependence on the source duration is linear or quadratic for stationary processes. This affects the amplitude of the GW spectrum, both in the causality tail and at the peak, showing that the assumption of stationarity is a very relevant one, as far as the GW spectral shape is concerned. Finally, we present a general semi-analytical template of the resulting GW spectrum, as a function of the parameters of the phase transition., Comment: 32 pages, 13 figures, 1 table, accepted in Phys. Rev. D

We study two possible cosmological consequences of a first-order phase transition in the temperature range of $1$ GeV to $10^3$ TeV: the generation of a stochastic gravitational wave background (SGWB) within the sensitivity of the Laser Interferometer Space Antenna (LISA) and, simultaneously, primordial magnetic fields that would evolve through the Universe's history and could be compatible with the lower bound from $\gamma$-ray telescopes on intergalactic magnetic fields (IGMF) at present time. We find that, if even a small fraction of the kinetic energy in sound waves is converted into MHD turbulence, a first order phase transition occurring at temperature between $1$ and $10^6$ GeV can give rise to an observable SGWB signal in LISA and, at the same time, an IGMF compatible with the lower bound from the $\gamma$-ray telescope MAGIC, for all proposed evolutionary paths of the magnetic fields throughout the radiation dominated era. For two values of the fraction of energy density converted into turbulence, $\varepsilon_{\rm turb}=0.1$ and $1$, we provide the range of first-order phase transition parameters (strength $\alpha$, duration $\beta^{-1}$, bubbles wall speed $v_w$, and temperature $T_*$), together with the corresponding range of magnetic field strength $B$ and correlation length $\lambda$, that would lead to the SGWB and IGMF observable with LISA and MAGIC. The resulting magnetic field strength at recombination can also correspond to the one that has been proposed to induce baryon clumping, previously suggested as a possible way to ease the Hubble tension. In the limiting case $\varepsilon_{\rm turb} \ll 1$, the SGWB is only sourced by sound waves, however, an IGMF is still generated. We find that values as small as $\varepsilon_{\rm turb} \sim O(10^{-13})$ (helical) and $O (10^{-9})$ (non-helical) can provide IGMF compatible with MAGIC's lower bound., Comment: 10 pages, 4 figures

The European Pulsar Timing Array (EPTA) and Indian Pulsar Timing Array (InPTA) collaborations have measured a low-frequency common signal in the combination of their second and first data releases respectively, with the correlation properties of a gravitational wave background (GWB). Such signal may have its origin in a number of physical processes including a cosmic population of inspiralling supermassive black hole binaries (SMBHBs); inflation, phase transitions, cosmic strings and tensor mode generation by non-linear evolution of scalar perturbations in the early Universe; oscillations of the Galactic potential in the presence of ultra-light dark matter (ULDM). At the current stage of emerging evidence, it is impossible to discriminate among the different origins. Therefore, in this paper, we consider each process separately, and investigate the implications of the signal under the hypothesis that it is generated by that specific process. We find that the signal is consistent with a cosmic population of inspiralling SMBHBs, and its relatively high amplitude can be used to place constraints on binary merger timescales and the SMBH-host galaxy scaling relations. If this origin is confirmed, this is the first direct evidence that SMBHBs merge in nature, adding an important observational piece to the puzzle of structure formation and galaxy evolution. As for early Universe processes, the measurement would place tight constraints on the cosmic string tension and on the level of turbulence developed by first-order phase transitions. Other processes would require non-standard scenarios, such as a blue-tilted inflationary spectrum or an excess in the primordial spectrum of scalar perturbations at large wavenumbers. Finally, a ULDM origin of the detected signal is disfavoured, which leads to direct constraints on the abundance of ULDM in our Galaxy., Comment: 30 pages, 23 figures, replaced to match the version published in Astronomy & Astrophysics, note the change in the numbering order in the series (now paper IV)

We study the propagation of cosmological gravitational wave (GW) backgrounds from the early radiation era until the present day in modified theories of gravity. Comparing to general relativity (GR), we study the effects that modified gravity parameters, such as the GW friction $\alpha_{\rm M}$ and the tensor speed excess $\alpha_{\rm T}$, have on the present-day GW spectrum. We use both the WKB estimate, which provides an analytical description but fails at superhorizon scales, and numerical simulations that allow us to go beyond the WKB approximation. We show that a constant $\alpha_{\rm T}$ makes relatively insignificant changes to the GR solution, especially taking into account the constraints on its value from GW observations by the LIGO--Virgo collaboration, while $\alpha_{\rm M}$ can introduce modifications to the spectral slopes of the GW energy spectrum in the low-frequency regime depending on the considered time evolution of $\alpha_{\rm M}$. The latter effect is additional to the damping or growth occurring equally at all scales that can be predicted by the WKB approximation. In light of the recent observations by pulsar timing array (PTA) collaborations, and the potential observations by future detectors such as SKA, LISA, DECIGO, BBO, or ET, we show that, in most of the cases, constraints cannot be placed on the effects of $\alpha_{\rm M}$ and the initial GW energy density $\mathcal{E}_{\rm GW}^*$ separately, but only on the combined effects of the two, unless the signal is observed at different frequency ranges. In particular, we provide some constraints on the combined effects from the reported PTA observations., Comment: 36 pages, 11 figures, 2 tables

We propose that the recent observations reported by the different Pulsar Timing Array (PTA) collaborations (i.e.~IPTA, EPTA, PPTA, and NANOGrav) of a common process over several pulsars could correspond to a stochastic gravitational wave background (SGWB) produced by turbulent sources in the early universe, in particular due to the magnetohydrodynamic (MHD) turbulence induced by primordial magnetic fields. I discuss recent results of numerical simulations of MHD turbulence and present an analytical template of the SGWB validated by the simulations. We use this template to constrain the magnetic field parameters using the results reported by the PTA collaborations. Finally, we compare the constraints on the primordial magnetic fields obtained from PTA with those from blazar signals observed by Fermi Large Area Telescope (LAT), from ultra high-energy cosmic rays, and from the cosmic microwave background. We show that a non-helical primordial magnetic field produced at the scale of the quantum chromodynamics phase transition is compatible with such constraints and it could additionally provide with a magnetic field at recombination that would help to alleviate the Hubble tension., Comment: Contribution to the 2022 Gravitation session of the 56th Rencontres de Moriond

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.

The NANOGrav, Parkes, European, and International Pulsar Timing Array (PTA) Collaborations have reported evidence for a common-spectrum process that can potentially correspond to a stochastic gravitational wave background (SGWB) in the 1--100 nHz frequency range. We consider the scenario in which this signal is produced by magnetohydrodynamic (MHD) turbulence in the early Universe, induced by a nonhelical primordial magnetic field at the energy scale corresponding to the quark confinement phase transition. We perform MHD simulations to study the dynamical evolution of the magnetic field and compute the resulting SGWB. We show that the SGWB output from the simulations can be very well approximated by assuming that the magnetic anisotropic stress is constant in time, over a time interval related to the eddy turnover time. The analytical spectrum that we derive under this assumption features a change of slope at a frequency corresponding to the GW source duration that we confirm with the numerical simulations. We compare the SGWB signal with the PTA data to constrain the temperature scale at which the SGWB is sourced, as well as the amplitude and characteristic scale of the initial magnetic field. We find that the generation temperature is constrained to be in the 1--200 MeV range, the magnetic field amplitude must be $>1$\% of the radiation energy density at that time, and the magnetic field characteristic scale is constrained to be $>10$\% of the horizon scale. We show that the turbulent decay of this magnetic field will lead to a field at recombination that can help to alleviate the Hubble tension and can be tested by measurements in the voids of the Large Scale Structure with gamma-ray telescopes like the Cherenkov Telescope Array., Comment: 28 pages, 9 figures, 1 table, accepted for publication in Phys. Rev. D, files and code available in https://github.com/AlbertoRoper/GW_turbulence

We study the leading-order nonlinear gravitational waves (GWs) produced by an electromagnetic (EM) stress in reheating magnetogenesis scenarios. Both nonhelical and helical magnetic fields are considered. By numerically solving the linear and leading-order nonlinear GW equations, we find that the GW energy from the latter is usually larger. We compare their differences in terms of the GW spectrum and parameterize the GW energy difference due to the nonlinear term, $\Delta\mathcal{E}_{\rm GW}$, in terms of EM energy $\mathcal{E}_{\rm EM}$ as $\Delta\mathcal{E}_{\rm GW}=(\tilde p\mathcal{E}_{\rm EM}/k_*)^3$, where $k_*$ is the characteristic wave number, $\tilde p=0.84$ and $0.88$ are found in the nonhelical and helical cases, respectively, with reheating around the QCD energy scale, while $\tilde p=0.45$ is found at the electroweak energy scale. We also compare the polarization spectrum of the linear and nonlinear cases and find that adding the nonlinear term usually yields a decrease in the polarization that is proportional to the EM energy density. We parameterize the fractional polarization suppression as $|\Delta \mathcal{P}_{\rm GW}/\mathcal{P}_{\rm GW}|=\tilde r \mathcal{E}_{\rm EM}/k_*$ and find $\tilde r = 1.2 \times 10^{-1}$, $7.2 \times 10^{-4}$, and $3.2 \times 10^{-2}$ for the helical cases with reheating temperatures $T_{\rm r} = 300 {\rm TeV}$, $8 {\rm GeV}$, and $120 {\rm MeV}$, respectively. Prospects of observation by pulsar timing arrays, space-based interferometers, and other novel detection proposals are also discussed., Comment: 17 pages, 7 figures, 3 tables, to be submitted to PRD

We use direct numerical simulations of decaying primordial hydromagnetic turbulence with helicity to compute the resulting gravitational wave (GW) production and its degree of circular polarization. We find a clear dependence of the polarization of the resulting GWs on the fractional helicity of the turbulent source and we show that driven magnetic fields produce GWs more efficiently than magnetic fields that are initially present, leading to larger spectral amplitudes. The helicity does not have a huge impact on the maximum spectral amplitude in any of the two types of turbulence considered. However, the GW spectrum at wave numbers away from the peak becomes smaller for larger values of the magnetic fractional helicity. The degree of circular polarization approaches zero at frequencies below the peak, and reaches its maximum at the peak. At higher frequencies, it stays finite if the magnetic field is initially present, and it approaches zero if it is driven. We predict that the spectral peak of the GW signal can be detected by LISA if the turbulent energy density is at least $\sim 3\%$ of the radiation energy density, and the characteristic scale is a hundredth of the horizon at the electroweak scale. We show that the resulting GW polarization is unlikely to be detectable by the anisotropies induced by our proper motion in the dipole response function of LISA. Such signals can, however, be detectable by cross-correlating data from the LISA-Taiji network for turbulent energy densities of $\sim 5\%$, and fractional helicity of 0.5 to 1. Second-generation space-based GW detectors, such as BBO and DECIGO, would allow for the detection of a larger range of the GW spectrum and smaller amplitudes of the magnetic field., Comment: 49 pages, 18 figures, 1 table, accepted in JCAP

I briefly discuss recent results of numerical simulations addressing the generation of a cosmological gravitational wave background produced by turbulence sources in the early universe. Contribution to the 2021 Gravitation session of the 55th Rencontres de Moriond., Comment: Contribution to the 2021 Gravitation session of the 55th Rencontres de Moriond

We study the gravitational wave (GW) signal sourced by primordial turbulence that is assumed to be present at cosmological phase transitions like the electroweak and quantum chromodynamics phase transitions. We consider various models of primordial turbulence, such as those with and without helicity, purely hydrodynamical turbulence induced by fluid motions, and magnetohydrodynamic turbulence whose energy can be dominated either by kinetic or magnetic energy, depending on the nature of the turbulence. We also study circularly polarized GWs generated by parity violating sources such as helical turbulence. Our ultimate goal is to determine the efficiency of GW production through different classes of turbulence. We find that the GW energy and strain tend to be large for acoustic or irrotational turbulence, even though its tensor mode amplitude is relatively small at most wave numbers. Only at very small wave numbers is the spectral tensor mode significant, which might explain the efficient GW production in that case., Comment: 25 pages, 8 figures

We perform direct numerical simulations to compute the net circular polarization of gravitational waves from helical (chiral) turbulent sources in the early Universe for a variety of initial conditions, including driven (stationary) and decaying turbulence. We investigate the resulting gravitational wave signal assuming different turbulent geneses such as magnetically or kinetically driven cases. Under realistic physical conditions in the early Universe we compute numerically the wave number-dependent polarization degree of the gravitational waves. We find that the spectral polarization degree strongly depends on the initial conditions. The peak of the spectral polarization degree occurs at twice the typical wavenumber of the source, as expected, and for fully helical decaying turbulence, it reaches its maximum of nearly 100\% {\it only} at the peak. We determine the temporal evolution of the turbulent sources as well as the resulting gravitational waves, showing that the dominant contribution to their spectral energy density happens shortly after the activation of the source. Only through an artificially prolonged decay of the turbulence can further increase of the gravitational wave amplitude be achieved. We estimate the detection prospects for the net polarization, arguing that its detection contains {\it clean} information (including the generation mechanisms, time, and strength) about the sources of possible parity violations in the early Universe., Comment: 9 pages, 5 figures; discussions added; references updated and added; conclusions unchanged; Phys. Rev. Research in press

The NANOGrav collaboration has recently reported evidence for the existence of a stochastic gravitational wave background in the 1-100 nHz frequency range. We argue that such background could have been produced by magneto-hydrodynamic (MHD) turbulence at the QCD scale. From the NANOGrav measurement one can infer the magnetic field parameters: comoving field strength close to microGauss and a correlation length close to 10\% of the Hubble radius at the QCD phase transition epoch. We point out that the turbulent decay of a non-helical magnetic field with such parameters leads to a magnetic field at the recombination epoch, which would be sufficiently strong to provide a solution to the Hubble tension problem, as recently proposed. We also show that the MHD turbulence interpretation of the NANOGrav signal can be tested via measurements of the relic magnetic field in the voids of the large scale structure, with gamma-ray telescopes like CTA., Comment: 6 pages, 2 figures

We perform direct numerical simulations of magnetohydrodynamic turbulence in the early universe and numerically compute the resulting stochastic background of gravitational waves and relic magnetic fields. These simulations do not make the simplifying assumptions of earlier analytic work. If the turbulence is assumed to have an energy-carrying scale that is about a hundredth of the Hubble radius at the time of generation, as expected in a first-order phase transition, the peak of gravitational wave power will be in the mHz frequency range for a signal produced at the electroweak scale. The efficiency of gravitational wave (GW) production varies significantly with how the turbulence is driven. Detectability of turbulence at the electroweak scale by the planned Laser Interferometer Space Antenna (LISA) requires anywhere from 0.1% to 10% of the thermal plasma energy density to be in plasma motions or magnetic fields, depending on the model of the driving process. Our results predict a new universal form below the spectral peak frequency that is shallower than previously thought. This implies larger values of the GW energy spectra in the low-frequency range. This extends the range where turbulence is detectable with LISA to lower frequencies, corresponding to higher energy scales than the assumed energy-carrying scale., Comment: 11 pages, 1 appendix, 7 + 1 (appendix) figures, 2 tables, published in PRD, updated figures and corrected normalization error

Blazar observations point toward the possible presence of magnetic fields over intergalactic scales of the order of up to $\sim1\,$Mpc, with strengths of at least $\sim10^{-16}\,$G. Understanding the origin of these large-scale magnetic fields is a challenge for modern astrophysics. Here we discuss the cosmological scenario, focussing on the following questions: (i) How and when was this magnetic field generated? (ii) How does it evolve during the expansion of the universe? (iii) Are the amplitude and statistical properties of this field such that they can explain the strengths and correlation lengths of observed magnetic fields? We also discuss the possibility of observing primordial turbulence through direct detection of stochastic gravitational waves in the mHz range accessible to LISA., Comment: Invited talk given at FM8 "New Insights in Extragalactic Magnetic Fields", XXX IAU GA, Vienna, Austria, August 29, 2018; 4 pages; 2 figures; accepted for publication

Gradient- and curl-type or E- and B-type polarizations have been routinely analyzed to study the physics contributing to the cosmic microwave background polarization and galactic foregrounds. They characterize the parity-even and parity-odd properties of the underlying physical mechanisms, for example hydromagnetic turbulence in the case of dust polarization. Here we study spectral correlation functions characterizing the parity-even and parity-odd parts of linear polarization for homogeneous and inhomogeneous turbulence to show that only the inhomogeneous helical case can give rise to a parity-odd polarization signal. We also study nonhelical turbulence and suggest that a strong nonvanishing (here negative) skewness of the E polarization is responsible for an enhanced ratio of the EE to the BB (quadratic) correlation in both helical and nonhelical cases. This could explain the enhanced EE/BB ratio observed recently for dust polarization. We close with a preliminary assessment of using linear polarization of the Sun to characterize its helical turbulence without being subjected to the pi ambiguity that magnetic inversion techniques have to address., Comment: 11 pages, 13 figures, 4 tables, ApJ, in press

Hydromagnetic turbulence produced during phase transitions in the early universe can be a powerful source of stochastic gravitational waves (GWs). GWs can be modelled by the linearised spatial part of the Einstein equations sourced by the Reynolds and Maxwell stresses. We have implemented two different GW solvers into the {\sc Pencil Code} -- a code which uses a third order timestep and sixth order finite differences. Using direct numerical integration of the GW equations, we study the appearance of a numerical degradation of the GW amplitude at the highest wavenumbers, which depends on the length of the timestep -- even when the Courant--Friedrichs--Lewy condition is ten times below the stability limit. This degradation leads to a numerical error, which is found to scale with the third power of the timestep. A similar degradation is not seen in the magnetic and velocity fields. To mitigate numerical degradation effects, we alternatively use the exact solution of the GW equations under the assumption that the source is constant between subsequent timesteps. This allows us to use a much longer timestep, which cuts the computational cost by a factor of about ten., Comment: 30 pages, 9 figures, 2 appendices. Published at GAFD, special edition on "Physics and Algorithms of the Pencil Code", minor typos corrected

We present new simulations of decaying hydromagnetic turbulence for a relativistic equation of state relevant to the early universe. We compare helical and nonhelical cases either with kinetically or magnetically dominated initial fields. Both kinetic and magnetic initial helicities lead to maximally helical magnetic fields after some time, but with different temporal decay laws. Both are relevant to the early universe, although no mechanisms have yet been identified that produce magnetic helicity with strengths comparable to the big bang nucleosynthesis limit at scales comparable to the Hubble horizon at the electroweak phase transition. Nonhelical magnetically dominated fields could still produce picoGauss magnetic fields under most optimistic conditions. Only helical magnetic fields can potentially have nanoGauss strengths at scales up to 30 kpc today., Comment: 17 pages, 13 figures; submitted to PRD

We show that in decaying hydromagnetic turbulence with initial kinetic helicity, a weak magnetic field eventually becomes fully helical. The sign of magnetic helicity is opposite to that of the kinetic helicity - regardless of whether or not the initial magnetic field was helical. The magnetic field undergoes inverse cascading with the magnetic energy decaying approximately like t^{-1/2}. This is even slower than in the fully helical case, where it decays like t^{-2/3}. In this parameter range, the product of magnetic energy and correlation length raised to a certain power slightly larger than unity, is approximately constant. This scaling of magnetic energy persists over long time scales. At very late times and for domain sizes large enough to accommodate the growing spatial scales, we expect a cross-over to the t^{-2/3} decay law that is commonly observed for fully helical magnetic fields. Regardless of the presence or absence of initial kinetic helicity, the magnetic field experiences exponential growth during the first few turnover times, which is suggestive of small-scale dynamo action. Our results have applications to a wide range of experimental dynamos and astrophysical time-dependent plasmas, including primordial turbulence in the early universe., Comment: 15 pages, 14 figures, submitted to Phys. Rev. Fluids