Rosetta Stone of Neural Mass Models

TL;DR

The paper addresses the fragmentation of neural mass modeling by proposing a unifying ladder that connects simple harmonic motion to full mean-field models through a push–pull motif between interacting populations. Starting from a linear harmonic oscillator and progressively adding damping, forcing, nonlinearities, and biologically grounded synaptic transfer, it derives principled links between phase- and amplitude-based formalisms (HO, WC, SL, NMM1, NMM2) and demonstrates how networks respond to external perturbations and stimulation. Key contributions include exact mean-field reductions for QIF-based NMM2, analytic tractable limits for WC/NMM1, and a coherent design path for translating between scales and modalities, enabling principled model selection and in silico perturbations. The framework clarifies how linear resonators underlie functional connectivity, how Hopf bifurcations govern regime transitions, and how synaptic kinetics shape rhythms, with practical implications for neuromodulation, clinical interventions, and multimodal data interpretation.

Abstract

Brain dynamics dominate every level of neural organization -- from single-neuron spiking to the macroscopic waves captured by fMRI, MEG, and 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 \emph{Rosetta Stone} for brain oscillation models -- one that lets the field speak a common dynamical language while preserving the dialectical richness that fuels discovery.

Figures (13)

• Figure 1.1: Oscillatory dynamics across increasing model complexity. Top row: phase-space portraits. Bottom row: qualitative one-parameter bifurcation diagrams. The horizontal arrow indicates increasing model complexity from left to right. Blue traces highlight attracting (or neutrally stable) periodic orbits and their stable cycle branches; gray traces denote unstable invariant sets and illustrative transients. Panels: (a) phase-only oscillator ( no amplitude dynamics); (b) damped oscillator with a globally attracting fixed point; (c) Stuart–Landau oscillator (SL) where a Hopf bifurcation creates a stable limit cycle; (d) Wilson–Cowan (WILCO) $E$--$I$ model with coexistence of equilibria and oscillations; (e–f) two neural-mass models (NMM1, NMM2) showing parameter-dependent onset and growth of oscillations and possible bistability. Sketches are schematic and not to scale.
• Figure 4.1: Post synaptic potentials in the first-order or the more realistic second-order synapse models. These plots are the response to an "impulse" at time zero, i.e., the solutions to $\hat{L}[x] = \delta (t)$.
• Figure 4.2: Simple recurrent population model. Here, the population receives an external current $I_\text{ext}$ and self-input via a synapse $s$, and outputs it firing rate $r$. Two time scales (delays) enter: the membrane response time constant $\tau_m$ for updating the firing rate $r$, and the synaptic current time constant $\tau_s$. See Eq. \ref{['eq:simpleNMM']} in the text for details.
• Figure 6.1: Four models using the second-order formalism: a) PING-like push–pull motif, b) Jansen–Rit model jansen_neurophysiologically-based_1993grimbert_analysis_2006, c) Wendling model wendling_computational_2016, d) Laminar model ruffiniP118BiophysicallyRealistic2020sanchez-todoPhysicalNeuralMass2023. Synapses are shown as arrowheads (excitatory) or buttons (inhibitory). Noise or external inputs are indicated on pyramidal cells, although other targets ($\mu$) are also possible.
• Figure 7.1: NMM2 diagrams. (A): Diagram for self-coupled population with connectivity $J$ receiving and external input $Q$. (B): generalization for multiple populations. (C): A generic two-population model.