Stochastic Background Research Articles

Gravitational wave spectrum of chain inflation

Chain inflation is an alternative to slow-roll inflation in which the inflaton tunnels along a large number of consecutive minima in its potential. In this work we perform the first comprehensive calculation of the gravitational wave (GW) spectrum of chain inflation. In contrast to slow-roll inflation the latter does not stem from quantum fluctuations of the gravitational field during inflation, but rather from the bubble collisions during the first-order phase transitions associated with vacuum tunneling. Our calculation is performed within an effective theory of chain inflation which builds on an expansion of the tunneling rate capturing most of the available model space. The effective theory can be seen as chain inflation’s analog of the slow-roll expansion in rolling models of inflation. The near scale-invariance of the scalar power spectrum translates to a quasiperiodic shape of the inflaton potential in chain inflation, with the tunneling rate changing very slowly during the e-folds leading to cosmic microwave background observables. We show that chain inflation produces a very characteristic double-peak GW spectrum: a faint high-frequency peak associated with the gravitational radiation emitted during inflation, and a strong low-frequency peak associated with the graceful exit from chain inflation (marking the transition to the radiation-dominated epoch). There exist very exciting prospects to test the gravitational wave signal from chain inflation at the aLIGO-aVIRGO-KAGRA network, at LISA and /or at pulsar timing array experiments. A particularly intriguing possibility we point out is that chain inflation could be the source of the stochastic gravitational wave background recently detected by NANOGrav, PPTA, EPTA, and CPTA. We also show that the gravitational wave signal of chain inflation is often accompanied by running/ higher running of the scalar spectral index to be tested at future cosmic microwave background experiments. Published by the American Physical Society 2024

Probing stochastic ultralight dark matter with space-based gravitational-wave interferometers

Ultralight particles are theoretically well-motivated dark matter candidates. In the vicinity of the Solar System, these ultralight particles can be described as a superposition of plane waves, resulting in a stochastic field with sizable amplitude fluctuations on scales determined by the velocity dispersion of dark matter. In this work, we systematically investigate the sensitivity of space-based gravitational-wave interferometers to the stochastic ultralight dark matter field within the frequentist framework. We derive the projected sensitivity of a single detector using the time-delay interferometry. Our results show that space-based gravitational-wave interferometers have the potential to probe unconstrained regions in parameter space and improve the current limit on coupling strengths. Furthermore, we explore the sensitivity of a detector network and investigate the optimal configuration for ultralight dark matter detection. We introduce the overlap reduction function for ultralight dark matter, which quantifies the degree of correlation between the signals observed by different detectors. We find that the configuration, where the signals observed by two detectors are uncorrelated, is the optimal choice for ultralight dark matter detection due to a smaller chance of missing signals. This contrasts with the detection of stochastic gravitational-wave background, where the correlated configuration is preferred. Our results may provide useful insights for potential joint observations involving space-based gravitational-wave detectors like LISA and Taiji, as well as other ultralight dark matter detection networks operating in the coherence limit. Published by the American Physical Society 2024

Measuring the circular polarization of gravitational waves with pulsar timing arrays

The circular polarization of the stochastic gravitational wave background (SGWB) is a key observable for characterizing the origin of the signal detected by Pulsar Timing Array (PTA) collaborations. Both the astrophysical and the cosmological SGWB can have a sizeable amount of circular polarization, due to Poisson fluctuations in the source properties for the former, and to parity violating processes in the early universe for the latter. Its measurement is challenging since PTA are blind to the circular polarization monopole, forcing us to turn to anisotropies for detection. We study the sensitivity of current and future PTA datasets to circular polarization anisotropies, focusing on realistic modelling of intrinsic and kinematic anisotropies for astrophysical and cosmological scenarios respectively. Our results indicate that the expected level of circular polarization for the astrophysical SGWB should be within the reach of near future datasets, while for cosmological SGWB circular polarization is a viable target for more advanced SKA-type experiments. Published by the American Physical Society 2024

Resonant conversion of gravitational waves in neutron star magnetospheres

High-frequency gravitational waves are the subject of rapidly growing interest in the theoretical and experimental community. In this work we calculate the resonant conversion of gravitational waves into photons in the magnetospheres of neutron stars via the inverse Gertsenshtein mechanism. The resonance occurs in regions where the vacuum birefringence effects cancel the classical plasma contribution to the photon dispersion relation, leading to a massless photon in the medium which becomes kinematically matched to the graviton. We set limits on the amplitude of a possible stochastic background of gravitational waves using X-ray and IR flux measurements of neutron stars. Using Chandra (2–8 keV) and NuSTAR (3–79 keV) observations of RX J1856.6-3754, we set strain limits hclim≃10−26–10−24 in the frequency range 5×1017 Hz≲f≲2×1019 Hz. Our limits are many orders of magnitude stronger than existing constraints from individual neutron stars at the same frequencies. We also use recent JWST observations of the Magnetar 4U 0142+61 in the range 2.7×1013 Hz≲f≲5.9×1013 Hz, setting a limit hclim≃5×10−19. These constraints are in complementary frequency ranges to laboratory searches with CAST, OSQAR and ALPS II. We expect these limits to be improved both in reach and breadth with a more exhaustive use of telescope data across the full spectrum of frequencies and targets. Published by the American Physical Society 2024

Beyond the Hellings–Downs curve: Non-Einsteinian gravitational waves in pulsar timing array correlations

Pulsar timing arrays (PTAs) have revealed galaxy-size gravitational waves (GWs) in the form of a stochastic gravitational wave background (SGWB), correlating the radio pulses emitted by millisecond pulsars. This discovery naturally leads to the question of the origin and the nature of the SGWB; the latter is synonymous to testing how quadrupolar the inter-pulsar spatial correlation is. In this paper, we investigate the nature of the SGWB by considering correlations beyond the Hellings–Downs (HD) curve of Einstein’s general relativity. We scrutinize the HD and non-Einsteinian GW correlations with the North American Nanohertz Observatory for Gravitational Waves and the Chinese PTA data, and find that both data sets allow a graviton mass of mg ≲ 1.04 × 10−22 eV/c2 and subluminal traveling waves. We discuss gravitational physics scenarios beyond general relativity that could host non-Einsteinian GW correlations in the SGWB and highlight the importance of the cosmic variance inherited from stochastic variations across realizations in interpreting PTA observations.

Constraining gravitational-wave backgrounds from conversions into photons in the Galactic magnetic field

High-frequency gravitational waves (f≳1 MHz) may provide a unique signature for the existence of exotic physics. The lack of current and future gravitational-wave experiments sensitive at those frequencies leads to the need of employing different indirect techniques. Notably, one of the most promising ones is constituted by graviton-photon conversions in magnetic fields. In this work, we focus on conversions of a gravitational-wave background into photons inside the Milky Way magnetic field, taking into account the state-of-the-art models for both regular and turbulent components. We discuss how graviton-to-photon conversions may lead to imprints in the cosmic photon background spectrum in the range of frequencies f∼109–1026 Hz, where the observed photon flux is widely explained by astrophysics emission models. Hence, the absence of any significant evidence for a diffuse photon flux induced by graviton-photon conversions allows us to set stringent constraints on the gravitational-wave strain hc, strengthening current astrophysical bounds by ∼1–2 orders of magnitude in the whole range of frequencies considered. Published by the American Physical Society 2024

Primordial gravitational waves of big bounce cosmology in light of stochastic gravitational wave background

Primordial gravitational waves of big bounce cosmology in light of stochastic gravitational wave background

Detection prospects of gravitational waves from SU(2) axion inflation

We study detection prospects of a gravitational-wave background (GWB) sourced by SU(2) gauge fields considering all possible observational constraints. More precisely, we consider bounds set by cosmic microwave background measurements, primordial black hole overproduction, as well as backreaction of the gauge fields on the background evolution. Gravitational-waves data from the first three observing runs of the LIGO-Virgo-KAGRA Collaboration show no evidence for a GWB contribution from axion inflation. However, we are able to place conservative constraints on the parameters of the SU(2) inflation with current data. We investigate conditions on the inflationary potential that would lead to a detectable signal that evades astrophysical and cosmological constraints and discuss detection prospects for third generation networks. Published by the American Physical Society 2024

Measuring the anisotropies in astrophysical and cosmological gravitational-wave backgrounds with Taiji and LISA networks

Measuring the anisotropies in astrophysical and cosmological gravitational-wave backgrounds with Taiji and LISA networks

Memory burden effect mimics reheating signatures on SGWB from ultra-low mass PBH domination

Ultra-low mass primordial black holes (PBH), briefly dominating the expansion of the universe, would leave detectable imprints in the secondary stochastic gravitational wave background (SGWB). Such a scenario leads to a characteristic doubly peaked spectrum of SGWB and strongly depends on the Hawking evaporation of such light PBHs. However, these observable signatures are significantly altered if the memory burden effect during the evaporation of PBHs is taken into account. We show that for the SGWB induced by PBH density fluctuations, the memory burden effects on the Hawking evaporation of ultra-low mass PBHs can mimic the signal arising due to the non-standard reheating epoch before PBH domination. Finally, we point out that this degeneracy can be broken by the simultaneous detection of the first peak in the SGWB, which is typically induced by the inflationary adiabatic perturbations.

Impact of a primordial gravitational wave background on LISA resolvable sources

Impact of a primordial gravitational wave background on LISA resolvable sources

Uncertainty of the white dwarf astrophysical gravitational wave background

The astrophysical gravitational wave background (AGWB) is a stochastic gravitational wave (GW) signal emitted by different populations of in-spiralling binary systems containing compact objects throughout the Universe . In the frequency range between 10$^ $ and 10$^ $ hertz (Hz), it will be detected by future space-based gravitational wave detectors, such as Laser Interferometer Space Antenna (LISA). In a recent work, we concluded that the white dwarf (WD) contribution to the AGWB dominates that of black holes (BHs) and neutron stars (NSs) We aim to investigate the uncertainties of the WD AGWB that arise from the use of different stellar metallicities, star formation rate density (SFRD) models, and binary evolution models. We used the code we previously developed to determine the WD component of the AGWB. We used a metallicity-dependent SFRD based on an earlier work to construct five different SFRD models. We used four different population models based on a range of common-envelope treatments and six different metallicities for each model. For all possible combinations, the WD component of the AGWB is dominant over other populations of compact objects. The effects of metallicity and population model are less significant than the effect of a (metallicity dependent) SFRD model. We find a range of about a factor of 5 in the level of the WD AGWB around a value of $ WD $ at 1 mHz and a shape that is weakly dependent on the model. We find the uncertainty for the WD component of the AGWB to be about a factor of 5. We note that there are other uncertainties that have an effect on this signal as well. We discuss whether the turnover of the WD AGWB at 10 mHz will be detectable by LISA and find it to be likely. We confirm our previous findings asserting that the WD component of the AGWB dominates over other populations, in particular, BHs.

Projections of the uncertainty on the compact binary population background using popstock

The LIGO-Virgo-KAGRA collaboration has announced the detection to date of almost 100 binary black holes that have been used in several studies to infer the features of the underlying binary black hole population. From these objects it is possible to predict the overall gravitational-wave (GW) fractional energy density contributed by black holes throughout the Universe, and thus estimate the gravitational-wave background (GWB) spectrum emitted in the current GW detector band. These predictions are fundamental in our forecasts for background detection and characterisation, with both present and future instruments. The uncertainties in the inferred population strongly impact the predicted energy spectrum, and in this paper we present a new flexible method to quickly calculate the energy spectrum for varying black hole population features, such as the mass spectrum and redshift distribution. We have implemented this method in an open-access package popstock and extensively tested its capabilities. Using popstock we investigated how uncertainties in these distributions impact our detection capabilities, and present several caveats for background estimation. In particular, we find that the standard assumption that the background signal follows a two-thirds power law at low frequencies is both waveform and mass-model dependent, and that the power-law signal is likely shallower than previously modelled, given the current waveform and population knowledge.

Efficient Computation of Overlap Reduction Functions for Pulsar Timing Arrays.

Pulsar timing arrays seek and study gravitational waves (GWs) through the angular two-point correlation function of timing residuals they induce in pulsars. The two-point correlation function induced by the standard transverse-traceless GWs is the famous Hellings-Downs curve, a function only of the angle between the two pulsars. Additional polarization modes (vector or scalar) that may arise in alternative-gravity theories have different angular correlation functions. Furthermore, anisotropy, linear, or circular polarization in the stochastic GW background gives rise to additional structure in the two-point correlation function that cannot be written simply in terms of the angular separation of the two pulsars. In this Letter, we provide a simple formula for the most general two-point correlation function-or overlap reduction function (ORF)-for a gravitational-wave background with an arbitrary polarization state, possibly containing anisotropies in its intensity and polarization (linear and/or circular). We provide specific expressions for the ORFs sourced by the general-relativistic transverse-traceless GW modes as well as vector (or spin-1) modes that may arise in alternative-gravity theories.

Early-type galaxy speciation: elliptical (E) and ellicular (ES) galaxies in the Mbh–M*,sph diagram, and a merger-driven explanation for the origin of ES galaxies, antitruncated stellar discs in lenticular (S0) galaxies, and the Sérsicification of E galaxy light profiles

ABSTRACT In a recent series of papers, supermassive black holes were used to discern pathways in galaxy evolution. By considering the black holes’ coupling with their host galaxy’s bulge/spheroid, the progression of mass within each component has shed light on the chronological sequence of galaxy speciation. Offsets between the galaxy-morphology-dependent $M_{\rm bh}$–$M_{\rm \star ,sph}$ scaling relations trace a pattern of ‘punctuated equilibrium’ arising from merger-driven transitions between galaxy types, such as from spirals to dust-rich lenticulars and further to ‘ellicular’ and elliptical galaxies. This study delves deeper into the distinction between the ellicular galaxies – characterised by their intermediate-scale discs – and elliptical galaxies. Along the way, it is shown how some antitruncated large-scale discs in lenticular galaxies can arise from the coexistence of a steep intermediate-scale disc and a relatively shallow large-scale disc. This observation undermines application of the popular exponential-disc plus Sérsic-bulge model for lenticular galaxies and suggests some past bulge mass measurements have been overestimated. Furthermore, it is discussed how merger-driven disc-heating and blending likely leads to the spheroidalisation of discs and the conglomeration of multiple discs leads to the (high-n) Sérsicification of light profiles. The ellicular and elliptical galaxy distribution in the $M_{\rm bh}$–$M_{\rm \star ,sph}$ diagram is explored relative to major-merger-built lenticular galaxies and brightest cluster galaxies. The super-quadratic $M_{\rm bh}$–$M_{\rm \star }$ relations, presented herein, for merger-built systems should aid studies of massive black hole collisions and the gravitational wave background. Finally, connections to dwarf compact elliptical and ultracompact dwarf galaxies, with their 100–1000 times higher $M_{\rm bh}/M_{\rm \star ,sph}$ ratios, are presented.

Electromagnetic signatures from accreting massive black hole binaries in time domain photometric surveys

We study spectral and time variability of accreting massive black hole binaries (MBHBs) at milli-parsec separations surrounded by a geometrically thin circumbinary disc. To this end, we present the first computation of the expected spectral energy distribution (SED) and light curves (LCs) from 3D hyper-Lagrangian resolution hydrodynamic simulations of these systems. We modelled binaries with a total mass of $10^6$ M$_ eccentricities of $e=0,\, 0.9$, and a mass ratio of $q=0.1,\, 1$. The circumbinary disc has an initial aspect ratio of 0.1, features an adiabatic equation of state, and evolves under the effect of viscous heating, black-body cooling, and self gravity. To construct the SED, we considered black-body emission from each element of the disc and we added a posteriori an X-ray corona with a luminosity proportional to that of the mini-discs that form around each individual black hole. We find significant variability of the SED, especially at high energies, which translates into LCs displaying distinctive modulations of a factor of $ 2$ in the optical and of $ 10$ in UV and X-rays. We analysed in detail the flux variability in the optical band that will be probed by the Vera Rubin Observatory (VRO). We find clear modulations on the orbital period and half of the orbital period in all systems. Only in equal-mass binaries, we find an additional, longer-timescale modulation, associated with an over-density forming at the inner edge of the circumbinary disc (commonly referred to as a lump). When considering the VRO flux limit and nominal survey duration, we find that equal-mass, circular binaries are unlikely to be identified, due to the lack of prominent peaks in their Fourier spectra. Conversely, unequal-mass and/or eccentric binaries can be singled out up to $z 0.5$ (for systems with bol $ erg s$^ $) and $z 2$ (for systems with bol $ erg s$^ $). Identifying electromagnetic signatures of MBHBs at separations of $ $ pc is of paramount importance to understand the physics of the gravitational wave (GW) sources of the future Laser Interferometer Space Antenna, and to pin down the origin of the GW background (GWB) observed in pulsar timing arrays.

Ultralight dark matter explanation of NANOGrav observations

The angular correlation of pulsar residuals observed by NANOGrav and other pulsar timing array collaborations show evidence in support of the Hellings-Downs correlation expected from stochastic gravitational wave background (SGWB). In this paper, we offer a nongravitational wave explanation of the observed pulsar timing correlations as caused by an ultralight Lμ−Lτ gauge boson dark matter (ULDM). ULDM can affect the pulsar correlations in two ways. The gravitational potential of vector ULDM gives rise to a Shapiro time delay of the pulsar signals and a nontrivial angular correlation (as compared to the scalar ULDM case). In addition, if the pulsars have a nonzero charge of the dark matter gauge group, then the electric field of the local dark matter causes an oscillation of the pulsar and a corresponding Doppler shift of the pulsar signal. We point out that pulsars carry a significant charge of muons, and thus the Lμ−Lτ vector dark matter contributes to both the Doppler oscillations and the time delay of the pulsar signals. The synergy between these two effects provides a better fit to the shape of the angular correlation function, as observed by the NANOGrav collaboration, compared to the standard SGWB explanation or the SGWB combined with time delay explanations. Our analysis shows that, in addition to the SGWB signal, there may potentially be excess timing residuals attributable to the Lμ−Lτ ULDM. Published by the American Physical Society 2024

Theoretical probes of Higgs boson-axion nonperturbative couplings

In this work we investigate the phenomenological implications of several nontrivial Higgs-boson- axion couplings, which cover most of the possible nonperturbative scenarios. Specifically we consider the combination of having higher-order nonrenormalizable Higgs-boson-axion couplings originating from a weakly coupled effective theory combined with nonperturbative couplings of the form ∼εΛc2|H|2cos(ϕfa). Since we consider the misalignment axion, the nonperturbative couplings can be expanded in the form of a perturbation expansion in powers of ϕ/fa, thus after the electroweak symmetry breaking, the effective potential of the axion is drastically affected by these terms. We investigate the phenomenological implications of these terms for various values of the mass scale Λc, and some scenarios are theoretically disfavored, while other scenarios with nonperturbative Higgs-boson-axion couplings of the form ∼εΛc2|H|2cos(ϕfa) with Λc∼10−10×ma and ma∼10−10 eV, lead to a characteristic pattern in the stochastic gravitational wave background via the deformation of the background equation of state parameter occurring at frequencies probed by the Einstein Telescope. We also consider loop effects from the Higgs sector caused by the term ∼εΛc2|H|2cos(ϕfa) if we close the Higgs in one loop. Published by the American Physical Society 2024

Signatures of circumbinary disc dynamics in multimessenger population studies of massive black hole binaries

ABSTRACT We investigate the effect of the cutting-edge circumbinary disc (CBD) evolution models on massive black hole binary (MBHB) populations and the gravitational wave background (GWB). We show that CBD-driven evolution leaves a tell-tale signature in MBHB populations, by driving binaries towards an equilibrium eccentricity that depends on the binary mass ratio. We find high orbital eccentricities ($e_{\rm b} \sim 0.5$) as MBHBs enter multimessenger observable frequency bands. The CBD-induced eccentricity distribution of MBHB populations in observable bands is independent of the initial eccentricity distribution at binary formation, erasing any memory of eccentricities induced in the large-scale dynamics of merging galaxies. Our results suggest that eccentric MBHBs are the rule rather than the exception in upcoming transient surveys, provided that CBDs regularly form in MBHB systems. We show that the GWB amplitude is sensitive to CBD-driven preferential accretion onto the secondary, resulting in an increase in GWB amplitude $A_{\rm yr^{-1}}$ by over 100 per cent with just 10 per cent Eddington accretion. As we self-consistently allow for binary hardening and softening, we show that CBD-driven orbital expansion does not diminish the GWB amplitude, and instead increases the amplitude by a small amount. We further present detection rates and population statistics of MBHBs with $M_{\rm b} \gtrsim 10^6 \, {\rm M}_{\odot }$ in Laser Interferometer Space Antenna, showing that most binaries have equal mass ratios and can retain residual eccentricities up to $e_{\rm b} \sim 10^{-3}$ due to CBD-driven evolution.

Constraining the Origin of the Nanohertz Gravitational-wave Background by Pulsar Timing Array Observations of Both the Background and Individual Supermassive Binary Black Holes

The gravitational waves (GWs) from supermassive binary black holes (BBHs) have long been sought by pulsar timing array (PTA) experiments, in the forms of both a stochastic GW background (GWB) and individual sources. Evidence for a GWB was reported recently by several PTAs with origins to be determined. Here we use a BBH population synthesis model to investigate the detection probability of individual BBHs by the Chinese PTA (CPTA) and the constraint on the GWB origin that may be obtained by PTA observations of both GWB and individual BBHs. If the detected GWB signal is entirely due to BBHs, a significantly positive redshift evolution (∝ (1 + z)2.07) of the mass scaling relation between supermassive black holes and their host galaxies is required. In this case, we find that the detection probability of individual BBHs is ∼85% or 64% if using a period of 3.4 yr of CPTA observation data, with an expectation of ∼1.9 or 1.0 BBHs detectable with a signal-to-noise ratio ≥3 or 5, and it is expected to increase to >95% if the observation period is extended to 5 yr or longer. Even if the contribution from BBHs to the GWB power signal is as small as ∼10%, a positive detection of individual BBHs can still be expected within an observation period of ∼10 yr. A nondetection of individual BBHs within several years from now jointly with the detected GWB signal can put a strong constraint on the upper limit of the BBH contribution to the GWB signal and help identify/falsify a cosmological origin.