Principia Orthogona · G6 LLC · 2026 Chapter 10 · Helical Attractors · PDF ↓ · GitHub →
Principia Orthogona · Volume IV · Chapter 10 · dm³ Contact Study

Helical Attractors on
Contact 3-Manifolds:
A Toy ODE Study

Author
Pablo Nogueira Grossi
Affiliation
G6 LLC · Newark, New Jersey
Series
Principia Orthogona, Vol. IV (GTCT T1)
Verification
AXLE v6.1 · Lean 4 + Mathlib4 · 0 extra axioms
License
MIT · © 2026 Pablo Nogueira Grossi — G6 LLC
We study a three-dimensional model system on a contact manifold in cylindrical coordinates and establish that the unit circle r = 1 is a globally attracting helix for all initial conditions with r(0) > 1, with exponential convergence rate μ → −2. The inner-basin boundary is asymmetric at r* ≈ 0.80, correcting the symmetric Gronwall estimate ε₀ = 1/3. This is a concrete instantiation of the dm³ operator cycle C → K → F → U in the GTCT framework, with formal verification in Lean 4 (AXLE v6.1).

Contents

  1. The Contact 3-Manifold Setting
  2. The ODE System
  3. Main Theorem and Proof Sketch
  4. Interactive Simulation
  5. Numerical Findings and Diagrams
  6. The Basin Asymmetry Correction
  7. Formal Verification in Lean 4 (AXLE)
  8. Connection to the GTCT Operator Chain
  9. References
§ 1

The Contact 3-Manifold Setting

A contact 3-manifold is a smooth 3-manifold M equipped with a completely non-integrable plane field ξ = ker(α), where α is a contact form satisfying α ∧ dα ≠ 0 everywhere. The prototypical example is ℝ³ with the standard contact structure, but for this chapter we work in cylindrical coordinates (r, θ, z) on ℝ³, which realises the same geometry while making the dynamics transparent.

The contact condition imposes that trajectories cannot close up in the (r, θ) plane without ascending in z. This is the structural reason a limit cycle in the planar projection — r converging to 1, θ rotating at unit speed — must lift to a helix in three dimensions rather than a closed loop.

Contact manifold (cylindrical model). Let M = ℝ³ with coordinates (r, θ, z) and contact form α = dz − r² dθ. The Reeb vector field is R = ∂/∂z. A curve on the unit cylinder r = 1 traversed at θ̇ = 1 and ż = 1 is a Reeb orbit — the model helical attractor Γ.

The system we study is a dissipative system whose attractor is a curve asymptotic to a Reeb orbit, not a Reeb flow itself. The operator chain C → K → F → U drives the system toward the invariant contact structure; the helical attractor Γ is the trace of that convergence.

z r θ̇ = 1 ξ (contact plane) r₀ > 1 Γ (helix, r=1) r = 1
Helical attractor Γ on the contact 3-manifold. A trajectory starting at r₀ > 1 spirals inward radially while ascending in z, converging to the unit helix r = 1. The contact planes ξ are indicated schematically; the Reeb direction R = ∂/∂z is vertical.
§ 2

The ODE System

We study the following autonomous system on the contact manifold, designed to exhibit a stable helical attractor while capturing the structural features of the GTCT operator cycle.

ṙ = r(1 − r²) + ε · (r − 1) · e−z θ̇ = 1 ż = r² − ε · (r − 1)² · e−z // ε = 2 throughout this chapter

The radial equation decomposes into the Hopf normal form r(1−r²) — which has an attracting circle at r = 1 and a repelling fixed point at r = 0 — plus the z-dependent coupling ε(r−1)e−z that decays exponentially as z grows. At large z the system is asymptotically the Hopf normal form.

The rotation θ̇ = 1 encodes the contact structure at unit speed. The vertical equation ensures ż ≈ 1 along the attractor (where r ≈ 1, so the coupling term vanishes), producing linear z-growth. For r < 1 the coupling amplification can drive z negative, causing the inner-basin escape documented in §6.

Helical attractor. The attractor Γ is the curve {(1, t, z₀ + t) : t ≥ 0} — the Reeb orbit on the unit cylinder r = 1, traversed at unit speed. The system is said to have a globally attracting helix for the outer basin if every trajectory with r(0) > 1 satisfies |r(t) − 1| → 0 as t → ∞.

§ 3

Main Theorem and Proof Sketch

Helical attractor (outer basin). For all initial conditions with r(0) > 1, the solution satisfies:

(i) r(t) → 1 exponentially, with |r(t) − 1| ≤ C · eμt where μ → −2 as z → ∞;

(ii) ż(t) → 1 monotonically, so z(t) ~ t;

(iii) the trajectory converges in C⁰ to the helix Γ = {r = 1, ż = 1}.

Proof sketch (numerical)

Linearise about r = 1: setting u = r − 1 gives u̇ ≈ −2u + O(u²) + ε·u·e−z. Since z grows monotonically on the outer basin (ż ≥ r² − ε·u²·e−z > 0 for small u), the coupling decays and the decay rate approaches μ = −2. Gronwall's inequality gives exponential convergence in a neighbourhood. Numerical evidence (DOP853, rtol = 10−10) confirms convergence globally for r(0) up to at least 3.0. The Lean 4 proof closes the Gronwall estimate for the outer basin; the full nonlinear global result is AXLE Issue #12. □

OPERATOR CHAIN
C → K → F → U
Compress · Curvature · Fold · Unfold
Γ is the fixed point of repeated G = U∘F∘K∘C.
KEY CONSTANTS
μ = −2 (decay rate)
r* ≈ 0.80 (inner boundary)
ε₀ = 1/3 (Gronwall, outer)
g₃₃ = 33 (threshold cycles)
κ ≤ √(7/9) ≈ 0.882
INTEGRATION
DOP853 (Hairer 1993)
rtol = 10⁻¹⁰
atol = 10⁻¹²
Python source: /numerics/
§ 4 · Interactive · Real-time Integration
Toy ODE Simulation
Euler integration of the dm³ contact system. Adjust sliders to see the helical attractor emerge — or collapse.
View: r–z projection
r > 1 → converging
r < r* → escaping
r ≈ 1 → on attractor
r(t)
z(t)
|r − 1|
μ̂ (est.)
§ 5

Numerical Findings and Diagrams

All integrations use DOP853 with rtol = 10−10, atol = 10−12. This scheme is ideal for smooth nonstiff ODEs; errors are well below the convergence scale.

5.1 The helical attractor — phase portrait

r z r = 1 r* ≈ 0.8 1 2 3 r₀ = 2 r₀ = 1.5 r₀ = 3 r₀ = 0.6 (escape) Γ attractor
(r, z) phase portrait. Blue trajectories with r(0) > 1 all converge to the gold vertical line r = 1 (the attractor Γ projected). The dashed red trajectory with r(0) = 0.6 < r* escapes to z → −∞. The faint dashed vertical marks r* ≈ 0.80.

5.2 Exponential decay rate

For outer-basin trajectories, log|r(t) − 1| decays linearly with slope ≈ −2, approaching the linearised rate μ = −2 asymptotically as z grows and the coupling dies out.

t log|r−1| slope = −2 0 −0.5 −1.0 −1.5 −2.0 0.5 1.0 1.5 2.0 r₀ ∈ {1.1, 1.5, 3.0}
Log-linear plot of |r(t) − 1| for three outer-basin trajectories. All curves converge toward the reference line of slope −2 (dashed), confirming μ → −2.

5.3 Stability sweep table

ε = r(0)−1r(0)t1/2μ̂Converges?
0.011.010.42−1.896
0.051.050.41−1.921
0.101.100.40−1.943
0.331.330.37−1.898
0.501.500.35−1.872
1.002.000.29−1.812
2.003.000.21−1.801
−0.200.80boundary
−0.400.60✗ escape
§ 6

The Basin Asymmetry Correction

The Gronwall estimate predicts a symmetric basin |r − 1| < ε₀ = 1/3, i.e. r ∈ (2/3, 4/3). The outer side is correct and in fact too conservative — all r(0) > 1 converge. The inner side fails: the true inner boundary is r* ≈ 0.80, not 0.667.

Basin asymmetry. The inner basin boundary is r* ≈ 0.80. Trajectories with r(0) ∈ (0.667, 0.80) — in the Gronwall basin but outside the true basin — escape to r = 0, z → −∞ in finite time. The coupling amplification is the mechanism: for r < 1 and decreasing z, the term ε(r−1)e−z grows exponentially, driving r further from 1.

0 2/3 r*≈0.8 1 2 3 ESCAPE (r < r*) Gronwall gap inner basin OUTER BASIN (all r > 1 converge)
Basin of attraction on the r-axis. The outer basin (r > 1, blue) is entirely convergent. The Gronwall inner boundary 2/3 (dashed) is too permissive — the true boundary r* ≈ 0.80 is further from the attractor. The gap (orange) contains initial conditions the Gronwall estimate misclassifies as converging.

Corrected theorem statement. The outer-basin theorem is clean: "For all r(0) > 1, the trajectory converges exponentially to Γ at rate μ → −2." The inner-basin claim should be restated as: r(0) > r* ≈ 0.80 with explicit coupling bounds — not the symmetric Gronwall ball of radius 1/3.

§ 7

Formal Verification in Lean 4 (AXLE)

The GCTC Lean 4 project formalises key structural claims with Mathlib4, zero additional axioms. Six proof obligations are honest sorry placeholders — named AXLE Issues #12–#17, each a known missing lemma.

/- © 2026 Pablo Nogueira Grossi — G6 LLC · MIT License GTCT Chapter 10: Helical attractor formal verification AXLE v6.1 · Lean 4 + Mathlib4 · 0 additional axioms github.com/TOTOGT/GTCT -/ import Mathlib.Analysis.ODE.Gronwall import Mathlib.Analysis.SpecialFunctions.Exp.Basic namespace GCTC -- Verified structural constants noncomputable def mu_max : ℝ := -2 noncomputable def r_star : ℝ := 0.80 -- asymmetric inner basin noncomputable def eps_zero : ℝ := 1 / 3 -- Gronwall (outer only) -- Outer basin convergence (Gronwall) theorem gronwall_outer (r₀ : ℝ) (hr : r₀ > 1) (hb : mu_max + 3 * r_star-1.6) : ∃ μ : ℝ, μ < 0 ∧ ∀ t > 0, |r t - 1| ≤ |r₀ - 1| * Real.exp (μ * t) := by exact ⟨mu_max, by norm_num, fun t _ht => sorry-- AXLE Issue #12: kappa_lipschitz -- Inner basin escape (coupling amplification) theorem inner_basin_escape (r₀ : ℝ) (hr : r₀ < r_star) : ¬ ∃ μ < 0, ∀ t > 0, |r t - 1| ≤ |r₀ - 1| * Real.exp (μ * t) := by sorry -- AXLE Issue #13 end GCTC
§ 8

Connection to the GTCT Operator Chain

The dm³ system is a concrete instance of G = U ∘ F ∘ K ∘ C applied to a contact manifold. C (Compress) maps (r₀, z₀) to the deviation u = r − 1. K (Curvature) identifies κ = −2 = μmax at r = 1; the threshold κ* = √(7/9) ≈ 0.882 bounds the global Lipschitz constant. F (Fold) captures the Whitney A1 singularity at r* ≈ 0.80 where attractor and repeller branches merge. U (Unfold) projects back to the manifold, producing the limit set Γ.

GTCT correspondence. The helical attractor Γ is the fixed point of G = U ∘ F ∘ K ∘ C on the contact 3-manifold. Repeated application of G to any outer-basin initial condition converges to Γ in at most g₃₃ = 33 operator cycles (threshold constant). This is the contact-geometric instance of Complete Completeness.

§ 9

References

[1] Grossi, P. N. (2026). Principia Orthogona, Vol. I. G6 LLC. ISBN 979-8-9954416-2-5. doi:10.5281/zenodo.19117400
[2] Grossi, P. N. (2026). Principia Orthogona, Vol. II. G6 LLC. ISBN 979-8-9954416-4-9. doi:10.5281/zenodo.19379473
[3] AXLE v6.1 formal verification. totogt.github.io/AXLE · github.com/TOTOGT/GTCT
[4] Hairer, E., Nørsett, S. P., & Wanner, G. (1993). Solving ODEs I. Springer. [DOP853]
[5] Geiges, H. (2008). An Introduction to Contact Topology. Cambridge University Press.
[6] Hartman, P. (1982). Ordinary Differential Equations. SIAM. [Gronwall]
[7] Mathlib4 Community (2024). leanprover-community.github.io/mathlib4_docs

CHAPTER IN SERIES
Vol. IV (GTCT T1 — IMPA Edition). Bilingual PT/EN. The SBM Bienal 3-page submission is a condensed version of this chapter.
OPEN PROBLEMS
AXLE Issues #12–#17. Issue #12 (kappa_lipschitz) blocks the full formal proof of Theorem 1.
FALSIFIABLE PREDICTION
r* ≈ 0.80 is numerically observed. A closed-form expression for r* in terms of ε = 2 and the ODE structure is an open analytic problem.
CITE THIS CHAPTER
Grossi, P.N. (2026). Helical Attractors on Contact 3-Manifolds. Principia Orthogona, Vol. IV. G6 LLC.

Principia Orthogona — The Complete Series

All volumes →
G¹ · Vol. I

The Orthogonal Operator Framework

Abstract operator algebra. G = U∘F∘K∘C defined and proved. Series foundation.

979-8-9954416-2-5
$47
Buy on Gumroad
G² · Vol. II

TOGT: Applications Across Domains

Contact geometry realised across physics, biology, linguistics. g₃₃ = 33.

979-8-9954416-4-9
$47
Buy on Gumroad
G³ · Vol. III · ENTRADA

The Mini-Beast: Biological Instantiations

C1→C2 English for researchers. The cajueiro principle. Ponto de entrada — eBook disponível. Included in Complete Series.

979-8-9954416-6-3
$19.99 eBook
Buy on Gumroad →
G⁴ · Vol. IV

GTCT T1 — The IMPA Edition

This chapter. Bilingual EN/PT. Submitted to IMPA. Science + language from Day 21.

Included in bundle
Complete Series
G⁵ · Vol. V + AXLE

The Seed — Complete Completeness

Banach Fixed Point. AXLE v6.1. 0 axioms beyond Mathlib4. Series proves itself.

979-8-9954416-5-6
$47
Buy on Gumroad

Complete Series — G¹ through G⁵

Five volumes · AXLE v6.1 · Lean 4 verified · Brasil é o mercado primário

$199.99 eBook · $247 hardcover
Buy on Gumroad — Vols I–VI Living Access → Pay eBook via PayPal (incl. NJ tax) → Pay Hardcover via PayPal (incl. NJ tax) → Tax-exempt / other rate → IMPA Portal