The reviewed record of science sign in
Pith

arxiv: 2606.12952 · v1 · pith:URZ2ZOXL · submitted 2026-06-11 · astro-ph.HE

Gravitational wave background from extreme-mass-ratio inspirals

Haoyu Zhao , Yuanhao Zhang , Xilong Fan , Wenbiao Han This is my paper

Reviewed by Pith2026-06-27 06:11 UTCgrok-4.3pith:URZ2ZOXLopen to challenge →

classification astro-ph.HE
keywords gravitational wavesextreme mass ratio inspiralsLISAsupermassive black holescompact objectscharacteristic straingalactic centersstochastic background
0
0 comments X

The pith

The masses of compact objects determine the strength of the gravitational wave background from extreme-mass-ratio inspirals by up to an order of magnitude.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper calculates the expected gravitational wave background signal from populations of extreme-mass-ratio inspirals. It examines how the signal's characteristic strain depends on the distribution of sources, the spin of the central black holes, and especially the type of compact object involved. The results show that black hole companions produce a much stronger background than neutron stars or white dwarfs. This matters because the background will be a major signal for the LISA detector, and understanding its possible range helps in planning observations and interpreting data about galactic centers.

Core claim

By analyzing the astrophysical distribution of EMRI sources and key parameters including the spin of supermassive black holes and the masses of compact objects, the distribution range of the characteristic strain of the GWB from EMRIs is determined. The final eccentricity distributions have negligible effect on the intensity due to rapid circularization. The spin of the SMBH enhances the GW characteristic strain by approximately 1%. The masses of COs significantly affect the characteristic strain, with BH as CO producing a GW signal intensity approximately one order of magnitude higher than NS or WD cases.

What carries the argument

The calculation of the characteristic strain of the GWB from EMRIs based on the statistical properties of source distributions and parameter variations.

Load-bearing premise

The astrophysical distribution of EMRI sources and the specific calculation methods accurately capture the statistical properties and dynamical processes in galactic centers.

What would settle it

Detection of a GWB characteristic strain level that is inconsistent with the predicted range across all tested compact object masses, or an absence of the expected difference between black hole and other compact object cases.

Figures

Figures reproduced from arXiv: 2606.12952 by Haoyu Zhao, Wenbiao Han, Xilong Fan, Yuanhao Zhang.

Figure 1
Figure 1. Figure 1: e − f relation. The horizontal axis represents the logarithm of the ratio between the orbital frequency forb and the initial frequency f0 , while the vertical axis represents the eccentricity e . Different colors correspond to the evolution in the frequency domain for initial eccentricities e0 = 0.2, 0.3, 0.5, 0.7, 0.9, 0.99 , respectively. This shows that for a given initial eccentricity e0 , e decreases … view at source ↗
Figure 2
Figure 2. Figure 2: Comparison of the fitting function and the analytical function when e0 = 0.9.The red solid line represents the analytical e − f function, while the blue dashed line represents the interpolated fitting formula. We can see that the fit is relatively good in the low-frequency region, but there is some deviation at higher frequencies. where β = χ 2/3/σ0 , χ = forb/f0 , forb is the orbital frequency. σ(e0) = e … view at source ↗
Figure 3
Figure 3. Figure 3: The function g(n, e) evolves with eccentricity e for harmonic numbers n = 1, 2, 3, 4, 5, 6, where g(2, 0) = 1 when n = 2 and e = 0, representing a circular orbit. The graph illustrates that higher￾order harmonics contribute more significantly to the total gravitational wave radiation as eccentricity increases. 6 [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: The solid lines represent the distribution of Φn as a function of eccentricity e for harmonic orders n = 1, 2, 3, 4, 5, 6. As Φ lacks an analytical solution, it must be determined through numerical methods or polynomial interpolation. In the range 0 < e < 0.9, the interpolation provides a good fit, allowing the use of a fitting function to approximate the values and thus significantly reduce computational … view at source ↗
Figure 5
Figure 5. Figure 5: This figure shows two different SMBH mass distribution models. The blue dashed line represents the B12 model, which corresponds to an optimistic scenario, while the orange dashed line corresponds to the G10 model, representing a lower threshold. The histogram represents the mass distribution we fitted in model 1. Plotting histograms can help us verify the accuracy of the fitted data by providing a visual c… view at source ↗
Figure 6
Figure 6. Figure 6: This figure shows the variation of the GWB characteristic strain with frequency for three different EMRI models. Model 1 is the fiducial model, model 11 represents the pessimistic scenario, and Model 12 is the optimistic scenario. The pink shaded region represents the distribution range of hc. The black dashed line represents the LISA’s sensitivity curve (Tobs = 4 yr). . 12 [PITH_FULL_IMAGE:figures/full_f… view at source ↗
Figure 7
Figure 7. Figure 7: This figure illustrates the comparison between the final eccentricity distributions and the corresponding initial eccentricity e0 distributions for EMRIs, with final eccentricities restricted to the range 0–0.2. The solid curves represent the final eccentricity distributions after the systems have evolved under GW emission. Overlaid on these curves are histograms representing the initial e0 distributions: … view at source ↗
Figure 8
Figure 8. Figure 8: Evolution of orbital eccentricity e(t) over 4 years for EMRI systems with different initial eccentricities. The blue, orange, and green curves correspond to initial eccentricities e0 ≈ 0.01, e0 = 0.5, and e0 ≈ 0.9, respectively. Systems with very high e0 exhibit rapid circularization within a short timescale (∼ 0.01 yr), while those with intermediate or low initial eccentricities evolve more gradually. Thi… view at source ↗
Figure 9
Figure 9. Figure 9: Characteristic strain spectra hc(f) for EMRI systems with different initial eccentricities e0 = 0.01, 0.5, and 0.9 for single sources at z = 4.5, shown against the LISA sensitivity curve (dashed line). Despite the variation in initial eccentricity, the resulting strain curves lie close to each other across the LISA-sensitive frequency band, with differences of at most a few factors. 15 [PITH_FULL_IMAGE:fi… view at source ↗
Figure 10
Figure 10. Figure 10: This figure shows the characteristic strain of the GWB corresponding to two different final eccentricity distributions. The blue solid line represents the case where the final eccentricity follows a uniform distribution, ranging from 0 to 0.2, while the red dashed line represents the case where the final eccentricity follows a power-law distribution with k = -0.9, also within the range of 0 to 0.2. It can… view at source ↗
Figure 11
Figure 11. Figure 11: The blue solid line represents the case where the central black hole has spin, with a spin parameter a = 0.98 , while the red dashed line represents the non-spinning case. It can be observed that in the low-frequency range, the effect of spin is not significant, but in the high-frequency range, a rapidly spinning black hole significantly enhances the GW radiation intensity. 17 [PITH_FULL_IMAGE:figures/fu… view at source ↗
Figure 12
Figure 12. Figure 12: This figure illustrates the impact of different compact objects on the gravitational wave signal. The term “CO” refers to all compact objects, which include black holes (BH), neutron stars (NS), and white dwarf-neutron star (WD-NS) systems, with masses in the range of 1–50 M⊙ . When the compact object is a BH , whose mass range is 5–50 M⊙, the gravitational wave signal is the strongest, especially in the … view at source ↗
read the original abstract

The gravitational wave background (GWB) produced by extreme-mass-ratio inspirals (EMRIs) serves as a powerful tool for probing the astrophysical and dynamical processes in galactic centers. EMRI systems are a primary target for the space-based detector LISA due to their long-lived signals and high signal-to-noise ratios. This study explores the statistical properties of the GWB from EMRI, focusing on the calculation methods for the GWB, the astrophysical distribution of EMRI sources, and the influence of key parameters, including the spin of supermassive black holes (SMBHs) and the masses of compact objects (COs). By analyzing these factors, we determine the distribution range of the characteristic strain of the GWB from EMRIs. We find that the final eccentricity distributions appear to have negligible effect on the intensity of the GWB due to rapid circularization before they become detectable and the spin of the SMBH enhances the GW characteristic strain by approximately 1$\%$ compared to cases without spin effects. The masses of COs can also significantly affect the characteristic strain of the GWB from EMRIs, with Black Hole (BH) as CO producing a GW signal intensity that is approximately one order of magnitude higher compared to cases where Neutron Star (NS) or White Dwarf (WD) are the COs.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper computes the stochastic gravitational-wave background (GWB) from extreme-mass-ratio inspirals (EMRIs) for LISA, focusing on the effects of source eccentricity distributions, supermassive black hole (SMBH) spin, and compact-object (CO) mass/type on the characteristic strain. It reports that final eccentricity distributions have negligible impact due to rapid circularization, SMBH spin increases the strain by ~1%, and black-hole COs produce ~10× higher strain than neutron-star or white-dwarf COs, with the range of strains determined by the adopted astrophysical source distributions and rates.

Significance. If the central numerical claims hold after correction for mass-dependent rates, the work would supply concrete guidance on the dominant astrophysical uncertainties in the EMRI GWB for LISA data analysis. The explicit comparison across CO types is a useful addition to the literature, but only if the underlying source-rate model is shown to be self-consistent with dynamical capture physics.

major comments (2)
  1. [Abstract / results on CO mass] Abstract (and presumably the results section presenting the order-of-magnitude claim): the statement that BH COs produce approximately one order of magnitude higher characteristic strain than NS/WD COs is presented as a robust outcome of the astrophysical distribution. However, the EMRI formation rate Γ must depend on m_CO through mass segregation, dynamical friction, and relaxation timescales; if the calculation normalizes the same total rate or density across CO types rather than using an m-dependent ρ(m) or τ_relax(m), the quoted factor of 10 does not reflect the integrated background. This directly undermines the central claim about CO-mass effects.
  2. [Methods] Methods section (as flagged by the absence of derivation steps, error bars, or exclusion criteria for the quoted 1% and order-of-magnitude numbers): the abstract supplies specific numerical outcomes without showing how the source distributions were constructed or how the GWB integral was evaluated. A full methods description is required to verify whether post-hoc choices in the rate normalization affect the reported differences.
minor comments (2)
  1. Notation for the characteristic strain h_c(f) and the precise definition of the GWB integral should be stated explicitly once in the text rather than assumed from prior literature.
  2. Figure captions should indicate whether the plotted curves correspond to fixed total rate or to m_CO-dependent capture rates.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript on the EMRI gravitational-wave background. The comments raise important points about the robustness of the CO-mass comparison and the level of methodological detail. We address each major comment below and have revised the manuscript to improve clarity and transparency.

read point-by-point responses

  1. Referee: [Abstract / results on CO mass] Abstract (and presumably the results section presenting the order-of-magnitude claim): the statement that BH COs produce approximately one order of magnitude higher characteristic strain than NS/WD COs is presented as a robust outcome of the astrophysical distribution. However, the EMRI formation rate Γ must depend on m_CO through mass segregation, dynamical friction, and relaxation timescales; if the calculation normalizes the same total rate or density across CO types rather than using an m-dependent ρ(m) or τ_relax(m), the quoted factor of 10 does not reflect the integrated background. This directly undermines the central claim about CO-mass effects.

    Authors: We agree that a mass-independent normalization would undermine the claim. Our adopted source distributions are taken from literature models that already incorporate mass segregation, dynamical friction, and relaxation timescales that differ by CO type, producing both higher per-source strain and higher effective rates for BH COs. The reported factor of ~10 therefore reflects the integrated background under those m_CO-dependent inputs. To make this explicit, we have added a dedicated paragraph in the methods section describing the rate model, citing the specific references for the m-dependent distributions, and stating that the same total rate was not imposed across CO types. revision: yes

  2. Referee: [Methods] Methods section (as flagged by the absence of derivation steps, error bars, or exclusion criteria for the quoted 1% and order-of-magnitude numbers): the abstract supplies specific numerical outcomes without showing how the source distributions were constructed or how the GWB integral was evaluated. A full methods description is required to verify whether post-hoc choices in the rate normalization affect the reported differences.

    Authors: We accept that the original methods section was insufficiently detailed for independent verification. We have expanded it to include: (i) the explicit form of the GWB integral and the numerical quadrature method used, (ii) the step-by-step construction of the eccentricity, spin, and mass distributions from the cited astrophysical models, and (iii) the precise normalization procedure together with any exclusion criteria applied to the Monte-Carlo realizations. These additions directly address concerns about post-hoc rate choices and allow readers to reproduce the quoted 1% and order-of-magnitude results. revision: yes

Circularity Check

0 steps flagged

No circularity: claims rest on external astrophysical inputs without reduction to self-fit

full rationale

The abstract presents results on GWB strain dependence on CO mass, SMBH spin, and eccentricity as outcomes of analyzing source distributions and calculation methods. No equations, self-citations, or derivations are provided that reduce a prediction to a fitted input by construction, nor any self-definitional loops or uniqueness theorems imported from the authors. The distribution of sources is treated as an input parameter whose effect is computed, not redefined from the output strain. This satisfies the default expectation of a non-circular paper when no load-bearing step can be quoted as equivalent to its own inputs.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 0 invented entities

Central claim depends on multiple unfixed astrophysical inputs (source rates, mass functions, spin distributions) and the assumption that standard EMRI waveform and population-synthesis models are adequate; no independent evidence for these inputs is supplied in the abstract.

free parameters (3)
  • EMRI source distribution and rate
    Astrophysical distribution of EMRI sources is varied to determine the range of characteristic strain.
  • Compact object mass and type
    Masses and types (BH, NS, WD) are treated as variable inputs that directly scale the reported signal intensity.
  • SMBH spin parameter
    Spin is introduced as a key parameter whose effect is quantified at the 1% level.
axioms (2)
  • domain assumption Rapid circularization occurs before EMRIs become detectable
    Invoked to conclude that final eccentricity distributions have negligible effect on GWB intensity.
  • domain assumption Standard general-relativistic waveform models apply to EMRIs
    Underlying all characteristic-strain calculations.

pith-pipeline@v0.9.1-grok · 5771 in / 1581 out tokens · 24905 ms · 2026-06-27T06:11:56.504175+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

29 extracted references · 3 canonical work pages · 1 internal anchor

  1. [1]

    2010 Class

    Amaro-Seoane P and Gair J R et al. 2010 Class. Quantum Grav. 24 R113

    2010

  2. [2]

    Quantum Grav

    Robson T and Cornish N J and Liu C 2019 Class. Quantum Grav. 36 105011

    2019

  3. [3]

    Laser Interferometer Space Antenna

    Amaro-Seoane P et al. 2017 arXiv 1702.00786

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  4. [4]

    Thorne K S 1980 Rev. Mod. Phys. 52 299

    1980

  5. [5]

    Rauch K P and Tremaine S 1996 New Astron. 1 149

    1996

  6. [6]

    Tagawa H and Haiman Z and Kocsis B 2020 Astrophys. J. 898 25

    2020

  7. [7]

    Maggiore M 2008 Gravitational Waves: Theory and Experiments (Oxford: Oxford University Press)

    2008

  8. [8]

    Renzini A I and Goncharov B and Jenkins A C and Meyers P M 2022 arXiv 2202.00178

    work page Pith review arXiv 2022

  9. [9]

    Christensen N 2019 Rep. Prog. Phys. 82 016903

    2019

  10. [10]

    Barausse E 2012 Mon. Not. R. Astron. Soc. 423 2533

    2012

  11. [11]

    Milosavljević M and Merritt D 2003 Astrophys. J. 596 860

    2003

  12. [12]

    Barack L and Cutler C 2007 Phys. Rev. D 75 042003

    2007

  13. [13]

    Quantum Grav

    Vallisneri M 2013 Class. Quantum Grav. 30 224010

    2013

  14. [14]

    Hughes S A 2001 Phys. Rev. D 64 064004

    2001

  15. [15]

    Quantum Grav

    Gair J R and Barack L and Creighton T and Cutler C and Larson S L and Phinney E S and Vallisneri M 2004 Class. Quantum Grav. 21 S1595

    2004

  16. [16]

    Enoki M and Nagashima M 2007 Prog. Theor. Phys. 117 241

    2007

  17. [17]

    Yunes N and Arun K G and Berti E and Will C M 2009 Phys. Rev. D 80 084001

    2009

  18. [18]

    Huerta E A and McWilliams S T and Gair J R and Taylor S R 2015 Phys. Rev. D 92 063010 20

    2015

  19. [19]

    Peters P C and Mathews J 1963 Phys. Rev. 131 435

    1963

  20. [20]

    Finn L S and Thorne K S 2000 Phys. Rev. D 62 124021

    2000

  21. [21]

    Babak S and Gair J and Sesana A and Barausse E and Sopuerta C F and Berry C P L and Berti E and Amaro- Seoane P and Petiteau A and Klein A 2017 Phys. Rev. D 95 103012

    2017

  22. [22]

    Quantum Grav

    Barack L 2009 Class. Quantum Grav. 26 213001

    2009

  23. [23]

    Phinney E S 2001 arXiv 0108028

    2001

  24. [24]

    Bonetti M and Sesana A 2020 Phys. Rev. D 102 103023

    2020

  25. [25]

    China Phys

    Li G L and Tang Y and Wu Y L 2022 Sci. China Phys. Mech. Astron. 65 109732

    2022

  26. [26]

    Qunbar I and Stone N C 2023 arXiv 2304.13062

    work page arXiv 2023

  27. [27]

    2009 Astrophys

    Gültekin K et al. 2009 Astrophys. J. 698 198221

    2009

  28. [28]

    Kormendy J and Ho L C 2013 Annu. Rev. Astron. Astrophys. 51 511653

    2013

  29. [29]

    Graham A W and Scott N 2013 Astrophys. J. 764 151 21

    2013