{"id":"1f68b264-49b1-4d2a-9a6d-37190a6457b8","arxiv_id":"2501.07998","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In scalar-tensor theories, the 1/r^3 term used to extract neutron star tidal Love numbers contains a Love-number-independent contamination, computed here for minimally coupled and DEF scalar fields.","lead":"Gravitational wave measurements infer neutron star structure through tidal effects, but in modified gravity theories the standard extraction of the tidal response from the far field is ambiguous. This paper identifies an extra contribution to the far-field term that mimics a tidal response, computes it for the DEF scalar-tensor model, and shows it shifts the extracted tidal coefficients by up to about 15 percent in an excluded corner of parameter space.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The even-power parity assumption in Eq. (4.43) is the load-bearing step and is asserted rather than proven; a direct diagram check is needed, though the Fig. 5 internal consistency check already mitigates the risk.","rationale":"This paper's core contribution is the identification and quantification of a non-tidal 1/r^3 contamination in scalar-tensor Love-number extraction. Reading the argument in good faith, the chain is: (i) the EFT point-particle action defines the Love numbers; (ii) the asymptotic expansions are complete up to the four free coefficients; (iii) the only non-EFT structural input is the even-power parity structure of the zero-Love-number response in isotropic coordinates, Eq. (4.43); and (iv) the subsequent matching is lengthy but internally consistent. The numerical GR validation (Fig. 3) and the mixed-Love-number consistency check (Fig. 5) provide independent support. The latter is especially strong because the corrected lambda_hphi from the two extraction channels agrees to 0.07%, while the uncorrected values differ by up to 13%; this makes it unlikely that the non-tidal subtraction is substantially wrong. The residual concern is that Eq. (4.43) is justified by a symmetry of the linearized point-particle action plus a dimensional-analysis claim, rather than by explicit diagrammatic computation, and that non-linear point-particle operators could in principle break the parity. Because this is the one unproven structural input, I would not change the accept verdict, but I would recommend the direct diagram check described above as a follow-up. I agree with the reader that this is the weakest assumption and do not see a stronger or more decisive vulnerability in the argument.","tokens_in":44190,"tokens_out":20825,"duration_ms":237411,"concrete_test":"Run an explicit tree-level EFT calculation of the classical response to an external quadrupole in isotropic coordinates using the full point-particle action (4.22) and the bulk action (4.20), keeping all diagrams with up to five worldline couplings (Fig. 1c) and all derivative vertices, and read off the coefficients of r^1, r^-1, r^-3 and r^-5 in the perturbed fields for lambda = 0. If any such odd-power coefficient is non-zero, Eq. (4.43) fails and the non-tidal formulas (4.65)-(4.66) must be revised; if all vanish, the parity structure is confirmed and the central claim stands. As a complementary check, repeat the order-by-order matching of Sec. IVC keeping the full point-particle action rather than the linearized form (4.23) and verify that H0,3^0 and delta_phi_3^0 are unchanged.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central decomposition (4.72)-(4.73), and its DEF analog (5.25)-(5.26), separates the 1/r^3 coefficient into a tidal part and a non-tidal part. The non-tidal part (4.65)-(4.66) is obtained by matching the Schwarzschild-coordinate expansion to the isotropic-coordinate expansion (4.43), which asserts that for zero Love numbers the perturbed fields contain only even powers of the radial coordinate. This parity statement is the load-bearing step. The inversion symmetry (4.42) is checked for the linearized point-particle action (4.23), not for the full point-particle action (4.22) or (3.2), and the transition from an even number of worldline couplings to even powers of r is not demonstrated at the diagram level. Odd-power terms could enter through derivative vertices or through non-linear point-particle operators such as p phi^2/M_Pl^2; if they do, Eqs. (4.65)-(4.66) and hence (5.32)-(5.35) would be incomplete. The paper's Fig. 5 check, where lambda_hphi obtained from c_hphi agrees with lambda_hphi obtained from c_phih to better than 0.07% after correction, is strong indirect evidence that the non-tidal subtraction is correct, because a wrong subtraction would generically spoil this agreement. Still, the parity assumption itself is asserted rather than proven.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper analyzes the extraction of quadrupolar tidal Love numbers for neutron stars in Horndeski scalar-tensor theories. Building on the perturbation equations of Ref. [57], the authors derive linear static even-parity perturbation equations and show that the asymptotic 1/r^3 coefficients of the metric and scalar perturbations contain, in addition to the tidal Love numbers, a non-tidal contribution that is independent of the Love numbers. Using an effective field theory (EFT) point-particle approach, they compute this contamination for a minimally coupled scalar field (Sec. IV) and for the Damour–Esposito–Farèse (DEF) model (Sec. V), obtaining closed-form expressions for the corrected Love numbers. They numerically compute the DEF Love numbers with the DD2 equation of state and demonstrate that neglecting the contamination changes the extracted Love numbers by up to 15% for β = −6. They also discuss why the same method cannot be straightforwardly applied to scalar-Gauss-Bonnet gravity (Sec. VI).","tokens_in":44450,"tokens_out":6129,"duration_ms":63392,"significance":"If the central derivation holds, the paper resolves a previously overlooked ambiguity in computing tidal Love numbers in scalar-tensor theories. The result that the 1/r^3 asymptotic coefficient is not purely tidal is conceptually important and has direct implications for gravitational-wave parameter estimation. The paper's strengths include the derivation of explicit, ready-to-use formulas (Eqs. (4.72)–(4.73), (5.32)–(5.35)), the validation of the numerical extraction against Hinderer's exact GR formula (Fig. 3), and the nontrivial internal consistency check in Fig. 5, where the mixed Love number λ_hϕ extracted via two independent channels agrees to better than 0.07% after correction. The presentation is detailed and the appendices provide extensive asymptotic expansions that would be useful for future work.","major_comments":[{"comment":"The even-power parity assumption in Eq. (4.43) is load-bearing for the central claim, but it is asserted rather than proven. The inversion symmetry (4.42) is verified only for the truncated linearized point-particle action (4.23), not for the full action (4.22) or for generic operators such as the p φ^2/M_Pl^2 term introduced in Eq. (3.2). Such operators break the symmetry and could in principle contribute to the 1/r^3 coefficient, altering the non-tidal parts (4.65)–(4.66) and the resulting Love-number formulas (5.32)–(5.35). The footnote justifying the truncation rules out some classes of diagrams but does not explicitly exclude p-type operators at the relevant order. I request a diagram-level proof, or at least a clear statement of the conditions under which Eq. (4.43) holds. The numerical agreement in Fig. 5 is reassuring, but it tests only the DEF model with the specific PP action used and does not constitute a general proof.","section":"Sec. IVB, Eq. (4.43)"},{"comment":"The universality claim stated after Eq. (4.73)—that the result applies to any scalar-tensor theory reducing to Einstein gravity with a minimally coupled scalar in vacuum—is conditional on the specific truncation of the point-particle action in Eq. (4.23). If higher-derivative or nonlinear worldline operators are present, the inversion symmetry (4.42) is generically broken and the decomposition into tidal and non-tidal parts may need to be revisited. The authors should either prove that such operators cannot contribute to the 1/r^3 coefficient (for example by explicit diagrammatic counting) or explicitly restrict the statement of universality to the class of models where those operators are absent. This clarification is necessary because the paper's main phenomenological conclusions for the DEF model rely on this decomposition.","section":"Sec. IVC, Eqs. (4.72)–(4.73)"}],"minor_comments":[{"comment":"There is a typo: \"neither neither G4,X nor G5 vanish\" should read \"neither G4,X nor G5 vanish\".","section":"Sec. II, first paragraph"},{"comment":"The footnote says the calculation is \"not directly comparable to the expansion in Eq. (2.7)\", but Eq. (2.7) is a background equation; the intended reference is likely Eq. (4.7).","section":"Sec. IVB, footnote 3"},{"comment":"The phrase \"O(1 ∼ 10) %\" is nonstandard notation; it should be written as \"O(1–10)%\".","section":"Abstract"},{"comment":"The text refers to \"the coefficients H0,3 and δϕ3 of Eqs. (4.13) and (4.14)\", but Eqs. (4.13)–(4.14) are expansions; the wording should be \"the coefficients H0,3 and δϕ3 in Eqs. (4.13) and (4.14)\".","section":"Sec. IVC, paragraph after Eq. (4.63)"},{"comment":"The axis label uses the symbol φ (Greek phi without subscript) in \"|λhφ/ ˆλhφ − 1|\"; for consistency with the text, the subscript should be the scalar-field symbol ϕ.","section":"Fig. 5"},{"comment":"A comma is missing: \"as long as one uses Schwarzschild coordinates combined with the Regge-Wheeler gauge — in scalar-tensor theory\" should have a comma after \"gauge\".","section":"Sec. I, last paragraph of introduction"}],"recommendation":"major_revision","confidential_remarks":"The paper is well within the scope of Physical Review D and addresses a timely topic with high practical relevance. The central claim is plausible and supported by nontrivial numerical cross-checks, but the parity assumption in Sec. IVB needs a more rigorous justification or a clear statement of its domain of validity. I believe this is fixable within the manuscript's scope, so I recommend major revision rather than rejection. The author self-citations to Ref. [57] and related work are appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know about this paper is that it settles a real ambiguity. In scalar-tensor theories, the 1/r^3 coefficient in the asymptotic expansion of the perturbations is not purely tidal; there is an extra Love-number-independent contribution at the same order. The authors compute it for a minimally coupled scalar and for the DEF model, and it shifts the extracted Love numbers by a few percent at the observationally allowed beta = -4.5 and by up to 15% at the excluded beta = -6. That matters for next-generation GW tests.\n\nThe contamination terms, Eqs. (4.65)-(4.66) and their DEF counterparts, are genuinely new relative to the earlier scalar-tensor Love number papers. The trick of using the inversion symmetry to avoid diagram-by-diagram computation is clever and works well here. The derivation is long but internally consistent. The numerical checks are strong: Fig. 3 recovers Hinderer's exact GR Love numbers to better than 0.01% with the r^-8 fit, and Fig. 5 shows that the mixed Love number extracted two independent ways agrees to 0.01-0.07% after correction, while the uncorrected values disagree by up to 13%. That agreement is a non-trivial test of the subtraction.\n\nThe soft spot is the parity assumption in Eq. (4.43): that for vanishing Love numbers the perturbed fields in isotropic coordinates contain only even powers of the radial coordinate. This is inferred from the inversion symmetry together with dimensional analysis, not proven by explicit diagram counting. The symmetry is checked for the linearized point-particle action (4.23), not the full action (4.22) or (3.2), so odd-power terms could in principle sneak in through derivative vertices or non-linear operators like p phi^2. If they did, the contamination formulas would be incomplete. But the stress-test concern is mitigated: a wrong subtraction would generically spoil the Fig. 5 agreement, and it does not. So this is a caveat, not a fatal flaw. A direct diagram check in a follow-up would close the gap.\n\nThis paper is for people doing precision GW tests of scalar-tensor gravity and for EFT practitioners working on tidal response. It deserves a serious referee. I would send it to review and, in the report, ask the authors to justify the even-power parity more carefully, perhaps with an appendix or a representative diagram check. But the central claim holds up.","headline":"Careful, important resolution of the 1/r^3 extraction ambiguity in scalar-tensor Love numbers; the parity assumption is the main soft spot, but the internal consistency check largely covers it.","tokens_in":45014,"tokens_out":2118,"would_cite":true,"duration_ms":22922,"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 shows that in scalar-tensor theories the $1/r^3$ coefficient in the asymptotic perturbation fields contains a contamination independent of the tidal Love numbers, and derives corrected formulas for the minimally coupled scalar…","keywords":["tidal Love numbers","neutron stars","scalar-tensor theories","Horndeski gravity","effective field theory","asymptotic expansion","spontaneous scalarization","gravitational waves"],"falsifier":"Compute the five-worldline-coupling diagrams of Fig. 1c explicitly in isotropic coordinates: the paper's parity argument predicts they vanish, so any nonzero result would falsify the even-power structure (4.43) and with it the contamination formulas (4.65)-(4.66).","tokens_in":43978,"feed_emoji":"🌌","tokens_out":15103,"duration_ms":121539,"temperature":0.7,"pith_summary":"This paper addresses how tidal Love numbers of neutron stars should be extracted from perturbation theory in Horndeski scalar-tensor theories, the most general scalar-tensor theories with second-order field equations. It establishes that the standard identification of the $1/r^3$ term in the asymptotic expansion of the perturbation fields with the quadrupolar tidal response is contaminated: an additional contribution, independent of the Love numbers, appears at the same order. The authors compute this contamination for a minimally coupled massless scalar field and for the DEF model (a scalar-tensor theory with nonminimal coupling $F(\\phi)=e^{-\\beta\\phi^2/(2M_{\\rm Pl}^2)}$), and derive corrected extraction formulas. They find that neglecting the contamination shifts the DEF Love numbers by up to about 15% for $\\beta = -6$ and by a few percent for $\\beta = -4.5$, which matters for gravitational-wave parameter estimation.","feed_headline":"Neutron-star Love numbers carry a hidden non-tidal term","feed_subtitle":"The 1/r^3 tail mixes tidal and scalar-charge effects, shifting extracted Love numbers by up to 15 percent.","key_machinery":"The load-bearing tool is an effective-field-theory point-particle action in isotropic coordinates. The bulk plus background worldline action is invariant under inversion of both fields and couplings, $\\hat\\varphi \\to -\\hat\\varphi$, $\\hat\\phi \\to -\\hat\\phi$, $M \\to -M$, $Q \\to -Q$, which together with dimensional analysis forces the perturbation fields $\\delta\\hat\\varphi$ and $\\delta\\hat\\phi$ of an object with zero Love numbers to contain only even powers of the radial coordinate (Eq. (4.43)). Transforming the EFT expansion back to Schwarzschild coordinates and matching through the gauge-invariant combination $\\Psi = H_0 - f'\\delta\\phi/(f\\phi')$ converts that even-power structure into explicit $1/r^3$ contamination terms, which are then separated from the tidal contributions $\\lambda_{hh}$, $\\lambda_{h\\phi}$, and $\\lambda_{\\phi\\phi}$ read off from the point-particle tidal action.","core_discovery":"The central claim is that in scalar-tensor theories the coefficient of the $1/r^3$ term in the asymptotic expansion of the metric and scalar perturbations is not purely tidal. Writing $H_0(r) = H_{0,-2} r^2 + \\dots + H_{0,3}/r^3 + \\dots$ and $\\delta\\phi(r) = \\delta\\phi_{-2} r^2 + \\dots + \\delta\\phi_3/r^3 + \\dots$, the paper splits $H_{0,3} = H^\\lambda_{0,3} + H^0_{0,3}$ and $\\delta\\phi_3 = \\delta\\phi^\\lambda_3 + \\delta\\phi^0_3$, where the superscript $\\lambda$ contains the tidal Love numbers and the superscript $0$ is a non-tidal part built only from the mass $M$ and scalar charge $\\phi_1$. For the minimally coupled scalar this non-tidal part is $H^0_{0,3} = H_{0,-2}(3 M^3 \\phi_1^2/(5 M_{\\rm Pl}^2) - 3 M \\phi_1^4/(80 M_{\\rm Pl}^4))$, with the same expression for $\\delta\\phi^0_3$ in terms of $\\delta\\phi_{-2}$. The DEF-model generalization is given in Eqs. (5.25)-(5.26), and the corrected Love-number extraction formulas in Eqs. (5.32)-(5.35).","pith_inferences":["The same even-power argument would apply to dipolar and higher-multipole Love numbers in massless scalar-tensor theories, suggesting the ambiguity is not specific to the quadrupole and could be resolved by the same EFT matching.","Because the contamination terms are built from the mass $M$ and the scalar charge $\\phi_1$, their fractional importance grows with compactness and scalar charge, so published DEF Love-number tables that ignored them may need revision for strongly scalarized stars.","A practical diagnostic for future work would be to check whether the two independent extractions of $\\lambda_{h\\phi}$ agree: if they disagree at the percent level or more, the non-tidal $1/r^3$ contamination is present and must be subtracted.","For massive scalar fields the contamination would be exponentially suppressed at infinity, so the correction is relevant mainly for massless or extremely light scalars, which could guide which scalar-tensor models need the corrected extraction."],"forward_implications":["In any scalar-tensor theory whose exterior reduces to Einstein gravity with a massless minimally coupled scalar, the corrected formulas (4.72)-(4.73) must be used to extract $\\lambda_{hh}$, $\\lambda_{h\\phi}$, and $\\lambda_{\\phi\\phi}$ from the $1/r^3$ coefficients.","For DEF neutron stars with the DD2 equation of state, using the raw $1/r^3$ coefficients instead of the corrected formulas shifts $\\Lambda_{hh}$ by up to about 3%, $\\Lambda_{\\phi\\phi}$ by up to about 5%, and $\\Lambda_{h\\phi}$ by up to 15% for $\\beta=-6$; for $\\beta=-4.5$ the shifts are a few percent or less.","The corrected extraction makes $\\lambda_{h\\phi}$ obtained from the gravitational channel agree with that from the scalar channel to within numerical accuracy, whereas the uncorrected extraction disagrees by up to 13%.","On the GR branch of DEF the mixed Love number $\\Lambda_{h\\phi}$ vanishes as parity requires, while $\\Lambda_{\\phi\\phi}$ stays nonzero, and spontaneous scalarization strongly enhances $\\Lambda_{\\phi\\phi}$ for $\\beta=-6$.","For scalar-Gauss-Bonnet gravity the inversion symmetry is broken, so the contamination cannot be fixed by the same parity argument and requires explicit diagrammatic computation."],"supporting_citations":[{"why":"supplies the exact general-relativistic quadrupolar Love number formula used as the validation baseline for the numerical matching procedure.","marker":"[27]"},{"why":"derives the gauge-ready linear perturbation equations for relativistic stars in Horndeski theories that this paper extends to static even-parity tides.","marker":"[57]"},{"why":"first noted the possible ambiguity in extracting Love numbers from asymptotic coefficients.","marker":"[60]"},{"why":"analyzes the ambiguity in relativistic tidal deformability that motivates the contamination computation.","marker":"[61]"},{"why":"is the prior scalar-tensor Love-number calculation whose extraction did not include the non-tidal $1/r^3$ contribution.","marker":"[64]"},{"why":"provides the effective-field-theory point-particle action for extended objects on which the EFT side is built.","marker":"[74]"},{"why":"supplies the isotropic-coordinate EFT treatment and the inversion-symmetry/even-power argument used to fix the contamination.","marker":"[78]"},{"why":"introduces the DEF spontaneous-scalarization model to which the corrected formulas are applied.","marker":"[51]"}],"fun_headline_variants":["Love numbers hide a scalar charge term","Neutron-star Love numbers off by up to 15%","Tidal Love numbers get a non-tidal twist","Scalar charge skews neutron-star Love numbers","New Love number formula for scalar-tensor gravity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation rests on the claim that for an object with vanishing Love numbers the perturbation fields in isotropic coordinates contain only even powers of the radial coordinate, which follows from inversion symmetry plus dimensional analysis.","fun_headline_variants_meta":{"raw":{"variants":["Love numbers hide a scalar charge term","Neutron-star Love numbers off by up to 15%","Tidal Love numbers get a non-tidal twist","Scalar charge skews neutron-star Love numbers","New Love number formula for scalar-tensor gravity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000711,"raw_usage":{"total_tokens":3291,"prompt_tokens":1126,"completion_tokens":2165,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":742,"completion_tokens_details":{"reasoning_tokens":2091}},"tokens_in":742,"tokens_out":2165,"duration_ms":14717,"temperature":1.0,"reasoning_tokens":2091,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:32:55.494006+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the five-worldline-coupling diagrams of Fig. 1c explicitly in isotropic coordinates: the paper's parity argument predicts they vanish, so any nonzero result would falsify the even-power structure (4.43) and with it the contamination formulas (4.65)-(4.66).","supporting_citations":[],"review_version":1}