{"id":"6616530c-d732-4140-b615-52faff4a42b7","arxiv_id":"2411.17247","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Eccentric intermediate-mass-ratio inspirals around vector gravitational atoms acquire faster decay, stronger circularization, and negative periastron precession, making the cloud visible to LISA-like detectors.","lead":"This paper computes how a cloud of ultralight vector particles around a black hole, formed by superradiance, changes the orbit and gravitational wave signal of a small inspiraling companion. It finds the cloud speeds up orbital decay and alters periastron precession strongly enough that future space-based detectors like LISA might identify such clouds.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The off-resonant-cloud assumption is not validated for the benchmark: Table I's LZ parameters use q=1e-4, beta=0.05, while the waveform benchmark uses q=1e-3, beta=0.024 (z scales as q/beta, ~20x larger), and the inspiral must have crossed the |1011> -> |3233> resonances before entering the…","rationale":"The reader's weakest assumption correctly identifies the unperturbed-cloud condition, and my concern lives in the same place. I sharpen it in two ways. First, the paper's demonstration of nonadiabaticity uses q=1e-4 and beta=0.05, not the benchmark q=1e-3 and beta=0.024; since Appendix F itself states z ~ q/beta, the benchmark sits about a factor of 20 closer to the adiabatic regime, and the paper does not print the benchmark LZ parameters. Second, the benchmark observation window begins below x(1), so it sidesteps rather than tests the resonance crossings that a realistic inspiral would have experienced earlier. This matters because the central claim is not merely that a stationary |1011> cloud would imprint a waveform, but that such a cloud survives the inspiral in that state. The paper's qualitative caveat that adiabatic transitions can occur for sufficiently small clouds or large companions is therefore not just a footnote; it bounds the parameter region in which the headline detectability conclusion holds. My proposed check is decisive and cheap: it reuses the same matrix elements and ionization fluxes already in the paper. I do not regard this as a rejection of the paper; the benchmark may well pass the check. But the conditionality of the reader's verdict is tied to this untested point, and the check should be reported before the quantitative LISA conclusions are relied upon. The reader's verdict of CONDITIONAL remains appropriate.","tokens_in":40098,"tokens_out":22187,"duration_ms":223813,"concrete_test":"Recompute the Landau-Zener parameter z for the |1011> <-> |3233> (and n'>=2) resonances for the benchmark parameters (q=1e-3, alpha=0.03, beta=0.024, e=0.01/0.3/0.7, i=0/pi/3) and for the full Fig. 11 grid, using Eqs. (F2)-(F3) and |Omega-dot| from Pgw+Pion. For every case with z>0.01, integrate the two-level system across the resonance coupled to the orbital evolution equations (26)-(28). If any case with SNR/SNRcrit>1 shows a change in |c_1011| greater than 10% or a dephasing larger than 1/SNR, the stationary-cloud waveform model is invalid in that region of parameter space.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the |1011> cloud remains effectively static and that first-order ionization fluxes describe the dissipative dynamics depends on all bound-bound Bohr resonances being nonadiabatic and nonresonant mixing being weak. The paper's only quantitative support is Table I and Fig. 14, both computed for q=1e-4, alpha=beta=0.05, and, for Fig. 14, equatorial circular orbits. The waveform benchmark instead has q=1e-3, beta=0.024 (beta=0.813 alpha), with initial eccentricities up to 0.6 and inclination i=2*pi/3. Appendix F states that for q<<1, z is proportional to q/beta, so the benchmark's z values are roughly 20 times larger than those in Table I. More importantly, the observed one-year window starts at xa about 0.54, below the first Bohr resonance x(1) about 1.6; an astrophysical IMRI formed outside the cloud must have crossed the |1011> <-> |3233> resonances at x about 1.71 and 2.72 earlier. If any of those crossings is adiabatic, which is possible for m2 significantly above 10 solar masses or beta appreciably below 0.024, the cloud would be transferred or depleted before the LISA band, invalidating the stationary potential and the ionization fluxes used to generate the waveforms. Nonresonant mixing is likewise only estimated for equatorial circular orbits, not for the benchmark's inclined eccentric orbit, so the quantitative claim rests on an extrapolation beyond the regime that was checked.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the orbital dynamics and gravitational-wave signals of intermediate mass-ratio inspirals (IMRIs) around a vector ultralight-boson cloud formed by superradiance, focusing on the fastest-growing |1011> state. Using a Newtonian perturbative treatment, the authors consistently combine the conservative effect of the stationary cloud potential and gravitomagnetic field with the dissipative effect of cloud ionization computed via Fermi's golden rule. They derive scaling relations for the ionization fluxes, compute secular changes of the orbital elements (including a negative periastron precession contribution), and evolve adiabatic inspirals with the analytical kludge waveform model. For a benchmark system with m1=10^4 solar masses, m2=10 solar masses, alpha=0.03, beta=0.0244, and d=1 Mpc, they find that the cloud accelerates orbital decay and circularization, and that one- and five-year LISA waveforms have nearly zero overlap with vacuum waveforms (faithfulness ~ -10^-3) with SNR above 100. They conclude that such vector gravitational atoms may be detectable by LISA and Taiji.","tokens_in":115,"tokens_out":14128,"duration_ms":187059,"significance":"The paper makes a solid step beyond earlier order-of-magnitude estimates of dynamical friction in gravitational atoms by computing ionization fluxes and their backreaction from first principles, and by including conservative precession effects. The scalar-reduction identity for the |1011> state (Eq. 29) is elegant and allows the use of established scalar-atom ionization results. The clear scaling laws (Eqs. D4-D6 and the q-scaling discussion) make the model easy to test, and the paper is careful to compare the magnitudes of all forces in Fig. 4. The central claim—that eccentric IMRI waveforms in a |1011> cloud are distinguishable from vacuum waveforms—is falsifiable and, if correct, would provide a new channel to probe ultralight vector bosons. However, the validity of the claim depends crucially on the cloud surviving the inspiral in its |1011> state, and this is where the validation is incomplete.","major_comments":[{"comment":"The central assumption that the |1011> cloud remains an unperturbed single eigenstate is checked in Table I and Fig. 14 only for parameters different from the waveform benchmark. Table I uses alpha=beta=0.05 and q=10^-4, and Fig. 14 uses q=10^-4, while the benchmark in Sec. IV uses q=10^-3 and beta=0.813 alpha=0.0244. Since Appendix F states that z is proportional to q/beta when q<<1 and P is approximately P_ion, the benchmark's Landau-Zener parameters are about 20 times larger than the tabulated values; for the specific |3233> entries this still gives z of order 10^-3 or smaller, but no LZ calculation is shown for the other allowed Bohr resonances (|2122>, |21j'0>, |n'011>) or for the eccentricity values of the benchmark. Moreover, the benchmark inspiral starts at x_a ~ 0.54, which is below the first Bohr resonance at x^(1) ~ 1.6; the paper's assumption that the inspiral starts with x_a < x^(1) is imposed, not derived. A companion formed outside the cloud must cross the resonances at x ~ 2.72 and 1.71 before reaching the LISA band. The absence of a direct validation at the benchmark parameters is a load-bearing gap: if any of these crossings is adiabatic, the cloud would be depleted or transferred, invalidating the potential and ionization fluxes used to generate the waveforms and the claimed detectability.","section":"Sec. IV and Appendix F"},{"comment":"The nonresonant bound-state mixing is estimated only for an equatorial circular orbit. The paragraph containing Eq. (F8) explicitly states 'In the case of equatorial corotating (i=0) circular orbit,' and Fig. 14 shows P_mix and (M_dot_c)_mix only for that case. The benchmark waveforms in Sec. IV, however, use eccentric orbits with e0 up to 0.85 (Fig. 9) and e0=0.6 (Fig. 10). The ionization flux is independent of inclination because of the scalar |100> reduction, but the bound-bound mixing is not, as Appendix F itself notes that a dedicated study would be required for a generic inclined eccentric orbit. Applying the stationary-cloud model to the full benchmark parameter space without quantifying nonresonant mixing there is an extrapolation, not a result. The authors should compute or bound this effect at the benchmark point, or restrict the claims to the equatorial circular case.","section":"Appendix F and Sec. IV"}],"minor_comments":[{"comment":"The abstract describes the cloud as being 'saturated in its SR ground state'; for vector fields the fastest-growing superradiant mode is |1011>, which is not the n=1, l=0 hydrogenic ground state. Please add a brief clarification to avoid confusion.","section":"Abstract and Sec. II"},{"comment":"The formula alpha = 0.1 x (mu/1.3 x 10^-11 eV)(M/M_sun) mixes different unit conventions; please state the conversion constant explicitly in natural units.","section":"Eq. (40)"},{"comment":"The upper-left panel of Fig. 6 is mentioned in the text as showing the magnitudes of the precession rates relative to BH-induced rates, but the plotted quantities and their color coding are not described in the caption. Please expand the caption.","section":"Fig. 6"},{"comment":"Please double-check the sign of the tau_ion,X term in the inclination evolution equation; in the limit of a spherically symmetric cloud the X-component should vanish, and the relation between Eqs. (23) and (24) should be stated consistently.","section":"Sec. III.A, Eq. (24)"},{"comment":"The term 'nondimensional' is used throughout; 'dimensionless' is the more common English term. No substantive issue.","section":"Throughout"},{"comment":"The statement that the paper corrects the treatment of the cloud depletion effect in Appendix E should explicitly note that this corrects the authors' own previous work [83]; as written, it is unclear which earlier treatment is being corrected.","section":"Introduction"}],"recommendation":"major_revision","confidential_remarks":"The paper is appropriate for a general-relativity or astrophysics journal. The main risk is that the detection claim may be overstated if the cloud is destroyed by resonances before entering the LISA band. I recommend requesting the missing Landau-Zener and nonresonant-mixing computations at the benchmark parameters, and a discussion of the formation-history assumption. The emphasis on the authors' own prior work is acceptable but should be moderated."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper does something useful: it takes the established ionization-flux machinery for gravitational atoms and applies it to the vector |1011> cloud in eccentric intermediate-mass-ratio inspirals, adding the conservative effect of the cloud's Newtonian potential. The new pieces are the eccentric-orbit ionization fluxes, the negative periastron precession from the spherical potential, and the waveform/Fisher forecasts for LISA and Taiji. The derivation is internally consistent: alpha and q are small, orbital radii stay above relativistic scales, and the magnitude comparisons in the appendices justify neglecting oscillatory potentials, cloud depletion, and resonant transitions in the regime they check. They also correct a mistake from their earlier circular-orbit paper regarding the cloud-depletion effect, which is good practice.\n\nThe central result—that the cloud accelerates decay and circularization and imprints a negative periastron precession that makes the waveform essentially orthogonal to the vacuum waveform—is likely correct under the stated assumptions.\n\nThe main soft spot is the off-resonance assumption for the benchmark. The LZ parameters in Table I are computed for q=1e-4, beta=0.05, while the waveform benchmark uses q=1e-3, beta=0.024, which pushes z up by roughly a factor of twenty. The benchmark's one-year window starts at xa~0.54, below the first Bohr resonance at ~1.6, so an astrophysical inspiral must have crossed the |1011>->|3233> resonances earlier. The paper states it focuses on off-resonant systems, but it does not demonstrate that the benchmark system stays off-resonant through those crossings. That is a genuine gap, though probably fixable with a short calculation. The nonresonant mixing estimate is also only for equatorial circular orbits, not the benchmark's inclined eccentric orbit. These are not fatal, but they mean the detectability claim is extrapolated slightly beyond the checked regime.\n\nA separate limitation is that no code or data are provided for the P functions and the ionization sums, and the Fisher forecast fixes several extrinsic parameters, which is standard but optimistic. There is no event-rate estimate, so the benchmark source at 1 Mpc is an existence proof, not a projection.\n\nFor a reader in gravitational-wave astronomy or ultralight-boson phenomenology, this is worth a careful read. I would send it to a serious referee, asking specifically for a check of the resonance adiabaticity at the benchmark parameters and for release of the numerical machinery. The paper deserves a round of revision, not a desk reject.","headline":"A credible extension of gravitational-atom physics to vector clouds in eccentric IMRIs; the detectability claim rests on an off-resonance assumption checked only for parameters different from the benchmark.","tokens_in":41018,"tokens_out":2918,"would_cite":true,"duration_ms":28125,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A vector boson cloud grown around a black hole would leave detectable, near-orthogonal imprints in the gravitational waves of eccentric intermediate-mass-ratio inspirals.","keywords":["vector gravitational atoms","ultralight vector bosons","superradiance","intermediate-mass-ratio inspirals","eccentric orbits","gravitational wave astronomy","cloud ionization","LISA"],"falsifier":"A concrete check: search LISA data for an eccentric IMRI with $m_1=10^4\\,M_\\odot$, $m_2=10\\,M_\\odot$, $\\alpha=0.03$, $\\beta=0.8\\alpha$, $f_0=2\\times10^{-4}\\,\\mathrm{Hz}$, $e_0=0.6$, and $d=1\\,\\mathrm{Mpc}$. If the recovered waveform matches a vacuum binary with faithfulness above $0.99$ — positive Schwarzschild periastron precession, braking index in $[4/3,11/3]$, and no ionization-accelerated decay or circularization — the claim that this vector cloud dominates the waveform in that band is falsified for that event. A companion numerical check would be to evolve the cloud self-consistently for parameters where the Landau-Zener parameter approaches unity and see whether resonances deplete the cloud before the predicted dephasing accumulates.","tokens_in":39840,"feed_emoji":"🛰️","tokens_out":8095,"duration_ms":73697,"temperature":0.7,"pith_summary":"The paper sets out to show that an ultralight vector boson cloud formed around a spinning black hole — a vector gravitational atom in its fastest-growing $|1011\\rangle$ ground state — is not a passive backdrop for an inspiraling companion. In the leading Newtonian regime with small mass ratio, the cloud's stationary gravity produces a negative secular periastron precession at separations near the gravitational Bohr radius, while ionization of the cloud by the companion's tidal field accelerates orbital decay and circularization. The paper computes the resulting adiabatic orbital evolution and gravitational waveforms, and finds that for a benchmark intermediate-mass-ratio inspiral ($m_1=10^4\\,M_\\odot$, $m_2=10\\,M_\\odot$, $\\alpha=0.03$, $\\beta=0.8\\alpha$, $d=1\\,\\mathrm{Mpc}$) the one-year and five-year overlaps between cloud and vacuum waveforms are essentially zero ($\\mathcal{F}\\sim -10^{-3}$) at signal-to-noise ratios of order $10^2$. The conclusion a sympathetic reader is asked to accept: future space-based gravitational-wave detectors could identify a vector gravitational atom through the orbital and waveform imprints of an eccentric intermediate-mass-ratio inspiral.","feed_headline":"A vector boson cloud rewrites eccentric IMRI waveforms","feed_subtitle":"For a benchmark 10^4-solar-mass system, LISA overlaps with vacuum fall to ~-10^-3 at SNR ~100.","key_machinery":"The load-bearing object is the vector gravitational atom in the $|1011\\rangle$ state, the fastest-growing superradiant mode, whose nonrelativistic wave function is proportional to the scalar $|100\\rangle$ hydrogenic state, $\\Psi_{1011}=-\\frac{1}{\\sqrt2}(1,i,0)^T\\Psi_{100}$. The paper's machinery is a first-order perturbative model of the companion's orbit: Gaussian osculating-element equations driven by the conservative acceleration $-\\nabla\\Phi + v\\times(\\nabla\\times\\Xi)$ from the cloud's stationary Newtonian and gravitomagnetic potentials, together with dissipative ionization fluxes obtained from Fermi golden rule bound-free transition rates under the companion's tidal potential decomposed into multipoles. The ionization fluxes carry sharp jumps at semimajor axes $x^{(g)}=(2g)^{2/3}(1+q)^{1/3}$, and the secular evolution is fed into the analytical-kludge quadrupole waveform with LISA noise to produce faithfulness, signal-to-noise, and Fisher-matrix statements.","core_discovery":"The central claim is that the $|1011\\rangle$ vector cloud leaves a distinctive, observationally separable mark on eccentric intermediate-mass-ratio inspirals. Treating the cloud as a nearly unperturbed single eigenstate, the paper combines two effects: a conservative one from the cloud's Newtonian potential, which yields a negative contribution to the secular periastron precession that can rival or exceed the Schwarzschild precession at $x_a\\sim 1$, and a dissipative one from ionization, whose backreaction produces energy and angular-momentum fluxes that accelerate orbital decay and circularization. These effects make the binary more eccentric than a vacuum inspiral at a given orbital frequency and can drive the braking index below $4/3$ or even negative. For the benchmark system the cloud and vacuum waveforms are almost orthogonal ($\\mathcal{F}\\approx -10^{-3}$ over one and five years) with $\\mathrm{SNR}\\sim 10^2$, and Fisher-matrix estimates give fractional errors below unity for the cloud parameters $\\alpha$ and $\\beta$, so a detection would both reveal the cloud and constrain the boson.","pith_inferences":["Beyond the paper: if such an event is observed, the detection would give a direct measurement of the boson mass from the Bohr-radius relation, independent of the cloud's own monochromatic gravitational-wave line.","Beyond the paper: the strong eccentricity dependence shown in the paper's parameter scans suggests that searches should prioritize the most eccentric events, for which the cloud-vs-vacuum overlap drops fastest and the SNR difference grows largest.","Beyond the paper: the Landau-Zener parameters in Table I indicate that the nonadiabaticity assumption can break down for larger companion masses or smaller cloud masses, so the predicted waveforms and detectability window may close in that region of parameter space.","Beyond the paper: a detected cloud signal would not by itself identify the boson spin, since the dominant ionization and potential effects are shared with scalar $|100\\rangle$ clouds; distinguishing the spin would require the gravitomagnetic precession or the cloud's own GW line."],"forward_implications":["For an eccentric IMRI with $m_1=10^4\\,M_\\odot$, $m_2=10\\,M_\\odot$, $\\alpha=0.03$, and $\\beta=0.8\\alpha$, one-year or five-year LISA observations can distinguish the vector-cloud waveform from a vacuum inspiral, since $\\mathrm{SNR}\\sim 10^2$ and the faithfulness $F$ is about $-10^{-3}$.","The cloud makes the binary more eccentric at a given orbital frequency and can produce a small or negative braking index, in contrast to the vacuum range $n_b\\in[4/3,11/3]$.","Ionization accelerates orbital decay and circularization, and the eccentricity decays more slowly relative to the semimajor axis, shifting the inspiral track away from the vacuum trajectory.","Fisher-matrix estimates show that $\\alpha$ and $\\beta$ can be measured with fractional errors below unity; in the absence of a cloud, a one-year observation would constrain $\\beta\\lesssim 3\\times10^{-6}$ for fixed $\\alpha=0.03$.","Because the ionization fluxes and stationary potential of the vector $|1011\\rangle$ state are degenerate with those of scalar $|100\\rangle$ and certain vector and tensor states, the same waveforms can be used to search for a broader class of ultralight boson clouds."],"supporting_citations":[{"why":"Supplies the superradiant growth rates and identifies the vector $|1011\\rangle$ state as the fastest-growing mode that saturates into the cloud.","marker":"[66]"},{"why":"Provides the selection rules, transition matrix elements, and Landau-Zener resonance treatment used to justify neglecting bound-bound transitions.","marker":"[87]"},{"why":"Establishes the ionization model for gravitational atoms whose Fermi golden rule bound-free fluxes are adopted here.","marker":"[97]"},{"why":"Computes dynamical friction and ionization fluxes for inclined eccentric orbits, extended in this paper to the vector $|1011\\rangle$ cloud.","marker":"[99]"},{"why":"Provides the earlier outspiral and cloud-depletion treatment that this paper corrects and replaces with ionization fluxes and a consistent potential.","marker":"[83]"},{"why":"Gives the relativistic and hyperfine spectrum of gravitational atoms used for level corrections and resonance estimates.","marker":"[114]"},{"why":"Supplies the accurate gravitational-wave depletion power of the $|1011\\rangle$ cloud used for cloud age and mass evolution.","marker":"[124]"},{"why":"Gives the periastron-advance signature and parameter-estimation relation used in the analytical-kludge waveform and Fisher estimates.","marker":"[126]"},{"why":"Provides the analytical-kludge waveform construction used to generate the binary's time-domain and frequency-domain signals.","marker":"[133]"},{"why":"Provides the LISA noise power spectral density used for SNR, faithfulness, and Fisher-matrix calculations.","marker":"[136]"}],"fun_headline_variants":["Eccentric IMRIs reveal vector cloud signatures","Vector boson clouds reshape eccentric inspiral waves","Gravitational atom imprints on eccentric IMRI orbits","Eccentric IMRIs offer new test of ultralight bosons"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes the small companion leaves the $|1011\\rangle$ cloud essentially unperturbed: all bound-bound Bohr resonances stay nonadiabatic (Landau-Zener parameter $z\\ll1$ as computed in Table I), nonresonant bound-state mixing is negligible (Appendix F), and the cloud mass changes by only about $10^{-2}$ during the inspiral, so the stationary single-eigenstate potential and first-order ionization fluxes describe the entire evolution.","fun_headline_variants_meta":{"raw":{"variants":["Eccentric IMRIs reveal vector cloud signatures","Vector boson clouds reshape eccentric inspiral waves","Gravitational atom imprints on eccentric IMRI orbits","Eccentric IMRIs offer new test of ultralight bosons"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000248,"raw_usage":{"total_tokens":1564,"prompt_tokens":978,"completion_tokens":586,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":594,"completion_tokens_details":{"reasoning_tokens":518}},"tokens_in":594,"tokens_out":586,"duration_ms":6044,"temperature":1.0,"reasoning_tokens":518,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:21:15.857609+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete check: search LISA data for an eccentric IMRI with $m_1=10^4\\,M_\\odot$, $m_2=10\\,M_\\odot$, $\\alpha=0.03$, $\\beta=0.8\\alpha$, $f_0=2\\times10^{-4}\\,\\mathrm{Hz}$, $e_0=0.6$, and $d=1\\,\\mathrm{Mpc}$. If the recovered waveform matches a vacuum binary with faithfulness above $0.99$ — positive Schwarzschild periastron precession, braking index in $[4/3,11/3]$, and no ionization-accelerated decay or circularization — the claim that this vector cloud dominates the waveform in that band is falsified for that event. A companion numerical check would be to evolve the cloud self-consistently for parameters where the Landau-Zener parameter approaches unity and see whether resonances deplete the cloud before the predicted dephasing accumulates.","supporting_citations":[{"cited_title":"Proposal for Determining the Total Masses of Eccentric Binaries Using Signature of Periastron Advance in Gravitational Waves","cited_arxiv_id":"astro-ph/0111107","evidence_quote":"Gives the periastron-advance signature and parameter-estimation relation used in the analytical-kludge waveform and Fisher estimates."}],"review_version":1}