{"id":"acadd6f1-a5a3-42d0-98c0-9bf4181c7d8c","arxiv_id":"2601.06598","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":1,"one_line_summary":"A shallow arch's response to patterned preferred curvature is governed by the overlap of the pattern with Euler buckling modes, yielding design rules for snapping, energy release, and switching paths.","lead":"This paper shows that the shape of the bending pattern stamped onto an elastic arch decides how and when it snaps, and that the pattern's overlap with the arch's natural buckling modes predicts the snapping threshold. It gives simple design rules for maximizing snap energy or minimizing stimulation, verified in liquid-crystal elastomer arches.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The load-bearing issue is scope: Eq. 2's shallow linearization and the 1D strip energy are only validated at ε≈1.2% for one width, yet the abstract claims the mode-decomposition spectrum 'completely describes any such system' — Fig. S7 already shows up to 8% threshold error at ε=0.05.","rationale":"The reader's weakest-assumption exactly identifies the same load-bearing concern: the shallow-angle, 1D elastica model is only quantitatively validated near ε≈1.2% and for the chosen narrow strip, while the abstract and several forward-looking statements claim complete predictive generality. The paper itself supplies the supporting evidence for the limitation: Fig. S7 reports up to 8% error at ε=0.05, and S2.3 concedes that intermediate widths require a 2D plate description. Because the central claim of the paper is built on this mode-decomposition spectrum, the domain restriction is not cosmetic: it controls whether the theory 'completely describes' or only describes a limited shallow-sender regime. I found no stronger internal inconsistency: the mode orthogonality proof in S3.3, the resonance form of Eq. 3, and the stability argument in S3.4 are coherent, and the predictions are corroborated by fully nonlinear numerics and by experiments at low ε with no target fitting. The appropriate disposition remains CONDITIONAL: accept the mechanics but require the authors to explicitly scope the 'complete' claim to shallow arches in the flat-clamped, narrow/very-wide strip regime, and to provide error statistics or code/data release. Since the reader already reached this conditional verdict, no verdict adjustment is needed.","tokens_in":20860,"tokens_out":10218,"duration_ms":110702,"concrete_test":"Take a set of representative κ̄(s) profiles (κ0, κ1, central/offset rectangles, κmaxE, κminThr, and random mode mixtures) and compare the Eq. 3 compression-spectrum predictions (threshold, snapping mode, center-height tracks) against the full nonlinear solve_bvp and discrete-rod solvers at ε=0.01, 0.05, 0.10, 0.15, and against 2D MorphoShell at w = 3.5, 10, 20 mm. If threshold or mode predictions deviate by more than ~5% anywhere inside the claimed 'any such system' domain, the abstract must be qualified to the shallow, slender/very-wide, low-ε regime.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Central claim: Eq. 3's compression spectrum built from κ̄(s) mode decomposition completely predicts morphology, snapping threshold, and snapping mode for flat-clamped patterned arches. This requires the linearized ODE (Eq. 2: θ''+fθ+v=κ̄', small θ) and the 1D bending energy to remain quantitatively valid along the entire loading path up to the snap. The paper's own SI gives the two limiting facts: Fig. S7 shows shallow-vs-nonlinear snapping-threshold deviations up to 8% already at ε=0.05 (and growing for larger ε), and S2.3/MorphoShell shows the 1D energy is reliable only for narrow strips (or very wide sheets), with intermediate widths requiring a full 2D plate treatment. All experiments quoted are at ε=1.2% on a single 3.5 mm-wide sample, i.e., inside the good-agreement pocket. The abstract's 'completely describes any such system', and the implied predictive generality for MEMS/robotics/meta-materials, is therefore stronger than the demonstrated envelope. This is a scope overclaim rather than an internal inconsistency; the theoretical construction itself is consistent and is independently supported by fully nonlinear numerics and by experiments at low compression.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a theoretical framework for flat-clamped elastic arches with patterned preferred curvature, backed by liquid-crystal-elastomer experiments and nonlinear numerics. The authors linearize the elastica equation (Eq. 2), decompose the preferred curvature into Euler-buckling modes, and derive a per-mode compression relation (Eq. 3). The resulting 'compression spectrum' predicts arch morphology, snapping threshold, and the mode through which snapping occurs. They use it to design binary stimulation profiles that maximize released energy or minimize the snapping threshold, to create symmetric or asymmetric switching pathways, and to revisit snapping of passive arches with clamped end angles. The theory is internally consistent and is checked against nonlinear solve_bvp, discrete-rod simulations, and experiments at ε ≈ 1.2%.","tokens_in":21223,"tokens_out":4596,"duration_ms":51369,"significance":"If the scope is stated accurately, this is a valuable contribution. The modal decomposition is elegant: it reduces a weakly nonlinear arch problem to a one-dimensional spectrum that yields simple design rules, and the paper demonstrates surprising and useful results (e.g., non-overlapping stimulation drives continuous higher-order transitions, and offset binary patterns optimize energy release and threshold). The paper is strengthened by complementary numerical methods, reproducible solve_bvp and discrete-rod implementations, and experimental validation of hysteresis, continuous transitions, and optimized profiles. The main weakness is that the abstract and conclusions claim more generality than the validated envelope: the shallow linearization (Eq. 2) is only quantified for small compression, and the 1D bending energy only applies to narrow or very wide strips. The paper itself documents these restrictions, so the fix is a matter of qualification rather than a fundamental flaw.","major_comments":[{"comment":"The abstract's claim that the shallow-arch theory 'completely describes any such system' is broader than the validated scope. Fig. S7 shows that shallow-vs-nonlinear snapping-threshold deviations reach up to 8% already at ε = 0.05 and grow for larger ε, while §S2.3 states that the 1D energy is reliable only for narrow or very wide strips. The experiments shown are at ε ≈ 1.2% for a single 3.5 mm wide strip. Since Eq. 3 is built on the linearization, the claims of 'complete' description and of direct predictive value for MEMS/robotics/metamaterials should be qualified by stating the quantitative validity envelope (small ε, 1D strip-width regime). This is a scope overclaim, not an internal inconsistency, but it is load-bearing for the central claim and should be corrected.","section":"Abstract, §S3, Fig. S7"},{"comment":"The theory treats the compression ε as a fixed constraint, but the experiments measure a variable compression along the loading path (Δε ≈ 0.075% at f0 and ≈ 0.15% at f1, §S1.5), caused by stretching of the strip under compressive force. The theoretical predictions shown in Figs. 1–3 are therefore not entirely closed-loop predictions from a prescribed clamp displacement and stimulus; they use the experimentally measured ε. The manuscript should state this explicitly and quantify when the fixed-ε approximation is predictive. This does not undermine the low-ε results reported, but it is important for the claimed predictive power of Eq. 3.","section":"Eq. (3), §S1.5"},{"comment":"The linear-stability analysis at the spectral minimum fm establishes the onset of a growing perturbation of the form ∂f w0. It is implicitly assumed that this loss of linear stability coincides with a snap to a different equilibrium, which is true when no other stable branch is available in the nonlinear problem. The authors do verify this with nonlinear solve_bvp for the specific cases in Fig. S7, but the paper should state that the identification of fm with the snapping threshold relies on this nonlinear branch structure, rather than on the linear analysis alone. This is a clarity issue, not a correctness error.","section":"§3.4, Fig. 2(b)"}],"minor_comments":[{"comment":"The title contains a typo: 'Pref erred' should be 'Preferred'. Also the abstract and main text use slightly different spacing/formatting for κ̄; please harmonize.","section":"Title / Abstract"},{"comment":"The axes and curves in the 'compression spectrum' plot are not defined in the caption. A reader needs to know what g(f) is, what the horizontal line represents, and what the dashed/solid curves mean; please expand the caption.","section":"Fig. 1(c)"},{"comment":"In the piecewise expression for θ(s), the third branch uses 'x' instead of 's' in 'γ/2 ≤ x'. Please correct.","section":"§S3.8, Eq. (S51)"},{"comment":"The Fourier decomposition uses k = 701 modes, but no convergence criterion is stated. Since the compression spectrum is central, please include a short statement on convergence with respect to k and the number of quadrature points.","section":"§S3.7"},{"comment":"The observation of a symmetric snap-back during edge stimulation is attributed to thermal diffusion (reduction in α0 and increase in α2). This is plausible but speculative; label it as a hypothesis or support it with a temperature-profile measurement.","section":"§4 (symmetry of snap-back)"},{"comment":"Please add a data/code availability statement, given the custom numerical routines (solve_bvp, discrete rod, optimization) that are essential to reproduce the figures.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"This is a strong paper with an elegant and useful theory, and the experimental work is impressive. The main reason for major revision is the scope overclaim in the abstract; the authors themselves provide the evidence for the restriction (Fig. S7, §S2.3), so I expect a straightforward revision. I would be comfortable with acceptance once the claims are qualified and the fixed-ε/variable-ε issue is clarified."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a genuinely useful paper and I'd send it to review without hesitation. The modal-decomposition picture — expanding arbitrary preferred curvature into Euler buckling modes and reading snapping off the resulting compression spectrum (Eq. 3) — is new and organises a lot of scattered empirical results on central-segment stimulation. The design rules (binary profiles that maximise energy release or minimise threshold, symmetric vs antisymmetric pathways, continuous switching around the critical point) are practical output that makes the theory worth having. The experiments at ε=1.2% are solid, and the SI is thorough: full nonlinear solve_bvp, discrete-rod dynamics, 2D MorphoShell checks, and mode-orthogonality proofs.\n\nNow the soft spots, in proportion. The main one is scope, and it's real. The abstract claims the shallow-arch theory 'completely describes any such system', but the theory is built on the linearised Eq. 2 and a 1D bending energy that, by the authors' own SI (S2.3), holds only for narrow strips or very wide sheets. Fig. S7 shows the shallow theory misses the nonlinear snap threshold by up to 8% already at ε=0.05, growing at larger compression. All headline experiments sit at ε≈1.2% on a single 3.5 mm-wide sample. The demonstrated envelope is much narrower than the claim. That is an overclaim, not an internal inconsistency — the construction is consistent and independently backed by nonlinear numerics — but the abstract needs qualifying to the flat-clamped, shallow, 1D regime.\n\nTwo smaller things. The experimental evidence rests on one sample, no repeated-sample error bars; convincing hysteresis loops, but I'd want statistics before the numbers harden into design tables. No code or raw data deposited — the Fourier-decomposition and solve_bvp pipeline is described carefully enough to reproduce, but shipping it would make the design rules properly usable.\n\nNone of this touches the central mechanism. The compression-spectrum idea is correct, the match to nonlinear numerics is good, and the circularity burden is low: measured κ(s) and ε go in, thresholds and energies come out. The self-cited [22] is used exactly where it belongs, to set the 1D/2D boundary.\n\nWho's it for: bistable-beam people, MEMS switch designers, LCE actuator folks; the pathway-symmetry results will interest the metamaterials crowd. Send it to a serious referee with a request to check the claim-envelope tightening lands. Accept after revision.","headline":"Solid, useful theory with good experiments — but the abstract's 'completely describes any such system' outruns the validated shallow, 1D envelope.","tokens_in":21671,"tokens_out":2824,"would_cite":true,"duration_ms":26281,"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":"When an elastic arch carries a patterned preferred curvature, its entire snapping behavior — morphology, threshold, and snapping mode — is predicted by how that pattern decomposes into Euler-buckling modes.","keywords":["elastic arch","bistable switching","snap-through","preferred curvature","Euler buckling modes","compression spectrum","liquid crystal elastomer","photo-thermal actuation"],"falsifier":"Measure the snapping threshold of a flat-clamped arch stimulated by a pure κ1 pattern (orthogonal to the fundamental mode) at compression ε = 0.1: the shallow spectrum predicts a continuous transition with no snap, whereas the paper's own nonlinear solver shows the threshold moves with compression; a measured discontinuous snap there would contradict the modal-decomposition claim.","tokens_in":20769,"feed_emoji":"🌀","tokens_out":6788,"duration_ms":70386,"temperature":0.7,"pith_summary":"This paper addresses a classic bistable structure — the compressed elastic arch — and asks what happens when the arch's preferred curvature is not uniform but patterned in space and time, as light-driven liquid-crystal elastomers allow. The authors claim that a shallow-arch elastica theory gives a complete description: decompose the preferred-curvature pattern into Euler-buckling modes, and the compression spectrum built from Equation 3 predicts the arch's shape, its snapping threshold, and which mode it snaps through. If the pattern overlaps the fundamental mode, the arch snaps up/down; if it is orthogonal, stimulation drives a smooth second-order transition into a higher-order shape. The same spectrum lets a designer pick binary patterns that maximize released energy or minimize the stimulation threshold, and design symmetric or asymmetric switching pathways, with the same machinery extending to passive arches snapped by boundary-angle control. A single experimental arch demonstrates these predictions under light, giving a direct design rule for switchable mechanical elements.","feed_headline":"Arch snapping reduces to a sum over buckling modes","feed_subtitle":"A pattern-based spectrum predicts when an arch snaps, how it snaps, and how to make it switch without snapping.","key_machinery":"The central object is the compression spectrum: the relation ε = Σ α_i² f_i²/(f−f_i)² obtained from the shallow-arch elastica (θ'' + fθ + v = κ') by decomposing preferred curvature into compression-normalized Euler-buckling modes with amplitudes α_i. Each mode contributes a divergence at its buckling force f_i; unsplit modes appear as degenerate spikes. The spectrum maps all possible arch states, and instability is identified where the horizontal line of fixed compression/stimulation can no longer intersect a solution branch — an endpoint or a minimum between divergences — reducing snap prediction to reading a spectral plot.","core_discovery":"Expressing preferred curvature as a sum of Euler-buckling modes, the paper derives an algebraic compression spectrum, ε = Σ α_i² f_i²/(f−f_i)², and reads arch states directly from it: solutions slide along branches as stimulation grows, and snapping happens when the branch ends — at the degenerate spike of an unsplit odd mode or at a minimum between two split divergences. Stimulation orthogonal to the fundamental mode drives a continuous second-order transition to a higher-order shape with vanishing central displacement. The same spectrum yields optimized binary patterns for maximum energy release and minimum threshold, and a phase diagram whose critical point can be circled to switch up/dow","pith_inferences":["Because the compression-spectrum construction depends only on the linear buckling modes of the base system, the same design logic should transfer to other slender bistable structures — pinned or mixed-clamped beams, rings, or shells — wherever a mode basis exists.","A testable extension would be temporal patterning: continuously morphing the pattern (rather than the amplitude) could navigate the phase diagram to produce non-reciprocal or sequential snapping sequences in a single arch, enabling multi-step mechanical logic.","The paper's robustness observation (a stimulation width of γ ≈ 2π/k1 ≈ 0.7 decouples from mode 1) suggests fabricating patterns with that width will make arch switches far more tolerant of misalignment in real devices.","The experiments sit at ε = 1.2% where the shallow theory is accurate; extending the design rules to strongly compressed arches will need the nonlinear solver, but the qualitative mode-selection picture may survive beyond the validated range."],"forward_implications":["A designer can prescribe preferred-curvature patterns to make an arch snap or switch smoothly, and to control the symmetry of the pathway: including κ1 drives antisymmetric snaps, while symmetric multi-mode patterns (κ0 + κ2 + ...) enable symmetric through-snaps.","Optimized binary patterns exist: a single offset segment releases up to 0.661 α² at threshold α = 10.33√ε, more than twice the 0.239 α² of pure fundamental-mode stimulation, while the low-threshold pattern snaps at α = 9.11√ε versus 26.28√ε.","Boundary-driven passive arches fit the same formalism: clamping-angle patterns act like preferred-curvature mode amplitudes, explaining why asymmetric clamping gives continuous transitions and symmetric clamping gives antisymmetric snap-through.","Stimulating orthogonal to the fundamental mode yields a reversible second-order morphological transition rather than a snap — a way to morph an arch into normally unstable higher-order shapes.","Switching without snapping is possible by looping around the critical point in the (α0, α1) phase diagram, returning to the same stimulation level in the opposite bistable well."],"fun_headline_variants":["Arch snapping predicted by buckling-mode spectrum","Patterned curvature tunes arch snap and switch","Euler modes decode arch bistability","Preferred curvature controls arch snapping","Buckling-mode sum dictates arch behavior"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the small-angle, one-dimensional elastica model remains quantitatively valid all the way up to the snap threshold; the paper's own comparison (Fig. S7) shows the shallow theory deviates by up to 8% in snapping thresholds once compression reaches ε ≈ 0.05, and the 1D bending energy holds only for narrow or very wide strips.","fun_headline_variants_meta":{"raw":{"variants":["Arch snapping predicted by buckling-mode spectrum","Patterned curvature tunes arch snap and switch","Euler modes decode arch bistability","Preferred curvature controls arch snapping","Buckling-mode sum dictates arch behavior"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000142,"raw_usage":{"total_tokens":991,"prompt_tokens":718,"completion_tokens":273,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":462,"completion_tokens_details":{"reasoning_tokens":227}},"tokens_in":462,"tokens_out":273,"duration_ms":3278,"temperature":1.0,"reasoning_tokens":227,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T11:21:53.436127+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the snapping threshold of a flat-clamped arch stimulated by a pure κ1 pattern (orthogonal to the fundamental mode) at compression ε = 0.1: the shallow spectrum predicts a continuous transition with no snap, whereas the paper's own nonlinear solver shows the threshold moves with compression; a measured discontinuous snap there would contradict the modal-decomposition claim.","supporting_citations":[],"review_version":1}