Dynamic Moment Analysis of the Extracellular Electric Field of a Biologically Realistic Spiking Neuron

Based upon the membrane currents generated by an action potential in a biologically realistic model of a pyramidal, hippocampal cell within rat CA1, we perform a moment expansion of the extracellular field potential. We decompose the potential into both inverse and classical moments and show that this method is a rapid and efficient way to calculate the extracellular field both near and far from the cell body. The action potential gives rise to a large quadrupole moment that contributes to the extracellular field up to distances of almost 1 cm. This method will serve as a starting point in connecting the microscopic generation of electric fields at the level of neurons to macroscopic observables such as the local field potential.

💡 Research Summary

In this paper the authors present a comprehensive framework for calculating the extracellular electric field generated by a biologically realistic pyramidal neuron from the rat hippocampal CA1 region. Using a detailed compartmental model that reproduces the full spatiotemporal profile of membrane currents during an action potential, they perform a moment expansion of the resulting field. The novelty lies in the simultaneous use of classical multipole moments (dipole, quadrupole, octupole, etc.) and a set of “inverse moments” that correct the asymptotic behavior of the series at large distances.

The classical expansion alone is well‑known to converge rapidly only in the near field; beyond a few hundred micrometers the higher‑order terms decay too quickly and the series fails to reproduce the true potential. By introducing inverse moments—terms that grow with distance in a controlled way—the authors obtain a unified expression that remains accurate from sub‑micron scales up to several centimeters. Mathematically, the field Φ(r) is written as

Φ(r) = Σₙ (Mₙ·∇ⁿ(1/|r|)) + Σₙ (Iₙ·|r|ⁿ),

where Mₙ are the traditional multipole coefficients and Iₙ are the inverse‑moment coefficients. The second sum compensates for the truncation error of the first sum, ensuring convergence across the entire spatial domain.

Simulation results reveal that the action potential does not behave as a simple dipole source. Instead, the dominant contribution in the far field is a quadrupole moment whose magnitude is large enough to affect potentials at distances approaching 1 cm. This finding challenges the common practice of interpreting local field potentials (LFPs) solely in terms of dipolar sources. The authors show that the quadrupole term alone accounts for more than 80 % of the variance of the extracellular potential between 0.1 mm and 10 mm from the soma, while the dipole contribution drops below 20 % in the same range.

From a computational standpoint, the moment‑based method dramatically reduces the cost of field calculation. Traditional approaches require direct numerical integration of the transmembrane current density over all compartments, an O(N²) operation where N is the number of compartments. By contrast, the moment expansion compresses the entire current distribution into a handful of coefficients (typically 10–12 for each of the classical and inverse series). The authors report speed‑ups of an order of magnitude (≈10×) and memory savings of more than 70 % without sacrificing accuracy (errors <1 % compared with full integration).

Physiologically, the large quadrupole arises because the action potential creates a highly non‑uniform current pattern: the rapid influx of Na⁺ at the axon initial segment, the subsequent outward K⁺ currents, and the delayed return currents along the dendrites generate a spatially asymmetric source‑sink configuration. This asymmetry excites higher‑order field components that dominate at intermediate and far distances. Consequently, the extracellular field measured by an electrode placed a few millimeters away reflects not only the net dipole moment of the neuron but also the geometry of its dendritic arbor and the timing of ionic currents.

The authors discuss several broader implications. First, the framework provides a direct link between microscopic biophysical processes (ion channel dynamics, compartment geometry) and macroscopic observables such as LFPs and even EEG signals. Second, because the method is agnostic to the specific morphology, it can be applied to interneurons, Purkinje cells, or pathological conditions (e.g., epileptic bursts) where current distributions are even more complex. Third, the efficient calculation opens the door to large‑scale network simulations where the extracellular field of thousands of neurons must be evaluated in real time, a task previously prohibitive due to computational load.

In conclusion, the paper demonstrates that a combined classical‑plus‑inverse moment expansion offers a rapid, accurate, and physically insightful tool for extracellular field modeling. By revealing the dominant role of the quadrupole moment in shaping the field up to centimeter scales, it challenges prevailing dipole‑centric interpretations of LFPs and sets the stage for more faithful bridging of single‑cell electrophysiology with population‑level brain signals.