{"id":"4236f959-e042-49d2-819e-40cc7850a44c","arxiv_id":"2608.09811","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"LISA spin measurements could exclude scalar and vector ultralight bosons over roughly four orders of magnitude in mass, while post-merger follow-up searches offer a narrower, model-dependent window for detecting vector boson clouds.","lead":"This paper calculates how often LISA, the planned space-based gravitational wave detector, could catch signs of ultralight bosons by measuring the spins of massive black hole mergers and by listening for gravitational wave signals from particle clouds around the merger remnants. It finds that spin measurements could rule out a wide range of boson masses, while follow-up searches could detect a narrower range, with the odds depending strongly on how black holes grow.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Sec. III C follow-up exclusion curves ignore the vector boson's spin-down of binary components, which the paper itself models in Sec. III D; this makes Table III and Fig. 3 follow-up exclusion ranges optimistic.","rationale":"The reader's verdict is CONDITIONAL and already flags, in its rationale, that follow-up exclusion curves are computed under a no-boson assumption while the detection analysis accounts for spin-down. However, the reader's stated weakest_assumption is the accretion timescale, not this internal inconsistency. I regard the Sec. III C / III D inconsistency as the single most load-bearing concern because it affects a central quantitative result (the vector follow-up exclusion ranges in Table III and Fig. 3) and is a self-consistency issue within the paper's own framework, rather than an external astrophysical assumption. The accretion-timescale concern is real but is explicitly acknowledged and partially stress-tested by the 0.01 T_S comparison; by contrast, the follow-up exclusion inconsistency is untested and the paper itself supplies the machinery needed to fix it. The spin-measurement exclusions, which are the dominant results, are not affected by this issue, so the overall verdict of conditional acceptance remains appropriate. The revision should either recompute the follow-up exclusion curves with the Sec. III D spin-down prescription or clearly label them as upper bounds on exclusion power.","tokens_in":24699,"tokens_out":7078,"duration_ms":67793,"concrete_test":"Recompute the vector follow-up exclusion probability for the Q3nodelays model across the Table III range (e.g., mu = 1e-16 eV) using the Sec. III D pipeline: apply superradiant spin-down to all binary components before merger, compute remnant mass and spin with Eqs. (54) and (57) of Ref. [50], then calculate the SNR >= 10 probability. If the corrected exclusion probability drops below 0.99 at any mass where Table III claims exclusion, the follow-up exclusion column is overstated and must be revised or reframed as an upper bound.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The follow-up exclusion analysis in Sec. III C is internally inconsistent with the detection analysis in Sec. III D. For each remnant, Sec. III C uses the mass and spin from the original no-boson population catalogs and asks whether a vector-boson cloud around that remnant would produce SNR >= 10; non-detection is then counted as excluding that boson mass. But if a vector boson of that mass exists, superradiance acts on the binary components before merger, reducing their spins and therefore the remnant spin. This is exactly what Sec. III D does when computing detection probabilities, and the paper states it explicitly: 'A key limitation arises from the fact that the same superradiant instability responsible for generating post-merger signals also acts on the binary components prior to merger... leading to remnants with significantly reduced spins' (Sec. IV). The right panel of Fig. 2 shows this spin reduction, and Fig. 6 shows that follow-up detection relies on remnants with spins 0.55-0.85, while the exclusion population in the left panels peaks at spins above 0.7-0.9. Consequently, absence of a follow-up signal from a hypothetical high-spin remnant cannot exclude a boson mass that would have prevented that remnant from being high-spin in the first place. The vector follow-up exclusion curves in Fig. 3 and the 'Vector via Follow-up' column of Table III are therefore optimistic and should be recomputed with the Sec. III D procedure or explicitly reframed as upper bounds. This is a correctness issue inside the paper's own modeling framework, not merely an astrophysical uncertainty.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper forecasts the ability of LISA to constrain or detect ultralight scalar and vector bosons through black-hole superradiance, using three massive black-hole population models (PopIII, Q3delays, Q3nodelays) and the open-source SuperRad package. Two observational channels are considered: spin measurements of merging massive black-hole binaries, and follow-up searches for quasi-monochromatic gravitational-wave signals from superradiant clouds around merger remnants. The paper reports that spin measurements can exclude scalar masses approximately in [5e-18, 1e-14] eV and vector masses approximately in [6e-19, 2e-14] eV, while follow-up searches are sensitive to a narrower vector mass range around [3e-17, 3e-15] eV, with detection probabilities that can reach near unity for Q3nodelays near 1e-16 eV. The spin-measurement exclusion logic is clearly described and tested against variations in SNR threshold, spin uncertainty, and spin-down timescale; however, the follow-up exclusion analysis is internally inconsistent with the paper's own treatment of pre-merger spin-down.","tokens_in":24994,"tokens_out":5799,"duration_ms":54795,"significance":"If the forecasts are correct, this is a useful and timely projection for LISA's new-physics reach. The paper's strengths are its use of a publicly available, physically detailed superradiance model (SuperRad), the adoption of updated population catalogs that incorporate pulsar-timing-array information, the inclusion of both scalar and vector bosons, and the explicit testing of SNR and spin-uncertainty assumptions. The spin-measurement constraints are coherent and largely consistent with earlier work in the literature. The main significance-limiting issue is the follow-up exclusion calculation, which currently overstates the constraining power of remnant searches because it does not account for the very spin-down effect the paper itself models elsewhere. This issue is fixable and should not obscure the value of the spin-based forecasts.","major_comments":[{"comment":"The follow-up exclusion calculation in Sec. III C is internally inconsistent with the detection calculation in Sec. III D. In Sec. III C, each remnant is taken directly from the no-boson catalogs, and a vector mass is counted as excluded if SuperRad predicts a follow-up SNR >= 10 for a cloud around that remnant. But Sec. III D and the limitation stated in Sec. IV show that the same instability acts on the binary components prior to merger, producing remnants with significantly reduced spins (right panel of Fig. 2 and right panels of Fig. 6). For the masses in the 'Vector via Follow-up' column of Table III, the high-spin remnants on which the exclusion is based would generally not exist if the boson were present, so the absence of a follow-up signal cannot exclude those masses. The correct procedure would be to apply the Sec. III D spin-down to the binary components, recompute the remnant mass and spin, and then evaluate the follow-up SNR; if the resulting signal is undetectable, the non-detection is consistent with the boson and should not be counted as an exclusion. As written, the vector follow-up curves in Fig. 3 and Table III are optimistic and should be recomputed or explicitly relabeled as applying only to the no-boson catalogs.","section":"Sec. III C and Sec. III D"},{"comment":"The spin-measurement exclusion ranges in Table III rest on the assumption that accretion re-spins massive black holes on the Salpeter timescale T_S = 4.5e7 yr, so that any instability slower than this leaves no observable spin-down signature. The 0.01 T_S curve in Fig. 3 is a useful robustness check, but it only changes the threshold on the superradiant growth time; it does not model faster or episodic coherent accretion, partial spin-down followed by re-spin, or the Blandford-Znajek spin-down discussed in Sec. II C. The abstract and Table III quote the T_S-based ranges without reporting the spread across these accretion assumptions, which makes the claimed exclusion ranges appear more robust than the modeling uncertainty supports. I recommend either reporting the 0.01 T_S ranges alongside the fiducial ranges in Table III and the abstract, or adding an explicit caveat that the quoted ranges are conditional on the Salpeter-time accretion assumption.","section":"Sec. II C"}],"minor_comments":[{"comment":"The abstract quotes the scalar range [5e-18, 1e-14] eV and the vector range [6e-19, 2e-14] eV, which correspond to the Q3nodelays model in Table III; the ranges for Q3delays and PopIII differ by almost an order of magnitude. The abstract should either identify these as the most optimistic model values or quote the model dependence explicitly.","section":"Abstract and Table III"},{"comment":"The catalog-realization convergence criterion is only stated as a 'standard error of order 0.1'; specifying the monitored quantity (e.g., mean number of SNR>8 mergers) and the acceptable tolerance would make the convergence test easier to reproduce.","section":"Sec. III A"},{"comment":"The detection-probability discussion states that Q3nodelays detection probability reaches about 80% in [2e-16, 1.5e-15] eV, but Fig. 4 appears to show the curve crossing 0.8 only near the upper end of that interval; a clearer reading of the threshold would help.","section":"Sec. IV and Fig. 4"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the spin-based forecasts are substantial. The follow-up exclusion inconsistency is the main blocking issue: it affects a headline claim in the abstract and Table III. Once that is addressed, the manuscript should be suitable for publication. I would not reject on the basis of the current issues, but the revision needs to be substantive rather than cosmetic."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing you should know: this is a solid forecasting paper for LISA's ability to probe ultralight bosons, but one of its headline numbers — the vector-boson follow-up exclusion ranges in Fig. 3 and Table III — is optimistic because the analysis in Sec. III C uses remnant spins from catalogs with no boson, ignoring the pre-merger spin-down that the paper's own detection analysis in Sec. III D accounts for. The paper explicitly acknowledges that the same instability spins down binary components and reduces remnant spins, so the follow-up exclusion curves are upper bounds at best. This should be fixed or reframed.\n\nWhat is genuinely new: the vector-boson follow-up forecasts, the self-consistent detection calculation that includes pre-merger spin-down, and the use of updated population catalogs calibrated to pulsar timing array data. The scalar spin-measurement constraints reproduce earlier results (Refs. [18, 64]) with the new catalogs — not a new result, but a useful consistency check. The authors ship an open-source code, SuperRad, so the spin-down and strain calculations are reproducible. The spin-measurement exclusion logic is tested against SNR thresholds, spin uncertainties, and two spin-down timescales (T_S and 0.01 T_S), and it is robust to those variations. That part is trustworthy.\n\nThe soft spots are the follow-up exclusion inconsistency and, secondarily, the accretion timescale assumption. The paper assumes T_S = 4.5e7 yr as the fiducial upper limit for spin-down; real accretion histories could be faster or more episodic. The 0.01 T_S test only partially covers that, but the paper is explicit about the limitation, so it is a caveat, not a hidden flaw.\n\nOne more thing: the detection probability calculation is done correctly with pre-merger spin-down, and the results show direct detection is harder than the exclusion curves imply. That asymmetry is a sign the authors understood the issue but didn't propagate it into the exclusion numbers. A referee should ask them to reconcile the two.\n\nBottom line: this deserves a serious referee. I would accept it conditionally, with the follow-up exclusion curves recomputed or explicitly labeled as upper bounds. The spin-measurement forecasts alone justify publication.","headline":"Solid LISA forecast paper, but the vector follow-up exclusion curves ignore the pre-merger spin-down that the paper itself models, making those exclusion ranges optimistic.","tokens_in":25568,"tokens_out":2569,"would_cite":true,"duration_ms":20979,"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":"LISA observations of merging massive black holes could exclude ultralight scalar and vector bosons over roughly four orders of magnitude in mass, and targeted follow-up of merger remnants could detect vector boson clouds near 1e-16 eV.","keywords":["ultralight bosons","black hole superradiance","LISA","gravitational waves","black hole spins","massive black hole binaries","boson clouds","follow-up searches"],"falsifier":"Take any LISA-detected massive binary whose confidently measured primary spin exceeds $a_{\\rm max}$ for a boson mass inside the claimed exclusion window under the fiducial 10% spin uncertainty and $T_S$; one such event would already break the exclusion for that mass. Conversely, if a vector boson near $10^{-16}$ eV exists and the Q3nodelays population is accurate, a four-year mission with no post-merger cloud signal at SNR $\\geq 10$ would contradict the paper's near-unity detection probability.","tokens_in":24479,"feed_emoji":"🛰️","tokens_out":7333,"duration_ms":59000,"temperature":0.7,"pith_summary":"This paper forecasts what the space-based gravitational-wave observatory LISA could learn about ultralight bosons, hypothetical particles motivated by dark-matter and string-theory scenarios, by watching superradiant instabilities drain angular momentum from rotating massive black holes. The authors argue that precise spin measurements of black holes in merging binaries can rule out scalar boson masses in roughly $[5\\times10^{-18},10^{-14}]$ eV and vector boson masses in $[6\\times10^{-19},2\\times10^{-14}]$ eV, with the exact window depending on the assumed massive-black-hole population. They further show that targeted follow-up searches for the quasi-monochromatic gravitational waves emitted as a boson cloud dissipates around a merger remnant can constrain vector masses in a narrower $[3\\times10^{-17},3\\times10^{-15}]$ eV band, and that if a vector boson near $10^{-16}$ eV exists, detection prospects range from negligible to near-certain across population models. This matters because a null result would sharpen constraints on physics beyond the Standard Model, while a detection would be direct evidence for new ultralight particles.","feed_headline":"LISA spin data can exclude ultralight bosons across four mass decades","feed_subtitle":"A four-year LISA mission could probe boson masses down to 1e-18 eV via black hole spins and remnant follow-up.","key_machinery":"The machinery is black-hole superradiance: a bosonic field with mass $\\mu$ and azimuthal mode $m$ extracts energy and angular momentum from a Kerr black hole when $\\omega_R < m\\Omega_H$, growing into a cloud until saturation. The controlling parameter is the gravitational fine-structure constant $\\alpha = M\\mu$ (black-hole mass times boson mass in geometric units), which sets growth rates, spin-down to a maximum allowed spin $a_{\\rm max}(M,\\mu,\\tau_{\\rm sd})$, and the quasi-monochromatic gravitational-wave strain from the dissipating cloud. The authors use the SuperRad model to evolve all unstable modes and compute these quantities, and set the spin-down timescale $\\tau_{\\rm sd}$ against the Salpeter accretion timescale $T_S = 4.5\\times10^7$ yr to decide whether a spin-down signature survives. The population catalogs then supply which binaries LISA would see and which remnants could host detectable clouds.","core_discovery":"The central claim is that LISA, through two complementary channels, can exclude and potentially detect ultralight bosons via black-hole superradiance. Using three massive-black-hole population models (light-seed PopIII and heavy-seed Q3 with and without delays), the authors find that spin measurements of binaries with SNR $\\geq 20$ can exclude, with probability greater than 0.99, scalar masses in $[5\\times10^{-18},10^{-14}]$ eV and vector masses in $[6\\times10^{-19},2\\times10^{-14}]$ eV, assuming a Salpeter spin-down timescale $T_S = 4.5\\times10^7$ yr. Follow-up searches of post-merger remnants, requiring SNR $\\geq 10$ cloud signals, exclude only vector masses in $[3\\times10^{-17},3\\times10^{-15}]$ eV; scalar clouds grow too slowly to be seen. When the backreaction of a real vector boson on binary spins is included, the detection probability exceeds roughly 80% for the Q3nodelays model in $[2\\times10^{-16},1.5\\times10^{-15}]$ eV, while other models give lower but non-negligible prospects. The key tension identified is that the same superradiance that creates the cloud spins down the progenitor black holes before merger, suppressing remnant spins and weakening the very signals follow-up searches target.","pith_inferences":["The paper's fiducial Salpeter-timescale assumption brackets only two spin-down timescales ($T_S$ and $0.01T_S$); if accretion is systematically slower than assumed, the excluded ranges would extend further, and if faster, they would shrink, so the quoted numbers should be read as model-dependent rather than hard physical bounds.","The same spin-down formalism could be folded into a joint Bayesian analysis of all LISA events to constrain boson mass and population parameters simultaneously, which the paper notes as future work but does not carry out.","Because the exclusions rely on pure gravitational coupling, adding self-interactions or dark-photon kinetic mixing above the rough thresholds the paper estimates ($f \\gtrsim 10^{15}{-}10^{16}$ GeV, $\\epsilon \\lesssim 0.001$) could suppress cloud growth and weaken both exclusions and detection prospects.","The mass-spin feature shown for a $5\\times10^{-16}$ eV vector boson suggests that even without a detected cloud signal, the statistical shape of the LISA spin distribution could itself be evidence for superradiance, a testable signature beyond single-event exclusion."],"forward_implications":["For all three population models, spin measurements exclude scalar and vector bosons over roughly four orders of magnitude in mass, with PopIII shifting to higher masses and Q3nodelays giving the broadest window.","Follow-up searches cannot constrain scalar bosons, whose cloud growth times exceed the LISA mission, but they add a complementary vector-boson exclusion window around $10^{-16}$ to $10^{-15}$ eV.","If an ultralight vector boson exists in $[10^{-16},2\\times10^{-15}]$ eV, the probability that LISA sees at least one post-merger cloud signal ranges from very small to near unity depending on the population; Q3nodelays is the most optimistic, reaching above 80%.","Superradiant spin-down of binary components before merger shifts the remnant spin distribution downward, which reduces follow-up detectability; this backreaction is essential to the detection forecast.","Constraints are robust to changes in SNR threshold and spin-measurement accuracy (1% vs 10%), but shift by up to an order of magnitude with the assumed spin-down timescale and population model."],"supporting_citations":[{"why":"Supplies the SuperRad waveform model used to evolve boson clouds and extract spin-down and gravitational-wave observables.","marker":"[57]"},{"why":"Provides the updated SuperRad treatment of relativistic vector and scalar instabilities used for the spin-down and signal calculations.","marker":"[58]"},{"why":"Defines the semi-analytical massive-black-hole population models (PopIII and Q3 variants) that supply merger rates and binary properties.","marker":"[49]"},{"why":"Provides the updated population catalogs, tuned to pulsar-timing-array observations, that set the detectable LISA merger populations.","marker":"[50]"},{"why":"Establishes the method for computing the maximum allowed black-hole spin $a_{\\rm max}$ as a function of boson mass and spin-down timescale.","marker":"[24]"},{"why":"Earlier LISA scalar-boson constraint forecast whose mass range the paper compares against its spin-measurement results.","marker":"[18]"},{"why":"Earlier LISA scalar-boson constraint forecast used as a baseline for the spin-measurement exclusion ranges.","marker":"[64]"},{"why":"Supplies lisabeta, used to compute sky- and orientation-averaged signal-to-noise ratios for the merging binaries.","marker":"[120]"},{"why":"Provides the IMRPhenomHM waveform model used to model the merger gravitational-wave signals for detectability.","marker":"[121]"}],"fun_headline_variants":["LISA's twin probes shrink ultralight boson search space","Black hole spins seen by LISA could rule out ultralight bosons","LISA spin data alone can exclude ultralight bosons across decades","Remnant follow-ups add boson detection promise for LISA"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The forecasts assume that accretion re-spins massive black holes on a Salpeter timescale of about $4.5\\times10^7$ years, so any superradiant spin-down slower than that leaves no observable signature; if real accretion histories are faster, episodic, or less coherent, the excluded mass ranges shrink or shift.","fun_headline_variants_meta":{"raw":{"variants":["LISA's twin probes shrink ultralight boson search space","Black hole spins seen by LISA could rule out ultralight bosons","LISA spin data alone can exclude ultralight bosons across decades","Remnant follow-ups add boson detection promise for LISA"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000264,"raw_usage":{"total_tokens":1723,"prompt_tokens":1186,"completion_tokens":537,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":802,"completion_tokens_details":{"reasoning_tokens":460}},"tokens_in":802,"tokens_out":537,"duration_ms":5441,"temperature":1.0,"reasoning_tokens":460,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T10:20:55.458867+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take any LISA-detected massive binary whose confidently measured primary spin exceeds $a_{\\rm max}$ for a boson mass inside the claimed exclusion window under the fiducial 10% spin uncertainty and $T_S$; one such event would already break the exclusion for that mass. Conversely, if a vector boson near $10^{-16}$ eV exists and the Q3nodelays population is accurate, a four-year mission with no post-merger cloud signal at SNR $\\geq 10$ would contradict the paper's near-unity detection probability.","supporting_citations":[{"cited_title":"Baryakhtar, R","cited_arxiv_id":null,"evidence_quote":"Establishes the method for computing the maximum allowed black-hole spin $a_{\\rm max}$ as a function of boson mass and spin-down timescale."}],"review_version":1}