From Maxwell's equations to the theory of current-source density analysis

Electrocardiographic approach strategies in patients with Parkinson disease treated with deep brain stimulation

Deep brain stimulation (DBS) is an interdisciplinary and reversible therapy that uses high-frequency electrical stimulation to correct aberrant neural pathways in motor and cognitive neurological disorders. However, the high frequency of the waves used in DBS can interfere with electrical recording devices (e.g., electrocardiogram, electroencephalogram, cardiac monitor), creating artifacts that hinder their interpretation. The compatibility of DBS with these devices varies and depends on factors such as the underlying disease and the configuration of the neurostimulator. In emergencies where obtaining an electrocardiogram is crucial, the need for more consensus on reducing electrical artifacts in patients with DBS becomes a significant challenge. Various strategies have been proposed to attenuate the artifact generated by DBS, such as changing the DBS configuration from monopolar to bipolar, temporarily deactivating DBS during electrocardiographic recording, applying frequency filters both lower and higher than those used by DBS, and using non-standard leads. However, the inexperience of medical personnel, variability in DBS models, or the lack of a controller at the time of approach limit the application of these strategies. Current evidence on their reproducibility and efficacy is limited. Due to the growing elderly population and the rising utilization of DBS, it is imperative to create electrocardiographic methods that are easily accessible and reproducible for general physicians and emergency services.

kCSD-python, reliable current source density estimation with quality control

Interpretation of extracellular recordings can be challenging due to the long range of electric field. This challenge can be mitigated by estimating the current source density (CSD). Here we introduce kCSD-python, an open Python package implementing Kernel Current Source Density (kCSD) method and related tools to facilitate CSD analysis of experimental data and the interpretation of results. We show how to counter the limitations imposed by noise and assumptions in the method itself. kCSD-python allows CSD estimation for an arbitrary distribution of electrodes in 1D, 2D, and 3D, assuming distributions of sources in tissue, a slice, or in a single cell, and includes a range of diagnostic aids. We demonstrate its features in a Jupyter Notebook tutorial which illustrates a typical analytical workflow and main functionalities useful in validating analysis results.

Spontaneous slow wave oscillations in extracellular field potential recordings reflect the alternating dominance of excitation and inhibition

In the resting state, cortical neurons can fire action potentials spontaneously but synchronously (Up state), followed by a quiescent period (Down state) before the cycle repeats. Extracellular recordings in the infragranular layer of cortex with a micro-electrode display a negative deflection (depth-negative) during Up states and a positive deflection (depth-positive) during Down states. The resulting slow wave oscillation (SWO) has been studied extensively during sleep and under anaesthesia. However, recent research on the balanced nature of synaptic excitation and inhibition has highlighted our limited understanding of its genesis. Specifically, are excitation and inhibition balanced during SWOs? We analyse spontaneous local field potentials (LFPs) during SWOs recorded from anaesthetised rats via a multi-channel laminar micro-electrode and show that the Down state consists of two distinct synaptic states: a Dynamic Down state associated with depth-positive LFPs and a prominent dipole in the extracellular field, and a Static Down state with negligible (≈ 0 mV $ \approx 0{\mathrm{\;mV}}$ ) LFPs and a lack of dipoles extracellularly. We demonstrate that depth-negative and -positive LFPs are generated by a shift in the balance of synaptic excitation and inhibition from excitation dominance (depth-negative) to inhibition dominance (depth-positive) in the infragranular layer neurons. Thus, although excitation and inhibition co-tune overall, differences in their timing lead to an alternation of dominance, manifesting as SWOs. We further show that Up state initiation is significantly faster if the preceding Down state is dynamic rather than static. Our findings provide a coherent picture of the dependence of SWOs on synaptic activity. KEY POINTS: Cortical neurons can exhibit repeated cycles of spontaneous activity interleaved with periods of relative silence, a phenomenon known as 'slow wave oscillation' (SWO). During SWOs, recordings of local field potentials (LFPs) in the neocortex show depth-negative deflection during the active period (Up state) and depth-positive deflection during the silent period (Down state). Here we further classified the Down state into a dynamic phase and a static phase based on a novel method of classification and revealed non-random, stereotypical sequences of the three states occurring with significantly different transitional kinetics. Our results suggest that the positive and negative deflections in the LFP reflect the shift of the instantaneous balance between excitatory and inhibitory synaptic activity of the local cortical neurons. The differences in transitional kinetics may imply distinct synaptic mechanisms for Up state initiation. The study may provide a new approach for investigating spontaneous brain rhythms.

An electrodiffusive neuron-extracellular-glia model for exploring the genesis of slow potentials in the brain

Within the computational neuroscience community, there has been a focus on simulating the electrical activity of neurons, while other components of brain tissue, such as glia cells and the extracellular space, are often neglected. Standard models of extracellular potentials are based on a combination of multicompartmental models describing neural electrodynamics and volume conductor theory. Such models cannot be used to simulate the slow components of extracellular potentials, which depend on ion concentration dynamics, and the effect that this has on extracellular diffusion potentials and glial buffering currents. We here present the electrodiffusive neuron-extracellular-glia (edNEG) model, which we believe is the first model to combine compartmental neuron modeling with an electrodiffusive framework for intra- and extracellular ion concentration dynamics in a local piece of neuro-glial brain tissue. The edNEG model (i) keeps track of all intraneuronal, intraglial, and extracellular ion concentrations and electrical potentials, (ii) accounts for action potentials and dendritic calcium spikes in neurons, (iii) contains a neuronal and glial homeostatic machinery that gives physiologically realistic ion concentration dynamics, (iv) accounts for electrodiffusive transmembrane, intracellular, and extracellular ionic movements, and (v) accounts for glial and neuronal swelling caused by osmotic transmembrane pressure gradients. The edNEG model accounts for the concentration-dependent effects on ECS potentials that the standard models neglect. Using the edNEG model, we analyze these effects by splitting the extracellular potential into three components: one due to neural sink/source configurations, one due to glial sink/source configurations, and one due to extracellular diffusive currents. Through a series of simulations, we analyze the roles played by the various components and how they interact in generating the total slow potential. We conclude that the three components are of comparable magnitude and that the stimulus conditions determine which of the components that dominate.

What we can and what we cannot see with extracellular multielectrodes

Extracellular recording is an accessible technique used in animals and humans to study the brain physiology and pathology. As the number of recording channels and their density grows it is natural to ask how much improvement the additional channels bring in and how we can optimally use the new capabilities for monitoring the brain. Here we show that for any given distribution of electrodes we can establish exactly what information about current sources in the brain can be recovered and what information is strictly unobservable. We demonstrate this in the general setting of previously proposed kernel Current Source Density method and illustrate it with simplified examples as well as using evoked potentials from the barrel cortex obtained with a Neuropixels probe and with compatible model data. We show that with conceptual separation of the estimation space from experimental setup one can recover sources not accessible to standard methods.

Every set of measurements is a window into reality rendering its incomplete or distorted picture. It is often difficult to relate the obtained representation of the world to underlying ground truth. Here we show, for brain electrophysiology, for arbitrary experimental setup (distribution of electrodes), and arbitrary analytical setup (function space of current source densities), that one can identify distributions of current sources which can be recovered precisely, and those which are invisible in the system. This shows what is and what is not observable in the studied system for a given setup, allows to improve the analysis results by modifying analytical setup, and facilitates interpretation of the measured sets of LFP, ECoG and EEG recordings. In particular we show that with conceptual separation of the estimation space from experimental setup one can recover source distributions not accessible to standard methods.

Lag-invariant detection of interactions in spatially-extended systems using linear inverse modeling

Measurements on physical systems result from the systems' activity being converted into sensor measurements by a forward model. In a number of cases, inversion of the forward model is extremely sensitive to perturbations such as sensor noise or numerical errors in the forward model. Regularization is then required, which introduces bias in the reconstruction of the systems' activity. One domain in which this is particularly problematic is the reconstruction of interactions in spatially-extended complex systems such as the human brain. Brain interactions can be reconstructed from non-invasive measurements such as electroencephalography (EEG) or magnetoencephalography (MEG), whose forward models are linear and instantaneous, but have large null-spaces and high condition numbers. This leads to incomplete unmixing of the forward models and hence to spurious interactions. This motivated the development of interaction measures that are exclusively sensitive to lagged, i.e. delayed interactions. The drawback of such measures is that they only detect interactions that have sufficiently large lags and this introduces bias in reconstructed brain networks. We introduce three estimators for linear interactions in spatially-extended systems that are uniformly sensitive to all lags. We derive some basic properties of and relationships between the estimators and evaluate their performance using numerical simulations from a simple benchmark model.

Local and Volume-Conducted Contributions to Cortical Field Potentials

Brain field potentials (FPs) can reach far from their sources, making difficult to know which waves come from where. We show that modern algorithms efficiently segregate the local and remote contributions to cortical FPs by recovering the generator-specific spatial voltage profiles. We investigated experimentally and numerically the local and remote origin of FPs in different cortical areas in anesthetized rats. All cortices examined show significant state, layer, and region dependent contribution of remote activity, while the voltage profiles help identify their subcortical or remote cortical origin. Co-activation of different cortical modules can be discriminated by the distinctive spatial features of the corresponding profiles. All frequency bands contain remote activity, thus influencing the FP time course, in cases drastically. The reach of different FP patterns is boosted by spatial coherence and curved geometry of the sources. For instance, slow cortical oscillations reached the entire brain, while hippocampal theta reached only some portions of the cortex. In anterior cortices, most alpha oscillations have a remote origin, while in the visual cortex the remote theta and gamma even surpass the local contribution. The quantitative approach to local and distant FP contributions helps to refine functional connectivity among cortical regions, and their relation to behavior.

An electrodiffusive, ion conserving Pinsky-Rinzel model with homeostatic mechanisms

In most neuronal models, ion concentrations are assumed to be constant, and effects of concentration variations on ionic reversal potentials, or of ionic diffusion on electrical potentials are not accounted for. Here, we present the electrodiffusive Pinsky-Rinzel (edPR) model, which we believe is the first multicompartmental neuron model that accounts for electrodiffusive ion concentration dynamics in a way that ensures a biophysically consistent relationship between ion concentrations, electrical charge, and electrical potentials in both the intra- and extracellular space. The edPR model is an expanded version of the two-compartment Pinsky-Rinzel (PR) model of a hippocampal CA3 neuron. Unlike the PR model, the edPR model includes homeostatic mechanisms and ion-specific leakage currents, and keeps track of all ion concentrations (Na⁺, K⁺, Ca²⁺, and Cl⁻), electrical potentials, and electrical conductivities in the intra- and extracellular space. The edPR model reproduces the membrane potential dynamics of the PR model for moderate firing activity. For higher activity levels, or when homeostatic mechanisms are impaired, the homeostatic mechanisms fail in maintaining ion concentrations close to baseline, and the edPR model diverges from the PR model as it accounts for effects of concentration changes on neuronal firing. We envision that the edPR model will be useful for the field in three main ways. Firstly, as it relaxes commonly made modeling assumptions, the edPR model can be used to test the validity of these assumptions under various firing conditions, as we show here for a few selected cases. Secondly, the edPR model should supplement the PR model when simulating scenarios where ion concentrations are expected to vary over time. Thirdly, being applicable to conditions with failed homeostasis, the edPR model opens up for simulating a range of pathological conditions, such as spreading depression or epilepsy.

Neurons generate their electrical signals by letting ions pass through their membranes. Despite this fact, most models of neurons apply the simplifying assumption that ion concentrations remain effectively constant during neural activity. This assumption is often quite good, as neurons contain a set of homeostatic mechanisms that make sure that ion concentrations vary quite little under normal circumstances. However, under some conditions, these mechanisms can fail, and ion concentrations can vary quite dramatically. Standard models are thus not able to simulate such conditions. Here, we present what to our knowledge is the first multicompartmental neuron model that accounts for ion concentration variations in a way that ensures complete and consistent ion concentration and charge conservation. In this work, we use the model to explore under which activity conditions the ion concentration variations become important for predicting the neurodynamics. We expect the model to be of great value for the field of neuroscience, as it can be used to simulate a range of pathological conditions, such as spreading depression or epilepsy, which are associated with large changes in extracellular ion concentrations.

Finite Element Simulation of Ionic Electrodiffusion in Cellular Geometries

Mathematical models for excitable cells are commonly based on cable theory, which considers a homogenized domain and spatially constant ionic concentrations. Although such models provide valuable insight, the effect of altered ion concentrations or detailed cell morphology on the electrical potentials cannot be captured. In this paper, we discuss an alternative approach to detailed modeling of electrodiffusion in neural tissue. The mathematical model describes the distribution and evolution of ion concentrations in a geometrically-explicit representation of the intra- and extracellular domains. As a combination of the electroneutral Kirchhoff-Nernst-Planck (KNP) model and the Extracellular-Membrane-Intracellular (EMI) framework, we refer to this model as the KNP-EMI model. Here, we introduce and numerically evaluate a new, finite element-based numerical scheme for the KNP-EMI model, capable of efficiently and flexibly handling geometries of arbitrary dimension and arbitrary polynomial degree. Moreover, we compare the electrical potentials predicted by the KNP-EMI and EMI models. Finally, we study ephaptic coupling induced in an unmyelinated axon bundle and demonstrate how the KNP-EMI framework can give new insights in this setting.

ELFENN: A Generalized Platform for Modeling Ephaptic Coupling in Spiking Neuron Models

The transmembrane ionic currents that underlie changes in a cell's membrane potential give rise to electric fields in the extracellular space. In the context of brain activity, these electric fields form the basis for extracellularly recorded signals, such as multiunit activity, local field potentials and electroencephalograms. Understanding the underlying neuronal dynamics and localizing current sources using these signals is often challenging, and therefore effective computational modeling approaches are critical. Typically, the electric fields from neural activity are modeled in a post-hoc form, i.e., a traditional neuronal model is used to first generate the membrane currents, which in turn are then used to calculate the electric fields. When the conductivity of the extracellular space is high, the electric fields are weak, and therefore treating membrane currents and electric fields separately is justified. However, in brain regions of lower conductivity, extracellular fields can feed back and significantly influence the underlying transmembrane currents and dynamics of nearby neurons—this is often referred to as ephaptic coupling. The closed-loop nature of ephaptic coupling cannot be modeled using the post-hoc approaches implemented by existing software tools; instead, electric fields and neuronal dynamics must be solved simultaneously. To this end, we have developed a generalized modeling toolbox for studying ephaptic coupling in compartmental neuron models: E_LFENN (E_Lectric Field Effects in Neural Networks). In open loop conditions, we validate the separate components of E_LFENN for modeling membrane dynamics and associated field potentials against standard approaches (NEURON and LFPy). Unlike standard approaches however, E_LFENN enables the closed-loop condition to be modeled as well, in that the field potentials can feed back and influence membrane dynamics. As an example closed-loop case, we use E_LFENN to study phase-locking of action potentials generated by a population of axons running parallel in a bundle. Being able to efficiently explore ephaptic coupling from a computational perspective using tools, such as E_LFENN will allow us to better understand the physical basis of electric fields in the brain, as well as the conditions in which these fields may influence neuronal dynamics in general.

V1 microcircuit dynamics: altered signal propagation suggests intracortical origins for adaptation in response to visual repetition

Repetitive visual stimulation profoundly changes sensory processing in the primary visual cortex (V1). We show how the associated adaptive changes are linked to an altered flow of synaptic activation across the V1 laminar microcircuit. Using repeated visual stimulation, we recorded layer-specific responses in V1 of two fixating monkeys. We found that repetition-related spiking suppression was most pronounced outside granular V1 layers that receive the main retinogeniculate input. This repetition-related response suppression was robust to alternating stimuli between the eyes, in line with the notion that repetition-related adaptation is predominantly of cortical origin. Most importantly, current source density (CSD) analysis, which provides an estimate of local net depolarization, revealed that synaptic processing during repeated stimulation was most profoundly affected within supragranular layers, which harbor the bulk of cortico-cortical connections. Direct comparison of the temporal evolution of laminar CSD and spiking activity showed that stimulus repetition first affected supragranular synaptic currents, which translated into a reduction of stimulus-evoked spiking across layers. Together, these results suggest that repetition induces an altered state of intracortical processing that underpins visual adaptation. NEW & NOTEWORTHY Our survival depends on our brains rapidly adapting to ever changing environments. A well-studied form of adaptation occurs whenever we encounter the same or similar stimuli repeatedly. We show that this repetition-related adaptation is supported by systematic changes in the flow of sensory activation across the laminar cortical microcircuitry of primary visual cortex. These results demonstrate how adaptation impacts neuronal interactions across cortical circuits.

A Kirchhoff-Nernst-Planck framework for modeling large scale extracellular electrodiffusion surrounding morphologically detailed neurons

Many pathological conditions, such as seizures, stroke, and spreading depression, are associated with substantial changes in ion concentrations in the extracellular space (ECS) of the brain. An understanding of the mechanisms that govern ECS concentration dynamics may be a prerequisite for understanding such pathologies. To estimate the transport of ions due to electrodiffusive effects, one must keep track of both the ion concentrations and the electric potential simultaneously in the relevant regions of the brain. Although this is currently unfeasible experimentally, it is in principle achievable with computational models based on biophysical principles and constraints. Previous computational models of extracellular ion-concentration dynamics have required extensive computing power, and therefore have been limited to either phenomena on very small spatiotemporal scales (micrometers and milliseconds), or simplified and idealized 1-dimensional (1-D) transport processes on a larger scale. Here, we present the 3-D Kirchhoff-Nernst-Planck (KNP) framework, tailored to explore electrodiffusive effects on large spatiotemporal scales. By assuming electroneutrality, the KNP-framework circumvents charge-relaxation processes on the spatiotemporal scales of nanometers and nanoseconds, and makes it feasible to run simulations on the spatiotemporal scales of millimeters and seconds on a standard desktop computer. In the present work, we use the 3-D KNP framework to simulate the dynamics of ion concentrations and the electrical potential surrounding a morphologically detailed pyramidal cell. In addition to elucidating the single neuron contribution to electrodiffusive effects in the ECS, the simulation demonstrates the efficiency of the 3-D KNP framework. We envision that future applications of the framework to more complex and biologically realistic systems will be useful in exploring pathological conditions associated with large concentration variations in the ECS.

Investigating large-scale brain dynamics using field potential recordings: analysis and interpretation

New technologies to record electrical activity from the brain on a massive scale offer tremendous opportunities for discovery. Electrical measurements of large-scale brain dynamics, termed field potentials, are especially important to understanding and treating the human brain. Here, our goal is to provide best practices on how field potential recordings (electroencephalograms, magnetoencephalograms, electrocorticograms and local field potentials) can be analyzed to identify large-scale brain dynamics, and to highlight critical issues and limitations of interpretation in current work. We focus our discussion of analyses around the broad themes of activation, correlation, communication and coding. We provide recommendations for interpreting the data using forward and inverse models. The forward model describes how field potentials are generated by the activity of populations of neurons. The inverse model describes how to infer the activity of populations of neurons from field potential recordings. A recurring theme is the challenge of understanding how field potentials reflect neuronal population activity given the complexity of the underlying brain systems.

Long-term effects of direct current are reproduced by intermittent depolarization of myelinated nerve fibers

Direct current (DC) potently increases the excitability of myelinated afferent fibers in the dorsal columns, both during DC polarization of these fibers and during a considerable (>1 h) postpolarization period. The aim of the present study was to investigate whether similarly long-lasting changes in the excitability of myelinated nerve fibers in the dorsal columns may be evoked by field potentials following stimulation of peripheral afferents and by subthreshold epidurally applied current pulses. The experiments were performed in deeply anesthetized rats. The effects were monitored by changes in nerve volleys evoked in epidurally stimulated hindlimb afferents and in the synaptic actions of these afferents. Both were found to be facilitated during as well as following stimulation of a skin nerve and during as well as following epidurally applied current pulses of 5- to 10-ms duration. The facilitation occurring ≤2 min after skin nerve stimulation could be linked to both primary afferent depolarization and large dorsal horn field potentials, whereas the subsequent changes (up to 1 h) were attributable to effects of the field potentials. The findings lead to the conclusion that the modulation of spinal activity evoked by DC does not require long-lasting polarization and that relatively short current pulses and intrinsic field potentials may contribute to plasticity in spinal activity. These results suggest the possibility of enhancing the effects of epidural stimulation in human subjects by combining it with polarizing current pulses and peripheral afferent stimulation and not only with continuous DC. NEW & NOTEWORTHY The aim of this study was to define conditions under which a long-term increase is evoked in the excitability of myelinated nerve fibers. The results demonstrate that a potent and long-lasting increase in the excitability of afferent fibers traversing the dorsal columns may be induced by synaptically evoked intrinsic field as well as by epidurally applied intermittent current pulses. They thus provide a new means for the facilitation of the effects of epidural stimulation.

Somatotopic Representation of Second Pain in the Primary Somatosensory Cortex of Humans and Rodents

There is now compelling evidence that selective stimulation of Aδ nociceptors eliciting first pain evokes robust responses in the primary somatosensory cortex (S1). In contrast, whether the C-fiber nociceptive input eliciting second pain has an organized projection to S1 remains an open question. Here, we recorded the electrocortical responses elicited by nociceptive-specific laser stimulation of the four limbs in 202 humans (both males and females, using EEG) and 12 freely moving rats (all males, using ECoG). Topographical analysis and source modeling revealed in both species, a clear gross somatotopy of the unmyelinated C-fiber input within the S1 contralateral to the stimulated side. In the human EEG, S1 activity could be isolated as an early-latency negative deflection (C-N1 wave peaking at 710–730 ms) after hand stimulation, but not after foot stimulation because of the spatiotemporal overlap with the subsequent large-amplitude supramodal vertex waves (C-N2/P2). In contrast, because of the across-species difference in the representation of the body surface within S1, S1 activity could be isolated in rat ECoG as a C-N1 after both forepaw and hindpaw stimulation. Finally, we observed a functional dissociation between the generators of the somatosensory-specific lateralized waves (C-N1) and those of the supramodal vertex waves (C-N2/P2), indicating that C-fiber unmyelinated input is processed in functionally distinct somatosensory and multimodal cortical areas. These findings demonstrated that C-fiber input conveys information about the spatial location of noxious stimulation across the body surface, a prerequisite for deploying an appropriate defensive motor repertoire.

SIGNIFICANCE STATEMENT Unmyelinated C-fibers are the evolutionarily oldest peripheral afferents responding to noxious environmental stimuli. Whether C-fiber input conveys information about the spatial location of the noxious stimulation to the primary somatosensory cortex (S1) remains an open issue. In this study, C-fibers were activated by radiant heat stimuli delivered to different parts of the body in both humans and rodents while electrical brain activity was recorded. In both species, the C-fiber peripheral input projects to different parts of the contralateral S1, coherently with the representation of the body surface within this brain region. These findings demonstrate that C-fiber input conveys information about the spatial location of noxious stimulation across the body surface, a prerequisite for deploying an appropriate defensive motor repertoire.

EEG microstates as a tool for studying the temporal dynamics of whole-brain neuronal networks: A review

The present review discusses a well-established method for characterizing resting-state activity of the human brain using multichannel electroencephalography (EEG). This method involves the examination of electrical microstates in the brain, which are defined as successive short time periods during which the configuration of the scalp potential field remains semi-stable, suggesting quasi-simultaneity of activity among the nodes of large-scale networks. A few prototypic microstates, which occur in a repetitive sequence across time, can be reliably identified across participants. Researchers have proposed that these microstates represent the basic building blocks of the chain of spontaneous conscious mental processes, and that their occurrence and temporal dynamics determine the quality of mentation. Several studies have further demonstrated that disturbances of mental processes associated with neurological and psychiatric conditions manifest as changes in the temporal dynamics of specific microstates. Combined EEG-fMRI studies and EEG source imaging studies have indicated that EEG microstates are closely associated with resting-state networks as identified using fMRI. The scale-free properties of the time series of EEG microstates explain why similar networks can be observed at such different time scales. The present review will provide an overview of these EEG microstates, available methods for analysis, the functional interpretations of findings regarding these microstates, and their behavioral and clinical correlates.

Ion diffusion may introduce spurious current sources in current-source density (CSD) analysis

Standard CSD analysis does not account for ionic diffusion. Using biophysically realistic computer simulations, we show that unaccounted-for diffusive currents can lead to the prediction of spurious current sources. This finding may be of strong interest for in vivo electrophysiologists doing extracellular recordings in general, and CSD analysis in particular.

Current-source density (CSD) analysis is a well-established method for analyzing recorded local field potentials (LFPs), that is, the low-frequency part of extracellular potentials. Standard CSD theory is based on the assumption that all extracellular currents are purely ohmic, and thus neglects the possible impact from ionic diffusion on recorded potentials. However, it has previously been shown that in physiological conditions with large ion-concentration gradients, diffusive currents can evoke slow shifts in extracellular potentials. Using computer simulations, we here show that diffusion-evoked potential shifts can introduce errors in standard CSD analysis, and can lead to prediction of spurious current sources. Further, we here show that the diffusion-evoked prediction errors can be removed by using an improved CSD estimator which accounts for concentration-dependent effects.

NEW & NOTEWORTHY Standard CSD analysis does not account for ionic diffusion. Using biophysically realistic computer simulations, we show that unaccounted-for diffusive currents can lead to the prediction of spurious current sources. This finding may be of strong interest for in vivo electrophysiologists doing extracellular recordings in general, and CSD analysis in particular.