| name | multi-ensemble-mean-field-oscillators |
| category | neuroscience |
| description | Multi-ensemble mean-field reduction for networks of globally coupled phase oscillators with arbitrary (empirical) frequency distributions. Extends Ott-Antonsen to heterogeneity beyond Lorentzians. arXiv:2607.09516 |
| created | 2026-07-13T00:00:00.000Z |
| source | arXiv:2607.09516v1 (Gast, Takasu, Schmidt, Kennedy, 2026-07-10) |
Multi-Ensemble Mean-Field Reduction for Coupled Phase Oscillators
Overview
The Ott-Antonsen (OA) ansatz gives an exact, drastic dimensionality reduction for
globally coupled Kuramoto-style phase oscillators — but only when the oscillator
natural-frequency distribution belongs to a small closed family (e.g. Lorentzian,
Gaussian, Gamma). Real biological/physical systems have empirical frequency
distributions that violate this closure. This paper introduces a data-driven
multi-ensemble method that recovers OA-level dimensionality reduction for
arbitrary frequency distributions by decomposing the distribution into a mixture
of OA-admissible components.
Why It Matters for Neuroscience
Neural populations (cortical oscillators, central pattern generators, hippocampal
theta/gamma mixtures) are heterogeneous. Being able to apply OA stability,
sensitivity, and bifurcation analysis to empirically measured frequency histograms
— rather than forcing a Lorentzian fit — lets you study real circuits with the
full toolkit of low-dimensional mean-field theory.
Core Methodology
1. Mixture Decomposition of the Frequency Distribution
- Given an empirical (or measured) natural-frequency density
g(ω), fit it as a
convex mixture of K OA-admissible base distributions g_k(ω):
g(ω) ≈ Σ_k π_k g_k(ω), Σ_k π_k = 1, π_k ≥ 0.
- Lorentzians are the canonical choice for
g_k (each yields a closed OA manifold),
but any distribution in the OA-closed family works.
- Use a standard mixture fit (e.g. expectation-maximization, or least-squares on the
empirical histogram / kernel density).
2. Per-Ensemble Ott-Antonsen Reduction
- For each component
k, apply the OA ansatz to obtain its low-dimensional order
parameter z_k(t) (complex mean field) governed by a closed ODE:
dz_k/dt = f(z_k, {z_j}, coupling).
- Each ensemble contributes its own order parameter and its coupling to the global
mean field (weighted by
π_k).
3. Coupled Low-Dimensional System
- The full population is now described by
K complex variables instead of N phases
(or an infinite density). This is the multi-ensemble mean-field system.
- Bifurcation, stability, and sensitivity analyses (Hopf, saddle-node, Lyapunov
exponents) run on the
K-dimensional system — cheap even for N → ∞.
Implementation Steps
- Collect/measure the natural-frequency distribution
g(ω) of the oscillator
population (from data or a specified PDF).
- Fit a mixture of
K Lorentzians (or other OA-closed bases) to g(ω). Choose
K by a model-selection criterion (BIC/AIC) — often small (K = 3–8).
- Write OA equations per ensemble: for Kuramoto with global coupling
K_c and
mean field Z = Σ_k π_k z_k:
dz_k/dt = -iω̄_k z_k - Δ_k/2 (z_k - z_k*) + (K_c/2)(Z - Z* z_k)(1 - z_k²)/...
(use the standard OA Lorentzian reduction for each component).
- Integrate / analyze the coupled
z_k system; recover macroscopic observables
(sync order parameter |Z|, incoherent fraction) and run bifurcation scans over
K_c, Δ_k, π_k.
- Validate against direct
N-oscillator simulation (e.g. N = 10^4) on the
empirical g(ω) to confirm the reduction matches.
Key Innovation
Renders the Ott-Antonsen equations directly applicable to empirical frequency
distributions via a data-driven multi-ensemble decomposition, achieving drastic
dimensionality reduction and enabling stability/sensitivity/bifurcation analysis of
real-world physical and biological oscillator systems previously outside OA's closure.
Activation / Triggers
multi-ensemble mean-field, Ott-Antonsen extension, heterogeneous phase oscillators,
arbitrary frequency distribution, Kuramoto reduction, coupled oscillator stability,
bifurcation analysis heterogeneous populations, empirical frequency distribution
Verification
- Reduction matches direct simulation on the empirical distribution within tolerance.
- Mixture fit captures the empirical histogram (residuals small; BIC justifies K).
- Bifurcation diagram of the K-dim system reproduces the full-population transition.
