flexible-phase-locking-cortical-theta

# Flexible Phase-Locking in Cortical Theta Oscillators > Multi-timescale inhibitory current interactions generate flexible phase-locking via delayed Hopf bifurcation, enabling cortical oscillators to entrain to rhythms substantially slower than their intrinsic frequency. ## Metadata - **Source**: arXiv:2605.08014 - **Authors**: Yangyang Wang, Benjamin R. Pittman-Polletta - **Published**: 2026-05-08 - **Categories**: q-bio.NC, math.DS ## Core Methodology ### Key Innovation Cortical oscillators can flexibly phase-lock to inputs spanning a wide range of timescales (including rhythms substantially slower than intrinsic frequency) through a **delayed Hopf bifurcation (DHB)** mechanism driven by multi-timescale inhibitory current interactions. This expands the entrainment frequency range far beyond what single-timescale models can achieve. ### Technical Framework **Three-Timescale Structure:** 1. **Fast timescale** — Intrinsic theta oscillation (4-8 Hz) 2. **Intermediate timescale** — Theta-timescale inhibitory current I_m, expands phase-locking range by prolonging delayed recovery along superslow manifold 3. **Superslow timescale** — Delta-timescale (1-4 Hz) inhibitory potassium current I_{K_{SS}}, critical for entrainment flexibility under external forcing **Delayed Hopf Bifurcation Mechanism:** - Slow and superslow inhibitory processes interact to generate **prolonged post-input recovery delays** - The DHB creates a memory-like effect where the system "remembers" previous inputs, allowing it to track slower rhythms - I_{K_{SS}} plays minimal role in unforced oscillatory dynamics but is recruited specifically for phase-locking under external forcing - I_m is not redundant but cooperatively expands the entrainment range **Dynamical Systems Analysis:** - Use geometric singular perturbation theory to analyze the three-timescale system - Identify the critical manifold structure (fast, slow, superslow manifolds) - Analyze the delayed Hopf bifurcation phenomenon: trajectory stays near repelling slow manifold after bifurcation point - Compute the delay duration and its dependence on timescale separation parameters ### Mathematical Framework ```python # Conceptual model structure # Three-timescale ODE system: # dx/dt = f(x, y, z) # Fast variables (membrane potential) # dy/dt = ε * g(x, y, z) # Slow variables (I_m dynamics) # dz/dt = ε² * h(x, y, z) # Superslow variables (I_{K_{SS}} dynamics) # where 0 < ε << 1 represents timescale separation ``` ## Applications - **Speech processing** — Explain how auditory cortex tracks speech rhythms at variable rates - **Neural coding** — Understand how oscillatory neurons encode temporal information across scales - **Brain-computer interfaces** — Design oscillatory decoders that adapt to variable input rates - **Computational psychiatry** — Model disruptions in temporal processing (e.g., dyslexia, aphasia) ## Pitfalls - DHB effects are sensitive to noise level — biological noise may disrupt the delayed recovery - Three-timescale separation requires careful parameter tuning; insufficient separation eliminates the effect - The mechanism is specific to forced oscillators — spontaneous dynamics may not exhibit the same properties - Validation requires biophysically grounded models; simplified phase oscillator models may miss key dynamics ## Related Skills - kuramoto-brain-network - neural-dynamics-universal-translator - attractor-metadynamics-neural - neuromodulation-rhythmic-pattern-control