{"id":"3ea3ec6e-85b0-4f8d-8b63-4443223568ea","arxiv_id":"1908.03352","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":3.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":1,"one_line_summary":"Explicit sub-Riemannian geodesics for the Heisenberg approximation of the rolling disc, with the first cut point derived from the isotropy symmetry.","lead":"This paper applies geometric control theory and symmetry analysis to the classical vertical rolling disc, deriving explicit geodesic formulas for the Heisenberg-group nilpotent approximation. It shows how a rotational symmetry locates the first cut point and demonstrates all computations in Maple, making a useful worked example for engineers and students.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Cut-point claim is proven only for the Heisenberg model and only as an upper bound; transfer to the rolling disc lacks estimates.","rationale":"The reader's conditional verdict is appropriate. The algebra and symmetry derivations check out: the vertical systems are correct for the stated Lie brackets, the nilpotent approximation is the standard Heisenberg model, and the sl(3,R) symmetry algebra with the no-invariant-metric lemma is a known flat parabolic-geometry fact, supported here by a Tanaka prolongation computation in Maple. The weak spot is the interpretation of the cut-point result. The paper itself hedges by saying the observation 'reflects known results' and by comparing only numerically, but Section 4 presents this as determining the first cut point. The missing pieces are, first, an equality argument in the nilpotent model showing that no earlier cut time occurs, and second, any rigorous bound relating the original and nilpotent cut loci. Both are checkable. Because these are gaps in justification rather than demonstrated errors, and because the main theorems are explicitly about the nilpotent and Lagrangian structures, the conditional verdict remains appropriate and no verdict change is needed.","tokens_in":27210,"tokens_out":10545,"duration_ms":120120,"concrete_test":"Integrate the original system (17)-(18) for Example 1 (h1(0)=1/2, h2(0)=sqrt(3)/2, h3(0)=2, origin) together with the variational equation of the endpoint map. Find the first time T* where a nontrivial Jacobi field vanishes at 0 and T*. If T* differs from 2*pi/2 = pi by more than a few percent, or the endpoint has x^2+theta^2 not close to 0, the nilpotent cut-point prediction is not quantitatively reliable for the rolling disc. Repeat for h3(0)=20 with the corresponding interval.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the passage from Prop. 4.4/Cor. 4.5 to a statement about the vertical rolling disc. In the nilpotent model, the isotropy flow of t0 gives, for each generic geodesic (30), a one-parameter family of distinct local minimizers with the same endpoint at t = 2*pi/C1. That proves the geodesic cannot be optimal after that time, i.e. the cut time is at most 2*pi/C1. It does not prove equality: there could be another minimizer reaching some earlier point of the same geodesic, and the paper supplies no conjugate-point or distance comparison ruling this out before invoking known results. More importantly, t0 and the set S are symmetries of the nilpotent structure only; Remark 4.2 states that the original rolling-disc structure has no such isotropy symmetry. Section 3.4 compares two trajectories numerically but gives no error estimate, and the comparison interval ends exactly at the predicted cut time. Since the cut locus is a global object, the standard nilpotent approximation does not automatically transfer cut times from the Heisenberg model to the disc. Without a quantitative estimate or a statement that the results are only about the nilpotent model, the central geometric conclusion is not established for the disc.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the vertical rolling disc as a nonholonomic control system, derives the Pontryagin extremal equations for the original system and for its homogeneous nilpotent approximation, and then focuses on the Heisenberg nilpotent model. It gives explicit solutions of the extremal equations (Propositions 3.7 and 3.8), computes the full infinitesimal symmetry algebra of the sub-Riemannian structure (Proposition 4.1), and uses the isotropy symmetry t0 to construct a one-parameter family of local minimizers with a common endpoint in the fixed-point set S (Proposition 4.4 and Corollary 4.5). The paper then studies the Lagrangian contact structure of the nilpotent model, proves that its symmetry algebra is sl(3,R) (Proposition 5.2), and shows that no positive definite sub-Riemannian metric is invariant under this full symmetry algebra (Lemma 5.4). The text is accompanied by extensive Maple/DifferentialGeometry code and several numerical comparisons with the original rolling-disc system.","tokens_in":27423,"tokens_out":30230,"duration_ms":273342,"significance":"The paper is a careful, self-contained computation in the geometric control of the Heisenberg nilpotent model associated with the vertical rolling disc. The explicit solutions in Proposition 3.8, the symmetry classification in Propositions 4.1 and 5.2, and the non-existence result in Lemma 5.4 are cleanly stated and are supported by reproducible CAS code and direct verification of the differential equations. The symmetry-based construction in Section 4.3 gives a transparent geometric explanation of the known Heisenberg conjugate/cut-locus results, and the discussion of the metric class attached to the Lagrangian contact structure is a useful contribution. The main limitation is that the cut-point and symmetry results are proven for the nilpotent approximation only; the connection to the original rolling disc is made through numerical examples and informal statements, not through quantitative estimates.","major_comments":[{"comment":"The proof establishes that, for each generic extremal (30), the isotropy flow produces a one-parameter family of distinct local minimizers with the same endpoint at time t=2*pi/C1, which shows that the geodesic cannot be optimal after that time. This is an upper bound on the cut time. The phrase \"and it is the first point with this property\" refers to the first intersection with the fixed-point set S, not the first cut point. If the authors intend Proposition 4.4 and Corollary 4.5 to assert equality with the cut point, an additional argument (for example, a conjugate-point or distance comparison) is required; as written, the equality is only supported by the citations [2,26] mentioned after Corollary 4.5. Please make the exact logical status of the cut-point claim explicit.","section":"Section 4.3, Proposition 4.4"},{"comment":"The numerical comparison in Section 3.4 is purely visual and contains no error estimate, and Remark 4.2 explicitly states that the original rolling-disc structure has no isotropy symmetry analogous to t0. Consequently, the cut-point and symmetry results in Section 4 are established only for the nilpotent Heisenberg model, not for the vertical rolling disc itself. The title, the abstract, and the sentence in Section 3 saying that the nilpotent approximation \"still describes the system appropriately\" should be qualified to reflect this, or a quantitative transfer statement must be supplied. As written, a reader could reasonably infer a cut-locus statement for the original system that is not proved.","section":"Sections 3.4 and 4.3 (scope)"}],"minor_comments":[{"comment":"The condition C2*C3 different from 0 is stronger than needed: the first-intersection computation with S also works when exactly one of C2, C3 vanishes, and the genuinely degenerate case to exclude is C2=C3=0. In addition, the time 2*pi/C1 should be written as 2*pi/|C1|, or the authors should state explicitly that C1>0 is assumed.","section":"Proposition 4.4"},{"comment":"The displayed formula for tau_s contains (dx - s/2 dy)^2, but the flow is x maps to x + s y/2, so the pullback of r=dx^2+d theta^2 gives (dx + s/2 dy)^2. Please check this sign and also clarify the convention used for f_s^* on vector fields, since the factors w^2/4 and 4/w^2 appear inverted relative to the usual pushforward convention.","section":"Section 5.4"},{"comment":"The displayed formula for Ye7 is missing the final basis vector: it should read Ye7 = 2x^2*dx + 2xy*dy + (2y - x*theta)*d_theta.","section":"Appendix A"},{"comment":"The bracket table appears to use the convention that the entry in row i and column j represents [e_j, e_i] rather than the more common [e_i, e_j]. Please state the convention explicitly or transpose the table to avoid confusion.","section":"Table 2"},{"comment":"The text contains the typo \"Im other words\"; it should read \"In other words\".","section":"Section 1.1"}],"recommendation":"minor_revision","confidential_remarks":"The paper is a competent and useful computation on the Heisenberg nilpotent model, but its packaging overstates the connection to the original vertical rolling disc. The cut-point part largely recovers known results for the Heisenberg group, so the novelty lies mainly in the symmetry computations and the metric discussion in Section 5. The authors should be asked to clearly delimit the scope of the cut-point claims and to fix the sign/convention typos; with those changes the paper would be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear X,\n\nRead the Hrdina–Navrat–Zalabova paper. Quick take: it's a clean, honest methods paper about the Heisenberg nilpotent approximation of the vertical rolling disc, not a new-theorem paper. Most of the mathematical content is a re-derivation of known results — the nilpotent approximation, the explicit geodesics, the isotropy symmetry and the cut point are all in [2] and [26], and the sl(3,R) symmetry algebra of the Lagrangian contact structure is standard parabolic geometry. The paper says so. Its value is in the explicit, self-contained derivations and the Maple/DifferentialGeometry code, which is real reproducible evidence.\n\nWhat it does well: PMP setup is correct, vertical and horizontal systems check out, solutions (29)-(30) are explicit and satisfy the equations, the symmetry computation via Tanaka prolongation in Appendix A is detailed, and Lemma 5.4 is neatly proven. The numerical comparison in Section 3.4 is honest — two examples, no error estimates, but presented as demonstration. The attribution of the cut point to [2,26] is exactly right.\n\nSoft spots: Two. First, the cut-point discussion in Section 4 is only about the nilpotent model, and Prop 4.4 alone establishes at most an upper bound on the cut time: the orbit of the isotropy symmetry gives a distinct local minimizer with the same endpoint at t=2π/C1, which means optimality ceases at or before that time; equality requires the known conjugate-locus results, which the paper cites but could be more explicit about when presenting Prop 4.4 as a 'recovery' of the cut point. Second, the transfer from the nilpotent model to the actual rolling disc is informal. The paper never proves that the cut locus of the approximation controls the cut locus of the disc, and the two numerical comparisons in Section 3.4 don't amount to an estimate. The paper is careful not to overclaim — Remark 4.2 notes the original structure lacks the isotropy symmetry — but the title and framing might lead a reader to think the cut-point result applies to the disc. One sentence of scope limitation would fix it.\n\nWho is this for: someone wanting a worked, executable example of symmetry methods in sub-Riemannian control, or checking Maple implementations. It deserves peer review — it's sound, self-contained, and reproducible — but expect the referee to ask for a clearer statement of what is genuinely new and a tighter discussion of the nilpotent approximation's limits.\n\nFWIW, I'd send it out rather than desk-reject.","headline":"A sound, honest methods paper that re-derives known Heisenberg results for the rolling disc; the useful contribution is the reproducible Maple computation, not a new theorem.","tokens_in":27971,"tokens_out":6871,"would_cite":false,"duration_ms":69263,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["53C17","70Q05","22E60","37J60"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper proves that each generic geodesic in the nilpotent rolling-disc model first stops being optimal exactly where a rotational symmetry fixes it, and gives the location explicitly.","keywords":["vertical rolling disc","nilpotent approximation","Heisenberg group","sub-Riemannian geometry","cut point","isotropy symmetry","Lagrangian contact structure","maximum principle"],"falsifier":"For the nilpotent system with $C_1=2$ and $C_2^2+C_3^2=1$ (as in Example 1), the proposition predicts the first cut point $(0,\\pi/4,0)$ at time $\\pi$; solve the original rolling-disc extremal system numerically with the same initial covector and test whether any admissible curve shorter than the extremal connects the origin to that point before time $\\pi$. If such a shorter curve exists, or if the first intersection of the true extremal with $x=\\theta=0$ occurs at a value of $y$ different from the nilpotent prediction, the approximation does not control the original cut locus.","tokens_in":26975,"feed_emoji":"🎯","tokens_out":12556,"duration_ms":110332,"temperature":0.7,"pith_summary":"The paper studies optimal movement of a vertical rolling disc as a sub-Riemannian control problem and replaces the intractable extremal equations by a homogeneous nilpotent approximation modelled on the three-dimensional Heisenberg group. In that model every generic length-minimizing curve can be written explicitly, and the paper proves that a one-dimensional rotational symmetry fixes the exact point where each such curve first ceases to be optimal: the first cut point. The argument identifies this point directly from the symmetry's fixed-point set, without solving a separate optimality problem. The paper further proves that the Lagrangian contact structure, which keeps angular and plane velocities distinguished, has symmetry algebra $\\mathfrak{sl}(3,\\mathbb{R})$, and that no positive-definite sub-Riemannian metric is invariant under that full symmetry group. Explicit cut points matter because they are what make sub-Riemannian distance and motion planning tractable on the model that approximates the rolling disc.","feed_headline":"Symmetries locate first cut point of rolling-disc geodesics","feed_subtitle":"Every generic minimizer in the Heisenberg model first stops being optimal at a point fixed by a rotational symmetry.","key_machinery":"The load-bearing object is the isotropy symmetry $t_0$, a vector field whose flow rotates the horizontal distribution around the origin and whose fixed-point set is the line $S=\\{(0,y,0)\\}$. Applying $t_0$ to a local minimizer from (30) maps it to other local minimizers from the origin with the same parametrization and the same endpoint at time $2\\pi/C_1$; a symmetry with this property gives a direct certificate that the curve stops being optimal at the meeting point. Around this, the machinery is the left-invariant Hamiltonian description of the Heisenberg-group control problem, which splits extremals into the vertical system (27) and horizontal system (28), and the translation symmetries $t_1,t_2,t_3$ that move solutions between different initial conditions. For the Lagrangian contact structure, the algebraic prolongation of the grading $g_{-2}\\oplus g_{-1}$ yields the full symmetry algebra $\\mathfrak{sl}(3,\\mathbb{R})$, and the action of its degree-zero part on $g_{-1}$ rules out any invariant positive-definite metric.","core_discovery":"The central result is Proposition 4.4. In the nilpotent approximation with coordinates $(\\theta,x,y)$, every local minimizer from (30) with $h_3=C_1\\neq 0$ and $C_2C_3\\neq 0$ intersects the fixed-point set $S=\\{(0,y,0)\\}$ of the isotropy symmetry $t_0 = \\theta\\,\\partial_x + \\frac{\\theta^2-x^2}{2}\\,\\partial_y - x\\,\\partial_\\theta$ at the point $(0,\\pi(C_2^2+C_3^2)/C_1^2,0)$ at time $2\\pi/C_1$, and this is the first such intersection. Since two equal-length extremals meeting at the same point and time cannot be optimal past the meeting point, this intersection is the first cut point along the minimizer; for arc-length parametrized curves, where $C_2^2+C_3^2=1$, it simplifies to $(0,\\pi/C_1^2,0)$. The isotropy flow simultaneously produces a one-parameter family of equal-length minimizers from the origin to that point. The paper also establishes, by algebraic prolongation of the contact grading, that the symmetry algebra of the Lagrangian contact structure is $\\mathfrak{sl}(3,\\mathbb{R})$, and that no positive-definite sub-Riemannian metric is invariant under that algebra.","pith_inferences":["A natural but unproved conjecture is that the first cut time of the original rolling disc converges to $2\\pi/C_1$ as the endpoint approaches the origin; the paper's two numerical examples are consistent with this, but Section 3.4 gives no error estimate.","The same symmetry argument should apply to any control system with the same three-dimensional solvable controllability algebra, such as the reduced kinematic-car system (11), giving the same first cut-point formula after a coordinate change.","The $\\mathfrak{sl}(3,\\mathbb{R})$ symmetry suggests a motion-planning shortcut: instead of re-solving the extremal equations for each choice of velocity-unit ratio, one can pull known Heisenberg minimizers back through the symmetry flow and compare their cost under the family of metrics $\\tau_s$.","A sharper structural test would be to reduce the Heisenberg geodesic flow by the isotropy symmetry $t_0$; the fixed-point set $S$ should then appear as a caustic or conjugate locus in the reduced space, giving a symplectic explanation of the first cut time."],"forward_implications":["For the nilpotent approximation, every generic arc-length minimizer has its first cut point on the line $S$ at $(0,\\pi/C_1^2,0)$, so the first cut locus is explicitly parametrized by $C_1$.","Because the isotropy flow generates a one-parameter family of equal-length minimizers between the origin and each cut point, the cut locus near the origin is the rotational image of a single generic minimizer.","The translation symmetries mean a solution starting anywhere in a neighbourhood can be obtained from the origin solution by the group law with the same vertical covector, so the explicit formulas (30) solve every nearby optimal-control problem, not only those starting at the origin.","The symmetry algebra $\\mathfrak{sl}(3,\\mathbb{R})$ has maximal possible dimension for a three-dimensional Lagrangian contact structure, so the nilpotent model is flat in the parabolic-geometry sense; all metrics that preserve the angular/plane velocity splitting but vary the unit ratio are related by these symmetries.","Since no positive-definite metric is invariant under $\\mathfrak{sl}(3,\\mathbb{R})$, any choice of units for angular versus plane velocity breaks the full symmetry, leaving exactly the three translations as the metric-preserving symmetry group."],"supporting_citations":[{"why":"Supplies the maximum principle and the Hamiltonian system used to derive the extremal equations.","marker":"[1]"},{"why":"Supplies the sub-Riemannian framework, left-invariant Hamiltonians, and the definition of cut points used in Section 4.","marker":"[2]"},{"why":"Provides the algebraic construction of homogeneous nilpotent approximation used to build the Heisenberg model.","marker":"[4]"},{"why":"Supplies the nilpotent approximation method for control distributions expressed in privileged coordinates.","marker":"[17]"},{"why":"Defines Lagrangian contact structures and the flat model whose symmetry bound Proposition 5.2 matches.","marker":"[8]"},{"why":"Provides the explicit Jacobi-elliptic solutions of the original extremal system that motivate the nilpotent approximation.","marker":"[30]"},{"why":"Introduces the algebraic prolongation of graded Lie algebras used in the appendix to compute the symmetry algebra.","marker":"[32]"},{"why":"Supplies the classification of differential systems associated with simple graded Lie algebras, used in identifying the prolonged algebra as $\\mathfrak{sl}(3,\\mathbb{R})$.","marker":"[33]"},{"why":"Provides the prolongation procedure for filtered structures of constant type that yields the full geometric symmetry algebra.","marker":"[34]"}],"fun_headline_variants":["Symmetry fixes first cut point of rolling-disc minimizers","Isotropy symmetry pinpoints rolling-disc cut point","Geometric control: symmetry locates first cut point","Rolling disc: symmetry reveals exact cut point","Symmetry algebra locates first cut point in rolling disc"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The first-cut-point theorem is proved in the nilpotent approximation; the paper assumes, without a rigorous error estimate, that this approximation faithfully represents where geodesics of the original rolling disc stop being optimal.","fun_headline_variants_meta":{"raw":{"variants":["Symmetry fixes first cut point of rolling-disc minimizers","Isotropy symmetry pinpoints rolling-disc cut point","Geometric control: symmetry locates first cut point","Rolling disc: symmetry reveals exact cut point","Symmetry algebra locates first cut point in rolling disc"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000193,"raw_usage":{"total_tokens":1302,"prompt_tokens":850,"completion_tokens":452,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":466,"completion_tokens_details":{"reasoning_tokens":375}},"tokens_in":466,"tokens_out":452,"duration_ms":4677,"temperature":1.0,"reasoning_tokens":375,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:15:56.370583+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For the nilpotent system with $C_1=2$ and $C_2^2+C_3^2=1$ (as in Example 1), the proposition predicts the first cut point $(0,\\pi/4,0)$ at time $\\pi$; solve the original rolling-disc extremal system numerically with the same initial covector and test whether any admissible curve shorter than the extremal connects the origin to that point before time $\\pi$. If such a shorter curve exists, or if the first intersection of the true extremal with $x=\\theta=0$ occurs at a value of $y$ different from the nilpotent prediction, the approximation does not control the original cut locus.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the maximum principle and the Hamiltonian system used to derive the extremal equations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the sub-Riemannian framework, left-invariant Hamiltonians, and the definition of cut points used in Section 4."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the algebraic construction of homogeneous nilpotent approximation used to build the Heisenberg model."},{"cited_title":"Hermes , Nilpotent approximations of control systems and distributions","cited_arxiv_id":null,"evidence_quote":"Supplies the nilpotent approximation method for control distributions expressed in privileged coordinates."},{"cited_title":"ˇCap and J","cited_arxiv_id":null,"evidence_quote":"Defines Lagrangian contact structures and the flat model whose symmetry bound Proposition 5.2 matches."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the explicit Jacobi-elliptic solutions of the original extremal system that motivate the nilpotent approximation."},{"cited_title":"Tanaka, On differential systems, graded Lie algebras and pseudo-groups, J","cited_arxiv_id":null,"evidence_quote":"Introduces the algebraic prolongation of graded Lie algebras used in the appendix to compute the symmetry algebra."},{"cited_title":"Yamaguchi, Differential systems associated with simple graded Lie algebras , in Progress in Differential Geometry, Adv","cited_arxiv_id":null,"evidence_quote":"Supplies the classification of differential systems associated with simple graded Lie algebras, used in identifying the prolonged algebra as $\\mathfrak{sl}(3,\\mathbb{R})$."},{"cited_title":"Zelenko, On Tanakas Prolongation Procedure for Filtered Structures of Constant Type, Symmetry, Integrability and Geometry: Methods and Applications SIGMA 5, 094, pp","cited_arxiv_id":null,"evidence_quote":"Provides the prolongation procedure for filtered structures of constant type that yields the full geometric symmetry algebra."}],"review_version":1}