Recognition: unknown
Chaotic migration of LISA Extreme Mass Ratio Inspirals in a turbulent accretion disk: effect on waveform de-phasing
Pith reviewed 2026-05-09 23:14 UTC · model grok-4.3
The pith
Turbulent gas torques on extreme mass ratio inspirals can produce observable dephasing in LISA gravitational wave signals when mean linear torques alone do not.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that when the turbulent torque is modeled as a Gaussian distribution around the linear torque T_lin with normalization C, the resulting chaotic migration of the EMRI produces a gas-induced dephasing Delta psi_gas that exceeds the detectability threshold of 8/SNR for f_Edd ≳ 0.3, C ≳ 300, h0 ≳ 0.03, and α ≳ 0.1 in the final four years of a golden circular EMRI with M=10^6 M_sun, q=5e-5 at z=0.276 and SNR=50, whereas considering only T_lin yields unobservable dephasings.
What carries the argument
The turbulent torque T_turb, represented as a Gaussian random variable centered on the linear torque T_lin with a normalization factor C drawn from hydrodynamical simulations, which introduces stochasticity that drives chaotic orbital migration and accumulates waveform phase shifts.
Load-bearing premise
The assumption that turbulent torques act as Gaussian fluctuations around the mean linear torque, with strength fixed from prior simulations, without self-consistently evolving the coupled disk-EMRI system.
What would settle it
Detection or non-detection of dephasing signals larger than 8/SNR in actual LISA EMRI events from high-accretion turbulent disks would confirm or refute the model predictions.
Figures
read the original abstract
Gravitational wave (GW) detector LISA will observe near-coalescence extreme mass ratio inspirals (EMRIs), which typically form in galactic central accretion disks. Gas torques on EMRI will alter its GW-driven inspiral trajectory from the vacuum expectation, leading to potentially LISA-observable GW dephasing ($\Delta\psi_{\rm gas}$). Most studies compute $\Delta\psi_{\rm gas}$ for a thin, laminar disk, with negligible flow turbulence, where the disk exerts a fairly well-understood linear torque ($T_{\rm lin}$). However, these disks must be turbulent due to magneto-rotational instability in the inner regions. Hence, we present a proof-of-concept general, agnostic prescription for the turbulent torque ($T_{\rm turb}$) acting on an EMRI by modeling it as a Gaussian distribution around $T_{\rm lin}$, based on recent advances from a global hydrodynamical (HD) study. We compute $\Delta\psi_{\rm gas}$ for the ``golden'' circular EMRI with total source mass $M=10^6~{\rm M}_\odot$ and mass ratio $q=5\times10^{-5}$ in its final four-year evolution at redshift $z=0.276$ and signal-to-noise ratio (SNR) $=50$ by varying Eddington ratio ${\rm f}_{\rm Edd}$, turbulence normalization $C$ ($=~360$ in the aforementioned HD study), disk aspect ratio $h_0$, and turbo-viscous coefficient $\alpha$ in a reasonable parameters space. We find that for ${\rm f}_{\rm Edd}\gtrsim0.3$, $C\gtrsim300$, $h_0\gtrsim0.03$, and $\alpha\gtrsim0.1$, gas-induced dephasings are unobservable if only considering $T_{\rm lin}$ but could become detectable ($\Delta\psi_{\rm gas}>8/$SNR) if EMRIs exhibit chaotic migration due to turbulent gas flow. Hence, this work motivates running MHD simulations of accretion disks with embedded LISA EMRIs in the early in-spiral phase over long enough timescales to understand the evolution of their orbital elements and the imprint of the turbulent environment on their gravitational waveforms.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a proof-of-concept model in which the turbulent torque T_turb on a LISA EMRI is drawn from a Gaussian distribution centered on the linear torque T_lin, with width set by a normalization C drawn from a prior hydrodynamical study. For a fiducial circular EMRI (M=10^6 M_⊙, q=5×10^{-5}, z=0.276, SNR=50), the authors integrate the orbital evolution over the final four years and report that gas-induced dephasing Δψ_gas exceeds the detectability threshold 8/SNR for f_Edd ≳ 0.3, C ≳ 300, h0 ≳ 0.03 and α ≳ 0.1, whereas the same parameters yield undetectable dephasing when only T_lin is used. The work concludes by motivating self-consistent MHD simulations of disks with embedded EMRIs.
Significance. If the Gaussian turbulence prescription can be shown to produce a net random walk that survives orbital and viscous correlations, the result would indicate that disk turbulence can push otherwise undetectable gas effects into the LISA band, directly affecting waveform modeling and parameter estimation. The numerical integration itself is straightforward and the parameter survey is systematic; the paper also correctly flags the need for future coupled simulations. However, the claimed detectability window is defined by the same free parameters that set the torque variance, limiting the strength of the prediction.
major comments (3)
- [§2] §2 (Turbulent torque prescription): The Gaussian model T_turb ~ N(T_lin, C) is introduced without derivation from the MHD equations or direct validation against simulation data that include an embedded point mass; C is taken from a single prior global HD run that does not contain the EMRI, leaving open whether the quoted variance remains appropriate once the secondary perturbs the local flow.
- [§3] §3 (Numerical integration): The torque time series is sampled independently at each timestep, implicitly assuming white, uncorrelated fluctuations; the manuscript does not evolve the disk and EMRI self-consistently, so possible correlations on orbital or viscous timescales are not tested and could cause the random walk in semi-major axis (and thus Δψ_gas) to saturate rather than accumulate to detectable levels.
- [§4] §4 (Results and detectability): The statement that Δψ_gas > 8/SNR holds for f_Edd ≳ 0.3, C ≳ 300, etc., is obtained after scanning the four-dimensional parameter space; the threshold is therefore applied post hoc rather than as a blind, a-priori prediction, making the central claim dependent on the chosen values of the free parameters.
minor comments (3)
- [§3] The definition of Δψ_gas and the precise numerical implementation of the phase integral should be stated explicitly in the text (currently referenced only to prior work).
- [Figures] Figure captions and axis labels could more clearly indicate which curves correspond to the laminar (T_lin only) versus turbulent cases and whether error bands reflect Monte-Carlo realizations.
- [§3] A short convergence test with respect to timestep or number of turbulence realizations would strengthen the numerical results section.
Simulated Author's Rebuttal
We thank the referee for the thoughtful and constructive report. We address each major comment below and have revised the manuscript to incorporate clarifications and additional caveats where appropriate. Our responses aim to strengthen the presentation of this proof-of-concept study while honestly acknowledging its limitations.
read point-by-point responses
-
Referee: §2 (Turbulent torque prescription): The Gaussian model T_turb ~ N(T_lin, C) is introduced without derivation from the MHD equations or direct validation against simulation data that include an embedded point mass; C is taken from a single prior global HD run that does not contain the EMRI, leaving open whether the quoted variance remains appropriate once the secondary perturbs the local flow.
Authors: We agree that the Gaussian prescription is a phenomenological model rather than a first-principles derivation from the MHD equations with an embedded secondary. The value of C is drawn from a prior global hydrodynamical simulation that did not include an EMRI, so the local flow perturbation by the secondary is not accounted for. This is an inherent limitation of the current approach. In the revised manuscript we have expanded the discussion in §2 to explicitly state these caveats, to clarify that the model is intended only as a proof-of-concept, and to reiterate the motivation for future self-consistent MHD simulations that include the embedded EMRI. revision: partial
-
Referee: §3 (Numerical integration): The torque time series is sampled independently at each timestep, implicitly assuming white, uncorrelated fluctuations; the manuscript does not evolve the disk and EMRI self-consistently, so possible correlations on orbital or viscous timescales are not tested and could cause the random walk in semi-major axis (and thus Δψ_gas) to saturate rather than accumulate to detectable levels.
Authors: The independent sampling of the torque at each timestep is indeed an assumption that corresponds to white-noise fluctuations. We recognize that real disk turbulence is expected to exhibit correlations on orbital and viscous timescales, which could in principle cause the random walk to saturate. Because the present work does not perform self-consistent disk+EMRI evolution, we cannot quantify the impact of such correlations. We have added a dedicated paragraph in the revised §3 and in the conclusions that highlights this assumption and its possible consequences, while emphasizing that the study is designed to motivate precisely the long-term, self-consistent simulations needed to test these effects. revision: partial
-
Referee: §4 (Results and detectability): The statement that Δψ_gas > 8/SNR holds for f_Edd ≳ 0.3, C ≳ 300, etc., is obtained after scanning the four-dimensional parameter space; the threshold is therefore applied post hoc rather than as a blind, a-priori prediction, making the central claim dependent on the chosen values of the free parameters.
Authors: We accept that the reported detectability thresholds emerge from a systematic scan of the (f_Edd, C, h0, α) parameter space rather than from a blind, a-priori forecast. The purpose of the survey is exploratory: to identify the regions of parameter space in which turbulent torques could push gas-induced dephasing above the LISA detection threshold. In the revised manuscript we have rephrased the abstract, §4, and the conclusions to present the results explicitly as an exploration of the viable parameter regime, to stress the dependence on the free parameters, and to avoid any implication of a definitive prediction. revision: yes
Circularity Check
No significant circularity in derivation chain
full rationale
The paper introduces an agnostic Gaussian model for T_turb centered on T_lin with free normalization C drawn from an external HD study, then performs a forward parameter sweep over f_Edd, C, h0, and α to compute Δψ_gas for a fixed EMRI. This is a model-based exploration of detectability thresholds rather than a closed derivation; the output dephasing scales with the input variance by design of the random-walk integration, but no step equates a claimed prediction to its own fitted inputs or reduces via self-citation to an unverified premise. The central result is conditional on the assumed model parameters and does not invoke uniqueness theorems or rename prior results. The derivation remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
free parameters (4)
- C (turbulence normalization)
- f_Edd (Eddington ratio)
- h0 (disk aspect ratio)
- α (turbo-viscous coefficient)
axioms (2)
- domain assumption The disk is turbulent due to magneto-rotational instability and the torque fluctuations can be modeled as a stationary Gaussian process around the mean linear torque.
- domain assumption The EMRI remains on a circular orbit while the torque is applied over the final four years.
invented entities (1)
-
Gaussian turbulent torque T_turb
no independent evidence
Reference graph
Works this paper leans on
-
[1]
2023, Living Reviews in Relativity, 26, 2, doi: 10.1007/s41114-022-00041-y
Amaro-Seoane, P., Andrews, J., Arca Sedda, M., et al. 2023, Living Reviews in Relativity, 26, 2, doi: 10.1007/s41114-022-00041-y Arca Sedda, M., Amaro Seoane, P., & Chen, X. 2021, A&A, 652, A54, doi: 10.1051/0004-6361/202037785
-
[2]
doi:10.1103/PhysRevD.89.104059, https://arxiv.org/abs/1404.7149 arXiv:1404.7149 [gr-qc]
Barausse, E., Cardoso, V., & Pani, P. 2014, PhRvD, 89, 104059, doi: 10.1103/PhysRevD.89.104059
-
[3]
Bardeen, J. M., & Petterson, J. A. 1975, ApJL, 195, L65, doi: 10.1086/181711
-
[4]
Baruteau, C., & Lin, D. N. C. 2010, ApJ, 709, 759, doi: 10.1088/0004-637X/709/2/759
-
[5]
M., Mac Low, M.-M., McKernan, B., & Ford, K
Bellovary, J. M., Mac Low, M.-M., McKernan, B., & Ford, K. E. S. 2016, ApJL, 819, L17, doi: 10.3847/2041-8205/819/2/L17
-
[6]
Sathyaprakash, B. S. 2009, PhRvD, 80, 084043, doi: 10.1103/PhysRevD.80.084043 Cevallos Soto, A., & Zhu, Z. 2026, ApJ, 998, 209, doi: 10.3847/1538-4357/ae355c
-
[7]
2025, MNRAS, 540, 1998, doi: 10.1093/mnras/staf867
Chen, Y.-X., Wu, Y., Li, Y.-P., et al. 2025, MNRAS, 540, 1998, doi: 10.1093/mnras/staf867
-
[8]
Colpi, M., Danzmann, K., Hewitson, M., et al. 2024, arXiv e-prints, arXiv:2402.07571, doi: 10.48550/arXiv.2402.07571 9
work page internal anchor Pith review doi:10.48550/arxiv.2402.07571 2024
-
[9]
The Astrophysical Journal891(2), 154 (2020)
Deng, H., Mayer, L., & Latter, H. 2020, ApJ, 891, 154, doi: 10.3847/1538-4357/ab77b2
-
[10]
2021, MNRAS, 501, 3540, doi: 10.1093/mnras/staa3976
MacFadyen, A. 2021, MNRAS, 501, 3540, doi: 10.1093/mnras/staa3976
-
[11]
Derdzinski, A., Dittmann, A. J., Franchini, A., et al. 2025, arXiv e-prints, arXiv:2512.10893, doi: 10.48550/arXiv.2512.10893
-
[12]
2025a, PhRvD, 111, 084006, doi: 10.1103/PhysRevD.111.084006
Duque, F., Kejriwal, S., Sberna, L., Speri, L., & Gair, J. 2025a, PhRvD, 111, 084006, doi: 10.1103/PhysRevD.111.084006
-
[13]
Duque, F., Sberna, L., Spiers, A., & Vicente, R. 2025b, arXiv e-prints, arXiv:2510.02433, doi: 10.48550/arXiv.2510.02433
-
[14]
2024a, MNRAS, 532, 4060, doi: 10.1093/mnras/stae1764
Garg, M., Derdzinski, A., Tiwari, S., Gair, J., & Mayer, L. 2024a, MNRAS, 532, 4060, doi: 10.1093/mnras/stae1764
-
[15]
2022, MNRAS, 517, 1339, doi: 10.1093/mnras/stac2711
Mayer, L. 2022, MNRAS, 517, 1339, doi: 10.1093/mnras/stac2711
-
[16]
2025, ApJ, 993, 145, doi: 10.3847/1538-4357/ae10b4
Garg, M., Franchini, A., Lupi, A., Bonetti, M., & Mayer, L. 2025, ApJ, 993, 145, doi: 10.3847/1538-4357/ae10b4
-
[17]
2024b, MNRAS, 535, 3283, doi: 10.1093/mnras/stae2605
Garg, M., Sberna, L., Speri, L., Duque, F., & Gair, J. 2024b, MNRAS, 535, 3283, doi: 10.1093/mnras/stae2605
-
[18]
Garg, M., Tiede, C., & D’Orazio, D. J. 2024c, MNRAS, 534, 3705, doi: 10.1093/mnras/stae2357
-
[19]
2024d, MNRAS, 528, 4176, doi: 10.1093/mnras/stad3477
Garg, M., Tiwari, S., Derdzinski, A., et al. 2024d, MNRAS, 528, 4176, doi: 10.1093/mnras/stad3477
-
[20]
Gong, X., Xu, S., Gui, S., Huang, S., & Lau, Y.-K. 2021, in Handbook of Gravitational Wave Astronomy (Springer Singapore), 24, doi: 10.1007/978-981-15-4702-7 24-1
-
[21]
Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357, doi: 10.1038/s41586-020-2649-2
-
[22]
2012, ApJ, 750, 34, doi: 10.1088/0004-637X/750/1/34 Hossein Nouri, F., & Janiuk, A
Horn, B., Lyra, W., Mac Low, M.-M., & S´ andor, Z. 2012, ApJ, 750, 34, doi: 10.1088/0004-637X/750/1/34 Hossein Nouri, F., & Janiuk, A. 2024, A&A, 687, A184, doi: 10.1051/0004-6361/202348796
-
[23]
Hunter, J. D. 2007, Computing in Science & Engineering, 9, 90, doi: 10.1109/MCSE.2007.55
-
[24]
2011, PhRvD, 84, 024032, doi: 10.1103/PhysRevD.84.024032
Kocsis, B., Yunes, N., & Loeb, A. 2011, PhRvD, 84, 024032, doi: 10.1103/PhysRevD.84.024032
-
[25]
Kormendy, J., & Ho, L. C. 2013, ARA&A, 51, 511, doi: 10.1146/annurev-astro-082708-101811
work page internal anchor Pith review doi:10.1146/annurev-astro-082708-101811 2013
-
[26]
Kubli, N., Mayer, L., Deng, H., & Lin, D. N. C. 2026, MNRAS, 546, stag024, doi: 10.1093/mnras/stag024
-
[27]
2025a, Reports on Progress in Physics, 88, 056901, doi: 10.1088/1361-6633/adc9be
Li, E.-K., Liu, S., Torres-Orjuela, A., et al. 2025, Reports on Progress in Physics, 88, 056901, doi: 10.1088/1361-6633/adc9be
-
[28]
McKernan, B., Ford, K. E. S., Cook, H. E., et al. 2025, ApJ, 990, 217, doi: 10.3847/1538-4357/adf114
-
[29]
McPike, E., Perna, R., Ford, K. E. S., et al. 2026, arXiv e-prints, arXiv:2602.04135, doi: 10.48550/arXiv.2602.04135
-
[30]
Miniutti, G., Saxton, R. D., Giustini, M., et al. 2019, Nature, 573, 381, doi: 10.1038/s41586-019-1556-x
-
[31]
Murphy, G. C., & Pessah, M. E. 2015, ApJ, 802, 139, doi: 10.1088/0004-637X/802/2/139
-
[32]
Nelson, R. P. 2005, A&A, 443, 1067, doi: 10.1051/0004-6361:20042605
-
[33]
Paardekooper, S.-J., Baruteau, C., Crida, A., & Kley, W. 2010, MNRAS, 401, 1950, doi: 10.1111/j.1365-2966.2009.15782.x pandas development team, T. 2020, pandas-dev/pandas:
-
[34]
Pandas, Zenodo, doi: 10.5281/zenodo.3509134
-
[35]
Peters, P. C. 1964, Physical Review, 136, 1224, doi: 10.1103/PhysRev.136.B1224
-
[36]
Monthly Notices of the Royal Astronomical Society474(2), 2212–2232 (2018)
Riols, A., & Latter, H. 2018, MNRAS, 474, 2212, doi: 10.1093/mnras/stx2455 —. 2019, MNRAS, 482, 3989, doi: 10.1093/mnras/sty2804
-
[37]
2020, ApJ, 903, 133, doi: 10.3847/1538-4357/abbc1d
Secunda, A., Bellovary, J., Mac Low, M.-M., et al. 2020, ApJ, 903, 133, doi: 10.3847/1538-4357/abbc1d
-
[38]
I., & Sunyaev, R
Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337
1973
-
[39]
2023, Physical Review X, 13, 021035, doi: 10.1103/PhysRevX.13.021035
Speri, L., Antonelli, A., Sberna, L., et al. 2023, Physical Review X, 13, 021035, doi: 10.1103/PhysRevX.13.021035
-
[40]
Tanaka, H., Takeuchi, T., & Ward, W. R. 2002, ApJ, 565, 1257, doi: 10.1086/324713
-
[41]
Trani, A. A., & Di Cintio, P. 2025, A&A, 703, A6, doi: 10.1051/0004-6361/202555879
-
[42]
Begelman, M. C. 2025, arXiv e-prints, arXiv:2510.21187, doi: 10.48550/arXiv.2510.21187
-
[43]
Waskom, M. L. 2021, Journal of Open Source Software, 6, 3021, doi: 10.21105/joss.03021
-
[44]
Wu, Y., Chen, Y.-X., & Lin, D. N. C. 2024, MNRAS, 528, L127, doi: 10.1093/mnrasl/slad183
-
[45]
2011, PhRvL, 107, 171103, doi: 10.1103/PhysRevLett.107.171103
Yunes, N., Kocsis, B., Loeb, A., & Haiman, Z. 2011, PhRvL, 107, 171103, doi: 10.1103/PhysRevLett.107.171103
-
[46]
2022, MNRAS, 511, 6143, doi: 10.1093/mnras/stac299
Mayer, L. 2022, MNRAS, 511, 6143, doi: 10.1093/mnras/stac299
-
[47]
Zwick, L., Tiede, C., Trani, A. A., et al. 2024, PhRvD, 110, 103005, doi: 10.1103/PhysRevD.110.103005 10 APPENDIX A.ANALYTICAL APPROXIMATION OF THE TURBULENT DEPHASING GAUSSIAN DISTRIBUTION As given in Eq. (12), the overall turbulent dephasing is ∆ψturb = ∆0 KX j=0 ξj a7 j −a 7 j+1 .(A1) The mean of ∆ψ turb over realizations of the{ξ i}is ⟨∆ψtrub⟩ξi = ∆0 ...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.