{"id":"2e9f9b08-93c5-44a9-9b14-b58ba3bef54b","arxiv_id":"2502.07929","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"The dephasing of gravitational waves from an EMRI encodes the equation of state of an accreting dark fluid, mainly through the fluid's global gravitational pull on the orbiting black hole.","lead":"An EMRI black hole binary surrounded by a steadily accreting dark fluid accumulates a gravitational-wave phase shift that depends on the fluid's equation of state. The effect is negligible at cosmic densities but grows sharply with density, offering a possible route to probing dark fluid properties with future space-based detectors.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Omitted 2PN and spin phase errors near 70 r_ISCO are ~0.1 rad over 10 yr, far above the claimed 1e-6 rad fluid dephasing, and the paper gives no template-error control.","rationale":"The reader's conditional verdict is appropriate. I do not see a fatal internal contradiction in the central dephasing mechanism: the test-fluid-to-global-force switch can be interpreted as a linearized perturbation calculation, and the reader's weakest assumption on self-gravity is less decisive than the omitted-PN systematic. The apparent factor or sign issue in Eq. (75) affects mainly the sudden-singularity extension and does not by itself overturn the Section 5.1 w-dephasing claim. The most load-bearing problem for the central claim is that the quoted 1e-6 rad dephasing is many orders of magnitude smaller than the phase error expected from terms the paper explicitly drops, namely 2PN conservative corrections and spin effects. Because the paper compares two equally truncated trajectories, the leading 2PN term cancels in the formal difference, but the stated scientific conclusion requires comparing against a phase model of the vacuum binary that is accurate to at least the signal level. No evidence is given that 2PN, spin, or parameter-estimation degeneracies are controlled at that level. This does not invalidate the calculation as a preliminary estimate, but it does prevent the stronger interpretation that GW phase tracking can in practice constrain the dark-fluid EoS. A concrete numerical experiment with 2PN and spin terms added would settle whether this concern is real or whether the cancellation is sufficiently clean.","tokens_in":30209,"tokens_out":31917,"duration_ms":321133,"concrete_test":"Add the standard 2PN two-body conservative acceleration A_2PN (from the references in Eqs. (2)-(4)) to Equations (13)-(14), first with zero spins, and rerun the 10-year r_init = 70 r_ISCO fluid and no-fluid cases of Fig. 5. If the recomputed Delta-phi differs from the 1PN+2.5PN result by more than ~1e-6 rad, the EoS-dependent dephasing is not robust below the omitted-PN systematic. Then add leading spin-orbit (1.5PN) and spin-spin (2PN) terms for representative aligned spins chi_1 = chi_2 = 0.9 and recompute Delta-phi; this quantifies whether the 1e-6 rad signal can be separated from spin degeneracy. A second, independent check is to compare the same Delta-phi against a full GR/numerical-relativity-informed EMRI waveform to confirm that omitted conservative and dissipative higher-PN terms do not shift the phase difference by more than the quoted signal.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2 retains only 1PN and 2.5PN and explicitly neglects 2PN and all spin terms. At r_init = 70 r_ISCO, v^2/c^2 = GM/(c^2 r) = 1/420 ~ 2.4e-3, so an order-of-magnitude estimate of the dropped 2PN phase over the ~2360 orbits in 10 yr is (v/c)^4 * 2pi N_orbit ~ 5.7e-6 * 1.5e4 ~ 0.08 rad. This is roughly 1e5 times larger than the claimed Delta-phi ~ 1e-6 rad in Fig. 5. The authors compute Delta-phi as a difference between fluid and no-fluid trajectories that both omit 2PN, so the leading 2PN contribution cancels in that idealized comparison. However, the central conclusion that gravitational-wave phase tracking can constrain the fluid EoS presupposes that the vacuum phase model is accurate below the fluid signal; otherwise the 1e-6 rad w-dependence is degenerate with omitted 2PN, spin-orbit, spin-spin, and parameter-estimation systematics. The manuscript provides no accuracy estimate, no spin-addition test, and no Fisher or systematic-error computation. The problem worsens for the 10 r_ISCO examples in the lower panel of Fig. 5, where v^2/c^2 ~ 0.033 and the omitted 2PN phase over the inspiral is of order 0.1 rad. Thus the physical mechanism may be correct, but the claimed ability to infer w from phase tracking is not yet supported without controlling these larger omitted phase terms.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies an extreme-mass-ratio inspiral (a 10^6 M_sun Schwarzschild black hole with a 10 M_sun companion) embedded in a spherically symmetric, steadily accreting dark fluid with a linear equation of state p = alpha (rho - rho0) c^2. The authors combine 1PN and 2.5PN equations of motion with local dynamical friction from the fluid and a global gravitational force from the fluid enclosed inside the perturber's orbit. They compute gravitational-wave dephasing relative to a no-fluid binary, report a dependence of the dephasing on the equation-of-state parameter w (of order 10^-6 rad over 10 years at rho_inf = rho_crit and r_init = 70 r_ISCO), and extend the formalism to sudden cosmological singularities, finding those give negligible dephasing at cosmological densities. The paper validates its numerical setup against the dark-matter-spike results of Montalvo et al. (2024).","tokens_in":30634,"tokens_out":22413,"duration_ms":200385,"significance":"If the reported dephasing numbers were correct, the paper would provide a useful framework for studying environmental effects on EMRI waveforms and a concrete, falsifiable prediction connecting a dark fluid's equation of state to gravitational-wave phase evolution. The manuscript has clear strengths: the equations of motion are explicit, the accretion profiles are derived from the standard Babichev et al. formalism, and the validation against an earlier dark-matter-spike calculation is a good check of the numerical pipeline. However, the central quantitative claim at cosmologically relevant densities is called into question by a simple order-of-magnitude estimate given in the major comments below, and the paper does not provide the systematic-error analysis needed to support the inference statement in Section 5.1.","major_comments":[{"comment":"The magnitude of the reported dephasing is inconsistent with the input physics by many orders of magnitude. For rho_inf = rho_crit, m1 = 10^6 M_sun, and r = 70 r_ISCO = 6.2e11 m, the enclosed active mass is |M_enc| ~ (4 pi / 3) |rho + 3p/c^2| r^3 ~ 1.7e10 kg (for w = -1), so the fluid acceleration on the perturber is G |M_enc| / r^2 ~ 3e-24 m/s^2 and the fractional change of the central acceleration is |M_enc|/m1 ~ 8e-27. With about 1.2e3 orbits in 10 years, the accumulated dephasing is Delta-phi ~ pi N |M_enc|/m1 ~ 1e-23 rad, not the plotted 10^-6 rad. This 17-order-of-magnitude discrepancy suggests a unit conversion or rescaling error in the numerical implementation of Eq. (60) (for example, using the SI density without the G/c^2 conversion to active gravitational mass). The authors should verify their code against the enclosed-mass benchmark and report the physical acceleration and dephasing for a test case.","section":"5.1, Eq. (60), Fig. 5 upper panel"},{"comment":"The inference claim in Section 5.1 ('the equation of state of the accreting energy density can be inferred from its impact on the binary system's dynamics') is not supported by a systematic-error budget. The model omits 2PN and all spin-orbit and spin-spin terms; at r_init = 70 r_ISCO, (v/c)^2 ~ 2.4e-3, so the omitted 2PN phase over 10 years is of order (v/c)^4 * 2 pi N ~ 0.08 rad, roughly five orders of magnitude above the claimed fluid dephasing. Computing Delta-phi as a fluid-minus-no-fluid difference cancels the leading omitted PN terms in that idealized comparison, but a real measurement requires a vacuum waveform model accurate below the fluid signal, and the paper provides no Fisher-matrix, mismatch, or systematic-error estimate. The Conclusions acknowledge that such effects may be overshadowed, but the stronger statement in Section 5.1 should either be removed or backed by a concrete template-accuracy calculation.","section":"2, 5.1, and Conclusions"},{"comment":"The treatment of the fluid's self-gravity is asymmetric and its domain of validity is not quantified. The accretion profiles are derived under the explicit assumption that 'the fluid's energy density is sufficiently low, such that its self-gravity can be neglected' (Section 4), yet the same fluid is used to produce a global gravitational force on the perturber of the form G M_enc(r)/r^2 in Eq. (60). This is consistent only to leading order in M_enc(r)/m1, and the manuscript does not state this small parameter or give its value for the plotted examples. The extrapolation to 'higher energy densities' where the effect becomes significant is precisely the regime where M_enc(r)/m1 may no longer be small (for example, at 10^25 rho_crit and 70 r_ISCO, M_enc/m1 ~ 4e-2). The paper should specify the validity condition and report M_enc/m1 for the cases shown.","section":"4 and 5, Eq. (60)"}],"minor_comments":[{"comment":"The sentence 'As long as m1 << m2, such that the center of mass nearly coincides with the center of the supermassive black hole' has the mass ratio reversed; it should read m2 << m1.","section":"5"},{"comment":"The quadrupole moment in Eq. (62) includes the static, spherically symmetric T00_fluid term. Such a term contributes a constant to Mij and therefore has zero second time derivative, so it does not affect the waveform polarizations in Eqs. (18)-(19); the text should clarify why it is included or remove it.","section":"5, Eq. (62)"},{"comment":"The notation in Eq. (11) appears garbled ('tyre', 'ctyre2'); the definitions of ar r and ar t should be typeset as ar r = r/(c t_yr e) and ar t = t/t_yr, or similar, and the text should define all symbols.","section":"2, Eq. (11)"},{"comment":"The caption refers to 'Left Panel' and 'Right Panel', while the text refers to 'upper panel' and 'lower panel'; the panel labeling should be made consistent.","section":"Figure 7 caption"}],"recommendation":"reject","confidential_remarks":"The order-of-magnitude estimate in major comment 1 suggests the central numerical result of the paper is not physically correct as presented; if it is a unit conversion error, correcting it would reduce the predicted dephasing at rho_crit by roughly 17 orders of magnitude and eliminate the paper's main quantitative motivation. The paper also lacks the systematic-error analysis needed to support the inference claim. I would encourage the authors to re-examine the numerical implementation and, if a corrected version still yields a meaningful effect, to resubmit with explicit benchmarks and a template-accuracy estimate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The new thing here is the unified treatment of steady-state fluid accretion around a Schwarzschild black hole plus the enclosed-fluid global force on the orbiting perturber, extending Montalvo et al.'s static-spike analysis. The authors validate against Montalvo's waveforms, which is good, and they compute the EoS-dependent dephasing from first principles rather than fitting it. The sudden-singularity extension is a minor but legitimate bonus, and the paper is refreshingly explicit that this is a preliminary framework, not a detection claim.\n\nI read the central calculation as internally coherent. The dephasing is a difference between fluid and no-fluid trajectories, so the omitted 2PN and spin terms largely cancel in that idealized comparison. The reader's stress-test estimate of ~0.1 rad from 2PN is relevant to the paper's aspirational claim that 'the equation of state of the accreting energy density can be inferred,' but it does not invalidate the theoretical result that the fluid induces a w-dependent phase shift. What is missing is the error budget that would connect the two: no Fisher analysis, no spin-addition test, no systematic-error estimate. At cosmological densities the predicted shift is ~1e-6 rad over 10 years, which is many orders below LISA sensitivity and also below the omitted conservative terms unless those are controlled. The paper acknowledges this in words but does not quantify it.\n\nThe test-fluid-to-global-force switch—neglecting fluid self-gravity in the accretion profile but then using that same fluid as a source of the Newtonian force in Eq. (60)—is a standard approximation but deserves a sentence justifying consistency. Also, the sign of the integrand in Eq. (75) looked off to me; that is worth a careful check but is not load-bearing for the main result.\n\nWho should read this: people working on environmental effects in EMRI waveforms and dark-matter/fluid distributions around black holes. They will find the derivation clear and the validation useful. The detectability claim needs serious follow-up work before it can be advertised.\n\nIt deserves peer review, not desk rejection. A good referee will ask for a systematic-error estimate and a quantitative comparison between the fluid-induced phase and the omitted 2PN/spin phase for the same orbital setup. The paper is preliminary in the good sense: it sets up a framework and shows the physics, without pretending to have the final answer.","headline":"A coherent framework for accreting-fluid effects on EMRI dephasing, honestly labeled preliminary, but the detectability claim rests on systematics control the paper does not provide.","tokens_in":31159,"tokens_out":1614,"would_cite":true,"duration_ms":18336,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C35","83C57","83F05"],"pacs":["04.30.-w","04.70.Bw","95.35.+d","98.80.-k"],"model":"deepseek-v4-flash","headline":"An accreting dark fluid around a supermassive black hole imprints its equation of state on the gravitational-wave phase of an inspiraling stellar-mass companion.","keywords":["gravitational waves","dark fluid accretion","extreme mass ratio inspirals","post-Newtonian approximation","spherical accretion","equation of state","de-phasing","sudden cosmological singularities"],"falsifier":"Run a fully relativistic hydrodynamical simulation of a test fluid accreting onto a Schwarzschild black hole with a companion moving through it, including the fluid's self-gravity, and compare the resulting radial force on the companion with Eq. (60) at radii of $10$-$250\\,r_{\\rm ISCO}$; if the sign or magnitude of the force changes, the claimed $w$-dependence of the de-phasing would not survive.","tokens_in":29950,"feed_emoji":"🕳️","tokens_out":7100,"duration_ms":60413,"temperature":0.7,"pith_summary":"This paper argues that a dark fluid steadily accreting onto a supermassive black hole can be read off from the gravitational-wave phase of an extreme-mass-ratio inspiral, provided the ambient density is high enough. The central prediction is that the accumulated de-phasing between the fluid-dressed and the vacuum binary is controlled by the fluid's equation-of-state parameter $w$, with the sign of the de-phasing flipping near $w\\simeq 0$ at an initial separation of $70\\,r_{\\rm ISCO}$. At cosmologically critical densities the ten-year de-phasing is tiny, about $10^{-6}$ rad, but it grows with density and reaches $\\simeq -0.01$ rad in four years for $\\rho_\\infty=10^6\\rho_{\\rm crit}$, a range where future space-based detectors could see it. The paper also studies sudden cosmological singularities, finding that their effect on the waveform is negligible at cosmologically relevant densities even though they can make an initially circular orbit slightly eccentric.","feed_headline":"Dark fluid accretion flips the sign of a binary's wave de-phasing","feed_subtitle":"The predicted phase shift is tiny at cosmic densities but grows with density, and its sign tracks the fluid's equation of state.","key_machinery":"The machinery is a post-Newtonian Lagrangian two-body model in which the 1PN and 2.5PN corrections are inserted as generalized forces (Eqs. (13)-(14)), supplemented by two fluid effects. The accretion profile comes from the steady-state, spherically symmetric Michel solution for a perfect fluid with linear equation of state $p=\\alpha(\\rho-\\rho_0)c^2$, giving radial velocity and density profiles $u^r(r)$ and $\\rho(r)$ from Eqs. (55)-(56). The load-bearing object is the global radial force $f_r^{\\rm fluid}=-4\\pi G r^{-2}\\int_{r_S}^{r}\\left(\\rho(r')+3p(r')/c^2\\right) r'^2\\,dr'$ (Eq. (60)), the spherical-shell gravitational pull of the fluid on the perturber, which dominates over dynamical friction at large separations and whose zero controls the sign flip of the de-phasing. Waveforms are then computed from the quadrupole formula using the fluid-corrected trajectories, and phase shifts are read off as $\\Delta\\phi=\\phi_{\\rm fluid}-\\phi$.","core_discovery":"The paper claims that for an extreme-mass-ratio inspiral composed of a $10^6\\,M_\\odot$ Schwarzschild black hole and a $10\\,M_\\odot$ perturber, the dominant environmental effect of a spherically accreting dark fluid is not the local dynamical friction but the global gravitational pull of the fluid shell between the horizon and the perturber. The de-phasing $\\Delta\\phi\\equiv\\phi_{\\rm fluid}-\\phi$ at ten years is negative for $w\\lesssim 0$ and positive for $w\\gtrsim 0$ at $r_{\\rm init}=70\\,r_{\\rm ISCO}$ when $\\rho_\\infty=\\rho_{\\rm crit}$, passing through zero where the integrated force in Eq. (60) vanishes, and the precise crossing value shifts with distance. Quantitatively, for $\\rho_\\infty=\\rho_{\\rm crit}$ the ten-year shift is of order $10^{-6}$ rad, while at $\\rho_\\infty=10^6\\rho_{\\rm crit}$ and $r_{\\rm init}=10\\,r_{\\rm ISCO}$ the four-year de-phasing is $\\sim -0.01$ rad for $w\\simeq -1$. The paper further claims that sudden cosmological singularities produce only a negligible waveform de-phasing at cosmologically relevant densities, although they kick the perturber's radial velocity by an amount that grows with distance.","pith_inferences":["The paper leaves implicit that the sign flip near $w\\simeq 0$ could serve as a model-independent diagnostic: if an observed EMRI de-phasing tracks the enclosed-mass integral rather than the local density, that favors a pressure-supporting dark fluid over a collisionless spike, whose global and local effects scale differently.","At densities far above $\\rho_{\\rm crit}$, non-spherically symmetric accretion may make dynamical friction dominate again, as the paper itself notes, so the sign-flip diagnostic would need calibration on inflow geometry before being used as a precise equation-of-state probe.","A natural test of the framework would be a fully relativistic hydrodynamical simulation of Bondi-Michel accretion onto a Schwarzschild black hole with a perturbing companion, checking whether the integrated shell force matches Eq. (60) to the claimed precision at radii of $10$-$250\\,r_{\\rm ISCO}$.","The sudden-singularity velocity kick in Eq. (76) could be confronted with cosmological simulations of type-II singularities; if future observations ever tie an EMRI eccentricity to a singularity time, that would give an independent constraint on the jump parameter $\\eta$ and the background equation-of-state parameter $w_\\infty$."],"forward_implications":["At large separations of tens to hundreds of $r_{\\rm ISCO}$, de-phasing from spherical accretion is dominated by the fluid's global gravity rather than dynamical friction, so measurements of phase can constrain the enclosed fluid mass instead of only local dissipation.","The sign of $\\Delta\\phi$ at a given distance is a function of $w$: stiff fluids with $w\\gtrsim 0$ produce a positive shift while dark-energy-like fluids with $w\\lesssim 0$ produce a negative shift, with the transition point depending on orbital radius.","For cosmologically relevant densities the predicted shift is too small to observe, but for densities of $10^6\\rho_{\\rm crit}$ the four-year shift reaches $\\sim -0.01$ rad at $r_{\\rm init}=10\\,r_{\\rm ISCO}$, potentially within reach of next-generation space-based gravitational-wave detectors.","Sudden cosmological singularities can deform an initially circular orbit into a mildly eccentric one through a velocity kick that grows with distance, but the resulting gravitational-wave de-phasing is negligible at $\\rho_\\infty=\\rho_{\\rm crit}$.","The Lagrangian generalized-force framework applies to any dissipative or conservative environmental effect expressible as a force, so the same method can be reused for other dark-fluid equations of state and for more general fluid geometries."],"supporting_citations":[{"why":"Supplies the Lagrangian generalized-force framework with 1PN and 2.5PN corrections and a static dark-matter spike, which the paper validates and extends.","marker":"[43]"},{"why":"Provides the test-fluid spherical accretion solution and the linear equation of state $p=\\alpha(\\rho-\\rho_0)c^2$ used to derive the density and velocity profiles.","marker":"[49]"},{"why":"Gives the steady-state accretion of dark energy onto a Schwarzschild black hole, including the integrals of motion and the flux constant $A$ used in Eqs. (55)-(57).","marker":"[53]"},{"why":"Introduces the steady-state spherically symmetric accretion solution and the critical-point (sonic-point) regularity condition that fixes the accretion constant.","marker":"[51]"},{"why":"Provides the relativistic dynamical-friction force on a perturber moving through a collisional fluid, giving the components $F_r$ and $F_\\phi$ in Eqs. (63)-(64).","marker":"[77]"},{"why":"Supplies the Mach-number-dependent fits $I_r$ and $I_\\phi$ for the gravitational wake drag used in the dynamical-friction terms.","marker":"[82]"},{"why":"Supplies the relativistic dark-matter spike profile and effective scaling parameters used in the validation example with a static spike.","marker":"[46]"},{"why":"Defines sudden cosmological singularities and their pressure divergence, which the paper adapts to model the singularity-induced velocity kick.","marker":"[66]"},{"why":"Provides the parametrized scale factor for a past sudden singularity, used to compute the time-dependent background density and pressure in Section 5.5.","marker":"[93]"}],"fun_headline_variants":["Global dark fluid pull dominates local friction in binaries","Tiny dark fluid effect flips sign with pressure","Gravitational wave de-phasing tracks dark fluid's equation of state","Dark fluid shifts binary orbit enough to see in waves"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that the fluid's own gravity is too weak to affect its own density profile, yet the fluid's gravity is strong enough to pull the orbiting black hole; if that combination is internally inconsistent, the predicted phase shifts would change.","fun_headline_variants_meta":{"raw":{"variants":["Global dark fluid pull dominates local friction in binaries","Tiny dark fluid effect flips sign with pressure","Gravitational wave de-phasing tracks dark fluid's equation of state","Dark fluid shifts binary orbit enough to see in waves"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000727,"raw_usage":{"total_tokens":3331,"prompt_tokens":1094,"completion_tokens":2237,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":710,"completion_tokens_details":{"reasoning_tokens":2171}},"tokens_in":710,"tokens_out":2237,"duration_ms":17520,"temperature":1.0,"reasoning_tokens":2171,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T11:25:34.958260+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a fully relativistic hydrodynamical simulation of a test fluid accreting onto a Schwarzschild black hole with a companion moving through it, including the fluid's self-gravity, and compare the resulting radial force on the companion with Eq. (60) at radii of $10$-$250\\,r_{\\rm ISCO}$; if the sign or magnitude of the force changes, the claimed $w$-dependence of the de-phasing would not survive.","supporting_citations":[{"cited_title":"Post-Newtonian Effects in Compact Binaries with a Dark Matter Spike: A Lagrangian Approach","cited_arxiv_id":null,"evidence_quote":"Supplies the Lagrangian generalized-force framework with 1PN and 2.5PN corrections and a static dark-matter spike, which the paper validates and extends."},{"cited_title":"Accretion of matter by condensed objects","cited_arxiv_id":null,"evidence_quote":"Introduces the steady-state spherically symmetric accretion solution and the critical-point (sonic-point) regularity condition that fixes the accretion constant."},{"cited_title":"Relativistic dynamical friction in a collisional fluid","cited_arxiv_id":null,"evidence_quote":"Provides the relativistic dynamical-friction force on a perturber moving through a collisional fluid, giving the components $F_r$ and $F_\\phi$ in Eqs. (63)-(64)."},{"cited_title":"Sudden future singularities","cited_arxiv_id":null,"evidence_quote":"Defines sudden cosmological singularities and their pressure divergence, which the paper adapts to model the singularity-induced velocity kick."},{"cited_title":"Effects of a Late Gravitational Transition on Gravitational Waves and Anticipated Constraints","cited_arxiv_id":null,"evidence_quote":"Provides the parametrized scale factor for a past sudden singularity, used to compute the time-dependent background density and pressure in Section 5.5."}],"review_version":1}