Paper I
Operational Regime Classification
From Thermodynamic Axioms to Sleep Architecture
Introduction
How do we classify the operational state of a brain? Existing approaches rely on spectral power bands or heuristic staging rules. We show that three thermodynamic axioms—bounded phase space, no null state, and finite observational resolution—are sufficient to derive a complete classification of neural regimes.
Bounded Phase Space
The neural state space has finite measure: μ(Ω) < ∞. No infinite-energy configurations are accessible.
No Null State
The system is never in a state of zero activity. Even under deep anesthesia, residual oscillations persist.
Finite Resolution
Observational precision is bounded: δ > 0. This imposes a natural coarse-graining on the partition.
These axioms uniquely determine the Kuramoto model as the appropriate mathematical backbone for regime classification. The order parameter serves as the single sufficient statistic that partitions neural dynamics into five operationally distinct regimes, each with characteristic synchronization profiles, spectral signatures, and functional correlates.
The Kuramoto Model
A population of coupled oscillators with natural frequencies drawn from a distribution with spread .
Mean-Field Dynamics
Each oscillator evolves according to the mean-field coupling:
where is the natural frequency of oscillator and is the global coupling strength.
Order Parameter
The degree of phase coherence is captured by the Kuramoto order parameter:
indicates complete incoherence, indicates perfect synchrony. All intermediate values define the regime boundaries.
Critical Coupling
The system undergoes a phase transition at the critical coupling strength:
Below , the oscillators remain incoherent. Above , a macroscopic fraction synchronizes spontaneously. This bifurcation separates pathological (turbulent) from functional (cascade/coherent) operation.
The Five Operational Regimes
Every neural state maps to exactly one regime, classified by boundaries and structural factor values.
Complete desynchronization. No macroscopic order emerges. Individual oscillators run at their natural frequencies with no mutual entrainment. The system dissipates maximal energy with no coherent output.
Seizure prodrome, severe delirium, deep anesthesia washout
Fully developed turbulence, uncorrelated spin glass
Selective gating begins. Small clusters of oscillators lock transiently, creating apertures through which specific frequency bands pass. This regime enables flexible filtering without global coherence.
REM sleep, psychedelic states, creative divergent thinking
Partial synchrony onset, cluster formation in Josephson arrays
Cooperative synchronization propagates through the network. Macroscopic clusters form and dissolve on intermediate timescales. The system operates near criticality, maximizing dynamic range and information transmission.
N1/N2 sleep, relaxed wakefulness, default mode network
Cooperative synchronization, critical cascades, power-law avalanches
Healthy baseline operation. The majority of oscillators are entrained to a common frequency. Stable macroscopic rhythms support reliable information processing and motor coordination.
Alert wakefulness, focused attention, N3 deep sleep
Laser above threshold, superfluid phase, Bose-Einstein condensate
Hypersynchrony. Nearly all oscillators are phase-locked. The system loses flexibility and cannot modulate its output. Pathological in neural context: epileptic seizures represent this regime.
Tonic-clonic seizure, catatonia, status epilepticus
Rigid body rotation, ferromagnetic saturation
Sleep Architecture
Sleep stages map systematically onto the regime classification. The 90-minute ultradian cycle traces a periodic orbit through regime space.
Stage-to-Regime Mapping
| Sleep Stage | Regime | R Range | Dominant Band |
|---|---|---|---|
| N3 (Deep Sleep) | Phase-Locked | R > 0.9 | Delta (0.5–4 Hz) |
| W (Wakefulness) | coherent | 0.8 < R < 0.9 | Alpha/Beta (8–30 Hz) |
| N1 (Light Sleep) | cascade | 0.6 < R < 0.8 | Theta (4–8 Hz) |
| N2 (Spindle Sleep) | cascade | 0.5 < R < 0.7 | Sigma (12–15 Hz) |
| REM | turbulent | R < 0.5 | Theta (4–8 Hz) |
R Ordering
The order parameter preserves a strict ordering across all epochs:
This ordering is validated across all 500 epochs in the test dataset. N3 achieves near-unity synchronization (delta dominance), while REM shows maximal desynchronization.
Ultradian Cycling
The ~90-minute sleep cycle is a periodic orbit in regime space:
The trajectory descends from coherent (W) through cascade (N1/N2) to locked (N3), then jumps to turbulent/aperture (REM) before returning. Band power decomposition confirms delta dominance in N3 and theta dominance in REM.
Critical Coupling and Bifurcation
The transition from incoherence to synchrony occurs at a sharp bifurcation point determined analytically.
Bifurcation Point
For a Lorentzian frequency distribution with half-width , the critical coupling is:
Below , the only stable state is (incoherent). Above , a nonzero branch appears via supercritical pitchfork bifurcation.
Finite-Size Scaling
For finite , the transition is smoothed. The order parameter scales as:
The analytical prediction is validated against numerical simulation across system sizes to .
The Consciousness Window
Consciousness requires temporal overlap between perception decay and thought decay.
The consciousness window is defined as the temporal intersection of perceptual decay (the fading of sensory input) and thought decay (the fading of internal representation). When these two processes overlap in time, conscious experience arises. Outside this window, the system operates in either purely sensory (reflexive) or purely internal (unconscious processing) modes.
This window is maximized in the coherent regime (alert wakefulness) and collapses in both the turbulent regime (where perceptual integration fails) and the phase-locked regime (where internal dynamics are frozen).
Variance Minimization
The free energy principle yields exact predictions for variance scaling.
Free Energy
The free energy of the partition is proportional to the variance of the order parameter fluctuations. The system minimizes by reducing , driving toward increased synchronization.
Variance Floor
The variance cannot be reduced below the thermal floor set by the ratio of temperature to coupling. The scaling slope is exactly on a log-log plot of vs .
Figures
Six panels summarizing the key results of the regime classification framework.
Kuramoto order parameter R vs coupling strength K, showing the supercritical bifurcation at K_c
Five-regime phase diagram with R boundaries and representative time series for each regime
Sleep hypnogram overlaid with R trajectory, showing ultradian cycling through regime space
Band power decomposition across sleep stages: delta, theta, alpha, sigma, beta
Variance scaling: log-log plot of sigma^2 vs K with slope = -1.0
Consciousness window C = P_decay intersect T_decay as a function of regime