slow-rhythms-delay-coupled-oscillators

name	slow-rhythms-delay-coupled-oscillators
description	Systematic bifurcation analysis framework for discovering delay-induced slow rhythms in neural oscillator networks. Phase reduction + numerical continuation reveals Hopf/heteroclinic/saddle-node bifurcations organizing slow-fast dynamics. Applicable to FHN, ML, QIF models. Activation: delay-induced rhythms, slow-fast dynamics, phase reduction, bifurcation analysis, neural oscillators, numerical continuation.
tags	["neural-dynamics","bifurcation-analysis","phase-reduction","delay-coupling","oscillators","slow-rhythms","computational-neuroscience"]
source	arXiv:2606.20733
date	2026-06-17T00:00:00.000Z

Dissecting Emerging Slow Rhythms in Delay-Coupled Neural Oscillators

arXiv:2606.20733 | Published: 2026-06-17 Authors: Xinxin Qie, Matteo Martin, Shenquan Liu, Morten Gram Pedersen Subjects: nlin.CD, math.DS, q-bio.NC, q-bio.QM

Core Discovery

Synaptic transmission delays in inhibitory neural networks create an effective slow-fast structure in phase-difference dynamics, generating low-frequency components that are NOT intrinsic cellular properties but emerge from network coupling. This is a generic phenomenon not specific to any particular model.

Methodology Framework

Phase Reduction with Delays

1. Phase Response Curves (PRCs): Compute PRCs for individual oscillators
2. Phase-Difference Model: Derive delayed phase-difference equations for mutually inhibitory coupled pairs
3. Model Generality: Validated across FitzHugh-Nagumo, Morris-Lecar, and QIF-derived neural mass models

Bifurcation Analysis Pipeline

1. Identify synaptic delay τ as bifurcation parameter
2. Construct phase-plane for phase-difference dynamics
3. Apply numerical continuation (e.g., AUTO, MatCont)
4. Map multistability and limit cycles
5. Identify bifurcations: Hopf, heteroclinic, saddle-node-of-periodics
6. Correlate with slow modulating rhythms in full model

Key Insight

• Limit cycles in phase-reduced model → slow amplitude modulation in full model
• Delay creates effective timescale separation even when oscillators are identical
• Slow rhythms arise from phase-difference dynamics, not amplitude dynamics

Bifurcation Types Observed

1. Hopf Bifurcation: Transition from fixed point to limit cycle (onset of slow modulation)
2. Heteroclinic Bifurcation: Connection between saddle points (slow passage near saddles)
3. Saddle-Node-of-Periodics: Creation/destruction of limit cycles (amplitude modulation onset/offset)

Applications

Neuroscience

• Understanding theta/gamma coupling mechanisms
• Explaining slow oscillations in cortical networks without intrinsic slow currents
• Modeling delay-dependent rhythm generation in thalamocortical circuits

Computational Modeling

• Predicting delay-induced dynamics in large-scale brain models
• Designing neuromorphic circuits with controllable timescales
• Optimizing coupling delays for desired oscillation patterns

Clinical Relevance

• Understanding pathological slow rhythms in epilepsy
• Modeling tremor generation in basal ganglia networks
• Investigating delay alterations in neurodegenerative diseases

Implementation Guide

Phase Reduction Steps

# 1. Compute PRC for isolated oscillator
# Z(θ) = dθ/dI at phase θ for small perturbation I

# 2. Derive phase-difference equation with delay
# dφ/dt = ω + ε * [Z(φ) * f(φ(t-τ)) - Z(-φ) * f(-φ(t-τ))]
# where φ is phase difference, τ is synaptic delay

# 3. Analyze fixed points and limit cycles
# Fixed points: dφ/dt = 0 → phase-locked states
# Limit cycles: slow modulation of phase difference

Numerical Continuation

• Software: AUTO-07p, MatCont, or PyDSTool
• Parameters: Track solutions as τ varies
• Detection: Identify bifurcation points via test functions
• Output: Bifurcation diagrams showing solution branches

Validation Across Models

Model	Type	Delay Effect
FitzHugh-Nagumo	Relaxation oscillator	Slow-fast separation enhanced
Morris-Lecar	Conductance-based	Gamma-theta coupling
QIF Neural Mass	Population model	Macroscopic slow rhythms

All models show generic delay-induced slow rhythms via same bifurcation mechanisms.

Key Equations

Phase-Reduced Dynamics

dφ/dt = ω - ε * H(φ, φ(t-τ))
where H is coupling function derived from PRC

Bifurcation Condition

At Hopf: Re(λ) = 0, Im(λ) ≠ 0
λ = eigenvalue of linearized delayed system

Pitfalls & Considerations

1. Weak Coupling Assumption: Phase reduction valid only for ε << 1
2. Delay Estimation: τ must be known or estimated from data
3. Multiple Delays: Network with heterogeneous delays requires extension
4. Noise Effects: Stochastic perturbations can shift bifurcation points

Connections to Existing Skills

• [[adaptive-conduction-delays-haken-lighthouse]]: Complementary framework for spiking networks with plastic delays
• [[kuramoto-brain-network]]: Phase oscillator framework for brain networks
• [[bipartite-oscillator-synchronization]]: Synchronization in inhibitory-excitatory networks

References

@article{qie2026slow,
  title={Dissecting emerging slow rhythms in delay-coupled neural oscillators},
  author={Qie, Xinxin and Martin, Matteo and Liu, Shenquan and Pedersen, Morten Gram},
  journal={arXiv preprint arXiv:2606.20733},
  year={2026}
}