{"id":"70739a51-3172-4161-a532-c90a6e65ac2c","arxiv_id":"2504.18403","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":4,"one_line_summary":"Adding a quadratic radial correction changes the predicted lower QPO frequency near the ISCO, but fits to eight neutron stars show it is insufficient to rescue the relativistic precession model.","lead":"The authors extend the relativistic precession model of quasi-periodic oscillations by adding a quadratic anharmonic term to the radial oscillation equation, which shifts the predicted lower frequency near the innermost stable circular orbit. They fit eight neutron star sources and find the correction is real but does not resolve the model's tension with the Kerr hypothesis.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Fitted oscillation amplitudes violate the small-displacement condition on which the anharmonic expansion is based, so the claim that anharmonic corrections are physically necessary is not supported.","rationale":"The reader's weakest assumption correctly identifies that the fitted values of A0 are too large for the perturbative expansion to be valid. This is the most load-bearing concern because the paper's novel quantitative claim is that anharmonic corrections are physically necessary to describe eight NS QPO sources; that claim is supported only by MCMC fits in which A0 reaches several kilometers, comparable to or larger than the orbital radius and to the ISCO. The tests proposed here directly probe the internal consistency of the model: if the conditions on A0 are imposed, the anharmonic models may lose their statistical advantage, in which case the central conclusion would collapse. The theoretical derivation of the quadratic term itself is a straightforward Taylor expansion and is not in question, but the step from Eq. (15) to Eq. (16) is a secondary logical gap that warrants attention. Overall, the paper remains worth publishing conditionally, provided the authors restrict claims to the regime where the anharmonic expansion is valid and make the code and data available, as the reader already required.","tokens_in":23770,"tokens_out":16870,"duration_ms":162788,"concrete_test":"For every best-fit row in Table I, solve the model relation f_U = f_phi(r0; M, j) for r0 at each observed f_U^k, then compute alpha0(r0), alpha1(r0) from Eq. (27) and epsilon = alpha1 * A(r0). Verify the consistency condition |epsilon| <= 0.1 * alpha0 and A(r0) < r0 - r_ISCO for all rows; count how many rows violate. Then rerun the MCMC for the A1/A2/A3 models with a hard prior enforcing |alpha1 A0| <= 0.1 alpha0, and compare Delta AIC and Delta BIC against the harmonic K and SdS fits. If the anharmonic models no longer come within 2-3 units of the best harmonic model for any source, the claim that anharmonic corrections are necessary is falsified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The anharmonic model in Eqs. (19) and (20) is a Helmholtz oscillator obtained by truncating the effective potential at third order in δr. Its frequency, Eq. (24), is only the true radial epicyclic frequency if δr remains small enough that quartic and higher terms are negligible and, near the ISCO, if the particle does not cross the ISCO. The algebraic solution itself requires |α1 A0| ≤ α0/2 for the square roots in Eq. (23) to be real, with ϵ = α1 A0. In Table I, the MCMC posteriors return A0 = 5.05 km (GX 17+2, A1), 9.18 and 13.45 km (GX 17+2, A2 and A3), 6.77 km (4U0614+091, A1), and comparable values elsewhere. For the best-fit masses, the orbital radius set by f_U ≈ 0.8–1.2 kHz is only a few gravitational radii, typically 10–25 km; several A0 values are therefore comparable to or larger than r0, and near the ISCO α0 → 0 so |α1 A0|/α0 ≫ 1. At such amplitudes the cubic truncation is invalid, the particle would plunge inside the ISCO, and Eq. (24) is not the oscillation frequency of the underlying geodesic. Since the statistical conclusion that anharmonic corrections are 'physically necessary' rests on the improved MCMC fits for sources like GX 340+0 and 4U0614+091, and those improvements are obtained only by allowing amplitudes that break the perturbative expansion, the central phenomenological claim is not established. A secondary internal issue is the step from Eq. (15) to Eq. (16): setting each bracketed factor to zero is sufficient but not necessary for Eq. (15) to hold, and the true second-order geodesic equations should be verified independently.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes an extension of the relativistic precession model (RPM) for quasi-periodic oscillations (QPOs) by including third-order terms in the Taylor expansion of the effective potential. The authors derive a cubic (Helmholtz-like) radial oscillator, compute the resulting amplitude-dependent radial epicyclic frequency using Jacobi elliptic functions, and apply the model to eight neutron-star QPO sources with MCMC fits, comparing against harmonic Schwarzschild, Schwarzschild–de Sitter, and Kerr models. They conclude that anharmonic corrections are physically necessary but insufficient to fully reconcile the RPM with observations, and that the standard RPM needs a more fundamental revision.","tokens_in":24243,"tokens_out":10088,"duration_ms":107064,"significance":"If the central claim were established, the paper would make a useful contribution: it identifies a simple, geometry-independent feature of the perturbed geodesic equations—the leading anharmonic radial correction is quadratic in the radial displacement—and shows that this term grows relative to the harmonic term near the ISCO. The explicit elliptic-function solution and the full Table I of MCMC results are useful technical resources, and the paper is honest in reporting that anharmonic corrections do not resolve all tensions. The main significance is therefore conditional: the phenomenological evidence for the necessity of anharmonic terms is weakened by the fact that the amplitude A0 is a free parameter fitted to each source, and the fitted amplitudes often violate the small-oscillation assumption on which the expansion is based.","major_comments":[{"comment":"The derivation of the decoupled radial equation is internally inconsistent. The bullet list following Fig. 1 states that β1 δr δθ is about 20% of β0 δθ at every r and therefore “cannot be neglected,” yet the next paragraph concludes that “the above considerations enable us to remove the couplings between the two ODEs” and drops the polar equation. Dropping α2 δθ² from Eq. (16a) makes the radial equation independent of θ, but the polar equation (16b) still contains β1 δr δθ, so the two coordinates are not decoupled. The abstract’s claim that polar oscillations “remain approximately decoupled from radial ones” is contradicted by the paper’s own estimate. The model predictions use only the radial equation, so this might be harmless if the polar motion is treated as a spectator, but that should be stated explicitly and the justification for discarding β1 must be provided.","section":"§III.B, Eqs. (15)–(19)"},{"comment":"The anharmonic frequency shift, Eq. (24) and its small-amplitude expansion (26), is valid only when δr is small enough that the cubic truncation of the effective potential is accurate and the elliptic-solution reality condition is satisfied. The MCMC fits infer A0 values of 5.05, 9.18, and 13.45 km for GX 17+2 (A1–A3), 6.77 km for 4U0614+091 (A1), and several km for other sources. For M ≈ 1.4–2 M⊙, these amplitudes are comparable to the orbital radii inferred from the observed f_U ≈ 0.8–1.2 kHz (roughly 10–25 km), and near the ISCO α0 → 0 while α1 does not, so the controlling ratio α1 A0/α0 can become of order unity or larger. In that regime Eq. (24) is not the radial epicyclic frequency of the underlying geodesic, the particle would leave the perturbative region or cross the ISCO, and the fitted A0 cannot be interpreted as the amplitude of a small oscillation. The paper should impose a validity constraint such as |α1 A0| ≪ α0 at every sampled point, or demonstrate explicitly that the best-fit regions satisfy this condition.","section":"§III.C and Table I"},{"comment":"The statistical evidence does not support the qualitative claim that anharmonic corrections are “physically necessary” as stated in the abstract and Sect. V. The corrections introduce an extra free parameter A0 per source with no independent constraint, so part of the improvement in ln L0 is necessarily a fit. In Table I, the anharmonic models improve the fit decisively only for 4U0614+091 (A1) and mildly for GX 340+0 (A1/A2), while for GX 17+2 they are far worse than the harmonic SdS fit, and for GX 5-1, Sco X1, and 4U1608-52 they barely change ln L0 despite nonzero A0. Given the validity problem identified above, the improved likelihood for 4U0614+091 and GX 340+0 does not establish physical necessity; it shows that a model with an unconstrained amplitude parameter can absorb part of the tension. An independent estimate of A0, or a prior derived from the small-oscillation condition, is needed before this conclusion can be drawn.","section":"§IV and Table I"}],"minor_comments":[{"comment":"The step from Eq. (15) to Eqs. (16) is sufficient but not necessary: setting each bracketed factor to zero guarantees that the sum vanishes, but other solutions of the sum equation exist. The text should clarify that this is an ansatz restricting to solutions for which both brackets vanish separately.","section":"§III.B, Eq. (15)"},{"comment":"There is a typo in the line after Eq. (24): “Ωr = Ωr =√α0” should read “Ωr = √α0”.","section":"§III.C, Eq. (24)"},{"comment":"The heading “Anharmonicity vs radial epyciclic frequency” contains a typo: “epyciclic” should be “epicyclic”.","section":"§III.C heading"},{"comment":"The caption of Fig. 2 contains “ad A3”; this should be “and A3”.","section":"Fig. 2 caption"},{"comment":"The anharmonic corrections are applied only to the Kerr spacetime, while the harmonic SdS results are taken from previous work. The statement that SdS remains a viable replacement is therefore not a model comparison on equal footing; an anharmonic SdS analysis would be needed to make that comparison symmetric.","section":"§IV, Table I"}],"recommendation":"major_revision","confidential_remarks":"The central theoretical observation is interesting, but the statistical and phenomenological claims are currently stronger than the evidence. The authors should be pushed to add explicit validity checks on A0 and to either impose a physically motivated prior or present the anharmonic model as a purely phenomenological fit. A Bayesian evidence comparison, rather than AIC/BIC point estimates, would also strengthen the paper. The decoupling contradiction in §III.B should be fixed before considering publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the theoretical point is real, the statistical work is careful, and the authors are appropriately cautious in their conclusions, but the two main support beams for 'anharmonic corrections are physically necessary' are shaky: the decoupling step is inconsistent with their own Figure 1, and the fitted amplitudes break the perturbative regime in several sources. I would still referee it—the idea deserves a serious look—but the quantitative claim needs rework.\n\nWhat's new: the paper takes the RPM beyond the harmonic approximation by expanding the effective potential to third order, solves the resulting Helmholtz-type radial oscillator with Jacobi elliptic functions, and gets an amplitude-dependent epicyclic frequency. That frequency shift is the kernel. The MCMC comparison over eight sources with three amplitude profiles is a legitimate extension of the authors' earlier harmonic analyses. Credit is due: they separate statistical from physical selection criteria, and they openly report that the anharmonic corrections are insufficient to resolve the preference for SdS in most sources. That honesty matters.\n\nSoft spots, in order of size:\n\n1. The amplitude problem is the load-bearing one. A0 is a free parameter with a broad prior, and the posteriors in Table I give A0 = 9–13 km for GX 17+2 and ~7 km for 4U0614+091, with M around 1.4–2.7 solar masses. For those masses the orbital radius set by f_U ~ 1 kHz is only tens of km, so the displacement is a large fraction of r0, not a small oscillation. The frequency formula in Eq. (24) is derived from a cubic truncation and the elliptic solution requires |α1 A0| ≤ α0/2 for real roots; near the ISCO α0 → 0. The authors never check whether the best-fit parameters satisfy those consistency conditions. Since the improved fits are what motivate 'physically necessary,' this is the central weakness.\n\n2. The decoupling step is internally inconsistent. The text says β1δrδθ is ~20% of β0δθ and cannot be neglected, then removes the coupling anyway. For the radial frequency you only need the radial equation, and α2δθ² is negligible, so the fits may be unaffected; but the abstract's claim that polar oscillations remain approximately decoupled is contradicted by their own Figure 1. This is fixable with a cleaner derivation from the full second-order geodesic equations.\n\n3. Minor: Eq. (15) to Eq. (16) sets each bracketed factor to zero, which is sufficient but not necessary; the actual geodesic equations should be checked independently. Code and data are not provided, which makes the MCMC harder to audit.\n\nWho it's for: researchers working on QPO models, the RPM, or mass/spin estimates of neutron stars from kHz quasi-periodic oscillations. The theoretical part is worth knowing even if the fits are not adopted.\n\nRecommendation: accept for peer review with major revision. The idea is significant enough and the analysis is serious enough that a good referee's time is justified. If the amplitude consistency issue is fixed and the decoupling claim is corrected, this could be a solid contribution.","headline":"A real theoretical extension of the RPM with honest statistics, but the fitted amplitudes violate the perturbative regime and the decoupling step is inconsistent, so the 'necessary' claim isn't yet established.","tokens_in":24722,"tokens_out":5659,"would_cite":true,"duration_ms":60511,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.70.-s","04.70.Bw","04.40.Dg","02.70.Uu"],"model":"deepseek-v4-flash","headline":"The standard harmonic treatment of quasi-periodic oscillations omits a quadratic radial correction that becomes important near the innermost stable orbit.","keywords":["quasi-periodic oscillations","relativistic precession model","anharmonic corrections","radial epicyclic frequency","Kerr spacetime","neutron stars","Helmholtz oscillator","MCMC parameter estimation"],"falsifier":"Integrate the full geodesic equations numerically for the best-fit anharmonic Kerr parameters of one source, compute the resulting $f_L$--$f_U$ curve without the small-amplitude expansion, and compare it with the paper's Eq. (24) against the same data; a disagreement larger than the quoted uncertainties would show the anharmonic frequency formula, not the overall model, is responsible for the residual discrepancy.","tokens_in":23566,"feed_emoji":"🔭","tokens_out":13234,"duration_ms":118231,"temperature":0.7,"pith_summary":"This paper claims that the standard relativistic precession model of quasi-periodic oscillations, which treats orbit perturbations as harmonic oscillators, is missing a term that cannot be discarded in advance: expanding the effective potential to third order produces a radial restoring force quadratic in the displacement, $\\propto \\delta r^2$, in any stationary, axisymmetric, asymptotically flat spacetime. Near the innermost stable circular orbit this quadratic term rivals the harmonic one and shifts the radial epicyclic frequency, while the polar mode stays approximately decoupled. The authors test this by fitting eight neutron star sources first in the harmonic approximation in Schwarzschild, Schwarzschild--de Sitter, and Kerr spacetimes, then with the anharmonic correction in Kerr. They conclude that the anharmonic correction is physically necessary and statistically relevant but still insufficient to fully account for the observed quasi-periodic oscillations, so the relativistic precession model itself needs a more fundamental revision.","feed_headline":"Quasi-periodic oscillation fits demand a quadratic radial correction","feed_subtitle":"Displacement-squared shifts oscillation frequencies near the innermost orbit; Kerr still misses eight neutron stars.","key_machinery":"The load-bearing object is the third-order Taylor expansion of the effective potential for a test particle orbiting in a generic stationary, axisymmetric spacetime. Stability kills the linear terms at the equilibrium orbit, and the surviving third-order terms reduce the radial perturbation to an undamped Helmholtz oscillator, $\\delta r'' + \\alpha_0 \\delta r + \\alpha_1 \\delta r^2 = 0$, with $\\alpha_0$ the harmonic radial epicyclic coefficient and $\\alpha_1$ the quadratic coefficient. The exact solution is $\\delta r(t) = A + B\\,\\mathrm{sn}^2(\\omega t, k)$, where $\\mathrm{sn}$ is the Jacobi elliptic sine, and the angular frequency is built from the complete elliptic integral of the first kind; its small-$\\alpha_1$ expansion gives the amplitude-dependent frequency shift above. This object carries the argument because it converts the geometric statement that the quadratic correction is universal into a concrete shift of the predicted periastron-precession frequency that can be compared with the observed $f_L$--$f_U$ relation.","core_discovery":"Regardless of the spacetime geometry, the leading anharmonic correction to the radial epicyclic oscillation is quadratic in the radial displacement, $\\propto \\delta r^2$, and it is not negligible near the innermost stable circular orbit. Expanding the effective potential about a circular equatorial orbit to third order, the paper shows that the coupling terms between the radial and polar perturbations are subdominant, leaving the radial equation as an undamped Helmholtz oscillator, $\\delta r'' + \\alpha_0 \\delta r + \\alpha_1 \\delta r^2 = 0$. The exact Jacobi-elliptic solution makes the radial epicyclic frequency amplitude dependent, with leading correction $\\Omega_r = \\sqrt{\\alpha_0}\\,\\left(1 - \\frac{5A^2\\alpha_1^2}{12\\alpha_0^2}\\right) + \\mathcal{O}(\\alpha_1^3)$ for small $\\alpha_1$. Fitting the lower and upper quasi-periodic frequencies of eight neutron stars shows that anharmonic Kerr models become viable for some sources, but not for all; in the harmonic limit six of the eight sources prefer Schwarzschild--de Sitter over Kerr, which the authors interpret as a failure of the harmonic approximation and of the Kerr hypothesis within the standard relativistic precession model.","pith_inferences":["Editorial inference: if the quadratic radial correction is genuinely spacetime-independent, the same frequency shift should show up in other epicyclic models of quasi-periodic oscillations, such as those applied to black holes or white dwarfs, giving an independent test of the effective-potential expansion.","Testable extension: the amplitude dependence of the radial frequency turns the oscillation amplitude into a measurable parameter, and the posterior values of several kilometers suggest a fully non-perturbative treatment of the geodesics may be needed before the residual tension is judged real.","Possible connection: comparing the fitted amplitude parameter for each source with observed variability or flux modulation could separate anharmonic orbital effects from non-geodesic disk physics such as pressure gradients or magnetic fields."],"forward_implications":["Any analysis that stops at the harmonic order will misestimate the radial epicyclic frequency for orbits near the innermost stable circular orbit, so masses and spins inferred from such fits carry a systematic error.","In the anharmonic model the radial frequency depends on the oscillation amplitude, so fitting quasi-periodic oscillation data requires choosing an amplitude function rather than a purely geometric frequency relation.","Polar oscillations remain approximately decoupled from radial ones, so the two-frequency structure of the relativistic precession model survives the extension.","The anharmonic correction is not enough: Kerr remains physically excluded for most of the eight sources, and the harmonic preference for Schwarzschild--de Sitter does not disappear, so the model must be modified beyond this term."],"supporting_citations":[{"why":"introduces the relativistic precession model and the identification of the lower quasi-periodic frequency with periastron precession and the upper one with the Keplerian frequency.","marker":"[65]"},{"why":"supplies the innermost stable circular orbit radius and the coefficients $\\alpha_0$ and $\\alpha_1$ for radial epicyclic motion in Kerr spacetime.","marker":"[73]"},{"why":"provides the exact Helmholtz oscillator solution via Jacobi elliptic functions used to derive the amplitude-dependent radial frequency.","marker":"[74]"},{"why":"gives the harmonic Schwarzschild and Schwarzschild--de Sitter best fits whose statistical tensions motivate the anharmonic extension.","marker":"[35]"},{"why":"provides the harmonic Kerr spacetime best fits used as the baseline for comparing the anharmonic models.","marker":"[75]"},{"why":"defines the Akaike and Bayesian information criteria used to rank the harmonic and anharmonic models.","marker":"[76]"}],"fun_headline_variants":["Anharmonic radial term shifts QPO frequencies near ISCO","Quadratic radial correction challenges harmonic QPO model","Relativistic precession needs anharmonic radial kick","QPO fits demand quadratic radial anharmonicity","Kerr and Schwarzschild fail harmonic QPO test"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the radial oscillation amplitude is small and is captured by a single parameter per source; posterior amplitudes that reach several kilometers are comparable to the orbital radius of a neutron star, and at such amplitudes the small-displacement expansion used to derive the frequency shift is outside its validity range.","fun_headline_variants_meta":{"raw":{"variants":["Anharmonic radial term shifts QPO frequencies near ISCO","Quadratic radial correction challenges harmonic QPO model","Relativistic precession needs anharmonic radial kick","QPO fits demand quadratic radial anharmonicity","Kerr and Schwarzschild fail harmonic QPO test"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000164,"raw_usage":{"total_tokens":1316,"prompt_tokens":1083,"completion_tokens":233,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":699,"completion_tokens_details":{"reasoning_tokens":158}},"tokens_in":699,"tokens_out":233,"duration_ms":2729,"temperature":1.0,"reasoning_tokens":158,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T10:19:47.148716+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Integrate the full geodesic equations numerically for the best-fit anharmonic Kerr parameters of one source, compute the resulting $f_L$--$f_U$ curve without the small-amplitude expansion, and compare it with the paper's Eq. (24) against the same data; a disagreement larger than the quoted uncertainties would show the anharmonic frequency formula, not the overall model, is responsible for the residual discrepancy.","supporting_citations":[{"cited_title":"Stella and M","cited_arxiv_id":null,"evidence_quote":"introduces the relativistic precession model and the identification of the lower quasi-periodic frequency with periastron precession and the upper one with the Keplerian frequency."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the innermost stable circular orbit radius and the coefficients $\\alpha_0$ and $\\alpha_1$ for radial epicyclic motion in Kerr spacetime."},{"cited_title":"Hu, Journal of Sound and Vibration 293, 462 (2006), ISSN 0022-460X, URL https://www.sciencedirect","cited_arxiv_id":null,"evidence_quote":"provides the exact Helmholtz oscillator solution via Jacobi elliptic functions used to derive the amplitude-dependent radial frequency."},{"cited_title":"Numerical analysis of quasi-periodic oscillations with spherical spacetimes","cited_arxiv_id":"2212.10186","evidence_quote":"gives the harmonic Schwarzschild and Schwarzschild--de Sitter best fits whose statistical tensions motivate the anharmonic extension."}],"review_version":1}