{"id":"a7312a90-18ec-4412-bd65-48e66dc01e65","arxiv_id":"2412.03725","paper_version":1,"verdict":"REJECT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"The paper claims sawtooth body shapes are analytically optimal for snake locomotion on flat surfaces, but the variational proof is flawed.","lead":"This paper presents an analytical framework for snake locomotion, claiming that sawtooth waveforms optimize efficiency and speed on flat surfaces under anisotropic friction. The central variational derivation has a critical gap, so the main claim is not established.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The variational core of §3 is unsound: Eqs. 13 and 17 differentiate integrands with respect to θ while treating U, W, ϕ, ψ, and the closure constraint as independent of θ, although the paper itself states that all kinematic variables are functions of θ.","rationale":"The reader’s weakest assumption identifies precisely the load-bearing defect: the variational derivations in §3.1 and §3.2 treat U and other global kinematic variables as constants under shape variation, even though Eq. 9 defines U as an integral over θ(s) and the text explicitly acknowledges that all kinematic variables are functions of θ(s). The Euler–Lagrange equations (13) and (17) are therefore not valid stationarity conditions for the ratio η^{-1} or for the speed functional. This is a mathematical error in the proof of the central claim, not merely a disagreement with numerical consensus. The paper’s restricted sawtooth-family formulas, and the asymptotic agreement with Alben’s numerical results, are real and may be salvageable, but they do not rescue the variational argument. Because this concern is the same one the reader identified, my analysis agrees with the reader’s assessment. I recommend no change to the REJECT verdict: the strongest claim is not proven by the manuscript, although the underlying optimization problem and the restricted-form results remain worthy of further study.","tokens_in":9108,"tokens_out":4474,"duration_ms":50197,"concrete_test":"Independently re-derive Eq. 13 from Eq. 12 by computing the Gateaux derivative of η^{-1} with respect to θ. Take a candidate sawtooth θ0(s) = α on (−0.5,0) and −α on (0,0.5), with α from Eq. 19 at µ = 2, and perturb it by θ_ε(s) = θ0(s) + ε sin(2πs) (adjusting a constant to preserve ∫ sinθ ds = 0). Solve the force-balance equations to first order in ε, compute δU from Eq. 9 and δ∫cosψ ds, and evaluate δ(η^{-1}) = (1/U)δ∫cosψ ds − (η^{-1}/U)δU. If this derivative is nonzero, the claimed Euler–Lagrange stationarity is false and the sawtooth is not the variational optimum of the cost functional. If it is zero, the paper must still supply the missing argument that δU = 0 for all admissible perturbations.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation 12 defines η^{-1} = (1/U[θ]) ∫ cosψ(s;θ) ds, with U[θ] itself an integral over θ(s) in Eq. 9. The paper sets L = cosψ/U and, observing that θ′ does not appear, concludes from the Euler–Lagrange equation that ∂(cosψ/U)/∂θ = 0. This treats U as a constant independent of θ. Correctly, the first variation of a ratio is δ(η^{-1}) = (1/U) δ∫cosψ ds − (η^{-1}/U) δU, and δU is not zero for general admissible perturbations. The text even concedes, immediately before this step, that “we do not know the explicit form of U in terms of θ(s); we only know of its existence,” which makes the subsequent treatment of U as a constant especially problematic. The speed-maximization argument has the same defect: Eq. 17 treats v cos(θ+ψ−ϕ) as a local function of θ, while W, θ0, ϕ, and ψ are global functionals determined by the force-balance equations. Equation 18 itself contains W(θ) and θ0(θ) as s-independent but shape-dependent quantities, so setting ∂/∂θ on the integrand is not a valid stationarity condition. In addition, the closure constraint ∫ sinθ ds = 0 is not incorporated with a Lagrange multiplier, so the variation is not restricted to the admissible class. Thus the claim that θ = ±constant, and hence sawtooth optimality, does not follow from the variational argument as written. This is a load-bearing mathematical gap, not a cosmetic one: removing it removes the central proof, even though the restricted sawtooth-family formulas may still be salvageable.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops an analytical variational framework for optimal locomotion of a snake modeled as an inextensible, flexible, periodic body moving on a flat surface under Coulomb friction with anisotropic coefficients. It claims that minimizing the cost of transport and maximizing forward speed both lead to sawtooth body shapes, i.e., θ(s)=±constant, with the optimal angle determined by the friction ratio μ. The paper derives closed-form expressions for the optimal speed and efficiency as functions of wavelength λ (Eqs. 14–16, 19), compares the results with Alben's numerical optimization, and extends the same variational argument to sand-swimming and viscous-fluid propulsion. The central mathematical step in Section 3 is a Euler–Lagrange treatment of integrands that contain no derivative of θ, leading to the conclusion that cos ψ and then cos θ must be constant.","tokens_in":9447,"tokens_out":3161,"duration_ms":34180,"significance":"If the main derivation were correct, the paper would provide a compact analytical explanation for a range of numerical observations, including the emergence of sawtooth-like optimal shapes at low transverse friction, and would unify Lighthill's flagellar result with snake locomotion and other RFT media. The manuscript has clear strengths: it presents an explicit non-dimensional model, gives closed-form expressions for the restricted sawtooth family, reports no fitted parameters, and includes a direct comparison with Alben's independent numerical optimization. However, the purported variational proof of global optimality is invalid as written, because it treats the average speed and other global kinematic quantities as constants when they are functionals of the shape. Consequently, the central claim that the optimal shape must be a sawtooth is not established by the derivation, and the paper's main contribution reduces to an analysis of a one-parameter family of candidate shapes.","major_comments":[{"comment":"The derivation of the optimal efficiency condition is not a valid variational step. Equation (12) defines η^{-1} = (1/U) ∫ cos ψ ds, and Eq. (9) defines U as an integral over v(s), which itself depends on θ(s) through Eqs. (1)–(3) and the equilibrium equations. When the shape is varied, δU is generally nonzero, so the first variation contains the term −(η^{-1}/U) δU. The paper sets ∂/∂θ(cos ψ/U)=0, which treats U as a constant independent of θ. This is explicitly contradicted by the text immediately before Eq. (13), which states that all kinematic variables, including U, are functions of θ(s), and that the explicit form of U in terms of θ(s) is unknown. Since the constancy of cos ψ is the premise for the entire sawtooth conclusion, this is a load-bearing error.","section":"§3.1, Eqs. (12)–(13)"},{"comment":"The same functional-derivative defect appears in the speed maximization. Equation (17) applies ∂/∂θ to the integrand v cos(θ+ψ−ϕ) and concludes that this quantity is constant. But v, ψ, and ϕ are not local functions of θ(s); they are determined globally by the force-balance equations (8)–(10). Equation (18) makes the dependence explicit: the right-hand side contains W and θ0, which are s-independent but shape-dependent. A correct stationarity condition must account for variations of W, θ0, ϕ, and ψ induced by a change in θ(s). The paper does not do this, so the derivation of Eq. (19) is unsupported.","section":"§3.2, Eqs. (17)–(18)"},{"comment":"The admissible class of body shapes is restricted by the continuity/closure condition ∫ sin θ ds=0, as stated in Section 3. However, the Euler–Lagrange equations (13) and (17) are derived without introducing this constraint via a Lagrange multiplier or otherwise restricting the variation space. As a result, the derived necessary conditions apply to a larger set of shapes than the admissible ones. Even if the rest of the variational reasoning were correct, optimality among admissible periodic shapes would not follow without enforcing this constraint.","section":"§3, closure constraint"},{"comment":"The comparison with Alben's numerical optimization and the asymptotic match in Eq. (20) are not sufficient to validate the optimality claim. Because the sawtooth conclusion is derived from the invalid variational steps, the agreement in Figure 3 only shows that the one-parameter sawtooth family contains shapes with favorable efficiency; it does not demonstrate that these shapes are global optimizers. The statement that the numerical and analytical results are 'closely aligned' therefore overstates what the analytical framework actually proves.","section":"§4 and Figure 3"}],"minor_comments":[{"comment":"There are several typographical errors: 'Rwynolds number' should be 'Reynolds number', 'viscus' should be 'viscous', 'Heavyside' should be 'Heaviside', and 'non-dimensionalion' should be 'non-dimensionalization'.","section":"§1 and §2, typos"},{"comment":"The notation θ=±constant is used interchangeably with 'sawtooth shape,' but a periodic function composed of segments with slopes ±θ is not uniquely a sawtooth; the paper acknowledges this, yet the phrasing 'a sawtooth shape is the most natural shape' is informal. Clarifying the exact class of shapes and how the optimal one is selected within that class would improve precision.","section":"§3.1 and §3.2, notation"},{"comment":"The caption states that the numerical results are 'reconstructed from [7]' but does not describe the reconstruction procedure or its accuracy. Since this comparison is used to support the agreement claim, more detail is needed.","section":"§4, Figure 3"},{"comment":"In Section 1, the text 'David Hu et al. [5] also demonstrated' cites reference [5], which is Lighthill's book, not the experimental study by Hu et al. The intended citation appears to be [1] or [2].","section":"References"}],"recommendation":"reject","confidential_remarks":"The manuscript contains a clear and central mathematical error: the variational derivation in Section 3 treats U, W, θ0, ϕ, and ψ as constants when they are functionals of θ(s). This is not a matter of presentation or a missing reference; it invalidates the proof of the paper's main claim. The restricted sawtooth-family formulas may be salvageable as a separate, weaker contribution, but the paper as written does not establish global optimality. I see no sign of circular reasoning or parameter fitting; the issue is a derivation error, not an empirical overreach."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper has a real contribution in the closed-form wavelength formulas and the asymptotic match to Alben, but the main proof is unsound. The Euler–Lagrange step in §3 treats U and the other global kinematic variables as independent of the shape θ(s), even though the paper itself states they are functionals of θ. That is not a cosmetic slip; it is the load-bearing step that produces the sawtooth conclusion.\n\nWhat's new: for the restricted family of sawtooth shapes, Eqs. 14, 16 and 19 are clean closed-form results, and the amplitude scaling A ~ (μ/2)^{-1/4} matches Alben's numerical asymptotics. That is valuable and non-obvious. The framework's extension to sand and viscous RFT is a nice unifying thought, though Lighthill already had the flagellar sawtooth.\n\nSoft spots: the variational argument. Eq. 12 writes η^{-1} = (1/U[θ]) ∫ cosψ ds. When minimizing over θ, δU is not zero. The first variation has an extra term −(η^{-1}/U) δU. The paper acknowledges U is unknown as a functional of θ and then promptly ignores that by treating it as constant in the Euler–Lagrange equation. Same problem in §3.2 for speed: ∂/∂θ of the integrand is not the correct stationarity condition when W, θ0, ϕ are shape-dependent. The closure constraint ∫ sinθ ds = 0 is also never enforced with a multiplier, so the variation explores inadmissible shapes. Removing this step removes the proof that θ must be piecewise constant. The restricted sawtooth optimization still stands, but the claim that the full variational problem forces a sawtooth does not.\n\nAlso, the comparison to Alben in Fig 3 is reconstructive and qualitative, not a head-to-head on the same cost function; fine as a sanity check, but not independent validation of the proof.\n\nWho it's for: researchers in snake robotics and resistive-force-theory locomotion will find the explicit formulas useful, and the failure mode is an instructive example of why one cannot minimize a ratio of functionals by treating the denominator as fixed. It deserves a serious referee, but the current version should not be accepted as is. A corrected derivation, or an explicit restriction to the sawtooth family as a conjecture supported by numerics, would make this publishable.","headline":"Useful closed-form formulas for sawtooth snake gaits, but the central variational proof treats U—a functional of θ(s)—as a constant, so the sawtooth-optimality claim does not follow as written.","tokens_in":9976,"tokens_out":1856,"would_cite":false,"duration_ms":17494,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"For snakes on flat ground, the optimal body shape for both efficiency and speed is a sawtooth, with the tooth angle set by the friction ratio.","keywords":["snake locomotion","friction anisotropy","resistive force theory","variational optimization","cost of transport","sawtooth waveform","undulatory locomotion","low Reynolds number"],"falsifier":"Pick a friction ratio such as μ=1.5, solve the full force-balance equations (Eqs. 8–10) exactly for a family of smooth sinusoidal body shapes, and compare their cost of transport and speed with the sawtooth values from Eqs. 16 and 19; if any non-sawtooth shape beats the sawtooth on either objective, the variational reduction is incorrect. The same check can be done experimentally with a snake robot whose ventral friction is tuned to μ=1.5.","tokens_in":8831,"feed_emoji":"🐍","tokens_out":7945,"duration_ms":72017,"temperature":0.7,"pith_summary":"The paper derives an analytical answer to an old question: how should a snake shape its body to slither efficiently, or as fast as possible, on flat ground? Treating the snake as a uniform, flexible, rigid-free body whose interaction with the ground is anisotropic Coulomb friction, it writes both the cost of transport and the forward speed as integrals over the body-shape function and applies the Euler-Lagrange equations. The conclusion is that the optimal body shape is a sawtooth, θ(s)=±constant, with the tooth angle fixed by the friction ratio μ=μ_n/μ_f through closed-form expressions. If correct, the same sawtooth conclusion holds for any slender-body propulsion governed by resistive force theory, including sand-swimming and viscous-fluid swimming.","feed_headline":"Optimal snake slither is a sawtooth, not a wave","feed_subtitle":"A friction-only model gives closed-form angles for max efficiency and speed—same shape for sand and fluid swimmers.","key_machinery":"The central object is the body shape function θ(s), the local tangent angle relative to the direction of travel, defined over one period of the undulation. The paper expresses both the cost of transport (Eq. 12) and the average speed (Eq. 9) as functionals of θ(s) with no θ′ dependence, so the Euler-Lagrange equations collapse to algebraic stationarity conditions rather than differential equations. For efficiency the condition is that cosψ, the component of local velocity tangent to the body, be constant; for speed it is that cos(θ−φ) be constant. The force-balance equations then select θ=±constant, i.e. a sawtooth, and the closed-form relations (Eqs. 14, 16, 19) connect the sawtooth angle and wavelength to the friction ratio μ.","core_discovery":"On the paper's own terms, the central discovery is that the efficiency-optimal and speed-optimal gaits coincide in form: both are piecewise-linear sawtooth waves with θ(s)=±constant, not smooth sinusoidal curves. The argument starts from the cost of transport $η^{{-1}}$ = (1/U)∫ cosψ ds and the speed U = ∫ v cos(θ+ψ−φ) ds (Eqs. 12 and 9), assumes the average speed U can be held fixed during shape variation, and reduces the Euler-Lagrange condition to a pointwise constancy condition on cosψ (efficiency) or on cos(θ−φ) (speed). Combined with the transverse force balance, this forces θ to be constant in magnitude. The paper then derives U = λ − λ/(λ²+μ(1−λ²)) (Eq. 14), a sixth-order polynomial in λ for the efficiency optimum (Eq. 16), and the closed-form speed-optimal wavelength λ = $\\sqrt$((2μ+1 − $\\sqrt$(8μ+1))/(2(μ−1))) (Eq. 19). The same variational step applied to granular and viscous resistive-force laws yields the same sawtooth optimum, suggesting the result is a general feature of RFT-based undulatory locomotion.","pith_inferences":["If the sawtooth claim survives scrutiny, it implies that the smooth, nearly sinusoidal body waves observed in many real snakes are not optimal in the friction-only model; the discrepancy is a plausible signature of body stiffness, actuation constraints, or finite-length effects that the model deliberately omits.","A testable extension: measure the local tangent angle of an actual snake or snake robot on surfaces with engineered friction ratios, and compare the distribution of θ to the predicted sawtooth angle; a clear concentration near the predicted angle would support the model, while broad smooth distributions would indicate that omitted costs matter.","The variational reduction to a pointwise constancy condition may generalize to other RFT settings, e.g., swimming near a wall or in a stratified fluid, where the force law changes the stationarity condition and possibly selects different optimal shapes.","Because the paper treats U as fixed in the variation, a corrected variational treatment that accounts for δU could either confirm the sawtooth or shift the optimum; re-deriving Eq. 16 with the full functional derivative is a direct mathematical check."],"forward_implications":["An optimally efficient snake on a flat surface should adopt a sawtooth waveform; for the experimentally observed range 1<μ<2, the optimal efficiency angle stays near 45°, while the speed-optimal angle is steeper at low μ.","The speed-optimal wavelength has a closed form, λ = sqrt((2μ+1 − sqrt(8μ+1))/(2(μ−1))), and its wave amplitude scales as A ∼ (1/2)(μ/2)^(−1/4) at large μ, matching the asymptotic behavior of prior numerical optimization.","Because the same stationarity condition arises for granular and viscous resistive-force laws, the sawtooth conclusion extends to sand-swimming lizards and flagellar swimming, not just snake slithering.","The analytical results provide a parameter-free upper bound for efficiency and speed in more realistic finite-length snakes, where body rigidity and unsteady average speed add dissipation."],"supporting_citations":[{"why":"Supplies the steady-wave variational method for flagellar propulsion and the viscous resistive-force model that the paper extends to snake locomotion.","marker":"[5]"},{"why":"Provides the numerical optimization baseline for snake gaits and the asymptotic wave-amplitude scaling that the analytical results are compared against.","marker":"[7]"},{"why":"Defines the anisotropic Coulomb friction model (μf, μb, μn) used in the force balance.","marker":"[2]"},{"why":"Provides the experimental friction measurements and the negligible-inertia assumption underpinning the model.","marker":"[1]"},{"why":"Supplies the granular resistive-force law for sand-swimming that the paper uses to extend the sawtooth result.","marker":"[21]"}],"fun_headline_variants":["Sawtooth, not sine: optimal snake slither solved analytically","Snakes' efficient slither is a sawtooth, new analysis shows","Optimal snake motion: piecewise-linear sawtooth beats smooth wave","Analytical model: best snake slither is a sawtooth shape"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation treats the average forward speed U as a fixed number when the body shape is varied, even though U is itself an integral over that shape (Eq. 9); if a small shape change changes U, the Euler-Lagrange step forcing cosψ to be constant does not follow.","fun_headline_variants_meta":{"raw":{"variants":["Sawtooth, not sine: optimal snake slither solved analytically","Snakes' efficient slither is a sawtooth, new analysis shows","Optimal snake motion: piecewise-linear sawtooth beats smooth wave","Analytical model: best snake slither is a sawtooth shape"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000519,"raw_usage":{"total_tokens":2507,"prompt_tokens":929,"completion_tokens":1578,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":545,"completion_tokens_details":{"reasoning_tokens":1498}},"tokens_in":545,"tokens_out":1578,"duration_ms":12719,"temperature":1.0,"reasoning_tokens":1498,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T22:09:21.123608+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Pick a friction ratio such as μ=1.5, solve the full force-balance equations (Eqs. 8–10) exactly for a family of smooth sinusoidal body shapes, and compare their cost of transport and speed with the sawtooth values from Eqs. 16 and 19; if any non-sawtooth shape beats the sawtooth on either objective, the variational reduction is incorrect. The same check can be done experimentally with a snake robot whose ventral friction is tuned to μ=1.5.","supporting_citations":[{"cited_title":"Mathematical Biofluiddynamics","cited_arxiv_id":null,"evidence_quote":"Supplies the steady-wave variational method for flagellar propulsion and the viscous resistive-force model that the paper extends to snake locomotion."},{"cited_title":"Optimizing snake locomotion in the plane","cited_arxiv_id":null,"evidence_quote":"Provides the numerical optimization baseline for snake gaits and the asymptotic wave-amplitude scaling that the analytical results are compared against."},{"cited_title":"Slithering locomotion","cited_arxiv_id":null,"evidence_quote":"Defines the anisotropic Coulomb friction model (μf, μb, μn) used in the force balance."},{"cited_title":"Hu, Jasmine Nirody, Terri Scott, and Michael J","cited_arxiv_id":null,"evidence_quote":"Provides the experimental friction measurements and the negligible-inertia assumption underpinning the model."},{"cited_title":"Maladen, Yang Ding, Chen Li, and Daniel I","cited_arxiv_id":null,"evidence_quote":"Supplies the granular resistive-force law for sand-swimming that the paper uses to extend the sawtooth result."}],"review_version":1}