{"id":"74931422-9552-4810-894f-9a9c1d760b5a","arxiv_id":"2512.14943","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Near a fold point, the slow flow of an unstable oscillator coupled to a nonlinear energy sink obeys the dynamic saddle-node normal form, giving a ε^(1/3)-ε^(2/3) scaling law and an improved mitigation-limit prediction.","lead":"A single-author study shows that the slow oscillations of an unstable mechanical system with a nonlinear energy sink follow a universal fold-point scaling law, with distances shrinking like the 2/3 power of the mass-ratio parameter. The result yields a more accurate prediction of when the energy sink stops suppressing limit-cycle oscillations, validated on an aeroelastic wing model.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Tangent-space center-manifold approximation in §4.1 is the load-bearing weak point: it sets the coefficients a1,a2 (Eq. 47) that determine the quantitative constants in the scaling law (Eqs. 55,57) and the epsilon^{2/3} mitigation-limit correction (Eq. 64).","rationale":"The central claim of the paper is that the slow flow near the fold reduces to a dynamic saddle-node normal form, giving quantitative scaling constants and a mitigation-limit correction. The most load-bearing step is the center manifold reduction. The tangent-space approximation is uncontrolled: the stable fiber coordinate q_b is set to zero, but its curvature influences the quadratic coefficient a1 of the reduced normal form. Since the normal form is rescaled to q_a^2+v, the constants in the scaling law depend on a1/a2. The paper's numerical examples show good agreement, but they do not constitute a parametric check of a1/a2, and the approximation is acknowledged but not quantified. The Airy-branch selection and r∞-arrival assumptions are secondary: the Ai-branch is the one that matches the attracting branch, and the r∞-arrival is a plausible horizontal-jump approximation. Thus the conditional verdict is appropriate, and the concern should be resolved by a higher-order center-manifold computation or a systematic numerical scaling-law fit.","tokens_in":24141,"tokens_out":10033,"duration_ms":93436,"concrete_test":"Compute the center manifold correction ℓ(q_a,u) to second order for the fast subsystem (38) at the left fold, and re-evaluate a1 and a2 from G1(q_a, ℓ(q_a,u), u) + O(u^2). Then recompute K∞ in Eq. (57) and the first-order correction in Eq. (64). If the corrected values differ from the tangent-space Eqs. (47)/(57)/(64) by more than 10% for typical parameters (e.g., μ=0.25, α=5), the central quantitative predictions are not established without the full center-manifold calculation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Sec. 4.1, after diagonalizing the fast subsystem at the left fold, the paper invokes the center manifold theorem but then sets q_b = ℓ(q_a) = 0 ('tangent space approximation'). This is not justified: the true center manifold ℓ(q_a,u) has curvature, and substituting ℓ into G1 contributes to the coefficient of q_a^2 at the same order as the tangent-space term. Equation (46)/(47) therefore gives a1 and a2 only approximately. Since the normal form (48) is rescaled to q_a^2 + v, the constants K0 and K∞ in Eqs. (56)-(57) scale as a1/a2 times a known Airy-zero factor, and the mitigation-limit correction in Eq. (64) inherits K∞. Thus an unquantified error in a1/a2 directly shifts the predicted jump point, the predicted r∞, and the predicted epsilon^{2/3} correction to ρ_U. The exponents 1/3 and 2/3 are robust because they follow from the normal-form scaling, but the paper's quantitative claims (the 'scaling law' with Airy constants and the improved mitigation limit) depend on the constants being correct. The other two listed assumptions — Ai-branch selection and r∞-arrival — are less load-bearing: the Ai-branch is the one matching the attracting branch (Ai'/Ai → -√s), and the r∞-arrival is a reasonable approximation in the horizontal-jump limit.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes a mechanical system with one unstable mode coupled to a cubic nonlinear energy sink (NES). After modal reduction and complexification-averaging, the slow flow is a (2,1)-fast-slow system whose critical manifold has two fold points. The zeroth-order analysis recovers earlier predictions of the mitigation limit. The claimed novelty is a center-manifold reduction near the left fold to the normal form of a dynamic saddle-node bifurcation, an exact Airy-function solution of that normal form, and the resulting scaling law with 1/3 and 2/3 fractional powers of the mass-ratio parameter ε. From this scaling law the paper derives corrected jump/arrival points and an ε^{2/3} correction to the mitigation limit and to the optimal NES damping. The method is tested against direct numerical simulations of both the slow flow and the full aeroelastic wing model.","tokens_in":24518,"tokens_out":7813,"duration_ms":85149,"significance":"If established, the result is significant: it provides an explicit analytic description of the fold escape in a mechanical NES system, going beyond the usual zeroth-order analysis, and it makes falsifiable predictions involving Airy zeros and ε^{1/3}, ε^{2/3} scalings. A notable strength is that no parameter is fitted to the numerical results: the constants a1, a2, f_LF are evaluated from physical inputs, and 1.01879 and 2.33810 are tabulated Airy zeros. The numerical comparisons in Fig. 6 give nontrivial evidence that the predicted mitigation limit tracks the simulations, including for moderately large ε. The work is therefore a useful contribution to the NES literature, provided the three explicit approximations in Section 4 are properly justified or their error is quantified.","major_comments":[{"comment":"The center manifold reduction is invoked, but then the manifold is replaced by its tangent, q_b = ℓ(q_a)=0, 'for sake of simplicity'. This truncation sets the coefficients a1 and a2, which enter the normal form (48) and hence the quantitative constants in Eqs. (55)–(57) and in the mitigation-limit correction (64). The manuscript gives no estimate of the error made by omitting the curvature of ℓ, nor a proof that ℓ contributes only at higher order than the retained q_a^2 and u terms. Since the central quantitative claims depend on a1 and a2, the authors should either justify the truncation by an explicit normal-form calculation or provide a numerical convergence study showing that the omitted terms do not change the predicted constants at the claimed accuracy.","section":"§4.1, Eqs. (46)–(47)"},{"comment":"The text states: 'Assuming that x(−∞) = −√(−y), only the contribution of Ai(s) is kept (this is not proved here).' This is an explicit admission of an omitted proof at a load-bearing point: the choice of Ai instead of Bi determines the numerical constants in Eqs. (55)–(57) and in Eq. (64). The branch selection is in fact standard—Ai(s) is selected by matching to the attracting branch because Ai'/Ai → −√s as s→∞, whereas Bi'/Bi → +√s—but the paper should contain that argument rather than leaving the choice unproved.","section":"§4.2, Eq. (52)"},{"comment":"The new mitigation limit assumes that at the arrival point on the right attracting branch the limit value r∞ has already been reached. This is an uncontrolled approximation: the finite-ε jump has a nonzero duration during which r changes, and no error estimate is given. The numerical agreement in Fig. 5 is encouraging, but Eq. (58) is the basis of the central quantitative prediction (64), so the authors should state the expected order of the error and, ideally, verify it numerically over a range of ε and μ rather than for one parameter set.","section":"§4.3, Eq. (58)"}],"minor_comments":[{"comment":"The caption states 'ξ_h = 4' in one place while Eq. (71) fixes ξ_h = 5 and the text says 'ξ_x = 4 and ξ_φ = 8'. This inconsistency should be corrected.","section":"Fig. 4 caption"},{"comment":"The panel labels are garbled: '((a)ε=0.001, (b),ε=0.005, ε=0.02 and (d)ε=0.1' should read '(a) ε=0.001, (b) ε=0.005, (c) ε=0.02, (d) ε=0.1'.","section":"Fig. 6 caption"},{"comment":"The assumption f(r,s,Δ)=f_LF is introduced without comment. It is probably legitimate at leading order near the fold, but it should be stated explicitly as part of the truncation so that the reader can track all neglected terms.","section":"§4.1, text near Eq. (46)"},{"comment":"Typo: 'right attracting par of M0' should be 'right attracting part of M0'.","section":"§4.3, first paragraph"},{"comment":"The sentence 'using Eq. (40b)' should probably refer to both (40a) and (40b); as written it is slightly confusing.","section":"Eq. (40)"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper's contribution is real: it derives an explicit ε^{2/3} correction to the known zeroth-order mitigation limit for an NES attached to an unstable mechanical system, and it does so without fitting anything to numerics. The reduction of the slow flow near the fold to a dynamic saddle-node normal form is standard in spirit but nontrivial in application, and the resulting formulas (64) and (65), with Airy-zero constants, are new. The validation on both the slow flow and the full aeroelastic wing model is honestly presented, including the breakdown of the asymptotic expansion at larger ε.\n\nThe zeroth-order chain through Section 3 is clean. The scaling exponents 1/3 and 2/3 are robust because they follow from the normal-form scaling, not from any curve fitting. The paper also deserves credit for explicitly flagging its own shortcuts: the unproved Airy-branch selection and the r∞-arrival assumption are stated as such.\n\nThe real soft spot is the one the stress-test flags: Section 4.1 replaces the center manifold by its tangent space, q_b = ℓ(q_a) = 0, and this directly sets the coefficients a1, a2 in Eq. (47). Those coefficients determine K∞ and therefore the quantitative content of the scaling law and the corrections in Eqs. (64) and (65). The text gives no justification that the curvature of the true center manifold is negligible at the same order. For the specific wing system tested, the numerics suggest the approximation works, but that is empirical support, not a derivation. A referee should ask the author to either justify the tangent-space approximation for this class of systems or estimate the induced error, perhaps by computing the next term in the center-manifold expansion.\n\nThe other two assumptions are less load-bearing. The Ai-branch choice is physically the correct one for matching the attracting branch, and the r∞-arrival assumption is reasonable in the horizontal-jump limit. The divergence of the perturbation series at μ = μ* is acknowledged, and the numerical plots show exactly where the approximation degrades.\n\nWho should read this: anyone working on NES design for unstable structures, or on slow-fast systems with fold points. It is not a field reshuffling result, but it is a genuine improvement over the zeroth-order predictions. I would send it to peer review rather than desk reject, with the main request being a sharper treatment of the center-manifold approximation.","headline":"Real analytical result with no fitted parameters: the ε^{2/3} correction to the mitigation limit is credible, but the center-manifold tangent-space approximation leaves the quantitative constants unproved.","tokens_in":25023,"tokens_out":1825,"would_cite":true,"duration_ms":22199,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["34E15","34C23","70K50"],"pacs":[],"model":"deepseek-v4-flash","headline":"Near the left fold point of the critical manifold, the slow flow of a mechanically unstable system with a nonlinear energy sink reduces to the dynamic saddle-node normal form, giving a scaling law with exponents 1/3 and 2/3 and a corrected","keywords":["nonlinear energy sink","slow flow","critical manifold","fold point","dynamic saddle-node bifurcation","Airy function","scaling law","aeroelastic wing"],"falsifier":"Compute the first jump ordinate r_jump from direct numerical integration of the slow flow (38) for several values of epsilon, e.g., 0.001, 0.005, 0.02, and 0.1, at fixed mu and alpha, and plot r_jump - r_LF against epsilon^{2/3}. If the points do not fall on a straight line through the origin with slope 2.33810 a1 (f_LF a2)^{2/3} / a2, or if the intercept is not r_LF, the claimed scaling law fails.","tokens_in":23915,"feed_emoji":"🌀","tokens_out":2424,"duration_ms":26982,"temperature":0.7,"pith_summary":"The paper claims that when a primary structure with one unstable mode is coupled to a nonlinear energy sink, the slow dynamics near a fold point of its critical manifold can be reduced to the normal form of a dynamic saddle-node bifurcation. Solving this normal form yields an exact scaling law: the distance between the critical manifold and the actual trajectory scales with the mass-ratio parameter epsilon to the 1/3 and 2/3 powers, not with the first power. This law predicts where the trajectory leaves the manifold (the jump point) and where it lands (the arrival point), and from it the paper derives a corrected NES mitigation limit that depends on epsilon^{2/3}. If true, this replaces the zeroth-order prediction, which ignores epsilon entirely and is inaccurate even at small epsilon.","feed_headline":"NES mitigation limit gains an epsilon^{2/3} correction","feed_subtitle":"Near the fold point, the slow flow jumps earlier than predicted; Airy-function zeros set the landing point.","key_machinery":"The central object is the reduction of the slow flow to the normal form of the dynamic saddle-node bifurcation. Near a fold point, the center manifold theorem is used to eliminate one fast variable via a tangent-space approximation, leaving a single equation with a slowly varying bifurcation parameter. The exact solution of this normal form is expressed through Airy functions and their zeros, which determine the epsilon^{2/3} location of the departure and arrival points and the epsilon^{2/3} correction to the mitigation limit.","core_discovery":"The central claim is that, in a neighborhood of the left fold point of the critical manifold, the slow flow of an unstable mechanical system coupled to an NES is described by the normal form epsilon q' = q^2 + v, where v is a slowly varying parameter. This dynamic saddle-node normal form has an exact solution in terms of Airy functions: q(v) = epsilon^{1/3} Ai'(-epsilon^{-2/3} v)/Ai(-epsilon^{-2/3} v). Consequently, the jump point and arrival point of a relaxation cycle are set by the first zeros of the Airy function derivative and the Airy function respectively, producing the fractional exponents 1/3 and 2/3. The paper then uses this scaling law to correct the mitigation limit, replacing th","pith_inferences":["The same normal-form reduction should apply at the right fold point, yielding an analogous scaling law with Airy-function constants; this is not developed in the paper but follows directly from the symmetry of the fold geometry.","The exact Airy solution of the normal form is not an asymptotic approximation, so the scaling law may remain valid beyond the perturbative regime, provided the tangent-space and single-branch assumptions hold.","A testable extension is to measure the actual jump ordinates from direct numerical simulations of the slow flow for several epsilon values and fit K_infinity; a consistent epsilon^{2/3} slope would confirm the law and calibrate the correction term."],"forward_implications":["The finite-epsilon jump and arrival points of the slow flow are determined by Airy zeros, so they can be computed once the parameters a1, a2, and f_LF are known, without simulating the full system.","The NES mitigation limit depends on the mass-ratio parameter epsilon through the 2/3 power, meaning even small epsilon produces corrections visible at order epsilon^{2/3} rather than order epsilon.","The optimal NES damping coefficient shifts by a term proportional to epsilon^{2/3}, so NES design can explicitly account for the finite mass of the absorber.","The method applies to a full aeroelastic wing model, where the corrected mitigation limit matches numerical simulations of the full order system better than the zeroth-order prediction for a range of epsilon values.","The zeroth-order mitigation limit is shown to be insufficient: it deviates from numerical results even at epsilon = 0.001, while the new prediction remains accurate up to epsilon = 0.02 and qualitatively acceptable at epsilon = 0.1."],"fun_headline_variants":["Airy functions set the jump point in NES slow flow","NES slow flow jump: Airy zeros set mitigation limit","Slow flow scaling law: 1/3 and 2/3 exponents from Airy","Airy solution predicts earlier NES jump in slow flow","Saddle-node normal form dictates NES slow flow"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The quantitative constants in the scaling law and mitigation-limit correction rely on three unchecked steps: replacing the center manifold by its tangent plane, keeping only the Airy branch that satisfies the initial condition at negative infinity, and assuming the limit value r_infinity has already been reached when the trajectory lands on the right attracting branch.","fun_headline_variants_meta":{"raw":{"variants":["Airy functions set the jump point in NES slow flow","NES slow flow jump: Airy zeros set mitigation limit","Slow flow scaling law: 1/3 and 2/3 exponents from Airy","Airy solution predicts earlier NES jump in slow flow","Saddle-node normal form dictates NES slow flow"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000889,"raw_usage":{"total_tokens":3708,"prompt_tokens":813,"completion_tokens":2895,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":557,"completion_tokens_details":{"reasoning_tokens":2807}},"tokens_in":557,"tokens_out":2895,"duration_ms":18652,"temperature":1.0,"reasoning_tokens":2807,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T15:54:54.508548+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the first jump ordinate r_jump from direct numerical integration of the slow flow (38) for several values of epsilon, e.g., 0.001, 0.005, 0.02, and 0.1, at fixed mu and alpha, and plot r_jump - r_LF against epsilon^{2/3}. If the points do not fall on a straight line through the origin with slope 2.33810 a1 (f_LF a2)^{2/3} / a2, or if the intercept is not r_LF, the claimed scaling law fails.","supporting_citations":[],"review_version":1}