title: Computationally efficient simulation of electrical activity at cell membranes interacting with self-generated and externally imposed electric fields.
abstract: OBJECTIVE
  We present a computational method that implements a reduced set of Maxwell's equations to allow simulation of cells under realistic conditions: sub-micron cell morphology, a conductive non-homogeneous space and various ion channel properties and distributions.
  
  
  METHODS
  While a reduced set of Maxwell's equations can be used to couple membrane currents to extra- and intracellular potentials, this approach is rarely taken, most likely because adequate computational tools are missing. By using these equations, and introducing an implicit solver, numerical stability is attained even with large time steps. The time steps are limited only by the time development of the membrane potentials.
  
  
  RESULTS
  This method allows simulation times of tens of minutes instead of weeks, even for complex problems. The extracellular fields are accurately represented, including secondary fields, which originate at inhomogeneities of the extracellular space and can reach several millivolts. We present a set of instructive examples that show how this method can be used to obtain reference solutions for problems, which might not be accurately captured by the traditional approaches. This includes the simulation of realistic magnitudes of extracellular action potential signals in restricted extracellular space.
  
  
  CONCLUSIONS
  The electric activity of neurons creates extracellular potentials. Recent findings show that these endogenous fields act back onto the neurons, contributing to the synchronization of population activity. The influence of endogenous fields is also relevant for understanding therapeutic approaches such as transcranial direct current, transcranial magnetic and deep brain stimulation. The mutual interaction between fields and membrane currents is not captured by today's concepts of cellular electrophysiology, including the commonly used activation function, as those concepts are based on isolated membranes in an infinite, isopotential extracellular space. The presented tool makes simulations with detailed morphology and implicit interactions of currents and fields available to the electrophysiology community.
authors.affiliation: Max Planck Institute for Dynamics and Self-Organization, Göttingen, Germany. agdtoro@nld.ds.mpg.de
authors.forename: Andres
authors.initials: A
authors.lastname: Agudelo-Toro
authors.forename: Andreas
authors.initials: A
authors.lastname: Neef
mesh_terms: D000200:Action Potentials / Q000502:physiology; D000465:Algorithms; D002462:Cell Membrane / Q000502:physiology*; D003198:Computer Simulation; D046690:Deep Brain Stimulation; D004574:Electromagnetic Fields; D055724:Electrophysiological Phenomena; D005071:Evoked Potentials / Q000502:physiology; D005110:Extracellular Space / Q000502:physiology; D006801:Humans; D008959:Models, Neurological; D009431:Neural Conduction / Q000502:physiology; D009474:Neurons / Q000502:physiology / Q000648:ultrastructure; D017729:Presynaptic Terminals / Q000502:physiology / Q000648:ultrastructure; D012984:Software; D050781:Transcranial Magnetic Stimulation