Pith. sign in

REVIEW 4 major objections 6 minor 1 cited by

Searching for extreme mass ratio inspirals in LISA: from identification to parameter estimation

T0 review · 4 major / 6 minor · reviewed 2026-08-07 · deepseek-v4-flash

Pith's one-line read A phase-blind time-frequency search statistic takes LISA EMRI signals from wide-prior detection to parameter estimation.

desk verdict A genuinely new EMRI search statistic and a working pipeline, but the validation is too narrow to support the 'wide priors' claim. read the letter →

arxiv 2505.17814 v1 pith:FOTJZGL6 submitted 2025-05-23 gr-qc astro-ph.IMphysics.data-anphysics.space-ph

classification gr-qcastro-ph.IMphysics.data-anphysics.space-ph
keywords extreme-mass-ratioinspiralsLISAgravitational-wavesearchtime-frequencyanalysisgeneralizedmeanparameterestimationMarkovchainMonteCarlodifferentialevolution
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper tries to establish that a single search statistic can take an extreme-mass-ratio inspiral (EMRI) signal in LISA data from a blind search over wide priors all the way to a Bayesian posterior. The statistic $S_\lambda$ compares the short-time Fourier power spectra of data and template with a generalized-mean exponent $\lambda$, discarding phase coherence between time-frequency pixels while preserving the coherence of the frequency evolution. This choice widens the primary peak of the search landscape, so a stochastic optimizer can land near the true intrinsic parameters instead of getting stuck on the nonlocal secondary modes that plague standard matched-filter searches. The authors demonstrate the full chain—initial power-match search, iterative prior narrowing, maximization of the ordinary signal-to-noise ratio, and MCMC sampling—on two simulated single-EMRI datasets with stationary Gaussian noise and a galactic foreground, recovering the injections within their 68% credible intervals except for the initial phase. If the claim holds, extracting and characterizing a moderately loud EMRI becomes a laptop-scale job of about a day, which matters because LISA will see many overlapping sources that must be identified and subtracted in a global analysis.

What carries the argument

The load-bearing object is the statistic $S_\lambda$ of Eq. 8 with the time-frequency inner product $\langle x,y\rangle_{\mathrm{tf},\lambda}$ of Eq. 9. This is a generalized mean (a type of average controlled by $\lambda$) of per-window spectral matches between data and template, computed from short-time Fourier transforms; it deliberately discards phase differences between time-frequency pixels while preserving the track of the signal's frequency evolution. The exponent $\lambda$ controls the tradeoff: small $\lambda$ gives a wide primary peak but lets noise or secondary modes win, while large $\lambda$ suppresses those but narrows the peak. The paper sets $\lambda=3$ as the smallest exponent for which, on the single test injection, secondary modes stay below the injected match. The rest of the machinery is the staged optimizer: initial search on intrinsic parameters without the LISA response, iterative prior narrowing, a final maximization of the ordinary SNR $\rho$, and an MCMC chain for posteriors.

What would settle it

Run the full pipeline on a set of EMRI injections drawn from the same wide prior without retuning $\lambda$; if the fraction of runs that converge to the primary peak and recover the injection within credible intervals is far below the two showcased cases, the general search claim is falsified. A cheaper check is to draw 1000 random parameter sets around a second, different injection and test whether $\lambda=3$ still keeps all secondary modes below the injected match.

Watch

Extended reading notes

Core claim

The central discovery is that maximizing a phase-blind, time-frequency power match is enough to pull an EMRI out of noisy LISA data. The statistic, defined as $S_\lambda=\langle d,s(\theta)\rangle_{\mathrm{tf},\lambda}/\sqrt{\langle s(\theta),s(\theta)\rangle_{\mathrm{tf},\lambda}}$, with the time-frequency inner product summing generalized-mean powers of per-window spectra, sacrifices phase coherence between short-time windows but keeps the coherence of the frequency evolution. This widens the primary mode in the mass, eccentricity, and semi-latus-rectum landscape relative to the standard matched-filter SNR, while a sufficiently large $\lambda$ suppresses spurious secondary modes. The paper's pipeline first optimizes $S_\lambda$ without the LISA response to locate the intrinsic parameters, then narrows the priors, refines with the full response, switches to maximizing the ordinary SNR $\rho$, and finally runs MCMC to obtain posteriors. In two demonstrations (optimal SNRs 56 and 82), the recovered maximum-likelihood signal matches the injection closely enough that the residual is an order of magnitude smaller than the signal, and the posteriors contain the injected parameters except for the initial phase.

Load-bearing premise

The method's key tuning choice—the exponent $\lambda=3$—is selected by looking at the known injected signal, so the demonstrations do not show that the same setting works for unknown signals; if $\lambda$ must be retuned per source, the claimed blind search capability collapses.

Editorial extensions

If this is right

  • A single loud EMRI can be carried from a wide-prior search to MCMC posteriors in about 10 hours of optimization plus 8 hours of sampling on a mobile GPU, making large template banks unnecessary in this regime.
  • The residual after subtracting the recovered signal is an order of magnitude smaller than the signal itself, so the output can be fed into a global fit that then searches for the remaining LISA sources.
  • Because the statistic only needs a time-domain waveform, it can be re-run with more accurate waveform models, including spinning black holes, without changing the search architecture.
  • The demonstrated SNR regime (optimal SNR 56 and 82) sets a floor: fainter EMRIs will need a modified statistic or a longer integration to be found by this pipeline.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The choice of $\lambda=3$ is made by inspecting the injected signal, so the paper does not establish a principled way to set this exponent for an unknown source; an adaptive or marginalized $\lambda$ would be a natural extension to test.
  • Since $S_\lambda$ ignores phase, the final coherent SNR step must carry the phase information, and the two demonstrations already show the initial phase $\Phi_{\varphi_0}$ falling outside the 68% region, suggesting phase parameters are the first to degrade at lower SNR.
  • The two-week prior window on the time to plunge means the method already assumes a rough end time for the inspiral; a fully blind search would need to scan such windows, multiplying the cost.
  • The paper leaves the behavior of $S_\lambda$ across the full EMRI parameter space unexplored, so the immediate stress test is an ensemble of injections with varied masses, eccentricities, sky locations, and SNRs using a fixed $\lambda$.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 6 minor

Summary. The paper introduces a new search statistic for EMRI signals in LISA data, combining a time-frequency power-matched generalized mean (Eq. 8) with a staged optimization pipeline: first optimizing S without the LISA response on intrinsic parameters, then refining with S and finally maximizing the standard SNR, followed by MCMC parameter estimation. The method is demonstrated on two injected signals (optimal SNR 56 and 81) using the FastEMRIWaveforms package and the fast LISA response, with reported end-to-end runtimes of about a day on a laptop. The abstract claims that the strategy allows EMRI signals to be found 'from wide priors all the way to performing parameter estimation'.

Significance. If the central claim were fully supported, this would be a practical advance: it is one of the first demonstrations of an end-to-end chain from a coarse search to posterior sampling for a single EMRI, and it makes concrete use of public waveform and response tools. The paper is honest about its computational cost and provides reproducible ingredients (code, waveform models, and detailed pseudocode). However, the evidence presented is currently limited to two high-SNR injections, and the key hyperparameter λ is selected using the same injection used for validation. The 'wide priors' claim is also weakened by the narrow two-week tp prior. The method is therefore best viewed as a promising pipeline with an existence proof, not a validated search-from-wide-priors capability.

major comments (4)
  1. [Sec. II.C, Eq. (8), Fig. 3] The selection of λ = 3 is post-hoc: the paper states that λ is chosen as 'the smallest value for which secondary modes do not exceed the match statistic of the injected parameters,' using the same injection (Table I) that is later used for the demonstrations in Sec. III. This is circular: the success of the two recoveries does not establish that λ = 3 is appropriate for other EMRI parameters, other noise realizations, or other prior volumes. Please either (a) calibrate λ on a set of injections that are not used in the final demonstrations, including varied intrinsic parameters and noise draws, (b) perform a sensitivity study showing the dependence of the results on λ, or (c) explicitly reframe the claim as an existence proof for a fixed, injection-tuned choice. As it stands, the abstract's general claim of a new search strategy that works from wide priors is not supported.
  2. [Sec. II.E, Table II, Sec. III] The claim of searching 'from wide priors' is only partially demonstrated. The prior on the time to plunge tp is restricted to [0.42, 0.46] yr, a two-week window centered on the true injected value tp = 0.44 yr (Table II). The text in Sec. II.E states that 'the first iteration of the search restricts the prior on the time to plunge, tp, to a two-week window,' and Sec. III explains that either the pipeline is repeated with a moving window or a previous detection algorithm has already localized the signal in time. This means the search is not blind over the full mission duration. Please specify how the moving-window scan is implemented (step size, number of windows, total cost) and, ideally, demonstrate a search over a longer time interval without prior knowledge to substantiate the 'wide priors' claim.
  3. [Sec. III, Sec. II.C] There is no detection statistic with a calibrated threshold or false-alarm analysis. The text in Sec. III says that 'if the identified signal has an SNR above a certain threshold the signal is determined as a detection,' but no threshold is defined, no false-alarm probability is computed, and no noise-only trials are reported. The histograms in Fig. 3 show the distribution of the statistic for random prior samples, but they do not provide a detection threshold over many noise realizations. Without such an analysis, the pipeline cannot be used to claim that it can 'find' signals in unknown data rather than recover a known injection. Please add a noise-only study and define a detection criterion.
  4. [Sec. III, Conclusions] The two demonstrations use only high-SNR injections (ρoptimal = 56 and 81), and the paper's own Conclusions state that 'a current limitation of the method is its reliance on relatively high-SNR signals that stand out within a noisy power spectrum.' This is a load-bearing caveat for the central claim: if the method only works for very loud signals, then it is not a general EMRI search strategy. Please characterize the minimum SNR at which the pipeline can still recover a signal, or at least report the behavior of the statistic as a function of SNR, to place the two examples in context. Alternatively, temper the abstract's claim to 'high-SNR EMRI signals.'
minor comments (6)
  1. [Introduction] In the sentence 'Previous works have adopted ... with phenomenological waveform seachers [17],' the word 'seachers' should be 'searchers.'
  2. [Acknowledgements] The phrase 'previous work on ERMI signal simulation' should read 'EMRI signal simulation'.
  3. [Fig. 2] The caption of Fig. 2 does not specify which value of λ corresponds to each panel, making it difficult to verify the statements in the text that 'S λ=1 exhibits secondary modes' in the bottom-left panel and that 'S λ=4 results in a narrower primary mode compared to S λ=3.' Please label each panel with its λ value.
  4. [Algorithm 1] The pseudocode reads 'F unctionEMRI search(d)' with an apparent space inserted; this should be 'Function EMRI_search(d).' Also, the comment about the two-week window in Algorithm 1 is redundant with the main text and could be clarified.
  5. [Sec. III.A] The sentence 'Notably, the detected SNR of the recovered signal slightly surpasses that of the injected signal, which can be attributed to noise fluctuations. The condition ρ(θMLE) > ρ(θinjected) further validates the success of the optimization procedure' is somewhat overstated: for a noisy realization, the maximum matched-filter SNR is generally expected to exceed the injection SNR even for a correct template, so this condition does not by itself validate the optimization. It would be more informative to compare the recovered parameters to the injected values.
  6. [Appendix A] The sentence 'The full posterior Figure 11 with all parameters can be found in the appendix A' is grammatically awkward and should be rephrased, e.g., 'The full posterior, including all parameters, is shown in Figure 11 in Appendix A.'

Circularity Check

1 steps flagged · score 5.0 of 10

λ=3 is selected using the same injected signal used for the recovery demonstrations, making the headline results in-sample rather than independent predictions.

  1. fitted input called prediction [Sec. II.C (Matching spectrograms), Fig. 3 and following text; applied in Sec. III.A, Table III]
    "To balance these effects, we set λ = 3 for the remainder of this study, as it is the smallest value for which secondary modes do not exceed the match statistic of the injected parameters."

    The free exponent λ in the new search statistic S_λ is chosen by evaluating S_λ on 1000 random prior draws against the single injected EMRI used for all subsequent validation, selecting the smallest λ such that no random draw beats the match statistic of the injected parameters. The pipeline in Sec. III.A is then run on that same injection and reported as a successful extraction. The recovery therefore demonstrates that the optimizer can climb the landscape after the landscape's key hyperparameter was tuned to make the injected point stand out; it is not an out-of-sample test of the search statistic. The second example keeps the same λ=3, so it inherits the tuning choice, and the paper notes that an extensive investigation over the EMRI parameter space is left for future work.

full rationale

The central circularity is the selection of λ=3 in Sec. II.C using the same injected signal whose recovery is then presented as the main result (Sec. III.A, Table III). Because λ controls the width and prominence of the primary peak of S_λ, choosing it so that the injected parameters have the highest match among random prior samples builds the target recovery into the search statistic's tuning. This is a fitted-input-called-prediction pattern and justifies a score around 5. No other circular steps were found: the likelihood, SNR, time-frequency inner product, noise estimation, and MCMC sampling are defined independently and are not reduced to the paper's own conclusions; the waveform and response models are external tools. The paper's 'wide priors' claim is weakened by the two-week tp window (Sec. II.E, Table II) and the absence of a false-alarm or detection-threshold study, but that is a scope/overclaim issue rather than a circularity. The second injection (Table IV) provides some independent evidence, but it still uses the λ chosen for the first signal, so it does not remove the in-sample tuning concern.

Assumptions & free parameters 4 free parameters · 4 assumptions · 0 invented entities

The central claim rests on one tuned hyperparameter (λ), a narrow time-to-plunge prior, and several domain assumptions about noise and waveform fidelity. No new physical entities are introduced.

free parameters (4)
  • λ (generalized mean exponent) = 3
    Chosen as the smallest exponent such that, for 1000 random prior samples, no secondary mode exceeds the match statistic of the injected parameters (Sec. II.C, Fig. 3).
  • tp prior window = [0.42, 0.46] yr
    Two-week window centered on the injected plunge time 0.44 yr; restricts the search in time, undermining the 'wide priors' claim.
  • STFT window and overlap = 50000/dt samples, 50% overlap, Hann window
    Analysis settings chosen manually to optimize time-frequency resolution; no sensitivity study is provided.
  • Noise smoothing parameters = 30-bin moving median, Savitzky-Golay filter
    Chosen to smooth the PSD estimate; the values are fixed by hand and not varied.
assumptions (4)
  • domain assumption TDI channels A and E are independent with stationary Gaussian noise
    Used for the likelihood (Eq. 3) and the search statistic; real LISA noise may be non-stationary and correlated.
  • domain assumption The FastEMRIWaveforms (few) model accurately represents EMRI signals for injection and recovery
    Both injection and recovery use the same approximate waveform model, so systematic waveform errors are not tested.
  • domain assumption The PSD estimated from the data (including the signal) is a good proxy for the true noise PSD
    Noise estimation is performed on the same data containing the loud EMRI; the paper reports 3% relative error but uses this estimate in the likelihood.
  • standard math Differential evolution converges to the global maximum of the search statistics
    The optimization is stochastic; no guarantee of global convergence is provided.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Searching for extreme mass ratio inspirals in LISA: from identification to parameter estimation." pith.science (2026). https://pith.science/paper/FOTJZGL6

@misc{pith2026250517814,
  author       = {Pith},
  title        = {Pith review of: Searching for extreme mass ratio inspirals in LISA: from identification to parameter estimation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/FOTJZGL6}},
  note         = {Machine review of arXiv:2505.17814}
}
read the original abstract

The Laser Interferometer Space Antenna (LISA) is a planned space-based observatory designed to detect gravitational waves (GWs) within the millihertz frequency range. LISA is anticipated to observe the inspiral of compact objects into black holes at the centers of galaxies, so called extreme-mass-ratio inspirals (EMRIs). However, the extraction of these long-lived complex signals is challenging due to the large size and multimodality of the search space. In this study, we introduce a new search strategy that allows us to find EMRI signals in noisy data from wide priors all the way to performing parameter estimation. This work is an important step in understanding how to extract EMRIs from future LISA data.

Figures

Figures reproduced from arXiv: 2505.17814 by the authors.

Figure 1
Figure 1. FIG. 1: Spectrogram of the TDI A channel of the noisy [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: Search functions computed varying only two parameters at a time. The left panels show the SNR function, [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3: Histogram of different detection statistics for [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figures from the paper (9 more)
Figure 4
Figure 4. Figure 4: FIG. 4: Noise estimation for the data shown in Fig. [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5: The data, injected and recovered signal in the [PITH_FULL_IMAGE:figures/full_fig_p006_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6: The search chain of optimizing [PITH_FULL_IMAGE:figures/full_fig_p007_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7: The search chain of optimizing the SNR [PITH_FULL_IMAGE:figures/full_fig_p007_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8: Green solid lines denote the injected [PITH_FULL_IMAGE:figures/full_fig_p008_8.png]
Figure 9
Figure 9. Figure 9: FIG. 9: Search for the increased mass signal. The resulting chain optimizes [PITH_FULL_IMAGE:figures/full_fig_p009_9.png]
Figure 10
Figure 10. Figure 10: FIG. 10: Search for the increased mass signal, performed by optimizing the SNR [PITH_FULL_IMAGE:figures/full_fig_p009_10.png]
Figure 11
Figure 11. Figure 11: FIG. 11: The green solid lines indicate the injected parameters, while the orange dashed lines represent the recovered [PITH_FULL_IMAGE:figures/full_fig_p011_11.png]
Figure 12
Figure 12. Figure 12: FIG. 12: Corner plot of the increased mass signal. The injected parameters are shown as green solid lines, while the [PITH_FULL_IMAGE:figures/full_fig_p012_12.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Coherent End-to-End Search for Generic Extreme-Mass-Ratio Inspirals

    gr-qc 2026-08 conditional novelty 7.0 of 10

    A hierarchical search that uses the clustering of high-likelihood secondary maxima recovers generic EMRI signals in 14-dimensional parameter space from simulated LISA data.

Reference graph

Works this paper leans on

46 extracted references · 23 canonical work pages · cited by 1 Pith paper

  1. [21]

    We introduce a novel search statistic and employ it to achieve an accurate identification of EMRI signals, leading to the successful extraction of such sig- nals

    by reaching the primary peak and performing pa- rameter estimation, simulating a full search and param- eter estimation pipeline for a single EMRI in stationary arXiv:2505.17814v1 [gr-qc] 23 May 2025 2 Gaussian noise. We introduce a novel search statistic and employ it to achieve an accurate identification of EMRI signals, leading to the successful extrac...

  2. [1]

    Babak, J

    S. Babak, J. Gair, A. Sesana, E. Barausse, C. F. Sop- uerta, C. P. L. Berry, E. Berti, P. Amaro-Seoane, A. Pe- titeau, and A. Klein, Science with the space-based inter- ferometer LISA. V: Extreme mass-ratio inspirals, Phys. Rev. D 95, 103012 (2017), arXiv:1703.09722 [gr-qc]

  3. [2]

    Colpi, K

    M. Colpi, K. Danzmann, M. Hewitson, K. Holley- Bockelmann, P. Jetzer, G. Nelemans, A. Petiteau, D. Shoemaker, C. Sopuerta, R. Stebbins, et al., Lisa definition study report, arXiv preprint arXiv:2402.07571 (2024)

  4. [3]

    Amaro-Seoane, Relativistic dynamics and extreme mass ratio inspirals, Living Rev

    P. Amaro-Seoane, Relativistic dynamics and extreme mass ratio inspirals, Living Rev. Rel. 21, 4 (2018), arXiv:1205.5240 [astro-ph.CO]

  5. [4]

    C. P. L. Berry, S. A. Hughes, C. F. Sopuerta, A. J. K. Chua, A. Heffernan, K. Holley-Bockelmann, D. P. Mi- haylov, M. C. Miller, and A. Sesana, The unique poten- tial of extreme mass-ratio inspirals for gravitational-wave astronomy (2019), arXiv:1903.03686 [astro-ph.HE]

  6. [5]

    J. R. Gair, L. Barack, T. Creighton, C. Cutler, S. L. Larson, E. S. Phinney, and M. Vallisneri, Event rate estimates for LISA extreme mass ratio capture sources, Class. Quant. Grav. 21, S1595 (2004), arXiv:gr- qc/0405137

  7. [6]

    Speri, S

    L. Speri, S. Barsanti, A. Maselli, T. P. Sotiriou, N. War- burton, M. van de Meent, A. J. K. Chua, O. Burke, and J. Gair, Probing fundamental physics with Extreme Mass Ratio Inspirals: a full Bayesian inference for scalar charge (2024), arXiv:2406.07607 [gr-qc]

  8. [7]

    Barack and C

    L. Barack and C. Cutler, Using LISA EMRI sources to test off-Kerr deviations in the geometry of massive black holes, Phys. Rev. D 75, 042003 (2007), arXiv:gr- qc/0612029

Show all 46 references
  1. [8]

    C. F. Sopuerta and N. Yunes, Extreme and Intermediate- Mass Ratio Inspirals in Dynamical Chern-Simons Modified Gravity, Phys. Rev. D 80, 064006 (2009), arXiv:0904.4501 [gr-qc]

  2. [9]

    J. R. Gair, C. Tang, and M. Volonteri, LISA extreme- mass-ratio inspiral events as probes of the black hole mass function, Phys. Rev. D 81, 104014 (2010), arXiv:1004.1921 [astro-ph.GA]

  3. [10]

    Speri, A

    L. Speri, A. Antonelli, L. Sberna, S. Babak, E. Ba- rausse, J. R. Gair, and M. L. Katz, Probing Accretion Physics with Gravitational Waves, Phys. Rev. X 13, 021035 (2023), arXiv:2207.10086 [gr-qc]

  4. [11]

    P. S. Cole, G. Bertone, A. Coogan, D. Gaggero, T. Kary- das, B. J. Kavanagh, T. F. M. Spieksma, and G. M. Tomaselli, Distinguishing environmental effects on binary black hole gravitational waveforms, Nature Astron. 7, 943 (2023), arXiv:2211.01362 [gr-qc]

  5. [12]

    A. J. K. Chua and C. J. Cutler, Nonlocal parameter degeneracy in the intrinsic space of gravitational-wave signals from extreme-mass-ratio inspirals, Phys. Rev. D 106, 124046 (2022), arXiv:2109.14254 [gr-qc]

  6. [13]

    Babak, J

    S. Babak, J. G. Baker, M. J. Benacquista, N. J. Cornish, S. L. Larson, I. Mandel, S. T. McWilliams, A. Petiteau, E. K. Porter, E. L. Robinson, et al., The mock lisa data challenges: from challenge 3 to challenge 4, Classical and Quantum Gravity 27, 084009 (2010)

  7. [14]

    Babak, J

    S. Babak, J. R. Gair, and E. K. Porter, An Algorithm for detection of extreme mass ratio inspirals in LISA data, Class. Quant. Grav. 26, 135004 (2009), arXiv:0902.4133 [gr-qc]

  8. [15]

    N. J. Cornish, Detection Strategies for Extreme Mass Ratio Inspirals, Class. Quant. Grav. 28, 094016 (2011), arXiv:0804.3323 [gr-qc]

  9. [16]

    A. J. K. Chua, One-stop function for gravitational-wave detection, identification, and inference, Phys. Rev. D 106, 104051 (2022), arXiv:2205.08702 [gr-qc]. 11 FIG. 11: The green solid lines indicate the injected parameters, while the orange dashed lines represent the recovere...

  10. [17]

    Y. Wang, Y. Shang, S. Babak, Y. Shang, and S. Babak, EMRI data analysis with a phenomenological waveform, Phys. Rev. D 86, 104050 (2012), arXiv:1207.4956 [gr-qc]

  11. [18]

    Badger, J

    C. Badger, J. A. Font, M. Sakellariadou, and A. Torres- Forn´ e, High-speed reconstruction of long-duration grav- itational waves from extreme-mass-ratio inspirals using sparse dictionary learning, Phys. Rev. D 110, 064074 (2024), arXiv:2407.02908 [gr-qc]

  12. [19]

    Yun, W.-B

    Q. Yun, W.-B. Han, Y.-Y. Guo, H. Wang, and M. Du, The detection, extraction and parameter esti- mation of extreme-mass-ratio inspirals with deep learn- ing, Sci. China Phys. Mech. Astron. 68, 210413 (2025), arXiv:2311.18640 [gr-qc]

  13. [20]

    Zhang, C

    X.-T. Zhang, C. Messenger, N. Korsakova, M. L. Chan, Y.-M. Hu, and J.-d. Zhang, Detecting gravitational waves from extreme mass ratio inspirals using convolutional neural networks, Phys. Rev. D 105, 123027 (2022), arXiv:2202.07158 [astro-ph.HE]

  14. [22]

    Ye, H.-M

    C.-Q. Ye, H.-M. Fan, A. Torres-Orjuela, J.-d. Zhang, and Y.-M. Hu, Identification of gravitational waves from extreme-mass-ratio inspirals, Phys. Rev. D 109, 124034 (2024), arXiv:2310.03520 [gr-qc]

  15. [23]

    T. A. Prince, M. Tinto, S. L. Larson, and J. W. Arm- strong, Lisa optimal sensitivity, Physical Review D 66, 122002 (2002)

  16. [24]

    Vallisneri, Synthetic lisa: Simulating time delay inter- ferometry in a model lisa, Physical Review D 71, 022001 12 FIG

    M. Vallisneri, Synthetic lisa: Simulating time delay inter- ferometry in a model lisa, Physical Review D 71, 022001 12 FIG. 12: Corner plot of the increased mass signal. The injected parameters are shown as green solid lines, while the recovered parameters, θMLE, are indicated...

  17. [25]

    T. B. Littenberg, N. J. Cornish, K. Lackeos, and T. Rob- son, Global analysis of the gravitational wave signal from galactic binaries, Phys. Rev. D 101, 123021 (2020)

  18. [26]

    A. J. K. Chua, M. L. Katz, N. Warburton, and S. A. Hughes, Rapid generation of fully relativis- tic extreme-mass-ratio-inspiral waveform templates for LISA data analysis, Phys. Rev. Lett. 126, 051102 (2021), arXiv:2008.06071 [gr-qc]

  19. [27]

    M. L. Katz, L. Speri, A. J. K. Chua, C. E. A. Chapman- Bird, N. Warburton, and S. A. Hughes, BlackHolePertur- bationToolkit/FastEMRIWaveforms: Frequency Domain Waveform Added! (2023)

  20. [28]

    A. J. Chua, C. R. Galley, and M. Vallisneri, Reduced- order modeling with artificial neurons for gravitational- wave inference, Phys. Rev. Lett. 122, 211101 (2019), arXiv:1811.05491 [astro-ph.IM]. 13

  21. [29]

    Fujita and M

    R. Fujita and M. Shibata, Extreme mass ratio inspirals on the equatorial plane in the adiabatic order, Phys. Rev. D 102, 064005 (2020), arXiv:2008.13554 [gr-qc]

  22. [30]

    L. C. Stein and N. Warburton, Location of the last sta- ble orbit in Kerr spacetime, Phys. Rev. D 101, 064007 (2020), arXiv:1912.07609 [gr-qc]

  23. [31]

    A. J. Chua and J. R. Gair, Improved analytic extreme- mass-ratio inspiral model for scoping out eLISA data analysis, Class. Quant. Grav. 32, 232002 (2015), arXiv:1510.06245 [gr-qc]

  24. [32]

    A. J. Chua, C. J. Moore, and J. R. Gair, Augmented kludge waveforms for detecting extreme-mass-ratio inspi- rals, Phys. Rev. D 96, 044005 (2017), arXiv:1705.04259 [gr-qc]

  25. [33]

    Barack and C

    L. Barack and C. Cutler, LISA capture sources: Approx- imate waveforms, signal-to-noise ratios, and parameter estimation accuracy, Phys. Rev. D 69, 082005 (2004), arXiv:gr-qc/0310125

  26. [34]

    Speri, M

    L. Speri, M. L. Katz, A. J. Chua, S. A. Hughes, N. War- burton, J. E. Thompson, C. E. Chapman-Bird, and J. R. Gair, Fast and fourier: extreme mass ratio inspiral wave- forms in the frequency domain, Frontiers in Applied Mathematics and Statistics 9, 1266739 (2024)

  27. [35]

    A. J. Chua, M. L. Katz, N. Warburton, and S. A. Hughes, Data for fast emri waveforms, 10.5281/zenodo.3981654 (2020)

  28. [36]

    Ye, H.-M

    C.-Q. Ye, H.-M. Fan, A. Torres-Orjuela, J.-d. Zhang, and Y.-M. Hu, Identification of gravitational waves from extreme-mass-ratio inspirals, Physical Review D 109, 124034 (2024)

  29. [37]

    Katz, CChapmanbird, L

    M. Katz, CChapmanbird, L. Speri, N. Karnesis, and N. Korsakova, mikekatz04/lisaanalysistools: First main release. (2024)

  30. [38]

    de Carvalho, Mean, what do you mean? (2016)

    M. de Carvalho, Mean, what do you mean? (2016)

  31. [39]

    P. S. Bullen, The power means, in Handbook of Means and Their Inequalities(Springer Netherlands, Dordrecht,

  32. [40]

    S. H. Strub, L. Ferraioli, C. Schmelzbach, S. C. St¨ ahler, and D. Giardini, Global analysis of lisa data with galactic binaries and massive black hole binaries, Physical Review D 110, 024005 (2024)

  33. [41]

    S. H. Strub, L. Ferraioli, C. Schmelzbach, S. C. St¨ ahler, and D. Giardini, Accelerating global parameter estima- tion of gravitational waves from galactic binaries using a genetic algorithm and gpus, Physical Review D 108, 103018 (2023)

  34. [42]

    Karnesis, M

    N. Karnesis, M. L. Katz, N. Korsakova, J. R. Gair, and N. Stergioulas, Eryn: a multipurpose sampler for bayesian inference, Monthly Notices of the Royal Astro- nomical Society 526, 4814 (2023)

  35. [43]

    M. Katz, N. Karnesis, and N. Korsakova, mikekatz04/eryn: first full release (2023)

  36. [44]

    Foreman-Mackey, D

    D. Foreman-Mackey, D. W. Hogg, D. Lang, and J. Good- man, emcee: the mcmc hammer, Publications of the As- tronomical Society of the Pacific 125, 306 (2013)

  37. [45]

    Virtanen, R

    P. Virtanen, R. Gommers, T. E. Oliphant, M. Haber- land, T. Reddy, D. Cournapeau, E. Burovski, P. Peter- son, W. Weckesser, J. Bright, et al., SciPy 1.0: Fun- damental algorithms for scientific computing in Python, Nature Methods 17, 261 (2020)

  38. [46]

    Storn and K

    R. Storn and K. Price, Differential evolution–a simple and efficient heuristic for global optimization over con- tinuous spaces, Journal of global optimization 11, 341 (1997)

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.