{"id":"0847a1a2-c426-4ff5-a598-a545dd55e326","arxiv_id":"1908.00564","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":8,"one_line_summary":"For large system size, the rescaled beta-FPUT recurrence time is a function only of S = E beta (N+1), linear for small |S| and scaling as |S|^{-1/2} for large |S|, in both positive and negative beta chains.","lead":"Numerical and analytical work maps the first Fermi-Pasta-Ulam-Tsingou recurrence time in the quartic beta-FPUT chain to a single parameter S = E beta (N+1), for both attractive and repulsive nonlinearities. The paper explains the large-S |S|^{-1/2} scaling from solitons in the mKdV continuum limit, with a kink mechanism unique to negative beta.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Large-|S| collapse is not demonstrated at fixed S: data in the power-law fit come mostly from one N, so the 'depends only on S' claim is weakest exactly where it matters.","rationale":"The paper does several things well: the small-|S| shifted-frequency perturbation theory yields a closed expression (19) that matches the numerics across both signs without fitting parameters, and the large-|S| |S|^{-1/2} behavior is independently reproduced by a soliton-velocity argument in the continuum. Those are genuine independent supports. The reason I do not simply accept the headline scaling is that the numerical collapse in the large-|S| regime is the only evidence that the lattice recurrence time, not just the continuum model, is a function of S alone, and the structure of the dataset makes that evidence weaker than it appears. Because recurrences stop for β>0 above Eβ_c ≈ 0.53, smaller N can never reach large S; because of the Appendix B threshold, the same is true for β<0. Thus each large-S point is generated at a different (and generally the largest) N, and the power-law fit cannot separate S-dependence from N-dependence. A fixed-S, varying-N check would settle this directly. The kink-antikink interpretation (Sec. V, Eq. 26) is indeed inferred from heat maps rather than a scattering calculation, but I regard it as a secondary concern: even if that interpretation were wrong, the numerically observed asymmetry in prefactors would remain, and the analytic prefactor is already approximate. The main numerical scaling, if the fixed-S test fails, would not survive. Since the paper is otherwise carefully argued and the required test is straightforward, the appropriate disposition remains conditional pending that test.","tokens_in":20574,"tokens_out":9507,"duration_ms":102350,"concrete_test":"Run fixed-S sweeps with S=60, and also S=20 and S=100 if recurrences survive: integrate the β-FPUT chain for N=127, 255, and 511 with Eβ = S/(N+1), using the same SABA2C integrator, dt, and recurrence-time definition as Sec. III. For β>0, S=60 gives Eβ=0.47 at N=127, below Eβ_c^+=0.53, so all three sizes should show recurrences; for β<0, require N > 4|S|-3 to stay in the large-lattice regime. If T_r at fixed S shifts by more than a few percent when N doubles, the large-|S| power-law fit is a finite-size artifact. Report the tabulated T_r values with error estimates, and additionally list the (N,E,β) points included in the current Figure 1 and 2 fits.","verdict_should_be":"UNCHANGED","load_bearing_attack":"For the central claim as stated — T_r depends only on S for large N, with T_r ~ C_\\pm |S|^{-1/2} at large |S| — the load-bearing point is the numerical demonstration of the large-|S| collapse itself. In Sec. III, the text notes that for β>0 the smaller-N data 'stop at certain values of S' because recurrences cease, and for β<0 Appendix B shows that the lattice is 'large' only when N > 4|S|-3. Hence the points that determine the power-law fits (13) and (15) at S ≳ 30 are attained only at the largest N; the large-S ends of Figures 1 and 2 are essentially a single-N trajectory, and there is no stated test of two different N at the same S in that regime. The fits also have no error bars. If that trajectory is a finite-N curve rather than the true N→∞ locus, the S-only claim in the regime where the |S|^{-1/2} law is asserted is not established. This is a gap in evidence, not an internal inconsistency: the exponent is independently supported by the mKdV soliton estimate, and the authors are candid about the ~30% prefactor mismatch. The reader's kink-antikink concern is real but secondary, since it affects only the explanation of the negative-β prefactor, not the numerical fact of the difference.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the first FPUT recurrence time in the β-FPUT chain for both positive and negative β. The central numerical claim is that the rescaled recurrence time T_r = t_r/(N+1)^3 depends, for large N, only on the combination S = Eβ(N+1). For small |S| the authors report a linear behavior with nearly β-independent slope, and for large |S| they report T_r ∝ |S|^{-1/2} with different multiplicative constants for β>0 and β<0. The paper also derives a closed-form small-|S| expression from shifted-frequency perturbation theory (Eq. (19)), studies the mKdV continuum limit numerically, interprets the β<0 case in terms of soliton-kink interactions, and derives the large-|S| power law from soliton velocities (Eqs. (33) and (36)). Additional sections quantify the parameter range in which recurrences exist and discuss the different energy-mixing behavior for the two signs of β.","tokens_in":20904,"tokens_out":5892,"duration_ms":64540,"significance":"If the central scaling claim holds, this is a valuable contribution: it provides a clean one-parameter description of β-FPUT recurrence times analogous to the R-scaling in the α-FPUT chain, and it extends the soliton-based recurrence picture to both signs of β, including the less-studied negative-β case. The paper's analytic small-|S| expression is derived from established perturbation theory rather than fitted to the recurrence data, and the large-|S| exponent is independently obtained from the mKdV soliton spectrum, which is a genuine strength. The authors are also commendably explicit about the ~30% discrepancy between their analytic prefactors and the numerical fits, and about the inferential character of the kink-antikink identification. The main weakness is that the numerical demonstration of the S-only collapse is incomplete precisely in the large-|S| regime where the power law is asserted; this is a load-bearing point for the paper's central claim, though not an internal inconsistency.","major_comments":[{"comment":"The central claim that T_r depends only on S for large N is not fully demonstrated in the large-|S| regime where the power-law fits are made. By the paper's own criteria, the β<0 data in the fitted range require the lattice to be 'large,' which Appendix B states requires N > 4|S|-3, and the β>0 data terminate at each N because recurrences cease. Consequently, the high-|S| points that determine fits (13) and (15) are reached only by the largest N, and no two-N comparison at fixed S is reported in that range; the fits are also presented without error bars. Please provide a direct collapse test or equivalent evidence in the fitted regime (e.g., T_r(S) for several N at common S values with S ≳ 30 and S ≲ -12), state the N values and the number of points used in each fit, and report residuals or confidence intervals. This gap concerns evidence rather than internal consistency, and the exponent has independent analytic support, but the 'depends only on S' claim needs this support where it is asserted.","section":"Sec. III, Eqs. (13) and (15), Figs. 1 and 2"},{"comment":"The different β<0 prefactor is explained through the assertion that the negative-β mKdV background becomes a traveling kink-antikink, which leads to the recurrence spacing L rather than 2L in Eq. (26). This identification is inferred from the heat maps in Fig. 4 ('given the numerical results, we argue') and is not supported by a quantitative scattering or inverse-scattering calculation. Because Eq. (26) and hence Eq. (36) depend on this spacing choice, the analytic derivation of the β<0 prefactor is conditional on the kink-antikink identification. Please add a quantitative diagnostic (for example, a measurement of the asymptotic background amplitude or of soliton phase shifts in the mKdV simulations), or clearly present Eq. (36) as an ansatz motivated by the observed spacing rather than as a derivation.","section":"Secs. V and VI, Fig. 4, Eq. (26)"}],"minor_comments":[{"comment":"The thresholds S ≳ 30 and S ≲ -12 are stated without an explicit criterion for what counts as an 'accurate' fit; please specify how these ranges were chosen.","section":"Sec. III, Eqs. (12)-(15)"},{"comment":"The last term in the displayed expression is typographically ambiguous: it should be written as a fraction, e.g., 16π^5/(24π^2 - 27S), rather than in a form that reads as two separate terms.","section":"Eq. (18)"},{"comment":"The statement that the boundary conditions are 'the same as if the shift were a real number' is terse for the PT-symmetric, non-Hermitian oscillator; a brief justification or reference would help the reader.","section":"Sec. VI, Eq. (29)"},{"comment":"The parameter φ∞ is used before it is defined; please define it explicitly as the background amplitude and state the sign restrictions that follow from the soliton existence conditions.","section":"Appendix D, Eq. (D1)"},{"comment":"The 95% confidence intervals for Eβ_c^+ and Eβ_c^- are reported without describing the fitting or uncertainty procedure; please add a sentence explaining how these intervals were obtained.","section":"Sec. VII"},{"comment":"The captions would be more self-contained if they specified the color scale values and the meaning of white regions in the heat maps.","section":"Figs. 4 and 7"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The headline: this is the first systematic study of FPUT recurrence times in the beta-chain for both signs of beta, and the central scaling claim—rescaled T_r depending only on S = E beta (N+1), with a linear small-S regime and a |S|^{-1/2} large-S law—is very likely correct. The paper is worth a serious referee, but the evidence for the large-S branch is thinner than the prose suggests.\n\nWhat's genuinely new: the S-only collapse for both signs of beta, the closed-form small-S expression from shifted-frequency perturbation theory that matches numerics beautifully, and the mKdV soliton interpretation that explains the factor-of-two prefactor difference between beta>0 and beta<0. The analytic derivation of the |S|^{-1/2} exponent from soliton velocities is independent of the numerics and lands on the right scaling, even though the prefactors overshoot by 27–33 percent. The authors are honest about that mismatch and about the approximations behind it.\n\nSoft spots, in proportion. The stress-test note is on target: the large-|S| power-law fits in Figs. 1 and 2 are effectively single-N data. For beta>0, smaller N stop because recurrences vanish; for beta<0, the lattice is \"large\" only for N > 4|S|-3. So no two different N are shown to agree at the same large S. That means the \"depends only on S\" claim is not actually demonstrated in the regime where the headline scaling is asserted. The exponent is rescued by the soliton calculation, but the prefactor and the S-only universality at large S rest on a weaker numerical base than the small-S collapse. This is a fixable gap—run two larger N at the same S values—but as written it is a real hole. Also, the kink-antikink interpretation for beta<0 rests on heat-map inspection, not on a scattering calculation; that's acknowledged but it leaves the beta<0 prefactor explanation heuristic. Minor: no error bars on the fits and no deposited code or data, which would help a lot.\n\nVerdict: send to review. The core result is important enough and the analytic machinery is solid enough to deserve referee time. The referee should push for the fixed-S multi-N test and for data availability, but the paper's central message will survive. I'd bring it to a reading group and would cite it if I worked on lattices or soliton recurrences.","headline":"Convincing scaling study of beta-FPUT recurrence times with real analytic support, but the large-|S| collapse needs a fixed-S multi-N test before it is fully established.","tokens_in":21447,"tokens_out":2686,"would_cite":true,"duration_ms":21748,"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":"For the quartic (β) FPUT chain, the rescaled first recurrence time is controlled for large N by the single parameter S = E β (N+1), matching numerics, perturbation theory, and mKdV soliton theory for both signs of β.","keywords":["Fermi-Pasta-Ulam-Tsingou recurrence","beta-FPUT chain","quartic anharmonic chain","mKdV solitons","kink-antikink","recurrence time scaling","shifted-frequency perturbation theory","metastable state"],"falsifier":"Track one solitary wave in a numerically integrated β < 0 mKdV solution with cosine initial data: the paper's mechanism predicts the wave changes from soliton to antisoliton when it crosses the background and that the two overlap points at the recurrence are separated by L, not 2L; observing a different crossing behavior or spacing would falsify the prefactor explanation even if the $|S|^{-1/2}$ law survives.","tokens_in":20304,"feed_emoji":"🔄","tokens_out":6479,"duration_ms":66787,"temperature":0.7,"pith_summary":"This paper asks what controls the first Fermi-Pasta-Ulam-Tsingou (FPUT) recurrence in the quartic (β) chain, whose nonlinear potential is proportional to the fourth power of nearest-neighbor displacement. It claims that for large system size N, the rescaled recurrence time $T_r = t_r/(N+1)^3$ is a function of the single parameter $S \\equiv E\\beta(N+1)$ alone, with no separate dependence on energy, coupling, or N. In the nearly linear regime (small $|S|$), the paper derives a closed-form expression for $T_r$ from shifted-frequency perturbation theory, and in the highly nonlinear regime (large $|S|$), it shows numerically and analytically that $T_r \\propto |S|^{-1/2}$ for both signs of β, with different prefactors. A sympathetic reader would care because this collapses an apparently high-dimensional dynamical problem onto one universal curve and ties lattice recurrences to soliton and kink physics in the modified Korteweg-de Vries (mKdV) equation.","feed_headline":"One parameter predicts the quartic FPUT recurrence time","feed_subtitle":"Rescaled recurrence times collapse onto one curve set by S = E β (N+1), for both signs of β.","key_machinery":"The load-bearing quantity is the dimensionless parameter $S \\equiv E\\beta(N+1)$, which combines the total energy, the quartic coupling sign and strength, and the number of masses into one control parameter. The analytic work rests on two tools: shifted-frequency perturbation theory, in which nonlinear normal-mode frequencies $\\Omega_k$ are computed to second order and the recurrence time is $t_r = 2\\pi/(3\\Omega_1 - \\Omega_3)$, and the continuum mKdV description, in which the first-mode initial condition breaks into soliton-antisoliton pairs; the recurrence time is estimated from the velocity spacing between consecutive solitons through $\\tau_r = \\frac{b+3}{2}\\,\\frac{L}{\\Delta v}$, with $b = \\mathrm{sgn}(\\beta)$. For β < 0, kink solutions exist and solitons flip between soliton and antisoliton form when crossing the background, which changes the effective spatial separation appearing in the velocity-spacing estimate and produces the different prefactor.","core_discovery":"The central discovery is that the first FPUT recurrence time in the β chain is, for large N, fully described by $S \\equiv E\\beta(N+1)$: the rescaled time $T_r = t_r/(N+1)^3$ is a single function of S for β > 0 and β < 0. Numerically, $T_r$ is linear in S for small $|S|$ with nearly identical positive slope for both signs, reaches a maximum near $S \\approx 4.2$, and for large $|S|$ follows $T_r \\approx 0.5862\\,S^{-1/2}$ for β > 0 and $T_r \\approx 0.3078\\,|S|^{-1/2}$ for β < 0. The paper reproduces the small-$|S|$ branch with a closed formula from shifted-frequency perturbation theory and the large-$|S|$ branch from soliton velocities in the mKdV continuum limit; the smaller β < 0 prefactor is attributed to soliton-kink interactions that make recurrences occur sooner. The paper also finds that recurrences cease to form beyond critical values $E\\beta_c^+ \\approx 0.53$ and $E\\beta_c^- \\approx -2.43$ for large N, which it connects to the breakdown of the metastable state.","pith_inferences":["The S-collapse suggests the same one-parameter scaling may organize other diagnostics of the β chain, such as energy mixing or higher-order recurrences; the paper only checks the first recurrence.","Because the paper connects β < 0 dynamics to a one-dimensional Bose gas in the quantum rotor regime, the predicted difference in recurrence prefactors could be tested in ultracold bosons in optical lattices at energies below the blow-up threshold.","The PT-symmetric harmonic-oscillator structure that appears in the mKdV Schrödinger problem for β > 0 could be used to compute phase-shift corrections; the paper leaves those corrections for future work.","The existence of recurrences at $|E\\beta|$ values almost ten times the energy above which blow-up is possible for β < 0 implies a metastable regime worth studying separately from the blow-up instability; this separation is not explored beyond the empirical threshold."],"forward_implications":["For large N, any experiment or simulation that fixes S should see the same first recurrence time whether the quartic coupling is attractive or repulsive, provided blow-up is avoided; the sign of β enters only through the multiplicative constant at large $|S|$.","The closed-form small-$|S|$ expression lets one predict recurrence times near the harmonic limit without integrating the lattice, using only S and N.","The $|S|^{-1/2}$ scaling in the highly nonlinear regime links lattice recurrence to soliton and kink dynamics, so recurrence-time data can be used to probe the soliton content of the mKdV initial condition.","Beyond a critical value of $E\\beta$ that is independent of N for large N, recurrences stop forming; equivalently, in the thermodynamic limit with E growing extensively, recurrences will not occur at all."],"supporting_citations":[{"why":"Supplies the shifted-frequency perturbation treatment and the expression $t_r = 2\\pi/(3\\Omega_1 - \\Omega_3)$ that this paper extends to a closed form in S.","marker":"[7]"},{"why":"Establishes the soliton-based explanation of FPUT recurrence in the KdV equation, the template for the mKdV interpretation used here.","marker":"[8]"},{"why":"Provides the scaling of the recurrence time in the cubic α-FPUT chain and motivates the $(N+1)^3$ rescaling and one-parameter reduction.","marker":"[12]"},{"why":"Shows a limiting four-mode Hamiltonian depends on $6\\beta E(N+1)/\\pi^2$, the precedent for introducing S.","marker":"[18]"},{"why":"Uses soliton velocities to estimate recurrence times, the approach generalized here to the mKdV soliton-velocity calculation.","marker":"[35]"},{"why":"Connects metastable-state breakdown to a critical $E\\beta$ for β > 0, which the paper uses to interpret the recurrence-formation threshold.","marker":"[54]"}],"fun_headline_variants":["One parameter collapses FPUT recurrence times","Single S-curve sets FPUT recurrence time","β± share one recurrence-time law","Solitons set FPUT recurrence time in β chain"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the β < 0 mKdV background is a traveling kink-antikink and that solitary waves flip their form when crossing it; this identification is inferred from spacetime heat maps rather than from a quantitative scattering calculation, and it is what sets the different β < 0 prefactor.","fun_headline_variants_meta":{"raw":{"variants":["One parameter collapses FPUT recurrence times","Single S-curve sets FPUT recurrence time","β± share one recurrence-time law","Solitons set FPUT recurrence time in β chain"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000261,"raw_usage":{"total_tokens":1736,"prompt_tokens":1228,"completion_tokens":508,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":844,"completion_tokens_details":{"reasoning_tokens":451}},"tokens_in":844,"tokens_out":508,"duration_ms":4947,"temperature":1.0,"reasoning_tokens":451,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:48:11.557331+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Track one solitary wave in a numerically integrated β < 0 mKdV solution with cosine initial data: the paper's mechanism predicts the wave changes from soliton to antisoliton when it crosses the background and that the two overlap points at the recurrence are separated by L, not 2L; observing a different crossing behavior or spacing would falsify the prefactor explanation even if the $|S|^{-1/2}$ law survives.","supporting_citations":[{"cited_title":"small” and “large","cited_arxiv_id":null,"evidence_quote":"Supplies the shifted-frequency perturbation treatment and the expression $t_r = 2\\pi/(3\\Omega_1 - \\Omega_3)$ that this paper extends to a closed form in S."},{"cited_title":"(C25) We note that in equations (C9-C25), all of the nonlinear frequencies are calculated to first order in β","cited_arxiv_id":null,"evidence_quote":"Establishes the soliton-based explanation of FPUT recurrence in the KdV equation, the template for the mKdV interpretation used here."},{"cited_title":"Gallavotti, ed., The Fermi-Pasta-Ulam Problem: A Status Report","cited_arxiv_id":null,"evidence_quote":"Provides the scaling of the recurrence time in the cubic α-FPUT chain and motivates the $(N+1)^3$ rescaling and one-parameter reduction."},{"cited_title":"Phenomena Associated with the Oscil- lations of a Nonlinear Model String,","cited_arxiv_id":null,"evidence_quote":"Shows a limiting four-mode Hamiltonian depends on $6\\beta E(N+1)/\\pi^2$, the precedent for introducing S."},{"cited_title":"Lax pair","cited_arxiv_id":null,"evidence_quote":"Uses soliton velocities to estimate recurrence times, the approach generalized here to the mKdV soliton-velocity calculation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Connects metastable-state breakdown to a critical $E\\beta$ for β > 0, which the paper uses to interpret the recurrence-formation threshold."}],"review_version":1}