Each of N cortical neural-mass oscillators has its own natural alpha-band frequency ωi, drawn from a Lorentzian spread of width σ around a mean f₀ = 10 Hz. They interact through local (Kuramoto) coupling Kc and are all driven by the same external tACS electric field oscillating at fs with strength Ks:
dθᵢ/dt = ωᵢ + (Kc/N) Σⱼ sin(θⱼ − θᵢ) + Ks · sin(θs(t) − θᵢ)
θs(t) = 2π fs t (external stimulation phase)
R e^(iψ) = (1/N) Σᵢ e^(iθᵢ) (Kuramoto order parameter)
R (0–1) is the phase-locking value: R ≈ 0 means the population fires with random relative phase, R → 1 means the whole cortical patch has synchronized. Two distinct mechanisms raise R here:
- Endogenous synchrony — cortical coupling Kc alone can synchronize the population once it exceeds a critical value (the classic Kuramoto transition), independent of any stimulation.
- Entrainment — the external field pulls the whole ensemble to lock onto fs itself (θᵢ tracks θs). This only succeeds within an "Arnold tongue": entrainment requires Ks large enough relative to |fs − f₀| and to the dispersion σ — a stronger, more detuned stimulus needs more current to entrain the same tissue.
This is the accepted mechanistic model for non-invasive tACS neuromodulation used in sleep, memory and depression research: a weak (1–2 mA) sinusoidal scalp current cannot directly drive neurons to fire, but it can nudge the phase of an already-oscillating population, and if the drive is close enough to the tissue's natural rhythm and coupling is high enough, the whole patch locks to the stimulation frequency.