{"id":"eb550826-6b7d-45c2-9faf-98aa5d91957d","arxiv_id":"1908.00721","paper_version":3,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"A multi-degree-of-freedom Melnikov theory predicts which conservative backbone curves persist as forced-damped responses, when they become maximal-response ridges, and when isolas are born.","lead":"This paper derives an analytic Melnikov criterion that tells which periodic motions of an undamped mechanical system survive when small damping and periodic forcing are added, and how they form the peaks and isolated branches of the forced response. It justifies the energy-balance and phase-lag methods used in vibration testing and verifies the predictions on a six-degree-of-freedom example.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified.","rationale":"The paper's central claim is a local persistence theorem, and every step needed for that claim is either proved in the appendix or explicitly assumed. The simple-zero condition in Theorem 3.1 is the standard Melnikov nondegeneracy condition, and the proof via the scalar bifurcation function is sound: both the displacement map and the constraint equation are reduced by the implicit function theorem, and the leading-order term is the stated Melnikov function. The bounded-away-from-zero case is equally standard. The extensions in Theorem 3.2 and Proposition 4.2 follow from the same reduction plus standard singularity theory, and the phase-lag criterion is derived cleanly from the quadratic-zero condition. The numerical examples provide independent support for the ridge and isola predictions, including the otherwise difficult-to-find isola birth. The reader's weakest assumption was smooth, small-epsilon persistence; this is not a hidden flaw because the paper explicitly discusses it in Remark A.1 and in the concluding limitations, and it is the usual scope of Melnikov-type arguments. The only artifact I noticed is the duplicated text block near the start of Section 5 in the provided manuscript, which has no bearing on the mathematical content. I therefore see no reason to change the reader's ACCEPT verdict.","tokens_in":40006,"tokens_out":13750,"duration_ms":162737,"concrete_test":"As a verification step still worth running: repeat the nonlinear-damping example of Section 5.2 at smaller perturbation values, e.g., epsilon=0.02 and epsilon=0.01, and confirm that the numerically continued saddle-node/isola curve (black line in Fig. 8c) approaches the analytically predicted ridge R1 (green line) with the expected O(epsilon) gap; if the gap does not shrink linearly, the ridge prediction would need to be re-examined.","verdict_should_be":"UNCHANGED","load_bearing_attack":"I read the manuscript as establishing a Melnikov-type persistence and bifurcation criterion for m-normal conservative periodic orbits under small time-periodic perturbations, with O(epsilon)-close smooth continuation guaranteed by Theorem 3.1 and the implicit-function reduction in Theorem A.2. The argument is internally consistent: the reduction to a scalar bifurcation function, the L-independence of the leading-order term, and the energy-balance interpretation all check out, and the numerical continuation at epsilon=0.05 and 0.1 reproduces the predicted ridges and isola birth. The principal limitation is the one the authors themselves flag in Remark A.1: Theorem 3.1 proves smooth persistence only, so periodic orbits that are O(epsilon)-close but not smoothly connected to Z are outside the criterion, and near branch points where 1-normality fails the result does not apply. This is a genuine but explicitly scoped boundary of the central claim, not an internal inconsistency. I do not find a load-bearing gap that would invalidate the accepted verdict.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper develops a Melnikov-type criterion for the persistence and bifurcation of periodic orbits of a conservative mechanical system under small damping and time-periodic forcing. The main mathematical result (Theorem 3.1) reduces the persistence problem for an m-normal periodic orbit to the zeros of a scalar Melnikov function M_{m:l}(s); a simple zero guarantees two O(ε)-close smooth continuations, while a quadratic zero leads to saddle-node, isola birth, or simple bifurcation phenomena (Theorem 3.2). For monoharmonic forcing with arbitrary dissipation, the authors compute M explicitly and derive ridge curves Γ_l(λ) that locate maximal response and isola births from the conservative backbone curve alone, as well as a generalized phase-lag quadrature criterion. The theoretical predictions are validated on a six-degree-of-freedom chain with linear and nonlinear damping, including the birth and merging of an isola. The central limitation, explicitly stated in Remark A.1, is that Theorem 3.1 addresses smooth persistence only; periodic orbits that are O(ε)-close but not smoothly connected to the conservative orbit are outside the scope.","tokens_in":129,"tokens_out":7277,"duration_ms":374346,"significance":"The results are significant for nonlinear structural dynamics and for perturbation theory. They extend subharmonic Melnikov theory beyond planar/integrable systems, exploiting one-parameter families of periodic orbits instead of integrability. They also provide rigorous justification for energy balance and force-appropriation methods under broader conditions than before. The predictions are falsifiable: ridge curves and bifurcation thresholds are computed from the unperturbed conservative limit, with no parameters fitted to forced-response data. The numerical example supports the claimed accuracy even at ε=0.1. The main caveat—smooth persistence only—is a genuine but clearly scoped boundary of the method.","major_comments":[],"minor_comments":[{"comment":"The statement that DΓ_l>0 yields a 'maximal response' is derived from the existence of two orbits for λ<λ0 and none for λ>λ0; it would be clearer to state explicitly that this is a fold with respect to λ and to point out that in the examples λ plays the role of the frequency ω, which is what makes the ridge coincide with the peaks in Figure 7(a).","section":"§4.2 (Proposition 4.2, Eq. (25))"},{"comment":"The values α=0.2481, β=-1.085, γ=0.8314 are selected to break the monotonic trend of the resistance; it might be worth stressing that this choice does not introduce any fitted parameter in the ridge prediction, since Γ_1(λ) is computed solely from the conservative limit.","section":"§5.2"},{"comment":"The paper explicitly limits the result to smooth persistence; I suggest cross-referencing Remark A.1 in the main-text discussion around Fig. 1 so that readers do not infer that all O(ε)-close forced responses are captured.","section":"§3.1 (Theorem 3.1, Remark A.1)"},{"comment":"The duplicated paragraph and several typographical errors (e.g., 'detachement' in the Figure 4 caption, and the stray 'some' in the draft text) should be corrected in a final version.","section":"§5 (Eq. (34) and Figure 6)"}],"recommendation":"minor_revision","confidential_remarks":"This is a solid paper. The central results are internally consistent, the numerical tests support the claims, and the limitations are honestly stated. The only concerns are presentational, so I recommend minor revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper deserves a serious referee. The core result is real and the authors mostly respect its boundaries.\n\nWhat is actually new: they reduce the persistence of a family of conservative periodic orbits under small damping and periodic forcing to zeros of a scalar Melnikov-type function, without assuming the unperturbed limit is integrable or low-dimensional. That is a genuine step beyond Chicone/Rhouma–Chicone, who handled similar reductions but with exact resonance and more restrictive setups. The ridge criterion (Proposition 4.2) that turns energy balance into a predictive curve in the forcing-amplitude/frequency plane is useful, and the explicit formula for monoharmonic forcing is simple enough for practitioners to use. The six-DOF example is the right kind of evidence: the ridges are computed from the conservative limit only, then matched to independently continued forced responses, including isola birth and merging. That is not a fitted curve; it is a prediction.\n\nThe paper is also honest about the main limitation. Remark A.1 concedes what Theorem 3.1 actually proves: smooth persistence of orbits O(ε)-close to the unperturbed orbit. Orbits that are close but only continuously connected, or not connected at all, are outside the criterion. That is a real boundary, but the authors state it plainly and do not oversell. The numerical examples use ε = 0.05 and 0.1, which are small but not infinitesimal, so the match gives the theory some breathing room.\n\nSoft spots, in proportion: Proposition 4.2’s language about “maximal response with respect to λ” is a bit compressed—it is really a fold in the continuation parameter λ, and calling it a maximal response in the amplitude-frequency plot relies on the backbone being parametrizable by both frequency and an amplitude measure. That is fine for the examples but worth a careful referee checking the exact geometric claim. The nonlinear damping coefficients in Section 5.2 are hand-picked to create an isola; that is acceptable for a demonstration, though it does not test the theory’s predictive power in the way the linear-damping sweep does. There is also a visible duplicated-text artifact in the arXiv version around Section 5 (the same paragraph appears twice with slightly different formatting). It has no scientific impact but should be cleaned up before publication.\n\nOverall, the math is internally consistent, the assumptions are stated, and the numerical evidence supports the claims as scoped. The paper will be useful to nonlinear structural dynamicists who want a fast analytic route from conservative backbone curves to forced-response peaks and isolas, and to anyone who wants a rigorous justification for energy-balance or phase-lag heuristics. I would cite it and would bring it to a reading group.\n\nRecommendation: send to peer review. Expect careful checking of Proposition 4.2 and the nondegeneracy conditions, but this is not a desk-reject candidate.\n\nBest,\n[You]","headline":"Solid, honest generalization of Melnikov's idea to multi-DOF conservative orbit families, with real numerical verification and clearly scoped limits; worth refereeing seriously.","tokens_in":40646,"tokens_out":855,"would_cite":true,"duration_ms":11678,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["34C25","37C27","70K40"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper reduces the survival of conservative periodic orbits under small damping and forcing to the zeros of a single Melnikov-type integral, and uses those zeros to predict resonance maxima and isola births.","keywords":["Melnikov function","backbone curves","periodic orbit persistence","forced-damped response","nonlinear normal modes","isola bifurcation","phase-lag quadrature","multi-degree-of-freedom systems"],"falsifier":"Return to the six-mass example with $\\alpha=0.2481$, $\\beta=-1.085$, $\\gamma=0.8314$, $\\varepsilon=0.1$, and use the first-mode conservative family. The paper predicts an isola birth near $e\\approx 0.4$ and a reconnection to the main branch near $e\\approx 1$, both $O(\\varepsilon)$-close to the ridge $e=\\Gamma_1(\\bar\\omega)$. Continuing forced periodic orbits directly and tracking saddle-node curves: if the isola does not appear, or appears at amplitudes outside the predicted interval, or if any forced branch bifurcates continuously from a conservative orbit for which $W_{1:1}<R$, the criterion is contradicted.","tokens_in":39801,"feed_emoji":"⚙️","tokens_out":9139,"duration_ms":87552,"temperature":0.7,"pith_summary":"The paper asks which periodic orbits of an undamped, unforced mechanical system (the conservative backbone curves) actually show up as resonance peaks once small damping and periodic forcing are switched on. It answers with a Melnikov-type function: an integral over one period of the unforced orbit that measures the energy balance of the perturbation. If this function has a simple zero, the conservative orbit persists as a nearby forced-damped periodic response; if it stays away from zero, no smooth response is born from that orbit. For single-frequency forcing with arbitrary dissipation, the criterion collapses to comparing the work done by the forcing with the energy dissipated by damping: two response branches appear when work beats dissipation, none when it does not, and a saddle-node at equality. The same function's quadratic zeros locate ridges where the response is maximal or where isolated response branches (isolas) are born.","feed_headline":"A Melnikov integral decides which backbone curves survive forcing","feed_subtitle":"Its zeros mark exactly where conservative orbits reappear as damped forced responses; its folds mark isolas.","key_machinery":"The central object is the subharmonic Melnikov function $M_{m:l}(s)$ of Eq. (9), an integral of the perturbation against the energy gradient along the unperturbed orbit. It is the leading-order term of the one-cycle energy balance; its zeros are exactly the phase shifts $s$ at which the perturbation does zero net work. Because it is scalar, the implicit function theorem reduces an $(n+1)$-dimensional persistence problem to checking whether $M_{m:l}$ crosses zero transversally. The second necessary ingredient is the $m$-normal periodic orbit, a nondegeneracy condition (geometric multiplicity of the $+1$ Floquet multiplier at most two, with a tangency condition in the multiplicity-two case) that guarantees the conservative orbit belongs to a smooth one-parameter family; this replaces the integrability assumptions of classical Melnikov theory.","core_discovery":"The central claim is that the fate of a one-parameter family of conservative periodic orbits under small damping and periodic forcing is governed by the zeros of the scalar Melnikov function $$M_{m:l}(s)=\\$int_0^{{m\\tau}}$\\langle \\nabla H(x_0(t+s;p)), g(x_0(t+s;p),t;\\tau m/l,0)\\rangle\\,dt.$$ Theorem 3.1 asserts that if $M_{m:l}$ has a simple zero at $s_0$, the $m$-normal orbit $Z$ continues smoothly, with initial condition and period $O(\\varepsilon)$-close to $Z$ and $m\\tau$; if $M_{m:l}$ is bounded away from zero, there is no such smooth continuation. Theorem 3.2 classifies quadratic zeros as saddle-node bifurcations, isola births, or simple bifurcations. Under monoharmonic forcing and arbitrary dissipation, the function factors as $M_{1:l}=W_{1:l}(e)\\cos(l\\omega s-\\alpha_{l,e})-R$, and Proposition 4.2 shows that the ridge $e=\\Gamma_l(\\lambda)=R(\\lambda)/A_{l,e}(\\lambda)$ is an $O(\\varepsilon)$-close locus of maximal or minimal forced responses when $\\Gamma_l$ has nonzero slope, and of isola births or simple bifurcations where the slope vanishes. The paper argues this justifies the energy principle and a generalized phase-lag quadrature criterion, and demonstrates the predictions on a six-degree-of-freedom chain.","pith_inferences":["Extension: because the ridge uses only conservative-limit data, the same computation could screen many candidate modes and damping laws for isola-prone zones before any forced-damped continuation is run.","Extension: a natural next step would be to derive a generalized ridge for multi-harmonic forcing, where each harmonic should contribute its own phase and amplitude term; the paper only treats monoharmonic forcing.","Extension: near the branch point where the first mode stops being $1$-normal, the Melnikov criterion is silent, so a higher-dimensional bifurcation function would be needed to cover that amplitude regime; this is a concrete boundary of the present result."],"forward_implications":["For monoharmonic forcing, each conservative orbit produces two forced responses exactly when the forcing's work exceeds the damping's resistance, one response at equality, and none below; this gives a closed-form amplitude threshold.","Ridges computed from the conservative limit alone locate the peaks of frequency-response diagrams to within $O(\\varepsilon)$, so numerical continuation of the forced-damped system is not needed to find them.","Quadratic zeros of the Melnikov function are analytic early-warning signatures of isolas, which are otherwise hard to find by continuation.","The phase-lag quadrature criterion is valid for asynchronous, multi-harmonic periodic motions with arbitrary smooth damping, provided the phase lag is measured in co-location with the forcing.","For the monoharmonic forcing class considered, the criterion automatically rules out superharmonic and ultrasubharmonic resonances, while subharmonic resonances are captured through the $1:l$ case."],"supporting_citations":[{"why":"Supplies the reduction of periodic-orbit continuation to a scalar bifurcation function, which the paper extends to non-integrable multi-degree-of-freedom systems.","marker":"[57]"},{"why":"Defines $m$-normal periodic orbits and guarantees a smooth one-parameter family of such orbits.","marker":"[60]"},{"why":"Establishes continuation of periodic orbits in conservative systems, backing the normality and family-persistence assumptions.","marker":"[61]"},{"why":"Introduces the phase-lag quadrature criterion that Proposition 4.3 justifies under weaker assumptions.","marker":"[34]"},{"why":"Formulates the energy principle for locating forced responses that the Melnikov function makes rigorous.","marker":"[24]"},{"why":"Classifies saddle-node, isola-birth, and simple-bifurcation singularities used in Theorem 3.2.","marker":"[63]"},{"why":"Supplies the numerical continuation toolbox used to compute the example backbone curves and frequency responses.","marker":"[33]"},{"why":"Introduces the original scalar Melnikov integral for planar oscillators that this work generalizes.","marker":"[47]"}],"fun_headline_variants":["Melnikov zeros pick which backbones survive forcing","A scalar integral decides backbone fate under damping and drive","Melnikov rule: zeros mark surviving backbones, folds birth isolas","Zeros of Melnikov function predict forced responses near conservative orbits","Melnikov's integral: the gatekeeper for forced backbone curves"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claim presupposes that the forced-damped response is an $O(\\varepsilon)$ smooth continuation of a nondegenerate ($m$-normal) conservative periodic orbit, so near branch points where normality fails, or for perturbations that are not small, the Melnikov criterion does not apply.","fun_headline_variants_meta":{"raw":{"variants":["Melnikov zeros pick which backbones survive forcing","A scalar integral decides backbone fate under damping and drive","Melnikov rule: zeros mark surviving backbones, folds birth isolas","Zeros of Melnikov function predict forced responses near conservative orbits","Melnikov's integral: the gatekeeper for forced backbone curves"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00069,"raw_usage":{"total_tokens":3154,"prompt_tokens":1002,"completion_tokens":2152,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":618,"completion_tokens_details":{"reasoning_tokens":2066}},"tokens_in":618,"tokens_out":2152,"duration_ms":16985,"temperature":1.0,"reasoning_tokens":2066,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:34:35.697597+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Return to the six-mass example with $\\alpha=0.2481$, $\\beta=-1.085$, $\\gamma=0.8314$, $\\varepsilon=0.1$, and use the first-mode conservative family. The paper predicts an isola birth near $e\\approx 0.4$ and a reconnection to the main branch near $e\\approx 1$, both $O(\\varepsilon)$-close to the ridge $e=\\Gamma_1(\\bar\\omega)$. Continuing forced periodic orbits directly and tracking saddle-node curves: if the isola does not appear, or appears at amplitudes outside the predicted interval, or if any forced branch bifurcates continuously from a conservative orbit for which $W_{1:1}<R$, the criterion is contradicted.","supporting_citations":[{"cited_title":"Rhouma and C","cited_arxiv_id":null,"evidence_quote":"Supplies the reduction of periodic-orbit continuation to a scalar bifurcation function, which the paper extends to non-integrable multi-degree-of-freedom systems."},{"cited_title":"Sepulchre and R.S","cited_arxiv_id":null,"evidence_quote":"Defines $m$-normal periodic orbits and guarantees a smooth one-parameter family of such orbits."},{"cited_title":"Muñoz-Almaraz, E","cited_arxiv_id":null,"evidence_quote":"Establishes continuation of periodic orbits in conservative systems, backing the normality and family-persistence assumptions."},{"cited_title":"Peeters, G","cited_arxiv_id":null,"evidence_quote":"Introduces the phase-lag quadrature criterion that Proposition 4.3 justifies under weaker assumptions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Formulates the energy principle for locating forced responses that the Melnikov function makes rigorous."},{"cited_title":"Golubitsky and S","cited_arxiv_id":null,"evidence_quote":"Classifies saddle-node, isola-birth, and simple-bifurcation singularities used in Theorem 3.2."},{"cited_title":"Dankowicz and F","cited_arxiv_id":null,"evidence_quote":"Supplies the numerical continuation toolbox used to compute the example backbone curves and frequency responses."},{"cited_title":"Melnikov","cited_arxiv_id":null,"evidence_quote":"Introduces the original scalar Melnikov integral for planar oscillators that this work generalizes."}],"review_version":1}