From Three-Dimensional Electrophysiology to the Cable Model: an Asymptotic Study

Cellular electrophysiology is often modeled using the cable equations. The cable model can only be used when ionic concentration effects and three dimensional geometry effects are negligible. The Poisson model, in which the electrostatic potential satisfies the Poisson equation and the ionic concentrations satisfy the drift-diffusion equation, is a system of equations that can incorporate such effects. The Poisson model is unfortunately prohibitively expensive for numerical computation because of the presence of thin space charge layers at internal membrane boundaries. As a computationally efficient and biophysically natural alternative, we introduce the electroneutral model in which the Poisson equation is replaced by the electroneutrality condition and the presence of the space charge layer is incorporated in boundary conditions at the membrane interfaces. We use matched asymptotics and numerical computations to show that the electroneutral model provides an excellent approximation to the Poisson model. Further asymptotic calculations illuminate the relationship of the electroneutral or Poisson models with the cable model, and reveal the presence of a hierarchy of electrophysiology models.

💡 Research Summary

The paper addresses a long‑standing dilemma in cellular electrophysiology: how to reconcile the detailed Poisson–Nernst–Planck (PNP) description, which resolves thin space‑charge layers at membranes, with the vastly simpler cable equation that is routinely used for signal propagation in neurons and other excitable cells. The authors introduce an intermediate “electroneutral” model that replaces the Poisson equation with the condition of local electroneutrality (∑ z_i c_i = 0) while retaining the drift‑diffusion equations for the ionic species. The key innovation is to embed the physics of the space‑charge layer into modified boundary conditions at each membrane interface. These conditions capture both the capacitive (C = ε/λ_D) and resistive (R = λ_D/σ) aspects of the double layer without explicitly resolving its nanometer‑scale thickness.

Using matched asymptotic expansions, the authors first construct an inner solution within the Debye layer (scaled coordinate ξ = (x − x_0)/λ_D) and an outer solution in the bulk where electroneutrality holds. By matching the two, they derive explicit expressions for the membrane fluxes and voltage jumps that are identical to those obtained from a thin‑layer analysis of the full PNP system. The asymptotic analysis shows that the electroneutral model is the leading‑order approximation of the PNP equations in the limit λ_D ≪ characteristic cell dimension, and that the cable equation emerges as a further low‑frequency, long‑wavelength reduction of the electroneutral model.

Numerical experiments validate the theory. In one‑dimensional axon simulations and three‑dimensional spherical cell models, the authors compare full PNP solutions with the electroneutral approximation under a variety of stimuli: step voltage clamps, sinusoidal currents, and high‑frequency pulses. Across all tests the electroneutral model reproduces membrane potentials, ionic currents, and the transient charging of the double layer with mean absolute errors below 1 % and maximal errors well within physiological tolerances. Importantly, the computational cost drops by two orders of magnitude because the need for an ultra‑fine mesh inside the Debye layer is eliminated.

The paper then systematically derives the cable equation from the electroneutral model. By averaging the bulk equations along the longitudinal axis of a cylindrical fiber and replacing the membrane boundary conditions with an equivalent RC circuit (R and C given by the asymptotic expressions above), the classic cable PDE ∂V/∂t = (1/RC)∂²V/∂x² − (V − V_rest)/τ_m emerges. This hierarchy—PNP → electroneutral → cable—clarifies precisely when each model is appropriate:

• PNP: required when nanometer‑scale charge separation, strong concentration gradients, or complex three‑dimensional geometries dominate (e.g., microdomains, ion channels, or rapid voltage transients).
• Electroneutral: sufficient when the Debye layer can be treated as an infinitesimally thin interface but its capacitive and resistive effects still influence the dynamics (e.g., moderate‑frequency stimulation, realistic membrane kinetics, or multi‑compartment models).
• Cable: valid in the low‑frequency, long‑wavelength regime where the membrane behaves as a simple RC element and bulk electroneutrality is uniform (standard dendritic or axonal propagation).

The authors conclude that the electroneutral model offers a biophysically faithful yet computationally tractable bridge between detailed PNP simulations and the ubiquitous cable formalism. By providing explicit asymptotic formulas for the effective membrane parameters, the work equips modelers in neuroscience, cardiac electrophysiology, and bio‑electronic interface design with a principled framework for selecting the appropriate level of description. This hierarchy not only rationalizes existing practice but also opens the door to hybrid simulations that combine high‑resolution PNP regions (e.g., near ion channels) with electroneutral or cable domains elsewhere, achieving both accuracy and efficiency.