{"id":"3e48db4a-c81c-43d7-81a7-4e50cc79ba6c","arxiv_id":"2506.08722","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A general-relativistic calculation of the quadratic tidal constant p2 for polytropic neutron stars shows that nonlinear tides lower the effective dynamical resonance frequency by 10 to 15 percent, confirming Newtonian mode-based results without using normal modes.","lead":"The authors compute a new relativistic tidal constant for neutron star models and show that nonlinear tidal effects can shift the star's dynamical resonance frequency down by as much as 15 percent. This matters for gravitational wave astronomy, because an earlier onset of tidal resonance would measurably change the waveforms from neutron star inspirals and affect how the dense-matter equation of state is inferred.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 15% shift in omega_* and the enhanced tidal response are not consequences of the computed p_2; they rest on the unvalidated single-pole resummation of Eq. (8.9).","rationale":"The paper's genuine new result is the computation of the nonlinear tidal constant p_2 for relativistic polytropes, with a detailed numerical implementation and a comparison of independent computations agreeing to six significant digits. That part deserves credit. However, the central physical message—that nonlinearities lower the effective resonance frequency by as much as 15% and thereby enhance the tidal response during inspiral—does not follow from p_2 alone. It follows only after Eq. (8.8) is replaced by the single-pole form Eq. (8.9), a step the authors describe as pragmatic and inspired by the Newtonian mode picture. The reader's weakest assumption identified exactly this continuation as the load-bearing premise, and I concur. No independent validation exists for the nonlinear case; the cited validation in Ref. [35] is for p_2 = 0. The proposed concrete test, computing the fourth-order time-derivative coefficient, would directly probe whether the true response has the assumed analytic structure. If the omega^4 coefficient matches the geometric-series prediction, the ansatz is credible; if not, the headline claim is unsupported by the calculation. The Appendix A surface-layer divergence for n < 1 is a second concern, but it does not directly enter the headline figures (which use n >= 1) and is acknowledged by the authors; it reinforces the conditional nature of the results without replacing the continuation as the principal issue. Therefore the reader's CONDITIONAL verdict is appropriate, and no change is needed.","tokens_in":24773,"tokens_out":4840,"duration_ms":61096,"concrete_test":"Compute the next-order term in the linear dynamical response, i.e. the coefficient of omega^4 in tilde_k_2(omega), using the frequency-domain methods of Ref. [35] extended to fourth time derivatives. If the Taylor expansion of tilde_k_2(omega) obtained from the computed k_2, kdot_2, and this new coefficient matches the expansion of k_2 / (1 - omega^2 / omega_*^2) to within ~10% in the omega^4 term, the single-pole ansatz is supported and the earlier-resonance claim stands; if the coefficient deviates by more than that, Eq. (8.9) is not justified and the central physical conclusion becomes conditional on an unverified resummation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline quantity omega_* is defined in Eq. (1.10) from k_2, kdot_2, and p_2, but the physical claim that the approach to resonance occurs earlier is obtained only after the low-frequency expansion (8.8) is resummed to the single-pole form (8.9), tilde_k_2 = k_2 / (1 - omega^2 / omega_*^2). This resummation is an ansatz: the paper provides no relativistic derivation, and the existing validation in Ref. [35] covers only the linear case p_2 = 0, where omega_* is compared with the f-mode frequency. For p_2 != 0 there is no check that higher-order terms in the simultaneous time-derivative and nonlinearity expansion are subdominant up to omega approximately omega_*, or that the exact response is analytic with a simple pole. Consequently the up-to-15% reduction in omega_* and the large enhancement in Fig. 5 are properties of the ansatz, not of the computed p_2; the paper itself labels the step 'pragmatic' (Sec. I E) and an 'extension' (Sec. VIII). Appendix A raises a related but separate issue: for n < 1 polytropes the jump in extrinsic curvature diverges at the surface, which is not used in the headline figures but signals that the perturbative scheme is not fully under control for the softest model of Fig. 1.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes the nonlinear (quadratic) relativistic tidal constant p2, together with p4, for polytropic neutron-star models, using a second-order static perturbation scheme matched to the exterior metric of Ref. [36]. It combines these with the previously computed dynamic constant kdot2 to define a frequency parameter omega_* through Eq. (1.10), and reports that the nonlinear term lowers omega_* by up to 15% relative to the linear estimate. The paper interprets this as an earlier approach to resonance and an enhanced tidal response during binary inspiral, in qualitative agreement with the Newtonian mode-coupling results of Yu et al.","tokens_in":25152,"tokens_out":10145,"duration_ms":111856,"significance":"If the main physical claim holds, the paper is a valuable general-relativistic confirmation of the importance of nonlinear dynamical tides without relying on a mode decomposition. The strengths are substantial: the second-order perturbation equations are given explicitly (Sec. IV), the junction conditions and tidal-moment redefinition invariance are discussed carefully (Secs. VI-VII), and Table I reports agreement to six significant digits between two independent computations. The frequency parameter is constructed from independently computed constants rather than fitted, and the paper is transparent about the pragmatic nature of the response-function extension. However, the physical conclusion about earlier resonance depends on an additional, only partially validated ansatz, so the significance is conditional pending that validation.","major_comments":[{"comment":"The central physical conclusion is not a direct consequence of the computed p2 alone. Equation (8.8) is a low-frequency expansion valid to O(Omega^2), and replacing it by the single-pole form (8.9), then interpreting omega_* as a resonance frequency, is an additional ansatz. The only validation cited (Ref. [35]) is for p2 = 0, where omega_* agrees with the f-mode; no test is provided for p2 != 0. Since Figs. 4 and 5 and the statement that the approach to resonance occurs earlier all depend on this continuation, the manuscript should either supply an independent check (for example, higher-order terms in the simultaneous expansion, comparison with nonlinear mode calculations, or numerical-relativity data) or explicitly demote these results to model-dependent predictions of the one-pole ansatz.","section":"Sec. VIII, Eq. (8.9); Sec. I E, Eqs. (1.8)-(1.10)"},{"comment":"For polytropes with n < 1 the extrinsic-curvature jump diverges at the stellar surface, so the junction conditions [K_ab] = 0 used in Sec. VI are not satisfied for these models. The rebuttal in Appendix A is an analogy with the density expansion, not a demonstration that the divergence is an artifact of the perturbative expansion. Because Fig. 1 and Table I present p2 values for n = 0.5, those values are not supported unless the divergence is shown to be harmless; alternatively the n = 0.5 curve should be removed or explicitly labeled as tentative.","section":"Appendix A, Eqs. (A1)-(A2)"}],"minor_comments":[{"comment":"The heading contains a typo: \"quadropole\" should be \"quadrupole\".","section":"Sec. IV heading"},{"comment":"The phrase \"the joint occurs at the deformed surface\" should read \"the junction occurs at the deformed surface\".","section":"Sec. VI"},{"comment":"The phrase \"exiting developments\" should be \"exciting developments\".","section":"Sec. I A"},{"comment":"The n = 0.5 model appears in Fig. 1 and Table I but not in the frequency-parameter figures because kdot2 is unavailable for n = 0.5; the text should state this mismatch explicitly to avoid the impression that omega_* was computed for n = 0.5.","section":"Figs. 1 and 4"},{"comment":"The definition of omega_* mixes the physically static p2 term with the dynamical kdot2 term; although the text justifies this as natural, a one-sentence reminder that the p2 contribution is static would help readers distinguish the two effects.","section":"Sec. I E, Eq. (1.10)"}],"recommendation":"major_revision","confidential_remarks":"The paper's novel technical content is the computation of p2 and p4 and their extraction in a redefinition-invariant way. The resonance-related conclusions rest on the one-pole continuation that is validated only in the linear case. In my view the manuscript is publishable after either a genuine validation of the nonlinear continuation or a clear reframing of the 15% frequency shift and Fig. 5 as predictions of the ansatz rather than established consequences. The n < 1 junction-condition issue also needs to be resolved or the affected results explicitly qualified."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things. First, the quadratic tidal constants p2 and p4 for relativistic polytropes are genuinely new, and the calculation looks careful—two independent numerical routes agreeing to six digits is real evidence. Second, the flashy result—that nonlinearities shift the effective resonance frequency by up to 15%—is not a direct measurement from the perturbation equations; it comes from plugging the computed constants into a single-pole ansatz for the response function. The paper says so (\"pragmatic extension\"), but the abstract and intro state the 15% as if it were a settled prediction.\n\nThe genuinely new contribution is the p2 computation itself. The framework is inherited from Poisson 2021 and the linear dynamics from Pitre & Poisson 2024; the new step is solving the second-order perturbation equations and extracting p2 and p4. The equations are explicit, the handling of the tidal-moment redefinition ambiguity is transparent, and the table gives numbers. That is solid, reproducible work by the standards of this field. The 6-digit agreement between two computations is particularly reassuring.\n\nThe main soft spot is exactly what the stress-test flags. Eq. (8.9) is an ansatz. The low-frequency expansion (8.8) is valid for small omega; the one-pole resummation is validated only for p2 = 0, against the f-mode, in the previous paper. Nothing here checks that for p2 ≠ 0 the exact response remains one-pole up to omega ≈ omega_star. Without that check, Fig. 5 is an illustration of the ansatz, not a prediction. The paper is not hiding this, but a referee should push them to state the caveat more prominently, and ideally to test against Yu et al.'s Newtonian mode calculation or available NR results.\n\nSecond minor issue: Appendix A shows the junction conditions for n < 1 polytropes have a divergent surface-layer term. The authors argue it is an artifact of expanding at the unperturbed surface, and they keep n ≥ 1 for the headline plots. That is reasonable, but it does mean the perturbative scheme is not fully controlled at the softest end of their parameter range.\n\nThe citation pattern is fine: they lean on their own earlier work, but those papers are the actual foundations here, and they cite Yu et al. and the numerical-relativity literature.\n\nBottom line: this deserves serious peer review. The p2 values will be useful even if the resonance claim is eventually softened. The appropriate referee request is: keep the calculation, discuss the domain of validity of the single-pole form, and align the abstract with the conditional nature of the physical conclusion. I would bring it to reading group and would cite it for the p2 table.","headline":"The new GR computation of p2 is solid and worth citing; the flashier 'earlier resonance' claim is a conditional consequence of an admittedly pragmatic single-pole resummation, so the abstract overstates the physics.","tokens_in":25611,"tokens_out":3780,"would_cite":true,"duration_ms":43407,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C35","83C25","85A15"],"pacs":["04.30.-w","04.25.Nx","97.60.Jd"],"model":"deepseek-v4-flash","headline":"Nonlinear tidal effects lower a neutron star's resonant frequency parameter by as much as 15 percent, so the approach to resonance begins earlier during a binary inspiral and the tidal response is amplified.","keywords":["dynamical tides","neutron stars","tidal deformability","general relativity","nonlinear tides","polytropic equations of state","gravitational waves","tidal response function"],"falsifier":"Compare the response function predicted by Eq. (1.8) with a high-precision numerical-relativity simulation of a polytropic neutron star in a tidal field whose frequency is swept up to $\\omega_*$; if the simulated $\\tilde{k}_2(\\omega)$ deviates substantially from $k_2(1-\\omega^2/\\omega_*^2)^{-1}$ before the peak, the earlier-resonance conclusion fails.","tokens_in":24539,"feed_emoji":"🌊","tokens_out":5730,"duration_ms":62877,"temperature":0.7,"pith_summary":"This paper argues that nonlinear aspects of the tidal deformation of a neutron star are comparable in size to the dynamical (time-derivative) corrections, and that they cannot be ignored when modeling dynamical tides in a binary inspiral. Working in full general relativity and without calling on a normal-mode decomposition, it computes the quadratic tidal constant $p_2$ for relativistic polytropes and combines it with the static and dynamic constants $k_2$ and $\\ddot{k}_2$ into a single-pole response function. The result is a frequency parameter $\\omega_*$ that is up to 15 percent smaller than the purely linear estimate, meaning the approach to resonance starts earlier in the inspiral and the tidal response is enhanced. This confirms, in a relativistic setting, the Newtonian mode-based finding of Yu et al. and shows that dynamical and nonlinear tides come together.","feed_headline":"Nonlinear tides pull neutron-star resonance earlier by 15%","feed_subtitle":"General-relativistic computation shows quadratic tidal effects cut the resonance frequency, boosting the star's tidal response during…","key_machinery":"The machinery is the simultaneous time-derivative and nonlinear expansion of the spacetime metric of a tidally deformed body, expressed through perturbation variables for quadrupole, hexadecapole, and monopole sectors. The load-bearing objects are the three tidal constants $k_2$, $\\ddot{k}_2$, and $p_2$ appearing in the relation between the mass quadrupole moment and the tidal moment, together with the frequency-domain response function $\\tilde{k}_2(\\omega)$, which the paper pragmatically extends from the low-frequency expansion $\\propto 1+\\omega^2/\\omega_*^2$ to the one-pole form $(1-\\omega^2/\\omega_*^2)^{-1}$. This one-pole form is what converts the computed constants into the prediction that resonance is approached earlier.","core_discovery":"The central claim is that the nonlinear tidal constant $p_2$ is numerically comparable to the dynamic constant $\\ddot{k}_2$, so that the effective resonance frequency $\\omega_*$, defined by $\\omega_*^2 = k_2/(\\ddot{k}_2 + \\frac14 p_2 M'/(M+M') GM/R^3)$, is lower than a linear treatment would predict. For polytropic equations of state the ratio of nonlinear to linear frequency parameters ranges from about 0.85 for low compactness to about 0.94 near the maximum mass. Because the tidal response function $\\tilde{k}_2(\\omega) = k_2 (1-\\omega^2/\\omega_*^2)^{-1}$ grows as $\\omega$ approaches $\\omega_*$, the earlier resonance produces a strong magnification of the tidal response. The paper also computes the hexadecapole quadratic constant $p_4$, with similar qualitative behavior.","pith_inferences":["The one-pole continuation is only checked against the f-mode frequency in the linear ($p_2=0$) case; testing the same formula against a direct numerical-relativity computation of the tidal response for the same polytrope, with nonlinear terms included, would settle whether the 15 percent shift is physical.","The near-universality of $p_2/\\ddot{k}_2$ across polytropic indices suggests the 6–15 percent resonance shift may persist for realistic equations of state, but this paper does not compute those.","Because the nonlinear term also contributes to a time-independent piece of the quadrupole moment (Eq. 8.7a), nonlinear tides may leave a detectable imprint in the inspiral phasing even before the resonance regime; the paper does not develop that consequence."],"forward_implications":["For equal-mass binaries, the nonlinear correction lowers $\\omega_*$ by roughly 15 percent for low-compactness stars and by roughly 6 percent near the maximum-mass configuration.","At a given orbital frequency, the response function $\\tilde{k}_2(\\omega)$ is magnified relative to the linear prediction because the one-pole at $\\omega_*$ is closer (Fig. 5).","The reduction is larger when the companion is more massive, through the factor $M'/(M+M')$ in Eq. (1.10), so the effect is strongest for asymmetric mass ratios.","Gravitational waveforms from inspiralling neutron-star binaries should include a nonlinear tidal phase correction of order 10–20 percent relative to the linearized description, consistent with the Newtonian result of Yu et al.","The same expansion yields the hexadecapole quadratic constant $p_4$, whose behavior parallels $p_2$."],"supporting_citations":[{"why":"Supplies the Newtonian, mode-based nonlinear tidal result that this paper confirms and generalizes to relativity.","marker":"[23]"},{"why":"Provides the dynamic tidal constant $\\ddot{k}_2$, the prior no-mode relativistic framework, and the linear-case validation of the one-pole response against the f-mode frequency.","marker":"[35]"},{"why":"Constructs the exterior metric and the post-Newtonian quadrupole-moment definition used to extract the tidal constants and formulate Eq. (1.6).","marker":"[36]"},{"why":"Supplies the Israel junction conditions used to join interior and exterior solutions at the deformed stellar surface.","marker":"[39]"},{"why":"Provides the standard expression for the tidal quadrupole moment of a circular orbit used in Sec. VIII to derive the frequency-domain response.","marker":"[40]"}],"fun_headline_variants":["Nonlinear tides shrink neutron-star resonance by 15%","Tidal nonlinearity cuts neutron-star resonance frequency","Quadratic tides boost neutron-star inspiral response","Nonlinear effects shift neutron-star resonance earlier","Relativistic tides: nonlinearity lowers resonance by 15%"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that the approximate one-resonance formula remains accurate up to the resonance frequency even after nonlinear terms are included; the paper only checks this approximation in the simpler linear case, where it reproduces the f-mode frequency.","fun_headline_variants_meta":{"raw":{"variants":["Nonlinear tides shrink neutron-star resonance by 15%","Tidal nonlinearity cuts neutron-star resonance frequency","Quadratic tides boost neutron-star inspiral response","Nonlinear effects shift neutron-star resonance earlier","Relativistic tides: nonlinearity lowers resonance by 15%"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000603,"raw_usage":{"total_tokens":2876,"prompt_tokens":1070,"completion_tokens":1806,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":686,"completion_tokens_details":{"reasoning_tokens":1730}},"tokens_in":686,"tokens_out":1806,"duration_ms":15625,"temperature":1.0,"reasoning_tokens":1730,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T05:04:35.093163+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compare the response function predicted by Eq. (1.8) with a high-precision numerical-relativity simulation of a polytropic neutron star in a tidal field whose frequency is swept up to $\\omega_*$; if the simulated $\\tilde{k}_2(\\omega)$ deviates substantially from $k_2(1-\\omega^2/\\omega_*^2)^{-1}$ before the peak, the earlier-resonance conclusion fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Newtonian, mode-based nonlinear tidal result that this paper confirms and generalizes to relativity."},{"cited_title":"Poisson, Compact body in a tidal environment: New types of relativistic Love numbers, and a post-Newtonian opera- 22 tional definition for tidally induced multipole moments, Phys","cited_arxiv_id":null,"evidence_quote":"Constructs the exterior metric and the post-Newtonian quadrupole-moment definition used to extract the tidal constants and formulate Eq. (1.6)."},{"cited_title":"Poisson and C","cited_arxiv_id":null,"evidence_quote":"Provides the standard expression for the tidal quadrupole moment of a circular orbit used in Sec. VIII to derive the frequency-domain response."}],"review_version":1}