BCOMBCOM
CalliopeKnowledge Librarian
WP0123
working_paperongoinginternalcomplete

Rosetta Stone of Neural Mass Models

Francesca Castaldo, Jordi Garcia-Ojalvo, Raul Palma, Pau Clusella, Giulio Ruffini

P1·Computational Neuropsychiatry & NeurophenomenologyL2·MathematicsL6·Brains

Brain dynamics dominate every level of neural organization—from single-neuron spiking to the macroscopic waves captured by functional magnetic resonance imaging (fMRI), magnetoencephalography (MEG), and electroencephalography (EEG) —yet the mathematical tools used to interrogate those dynamics remain scattered across a patchwork of traditions. Neural mass models (NMMs) (aggregate neural models) provide one of the most popular gateways into this landscape, but their sheer variety---spanning lumped parameter models, firing‐rate equations, and multi‐layer generators---demands a unifying framework that situates diverse architectures along a continuum of abstraction and biological detail. Here, we start from the idea that oscillations originate from a simple push-pull interaction between two or more neural populations. We build from the undamped harmonic oscillator and, guided by a simple push–pull motif between excitatory and inhibitory populations, climb a systematic ladder of detail. Each rung is presented first in isolation, next under forcing, and then within a coupled network, reflecting the progression from single‐node to whole‐brain modeling. By transforming a repertoire of disparate formalisms into a navigable ladder, we hope to turn NMM choice from a subjective act into a principled design decision, helping both theorists and experimentalists translate between scales, modalities, and interventions. In doing so, we offer a Rosetta Stone for brain oscillation models—one that lets the field speak a common dynamical language while preserving the dialectical richness that fuels discovery.

Building the Rosetta Stone of Neural Mass Models

The field of brain oscillation modeling has a translation problem. One camp writes down phase oscillators. Another builds firing-rate equations with synaptic filters. A third derives exact mean-field limits from spiking neurons. They're often describing the same phenomenon, but in different notations, with different assumptions buried in different places — so nobody can easily tell you how a Jansen-Rit model relates to a Kuramoto network, or why a Wilson-Cowan model and a Stuart-Landau oscillator are secretly close cousins. This paper builds the dictionary.

The core move is almost embarrassingly simple: every model in the hierarchy is built on a "push-pull" motif — two variables that drive each other in a loop, one increasing the other while being decreased by it in return, like a mass and spring trading energy back and forth, or an excitatory population exciting an inhibitory population that turns around and suppresses it. Start with the cleanest version of this idea, the undamped harmonic oscillator (constant-amplitude, perpetual rotation), and then add complexity one knob at a time: damping (the oscillation decays), forcing (something pushes it from outside), a cubic nonlinearity that caps the amplitude and creates a genuine self-sustaining limit cycle (this gets you the Stuart-Landau oscillator, the normal form of a Hopf bifurcation), then swap in biologically-flavored variables — firing rates and E-I population interactions — to get Wilson-Cowan, then add realistic synaptic filtering to get Jansen-Rit-style models (the paper calls these NMM1), and finally derive the whole thing exactly from a spiking network mean-field limit (NMM2, based on the Montbrió-Pazó-Roxin reduction). Each rung is presented alone, then under forcing, then coupled into a network — mirroring how you'd actually go from a single brain region to a whole-brain model wired up over a connectome.

Why bother with this ladder instead of just picking a favorite model? Because the choice of neural mass model is usually made by convention or convenience, not by matching the model's assumptions to the question being asked. If you want closed-form predictions of functional connectivity and power spectra, the linear damped-oscillator regime is enough and it's fully analytic — no simulation needed. If you want self-sustained rhythms with amplitude regulation, multistability, or turbulence-like whole-brain dynamics (useful for modeling wakefulness vs. sedation, or psychedelic states), you need the nonlinear Stuart-Landau regime. If you want to connect all the way down to what a population of spiking neurons is actually doing, you need NMM2. The paper explicitly doesn't claim one grand unification — it calls these "local bridges," valid near specific bifurcations or in specific limits — but it makes the connective tissue between them explicit and puts them all in one notation.

One especially useful contribution is the treatment of "forcing" as a single unified concept covering everything from network coupling, to external stimulation (tACS, TMS, DBS), to background noise, to pharmacological modulation. Because everything reduces locally to the same linearized damped-oscillator equation near a stable state, the paper shows that stimulation and pharmacology are just two different ways of perturbing the same underlying linear response function — one perturbs via an external drive term, the other via a shift in the system's damping/coupling parameters. That's a nice unifying insight: your drug and your electrode are, mathematically, probing the same resolvent.

The paper is explicit about scope: it's not a full survey of every neural mass tradition. Conductance-

Zenodo
10.5281/zenodo.21008779
DOI
10.1016/j.physrep.2026.05.004
Preprint
https://arxiv.org/abs/2512.10982
Publication
https://doi.org/10.1016/j.physrep.2026.05.004
WP ID
WP0123
Lifecycle
ongoing
Visibility
internal
Access level
open
Embargo until
Priority
Collab
closed
Venue
Physics Reports
DOI
10.1016/j.physrep.2026.05.004
Deadline
Owner
Source
drive_legacy
Repo path
WP0123
  • v1.0.0 (publication) · external-source
    Auto-created by update_metadata to host current_venue / current_doi (the recompute_paper_denorm trigger reads these from paper_versions, not papers).
  • 0.1.0 (draft) · auto-run-placeholder · zenodo:21008780