{"id":"86685c4b-515a-4c44-8068-3ec590f77357","arxiv_id":"2608.02736","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"First-order scalarization of a neutron star triggered by accretion emits a long-lived inverse-chirp gravitational breathing mode that is orders of magnitude louder than the continuous second-order case.","lead":"This paper shows that neutron stars can switch to a strongly 'scalarized' state through a sudden, first-order phase transition, producing a loud gravitational wave signal. A generalist might read it because it turns a previously hard-to-see effect of modified gravity into a concrete, possibly detectable chirp-like signal.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Existing merger constraints on the scalar mass may exclude the paper's mφ=10^-14 eV example; the one-sentence dismissal in the Discussion is unsupported and directly threatens the observability claim.","rationale":"The reader's weakest assumption is the hand-picked quartic coupling γ=12β^2. That is a real concern: the first-order transition at high mass requires a positive γ of order β^2, and the paper's 'naturalness' argument (cancellation of the φ^4 term in A^4) is one convention, not a derived relation. However, the paper explicitly tests a second value (γ=18β^2), presents γ=12β^2 as a zero of a specific coefficient, and its central claim is conditional on first-order scalarization existing. The more decisive parameter for the headline observability claim is the scalar mass mφ: the inverse-chirp frequency scales linearly with mφ, and the chosen value 10^-14 eV sits three orders of magnitude below a tentative bound that the paper dismisses without calculation. If that bound applies, the signal is either outside the detector band (early times) or at frequencies where the source burst has negligible spectral power, so the computed loud signal would not be observable. This is a distinct, and in my view more load-bearing, weakness than the γ choice, because it connects the central example to existing data and the paper's rebuttal is a single unsubstantiated sentence. I therefore partially agree with the reader: both are parameter choices, but the mφ/constraint issue is the one whose failure most directly collapses the claimed observability. The verdict remains CONDITIONAL because the concern is addressable by a definite numerical check and does not invalidate the mechanism itself.","tokens_in":11840,"tokens_out":27072,"duration_ms":242163,"concrete_test":"Compute quasi-equilibrium sequences of a 1.35+1.35 M_sun neutron-star binary (same HB EOS) with the paper's parameters (mφ=10^-14 eV, α0=10^-2, β=-5, γ=12β^2), using the solver of Kuan et al. [49], down to separations corresponding to orbital frequencies of 10-1000 Hz. If the stars acquire a large scalar charge before merger and the resulting dephasing relative to GR exceeds the GW170817 measurement uncertainty, the mφ=10^-14 eV example is excluded and the predicted signal is not astrophysically observable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The detectability of the predicted inverse-chirp breathing mode is controlled by the scalar mass mφ=10^-14 eV. Equation (6) places the signal at ~310 Hz one year after the transition and ~10 Hz after 1000 years, inside the aLIGO/Virgo/KAGRA band. The paper itself notes that mφ≳10^-11 eV would move the signal out of the detector band and cites reference [49] for a tentative bound mφ≳10^-11 eV for α0=0=γ. Its entire rebuttal is one sentence: 'these limits are qualitative and do not yet consider the first-order scalarization that we study.' No analysis is offered. First-order scalarization could plausibly strengthen, not weaken, merger constraints: during late inspiral, dynamical scalarization would turn on abruptly and emit a scalar-radiation burst whose absence in GW170817 would be a direct constraint on β, α0, and γ at this mφ. If the mφ bound survives, the 1-year inverse-chirp frequency becomes ~300 kHz, where a ~ms source burst has negligible Fourier power, so the computed signal would not be observable with current detectors. The central claim therefore rests on an unexamined assumption that the existing constraint is evaded by first-order dynamics.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper argues that spontaneous scalarization in scalar-tensor theories, traditionally modeled as a second-order phase transition, is often first-order once the quartic scalar-matter coupling γ is included. It claims that first-order scalarization leads to abrupt transitions that emit strong monopolar scalar radiation, which at large distances appears as a long-lived inverse-chirp gravitational-wave breathing mode. The authors demonstrate this with spherically symmetric simulations of a neutron star driven across the phase transition by accretion, using a scalar-tensor extension of the GR1D code, with parameters α0=1e-2, β=-5, mφ=1e-14 eV, and γ=12β^2 (or 18β^2). They show that the resulting strain spectral density exceeds the sensitivity of aLIGO, aVirgo, and KAGRA at 10 kpc for a 60-day coherent observation, whereas the γ=0 (second-order) case yields no observable signal. They also discuss generalizations to black holes and binaries.","tokens_in":12029,"tokens_out":8854,"duration_ms":79188,"significance":"If correct, the paper establishes an important qualitative point: first-order scalarization produces qualitatively different, much louder gravitational-wave signatures than the continuously studied second-order case. The numerical pipeline is credible, with a demonstrated third-order convergence of the scalar constraint and a 1000-fold variation of the accretion amplitude leaving the signal unchanged; wave extraction and propagation follow published methods. The predicted inverse-chirp breathing mode is a falsifiable, detector-relevant prediction, and the general claim is likely to stimulate further work on first-order scalarization in binaries and black-hole systems. However, the specific observability claim rests on a scalar mass mφ=1e-14 eV that is in tension with existing bounds, and the paper does not quantitatively address this tension, leaving the central claim only conditionally supported.","major_comments":[{"comment":"The observability of the inverse-chirp signal is controlled by mφ, and the paper itself states that mφ ≳ 1e-11 eV would move the signal out of the detector band. The only response to the existing tentative bound mφ ≳ 1e-11 eV (Ref. [49]) is one sentence (\"these limits are qualitative and do not yet consider the first-order scalarization that we study\"). This is not adequate support for a central observability claim. The first-order scenario could plausibly strengthen, not weaken, merger constraints: a first-order dynamical scalarization during inspiral would produce a burst of scalar radiation whose absence in GW170817 would directly constrain the theory. The authors should either provide a quantitative estimate of the scalar radiation in a binary inspiral with first-order scalarization and compare with current bounds, or identify a region of parameter space (e.g., with mφ smaller than the bound) that is consistent with observations and still yields the predicted signal. Without this, the conclusion that the signal is observable is not established.","section":"Results and Discussion, around Eq. (6)"},{"comment":"The statement that \"a large part of the parameter space leads to the loud signals\" is not supported by the presented analysis. The paper varies only γ among three values (0, 12β^2, 18β^2), and qualitatively discusses β, mφ, and α0, but does not scan the parameter space. In particular, the strain is proportional to α0 (Eq. (5)), and the paper notes that mφ changes the transition and may push it to second order; both effects can significantly shrink the observable region. A systematic scan (e.g., a contour of signal-to-noise ratio in the β–γ or α0–mφ plane) would justify the \"large part\" claim, or the claim should be softened to the specific, representative region studied.","section":"Results, parameter-dependence paragraph and Figure 3"}],"minor_comments":[{"comment":"In the first paragraph of Results, \"α = 10^{-2}\" should read \"α0 = 10^{-2}\" to match the notation introduced in Eq. (3).","section":"Results"},{"comment":"Reference [36] contains a typo: \"later discover\" should be \"later discovered\".","section":"Reference [36]"},{"comment":"In the End Matter, \"we take the fourier transform\" should be capitalized as \"Fourier transform\".","section":"End Matter, Wave Propagation"},{"comment":"Figure 3 would be clearer if the caption stated the assumed source distance (D=10 kpc) and coherent observation time (T=60 days), which currently appear only in the text.","section":"Figure 3"},{"comment":"The approximation h_B ≈ 2α0φ is presented without stating that it is valid in the far zone where φ is small; a brief clarification would avoid confusion because near the star the scalar field can be of order unity.","section":"Eq. (5)"}],"recommendation":"major_revision","confidential_remarks":"The mφ constraint issue is the main stumbling block. The authors should not be allowed to dismiss a constraint that directly controls the central prediction without quantitative analysis. Also, the parameter choices (α0=1e-2, β=-5) are aggressive, and the paper should at least cite or briefly discuss existing constraints for massive scalars. Otherwise, the paper's core claim is conditional on an unexamined assumption. The numerical work itself appears solid, and the paper is well written."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I read this paper carefully. The core result is real: when scalarization is a first-order transition, an accreting neutron star jumps discontinuously to a strongly scalarized state, and the massive scalar field generates a long-lived inverse-chirp breathing mode. That effect is new, and the numerical implementation looks solid. The simulation shows third-order convergence, the accretion trigger amplitude can be varied by a factor of 1000 without changing the signal, and the wave extraction follows published methods. The contrast between the loud first-order signal and the silent second-order case is clean. The waveform is a genuine dynamical output, not a fit. That's real evidence.\n\nThe soft spots are where the claims outrun the evidence. First, gamma = 12 beta^2 is presented as \"natural\" because it cancels the quartic term in the expansion of A^4. That's a reasonable heuristic, but it's still a free parameter. If gamma were near zero, the transition is second-order and the signal vanishes. They do show gamma = 18 beta^2 makes the signal louder, so within their chosen family the conclusion is robust, but it does not cover the theory space.\n\nSecond, and more worrying: the paper uses m_phi = 10^-14 eV and dismisses the tentative bound m_phi ≳ 10^-11 eV from binary merger studies in a single sentence. That dismissal does not hold up. During late inspiral, first-order scalarization would turn on abruptly and emit a scalar burst; the absence of such a burst in GW170817 is likely a direct constraint on beta, alpha0, and gamma at this mass, and it is not obvious the constraint weakens. If the bound survives, the signal frequency after one year becomes ~300 kHz, far outside the detector band. So the central observability result is threatened by an unexamined assumption. This is load-bearing, not a minor detail.\n\nMinor issues: no code or data released, and uncertainty quantification is thin. But those are addressable.\n\nThis paper is a serious contribution to the scalarization literature and deserves refereeing, but the current version is conditional. The authors need to show that existing bounds on the scalar mass do not already exclude their example. If that can be done, it is an important result.\n\nRecommendation: send to peer review.","headline":"A credible numerical study of a genuinely new effect—first-order scalarization producing an inverse-chirp breathing mode—but the observability claim leans on an untested parameter choice and an under-examined scalar-mass bound.","tokens_in":12626,"tokens_out":3203,"would_cite":true,"duration_ms":28052,"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":"When scalarization happens as a first-order phase transition, an accreting neutron star jumps discontinuously into a strongly scalarized state and emits a gravitational-wave breathing mode that current detectors could see at 10…","keywords":["spontaneous scalarization","first-order phase transition","scalar-tensor gravity","neutron stars","gravitational waves","breathing modes","massive scalar field","inverse chirp"],"falsifier":"Search existing LIGO/Virgo/KAGRA data for the predicted breathing-mode inverse chirp: a signal whose instantaneous frequency evolves as $f(t)\\approx f_*/\\sqrt{2(t-D)/D}$ with $f_*=2.42\\,\\mathrm{Hz}$ and $D\\approx10$ kpc, appearing within weeks after an accretion event on a nearby neutron star. A null result over a sizable sample, at the strain levels the paper computes for $\\gamma=12\\beta^2$, would rule out first-order scalarization for those parameters; a detection would confirm it.","tokens_in":11598,"feed_emoji":"🌊","tokens_out":11060,"duration_ms":87404,"temperature":0.7,"pith_summary":"Spontaneous scalarization is a phase transition in which a neutron star suddenly grows a macroscopic scalar-field cloud and departs from general relativity. This paper argues that in many astrophysically relevant scalar-tensor theories the transition is first-order rather than the smooth, second-order process usually assumed. It simulates an accreting neutron star pushed past the point of metastability: the star abruptly jumps to a strongly scalarized configuration, releases latent heat, and rings with coupled fluid and scalar oscillations. The resulting burst of massive-scalar radiation disperses into a long-lived inverse chirp, and its gravitational-wave breathing-mode strain is calculated to rise above the noise of current ground-based detectors for a source at 10 kiloparsecs. The same discontinuity, the paper argues, should create novel observable signatures in many other scalarization scenarios, turning first-order scalarization into a new test of deviations from general relativity.","feed_headline":"Abrupt scalarization could make neutron stars loud GW emitters","feed_subtitle":"A discontinuous transition emits a detectable inverse-chirp breathing mode; the smooth second-order version stays silent.","key_machinery":"The mechanism is a Landau-type energy expansion of the ADM mass, Eq. (4), in terms of a scalarization-strength measure $\\Phi$: $M_{\\rm ADM}=M_0+a\\Phi^2+\\tfrac12 b\\Phi^4+\\tfrac13 c\\Phi^6+\\cdots$. When the quartic coefficient $b$ is positive, scalarized minima emerge continuously (second-order transition); when $b<0$, weakly and strongly scalarized minima coexist behind an energy barrier and the globally preferred solution changes discontinuously (first-order). The sign of $b$ is controlled by the quartic matter–scalar coupling $\\gamma$ in $A(\\phi)$: positive $\\gamma$ lowers $b$ and favors first-order scalarization, and the paper adopts $\\gamma=12\\beta^2$ because that choice cancels the quartic term in the Taylor expansion of $A^4$, making it a natural representative point. The second ingredient is the massive scalar field, whose frequency-dependent group velocity converts the short source burst into a long-lived inverse chirp as it propagates to kiloparsec distances.","core_discovery":"The central discovery is the general fact that replacing second-order scalarization with its first-order version radically changes the observable outcome of scalarization events, because the transition is discontinuous. In the model studied—a massive scalar field with mass $m_\\phi = 10^{-14}\\,\\mathrm{eV}$ coupled to matter through $A(\\phi)=\\exp(\\alpha_0\\phi + \\tfrac12\\beta\\phi^2 - \\tfrac1{24}\\gamma\\phi^4)$ with $\\beta=-5$ and $\\gamma=12\\beta^2$—the equilibrium sequence of neutron stars has two locally stable branches separated by an energy barrier. When accretion raises the baryon mass past the local maximum, the star makes an abrupt, nonperturbative jump from a weakly scalarized to a strongly scalarized branch. This releases part of the binding-energy difference, excites stellar oscillations, and emits a monopolar scalar pulse. Because the scalar is massive, propagation disperses the pulse into an inverse chirp whose observed frequency slowly decreases on year timescales; the associated breathing-mode strain, scaled to 10 kpc and a 60-day observation, exceeds the noise amplitudes of aLIGO, aVirgo, and KAGRA. The same setup with $\\gamma=0$, where scalarization is second-order, produces no signal distinguishable from numerical error.","pith_inferences":["Because $\\gamma=12\\beta^2$ is a representative choice rather than a derived relation, the paper's observability prediction effectively turns the detection problem into a direct measurement of the quartic matter–scalar coupling: a loud inverse chirp from a known accreting neutron star would be evidence for a first-order transition, and its absence would push $\\gamma$ toward zero.","The inverse-chirp waveform's frequency evolution is fixed by the scalar mass and source distance alone (Eq. 6), so even without precise knowledge of the accretion trigger, a targeted narrow-band search in archival LIGO/Virgo/KAGRA data could test the prediction.","The same latent-heat release that powers the scalar burst may also heat the stellar material and produce electromagnetic or neutrino counterparts, which the paper does not compute but which would give independent confirmation if first-order scalarization actually occurs.","Eccentric binaries could cross the energy barrier repeatedly, producing a train of descalarization/scalarization bursts—a signature with no analogue in second-order scalarization and a potentially clean probe of the transition's hysteresis."],"forward_implications":["Accretion-induced first-order scalarization in massive scalar-tensor gravity produces gravitational-wave breathing modes whose amplitude spectral density exceeds the current detector noise curves at 10 kpc, giving a new channel to test deviations from general relativity.","A null search for such inverse-chirp signals in existing detector data would constrain the quartic coupling $\\gamma$ and thereby the order of scalarization, since $\\gamma=0$ predicts silence.","The same abrupt-transition logic extends to descalarization and to black-hole scalarization, including spin-induced transitions, where discontinuous jumps would create distinct observable bursts rather than smooth evolution.","In binary systems dynamical and induced scalarization would become abrupt at small separations, changing the expected waveform morphology relative to all second-order studies to date.","The signal frequency scales linearly with the scalar mass, so the effect is observable with current detectors only for $m_\\phi \\lesssim 10^{-11}\\,\\mathrm{eV}$; heavier scalar fields move the signal out of band."],"supporting_citations":[{"why":"Establishes that scalarization in the DEF model and simple extensions is generically first-order and supplies the Landau energy expansion used to classify transition order.","marker":"[13]"},{"why":"Shows that the quartic coupling gamma moves first-order scalarization to higher, astrophysically relevant stellar masses, motivating the gamma=12 beta^2 choice.","marker":"[21]"},{"why":"Provides the accretion prescription and the gravitational-wave signal from accretion-induced descalarization in massive scalar-tensor theory that this paper adapts to first-order scalarization.","marker":"[12]"},{"why":"Derives the long-lived inverse-chirp signal from massive-scalar dispersion and the method for propagating the scalar wave to kiloparsec distances.","marker":"[11]"},{"why":"Introduces the original Damour–Esposito-Farèse nonperturbative scalarization mechanism that the paper's coupling model extends.","marker":"[8]"},{"why":"Gives the massive scalar-tensor neutron-star structure and the flat-space wave propagation scheme used in the numerical evolutions.","marker":"[44]"},{"why":"Provides the aLIGO, aVirgo, and KAGRA noise amplitude spectral densities used to judge the observability of the predicted strain.","marker":"[46]"},{"why":"The modern review of spontaneous scalarization that frames the strong-field deviations from GR which first-order scalarization is claimed to make observable.","marker":"[7]"},{"why":"Supplies the flux-conservative hydrodynamic and scalar-tensor evolution equations on which the numerical code is based.","marker":"[62]"}],"fun_headline_variants":["First-order scalarization makes neutron stars chirp","Abrupt star transition emits detectable inverse chirp","Discontinuous scalarization turns neutron stars into beacons","Snap transition in neutron stars yields new gravity test","First-order phase shift in stars screams out signal"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The observability prediction rests on the theory having a positive quartic coupling $\\gamma$ of order $\\beta^2$; if the true coupling were zero or small, scalarization would be second-order and the predicted signal would vanish, and if it were smaller the transition would happen below about one solar mass where astrophysical neutron stars are rarer.","fun_headline_variants_meta":{"raw":{"variants":["First-order scalarization makes neutron stars chirp","Abrupt star transition emits detectable inverse chirp","Discontinuous scalarization turns neutron stars into beacons","Snap transition in neutron stars yields new gravity test","First-order phase shift in stars screams out signal"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000245,"raw_usage":{"total_tokens":1507,"prompt_tokens":888,"completion_tokens":619,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":504,"completion_tokens_details":{"reasoning_tokens":546}},"tokens_in":504,"tokens_out":619,"duration_ms":6477,"temperature":1.0,"reasoning_tokens":546,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T15:00:33.024185+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Search existing LIGO/Virgo/KAGRA data for the predicted breathing-mode inverse chirp: a signal whose instantaneous frequency evolves as $f(t)\\approx f_*/\\sqrt{2(t-D)/D}$ with $f_*=2.42\\,\\mathrm{Hz}$ and $D\\approx10$ kpc, appearing within weeks after an accretion event on a nearby neutron star. A null result over a sizable sample, at the strain levels the paper computes for $\\gamma=12\\beta^2$, would rule out first-order scalarization for those parameters; a detection would confirm it.","supporting_citations":[{"cited_title":"Evans, R","cited_arxiv_id":null,"evidence_quote":"Provides the aLIGO, aVirgo, and KAGRA noise amplitude spectral densities used to judge the observability of the predicted strain."},{"cited_title":"Rosca-Mead, U","cited_arxiv_id":null,"evidence_quote":"Supplies the flux-conservative hydrodynamic and scalar-tensor evolution equations on which the numerical code is based."}],"review_version":2}