pith. sign in

arxiv: 2603.09402 · v1 · submitted 2026-03-10 · 🧬 q-bio.NC · cond-mat.dis-nn· nlin.CD

Compact Dynamical Mean-Field Theory of Oscillator Networks

Kanishka Reddy This is my paper

Pith reviewed 2026-05-15 13:41 UTC · model grok-4.3

classification 🧬 q-bio.NC cond-mat.dis-nnnlin.CD
keywords dynamical mean-field theoryphase oscillatorssynchronizationphase response curveKuramoto modelneural mass modelsquenched disorderOtt-Antonsen reduction
0
0 comments X

The pith

A compact dynamical mean-field theory reduces large networks of phase oscillators to a single stochastic equation with self-consistent colored noise.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper constructs a dynamical mean-field theory for networks of coupled phase oscillators on the circle that incorporates both mean-field interactions and quenched disorder. The approach begins with wrapped Langevin dynamics and uses a path-integral formulation to average over the disorder in the thermodynamic limit, resulting in an effective single-oscillator model driven by a deterministic mean field and colored Gaussian noise. The noise covariance is determined self-consistently from the circular two-time correlator of the phases. This formalism recovers classical reductions like the Ott-Antonsen ansatz for vanishing disorder and extends to couplings derived from neuron phase response curves. It enables quantitative predictions of synchronization thresholds directly from single-neuron data for biophysical models such as adaptive exponential integrate-and-fire neurons.

Core claim

In the thermodynamic limit, averaging over the quenched randomness in the couplings reduces the network to a single-oscillator stochastic differential equation driven by a deterministic mean-field term and a self-consistent colored Gaussian noise. The noise's two-time covariance is fixed exactly by the circular correlator that respects the 2π periodicity of the phases, allowing the theory to handle arbitrary periodic coupling functions including those fitted from infinitesimal phase response curves.

What carries the argument

The self-consistent colored Gaussian noise whose covariance is determined by the circular two-time correlator of the oscillator phases.

If this is right

  • The theory exactly recovers the Ott-Antonsen reduction when disorder vanishes.
  • It reproduces standard Kuramoto and theta-neuron neural-mass equations in appropriate limits.
  • Inserting couplings fitted from iPRCs of adaptive exponential integrate-and-fire neurons produces quantitative predictions for synchronization thresholds.
  • The framework provides a direct mapping from single-neuron phase response data to network-level mean-field predictions for any phase-reducible oscillators.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Simulations of finite but large networks could be compared to the DMFT to quantify finite-size effects and the accuracy of the Gaussian noise approximation.
  • This method might be adapted to study other collective phenomena in oscillator networks beyond synchronization, such as chimera states or pattern formation.
  • Experimental measurements of phase response curves in neurons could be plugged into this theory to predict population behavior without simulating every cell.

Load-bearing premise

Disorder averaging in the thermodynamic limit yields precisely a colored Gaussian noise with no surviving higher-order or non-Gaussian corrections.

What would settle it

Simulate a large network of adaptive exponential integrate-and-fire neurons with iPRC-derived couplings and check whether the observed synchronization threshold matches the one predicted by solving the compact DMFT equations.

Figures

Figures reproduced from arXiv: 2603.09402 by Kanishka Reddy.

Figure 1
Figure 1. Figure 1: FIG. 1. iPRC-based parameterization of compact DMFT. [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. DMFT correlator closure and finite-size scaling. (a– [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
read the original abstract

We present a compact dynamical mean-field theory (DMFT) for large networks of coupled phase oscillators whose phases live on the circle $S^1$ and interact with both coherent mean-field coupling and quenched randomness. Starting from wrapped Langevin dynamics, we build a path-integral representation that keeps the $2\pi$-periodicity of the phases explicit. After averaging over the disorder in the thermodynamic limit, this construction reduces to a single-oscillator stochastic equation driven by a deterministic mean field and a self-consistent colored Gaussian noise, whose covariance is fixed by a circular two-time correlator. In the limit of vanishing disorder, the formalism reproduces the Ott--Antonsen reduction and recovers standard Kuramoto and theta-neuron neural-mass equations. The same framework accommodates arbitrary $2\pi$-periodic coupling functions, including those obtained from infinitesimal phase response curves (iPRCs) of biophysical neuron models. As an example, we show that for adaptive exponential integrate-and-fire neurons, inserting an iPRC-fitted coupling into the compact DMFT yields quantitative predictions for synchronization thresholds, providing a direct route from single-neuron phase response data to network-level mean-field predictions for arbitrary phase-reducible oscillators.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper develops a compact dynamical mean-field theory (DMFT) for large networks of coupled phase oscillators on S¹ with both coherent mean-field coupling and quenched randomness. Starting from wrapped Langevin dynamics, a path-integral representation is constructed that preserves 2π-periodicity. After disorder averaging in the thermodynamic limit, the network reduces to a single-oscillator SDE driven by deterministic mean field plus self-consistent colored Gaussian noise whose covariance is fixed by the circular two-time correlator. The formalism recovers the Ott-Antonsen reduction for vanishing disorder and accommodates arbitrary 2π-periodic couplings, including iPRC-derived functions from biophysical models; an aEIF neuron example is used to obtain quantitative synchronization-threshold predictions directly from single-neuron phase-response data.

Significance. If the central reduction is valid, the work supplies a concrete, data-driven route from measured iPRCs to network-level mean-field predictions for synchronization in arbitrary phase-reducible oscillator networks. This is potentially significant for neural modeling, as it links single-cell biophysical measurements to collective dynamics without requiring ad-hoc fitting at the population level.

major comments (3)
  1. [path-integral construction after quenched averaging] Path-integral construction and quenched averaging (the reduction to the effective single-oscillator SDE): the assertion that thermodynamic-limit disorder averaging produces exactly Gaussian colored noise with vanishing connected cumulants of order 3 and higher is load-bearing for the entire closure. For phases on S¹ and arbitrary 2π-periodic (including iPRC-derived) couplings under nonlinear self-consistency, the standard CLT argument does not automatically apply; an explicit demonstration or bound on residual non-Gaussianity and finite-N corrections is required, particularly for the parameter regimes of the aEIF example.
  2. [aEIF example] aEIF example and iPRC fitting: the quantitative synchronization-threshold predictions rely on inserting an iPRC-fitted coupling into the compact DMFT. The manuscript should report the explicit fitting procedure, the resulting coupling function, and direct numerical comparisons (e.g., to finite-N simulations) that quantify the accuracy of the predicted thresholds; without these, the claim of a 'direct route' from single-neuron data remains unverified.
  3. [vanishing disorder limit] Recovery of Ott-Antonsen limit: the statement that vanishing disorder reproduces the Ott-Antonsen reduction and standard Kuramoto/theta-neuron equations is a key consistency check. The explicit steps showing how the circular correlator and self-consistent noise reduce to the OA ansatz should be provided in detail, including any assumptions on the noise spectrum.
minor comments (2)
  1. Notation for the circular two-time correlator and the self-consistent noise covariance should be introduced with a clear equation number and distinguished from the ordinary two-time correlator to avoid ambiguity.
  2. The abstract and introduction would benefit from a brief statement of the precise form of the wrapped Langevin dynamics used as the starting point.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading of the manuscript and the constructive comments. We address each major point below and will incorporate the suggested clarifications and additions in the revised version.

read point-by-point responses

  1. Referee: Path-integral construction and quenched averaging (the reduction to the effective single-oscillator SDE): the assertion that thermodynamic-limit disorder averaging produces exactly Gaussian colored noise with vanishing connected cumulants of order 3 and higher is load-bearing for the entire closure. For phases on S¹ and arbitrary 2π-periodic (including iPRC-derived) couplings under nonlinear self-consistency, the standard CLT argument does not automatically apply; an explicit demonstration or bound on residual non-Gaussianity and finite-N corrections is required, particularly for the parameter regimes of the aEIF example.

    Authors: We agree that the Gaussian character of the effective noise after quenched averaging requires explicit justification for phases on the circle and nonlinear couplings. The manuscript invokes the central limit theorem for the sum of independent disorder terms in the N→∞ limit, but we will add a dedicated appendix subsection deriving the cumulant-generating function for the circular variables. This will demonstrate that connected cumulants of order ≥3 vanish as O(1/N) and provide explicit bounds evaluated at the aEIF parameters and coupling strengths used in the example. Finite-N corrections will be illustrated with a brief numerical check. revision: yes

  2. Referee: aEIF example and iPRC fitting: the quantitative synchronization-threshold predictions rely on inserting an iPRC-fitted coupling into the compact DMFT. The manuscript should report the explicit fitting procedure, the resulting coupling function, and direct numerical comparisons (e.g., to finite-N simulations) that quantify the accuracy of the predicted thresholds; without these, the claim of a 'direct route' from single-neuron data remains unverified.

    Authors: We thank the referee for highlighting this gap. The revised manuscript will expand the aEIF section to include: (i) the precise least-squares fitting procedure applied to the iPRC data, (ii) the explicit form of the resulting 2π-periodic coupling function (given both analytically as a Fourier series and plotted), and (iii) quantitative comparisons of the DMFT-predicted synchronization threshold against direct simulations of finite-N networks (N=500 and N=2000) across a range of parameters, reporting relative errors and standard deviations over multiple realizations. These additions will substantiate the data-driven route. revision: yes

  3. Referee: Recovery of Ott-Antonsen limit: the statement that vanishing disorder reproduces the Ott-Antonsen reduction and standard Kuramoto/theta-neuron equations is a key consistency check. The explicit steps showing how the circular correlator and self-consistent noise reduce to the OA ansatz should be provided in detail, including any assumptions on the noise spectrum.

    Authors: We agree that the consistency check deserves an explicit derivation. In the revision we will insert a new subsection (or short appendix) that walks through the reduction: setting the disorder variance to zero eliminates the colored-noise term, leaving a deterministic drive; the circular two-time correlator then satisfies the closed deterministic equations of the Ott–Antonsen ansatz; and the standard Kuramoto and theta-neuron mean-field equations are recovered for the corresponding choices of coupling. The assumption that the noise spectrum collapses to a deterministic (zero-variance) limit will be stated clearly. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained from explicit dynamics

full rationale

The paper begins from wrapped Langevin dynamics on S¹, constructs an explicit path-integral representation preserving 2π-periodicity, performs quenched averaging in the thermodynamic limit, and obtains a single-oscillator SDE whose colored Gaussian noise covariance is fixed by the two-time circular correlator. This self-consistency is the standard DMFT closure that must be solved for the order parameters; it is not tautological because the microscopic SDE plus averaging produces the effective equation whose solution yields the correlator. The zero-disorder limit recovers the Ott–Antonsen reduction and known Kuramoto/theta-neuron equations as an independent check. Insertion of an iPRC-fitted coupling function uses external single-neuron data to define the interaction term, after which synchronization thresholds are computed from the resulting self-consistent equations rather than being redefined by the fit. No load-bearing self-citations, uniqueness theorems, or ansatzes imported from prior work appear in the provided derivation chain.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The construction rests on the thermodynamic limit for disorder averaging and the validity of a path-integral representation for wrapped Langevin dynamics; no new particles or forces are introduced, and the only fitted elements are the coupling function parameters taken from external iPRC data.

free parameters (1)
  • iPRC coupling parameters
    The 2π-periodic coupling function is obtained by fitting to phase-response data from biophysical neuron models, introducing parameters that are not derived from first principles within the DMFT.
axioms (2)
  • domain assumption Thermodynamic limit allows exact reduction to self-consistent single-oscillator dynamics with colored Gaussian noise
    Invoked when averaging over quenched randomness to obtain the effective stochastic equation.
  • standard math Path-integral representation exactly preserves 2π-periodicity of the phase variables
    Starting assumption when mapping wrapped Langevin dynamics to the functional integral.

pith-pipeline@v0.9.0 · 5510 in / 1445 out tokens · 53140 ms · 2026-05-15T13:41:27.225948+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

47 extracted references · 47 canonical work pages

  1. [1]

    Y. Kuramoto, Self-entrainment of a population of cou- pled nonlinear oscillators, inInternational Symposium on Mathematical Problems in Theoretical Physics, Lecture Notes in Physics, Vol. 39, edited by H. Araki (Springer, Berlin, 1975) pp. 420–422

    work page 1975

  2. [2]

    S. H. Strogatz, From Kuramoto to Crawford: exploring the onset of synchronization in populations of coupled oscillators, Physica D143, 1 (2000)

    work page 2000

  3. [3]

    J. A. Acebr´ on, L. L. Bonilla, C. J. P. Vicente, F. Ritort, and R. Spigler, The Kuramoto model: A simple paradigm for synchronization phenomena, Rev. Mod. Phys.77, 137 (2005)

    work page 2005

  4. [4]

    F. A. Rodrigues, T. K. D. M. Peron, P. Ji, and J. Kurths, The Kuramoto model in complex networks, Phys. Rep. 610, 1 (2016)

    work page 2016

  5. [5]

    Wiesenfeld, P

    K. Wiesenfeld, P. Colet, and S. H. Strogatz, Frequency locking in Josephson arrays: Connection with the Ku- ramoto model, Phys. Rev. E57, 1563 (1998)

    work page 1998

  6. [6]

    Buzs´ aki and X.-J

    G. Buzs´ aki and X.-J. Wang, Mechanisms of gamma os- cillations, Annu. Rev. Neurosci.35, 203 (2012)

    work page 2012

  7. [7]

    Ott and T

    E. Ott and T. M. Antonsen, Low dimensional behavior of large systems of globally coupled oscillators, Chaos18, 037113 (2008)

    work page 2008

  8. [8]

    Watanabe and S

    S. Watanabe and S. H. Strogatz, Integrability of a glob- ally coupled oscillator array, Phys. Rev. Lett.70, 2391 (1993)

    work page 1993

  9. [9]

    Ott and T

    E. Ott and T. M. Antonsen, Long time evolution of phase oscillator systems, Chaos19, 023117 (2009)

    work page 2009

  10. [10]

    Montbri´ o, D

    E. Montbri´ o, D. Paz´ o, and A. Roxin, Macroscopic de- scription for networks of spiking neurons, Phys. Rev. X 5, 021028 (2015)

    work page 2015

  11. [11]

    Paz´ o and E

    D. Paz´ o and E. Montbri´ o, Low-dimensional dynamics of populations of pulse-coupled oscillators, Phys. Rev. X4, 011009 (2014)

    work page 2014

  12. [12]

    C. R. Laing, Derivation of a neural field model from a net- work of theta neurons, Phys. Rev. E90, 010901 (2014)

    work page 2014

  13. [13]

    C. Bick, M. Goodfellow, C. R. Laing, and E. A. Martens, Understanding the dynamics of biological and neural os- cillator networks through exact mean-field reductions: A review, J. Math. Neurosci.10, 9 (2020)

    work page 2020

  14. [14]

    Ottino-L¨ offler and S

    B. Ottino-L¨ offler and S. H. Strogatz, Volcano transition in a solvable model of frustrated oscillators, Physical Re- view Letters120, 264102 (2018)

    work page 2018

  15. [15]

    Paz´ o and R

    D. Paz´ o and R. Gallego, Volcano transition in popula- tions of phase oscillators with random nonreciprocal in- teractions, Physical Review E108, 014202 (2023)

    work page 2023

  16. [16]

    Y. Kati, J. Ranft, and B. Lindner, Self-consistent auto- correlation of a disordered kuramoto model in the asyn- chronous state, Physical Review E110, 054301 (2024)

    work page 2024

  17. [17]

    Hong and E

    H. Hong and E. A. Martens, First-order like phase tran- sition induced by quenched coupling disorder, Chaos32, 063125 (2022)

    work page 2022

  18. [18]

    Le´ on and D

    I. Le´ on and D. Paz´ o, Enlarged Kuramoto model: Sec- ondary instability and transition to collective chaos, Phys. Rev. E105, L042201 (2022)

    work page 2022

  19. [19]

    Le´ on and D

    I. Le´ on and D. Paz´ o, Dynamics and chaotic properties of the fully disordered Kuramoto model, Chaos35, 073140 (2025)

    work page 2025

  20. [20]

    C. Bick, M. J. Panaggio, and E. A. Martens, Chaos in Kuramoto oscillator networks, Chaos28, 071102 (2018)

    work page 2018

  21. [21]

    S. Olmi, A. Politi, and A. Torcini, Collective chaos in pulse-coupled neural networks, Europhys. Lett.92, 60007 (2010)

    work page 2010

  22. [22]

    Pr¨ user, S

    A. Pr¨ user, S. Rosmej, and A. Engel, Nature of the vol- cano transition in the fully disordered kuramoto model, Physical Review Letters132, 187201 (2024)

    work page 2024

  23. [23]

    Pr¨ user and A

    A. Pr¨ user and A. Engel, Role of coupling asymmetry in the fully disordered kuramoto model, Physical Review E 12 110, 064214 (2024)

    work page 2024

  24. [24]

    Villain, Theory of one- and two-dimensional magnets with an easy magnetization plane

    J. Villain, Theory of one- and two-dimensional magnets with an easy magnetization plane. ii. the planar, classical, two-dimensional magnet, Journal de Physique36, 581 (1975)

    work page 1975

  25. [25]

    L. M. Sieberer, G. Wachtel, E. Altman, and S. Diehl, Lat- tice duality for the compact kardar–parisi–zhang equa- tion, Physical Review B94, 104521 (2016)

    work page 2016

  26. [26]

    T. B. Luke, E. Barreto, and P. So, Complete classification of the macroscopic behavior of a heterogeneous network of theta neurons, Neural Computation25, 3207 (2013)

    work page 2013

  27. [27]

    Brown, J

    E. Brown, J. Moehlis, and P. Holmes, On the phase re- duction and response dynamics of neural oscillator pop- ulations, Neural Computation16, 673 (2004)

    work page 2004

  28. [28]

    Augustin, J

    M. Augustin, J. Ladenbauer, F. Baumann, and K. Ober- mayer, Low-dimensional spike rate models derived from networks of adaptive integrate-and-fire neurons: Com- parison and implementation, PLoS Computational Biol- ogy13, e1005545 (2017)

    work page 2017

  29. [29]

    Zerlaut, S

    Y. Zerlaut, S. Chemla, F. Chavane, and A. Destexhe, Modeling mesoscopic cortical dynamics using a mean- field model of conductance-based networks of adaptive exponential integrate-and-fire neurons, Journal of Com- putational Neuroscience44, 45 (2018)

    work page 2018

  30. [30]

    di Volo, A

    M. di Volo, A. Romagnoni, C. Capone, and A. Destexhe, Biologically realistic mean-field models of conductance- based networks of spiking neurons with adaptation, Neu- ral Computation31, 653 (2019)

    work page 2019

  31. [31]

    Carlu, O

    M. Carlu, O. Chehab, L. Dalla Porta, D. Depan- nemaecker, C. H´ eric´ e, M. Jedynak, E. K¨ oksal Ers¨ oz, P. Muratore, S. Souihel, C. Capone, Y. Zerlaut, A. Des- texhe, and M. di Volo, A mean-field approach to the dy- namics of networks of complex neurons, from nonlinear integrate-and-fire to hodgkin–huxley models, Journal of Neurophysiology123, 1042 (2020)

    work page 2020

  32. [32]

    Daido, Multibranch entrainment and scaling in large populations of coupled oscillators, Physical Review Let- ters77, 1406 (1996)

    H. Daido, Multibranch entrainment and scaling in large populations of coupled oscillators, Physical Review Let- ters77, 1406 (1996)

    work page 1996

  33. [33]

    P. S. Skardal, E. Ott, and J. G. Restrepo, Cluster synchrony in systems of coupled phase oscillators with higher-order coupling, Physical Review E84, 036208 (2011)

    work page 2011

  34. [34]

    P. C. Martin, E. D. Siggia, and H. A. Rose, Statistical dynamics of classical systems, Physical Review A8, 423 (1973)

    work page 1973

  35. [35]

    J. A. Hertz, Y. Roudi, and P. Sollich, Path integral meth- ods for the dynamics of stochastic and disordered sys- tems, J. Phys. A: Math. Theor.50, 033001 (2017)

    work page 2017

  36. [36]

    Sompolinsky, A

    H. Sompolinsky, A. Crisanti, and H.-J. Sommers, Chaos in random neural networks, Physical Review Letters61, 259 (1988)

    work page 1988

  37. [37]

    Crisanti and H

    A. Crisanti and H. Sompolinsky, Path-integral approach to random neural networks, Phys. Rev. E98, 062120 (2018)

    work page 2018

  38. [38]

    Crisanti, H

    A. Crisanti, H. Horner, and H.-J. Sommers, The spherical p-spin interaction spin-glass model, Z. Phys. B92, 257 (1993)

    work page 1993

  39. [39]

    G. B. Ermentrout and D. H. Terman,Mathematical Foundations of Neuroscience, Interdisciplinary Applied Mathematics, Vol. 35 (Springer, New York, 2010)

    work page 2010

  40. [40]

    A. T. Winfree, Biological rhythms and the behavior of populations of coupled oscillators, Journal of Theoretical Biology16, 15 (1967)

    work page 1967

  41. [41]

    T. I. Netoff, M. A. Schwemmer, and T. J. Lewis, Ex- perimentally estimating phase response curves of neu- rons: Theoretical and practical issues, inPhase Re- sponse Curves in Neuroscience: Theory, Experiment, and Analysis, Springer Series in Computational Neuro- science, Vol. 6, edited by N. W. Schultheiss, A. A. Prinz, and R. J. Butera (Springer, New Yor...

    work page 2012

  42. [42]

    N. W. Schultheiss, J. R. Edgerton, and D. Jaeger, Phase response curve analysis of a full morphological globus pallidus neuron model reveals distinct perisomatic and dendritic modes of synaptic integration, Journal of Neu- roscience30, 2767 (2010)

    work page 2010

  43. [43]

    G. B. Ermentrout and N. Kopell, Parabolic bursting in an excitable system coupled with a slow oscillation, SIAM Journal on Applied Mathematics46, 233 (1986)

    work page 1986

  44. [44]

    Brette and W

    R. Brette and W. Gerstner, Adaptive exponential integrate-and-fire model as an effective description of neuronal activity, Journal of Neurophysiology94, 3637 (2005)

    work page 2005

  45. [45]

    J. P. L. Hatchett and T. Uezu, Mean field and cavity analysis for coupled oscillator networks, Phys. Rev. E 78, 036106 (2008)

    work page 2008

  46. [46]

    F. L. Metz, Dynamical mean-field theory of complex sys- tems on sparse directed networks, Phys. Rev. Lett.134, 037401 (2025), arXiv:2406.06346 [cond-mat.dis-nn]

    work page arXiv 2025

  47. [47]

    compact dynamical mean-field theory of oscillator networks

    K. Reddy, Simulation code for “compact dynamical mean-field theory of oscillator networks” (2026)

    work page 2026