{"id":"0beaa1af-109e-4edb-832d-fb589105b9fc","arxiv_id":"2507.02045","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Scalar radiation from an EMRI in an ultralight scalar cloud is computed semi-analytically, showing dipole clouds decelerate and quadrupole clouds accelerate the inspiral, with up to about 100 rad dephasing after 18 months.","lead":"A small black hole spiraling into a supermassive black hole surrounded by an ultralight scalar cloud emits scalar radiation that shifts the gravitational-wave phase. This paper computes that radiation with a new semi-analytical method and shows the cloud's shape determines whether the inspiral speeds up or slows down.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (61) scales the scalar flux by the gravitational mode frequency m_g*omega_g instead of the scalar mode frequency omega; this could change flux magnitudes by a factor ~30 and potentially reverse the sign of the net dipole scalar flux.","rationale":"The paper is a careful and well-structured computation, with explicit truncation-error estimates and a cross-check against Dyson et al. However, the central quantitative results are the scalar energy fluxes and the derived dephasing, and those all pass through Eq. (61). The reader's identified weak point is the external Lorenz-gauge reconstruction code h(0,1); that is a reasonable concern, but the flux normalization in Eq. (61) is more directly load-bearing because it affects every reported flux value and the sign of the net dipole flux. A simple derivation from the paper's own definitions (Eqs. 55a, 59a, 60a) gives Edot_s,inf proportional to omega + omega_c, whereas Eq. (61a) uses m_g*omega_g. For the dominant modes these differ by a factor of roughly 28, so the published flux magnitudes, the sign of the net dipole flux at r0 < 26.7M, and the O(10^2) rad dephasing claim could all change substantially. The proposed test is definitive and cheap: evaluate one mode directly from the stress tensor using the paper's already-computed radial function. If the test confirms Eq. (61), the paper's central claims stand and my concern is withdrawn; if it confirms the alternative normalization, the paper needs a corrected flux formula and revised conclusions. Therefore the appropriate verdict is CONDITIONAL pending this check, rather than an unconditional accept or reject.","tokens_in":33352,"tokens_out":34428,"duration_ms":394038,"concrete_test":"Take the numerical Green's-function solution R_{ell m omega}(r) that the paper already produces for the ell = m = 3, m_g = 2 mode of the dipole cloud at r0 = 20M, insert it into Eq. (55a) using the asymptotic form R ~ A exp(i k r)/r, and compute Edot_s,inf = Edot_Phi,inf + omega_c * Qdot_inf directly from the stress tensor and Noether current. Compare the result with Eq. (61a). If the ratio equals (omega + omega_c)/(m_g*omega_g) ≈ 28, Eq. (61a) is mis-normalized and the central flux signs and dephasing estimates are unreliable; if the ratio is 1, the formula is correct and the concern is resolved. This check is a one-mode analytic exercise using the paper's own output and can settle the issue without new numerical infrastructure.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Starting from the paper's own definitions, for a scalar radiation mode R(r)S(theta) exp(-i omega t + i m phi), the Noether charge flux at infinity is Qdot_inf = lim r^2 integral j^r dOmega = 2k|A|^2, and Eq. (55a) gives Edot_Phi,inf = omega * Qdot_inf. Hence the total scalar flux defined in Eq. (60a) is Edot_s,inf = (omega + omega_c) * Qdot_inf. Eq. (61a), however, uses the prefactor m_g*omega_g, where m_g and omega_g are the azimuthal number and frequency of the sourcing metric perturbation h(0,1), not of the scalar mode. For the fiducial dipole cloud at r0 = 20M (a = 0.88M), the dominant infinity mode has ell = m = 3, so m_g = 2, omega_c ≈ 0.296, omega = omega_c + 2 Omega ≈ 0.318, and m_g*omega_g = 2 Omega ≈ 0.022; the normalization differs by roughly a factor of 28. The same issue appears in the horizon flux Eq. (61b) and in Eq. (61c), which sets Ldot/Edot = 1/Omega_g, whereas a single scalar mode carries Ldot/Edot = m/omega. The central claim that the dipole cloud produces a negative net flux inside r0 ≈ 26.7M relies on the negative horizon flux overwhelming the positive infinity flux; an O(28) change in the relative normalization of these fluxes can flip the sign of the total and invalidate the claim that the dipole cloud slows the inspiral. The cross-check with [105] does not resolve this, because [105] appears to adopt the same m_g*omega_g normalization. The paper does not derive Eq. (61) from the stress-energy tensor, and the text merely cites [46,105].","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a semi-analytical, modified-Teukolsky-based method to compute the scalar radiation emitted by an extreme mass-ratio inspiral embedded in an ultralight, superradiant scalar cloud around a Kerr black hole. The authors construct the scalar source from the background cloud profile and a Lorenz-gauge reconstructed metric perturbation, solve the radial equation by a Green's-function method, and evaluate energy and angular-momentum fluxes at the horizon and infinity for dipolar and quadrupolar clouds. They compare with the recent calculation of Dyson et al., find broad agreement, and use the fluxes to estimate a gravitational-wave dephasing of O(10^2) rad after 18 months for a representative system. The main physical claims are that the dipolar cloud produces a net negative scalar energy flux inside r0≈26.7M, slowing the inspiral, while the quadrupolar cloud produces an always-positive net flux, accelerating the inspiral.","tokens_in":33691,"tokens_out":31396,"duration_ms":365159,"significance":"If the calculation is correct, this is a valuable step toward waveform-accurate modeling of environmental effects in EMRIs: it provides a fully relativistic treatment of scalar radiation from a rotating black hole, extends previous work to quadrupolar clouds, quantifies truncation errors in the source and mode sums, and gives a falsifiable dephasing prediction for LISA-like detectors. The semi-analytical source decomposition and the Green's-function solution are presented in enough detail to be adapted to the gravitational sector in the announced follow-up work. The comparison with the independent numerical approach of [105] is a genuine validation of the implementation, and the paper is honest about the residual discrepancies, although the validation is weakened by a small spin mismatch between the two calculations.","major_comments":[{"comment":"The central validation is a comparison to [105], but the text states that [105] used a=0.877M while the present work uses a=0.88M. The quoted '<5%' and '<6%' agreement therefore includes a parameter mismatch. Please recompute the comparison at identical spin values, or, if that is impractical, quantify the effect of Delta a = 0.003M on each reported difference and restate the agreement accordingly. This is especially relevant for the 11.4% discrepancy at ell=8, m=2, which the text partially attributes to the spin mismatch; as written, one cannot separate numerical error from the parameter difference.","section":"Secs. III and V, Figs. 2, 4, 5, 9"},{"comment":"Equation (62) has a sign error in the denominator for modes with m_g<0. For the quadrupolar cloud with m=1 (so m_g=-1), a=0.88M and muM=0.3 give m omega_+ - omega_c ≈ -1.5e-4 / M, so the right-hand side of Eq. (62) is negative or complex and the printed inequality is not meaningful. The threshold quoted in the text, r0≈370M, follows instead from the condition omega_c - m omega_+ < |m_g| Omega, i.e., r0 < (|m_g|/(omega_c - m omega_+) - a)^{2/3}. Please correct Eq. (62) and the sentence 'as one can directly find from Eq. (62)'.","section":"Sec. V, Eq. (62)"},{"comment":"The scalar source and hence every flux in the paper depend on the Lorenz-gauge metric perturbation h^{(0,1)}_{mu nu} taken directly from the Mathematica notebooks of [98]. The comparison with [105] uses the same reconstructed metric, so it validates the present implementation of the contraction and angular projection but not the underlying reconstruction itself. Please add an explicit limitation statement to Sec. VI noting that the results inherit the accuracy of the Kerr-Lorenz-Circ code, or, if feasible, perform an independent test of the reconstructed h^{(0,1)} (for example, against the vacuum gravitational energy flux for circular equatorial orbits).","section":"Secs. III and VI"}],"minor_comments":[{"comment":"The captions say the radiation is plotted 'at the equilateral plane (i.e., theta = 0)'; since theta = 0 is the polar axis, this should presumably read 'equatorial plane (theta = pi/2)'.","section":"Figs. 7 and 8 captions"},{"comment":"There is a typo in the Discussion: 'comptuation' should be 'computation'.","section":"Sec. VI"},{"comment":"The dephasing delta_phi shown in Fig. 12 is never defined by an equation; please give the explicit formula used to evolve the orbital radius and to compute delta_phi from the scalar fluxes.","section":"Sec. V, Fig. 12"},{"comment":"The Mathematica notebook implementing the source construction is said to be provided 'upon request'; for reproducibility, please provide a persistent repository link alongside the existing supplementary notebook reference.","section":"Sec. III, end of the source evaluation"},{"comment":"The statement that the angular-momentum flux is completely determined by the energy flux through Eq. (61c) relies on the specific definitions in Eq. (60) and the sign convention for the Noether charge flux; a brief remark explaining that this is not the single-mode identity Ldot/Edot = m/omega would prevent confusion.","section":"Sec. V, Eq. (61c)"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern about Eq. (61) does not, in my reading, land: m_g omega_g equals omega - omega_c, and the Noether charge flux Qdot is negative for outgoing radiation, so Eqs. (61a)-(61c) are consistent with Eqs. (55), (59), and (60). The load-bearing issues are instead the spin mismatch in the validation comparison and the sign error in Eq. (62), both of which are fixable. I would ask the authors to address those before acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real content here is the scalar radiation fluxes for an EMRI in a superradiant scalar cloud, computed semi-analytically and cross-checked against Dyson et al. The dipolar cloud part is largely a replication, but the quadrupolar cloud analysis is new and the dephasing estimate gives the paper observational bite. The calculation looks sound to me: truncation errors are quantified, the source-projection method is a genuine technical improvement, and the agreement with [105] at the few-percent level is a credible independent check.\n\nOne thing I checked carefully: the stress-test claim about Eq. (61) does not hold up. It drops the minus sign in the paper’s definition of Qdot, Eq. (59). With that sign included, E_Phi,inf = -omega Qdot_inf, so E_s,inf = (omega_c - omega) Qdot_inf = -m_g Omega_g Qdot_inf. That is exactly the m_g omega_g prefactor in Eq. (61a), and similarly for the horizon. The suspected factor-of-28 error is a sign-convention artifact, not a bug. The paper could still do a better job deriving Eq. (61) instead of citing [46,105], but the formula is consistent with its own definitions.\n\nThe soft spots are real but manageable. The metric perturbation h(0,1) is taken from the Kerr-Lorenz-Circ notebooks without independent verification; if that code is wrong, the scalar source is wrong. The reproducibility story is also weaker than it should be: one line says the notebook is available on request, another gives a GitHub link with no commit hash. The claimed corrections to [105] are asserted but not substantiated. And the MTF framing overstates what is actually computed — the modified Teukolsky equation is written down but not solved, so this is really a scalar-flux paper with an MTF roadmap attached. The dephasing estimate is admittedly adiabatic and uses only the dominant gravitational quadrupole mode, which is fine for an estimate but not a waveform.\n\nThe 11% discrepancy in one subdominant infinity mode is minor; it is 10^12 smaller than the dominant mode, and the authors say so. The sign change of the horizon flux for the quadrupole cloud and the resulting acceleration of the inspiral is the most interesting new result, and I don't see anything that undermines it.\n\nWho is this for? Anyone working on environmental effects for LISA, especially EMRI dephasing in boson clouds. It deserves a serious referee. I would send it to review with requests to pin down the code availability and to add a short derivation or clear citation for Eq. (61), but I would not block on either.","headline":"Solid new flux calculations for EMRIs in scalar clouds; the suspected normalization bug in Eq. (61) is actually a sign-convention artifact — worth a careful referee, not a desk reject.","tokens_in":34306,"tokens_out":23451,"would_cite":true,"duration_ms":271207,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C57","83C35","83C25","83C10"],"pacs":["04.30.-w","04.70.-s","95.30.Sf"],"model":"deepseek-v4-flash","headline":"An EMRI inside a superradiant scalar cloud emits scalar radiation whose backreaction changes the inspiral rate in a cloud-mode-dependent way, producing gravitational-wave dephasing of order hundreds of radians.","keywords":["extreme mass-ratio inspiral","ultralight scalar cloud","superradiance","scalar radiation","modified Teukolsky formalism","energy flux","gravitational-wave dephasing","Kerr black hole"],"falsifier":"An independent, grid-based numerical integration of the scalar source term (without the Fourier-series separation) for the same parameters ($a=0.88M$, $\\mu M=0.3$, secondary at $r_0=10M$) should reproduce both the sign and magnitude of the net scalar energy flux: negative for the dipole cloud and positive for the quadrupole cloud. A null or opposite-sign net flux would refute the central claim.","tokens_in":33078,"feed_emoji":"🛰️","tokens_out":9927,"duration_ms":84989,"temperature":0.7,"pith_summary":"This paper studies an extreme-mass-ratio inspiral—a stellar-mass compact object spiraling into a supermassive Kerr black hole—when the black hole is dressed in an ultralight scalar cloud formed by superradiance. It computes, in full general relativity, the scalar radiation excited in the cloud by the secondary and the backreaction of that radiation on the orbit. The central result is that the sign of the net scalar energy flux depends on which cloud mode is populated: a dipole ($\\ell_c=m_c=1$) cloud gives a net negative flux inside $r_0\\approx26.7M$, slowing the inspiral, while a quadrupole ($\\ell_c=m_c=2$) cloud gives a positive flux at all radii, accelerating it. For a cloud mass $M_c=10^{-4}M$ and mass ratio $\\epsilon=10^{-5}$, the accumulated gravitational-wave phase shift reaches $O(10^2)$ radians after 18 months, well above the one-radian phase accuracy space-based detectors aim for, so the effect could be observable in future EMRI waveforms.","feed_headline":"A scalar cloud can brake or speed an EMRI","feed_subtitle":"Dipole clouds brake the inspiral inside 26.7M; quadrupole clouds speed it up—a dephasing of ~100 radians for LISA.","key_machinery":"The machinery is a two-parameter expansion in the mass ratio $\\epsilon$ and cloud amplitude $\\zeta$, organized with the modified Teukolsky formalism (a decoupling of the Newman-Penrose equations on non-vacuum black-hole backgrounds). The scalar sector reduces to the sourced Klein-Gordon equation $(\\Box-\\mu^2)\\Phi^{(1,1)}=S_\\Phi^{(1,1)}$, whose source is the contraction of the background cloud profile $\\Phi^{(1,0)}$ with the Lorenz-gauge reconstructed metric perturbation $h_{\\mu\\nu}^{(0,1)}$ of the secondary. The paper renders this source explicitly separable by projecting onto the Kinnersley tetrad, replacing directional derivatives with Chandrasekhar operators, and expanding the factors $\\Gamma(r,\\theta)^{-\\beta}\\bar\\Gamma(r,\\theta)^{-\\sigma}$ in Fourier series in $\\cos\\theta$ following [141]; the radial equation is then solved by a Green's function built from in/up homogeneous solutions. The selection rule $m=m_c+m_g$, $\\omega=\\omega_c+m_g\\Omega_g$ couples cloud modes to the secondary's orbital harmonics, and the transition of a radiation mode from oscillatory to exponentially decaying behavior near infinity (when $|\\omega|<\\mu$) produces the sharp flux feature at $r_0\\approx41.66M$ for the dipole cloud.","core_discovery":"The paper establishes that the scalar radiation field $\\Phi^{(1,1)}$ sourced by an EMRI in a quasi-bound scalar cloud around a Kerr black hole carries energy fluxes to infinity and into the horizon whose balance depends on the cloud's quantum numbers. For the fundamental dipolar cloud ($\\ell_c=m_c=1$, $n_c=0$) the horizon flux dominates below $r_0\\approx26.7M$ and is negative, so the total scalar flux removes orbital energy and slows the inspiral. For the quadrupolar cloud ($\\ell_c=m_c=2$, $n_c=0$) the total flux is positive at every radius, accelerating the inspiral, and the horizon flux itself becomes positive inside $r_0\\approx18.1M$, depositing energy into the black hole. Mode-by-mode energy fluxes agree with the independent calculation of [105] to within a few percent, with relative differences under $6\\%$ for the total horizon and infinity fluxes, and the resulting dephasing is $O(10^2)$ radians after 18 months for $M_c=10^{-4}M$ and $\\epsilon=10^{-5}$.","pith_inferences":["Because the sign of the net scalar flux is tied to the cloud's $(\\ell_c,m_c)$, a measured dephasing sign in a real EMRI would identify which superradiant mode dominates, effectively turning gravitational waves into a cloud-resolving probe.","The radius where the $m_g=1$ radiation mode crosses from oscillatory to evanescent behavior should appear as a sudden drop in the scalar energy flux, so a chirp 'kink' at the corresponding frequency could constrain the scalar mass even when the transition lies outside LISA's nominal band.","The same source-extraction machinery (Chandrasekhar operators plus Fourier decomposition of $\\Gamma$ factors) transfers directly to the modified Teukolsky source terms for gravitational radiation, so the follow-up gravitational calculation should inherit the same selection-rule structure and show corresponding mode-dependent flux signs."],"forward_implications":["Waveform templates for EMRIs in scalar clouds must include cloud-mode-dependent scalar backreaction: dipole clouds decelerate the secondary inside $r_0\\approx26.7M$, while quadrupole clouds accelerate it at all radii.","The dephasing can reach $O(10^2)$ radians after 18 months for $M_c=10^{-4}M$ and $\\epsilon=10^{-5}$, far exceeding the roughly one-radian phase-accuracy goal, making the environmental effect potentially measurable rather than a negligible correction.","Because the scalar-to-gravitational flux ratio is below $10^{-3}$ at $r_0\\leq10M$ for $M_c=10^{-4}M$, the effect is small but not negligible, and the few-percent differences between the two independent flux calculations still translate into $O(1)$ radian dephasing over a year.","The sharp decrease of the dipole-cloud infinity flux near $r_0\\approx41.66M$, while outside LISA's band for $\\mu M=0.3$, would become observable for scalar masses $\\mu M\\gtrsim0.4$, offering a possible diagnostic in a wider parameter region."],"supporting_citations":[{"why":"Independent fully numerical calculation of the same scalar radiation and fluxes; the baseline against which this paper's semi-analytical sources and fluxes are compared.","marker":"[105]"},{"why":"Supplies the Lorenz-gauge reconstructed metric perturbation $h_{\\mu\\nu}^{(0,1)}$ for circular equatorial Kerr orbits that feeds the scalar source term.","marker":"[98]"},{"why":"Introduces the modified Teukolsky formalism whose two-parameter expansion organizes the scalar and gravitational sectors for EMRIs in non-vacuum backgrounds.","marker":"[116]"},{"why":"Provides the quasi-bound scalar cloud solutions and Leaver's method used to compute the background profile $\\Phi^{(1,0)}$ and frequency $\\omega_c$.","marker":"[131]"},{"why":"Sets the flux and Noether-charge definitions, including $M_c=\\omega_c Q$ and the total scalar energy-flux formulas adopted here.","marker":"[46]"},{"why":"Supplies the Fourier-series decomposition of $\\Gamma^{-\\beta}\\bar\\Gamma^{-\\sigma}$ factors that makes the source separable in radius and angle.","marker":"[141]"}],"fun_headline_variants":["Cloud multipole decides EMRI inspiral speed","Dipole cloud slows EMRI, quadrupole speeds it","Scalar cloud's shape flips EMRI inspiral direction","EMRI in scalar cloud: dipole brake, quadrupole boost","Ultralight scalar cloud: dipole slows, quadrupole speeds EMRI"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation imports the Lorenz-gauge reconstructed metric perturbation $h_{\\mu\\nu}^{(0,1)}$ for circular equatorial Kerr orbits directly from the code of [98] without independently verifying that reconstruction, so if that external input is in error the scalar source $S_\\Phi^{(1,1)}$ and all fluxes derived from it would be affected.","fun_headline_variants_meta":{"raw":{"variants":["Cloud multipole decides EMRI inspiral speed","Dipole cloud slows EMRI, quadrupole speeds it","Scalar cloud's shape flips EMRI inspiral direction","EMRI in scalar cloud: dipole brake, quadrupole boost","Ultralight scalar cloud: dipole slows, quadrupole speeds EMRI"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001176,"raw_usage":{"total_tokens":4847,"prompt_tokens":919,"completion_tokens":3928,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":535,"completion_tokens_details":{"reasoning_tokens":3841}},"tokens_in":535,"tokens_out":3928,"duration_ms":93731,"temperature":1.0,"reasoning_tokens":3841,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T20:39:25.247328+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An independent, grid-based numerical integration of the scalar source term (without the Fourier-series separation) for the same parameters ($a=0.88M$, $\\mu M=0.3$, secondary at $r_0=10M$) should reproduce both the sign and magnitude of the net scalar energy flux: negative for the dipole cloud and positive for the quadrupole cloud. A null or opposite-sign net flux would refute the central claim.","supporting_citations":[],"review_version":1}