Stochastic Background Research Articles

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.

Explaining PTA results by metastable cosmic strings from SO(10) GUT

In a recent paper ( Phys. Rev. D 108 (2023) 095053), we have demonstrated that the 2023 PTA results, which hint at a stochastic gravitational wave (GW) background at nanohertz frequencies, point towards a promising model-building route for realizing SO(10) Grand Unification with embedded inflation. The proposed supersymmetric scenario solves the doublet-triplet splitting without fine-tuning, accounts for charged fermion and neutrino masses, avoids conflicts with current proton decay bounds, and includes only representations no larger than the adjoint. It features multi-step breaking of SO(10) to the Standard Model gauge symmetry, with inflation embedded such that metastable cosmic strings are produced at the end of inflation. This cosmic string network generates a stochastic GW background that can explain the PTA results. In this paper, we provide a detailed analysis of the singled out GUT model class, focusing on how the gauge coupling unification condition affects the scales of multi-step SO(10) breaking and the preferred GW spectra. The lowest breaking scale, linked to inflation, the generation of right-handed neutrino masses for the seesaw mechanism, and metastable cosmic string production, coincides with the range suggested by the PTA results.

A millisecond pulsar position determined to 0.2 mas precision with VLBI

Precise millisecond pulsar (MSP) positions determined with very long baseline interferometry (VLBI) hold the key to building the connection between the kinematic and dynamic reference frames respectively used for VLBI and pulsar timing. A frame connection would provide an important pathway to examining the planetary ephemerides used in pulsar timing, and would potentially enhance the sensitivities of the pulsar timing arrays used to detect stochastic gravitational-wave background in the nano-Hz regime. We aim to significantly improve the precision of the VLBI-based MSP position (from $ at present) by reducing the two dominant components in the positional uncertainty --- the propagation-related uncertainty and the uncertainty resulting from the frequency-dependent core shifts of the reference sources. We introduce a new differential astrometry strategy called PINPT (Phase-screen Interpolation plus frequeNcy-dePendent core shifT correction; pronounced ``pinpoint''), which entails the use of multiple calibrators observed at several widely separated frequencies. The strategy allows determination of the core shift and mitigates the impact of residual delay in the atmosphere. We implemented the strategy on an MSP that is well constrained astrometrically with VLBI and pulsar timing. Using the PINPT strategy, we determined core shifts for four AGNs around and derived a VLBI-based pulsar position with uncertainties of 0.17\,mas and 0.32\,mas in Right Ascension and Declination, respectively, approaching the uncertainty level of the best-determined timing-based MSP positions. Additionally, incorporating the new observations into historical ones, we refined the pulsar proper motion and the parallax-based distance to the $ level and the subparsec level, respectively. The realization of the PINPT strategy promises a factor-of-five positional precision enhancement (over conventional VLBI astrometry) for all kinds of compact radio sources observed at $ including most fast radio bursts.

Induced domain walls of QCD axion, and gravitational waves

We show that heavy axion domain walls induce domain walls of the QCD axion through a mixing between the heavy axion and the QCD axion, even when the pre-inflationary initial condition is assumed for the QCD axion. The induced domain walls arise because the effective θ parameter changes across the heavy axion domain walls, shifting the potential minimum of the QCD axion. When the heavy axion domain walls collapse, the induced QCD axion domain walls collapse as well. This novel mechanism for producing the QCD axions can explain dark matter even with the axion decay constant as small as 𝒪(109) GeV. In particular, this scenario requires domain wall collapse near the QCD crossover, potentially accounting for the stochastic gravitational wave background suggested by recent pulsar timing array observations, including NANOGrav. Using this mechanism, it is also possible to easily create induced domain walls for string axions or axions with a large decay constant, which would otherwise be challenging. We also comment on the implications for cosmic birefringence using induced axion domain walls.