Cyclically Symmetric Continuous ODE First Described: 1999 Featured Archetype

Thomas' Cyclically Symmetric Attractor

Discovered by René Thomas

Executive Summary

A mesmerizing, endlessly winding labyrinth with threefold cyclic symmetry and sine non-linearities.

System of Equations

Phase Space Formulation
dxdt=sin(y)bx\frac{dx}{dt} = \sin(y) - b x
dydt=sin(z)by\frac{dy}{dt} = \sin(z) - b y
dzdt=sin(x)bz\frac{dz}{dt} = \sin(x) - b z

Governing Parameters

bb
Dissipation Coefficient = 0.208186

Controls the frictional damping of the system. Chaos emerges for b < 0.208186; above this threshold, motion collapses to fixed points.

Overview

Proposed in 1999 by Belgian molecular biologist and cyberneticist René Thomas, this deceptively simple system of differential equations models friction-damped motion through a periodic three-dimensional force field.

Thomas was investigating feedback loops and cyclic networks in gene regulation, searching for the simplest differential equations that could sustain continuous deterministic chaos. The result is a mathematically minimal masterpiece characterized by pure cyclic symmetry: every equation is an exact cyclical permutation of the other two (xyzxx \to y \to z \to x).

The Mathematical Model

The Thomas system uses sinusoidal driving forces balanced against linear dissipation:

dxdt=sin(y)bx\frac{dx}{dt} = \sin(y) - bx

dydt=sin(z)by\frac{dy}{dt} = \sin(z) - by

dzdt=sin(x)bz\frac{dz}{dt} = \sin(x) - bz

Here, bb acts as a damping or friction parameter:

  • b>0.208186b > 0.208186: The system possesses stable equilibrium points; all trajectories eventually slow down and converge to rest.
  • b0.208186b \approx 0.208186: The system undergoes a bifurcation into deterministic chaos.
  • b0b \to 0: The system approaches a conservative (Hamiltonian) flow, generating an infinite, space-filling labyrinth.

Topological Characteristics

  • Infinite Labyrinth: Because the non-linear terms are periodic trigonometric functions (sin\sin), the vector field repeats across three-dimensional space, giving rise to interconnected tubes and chambers through which the particle wanders chaotically.
  • Cyclic Isometry: Rotating any state by 120120^\circ around the diagonal axis x=y=zx = y = z maps the attractor precisely onto itself.
  • Organic Weave: In long-exposure particle renderings, the Thomas attractor resembles a dense bundle of woven, glowing silk ribbons that curve gracefully around cubic cells without ever self-intersecting.
Previous Type
Aizawa Attractor
All Attractor Types
Next Type
Halvorsen Attractor