1. Overview & Biophysical Foundations
The Circular Motion Illusion (fundamentally rooted in the geometry of the Tusi Couple and the kinematics of Cardan Gears) demonstrates how an ensemble of purely one-dimensional rectilinear simple harmonic oscillations (SHM) synthesizes the compelling optical perception of a coherent, rigid two-dimensional rotating ring. When reduced to a single particle ($N=1$), the motion is unambiguously recognized as a straight 1D back-and-forth oscillation along a single diameter. When two orthogonal particles ($N=2$) oscillate $90^\circ$ out of phase, they trace the basic Cartesian components ($X = R \cos \theta$, $Y = R \sin \theta$). As the number of particles increases ($N \ge 8$) with equidistant phase distribution across radial spokes, the human visual system ceases to perceive discrete linear reciprocating tracks and instead perceives a continuous, rigid 2D ring revolving smoothly in a circle.
In 1247 AD, the Persian astronomer Nasir al-Din al-Tusi formulated this mathematical couple in his treatise Tadhkira fi ilm al-hay'a to resolve linear and circular planetary motions. Nicolaus Copernicus adopted the Tusi couple in De revolutionibus orbium coelestium (1543) to describe reciprocal planetary deviations. In 1570, Gerolamo Cardano demonstrated how planetary hypocycloidal gears with a 2:1 radius ratio convert pure rotary motion into straight-line reciprocating mechanical stroke without crank-slider linkages.
From a sensory psychophysics standpoint, the illusion highlights fundamental mechanisms of cortical visual processing. Neurons in the primary visual cortex (area V1) act as localized spatiotemporal motion energy filters, sensitive only to local component vectors (the classic "aperture problem"). When visual signals propagate along the dorsal visual stream to the middle temporal area (area MT / V5), higher-order visual circuits perform global spatial pooling and vector integration. Guided by Gestalt principles of common fate, grouping, and continuous contour closure, the cortical visual hierarchy favors the simpler, unified perceptual hypothesis of a single rigid rotating structure over the overwhelming cognitive complexity of tracking multiple independent 1D linear oscillators.
This principle of orthogonal and multiphasic harmonic synthesis is ubiquitous across biological systems: cardiac vectorcardiography maps myocardial electrical depolarization wavefronts as orthogonal dipole potentials, spinal central pattern generator (CPG) neural circuits coordinate rhythmic multi-limb locomotion through phase-coupled linear oscillators, and respiratory sinus arrhythmia couples cardiovascular and pulmonary rhythms across orthogonal harmonic phases.
2. Interactive Laboratory Operation Guide
This simulation environment provides real-time parametric manipulation of geometric harmonic oscillators, phase-space scrubbing, guided educational demonstration, and acoustic sonification:
- Baseline State (Pure Illusion): By default, all guide lines, rolling circles, and helpers are hidden. Observe the glowing ensemble of particles moving in a seamless, fluid circular formation against the dark oscilloscope grid.
- Starting from a Single Dot ($N = 1$): Drag the Number of Dots slider down to $1$. Observe that a single particle simply moves back and forth along a single horizontal diameter line. Increase to $N = 2$ to see two orthogonal oscillators, and then to $N = 4, 8, 16$ to watch the circular illusion emerge before your eyes.
- Guided Demonstration Tour (Start Demo): Clicking the
Start Demo button launches an automated 5-stage progressive tour:
- Stage 1 (Single Dot): Demonstrates pure 1D rectilinear harmonic motion on a single spoke.
- Stage 2 (Orthogonal Pair, $N=2$): Demonstrates quadrature sine and cosine linear decomposition.
- Stage 3 (Emergence, $N=4 \rightarrow 8$): Smoothly increases spoke density, showing the emergence of symmetry.
- Stage 4 (Pure Illusion, $N=16$): Hides all lines and guides, presenting the pure 2D revolving ring illusion.
- Stage 5 (Geometric Revelation): Activates the rolling Tusi hypocycloid circle and velocity vectors to explain the underlying planetary gear geometry.
Note: Interacting with any slider, button, or canvas element immediately interrupts the demo and restores your original baseline configuration.
- Revealing the Real Paths (Show Linear Spokes): In the Visual Guides & Diagnostics panel, toggle
Show Linear Spokes to reveal the fixed straight radial diameters. Fixate on any individual particle to verify that it never departs from its designated straight track.
- Revealing the Rolling Geometry (Show Rolling Tusi Guide): Enable
Show Rolling Tusi Guide to render the virtual hypocycloid circle of radius $r = R/2$ rolling inside the boundary. Every particle sits permanently on the circumference of this rolling disk.
- Direct Interactive Scrubbing: Click or touch and drag anywhere across the visualizer canvas. The engine calculates the polar angle $\theta = \operatorname{atan2}(y - c_y, x - c_x)$, freezing the simulation and allowing manual phase scrubbing.
- Harmonic Audio Sonification (Sound ON/OFF): Toggle the Web Audio synthesizer. As particles cross the central equilibrium point ($d = 0$) at peak velocity, the engine triggers pentatonic bell tones, translating the spatial phase distribution into acoustic polyrhythms.
- Reset Baseline: The
Reset Baseline button immediately restores the system to its initial pure illusion state with default velocity, particle count, and hidden guides.
3. Mathematical Derivation & Engineering Architecture
Let $R$ represent the outer boundary radius and $N$ represent the total number of linear spoke tracks ($N \ge 1$). The orientation angle of the $k$-th spoke ($k \in \{0, 1, \dots, N-1\}$) is defined by:
$$\phi_k = \frac{k \pi}{N}$$
Each particle executes simple harmonic motion along its spoke. Its signed displacement $d_k(t)$ from the coordinate origin $(c_x, c_y)$ as a function of phase $\theta(t) = \omega t$ is:
$$d_k(t) = R \cos(\theta(t) - \phi_k)$$
The Cartesian position vector $\mathbf{p}_k(t) = \begin{bmatrix} x_k(t) \\ y_k(t) \end{bmatrix}$ is given by projecting $d_k(t)$ along the unit direction vector $\begin{bmatrix} \cos(\phi_k) \\ \sin(\phi_k) \end{bmatrix}$:
$$\mathbf{p}_k(t) = \begin{bmatrix} c_x + R \cos(\theta(t) - \phi_k) \cos(\phi_k) \\ c_y + R \cos(\theta(t) - \phi_k) \sin(\phi_k) \end{bmatrix}$$
Applying the standard trigonometric product-to-sum identities:
$$\cos(A)\cos(B) = \frac{1}{2}[\cos(A+B) + \cos(A-B)]$$
$$\cos(A)\sin(B) = \frac{1}{2}[\sin(A+B) - \sin(A-B)]$$
With $A = \theta(t) - \phi_k$ and $B = \phi_k$, the expansion evaluates to:
$$x_k(t) - c_x = \frac{R}{2}\left[\cos(\theta(t)) + \cos(\theta(t) - 2\phi_k)\right]$$
$$y_k(t) - c_y = \frac{R}{2}\left[\sin(\theta(t)) - \sin(\theta(t) - 2\phi_k)\right] = \frac{R}{2}\left[\sin(\theta(t)) + \sin(-\theta(t) + 2\phi_k)\right]$$
In vector notation, this decomposes into the sum of two rotating vectors:
$$\mathbf{p}_k(t) - \mathbf{c} = \underbrace{\frac{R}{2} \begin{bmatrix} \cos(\theta(t)) \\ \sin(\theta(t)) \end{bmatrix}}_{\mathbf{C}_{\text{roll}}(t)\text{ (Center of Rolling Disk)}} + \underbrace{\frac{R}{2} \begin{bmatrix} \cos(-\theta(t) + 2\phi_k) \\ \sin(-\theta(t) + 2\phi_k) \end{bmatrix}}_{\text{Vector along Rolling Disk Perimeter}}$$
This algebraic identity provides the complete mathematical proof:
- The center $\mathbf{C}_{\text{roll}}(t)$ of an inner disk of radius $r = R/2$ orbits the global origin at radius $R/2$ with positive angular velocity $+\omega$.
- Every particle $\mathbf{p}_k(t)$ lies on the circumference of this inner disk, rotating relative to $\mathbf{C}_{\text{roll}}(t)$ with equal and opposite angular velocity $-\omega$.
- Because the rolling disk radius is exactly half the outer radius ($r = R/2$), the hypocycloidal path of every perimeter point degenerates into a pure straight line segment of length $2R$.
Instantaneous Velocity Vector: Differentiating $\mathbf{p}_k(t)$ with respect to time yields the instantaneous velocity vector $\mathbf{v}_k(t)$:
$$\mathbf{v}_k(t) = \frac{d\mathbf{p}_k}{dt} = -R \omega \sin(\theta(t) - \phi_k) \begin{bmatrix} \cos(\phi_k) \\ \sin(\phi_k) \end{bmatrix}$$
Velocity reaches its maximum magnitude $v_{\max} = R\omega$ precisely as the particle crosses the origin ($d_k = 0$), and drops to zero at the turning boundaries ($d_k = \pm R$), matching the kinetic-potential energy interchange of an ideal mechanical oscillator.
High-DPI Display Scaling & Performance Safeguards: The rendering pipeline operates on an unthrottled requestAnimationFrame loop using sub-pixel antialiasing. The canvas buffer dimensions are dynamically scaled by window.devicePixelRatio without mutating CSS layout styles, preventing layout shifts and infinite resizing loops. Floating-point coordinates are strictly bound to eliminate potential NaN deadlocks.