{"id":"0db322c4-f2b4-4a8f-892b-addd50059279","arxiv_id":"2412.06620","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Tidal contributions to neutron-star binary gravitational-wave phase in scalar-tensor gravity involve scalar, tensor, and mixed Love numbers, and generally combine to reduce the net tidal signal below the general-relativistic value.","lead":"This paper computes how tidal deformation of neutron stars changes the gravitational-wave signal from binary inspirals in scalar-tensor theories of gravity, working to first post-Newtonian and first finite-size order. It finds that the extra scalar and mixed tidal effects enter with opposite signs, so the net tidal imprint is often smaller than in general relativity.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (83) has a sign error: the GR-like x^5 term should be −(39/(2α^5ξ^2))Λ̃, not +(39/(2α^5ξ^2))Λ̃, so the printed formula fails to recover the standard GR tidal phase.","rationale":"The reader's weakest assumption concerned dynamical friction, an environmental effect the paper itself flags. The sign error is a more load-bearing, internal inconsistency: Eq. (83) is the ready-to-use formula quoted in the abstract and summary, and it does not reduce to the correct GR limit, contradicting the paper's own Eq. (82a) and Fig. 9. Regardless of the DF debate, a reader using Eq. (83) would get the wrong sign for the dominant x^5 tidal term, and the claimed sign cancellation underlying 'net ST tidal signatures smaller than GR' is not supported by the printed formula. This warrants a conditional acceptance pending correction and rechecking of the derived conclusions, rather than the DF-based conditional verdict alone.","tokens_in":45947,"tokens_out":40243,"duration_ms":376402,"concrete_test":"Evaluate Eq. (83) in the GR limit by setting β0=0 (q=0, S−=0, c3=c5=0, α=ξ=1) and compare the resulting x^5 coefficient to the standard leading-order tidal phase from Flanagan & Hinderer [42] and to the paper's own Eq. (82a) computed with (82d) and (84) for equal masses. If the coefficient is +39/2 Λ̃ rather than −39/2 Λ̃, the sign error is confirmed. More directly, substitute the equal-mass GR values ζ2=2λT, ˜ζ2=2mλT, M=2m, η=1/4 into (82d) and verify ρ^{nd,2}_tid = −(39/2)Λ̃, then compare with the bracket in (83).","verdict_should_be":"CONDITIONAL","load_bearing_attack":"In the ready-to-use QD phase formula, Eq. (83), the coefficient of x^5 is printed as +(39/(2α^5ξ^2))Λ̃ + c5, with Λ̃ defined in (84) as a positive combination of the tensor quadrupolar deformabilities. This conflicts with the derivation in Sec. 4.5.2. In the GR limit, q=0, S−=0, α=ξ=1, c3=c5=0, Eq. (82a) gives ψ_non-dip,tid = 3/(128ηx^{5/2}) ρ^{nd,2}_tid x^5. Using (82d), for equal masses ρ^{nd,2}_tid = -24 ˜ζ2/(M^6η) - 216 ζ2/M^5 - 48 ζ2/M^5. With ζ2=2λT, ˜ζ2=2mλT, M=2m, η=1/4, this equals -39λT/(2m^5). Since Λ̃ = λT/m^5, ρ^{nd,2}_tid = -(39/2)Λ̃, so the x^5 coefficient in the factored form should be −(39/2)Λ̃, not +. The printed sign would give a positive GR tidal phase, contradicting the standard result and the paper's own Fig. 9, where the λT contribution is negative. The sign determines the direction of the cancellation claimed to make ST tidal effects smaller than GR, so the central claim as presented rests on the wrong sign in Eq. (83).","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes tidal contributions to the gravitational-wave phase of inspiraling neutron-star binaries in massless scalar-tensor theories, working to first order in a combined post-Newtonian and small-finite-size expansion. The authors derive a skeletonized action with three types of tidal Love numbers, compute scalar and tensor energy fluxes and the associated Fourier-domain phase in both dipole-driven and quadrupole-driven regimes, and provide ready-to-use expressions for the tidal phase coefficients. They then apply these results to Gaussian scalar-tensor couplings with the WFF1, SLy, and H4 equations of state, surveying parameter space and presenting case studies that lead to the claim that the net tidal GW phase signature in scalar-tensor theories is smaller than in GR because contributions from different Love numbers enter with opposite signs and different frequency scalings.","tokens_in":46255,"tokens_out":8483,"duration_ms":81748,"significance":"If correct, this work provides the first systematic analytic treatment of quadrupolar tidal effects in scalar-tensor inspirals, showing that all three Love numbers enter the GW phase with distinct frequency scalings and cannot be absorbed into a single GR-like deformability. The derivation is careful and largely self-contained: the double expansion in PN and tidal parameters is performed explicitly, the SPA phase is built from the energy balance, and the authors verify the GR limit and cross-check their numerical evaluation of beta(phi) against a finite-difference derivative at the 0.37% level. The physical insight that sign cancellations can reduce the net tidal phase, together with the clear parameter-space study, makes the paper a useful methodological contribution for template-building and theory-specific tests beyond GR.","major_comments":[{"comment":"The coefficient of x^5 in the ready-to-use phase formula is printed with a positive sign: (39/(2 α^5 ξ^2)) Λ̃ + c5. This is inconsistent with the derivation in Sec. 4.5.2. In the GR limit (q=0, S_−=0, α=ξ=1, c3=c5=0), using ζ2=2 λ_T, ζ̃2=2 m λ_T, M=2m, η=1/4 in Eq. (82d) gives ρ^{nd,2}_tid = - (39/2) Λ̃, so the x^5 term in Eq. (83) should be negative to reproduce the standard GR tidal phase ψ_tid = -(3/(128 η))(39/2) Λ̃ x^5. The printed plus sign yields a positive GR tidal phase, contradicting the standard result cited in Refs. [42,86] and the negative λ_T contribution shown in the authors' own Fig. 9. Because the sign of this dominant term drives the cancellation argument that leads to the abstract's central claim of smaller net ST tidal effects, Eq. (83) and its restatement in Sec. 6 must be corrected, and the consequences for the discussion of sign structure should be re-examined.","section":"Sec. 4.5.3, Eq. (83)"},{"comment":"The authors estimate F_DF/F_GW ~ O(10^-1) at 200 Hz for a representative equal-mass SLy system, yet the phase evolution in Sec. 4.5 accounts only for scalar and tensor radiation reaction. The ST-vs-GR tidal phase differences presented in Fig. 8 are at the 0.1 rad level, so an environmental force of the estimated relative size is expected to produce a phase contribution of comparable order. The manuscript appropriately defers environmental effects to future work, but the abstract's claim that net tidal GW imprints in ST gravity are smaller than in GR applies to the modeled two-body system, not to the complete physical system. The authors should add a quantitative estimate of the environmental phase or otherwise explicitly qualify the scope of the central claim.","section":"Sec. 2.2, Eq. (10)"},{"comment":"The scalar tidal expansion is performed under the assumption (L_ϕ/r) ≪ 1 with L_ϕ ~ 6 R_body, but this condition is not well satisfied in the frequency band where the numerical results are presented. For a 1.4+1.4 M_sun SLy system at 200 Hz, the orbital separation is only a few times L_ϕ, and the QD-phase results in Figs. 8 and 9 are shown up to 500 Hz or 1 kHz. The paper should either quantify the error incurred by truncating the scalar-cloud multipole expansion at these separations or restrict the claimed quantitative accuracy to frequencies where the expansion parameter remains small.","section":"Sec. 2.2, Eq. (8)"}],"minor_comments":[{"comment":"There are several typographical errors: \"dubble\" should be \"double\", \"Tabel 1\" should be \"Table 1\", and \"induvidual\" should be \"individual\". In Sec. 3 the term \"Newtionian\" should be \"Newtonian\".","section":"Sec. 5.5.3"},{"comment":"The phrase \"lnΛ denotes the Coulomb algorithm\" appears to be a typo; it should read \"Coulomb logarithm\".","section":"Sec. 2.2, Eq. (9)"},{"comment":"The header of Table 1 is garbled in the current text (\"contribution to largest zero sign magnitude for phase coefficients for for NS-NS\") and should be reformatted for clarity.","section":"Table 1"},{"comment":"The dipolar-driven phase formula in Appendix C.2 is labeled Eq. (D.1), which conflicts with the appendix labeling (the formula appears in Appendix C.2). Please renumber or relabel for consistency.","section":"Appendix C.2, Eq. (D.1)"}],"recommendation":"major_revision","confidential_remarks":"The sign error in Eq. (83) is the most consequential issue: although the numerical figures appear to be generated from the correct expressions in Sec. 4.5.2, the central ready-to-use formula and its repetition in the conclusions convey the wrong sign for the dominant tidal term. This is easily fixable by changing the sign in Eq. (83) and in the quoted form in Sec. 6, and then re-verifying the statements that rely on the sign structure. The dynamical-friction issue is a known omission that the authors acknowledge; a careful quantitative caveat would strengthen the paper. The work is otherwise systematic and timely, and I expect that after these revisions it could be published."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a serious, careful calculation, but Eq. (83) is wrong as printed. I checked the GR limit: with q=0, S-=0, alpha=xi=1, equal masses, Eq. (82a) with (82d) gives rho^{nd,2}_tid = -39/2 Lambda-tilde, so the phase is negative and matches the standard GR tidal phase. Eq. (83) instead has +(39/(2 alpha^5 xi^2)) Lambda-tilde, which would give a positive GR tidal phase, contradicting both the standard result and the paper's own Fig. 9. This is not a minor typo: (83) is the central deliverable, repeated in the summary, and anyone implementing it would fail to recover GR. The underlying derivation appears sound, so the error is likely in the transcription to the factored form, but it must be fixed.\n\nWhat is actually new: they give the first complete leading quadrupolar tidal phase for ST inspirals with all three Love numbers (scalar, tensor, mixed), including the frequency scalings that make net ST tidal effects smaller than GR. The derivation is systematic: skeletonized action, double PN/tidal expansion, multipole waveforms, energy balance, SPA. They recover the known dipolar results and the GR limit in the derivation (though not in the printed formula). The numerical case studies are thorough, and the beta(phi) computation is cross-checked at 0.37% precision.\n\nSofter spots, in proportion: (1) The sign error is the elephant in the room. (2) The dynamical friction estimate in Sec. 2.2 gives F_DF/F_GW ~ 0.1 at 200 Hz, same order as the tidal phase effects they compute (~0.1 rad). The authors flag it as future work, but it makes the quantitative predictions conditional. (3) They use beta0=-6, which they acknowledge is ruled out by pulsar tests, only for qualitative illustration. No code is shipped, but they extend a public code from their prior paper.\n\nBottom line: the paper deserves a serious referee. The derivation is valuable and likely correct, but the manuscript as written cannot be used as-is. I would send it to peer review with a strong request to fix the sign in (83), re-check the figures, and add a caveat about dynamical friction. It is for people building ST inspiral templates or doing beyond-GR tests; they should not copy (83) until revised.","headline":"The derivation is sound but the headline formula (Eq. 83) has a sign error in the GR limit: the Lambda-tilde term should be negative, which makes the paper's central ready-to-use phase unusable until corrected.","tokens_in":46820,"tokens_out":9773,"would_cite":false,"duration_ms":86205,"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":"This paper claims that in scalar-tensor theories the net tidal contribution to the gravitational-wave phase of a neutron-star inspiral is generally smaller than in general relativity, even though it involves three distinct Love numbers.","keywords":["scalar-tensor gravity","neutron star tidal deformability","Love numbers","gravitational-wave phase","post-Newtonian approximation","scalarization","binary inspiral","tidal phase"],"falsifier":"A concrete check is to run a numerical-relativity simulation of an equal-mass, equal-charge scalarized neutron-star binary, where $S_-=0$ so Eq. (83) reduces to the $c_3 x^3$ and $c_5 x^5$ terms: if the extracted tidal phase has the opposite sign or the wrong frequency scaling, the sign structure collapses. A second check is to recompute the phase after adding the Sec. 2.2 dynamical-friction force; if the 0.1-rad tidal difference shifts by a comparable amount, the phase model is incomplete.","tokens_in":45653,"feed_emoji":"🌊","tokens_out":10672,"duration_ms":94743,"temperature":0.7,"pith_summary":"This paper aims to establish that tidal effects in gravitational waves from neutron-star binaries are more varied, and in a specific sense weaker, in scalar-tensor gravity than in general relativity. Working to first post-Newtonian order and first order in finite-size effects, the authors derive a ready-to-use Fourier-domain phase formula for the inspiral in which scalar, tensor, and mixed scalar-tensor tidal deformabilities enter at different powers of the frequency parameter. Because those terms carry opposite signs, the net tidal phase correction in scalar-tensor theory is generally smaller than the corresponding general-relativistic correction for the same binary. The result matters because template models built on a single GR-like tidal deformability would miss the scalar and mixed tidal signatures and could misread the neutron-star equation of state.","feed_headline":"Tidal gravitational waves shrink under scalar-tensor gravity","feed_subtitle":"Three Love numbers with opposite signs change the inspiral phase, so GR-style tidal templates may misread neutron stars.","key_machinery":"The load-bearing structure is a skeletonized two-body action in which each neutron star carries three families of adiabatic, even-parity tidal Love numbers $\\lambda_S^\\ell$, $\\lambda_T^\\ell$, and $\\lambda_{ST}^\\ell$, responding respectively to scalar tides, tensor tides, and the cross-response between the two. These enter through tidal fields built from derivatives of the near-zone scalar and gravitational potentials, produce effective tidal coefficients $\\zeta_\\ell$, $\\bar{\\zeta}_\\ell$, and $\\tilde{\\zeta}_\\ell$ in the Lagrangian and in the radiative multipoles, and are converted into phasing by adiabatic energy balance between the binding energy and the combined scalar-plus-tensor flux. The final object is the stationary-phase Fourier phase of Eq. (83), with the inspiral split into dipole-driven and quadrupole-driven frequency domains; the quadrupole-driven version is the one relevant for current ground-based gravitational-wave detectors.","core_discovery":"The paper's central claim is that the tidal part of the Fourier-domain gravitational-wave phase in the quadrupolar-driven inspiral regime is $$\\psi_{\\rm tid} = \\frac{3}{128\\eta $x^{{5/2}}$}\\left[c_2 S_- $x^{2}$ + c_3 $x^{3}$ + c_4 S_-^2\\left(\\log x - \\frac{2}{3}\\right)$x^{4}$ + \\left(\\frac{39}{2\\$alpha^{5}$\\$xi^{2}$}\\tilde{\\Lambda} + c_5\\right)$x^{5}$\\right],$$ where $x=(G\\alpha M\\omega)^{2/3}$, $S_-=(q_A-q_B)/(2\\sqrt{\\alpha})$, $\\alpha=1+q_A q_B$, $\\xi=1+S_+^2\\alpha/6$, and $\\tilde{\\Lambda}$ is the mass-weighted combination of quadrupolar tensor deformabilities with the same functional form as in GR. The coefficients $c_2$ and $c_4$ vanish for equal-mass, equal-charge binaries; $c_3$ involves scalar dipolar Love numbers; $c_5$ and $\\tilde{\\Lambda}$ involve quadrupolar tensor, scalar, and mixed Love numbers; and in the general-relativistic limit only the $\\tilde{\\Lambda}$ term survives. Because the scalar and scalar-tensor contributions enter with opposite sign from the tensor ones, the authors conclude that the net tidal phase in scalar-tensor gravity is generally smaller in magnitude than in GR, while still depending on all three types of Love numbers.","pith_inferences":["If the sign structure of Eq. (83) holds in more complete calculations, then parameter estimation with GR templates containing a single tidal deformability would map a scalar-tensor signal onto a biased equation-of-state measurement, and the size of that bias is a direct injection-study prediction.","Beyond the paper, the same sign logic suggests a null test for modified gravity: a measured tidal phase whose frequency dependence cannot be fit by any single $\\tilde{\\Lambda}$ value would indicate the presence of the $x^2$, $x^3$, or $x^4$ terms.","At deci-Hertz frequencies, where the dipole-driven regime becomes accessible, the tidal terms grow with different frequency powers, so future space-based detectors could separate the scalar and mixed Love numbers cleanly.","The paper's own dynamical-friction estimate implies that scalar-cloud back-reaction may shift the phase by an amount comparable to the tidal signal, so the next step is to model that environment and test whether the predicted ST-minus-GR tidal difference survives."],"forward_implications":["Inspiral templates in scalar-tensor gravity must include all three Love numbers, because the scalar and mixed tidal terms scale as $x^2$, $x^3$, $x^4$, and $x^5$ and cannot be absorbed into a single GR-like deformability.","Equal-mass, equal-charge binaries have $S_-=0$, which switches off the dipolar tidal terms $c_2S_-x^2$ and $c_4S_-^2x^4$, so systems with maximum scalar-charge asymmetry are the ones that maximize those terms.","For the surveyed equations of state and couplings, the dipole-driven regime ends below about 3 Hz, so ground-based detectors should be modelled with the quadrupole-driven phase formula (83).","The scalar-tensor Love number $\\lambda_{ST}$ contributes at quadrupolar order at a level comparable to or larger than the pure scalar Love number, so it cannot be neglected once scalar tides are included.","The smaller net tidal phase in scalar-tensor gravity creates a potential degeneracy with a softer equation of state in GR, but the much larger point-particle phase difference between the theories should resolve it in most cases."],"supporting_citations":[{"why":"Supplies the three scalar, tensor, and scalar-tensor Love numbers and their signs for scalarized neutron stars, which are the inputs to the phase formula.","marker":"[45]"},{"why":"Defines the leading-order GR tidal phase and the combination to which the scalar-tensor result must reduce in the GR limit.","marker":"[42]"},{"why":"Provides the mass-weighted tidal deformability $\\tilde{\\Lambda}$ used in the $x^5$ term.","marker":"[86]"},{"why":"Establishes the dipolar tidal phasing framework and the dipole-to-quadrupole transition estimate that this paper generalizes to quadrupolar, scalar, and mixed tides.","marker":"[49]"},{"why":"Models leading-order dipolar tidal effects in scalar-tensor theories, which the paper complements with quadrupolar effects.","marker":"[48]"},{"why":"Supplies the scalar-wave multipole and energy-flux machinery used for the scalar radiation sector.","marker":"[64]"},{"why":"Supplies the 2PN point-particle scalar-tensor waveforms and flux expressions used as the nontidal baseline.","marker":"[65]"}],"fun_headline_variants":["Three Love numbers, one twist: smaller tidal GW phase in scalar-tensor gravity","Scalar-tensor gravity trims tidal wave imprint","Tidal Love numbers cancel out in scalar-tensor inspiral","Net tidal phase shrinks with scalar-tensor Love numbers","Opposing Love numbers weaken tidal fingerprints in GWs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model assumes that the inspiral's orbital decay is driven only by scalar and tensor radiation reaction, with no comparable environmental force; the paper's own estimate in Sec. 2.2 finds that dynamical friction from the scalar cloud around each star is about one tenth of the gravitational-wave flux near 200 Hz, which is the same order of magnitude as the tidal phase differences it predicts.","fun_headline_variants_meta":{"raw":{"variants":["Three Love numbers, one twist: smaller tidal GW phase in scalar-tensor gravity","Scalar-tensor gravity trims tidal wave imprint","Tidal Love numbers cancel out in scalar-tensor inspiral","Net tidal phase shrinks with scalar-tensor Love numbers","Opposing Love numbers weaken tidal fingerprints in GWs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000428,"raw_usage":{"total_tokens":2211,"prompt_tokens":989,"completion_tokens":1222,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":605,"completion_tokens_details":{"reasoning_tokens":1137}},"tokens_in":605,"tokens_out":1222,"duration_ms":9396,"temperature":1.0,"reasoning_tokens":1137,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T19:39:29.803093+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete check is to run a numerical-relativity simulation of an equal-mass, equal-charge scalarized neutron-star binary, where $S_-=0$ so Eq. (83) reduces to the $c_3 x^3$ and $c_5 x^5$ terms: if the extracted tidal phase has the opposite sign or the wrong frequency scaling, the sign structure collapses. A second check is to recompute the phase after adding the Sec. 2.2 dynamical-friction force; if the 0.1-rad tidal difference shifts by a comparable amount, the phase model is incomplete.","supporting_citations":[],"review_version":1}