{"id":"5a074481-84f8-42e8-9b1a-82b641198b0b","arxiv_id":"2508.16704","paper_version":1,"verdict":"REJECT","confidence":"LOW","novelty_score":3.0,"correctness_risk":"high","formal_verification":"none","parameter_count":1,"one_line_summary":"The paper tabulates harmonic balance stability ranges for n-fold symmetric elastic ring buckling modes.","lead":"This short preprint applies harmonic balance to the classic problem of a circular elastic ring under uniform pressure and tabulates stability ranges for symmetric buckling modes n=2 through 12. The paper is a lightweight extension of the author's prior third-order harmonic balance work and does not derive its key stability results in the text.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Table 1's stability ranges are not derived in the paper: the third-order HBM equations and the stability criterion (second variation / eigenvalue) are absent, so the central claim is unverifiable from the text.","rationale":"The paper's title and abstract promise a stability analysis, but the quantitative content is Table 1: a list of stable p-ranges. For the central claim to hold, two things are required: (i) accurate equilibrium solutions from the harmonic balance method, and (ii) a correct stability test applied to those solutions. The manuscript supplies neither in reproducible form. Section 5 explicitly states that the stability diagrams are computed with the third-order HBM formulation from [15], but that formulation is not reproduced; the displayed Sec. 3 equations are only zeroth- and first-order and are garbled (undefined nn and M in the nu1 expression). More importantly, no stability criterion is defined. The only candidate in the text is self-intersection, which is a geometric property, not a Lyapunov or energetic stability criterion. Without specifying the second variation or an eigenvalue problem, the endpoints of Table 1 have no derivable meaning. The claimed agreement with 'bifurcation theory of linearization' is likewise asserted without equations or numerical comparison. This is not merely an exposition issue: different stability criteria yield different ranges, and the range is the entire result. The reader's weakest assumption identified the reliance on the prior work [15]; my concern sharpens this to the unstated stability criterion plus the missing third-order equations. The rejection is therefore justified, but it rests on unverifiability and incompleteness rather than on a demonstrated counterexample. No change to the reader's verdict is needed.","tokens_in":4706,"tokens_out":5494,"duration_ms":65850,"concrete_test":"Reproduce Table 1 from the exact ring BVP: numerically continue equilibrium branches of Eq. (1) with (2)-(3) for each n; at each p evaluate the second variation of the elastica energy on the symmetric subspace (equivalently, solve the linearized BVP for perturbations dnu satisfying the same boundary conditions and identify the smallest eigenvalue). Mark the branch unstable when that eigenvalue changes sign. Additionally, obtain from [15] the third-order HBM residual equations and re-derive Table 1 from them. If the sign-change p does not match the Table endpoints, or if the HBM reproduction differs, the claimed stability ranges are unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the mapping from the HBM equilibrium family to the 'stable range of p' in Table 1. Section 5 does not present that mapping. It cites [15] for a third-order HBM solution and then jumps to stability diagrams; the only criterion named is geometric: 'if we exceed the external loading, we would get intersecting shape... break-even point, where the stable shapes becomes unstable.' Self-intersection is a configuration-space event, not a stability boundary: an equilibrium with positive second variation remains Lyapunov stable even if the curve self-intersects, and contact changes the mechanics entirely. No second variation of the elastica energy, no linearized perturbation eigenvalue problem, and no Floquet/Jacobi analysis is given. Consequently the endpoints of Table 1 (e.g., 10.5 for n=2, 1191 for n=12) are asserted rather than computed. The zeroth- and first-order equations shown in Sec. 3 contain undefined symbols (nn, M in the nu1 formula) and do not determine the third-order solution used; the claimed agreement with 'bifurcation theory of linearization' (Sec. 5/6) has no equations, so it cannot be checked. Thus the central claim of accurate stability ranges rests entirely on unreproduced equations and an unspecified stability criterion.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the planar Euler elastica ring under uniform, n-fold symmetric pressure p. The author solves the nonlinear curvature-difference equation (Eq. 1) by harmonic balance: a zeroth-order approximation ν0(s) = A cos(ns) is derived in Section 3, first-order corrections are sketched (with poorly formatted formulas), and the third-order solution is taken from the author's prior work [15]. Section 4 presents equilibrium shapes and discusses contact and self-intersection transitions. Section 5 presents stability diagrams (p vs ν(0)) and Table 1, which lists claimed stable p-intervals for loading modes n = 2 through 12, asserted to begin at p = n² − 1 and to end at a 'break-even point' identified with the onset of self-intersecting shapes. The abstract claims that the harmonic balance method gives accurate stability calculations and agrees with the bifurcation theory of linearization.","tokens_in":5062,"tokens_out":7457,"duration_ms":78732,"significance":"If the stability ranges were rigorously derived, the paper would provide a useful systematic map of stable p-intervals for all symmetric buckling modes of the pressurized elastica ring, complementing the classical p > n² − 1 condition of Tadjbakhsh–Odeh and the elliptic-function descriptions of Djondjorov et al. The zeroth-order harmonic balance balance (Section 3, Eq. (4)) is algebraically consistent and worth acknowledging. However, the manuscript ships no code, no machine-checked proof, and reproduces none of the equations that determine Table 1; the third-order formulation of [15] is invoked but not given, and no independent numerical or analytical verification of the table is presented. The claimed agreement with the bifurcation theory of linearization is stated only in words. Consequently, the paper's central quantitative claim is not checkable from the text, and the significance is contingent on material the manuscript does not provide.","major_comments":[{"comment":"The central result—the stable p-intervals in Table 1—is not derived in the manuscript. The only stability criterion stated is geometric: 'if we exceed the external loading, we would get intersecting shape… break-even point, where the stable shapes becomes unstable' (end of Sec. 4; Sec. 5). Self-intersection of the centerline is not a stability boundary: an equilibrium with positive second variation of the elastica energy remains stable even when the curve self-intersects, and contact changes the mechanical system entirely. The paper gives no second-variation calculation, no linearized perturbation eigenvalue problem, and no Jacobi/Floquet condition. The endpoints p = 10.5 (n = 2), 70 (n = 3), 1191 (n = 12) are therefore asserted, not computed, and the abstract's claim of 'accurate stability calculations' is unsupported by the text.","section":"§5, Table 1; end of §4"},{"comment":"The displayed expressions for ν1, μ1, β1 are not well-formed: the symbols 'nn', 'M', and 'An' are never defined, and factors such as '8Ann − 32An2' and 'nn − 5n2 − 3' are illegible as mathematics. The relation of ν1 to ν0 is unexplained, and the iterative scheme (5) is not carried beyond the first order, although Section 5 states that the stability results come from the third-order formulation of [15], which is not reproduced. The only harmonic balance derivation actually present in this paper, the zeroth-order balance of Eq. (4), does not determine the quantities reported in Table 1.","section":"§3, first-order expressions after Eq. (5)"},{"comment":"The stability analysis explicitly delegates the solution to the author's own prior paper: 'In our previous work [15], we used the third-order harmonic balance formulation… Therefore, we employ the same formulation to compute stability diagrams.' Since neither that formulation nor the stability boundary is reproduced, and the agreement with 'bifurcation theory of linearization' is claimed without equations, Table 1 rests on an unreproduced, self-referential basis. The statement that 'NSolve gives the same results' is not accompanied by any numerical values or comparison plots, and the offered Mathematica notebooks (Sec. 5) are not shipped. A reader cannot reproduce or verify the paper's predictions from the text.","section":"§5, opening sentence"}],"minor_comments":[{"comment":"The sentence 'Troger and Steindl ref and Chaskalovic and Naili ref' contains unfilled placeholder 'ref' markers and the corresponding references are missing from the bibliography.","section":"§1, Introduction"},{"comment":"This paragraph is garbled: pcp and psi are undefined, and the condition 'x−2 = x+1' is illegible as typeset. The contact/self-intersection condition should be formulated explicitly with clear notation.","section":"§4, paragraph 'As external loading move to thepcp < psi'"},{"comment":"The paper says 'Section 5 concludes the paper', but Section 5 is 'Stability Analysis' and the conclusion is actually Section 6.","section":"§1, last paragraph"},{"comment":"The row for n = 11 is omitted without explanation, and the notation '(3, 10.5)' does not state whether endpoints are included or how the endpoints were extracted from the diagrams.","section":"Table 1"},{"comment":"The caption lists p-values (5.247, 21.65, 51.844, 97.834) that are not connected to Table 1 or to the known contact pressures from the elliptic-function or numerical literature [6, 8]; no comparison with those values is given.","section":"§4, Figure 2 caption"},{"comment":"The claim that 'Previous researcher have not computed the stability range and bifurcation diagrams … using elliptic functions' is an unsupported literature assertion; a comparison with at least one existing stability or eigenvalue analysis of elastic rings would be needed to substantiate it.","section":"§5, penultimate paragraph"},{"comment":"The inline fractions are ambiguous: '3/2ν(s)2 + 1/2ν(s)3' should be typeset as (3/2)ν(s)² + (1/2)ν(s)³ to avoid misinterpretation.","section":"Eqs. (1) and (4)"},{"comment":"The stability diagrams are never compared quantitatively with the numerical bifurcation diagram (Fig. 6), which is shown only for n = 5; the claimed agreement between HBM and lin- earization bifurcation analysis cannot be assessed from the figures.","section":"Figures 4–6"}],"recommendation":"reject","confidential_remarks":"For the editor: in its current form the manuscript does not meet the standard for publication. The paper contains unfilled reference placeholders, equations with undefined symbols, a misstated section structure, and—decisively—a central result (Table 1) whose derivation and stability criterion are absent. The stability boundary is identified with self-intersection, which is a configuration-space event rather than a condition on the second variation, so the criterion itself is questionable. The substantial reliance on the author's own prior work [15] without reproduction compounds the verifiability problem. These are load-bearing deficiencies in the core contribution, not local presentation issues; if the authors can supply the complete third-order HBM equations, a genuine second-variation or eigenvalue stability analysis, and a verified table, the work could be reconsidered as a substantially revised submission."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The headline is that this paper promises stability ranges for Euler elastica rings at modes n=2..12, but the table at the center of the paper is not derived anywhere in the manuscript. It is imported from the author's own prior paper [15] via a single sentence, and the only criterion named for 'stability' is the point of self-intersection. That is a geometric event, not a stability boundary, so the central claim is unverifiable as written.\n\nWhat is actually new: tabulating ranges for n up to 12 is a modest extension of the known harmonic-balance program (Wu et al. did n=2; the author's own paper did general n). The zeroth-order equations in Sec. 3 are algebraically consistent, and the ring model with closure conditions is standard and correctly stated. The mention of comparing against shooting/NSolve is fine, but no comparison data are shown.\n\nThe soft spots are several. First, the first-order expressions in Sec. 3 are garbled: undefined symbols (nn, M) and apparent typos ('8Ann','32An2') make them uncheckable. The third-order equations used to produce Table 1 are not given at all. Section 5 explicitly says 'we employ the same formulation' from [15], so Table 1 is an assertion, not a computation in this paper. Second, the stability criterion is missing. The text equates the 'break-even point' with where 'the stable shapes becomes unstable' due to self-intersection. Self-intersection in a 2D model is not the same as loss of Lyapunov stability; it is a configuration-space phenomenon, and contact would change the mechanics. No second variation of the energy, no linearized eigenvalue problem, no Jacobi/Floquet analysis appears. Third, the claimed agreement with 'bifurcation theory of linearization' has no equations, so it cannot be checked. Finally, the remark that 'previous researcher have not computed the stability range... using elliptic functions' is irrelevant and unsupported.\n\nOn the positive side, the author is clearly aware of the relevant literature and the zeroth-order derivation is correct. If the stability ranges in Table 1 are intended to be taken from [15], then the present paper adds little beyond a re-plot; if they are new, they need derivation. The phrase 'the author can provide the Mathematica notebooks upon request' does not replace that.\n\nWho is this for? A reader who already knows [15] and wants a quick lookup table might bookmark it, but you cannot cite the table independently. This is not a paper I would bring to a reading group except to illustrate how a stability claim should not be presented. It deserves a desk reject rather than a referee round, unless the author incorporates the missing formulation and a genuine stability analysis.\n\nRecommendation: reject/desk-reject; do not send to peer review as is.","headline":"Table 1's stability ranges are not derived in the paper; the 'stability' criterion is geometric self-intersection, so the central claim is unverifiable as written—though the zeroth-order algebra is fine.","tokens_in":5486,"tokens_out":4069,"would_cite":false,"duration_ms":43469,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that harmonic balance gives accurate stability ranges for the pressurized Euler elastica ring, yielding smooth stable-pressure intervals for symmetric modes up to n = 12.","keywords":["Euler's elastica ring","harmonic balance method","stability analysis","post-buckling","bifurcation","self-intersection","uniform pressure","symmetric loading"],"falsifier":"Recompute the equilibrium branch for, say, n = 2 and n = 5 by shooting or finite differences with fine arc-length continuation, and check whether a closed, non-self-intersecting ring still exists just above the upper bounds in Table 1 (p = 10.5 and p = 170). If a stable non-self-intersecting solution exists beyond the listed endpoint, the claimed stability boundary is not the true one; if continuation also terminates at exactly those pressures, the table is confirmed.","tokens_in":4631,"feed_emoji":"⭕","tokens_out":10496,"duration_ms":105610,"temperature":0.7,"pith_summary":"This paper claims that the harmonic balance method gives reliable stability ranges for a thin elastic ring compressed by uniform pressure, for symmetric buckling modes up to n = 12. The central result is Table 1, which lists the interval of the loading parameter p in which stable equilibrium shapes exist for each mode: for example, (3, 10.5) for n = 2 and (143, 1191) for n = 12. The lower end is the classical threshold n^2 − 1; the upper end is where the computed shape self-intersects and the branch becomes unstable. The paper checks this against a linearization/bifurcation analysis and against numerical shooting, finding agreement, and presents stability diagrams showing the closed-to-contact-to-over-contact-to-irregular progression. These ranges matter because they predict when a pressurized ring will collapse or self-intersect, which is relevant for biological membranes, soft actuators, and structural elements.","feed_headline":"Harmonic balance maps stable pressure windows for buckled rings","feed_subtitle":"New tables pinpoint, for modes 2 through 12, where pressurized ring shapes become self-intersecting.","key_machinery":"The load-bearing mechanism is the harmonic balance expansion of the curvature difference ν(s) as a cosine series in cos(nks), solved order by order; the paper uses the third-order version from its earlier work. This turns the nonlinear differential equation into algebraic equations for the Fourier coefficients, giving the loading parameters µ, β, and hence p = µ + β − 1, as functions of the amplitude A. Varying A traces the equilibrium branch p versus ν(0); the upper endpoint of the branch, where the computed ring self-intersects, defines the claimed stability boundary. A second, independent check is the linearization/bifurcation analysis, which marks the same unstable parameter region.","core_discovery":"The paper expands the curvature difference ν(s) as an n-fold cosine series and finds its coefficients by harmonic balance. Using the third-order approximation from the author's earlier work, it traces the equilibrium branch for each mode n as the pressure parameter p grows. The branch begins at p = n^2 − 1, passes through closed, contact, and over-contact shapes, and ends where the computed curve self-intersects; that endpoint is taken as the loss of stability. The endpoints form Table 1's stable-pressure intervals. Linearization about the branch and numerical shooting recover the same intervals. The claimed result: harmonic balance gives smooth, accurate stability ranges for all symmetric m","pith_inferences":["A pattern in Table 1's upper endpoints (close to 8n^2 for n = 8, 9, 10, 12) suggests p_max may follow a simple quadratic law; fitting this relation and testing it for more modes would give a closed-form design rule the paper does not state.","The identification of self-intersection with loss of stability is made numerically; proving that the self-intersection value of p coincides with the first unstable eigenmode of the linearized problem would turn the table into a theorem.","Since the harmonic balance ansatz keeps only symmetric cosine modes, the stability statement concerns symmetric perturbations; checking non-symmetric or dynamic perturbations would test whether the branch is stable in full function space."],"forward_implications":["For each symmetric mode n, the stable window (n^2 − 1, p_max) gives the full pressure range over which a ring keeps a non-self-intersecting equilibrium shape.","The two independent routes — harmonic balance and linearization — agree on the unstable regions, so either can be used to cross-check buckling predictions in related problems.","The harmonic balance stability data extend to higher modes such as n = 10 and n = 12, a range not tabulated by earlier elliptic-function solutions.","The closed-to-contact-to-over-contact-to-irregular sequence gives a geometric reading of how a ring approaches instability, visible directly in the computed shapes."],"supporting_citations":[{"why":"Supplies the third-order harmonic balance formulation on which all stability diagrams and Table 1 rest.","marker":"[15]"},{"why":"Derives the governing equation and the lower bound p > n^2 − 1 for the symmetric modes.","marker":"[21]"},{"why":"Provides the numerical post-buckling contact shapes that define the closed–contact–over-contact sequence used in the stability diagrams.","marker":"[8]"},{"why":"Introduces the harmonic balance approach for the n = 2 elastica ring, which this paper extends to all modes.","marker":"[24]"},{"why":"Provides the iterative harmonic balance procedure used to generate the higher-order approximations.","marker":"[12]"},{"why":"Supplies the equilibrium shapes for n = 2 through 5 reproduced as starting configurations in the stability plots.","marker":"[14]"}],"fun_headline_variants":["Harmonic balance pinpoints ring buckling limits for modes 2-12","Buckled ring stability windows revealed by harmonic balance","Self-intersection marks instability boundary in ring buckling","Two methods converge on ring stability pressure windows","Harmonic balance vs linearization: same ring buckling ranges"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The stability intervals in Table 1 are correct only if the third-order harmonic balance approximation from the author's earlier paper reproduces the true equilibrium shapes and if the pressure at which the computed shape self-intersects really is the point of instability.","fun_headline_variants_meta":{"raw":{"variants":["Harmonic balance pinpoints ring buckling limits for modes 2-12","Buckled ring stability windows revealed by harmonic balance","Self-intersection marks instability boundary in ring buckling","Two methods converge on ring stability pressure windows","Harmonic balance vs linearization: same ring buckling ranges"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000746,"raw_usage":{"total_tokens":3083,"prompt_tokens":583,"completion_tokens":2500,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":327,"completion_tokens_details":{"reasoning_tokens":2422}},"tokens_in":327,"tokens_out":2500,"duration_ms":20483,"temperature":1.0,"reasoning_tokens":2422,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T17:22:42.041518+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the equilibrium branch for, say, n = 2 and n = 5 by shooting or finite differences with fine arc-length continuation, and check whether a closed, non-self-intersecting ring still exists just above the upper bounds in Table 1 (p = 10.5 and p = 170). If a stable non-self-intersecting solution exists beyond the listed endpoint, the claimed stability boundary is not the true one; if continuation also terminates at exactly those pressures, the table is confirmed.","supporting_citations":[{"cited_title":"International Journal of Mechanical Sciences 2007;49:661-668","cited_arxiv_id":null,"evidence_quote":"Introduces the harmonic balance approach for the n = 2 elastica ring, which this paper extends to all modes."},{"cited_title":"and He, X., Newton-harmonic balancing approach for accurate solutions to nonlinear cubic-quintic Duffing oscillators","cited_arxiv_id":null,"evidence_quote":"Provides the iterative harmonic balance procedure used to generate the higher-order approximations."},{"cited_title":"Self-contact of a flexible loop under uniform hydrostatic pressure","cited_arxiv_id":null,"evidence_quote":"Supplies the third-order harmonic balance formulation on which all stability diagrams and Table 1 rest."},{"cited_title":"1997;18:59-74","cited_arxiv_id":null,"evidence_quote":"Derives the governing equation and the lower bound p > n^2 − 1 for the symmetric modes."},{"cited_title":"and Rubinow, S.I., Post buckling behavior of elastic tubes and rings with opposite sides in contact.SIAM J","cited_arxiv_id":null,"evidence_quote":"Provides the numerical post-buckling contact shapes that define the closed–contact–over-contact sequence used in the stability diagrams."},{"cited_title":"Simulation of a soap film spanning a flexible loop.International Journal of Non-Linear Mechanics","cited_arxiv_id":null,"evidence_quote":"Supplies the equilibrium shapes for n = 2 through 5 reproduced as starting configurations in the stability plots."}],"review_version":1}