{"id":"ee99ea57-08fe-4ed6-a86e-907e398a7c67","arxiv_id":"2505.09696","paper_version":2,"verdict":"REJECT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"OJ287 orbital timing excludes ultralight scalar bosons with masses around 0.85 to 2.2 x 10^-21 eV through superradiance-cloud dynamical friction, independent of dark matter assumptions.","lead":"This paper uses 120 years of orbital timing of the supermassive black hole binary OJ287 to constrain ultralight bosons by the drag their superradiant clouds would exert on the secondary black hole. The result is a dark-matter-independent excluded mass window and predictions for the nanohertz gravitational wave background.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Saturated-|211⟩ assumption conflicts with the measured primary spin: Eq. (3) gives χ_sat ≈ 0.24 for μ ≈ 8.5×10^-22 eV, so the lower edge of the claimed exclusion window is internally inconsistent.","rationale":"The reader's weakest assumption focused on the wave-like dynamical-friction formula and the quasi-static cloud picture. My stress test identifies a more fundamental self-consistency failure: even granting Eq. (9) and the 0.41σ-to-1% mapping, the paper's own Eq. (3) and Table I imply that a fully saturated |211⟩ cloud at the lower end of the claimed window would leave the primary BH with a spin far below the measured χ_P2. This is not a subtle model-dependence issue but an internal contradiction within the assumptions used to compute the headline bound. The statistical mapping from 0.41σ to ⟨P_DF⟩/⟨P_GW⟩ < 0.01 is also unexplained and would need derivation, but the spin inconsistency alone invalidates the lower edge of the excluded mass range and casts doubt on the whole saturated-cloud computation. I therefore agree with the REJECT verdict, though via a different path than the reader's primary concern. The proposed method remains interesting in principle, and a corrected treatment that self-consistently evolves cloud mass and spin, plus a properly derived power-ratio limit, could still yield a useful constraint.","tokens_in":19548,"tokens_out":17030,"duration_ms":181255,"concrete_test":"Recompute the exclusion window with an explicit spin-consistency constraint: for each μ in the claimed range, solve Eq. (3) for χ_sat and require the assumed saturated-cloud state to yield χ_sat equal to χ_P2 = 0.381 ± 0.004 (or, more conservatively, evolve the cloud mass using the spin-down budget Δχ M_BH² and the age τU instead of fixing β = α). Then recompute ⟨P_DF⟩ with the resulting β; since ⟨P_DF⟩ scales linearly in the cloud density, the excluded μ interval will shrink or shift whenever the saturated assumption is violated. This directly tests whether the drag power underlying Fig. 3 is compatible with the measured OJ287 primary spin.","verdict_should_be":"REJECT","load_bearing_attack":"The central bound assumes the |211⟩ superradiant cloud is fully saturated and computes ⟨P_DF⟩ from the saturated density profile (Eq. 4 with β ≃ α). Equation (3) then fixes the saturated spin for the m=2 mode as χ_sat = 2α/(1+α²). In the claimed window μ = (8.5–22)×10^-22 eV, α ranges from 0.121 to 0.31, giving χ_sat ≈ 0.24–0.53. The OJ287 primary spin used in the orbital fit is χ_P2 = 0.381 ± 0.004 (Table I). At the lower edge, μ = 8.5×10^-22 eV (α = 0.121), χ_sat ≈ 0.238, which is more than 35σ below χ_P2. Since the paper itself states that 170Γ^-1 < τU, the cloud would have saturated long ago and spun the primary down below the value used to determine OJ287's orbital parameters. Thus the drag power entering Eq. (15) and Fig. 3 is computed for a configuration that is inconsistent with the measured spin. The text's condition 'χ < χ_211,sat' is also false for these parameters when χ is read as the measured spin. If one instead abandons saturation, β is no longer ≃ α and ⟨P_DF⟩ ∝ β is overestimated, so the derived window is not supported.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a dark-matter-independent probe of ultralight bosons (ULBs) using the dynamical friction that a superradiant boson cloud exerts on the secondary black hole in the SMBH binary OJ287. The authors compute the orbit-averaged frictional power from a saturated |211> cloud, compare it with the observed orbital decay relative to the general-relativistic quadrupole prediction, and claim a new excluded ULB mass window mu ~ (8.5-22) x 10^-22 eV. They further argue that the same dynamical friction could alleviate the final-parsec problem for SMBH binaries and imprint a characteristic turnover in the nanohertz gravitational-wave background.","tokens_in":19800,"tokens_out":6728,"duration_ms":70295,"significance":"If the central bound is correct, this would be the first DM-independent constraint on ULBs derived from binary orbital decay rather than from black-hole spin statistics, and it would open a genuinely new observational channel in the mu ~ 10^-21 eV window. The extension to the final-parsec problem and the predicted PTA turnover are testable and give the framework predictive content beyond OJ287. The paper also contains original, genuinely useful estimates in Appendix B for the survival of the |211> state against resonances and ionization. The main strengths are that the orbital parameters are taken from external timing data rather than fitted, the null-excess assumption is deliberately conservative, and the derived PTA signatures are falsifiable. However, the analysis is currently built on an internally inconsistent saturated-cloud assumption and on an unexplained choice of the power-ratio threshold, so the claimed window is not yet supported.","major_comments":[{"comment":"The saturated-cloud assumption is inconsistent with the measured primary spin at the lower edge of the claimed window. For mu = 8.5 x 10^-22 eV, alpha = 0.121, and Eq. (3) gives chi_sat = 2 alpha/(1+alpha^2) = 0.238, which is more than 35 sigma below the measured chi_P2 = 0.381 +/- 0.004. Since the paper states that 170 Gamma^-1 < tau_U throughout this band, a saturated |211> cloud would have spun the primary down to chi_sat long ago, contradicting the spin value used in the orbital fit. Thus Eq. (15) and Fig. 3 are evaluated for a configuration that cannot simultaneously have a saturated cloud and the adopted OJ287 orbital parameters. Moreover, the text condition 'chi < chi_211,sat' appears to be the reverse of the superradiance condition: superradiant growth requires chi > chi_sat, not chi < chi_sat. This issue directly undermines the lower boundary of the claimed exclusion window and must be resolved by a self-consistent treatment of the cloud mass and the primary spin.","section":"Sec. III, Eq. (3) and Table I"},{"comment":"The mapping from a 0.41-sigma null excess to the power-ratio limit <P_DF>/<P_GW> < 0.01 is not derived. A 0.41-sigma deviation is a statement about the statistical significance of the difference between observed and GR-expected power, and it does not by itself imply a 1% upper limit on an additional power component. The value R_lim = 0.01 is load-bearing because it defines the boundary of the red excluded region in Fig. 3; changing it by a factor of a few would shift or eliminate the claimed mass window. The authors should either derive R_lim from the uncertainties in [27] or state explicitly which upper-tail probability and which error propagation they use.","section":"Sec. III, paragraph after Eq. (15)"},{"comment":"The central drag calculation applies the collisionless wave dynamical-friction formula (Eq. 9) to a coherent, quasi-static superradiant cloud. Appendix B gives careful order-of-magnitude estimates for resonance and ionization depletion, but it does not demonstrate that the coherent cloud's response to the companion produces the same Coulomb logarithm C_Lambda as the wave formula for an incoherent medium. A concrete test would be to compare Eq. (9) with the gravitational-atom dynamical-friction result of Ref. [24] for the same density profile and orbital parameters; if the results differ at O(1), the quoted excluded window and the PTA predictions would need to be revised accordingly.","section":"Sec. III, Eq. (9) and App. B"}],"minor_comments":[{"comment":"The sentence 'Throughout this band the superradiance growth time obeys 170 x Gamma^-1_211 < tau_U and the primary spin satisfies chi < chi_211,sat, with chi_P2 being more conservative compared to considering chi_P1' is duplicated verbatim in the manuscript.","section":"Sec. III"},{"comment":"The caption should clarify whether the light-blue region (saturated spin below chi_P2) is part of the excluded region or a separate consistency condition, since the text does not make this clear.","section":"Fig. 3 caption"},{"comment":"The definitions of x_p and x_97 are given only in the text following the equation; the equation would be easier to read if all symbols were defined immediately in the caption or with the equation.","section":"Eq. (13)"},{"comment":"The notation 'fb1 (0.9e-21eV, 9e9M_sun)' in Fig. 5 is not explained in the caption; the reader has to infer that fb is the turnover frequency from Eq. (22). Please define these labels explicitly.","section":"Sec. V"}],"recommendation":"major_revision","confidential_remarks":"The paper depends crucially on the unpublished reanalysis in Ref. [27] for the null-excess result and for the value R_lim = 0.01. The referee report requests an explicit derivation of this threshold, which the authors should be able to provide if the central claim is to be evaluated. The duplicated paragraph in Sec. III and the reversed inequality for chi_sat suggest that the manuscript needs a careful revision pass before it can be considered further. There is no indication of problematic citation practice; the main technical obstacle is internal consistency, not novelty."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this paper's idea is good, and the writing is clear, but the central constraint as stated does not hold together. The authors propose using the orbital decay of OJ287 to bound ultralight bosons via superradiant cloud dynamical friction, without any dark matter assumption. That is a genuinely new channel, and the paper does a service by laying out the formalism and connecting it to the final-parsec problem and PTA turnovers. The appendices on resonance depletion and ionization are thoughtful.\n\nThe problem is the saturated-cloud assumption. The drag power is computed for a saturated |211> cloud, with cloud mass beta ~ alpha. But saturation fixes the primary spin to chi_sat = 2 alpha / (1 + alpha^2) (their Eq. 3). For the lower part of the claimed excluded window, say mu ~ 8.5e-22 eV (alpha ~ 0.12), chi_sat ~ 0.24, while the measured spin used in the OJ287 orbital fit is chi_P2 = 0.381 +/- 0.004. If the cloud were saturated, the spin would be 0.24, not 0.38. The paper even states 'the primary spin satisfies chi < chi_211,sat', which is false for these parameters. For most of the window, the assumed configuration is inconsistent with the data used to define the orbit. The method could be salvaged by treating the cloud mass as a free parameter beta and self-consistently relating it to the measured spin, but that is not what is done.\n\nSecond, the statistical mapping from the 0.41 sigma null excess to the 1% power-ratio limit is never derived. A 0.41 sigma deviation from zero does not obviously correspond to <P_DF>/<P_GW> < 0.01; the paper needs to show the actual error propagation. As written it looks too aggressive. Third, the text says the limit 'forbids ULB masses larger than 6.4e-22 eV' and then claims an excluded window (8.5-22)e-22 eV. Those statements are irreconcilable. There is also a verbatim duplicated paragraph.\n\nThe wave-like drag formula for a coherent bound state is a concern, but it has literature support and is less central. The citation pattern is fine. This is not a desk reject: the idea is novel, the target window is interesting, and the flaws are fixable. But the headline bound should not be cited in its current form. A serious referee should engage: the authors need to resolve the spin-saturation inconsistency, properly derive the statistical limit, and state what is actually excluded. I'd send it to review with major-revision expectations.","headline":"A genuinely new DM-independent ULB probe from OJ287, but the headline excluded window is undercut by a spin-saturation inconsistency and an unexplained statistical mapping.","tokens_in":20394,"tokens_out":12021,"would_cite":false,"duration_ms":116595,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Century-long timing of the OJ287 black-hole binary rules out ultralight bosons with masses (8.5–22)×10⁻²² eV by excluding the orbital drag their superradiant clouds would cause.","keywords":["ultralight bosons","superradiance","gravitational atoms","dynamical friction","OJ287","supermassive black hole binary","gravitational wave background","final-parsec problem"],"falsifier":"Run a self-consistent simulation of a compact companion orbiting inside a saturated |211> gravitational-atom cloud with OJ287 parameters, including resonant transitions, ionization, and back-reaction, and compute the orbit-averaged drag power; the excluded window stands only if that power is below one percent of the quadrupole gravitational-wave power across mu=(8.5–22)×$10^{-22}$ eV. A future OJ287 timing measurement that resolves an excess within that window would also falsify the null-excess premise.","tokens_in":19296,"feed_emoji":"🕳️","tokens_out":12115,"duration_ms":115214,"temperature":0.7,"pith_summary":"Ultralight bosons around $10^{-21}$ eV, if they exist, can be amplified by a spinning black hole into a surrounding boson cloud. The paper's claim is that such a cloud around OJ287's primary black hole would drag on the companion during every 12-year orbit, adding an orbital-decay power that the data do not show. Taking OJ287's century of timing as consistent with general relativity to within $0.41\\sigma$, the paper limits any extra decay to less than about one percent of the gravitational-wave power and excludes boson masses $\\mu=(8.5\\text{--}22)\\times 10^{-22}$ eV. This matters because the bound uses orbital dynamics alone, with no assumption that the bosons constitute dark matter.","feed_headline":"OJ287 orbit rules out ultralight bosons from 8.5 to 22×10⁻²² eV","feed_subtitle":"Century of orbital timing shows no extra drag, closing a dark-matter-free probe window for new bosons.","key_machinery":"The engine is the gravitational atom: an ultralight scalar field bound to a spinning black hole in hydrogen-like eigenstates $|n\\ell m\\rangle$, with the fastest-growing $|211\\rangle$ state forming a dense boson cloud of Bohr radius $r_0=1/(\\mu\\alpha)$ and density profile $\\rho=A_{211}\\,\\beta M_{\\rm BH}/r_0^3$. The companion's drag is computed with the wave dynamical-friction force $F_{\\rm DF}=4\\pi M_*^2\\rho/v^2\\, C_\\Lambda(kr_\\Lambda)$, where $v$ is the relative speed and $C_\\Lambda$ a Coulomb logarithm; averaging this force over one OJ287 orbit gives the decay power compared with the standard quadrupole gravitational-wave formula. The bound follows from requiring that this average power stay below one percent of the gravitational-wave power while the cloud is both formed ($170\\,\\Gamma_{211}^{-1}<\\tau_U$) and stable (the primary spin stays below the saturated value).","core_discovery":"The central claim is that OJ287's orbital decay is a budget: the observed period decay is consistent with quadrupole gravitational-wave emission, so any superradiant cloud's dynamical friction must contribute no more than $\\langle P_{\\mathrm{DF}}\\rangle/\\langle P_{\\mathrm{GW}}\\rangle \\lesssim 0.01$. The paper evaluates that friction for a saturated $|211\\rangle$ gravitational-atom cloud around the primary, averaging the position-dependent density over the precessing eccentric orbit, and finds the condition is violated for boson masses $\\mu\\simeq(8.5\\text{--}22)\\times 10^{-22}$ eV, while also requiring the cloud to form within cosmic time and leave the primary's spin at its observed value. It further asserts that the same drag can push supermassive-black-hole binaries through the final-parsec bottleneck in under a gigayear, and that a cosmic population of such clouds would suppress the pulsar-timing-array gravitational-wave background by roughly ten to thirty percent and imprint a spectral turnover at a calculable frequency.","pith_inferences":["The same timing-ratio method transfers to any well-measured supermassive or intermediate-mass black-hole binary: wherever observed decay matches general relativity to percent level, the mass window whose cloud Bohr radius overlaps the orbit can be excluded, making OJ287 a template rather than an isolated case.","A null detection of the predicted gravitational-wave-background turnover would not immediately falsify the boson; it would first constrain the fraction of supermassive black holes that start with near-maximal spin, since the background calculation assumes maximal initial spins.","A laboratory or cosmological discovery of a boson inside the excluded window would create a sharp testable tension: the particle would have to avoid forming a superradiant cloud on OJ287's primary, for example through self-interactions or low initial spins, and that avoidance is exactly what the paper's conservative limits leave open."],"forward_implications":["If a boson mass in the excluded window $\\mu=(8.5\\text{--}22)\\times 10^{-22}$ eV exists, OJ287's primary must either have started spinning below the superradiance threshold or lost its cloud; otherwise the measured orbit would have decayed faster than observed.","Superradiant cloud drag can carry supermassive-black-hole binaries of roughly $10^8$–$10^{10}\\,M_\\odot$ through their final sub-parsec separations in under a gigayear, offering a dark-matter-free route past the final-parsec bottleneck.","The cosmic population of such clouds suppresses the nanohertz gravitational-wave background by about ten to thirty percent and produces a turnover at $f_b\\simeq 3.2\\,\\mathrm{nHz}\\,(\\mu/10^{-21}\\,\\mathrm{eV})^3(M_{\\rm BH}/10^{10}M_\\odot)^2$, a signature pulsar-timing arrays can test.","If the claimed OJ287 orbital-decay excess is confirmed, cloud drag for $\\mu\\simeq(8.8\\text{--}9.3)\\times 10^{-22}$ eV provides a new-physics explanation; under the current conservative null reading, the drag channel alone already forbids bosons heavier than $6.4\\times 10^{-22}$ eV."],"supporting_citations":[{"why":"Supplies the measured OJ287 orbital parameters—primary and secondary masses, eccentricity, period, precession, and spin—used throughout the decay-power integrals.","marker":"[39]"},{"why":"Provides the reanalysis of OJ287 timing that finds only a 0.41σ excess over general-relativity decay, which the paper adopts as its conservative null baseline.","marker":"[27]"},{"why":"Reports the larger claimed orbital-decay excess that the paper chooses not to rely on, and which its framework would match for mu near (8.8–9.3)×10^-22 eV if confirmed.","marker":"[29]"},{"why":"Gives the analytic superradiance growth-rate approximation used to determine whether a |211> cloud can form within cosmic time.","marker":"[35]"},{"why":"Derives the collisionless wave dynamical-friction formula, Eq. (9), that converts cloud density into drag on the companion.","marker":"[52]"},{"why":"Provides the effective impact parameter for eccentric orbits used in the Coulomb logarithm when integrating the drag over OJ287's orbit.","marker":"[53]"},{"why":"Supplies the Landau-Zener resonance formalism used in Appendix B to argue that bound-state resonances deplete at most a few percent of the OJ287 cloud.","marker":"[22]"},{"why":"Shows that ionization back-reaction can be described as dynamical friction and is subdominant for massive clouds, supporting the quasi-static cloud assumption.","marker":"[24]"},{"why":"Supplies the pulsar-timing-array gravitational-wave background measurements against which the predicted suppression and turnover are compared.","marker":"[56]"}],"fun_headline_variants":["OJ287 orbit rules out ultralight boson mass window","Century of OJ287 data excludes ultralight bosons","Binary OJ287 sets sharpest bounds on light bosons yet","No drag from ultralight bosons: OJ287 constrains masses","OJ287's precision timing kills ultralight boson range"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The bound assumes the boson halo around the primary keeps its shape while the companion passes through it, so the extra orbital decay is set only by the halo's smooth density; if the companion's gravity drains, ionizes, or reshapes the halo on orbital timescales, the computed drag and the excluded mass range would not follow.","fun_headline_variants_meta":{"raw":{"variants":["OJ287 orbit rules out ultralight boson mass window","Century of OJ287 data excludes ultralight bosons","Binary OJ287 sets sharpest bounds on light bosons yet","No drag from ultralight bosons: OJ287 constrains masses","OJ287's precision timing kills ultralight boson range"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000188,"raw_usage":{"total_tokens":1310,"prompt_tokens":900,"completion_tokens":410,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":516,"completion_tokens_details":{"reasoning_tokens":319}},"tokens_in":516,"tokens_out":410,"duration_ms":4316,"temperature":1.0,"reasoning_tokens":319,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:29:00.787740+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a self-consistent simulation of a compact companion orbiting inside a saturated |211> gravitational-atom cloud with OJ287 parameters, including resonant transitions, ionization, and back-reaction, and compute the orbit-averaged drag power; the excluded window stands only if that power is below one percent of the quadrupole gravitational-wave power across mu=(8.5–22)×$10^{-22}$ eV. A future OJ287 timing measurement that resolves an excess within that window would also falsify the null-excess premise.","supporting_citations":[{"cited_title":"Authenticating the Presence of a Relativistic Massive Black Hole Binary in OJ 287 Using its General Relativity Centenary Flare: Improved Orbital Parameters","cited_arxiv_id":"1808.09309","evidence_quote":"Supplies the measured OJ287 orbital parameters—primary and secondary masses, eccentricity, period, precession, and spin—used throughout the decay-power integrals."},{"cited_title":"KLEIN-GORDON EQUATION AND ROTATING BLACK HOLES,","cited_arxiv_id":null,"evidence_quote":"Gives the analytic superradiance growth-rate approximation used to determine whether a |211> cloud can form within cosmic time."}],"review_version":1}