{"id":"6f19de32-2013-47c2-94c0-6cb711918a8e","arxiv_id":"2509.10420","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Simulated low-mass-gap boson star binaries are severely mis-measured when analyzed with black hole waveform models, but a model including both spin-induced quadrupole and tidal effects recovers the correct masses.","lead":"This paper uses simulated gravitational wave signals to test whether observations can tell black holes from other compact objects in the 3 to 5 solar mass gap. It finds that analyzing a boson-star-like binary with an ordinary black hole model misidentifies a 4+4 solar mass pair as an 8+2 solar mass pair, so waveform model choices directly affect how such detections are interpreted.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The headline (4,4)-to-(8,2) misidentification is demonstrated only within TaylorF2 with a Kerr-ISCO cutoff; the paper's own caveats leave open that a real boson-star waveform would not produce this specific bias.","rationale":"The reader's weakest assumption is exactly the point I find most load-bearing: the entire demonstration relies on TaylorF2 for both injection and recovery, with a Kerr-ISCO cutoff that ignores boson-star matter effects. I agree because the paper's headline conclusion, that an incorrect model biases a (4,4) boson-star-like signal to (8,2) M_sun and that joint SIQM+tidal modeling fixes this, is only as strong as the faithfulness of the TaylorF2 approximation to a real spinning boson star binary. The authors themselves flag in Sec. IV that the ISCO calculation does not account for SIQM or tidal effects and that no complete inspiral-merger-ringdown waveforms exist for such objects. The successful recovery with the joint model is expected by construction, since injection and recovery use the same phase model, so it does not independently validate the recovery or the bias. This does not invalidate the paper; the model-mismatch study is a legitimate and useful demonstration within its stated approximations, and the paper is appropriately cautious about them. However, the specific quantitative claim, including the (8,2) M_sun misidentification and the assertion that these biases are fundamental and not cured by increased SNR, should be tested with at least a different injection cutoff or an independent waveform before being applied to real observations. That is precisely the condition the reader attached, so the CONDITIONAL verdict remains appropriate and I recommend no change.","tokens_in":20828,"tokens_out":4110,"duration_ms":39680,"concrete_test":"Rerun the Sec. III B (4,4) boson-star injection with the same TaylorF2 phase model but truncate the injected signal at the boson-star's own innermost circular orbit (derived from the self-interacting boson star relations of Pacilio et al. 2020, which the paper already cites for the SIQM and tidal values) instead of the Kerr ISCO, while keeping the BBH recovery truncation at the Kerr ISCO. If the recovered component masses no longer peak near (8,2) M_sun, the headline misidentification is an artifact of the Kerr-ISCO truncation and the paper's central claim needs to be qualified. A complementary check is to repeat the recovery using an independently constructed boson-star waveform (for example, one that adds spin-tidal or higher-PN terms to the TaylorF2 phase) as the injection.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative result is computed entirely within one waveform family: both injection and recovery use TaylorF2 inspiral-only phasing, with the SIQM and tidal terms entering identically, and the likelihood is truncated at the Kerr ISCO (Sec. II C, Sec. IV). The paper explicitly states that this ISCO calculation does not take into account the SIQM or tidal deformability effects and that no complete inspiral-merger-ringdown models exist for boson stars. Therefore, the bias that a (4,4) boson-star-like binary is recovered as (8,2) M_sun under a BBH model has only been shown for this approximate model against itself. The joint SIQM+tidal recovery finding the true masses is a self-consistency check, not an independent validation of the waveform's faithfulness. If a real spinning boson star binary terminates earlier (its own ISCO), includes spin-tidal or dynamical-tide terms, or emits a merger signal, the recovered masses under the BBH model could shift, shrink, or even disappear. Because the paper's main claim is that model error, not statistics, dominates parameter recovery for such signals, this external-faithfulness assumption is load-bearing.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a Bayesian parameter-estimation study of simulated compact-binary signals in the lower mass gap, using TaylorF2 inspiral-only waveforms augmented with spin-induced quadrupole moment (SIQM) and tidal-deformability terms. The analysis has three parts: (i) BBH injections recovered with BBH, SIQM, and tidal models to quantify null-test performance; (ii) boson-star-like injections with both SIQM and tidal effects recovered with BBH, SIQM-only, tides-only, and SIQM+tides models; and (iii) injections with each non-BH effect separately. The central result is that a simulated equal-mass (4,4) solar-mass boson-star-like binary recovered with a BBH model yields biased masses around (8,2) solar masses, which the authors argue would likely be interpreted as a neutron-star-black-hole binary, whereas the joint SIQM+tides recovery model recovers the true masses. The paper concludes that model error, not statistics, dominates parameter recovery for such signals and that joint modeling is required for reliable distinguishability tests in the low-mass gap.","tokens_in":21008,"tokens_out":6037,"duration_ms":46880,"significance":"If the result holds, it provides a concrete and quantitative warning for tests of black-hole nature in the lower mass gap: measuring SIQM and tidal parameters separately, or ignoring them entirely, can produce severe biases in inferred masses and spins. The paper's strengths are its use of full Bayesian inference with Bilby and dynesty, a systematic injection-recovery matrix, and explicit exploration of spin and SNR dependence. The headline (8,2) solar-mass misidentification is striking and easy to falsify or confirm with future tests. However, the demonstration is entirely internal to the TaylorF2 waveform family with a Kerr-ISCO cutoff; the paper's own Sec. IV caveats (no complete inspiral-merger-ringdown models for boson stars and an ISCO calculation that ignores the very effects being tested) mean that external waveform faithfulness remains an open question. The result is therefore best interpreted as a model-consistency demonstration until robustness to the cutoff and waveform family is shown.","major_comments":[{"comment":"The central quantitative claim—that a (4,4) solar-mass boson-star-like signal is recovered as (8,2) solar masses under a BBH model—is computed with the same TaylorF2 inspiral-only waveform family used for both injection and recovery, with the likelihood truncated at the Kerr ISCO frequency. The paper itself states in Sec. IV that this ISCO calculation \"does not take into account the SIQM or tidal deformability effects\" and that no complete inspiral-merger-ringdown models exist for boson stars. Consequently, the specific bias has been demonstrated only for this approximate model against itself; a real spinning boson-star binary with a different terminal frequency or additional matter effects could produce a different or absent bias. The authors should either (a) test the robustness of the bias to the choice of cutoff (e.g., truncating at a boson-star ISCO or at a frequency where the TaylorF2 phase terms saturate) and to an alternative waveform approximant, or (b) explicitly restrict the claims in the abstract and Sec. III B to the inspiral-only TaylorF2 model.","section":"Secs. II C, III B, and IV"},{"comment":"The statement that the binary \"will be identified\" as having masses (8,2) solar masses is made without quantitative uncertainty. Please report the posterior median and 90% credible interval for the component masses, and ideally the joint m1-m2 posterior, for the (4,4) injection recovered with the BBH model for the case shown in Fig. 5 (spins 0.6, 0.5, SNR 100). This is needed to judge whether (8,2) is a robust posterior mode rather than a local peak or an artifact of a strongly asymmetric tail, and to support the neutron-star--black-hole misidentification claim.","section":"Sec. III B (\"Effect of SNR\") and Fig. 5"},{"comment":"The injected and recovered waveforms use identical SIQM and tidal phase parametrizations (Eqs. 4 and 5). The (SIQM tides)inj:(SIQM tides)rec recovery therefore validates internal consistency rather than waveform faithfulness. This should be stated explicitly in the discussion; otherwise readers may over-interpret the successful recovery as evidence that the TaylorF2 boson-star model is an adequate representation of the true physical signal. The distinction matters because the paper's conclusion about the importance of model accuracy is only as strong as the external validity of the TaylorF2 parametrization.","section":"Secs. II A and III B"}],"minor_comments":[{"comment":"The color labeling is inconsistent between the body text and the figure caption: the text describes the correct-model recovery as \"red histograms (SIQM tides)inj : (SIQM tides)rec,\" while the caption assigns red to \"(SIQM tides)inj−(BBH) rec\" and blue to the joint model. Please reconcile the color code across text, caption, and legend.","section":"Sec. III B and Fig. 3"},{"comment":"The line-style description is contradictory: it first says \"SNR=50 (dashed line) and SNR=100 (solid line),\" then says \"the posteriors are broader with SNR=50 (solid line) and are more biased compared to SNR=100 (dashed line).\" Clarify which line style corresponds to which SNR.","section":"Fig. 6 caption"},{"comment":"The prior ranges used for the SIQM parameters delta_kappa_i and the tidal deformability parameters Lambda_i are not stated. Please specify these priors, as they directly affect the Bayes factors reported in Fig. 2 and the interpretation of the posteriors in Figs. 4 and 8.","section":"Sec. II B"},{"comment":"There are several typos and formatting errors: \"Baye's theorem\" should be \"Bayes' theorem\"; \"signal-to-ratio\" should be \"signal-to-noise ratio\"; \"chie f f\" appears in Sec. III C; \"paramete\" appears in the Fig. 2 caption; and the reference list contains nonstandard entries (e.g., Ref. [100]). These should be corrected in a final pass.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is internally consistent and the simulation framework is sound. My main concern is that the headline bias claim is being presented as a forecast for real observations without the waveform-faithfulness caveats appearing in the abstract. The authors' own Sec. IV admits that the ISCO truncation neglects the very effects being tested, and the joint recovery result is a self-consistency check rather than an external validation. A robustness test that varies the ISCO cutoff, or a clear rephrasing of the abstract to say \"within the TaylorF2 inspiral-only model,\" would substantially increase my confidence. If the authors add such a test and soften the abstract, the paper would be suitable for publication; without that, the central claim is overstated relative to the evidence presented."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. The paper runs clean full Bayesian injection-recovery tests for low-mass-gap compact binaries: a simulated (4,4) M_sun boson-star-like signal, built with TaylorF2 including spin-induced quadrupole moments (SIQM) and tidal deformability, is recovered as (8,2) M_sun when analyzed with a BBH model, while a recovery model that includes both effects returns the injected masses. That is a useful, concrete warning that single-effect tests (SIQM only or tidal only) are not enough in this mass range. The second thing: the headline numbers come from using the same inspiral-only waveform family for both injection and recovery, with the likelihood truncated at the Kerr ISCO, so the specific (4,4)->(8,2) misidentification is a property of TaylorF2 against itself, not yet a robust prediction for real boson-star signals.\n\nWhat is genuinely new: the framework of joint SIQM+tidal measurement already existed (Narikawa et al. 2021), and boson-star SIQM/tidal relations were already computed (Pacilio et al. 2020; Vaglio et al. 2022, 2023). The new quantitative content is the low-mass-gap bias study itself: the mass misidentification example, the comparison of single-effect versus joint recovery, and the result that increasing SNR from 50 to 100 does not cure the bias. Those are real additions.\n\nWhat the paper does well: the simulations are internally consistent. Full Bayesian inference with Bilby and dynesty, several spin configurations, two mass ratios, SNRs of 50 and 100, and posterior distributions are shown. The paper is also straightforward about its main caveats in Section IV: inspiral-only waveform, no complete IMR models for boson stars, and a Kerr ISCO cutoff that ignores SIQM and tidal effects.\n\nSoft spots, in proportion. The main one is the external faithfulness assumption. Because injection and recovery share the same TaylorF2 phase model, the correct-model-recovers-truth result is a self-consistency check, not a validation against a real boson-star signal. If a real spinning boson star has a different ISCO or emits merger radiation, the specific bias could shift or disappear. The paper acknowledges this, but the abstract and Section III B present the (8,2) result more firmly than the caveat warrants. The phrase \"first quantitative estimates\" in Section IV is also too strong, given the prior joint-measurement framework and boson-star relations in the cited literature. And there is no code release and no noise-realization study, so the numbers are not independently reproducible. Minor figure-caption inconsistencies add small friction but do not change the conclusions.\n\nOverall: the qualitative conclusion—model error, not statistics, can dominate parameter recovery in the low-mass gap, and joint SIQM+tidal modeling matters—holds within the model family studied. The paper deserves a serious referee. I would send it to peer review with the expectation of revision: release code, temper the priority claim, and make the waveform-faithfulness caveat more prominent in the abstract. This is the kind of paper that is useful for people interpreting events like GW230529 and planning O5 and third-generation tests.","headline":"Sound simulation-based warning about model-induced biases in low-mass-gap tests, but the headline numbers depend entirely on one inspiral-only waveform family and some priority claims are overstated.","tokens_in":21637,"tokens_out":3264,"would_cite":true,"duration_ms":27019,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Model error, not detector statistics, dominates parameter recovery for non-black-hole low-mass-gap binaries; joint spin-quadrupole and tidal modeling is needed for reliable identification.","keywords":["gravitational waves","low mass gap","boson stars","spin-induced quadrupole moments","tidal deformability","waveform systematics","parameter estimation","tests of black hole nature"],"falsifier":"Re-run the $(4,4)\\,M_\\odot$ equal-mass boson-star injection and recovery using a waveform that includes merger and ringdown, or that truncates at a boson-star-specific innermost circular orbit instead of the black-hole one. If the black-hole-binary recovery no longer returns a posterior near $(8,2)\\,M_\\odot$, or the joint SIQM+tidal recovery no longer returns $(4,4)\\,M_\\odot$, the paper’s central claim—that model error dominates and is cured by joint modeling—would be called into question. A complementary observational check is to take a future high-SNR low-mass-gap event and compare its black-hole-template posterior with its joint SIQM+tidal posterior; an $(8,2)$-like shift would confirm the bias mechanism.","tokens_in":20570,"feed_emoji":"🌊","tokens_out":14040,"duration_ms":104405,"temperature":0.7,"pith_summary":"Gravitational-wave observatories are beginning to see compact objects in the ~3–5 solar mass “low mass gap,” a regime where the same mass could be a black hole or an exotic star such as a boson star. This paper asks whether the two standard fingerprints of non-black-hole matter—spin-induced quadrupole moments (SIQM) and tidal deformability—can tell them apart, and what happens when the analysis model is wrong. Using simulated TaylorF2 inspiral signals and full Bayesian inference, the authors show that a boson-star-like binary with component masses $(4,4)\\,M_\\odot$, when analyzed with a black-hole-binary template, is recovered as $(8,2)\\,M_\\odot$ and would likely be classified as a neutron star–black hole binary. Only a recovery model that includes both SIQM and tidal effects returns the true masses; models with just one effect still leave biases. The paper’s point is that waveform-model accuracy, not detector sensitivity or signal strength, is the limiting factor for these distinguishability tests, and the biases are not cured by doubling the signal-to-noise ratio.","feed_headline":"4,4 becomes 8,2 under black-hole templates","feed_subtitle":"A simulated boson-star binary analyzed with black-hole templates is misidentified as a neutron star-black hole merger.","key_machinery":"The load-bearing object is the TaylorF2 inspiral-only post-Newtonian waveform, whose phase is written as $\\psi(f)=2\\pi f t_c+\\phi_c+\\frac{3}{128\\eta v^5}\\left(1+\\psi_{\\mathrm{PP}}+\\psi_{\\mathrm{SIQM}}+\\psi_{\\mathrm{tidal}}\\right)$. Spin-induced quadrupole moments enter through $\\kappa_i=1+\\delta\\kappa_i$ at 2PN and 3PN order, and tidal deformability enters through $\\tilde{\\Lambda}$ at 5PN and 6PN order. The same TaylorF2 family is used for both injection and recovery, and the likelihood is truncated at the innermost stable circular orbit frequency of a black hole. Because injection and recovery share the same phase model, any displacement between the true and recovered masses in the results is attributed to which physical effects are omitted from the recovery template, isolating model-error bias from waveform-model uncertainty.","core_discovery":"The paper’s central claim is that using an incomplete gravitational-wave model to analyze a non-black-hole low-mass-gap binary produces parameter biases large enough to misidentify the source’s nature, and that this bias is avoided only when the recovery model includes both spin-induced quadrupole moments and tidal deformability together. For a simulated equal-mass binary boson star with $\\delta\\kappa_1=\\delta\\kappa_2=10$ and $\\Lambda_1=\\Lambda_2=289$, the black-hole-binary recovery gives a posterior concentrated near $(8,2)\\,M_\\odot$—interpreted as a neutron star–black hole binary—while the SIQM-only, tides-only, and combined recoveries give approximately $(4.7,3.4)\\,M_\\odot$, $(4.4,3.5)\\,M_\\odot$, and the true $(4,4)\\,M_\\odot$ respectively. The authors also find that genuine black-hole-binary injections recovered with SIQM or tidal models show no significant bias and Bayes factors favor the black-hole hypothesis, so the misidentification is specific to analyzing non-black-hole signals with incomplete templates.","pith_inferences":["I infer that if such biases are present in real searches, a population of low-mass-gap non-BH binaries could be cataloged as neutron star–black hole binaries, changing inferred merger rates and formation-channel demographics; the paper does not quantify this population-level consequence.","The $(4,4)\\,M_\\odot \\to (8,2)\\,M_\\odot$ shift suggests that a single effective parameter (a combination of mass ratio, spin, and the omitted SIQM/tidal terms) absorbs the signal; a follow-up could build a diagnostic that flags when recovered parameters are driven by this degeneracy rather than by astrophysics.","Because the paper uses the same TaylorF2 family for injection and recovery, I would treat the specific $(8,2)$ numbers as a property of this model family; repeating with an independent inspiral-merger-ringdown waveform would test whether the misidentification is generic or specific to the inspiral-only approximation.","The persistence of the biases with higher SNR is established only within the inspiral-only model; with merger/ringdown and higher harmonics the degeneracies could break differently, which the paper itself flags as future work."],"forward_implications":["If an observed low-mass-gap event is a non-black-hole binary, standard black-hole-binary analyses will not just be uncertain; they can actively return the wrong masses and spins, potentially classifying the source as a neutron star–black hole merger.","Tidal-deformability tests are stronger than SIQM tests for confirming that a signal is a black hole, whereas SIQM tests are more sensitive to spin magnitudes; both effects need to be measured jointly to avoid mass-ratio bias.","Combining both effects is feasible: the paper recovers non-BH values of the spin-quadrupole parameter (ruling out the black-hole value) even in the joint parameter space, and using symmetric combinations reduces bias compared with estimating the individual parameters.","Doubling the signal-to-noise ratio from 50 to 100 does not remove the biases, so these systematic errors can dominate statistical uncertainties even in the era of next-generation detectors."],"supporting_citations":[{"why":"Supplies the self-interacting boson star parameters (SIQM deviation and tidal deformability) used for the non-BH injections.","marker":"[38]"},{"why":"Defines the spin-induced quadrupole moment test that flags departures from the black-hole value.","marker":"[21]"},{"why":"Computes the post-Newtonian spin-dependent phase terms in the TaylorF2 waveform used throughout.","marker":"[89]"},{"why":"Provides the 2PN and 3PN spin-induced quadrupole phase corrections added to the point-particle phase.","marker":"[91]"},{"why":"Establishes the tidal deformability formalism whose phase correction enters the waveform at leading order.","marker":"[41]"},{"why":"Shows a Bayesian framework for measuring SIQM and tidal effects simultaneously, which this work adapts to the low-mass gap.","marker":"[51]"},{"why":"Supplies the public waveform implementation containing both SIQM and tidal effects.","marker":"[100]"},{"why":"Supplies the Bayesian inference library used for posterior computation.","marker":"[101]"},{"why":"Supplies the nested-sampling algorithm used to compute posteriors and Bayes factors.","marker":"[103]"}],"fun_headline_variants":["Boson star binary misread as 8,2 by black-hole templates","Black-hole template turns 4,4 boson star into 8,2","Gravitational wave templates mislabel boson stars as NS-BH","Incomplete waveform model skews compact object ID in low mass gap","Waveform model choice changes inferred masses of exotic binaries"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The demonstration assumes that the TaylorF2 inspiral-only waveform faithfully represents both the injected boson-star signal and the recovery templates, with the SIQM and tidal phase terms entering identically, and cuts the analysis at the black-hole innermost stable circular orbit frequency without accounting for the boson star’s own matter effects; if a real spinning boson-star waveform differs in the late inspiral or merger, the measured biases may not transfer to actual observations.","fun_headline_variants_meta":{"raw":{"variants":["Boson star binary misread as 8,2 by black-hole templates","Black-hole template turns 4,4 boson star into 8,2","Gravitational wave templates mislabel boson stars as NS-BH","Incomplete waveform model skews compact object ID in low mass gap","Waveform model choice changes inferred masses of exotic binaries"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000908,"raw_usage":{"total_tokens":3993,"prompt_tokens":1127,"completion_tokens":2866,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":743,"completion_tokens_details":{"reasoning_tokens":2772}},"tokens_in":743,"tokens_out":2866,"duration_ms":16735,"temperature":1.0,"reasoning_tokens":2772,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T15:54:34.899641+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run the $(4,4)\\,M_\\odot$ equal-mass boson-star injection and recovery using a waveform that includes merger and ringdown, or that truncates at a boson-star-specific innermost circular orbit instead of the black-hole one. If the black-hole-binary recovery no longer returns a posterior near $(8,2)\\,M_\\odot$, or the joint SIQM+tidal recovery no longer returns $(4,4)\\,M_\\odot$, the paper’s central claim—that model error dominates and is cured by joint modeling—would be called into question. A complementary observational check is to take a future high-SNR low-mass-gap event and compare its black-hole-template posterior with its joint SIQM+tidal posterior; an $(8,2)$-like shift would confirm the bias mechanism.","supporting_citations":[],"review_version":1}