Stochastic Gravitational Wave Research Articles

Development of generic no-scale inflation

We develop generalized no-scale supergravity models of inflation, and then study the corresponding cosmological predictions as well as the formation of primordial black holes (PBHs) and scalar-induced gravitational waves (SIGWs). With a new parameter 0 < a ≤ 1, the generalized no-scale supergravity provides the continuous connections among the generic no-scale supergravity from string theory compactifications. The resulting prediction of the CMB, spectrum index ns , and tensor-to-scalar ratio r can be highly consistent with the latest Planck/BICEP/Keck Array observations. Notably, the models with a ≠ 1 give a smaller ratio r ≤ 10-3, which is flexible even under the anticipated tighter observational constraints at the future experiments. Additionally, these models have the potential to generate a broad-band stochastic gravitational wave background, and thus explain the NANOGrav 15yr signal. Furthermore, they predict the formation of PBHs with various mass scales, which could account for a significant portion of dark matter relic density in the Universe.

Stochastic gravitational wave background from binary black hole mergers dynamically assembled in dense star clusters

Stochastic gravitational wave background from binary black hole mergers dynamically assembled in dense star clusters

Impact of a local CP -odd domain in hot QCD on the axionic domain-wall interpretation of NANOGrav 15-year data

We argue that the axionic domain-wall with a QCD bias may be incompatible with the NANOGrav 15-year data on a stochastic gravitational wave (GW) background, when the domain wall network collapses in the hot-quantum chromodynamics (QCD) induced local CP-odd domain. This is due to the drastic suppression of the QCD bias set by the QCD topological susceptibility in the presence of the CP-odd domain with nonzero θ parameter of order one which the QCD sphaleron could generate. We quantify the effect on the GW signals by working on a low-energy effective model of Nambu-Jona-Lasinio type in the mean field approximation. We find that only at θ=π, the QCD bias tends to get significantly large enough due to the criticality of the thermal CP restoration, which would, however, give too big signal strengths to be consistent with the NANOGrav 15-year data and would also be subject to the strength of the phase transition at the criticality. Published by the American Physical Society 2024

Stochastic gravitational wave background: birth from string-wall death

We study a new source of stochastic gravitational wave background (SGWB) from the final collapse of a network of topological defects. Typically, the final collapse is considered negligible for generating gravitational waves (GWs) due to its subdominance compared with the network's long-term evolution in the scaling regime. However, in some cases, a network can be driven outside of horizon by inflation and later re-enter horizon. Then, the network's final collapse after re-entering horizon becomes the dominant GW source and therefore cannot be neglected. We demonstrate this phenomenon in the context of N DW = 1 string-wall networks which naturally arise in axion models, although the framework can be generalized to other types of topological networks. The final collapse of walls bounded by strings releases GWs. Our calculation of the corresponding GW spectrum suggests it could be related to the first few bins of the nano-Hertz SGWB signal possibly detected by various Pulsar Timing Array (PTA) collaborations. However, it is important to note that such GW spectrum falls within a relatively narrow frequency range, which may not completely account for the PTA signal that spans more than one order of magnitude in frequency. Furthermore, with different parameter choices, the resultant GWs generated in this mechanism could be probed by various GW interferometry experiments.

Probing the galactic and extragalactic gravitational wave backgrounds with space-based interferometers

We employ the formalism developed in [1] and [2] to study the prospect of detecting an anisotropic Stochastic Gravitational Wave Background (SGWB) with the Laser Interferometer Space Antenna (LISA) alone, and combined with the proposed space-based interferometer Taiji. Previous analyses have been performed in the frequency domain only. Here, we study the detectability of the individual coefficients of the expansion of the SGWB in spherical harmonics, by taking into account the specific motion of the satellites. This requires the use of time-dependent response functions, which we include in our analysis to obtain an optimal estimate of the anisotropic signal. We focus on two applications. Firstly, the reconstruction of the anisotropic galactic signal without assuming any prior knowledge of its spatial distribution. We find that both LISA and LISA with Taiji cannot put tight constraints on the harmonic coefficients for realistic models of the galactic SGWB. We then focus on the discrimination between a galactic signal of known morphology but unknown overall amplitude and an isotropic extragalactic SGWB component of astrophysical origin. In this case, we find that the two surveys can confirm, at a confidence level ≳ 3σ, the existence of both the galactic and extragalactic background if both have amplitudes as predicted in standard models. We also find that, in the LISA-only case, the analysis in the frequency domain (under the assumption of a time average of data taken homogeneously across the year) provides a nearly identical determination of the two amplitudes as compared to the optimal analysis.

Gravitational wave signals from early matter domination: interpolating between fast and slow transitions

An epoch of matter domination in the early universe can enhance the primordial stochastic gravitational wave signal, potentially making it detectable to upcoming gravitational wave experiments. However, the resulting gravitational wave signal is quite sensitive to the end of the early matter-dominated epoch. If matter domination ends gradually, a cancellation results in an extremely suppressed signal, while in the limit of an instantaneous transition, there is a resonant-like enhancement. The end of the matter dominated epoch cannot be instantaneous, however, and previous analyses have used a Gaussian smoothing technique to account for this, and consider only a limited regime around the fast transition limit. In this work, we present a study of the enhanced gravitational wave signal from early matter domination without making either approximation and show how the signal smoothly evolves from the strongly suppressed to strongly enhanced regimes.

Vector dark matter, inflation, and non-minimal couplings with gravity

We propose a cosmological dark matter production mechanism in the form of a longitudinal massive vector boson. We build upon the work [1] including non-minimal couplings of the massive vector with gravity, developing a well motivated set-up from an effective field theory perspective. We carefully track the dynamics of vector field in passing from inflation to radiation dominated universe to show that the late time abundance of longitudinal modes — excited initially by the quantum fluctuations during inflation — can provide the observed dark matter abundance for sufficiently weak non-minimal coupling and wide range of vector masses 5 × 10-7 ≲ m [eV] ≲ 5 × 103. The final abundance of dark matter depends on two parameter, the vector mass and its non-minimal coupling with gravity. We discuss experimental venues to probe this framework, including the production of a stochastic gravitational wave background. The latter is especially interesting, as the same mechanism that generates dark matter can potentially lead to the production of gravitational waves in the LISA frequency band, through the second-order effects of large dark matter iso-curvature perturbations at small scales. We take a first step in this direction, identifying the potential information that gravitational wave experiments can provide on the parameter space of dark matter within this scenario.

Connecting low-redshift LISA massive black hole mergers to the nHz stochastic gravitational wave background

Pulsar Timing Array (PTA) experiments worldwide recently reported evidence of a nHz stochastic gravitational wave background (sGWB) compatible with the existence of slowly inspiralling massive black hole (MBH) binaries (MBHBs). The shape of the signal contains valuable information about the evolution of z &lt; 1 MBHs above 108 M⊙, suggesting a faster dynamical evolution of MBHBs towards the gravitational-wave-driven inspiral or a larger MBH growth than usually assumed. In this work, we investigate if the nHz sGWB could also provide constraints on the population of merging lower-mass MBHBs (&lt; 107 M⊙) detectable by LISA. To this end, we use the L-Galaxies semi-analytical model applied to the Millennium suite of simulations. We generate a population of MBHs compatible simultaneously with current electromagnetic and nHz sGWB constraints by including the possibility that, in favourable environments, MBHs can accrete gas beyond the Eddington limit. The predictions of this new model for the sGWB show that the global (integrated up to high-z) LISA detection rate is not significantly affected when compared to a fiducial model whose nHz sGWB signal is ∼2 times smaller. In both cases, the global rate yields ∼12 yr−1 and is dominated by systems of 105 − 6 M⊙. The main differences are limited to low-z (z &lt; 3), high-mass (&gt; 106 M⊙) LISA MBHBs. The model compatible with the latest PTA results predicts up to ∼1.6 times more detections, with a rate of ∼1 yr−1. We find that these LISA MBHB systems have 50% probability of shining with bolometric luminosities &gt; 1043 erg s−1. Hence, in case PTA results are confirmed and given the current MBH modelling, our findings suggest there will be higher chances to perform multimessenger studies with LISA MBHB than previously expected.

Constraints on inflation with null energy condition violation from advanced LIGO and advanced Virgo's first three observing runs

The null energy condition (NEC) is a cornerstone of general relativity, and its violation could leave observable imprints in the cosmic gravitational wave spectrum. Theoretical models suggest that NEC violations during inflation can amplify the primordial tensor power spectrum, leading to distinct features in the stochastic gravitational wave background (SGWB). In this work, we search for these NEC-violating signatures in the SGWB using data from Advanced LIGO and Advanced Virgo's first three observing runs. Our analysis reveals no statistically significant evidence of such signals, allowing us to place stringent upper limits on the tensor power spectrum amplitude, P T,2, during the second inflationary stage. Specifically, we find that P T,2 ≲ 0.15 at a 95% confidence level. Notably, this upper limit is consistent with constraints derived from pulsar timing array observations, reinforcing the hypothesis that NEC violations during inflation could explain the signal detected by pulsar timing arrays. Our findings contribute to a deeper understanding of the early Universe and highlight the potential of current and future gravitational wave experiments in probing the physics of inflation and NEC violations.

Can Supercooled Phase Transitions Explain the Gravitational Wave Background Observed by Pulsar Timing Arrays?

Several pulsar timing array collaborations recently reported evidence of a stochastic gravitational wave background (SGWB) at nHz frequencies. While the SGWB could originate from the merger of supermassive black holes, it could be a signature of new physics near the 100MeV scale. Supercooled first-order phase transitions (FOPTs) that end at the 100MeV scale are intriguing explanations, because they could connect the nHz signal to new physics at the electroweak scale or beyond. Here, however, we provide a clear demonstration that it is not simple to create a nHz signal from a supercooled phase transition, due to two crucial issues that could rule out many proposed supercooled explanations and should be checked. As an example, we use a model based on nonlinearly realized electroweak symmetry that has been cited as evidence for a supercooled explanation. First, we show that a FOPT cannot complete for the required transition temperature of around 100MeV. Such supercooling implies a period of vacuum domination that hinders bubble percolation and transition completion. Second, we show that even if completion is not required or if this constraint is evaded, the Universe typically reheats to the scale of any physics driving the FOPT. The hierarchy between the transition and reheating temperature makes it challenging to compute the spectrum of the SGWB.

LIGO-Virgo constraints on dark matter and leptogenesis triggered by a first order phase transition at high scale

We study the possibility of constraining a scenario with high scale first order phase transition (FOPT) responsible for the cogenesis of baryon and dark matter using gravitational wave (GW) (non)observations. While the FOPT at high scale is responsible for generating baryon asymmetry through leptogenesis and dark matter via the mechanism, the resulting GW spectrum falls within the ongoing LIGO-Virgo experimental sensitivity. The dark matter is preferred to be in the nonthermal ballpark with sub-GeV masses and the criteria of successful dark matter relic rules out a large portion of the parameter space consistent with high scale FOPT and successful leptogenesis. Some part of the parameter space allowed from dark matter and leptogenesis criteria also gives rise to a large signal-to-noise ratio at ongoing experiments and hence can be disfavored in a conservative way from the nonobservation of such stochastic GW background. Future data from ongoing and planned experiments will offer a complementary and indirect probe of the remaining parameter space which is typically outside the reach of any direct experimental probe. Published by the American Physical Society 2024

A stochastic gravitational wave background in LISA from unresolved white dwarf binaries in the Large Magellanic Cloud

ABSTRACT The Laser Interferometer Space Antenna (LISA) is expected to detect a wide variety of gravitational wave sources in the mHz band. Some of these signals will elude individual detection, instead contributing as confusion noise to one of several stochastic gravitational-wave backgrounds (SGWBs) – notably including the ‘Galactic foreground’, a loud signal resulting from the superposition of millions of unresolved double white dwarf binaries (DWDs) in the Milky Way. It is possible that similar, weaker SGWBs will be detectable from other DWD populations in the local Universe, including the Large Magellanic Cloud (LMC). We use the Bayesian LISA Inference Package (blip) to investigate the possibility of an anisotropic SGWB generated by unresolved DWDs in the LMC. To do so, we compute the LMC SGWB from a realistic DWD population generated via binary population synthesis, simulate 4 years of time-domain data with blip comprised of stochastic contributions from the LMC SGWB and the LISA detector noise, and analyse this data with blip’s spherical harmonic anisotropic SGWB search. We also consider the case of spectral separation from the Galactic foreground. We present the results of these analyses and show, for the first time, that the unresolved DWDs in the LMC will comprise a significant SGWB for LISA.

Phases of Pseudo-Nambu-Goldstone bosons

We study the vacuum dynamics of pseudo-Nambu-Goldstone bosons (pNGBs) for SO(N + 1) → SO(N) spontaneous and explicit symmetry breaking. We determine the magnitude of explicit symmetry breaking consistent with an EFT description of the effective potential at zero and finite temperatures. We expose and clarify novel additional vacuum transitions that can arise for generic pNGBs below the initial scale of SO(N + 1) → SO(N) spontaneous symmetry breaking, which may have phenomenological relevance. In this respect, two phenomenological scenarios are analyzed: thermal and supercooled dark sector pNGBs. In the thermal scenario the vacuum transition is first-order but very weak. For a supercooled dark sector we find that, depending on the sign of the explicit symmetry breaking, one can have a symmetry-restoring vacuum transition SO(N – 1) → SO(N) which can be strongly first-order, with a detectable stochastic gravitational wave background signal.

Quantum tunneling in the early universe: stable magnetic monopoles from metastable cosmic strings

We present a novel mechanism for producing topologically stable monopoles (TSMs) from the quantum mechanical decay of metastable cosmic strings in the early universe. In an SO(10) model this mechanism yields TSMs that carry two units (4π/e) of Dirac magnetic charge as well as some color magnetic charge which is screened. For a dimensionless string tension parameter Gμ ≈ 10-9-10-5, the monopoles are superheavy with masses of order 1015-1017 GeV. Monopoles with masses of order 108-1014 GeV arise from metastable strings for Gμ values from ∼ 10-22 to 10-10. We identify the parameter space for producing these monopoles at an observable level with detectors such as IceCube and KM3NeT. For lower Gμ values the ultra-relativistic monopoles should be detectable at Pierre Auger and ANITA. The stochastic gravitational wave emission arise from metastable strings with Gμ ∼ 10-9-10-5 and should be accessible at HLVK and future detectors including the Einstein Telescope and Cosmic Explorer. An E 6 extension based on this framework would yield TSMs from the quantum mechanical decay of metastable strings that carry three units (6π/e) of Dirac magnetic charge.

Dissecting the stochastic gravitational wave background with astrometry

Astrometry, the precise measurement of star motions, offers an alternative avenue to investigate low-frequency gravitational waves through the spatial deflection of photons, complementing pulsar timing arrays reliant on timing residuals. Upcoming data from Gaia, Theia, and Roman can not only cross-check pulsar timing array findings but also explore the uncharted frequency range bridging pulsar timing arrays and LISA. We present an analytical framework to evaluate the feasibility of detecting a gravitational wave background, considering measurement noise and the intrinsic variability of the stochastic background. Furthermore, we highlight astrometry's crucial role in uncovering key properties of the gravitational wave background, such as spectral index and chirality, employing information-matrix analysis. Finally, we simulate the emergence of quadrupolar correlations, commonly referred to as the generalized Hellings-Downs curves.

Observing gravitational waves with solar system astrometry

The subtle influence of gravitational waves on the apparent positioning of celestial bodies offers novel observational windows [1,2,3,4]. We calculate the expected astrometric signal induced by an isotropic Stochastic Gravitational Wave Background (SGWB) in the short distance limit. Our focus is on the resultant proper motion of Solar System objects, a signal on the same time scales addressed by Pulsar Timing Arrays (PTA). We derive the corresponding astrometric deflection patterns, finding that they manifest as distinctive dipole and quadrupole correlations or, in some cases, may not be present. Our analysis encompasses both Einsteinian and non-Einsteinian polarisations. We estimate the upper limits for the amplitude of SGWBs that could be obtained by tracking the proper motions of large numbers of solar system objects such as asteroids. We find that for SGWBs with negative spectral indices, such as that generated by Super Massive Black Hole Binaries (SMBHB), the constraints from these observations could rival those from PTAs. With the Gaia satellite and the Vera C. Rubin Observatory poised to track an extensive sample of asteroids — ranging from 𝒪(105) to 𝒪(106), we highlight the significant future potential for similar surveys to contribute to our understanding of the SGWB.

Percolating cosmic string loops from evaporating primordial black holes

The Pulsar timing data from NANOGrav Collaboration has regenerated interest in the possibility of observing stochastic gravitational wave background arising from cosmic strings. Standard theory of formation of cosmic strings is based on a spontaneous symmetry breaking (SSB) phase transition in the early universe, with a string network forming via the so called Kibble mechanism. This string network rapidly evolves and reaches a scaling solution with a given spectrum of string loops and long strings at any given time. This scenario necessarily requires that the entire observable Universe goes through the SSB phase transition. This would not be possible, e.g., in models of low energy inflation, where the reheat temperature is much lower than the energy scale of cosmic strings. We point out a very different possibility, where a network of even high energy scale cosmic strings can form when the temperature of the Universe is much lower. We consider local heating of plasma in the early universe by evaporating primordial black holes (PBHs). It is known that for suitable masses of PBHs, Hawking radiation of evaporating primordial black holes may re-heat the surrounding plasma to high temperatures, restoring certain symmetries locally which are broken at the ambient temperature of the Universe at that stage. Expansion of the hot plasma cools it so that the locally restored symmetry is spontaneously broken again. If this SSB supports formation of cosmic strings, then string loops will form in this region around the PBH. Further, resulting temperature gradients lead to large pressure gradients such that plasma will develop radial flow with the string loops getting stretched as they get dragged by the flow. For a finite density of PBHs of suitable masses, one will get local hot spots, each one contributing to expanding cosmic string loops. For suitable PBH density, the loops from different regions may intersect. If that happens, then intercommutation of strings can lead to percolation, leading to the possibility of formation of infinite string network, even when the entire universe never goes through the respective SSB phase transition.

Fully non-Gaussian Scalar-Induced Gravitational Waves

Scalar-Induced Gravitational Waves (SIGWs) represent a particular class of primordial signals which are sourced at second-order in perturbation theory whenever a scalar fluctuation of the metric is present. They form a guaranteed Stochastic Gravitational Wave Background (SGWB) that, depending on the amplification of primordial scalar fluctuations, can be detected by GW detectors. The amplitude and the frequency shape of the scalar-induced SGWB can be influenced by the statistical properties of the scalar density perturbations. In this work we study the intuitive physics behind SIGWs and we analyze the imprints of local non-Gaussianity of the primordial curvature perturbation on the GW spectrum. We consider all the relevant non-Gaussian contributions up to fifth-order in the scalar seeds without any hierarchy, and we derive the related GW energy density ΩGW(f). We perform a Fisher matrix analysis to understand to which accuracy non-Gaussianity can be constrained with the LISA detector, which will be sensitive in the milli-Hertz frequency band. We find that LISA, neglecting the impact of astrophysical foregrounds, will be able to measure the amplitude, the width and the peak of the spectrum with an accuracy up to 𝒪(10-4), while non-Gaussianity can be measured up to 𝒪(10-3). Finally, we discuss the implications of our non-Gaussianity expansion on the fraction of Primordial Black Holes.

Axionic defects in the CMB: birefringence and gravitational waves

The evidence for a non-vanishing isotropic cosmic birefringence in recent analyses of the CMB data provides a tantalizing hint for new physics. Domain wall (DW) networks have recently been shown to generate an isotropic birefringence signal in the ballpark of the measured value when coupled to photons. In this work, we explore the axionic defects hypothesis in more detail and extending previous results to annihilating and late-forming networks, and by pointing out other smoking-gun signatures of the network in the CMB spectrum such as the anisotropic birefringent spectrum and B-modes. We also argue that the presence of cosmic strings in the network does not hinder a large isotropic birefringence signal because of an intrinsic environmental contribution coming from low redshifts thus leaving open the possibility that axionic defects can explain the signal. Regarding the remaining CMB signatures, with the help of dedicated 3D numerical simulations of DW networks, that we took as a proxy for the axionic defects, we show how the anisotropic birefringence spectrum combined with a tomographic approach can be used to infer the formation and annihilation time of the network. From the numerical simulations, we also computed the spectrum of gravitational waves (GWs) generated by the network in the post-recombination epoch and use previous searches for stochastic GW backgrounds in the CMB to derive for the first time a bound on the tension and abundance of networks with DWs that annihilate after recombination. Our bounds extend to the case where the network survives until the present time and improve over previous bounds by roughly one order of magnitude. Finally, we show the interesting prospects for detecting B-modes of DW origin with future CMB experiments.

Pulsar timing array source ensembles

The stochastic gravitational wave background for pulsar timing arrays is often modeled by a Gaussian ensemble which is isotropic and unpolarized. However, the Universe has a discrete set of polarized gravitational wave sources at specific sky locations. Can we trust that the Gaussian ensemble is an accurate description? To investigate this, we explicitly construct an ensemble containing N individual binary sources with circular orbits. The orbital inclination angles are randomly distributed, hence the individual sources are elliptically polarized. We then compute the first two moments of the Hellings and Downs correlation, as well as the pulsar-averaged correlation mean and (cosmic) variance. The first moments are the same as for a previously studied ensemble of circularly polarized sources. However, the second moments, and hence the variances, are different for the two ensembles. While neither discrete source model is exactly described by a Gaussian ensemble, we show that in the limit of large N, the differences are small. Published by the American Physical Society 2024