{"id":"978ac833-5e4d-4020-a06c-3e09958aaad0","arxiv_id":"2608.01491","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"For higher-derivative Lagrangians, Ostrogradsky's canonical coordinates and momenta are exactly the coefficients of the variation terms in the boundary term of delta L.","lead":"This paper shows undergraduates that the canonical momenta of Lagrangians with higher time derivatives can be read from the boundary term in the first variation of the action. It derives Ostrogradsky's Hamiltonian formulas from repeated integration by parts and provides code and worked examples.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Boundary-term uniqueness claim in Sec. III is false: rescaling Q and P preserves PδQ, so Eq. (56) does not force Eq. (57); Ostrogradsky formulas remain correct as a convention.","rationale":"The central derivation of the Ostrogradsky momenta in Section IV is algebraically correct: Eq. (59) is valid, the summation reordering is sound, and identifying Q_i=q^(i-1) with P_i as the coefficient in Θ reproduces Woodard's formulas. The non-degeneracy condition (12) is exactly what is needed for the Legendre transform. The genuine weakness is in Section III, where the paper claims Eq. (56) forces the canonical pair to be (q,∂L/∂qdot). This is false: for L=1/2 qdot^2, (Q,P)=(2q,qdot/2) gives PδQ=qdotδq and valid Hamilton equations with H=2P^2. The boundary term therefore determines the symplectic potential only up to exact forms and coordinate rescalings. This does not break Eq. (69)—if one chooses Q_i=q^(i-1) as the configuration variables, the coefficient in Θ is exactly Ostrogradsky's P_i. But the paper's language that the boundary term 'defines' or 'determines' the canonical pairs is too strong, and the derivation is better presented as a natural construction with an explicit convention. The reader's conditional verdict captures this appropriately; I would keep it.","tokens_in":8095,"tokens_out":17982,"duration_ms":205151,"concrete_test":"For L=1/2 qdot^2, compute the two pairs (Q,P)=(q,qdot) and (Q,P)=(2q,qdot/2). Verify that both satisfy Eq. (56) (both give P δQ = qdot δq) and that each yields valid Hamilton equations equivalent to qddot=0. If both pass, Section III's uniqueness claim is refuted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III's inference from Eq. (56) to Eq. (57) is invalid. Eq. (56) only equates two total time derivatives; the boundary one-forms may differ by an exact form, and the canonical coordinates themselves are not fixed. Concretely, for L=1/2 qdot^2, both (Q,P)=(q,qdot) and (Q,P)=(2q,qdot/2) satisfy Eqs. (53)-(56) and give the same boundary term P δQ = qdot δq, with Hamiltonians H=1/2 P^2 and H=2P^2, respectively. Thus the boundary term does not single out the identification (57). The same non-uniqueness (rescaling Q_i -> c_i Q_i, P_i -> P_i/c_i, or adding exact terms) carries to the higher-derivative case, so the statement that the boundary term 'defines' the canonical pairs overstates what Eq. (56) proves. The Ostrogradsky formulas (69) remain correct for the natural choice Q_i=q^(i-1), but that choice is a convention, not a consequence of the boundary term alone.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a derivation of Ostrogradsky's Hamiltonian formalism for higher-derivative Lagrangians starting from the boundary term in the first variation. Section IV computes the boundary term for L(q, q^(1), ..., q^(N)) by repeated integration by parts and reads off canonical pairs Q_i = q^(i-1), P_i = Σ_{j≥i} (-d/dt)^{j-i} ∂L/∂q^(j), exactly the Ostrogradsky momenta. Worked examples (perturbed harmonic oscillator, Pais–Uhlenbeck oscillator), multi-variable formulas, and supplementary computer code are provided. The central integration-by-parts computation is correct and reproduces the standard Ostrogradsky results. Section III attempts to justify why the boundary term 'defines' canonical pairs; this justification contains an overclaim that is not load-bearing for the algebraic formula but is load-bearing for the paper's stated pedagogical thesis.","tokens_in":8432,"tokens_out":8955,"duration_ms":87560,"significance":"If the presentation is corrected, the paper offers a genuinely useful undergraduate-level route to the Ostrogradsky construction and makes explicit a connection between boundary terms and canonical momenta that is often left implicit. The derivation of the boundary term (Eqs. (59)–(68)) is self-contained and correct, the examples are accurate, and the accompanying code is a valuable supplement. The main weakness is Section III's claim of uniqueness: the boundary term does not uniquely determine canonical pairs, since rescalings Q_i → c_i Q_i, P_i → P_i/c_i preserve P_i δQ_i. The Ostrogradsky formulas are correct as a standard and natural convention, but not as a forced consequence of Eq. (56). This is a local, fixable gap rather than an error in the main computation.","major_comments":[{"comment":"The assertion that Eq. (56) 'is only satisfied if' (Q_i, P_i) = (q_i, ∂L/∂q̇_i) is false. For L = ½ q̇², take Q = 2q, P = q̇/2. Then P δQ = q̇ δq, so Eq. (56) holds, and with H = 2P² Hamilton's equations are equivalent to q̈ = 0. More generally, any rescaling Q_i → c_i Q_i, P_i → P_i/c_i preserves P_i δQ_i, and adding an exact form changes the boundary term by a total derivative without affecting the variational principle. Thus Eq. (56) fixes the momenta only after one chooses Q_i = q^(i-1). That choice is a convention, not a consequence. Please either prove uniqueness under a stated normalization or reframe the claim as selecting the natural/standard canonical pair.","section":"Section III, Eqs. (56)–(57)"},{"comment":"The identification of Θ with P_i δQ_i and the subsequent reading off of Eq. (69) rely on the flawed uniqueness argument in Section III. In the higher-derivative case the same rescaling freedom exists: with Q_i = c_i q^(i-1) and P_i = P_i^O/c_i, the boundary term is unchanged. The paper should state explicitly that Eq. (69) is the canonical choice obtained by taking Q_i = q^(i-1), which is the standard Ostrogradsky convention. The boundary term alone does not select this pair; the derivation yields the correct formulas only after that coordinate choice is made.","section":"Section IV, Eq. (69)"}],"minor_comments":[{"comment":"The notation L on both sides is confusing: the left side is the configuration-space Lagrangian and the right side is the phase-space Lagrangian of Eq. (49). Also, if the goal is to discuss general canonical transformations, the equality in Eq. (53) is stronger than necessary; the two Lagrangians may differ by a total derivative. This is related to the overclaim in Eqs. (56)–(57).","section":"Section III, Eq. (53)"},{"comment":"The claim that 'no derivation of Ostrogradsky's construction exists in the modern literature which emphasizes the role of the boundary term' is too strong, given that Refs. [5–8] (Barnich–Henneaux–Schomblond, Lee–Wald, Crnkovic–Witten, Torre) already use boundary terms in covariant phase-space constructions. The paper's contribution is better framed as an elementary, self-contained exposition for mechanics rather than the first boundary-term-based derivation.","section":"Section I and Discussion"},{"comment":"The notation δF/δq^(s) for the variational derivative may be confused with the variation δF. Consider denoting this quantity by E_s(F) or a similar symbol, to avoid confusion with the variation of F.","section":"Section IV, Eq. (70)"},{"comment":"The manuscript mentions a computer program and a GitHub repository but gives little detail about the algorithm or its testing. If space permits, a short appendix describing the algorithm's handling of the non-degeneracy condition (12) and a few test cases would strengthen the reproducibility claim.","section":"Supplementary material"}],"recommendation":"major_revision","confidential_remarks":"This is a solid pedagogical paper for a journal like American Journal of Physics or European Journal of Physics. The main technical derivation is correct, and the examples are well chosen. The Section III uniqueness claim is the one substantive flaw; it is fixable by softening the claim and explicitly presenting the canonical-pair choice as a convention. I recommend major revision rather than rejection. The authors should also be asked to acknowledge the prior boundary-term-based approaches in Refs. [5–8] more accurately, since that affects the novelty assertion but not the correctness of the central computation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this paper is a solid teaching piece, not a research advance. It derives Ostrogradsky's momenta by reading them off the boundary term in the first variation, with worked examples and Python/Maple code. The main derivation in Sec. IV is correct and clearly explained. I'd be happy to assign it to an undergrad.\n\nWhat's new is the ordering and presentation. The boundary-term formula (63) and the momenta (69) are standard in the covariant phase space literature [5-8] and in Woodard's review [13]. The authors are upfront about that. Their contribution is making it elementary, plus the code. That's a legitimate contribution for physics.ed-ph.\n\nThe soft spot is Sec. III. Equation (56) only equates two total time derivatives; it does not force (57). For L = 1/2 qdot^2, both (Q,P) = (q,qdot) and (2q,qdot/2) give the same boundary term, with different Hamiltonians. So the claim that the boundary term 'defines' the canonical pair is too strong. The correct statement is that it picks a natural convention, and Ostrogradsky's formulas are that convention. This doesn't hurt the main derivation, but it should be fixed. The non-degeneracy condition (12) is explicitly assumed and the degenerate case excluded, which is fine.\n\nMinor: the claim that no modern derivation exists is a bit strong, given the covariant phase space references derive it at a more advanced level. Not a big deal. The citation pattern looks honest: Woodard, the covariant phase space papers, and recent PU-oscillator references are all relevant, and the self-citation [5] is appropriate.\n\nThis paper is for instructors and advanced undergrads, not for researchers. It deserves peer review; with the Sec. III wording fixed, I'd accept it. The code and examples make it reproducible. Send it to a referee.","headline":"A clean pedagogical derivation of Ostrogradsky's construction from the boundary term, worth teaching from, but the Sec. III uniqueness claim overreaches and should be softened.","tokens_in":8815,"tokens_out":1617,"would_cite":false,"duration_ms":15262,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["45.20.Jj"],"model":"deepseek-v4-flash","headline":"Reading the boundary term in the first variation of a higher-derivative Lagrangian reproduces Ostrogradsky's canonical momenta and Hamiltonian without any ad hoc construction.","keywords":["Ostrogradsky formalism","higher-derivative Lagrangians","first variational formula","boundary term","canonical momenta","Hamiltonian mechanics","Pais–Uhlenbeck oscillator","Ostrogradsky instability"],"falsifier":"Compute $\\Theta(\\delta q)$ for a non-degenerate higher-derivative Lagrangian such as $L = \\frac{1}{2}\\ddot q^2 - \\frac{1}{2}\\omega^2 q^2$, read off $Q_1=q$, $Q_2=\\dot q$, $P_1=-q^{(3)}$, $P_2=\\ddot q$, and build $H=P_1Q_2+P_2A-L$ with $A=P_2$. Then check whether Hamilton's equations are equivalent to the Euler-Lagrange equation $q^{(4)}+\\omega^2 q=0$. A single non-degenerate Lagrangian for which this equivalence fails would disprove the paper's central claim.","tokens_in":8048,"feed_emoji":"⚙️","tokens_out":14445,"duration_ms":110999,"temperature":0.7,"pith_summary":"This paper aims to show that the Ostrogradsky Hamiltonian formalism for higher-derivative Lagrangians follows from one structural fact: the boundary term in the first variation of the Lagrangian already contains the canonical pairs. For a Lagrangian $L(q, q^{(1)}, \\ldots, q^{(N)}, t)$, the variation splits into an Euler-Lagrange expression plus the time derivative of a boundary term $\\Theta(\\delta q)$. The coefficients of $\\delta q, \\delta \\dot{q}, \\ldots, \\delta q^{(N-1)}$ in $\\Theta$ are exactly the momenta conjugate to $Q_i = q^{(i-1)}$, and they match Ostrogradsky's definitions. This turns a recipe that usually appears clever and ad hoc into a mechanical consequence of integration by parts, and it gives students a direct route to Hamiltonians for systems such as the Pais–Uhlenbeck oscillator. If the derivation is right, the awkward-looking Ostrogradsky momenta need not be memorized; they can be read off the variation.","feed_headline":"Boundary term alone yields Ostrogradsky's canonical pairs","feed_subtitle":"In this paper, higher-derivative Hamiltonians come from the first variation's boundary term, not from ad hoc Ostrogradsky formulas.","key_machinery":"The central object is the boundary term $\\Theta(\\delta q)$ in the first variational formula. The paper isolates it with the integration-by-parts identity $F\\,\\delta q^{(k)} = (-1)^k \\frac{d^k F}{dt^k}\\delta q + \\frac{d}{dt}\\left(\\sum_{j=0}^{k-1} (-1)^j \\frac{d^j F}{dt^j}\\delta q^{(k-j-1)}\\right)$, applied to each term $\\frac{\\partial L}{\\partial q^{(k)}}\\delta q^{(k)}$. This produces $\\Theta$ as a sum already arranged like $\\sum_i P_i \\delta Q_i$, so the canonical coordinates and momenta are read off directly. Equivalently, the momenta are shifted variational derivatives $P_i = \\delta L/\\delta q^{(i)}$ with $\\delta F/\\delta $q^{{(s)}}$ = \\sum_{i=0}^{N-s}\\left(-\\frac{d}{dt}\\right)^i \\frac{\\partia","core_discovery":"The paper's central claim is that the first variational formula $\\delta L = E(L)\\delta q + \\frac{d}{dt}\\Theta(\\delta q)$ contains the entire canonical structure of a higher-derivative mechanical system. For a non-degenerate Lagrangian depending on derivatives up to order $N$, the boundary term is $\\Theta(\\delta q) = \\sum_{i=1}^N \\left(\\sum_{j=i}^N \\left(-\\frac{d}{dt}\\right)^{j-i} \\frac{\\partial L}{\\partial q^{(j)}}\\right)\\delta q^{(i-1)}$, which has the form $\\sum_i P_i \\delta Q_i$ with $Q_i = q^{(i-1)}$ and $P_i = \\sum_{j=i}^N \\left(-\\frac{d}{dt}\\right)^{j-i} \\frac{\\partial L}{\\partial q^{(j)}}$. These $P_i$ are precisely Ostrogradsky's momenta, so the Hamiltonian $H = \\sum_{i=1}^{N-1} P_i","pith_inferences":["The same boundary-term logic could be applied to degenerate higher-derivative Lagrangians: $\\Theta(\\delta q)$ is still defined even when $\\partial^2 L/\\partial (q^{(N)})^2 = 0$, so it may constrain the constrained Hamiltonian analysis, but the paper does not pursue that.","Reading momenta as shifted variational derivatives, $P_i = \\delta L/\\delta q^{(i)}$, suggests a uniform computational rule for multi-degree-of-freedom and field-theory generalizations that could be tested on textbook Lagrangians.","If the boundary term is taken as the primary object, the Legendre transform becomes secondary: one could in principle derive Hamilton's equations directly from $\\Theta$ without first constructing $H$, an extension the paper leaves implicit."],"forward_implications":["For $N=1$, $\\Theta = (\\partial L/\\partial \\dot q)\\delta q$, so the boundary-term rule reduces to the standard momentum definition and the usual Hamiltonian.","For $N=2$, it recovers $P_1 = \\partial L/\\partial \\dot q - \\frac{d}{dt}\\partial L/\\partial \\ddot q$ and $P_2 = \\partial L/\\partial \\ddot q$, the Ostrogradsky pairs used in the worked examples.","The boundary term organizes the $2N$ initial data of a generically $2N$-th-order Euler-Lagrange equation into $N$ canonical pairs, explaining the dimension of Ostrogradsky's phase space.","Because the Hamiltonian is linear in at least one canonical momentum for every non-degenerate higher-derivative Lagrangian, the construction reproduces the Ostrogradsky instability: the Pais–Uhlenbeck example has no ground state.","The accompanying code implements the derivation, so a Hamiltonian for a given non-degenerate higher-derivative Lagrangian can be produced by algorithm rather than by hand."],"supporting_citations":[{"why":"Supplies the standard Euler-Lagrange equation for higher derivatives and the textbook boundary-term context.","marker":"[1]"},{"why":"Gives the field-theory first-variational formula with which the paper's boundary term is consistent.","marker":"[6]"},{"why":"Provides the Pais–Uhlenbeck oscillator Lagrangian used as the main worked example.","marker":"[11]"},{"why":"Ostrogradsky's original construction, which the paper re-derives from the boundary term.","marker":"[12]"},{"why":"Modern review of the Ostrogradsky formulas and instability; the paper's Eq. (69) is checked against it.","marker":"[13]"}],"fun_headline_variants":["Variation's boundary term births Ostrogradsky's momenta","First variation alone builds higher-derivative Hamiltonians","Ostrogradsky from the boundary: a single formula","Undergrad-level derivation: Ostrogradsky from first variation"],"cache_read_input_tokens":2816,"weakest_assumption_plain":"The load-bearing premise is that the Lagrangian depends non-degenerately on its highest derivative, $\\partial^2 L/\\partial (q^{(N)})^2 \\neq 0$, so that $q^{(N)}$ can be solved for in terms of the canonical variables and the Hamiltonian can be built; degenerate higher-derivative Lagrangians are explicitly outside the paper's scope.","fun_headline_variants_meta":{"raw":{"variants":["Variation's boundary term births Ostrogradsky's momenta","First variation alone builds higher-derivative Hamiltonians","Ostrogradsky from the boundary: a single formula","Undergrad-level derivation: Ostrogradsky from first variation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000206,"raw_usage":{"total_tokens":1213,"prompt_tokens":705,"completion_tokens":508,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":449,"completion_tokens_details":{"reasoning_tokens":438}},"tokens_in":449,"tokens_out":508,"duration_ms":4987,"temperature":1.0,"reasoning_tokens":438,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T00:04:51.761371+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute $\\Theta(\\delta q)$ for a non-degenerate higher-derivative Lagrangian such as $L = \\frac{1}{2}\\ddot q^2 - \\frac{1}{2}\\omega^2 q^2$, read off $Q_1=q$, $Q_2=\\dot q$, $P_1=-q^{(3)}$, $P_2=\\ddot q$, and build $H=P_1Q_2+P_2A-L$ with $A=P_2$. Then check whether Hamilton's equations are equivalent to the Euler-Lagrange equation $q^{(4)}+\\omega^2 q=0$. A single non-degenerate Lagrangian for which this equivalence fails would disprove the paper's central claim.","supporting_citations":[{"cited_title":"Whittaker,A Treatise on the Analytical Dynamics of Particles and Rigid Bodies: With an Introduction to the Problem of Three Bodies, 2nd ed","cited_arxiv_id":null,"evidence_quote":"Gives the field-theory first-variational formula with which the paper's boundary term is consistent."},{"cited_title":"Woodard, Avoiding dark energy with 1/rmodifications of gravity, inThe Invisible Universe: Dark Matter and Dark Energy(Springer Berlin Heidelberg, 2007) pp","cited_arxiv_id":null,"evidence_quote":"Provides the Pais–Uhlenbeck oscillator Lagrangian used as the main worked example."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Ostrogradsky's original construction, which the paper re-derives from the boundary term."},{"cited_title":"Pavˇ siˇ c, International Journal of Geometric Methods in Modern Physics13, 1630015 (2016)","cited_arxiv_id":null,"evidence_quote":"Modern review of the Ostrogradsky formulas and instability; the paper's Eq. (69) is checked against it."}],"review_version":1}