A generalized theory for current-source density analysis in brain tissue

A generalized theory for current-source density analysis in brain tissue

Claude Bédard and Alain Destexhe111Corresponding author. Tel: +33 1 69 82 34 35, Fax: +33 1 69 82 34 35, Email: Destexhe@unic.cnrs-gif.fr Unité de Neurosciences, Information et Complexité (UNIC),
1 Avenue de la Terrasse (Bat 33),
91198 Gif-sur-Yvette, France
Physical Review E, in press, 2011.
July 15, 2019

The current-source density (CSD) analysis is a widely used method in brain electrophysiology, but this method rests on a series of assumptions, namely that the surrounding extracellular medium is resistive and uniform, and in some versions of the theory, that the current sources are exclusively made by dipoles. Because of these assumptions, this standard model does not correctly describe the contributions of monopolar sources or of non-resistive aspects of the extracellular medium. We propose here a general framework to model electric fields and potentials resulting from current source densities, without relying on the above assumptions. We develop a mean-field formalism which is a generalization of the standard model, and which can directly incorporate non-resistive (non-ohmic) properties of the extracellular medium, such as ionic diffusion effects. This formalism recovers the classic results of the standard model such as the CSD analysis, but in addition, we provide expressions to generalize the CSD approach to situations with non-resistive media and arbitrarily complex multipolar configurations of current sources. We found that the power spectrum of the signal contains the signature of the nature of current sources and extracellular medium, which provides a direct way to estimate those properties from experimental data, and in particular, estimate the possible contribution of electric monopoles.

87.19.La, 87.17.Aa

I Introduction

The current-source density (CSD) analysis Mitzdorf (1985); Pettersen (2006) is a method consisting of estimating the underlying current sources from a series of recordings of the extracellular electric potential. This method is widely used in neuroscience, and applies well to layered structures of the brain, such as cerebral cortex, hippocampus or cerebellum Mitzdorf (1985). The CSD analysis is based on the “standard” model of electric potentials in biological tissue Mitzdorf (1985); Pettersen (2006); Plonsey (2007), which rests on the hypothesis that the extracellular medium is resistive (ohmic) and uniform. Other influences, such as ionic diffusion, are assumed to play a negligible role on the propagation of the electric field.

Based on this set of hypotheses, the equation that determines the electric potential at macroscopic scales ( or more) is given by:


where is the electric conductivity of the extracellular medium. This expression can be obtained by applying the differential law of charge conservation and Ohm’s law. The term is interpreted as the volumic density of current sources. This equation forms the basis of the CSD analysis method Mitzdorf (1985); Pitts (1952); Pettersen (2006); Nicholson (1975).

According to Eq 1, the electric potential would only depend on electric conductivity and not at all on electric permittivity. However, Poisson’s law in a homogeneous medium implies that V will be twice smaller for twice larger with the same charge distribution, so it is paradoxical that permittivity is not taken into account in CSD analysis. Moreover, according to Eq 1, the electric potential is determined solely by the charge conservation law, and independently of Poisson’s law, which is contradictory with Gauss’ law in Maxwell equations.

If we take Gauss’ law into account, we can write

The general expression for is:

where is the Maxwell-Wagner time of the medium. However, is usually considered as negligible (typical values for biological tissue are , with and ) such that the current source density must be approximately zero, which is paradoxical. One way to resolve this paradox is to consider that the Maxwell-Wagner time is not negligible, or that electric parameters display strong spatial variations. However, such conditions contradict the hypothesis that the medium is resistive, and lay outside the domain of validity of Eq. 1 because in this case the impedance of the extracellular medium is complex in Fourier space Bédard et al.  (2004); Bédard and Destexhe (2009).

Despite this paradox, the standard model seems to apply relatively well to media such as brain tissue Mitzdorf (1985); Pettersen (2006). This model has, however, the drawback that it cannot be used to determine the validity of the hypotheses it is based on. Moreover, there is no clear definition of microscopic or macroscopic levels, and consequently, it is difficult to include possible frequency dependencies that could result from different physical phenomena at intermediate (mesoscopic) scales, such as ionic diffusion or membrane polarization Bédard and Destexhe (2009).

In the present paper, we introduce a more general formalism which does not rest on the hypotheses of classic CSD where the medium is hypothesized to be uniform and resistive, which also supposes that the electric parameters are constant and independent of frequency. The goal of this new formalism is to provide generalized expressions for CSD analysis in non-resistive media. Our aim is also to provide a theory which is general enough to enable testing different hypotheses concerning the nature of current sources and the electrical properties of the surrounding extracellular medium, which could then be directly estimated from experimental data.

Ii General Theory

In this section, we derive a mean-field theory of the electric field and potential resulting from current densities in biological tissue, by staying as general as possible.

ii.1 Definitions and scales

In the generalized formalism presented below, we will define the current sources from conductance variations. We will assume that the differential law of charge conservation holds in a given domain D, without defining a current source density per unit volume. We assume that the conductance variations in cellular membranes (especially around synapses), within domain D, are the principal origin of the extracellular electric field. This assumption is more realistic and biological compared to the “classic” assumption which is based on current source densities per unit volume because real current sources are caused by the opening or closing of membrane conductances in neurons. Note that these two different points of view can be complementary if we assume that the volumic density of current source is

In this case, the two points of view are mathematically equivalent. This is the reason why we wrote the source term as in Eq. 1.

An important assumption of the present formalism is that all observable phenomena can be modeled by fields which are twice differentiable (class ). While most fields will obey this criterion, it will considerably simplify the mathematical analysis to its simplest expression (commuting spatial and temporal first-order derivatives). In mean-field physics, by virtue of the Stone-Weierstrass theorem, it is always possible to make a uniform-convergence approximation of the observable phenomena by a mean-field model of class . Indeed, because the mean-field of a discontinuous field of first kind is necessarily a continuous field (the primitive of a discontinuous function of first kind is continuous), this restriction to class fields will not limit the applications of the theory developed here. Moreover, the theory will remain general because most fields are class in practice. It will not apply to very particular models, such as fields involving surfaces with infinitely small thickness, current sources without volume, fields that necessitate infinite energies222In such particular cases, the partial and first-order derivatives of fields , , and are not defined for every point of space and time.

Finally, the formalism developed below is only valid for well-defined ranges of spatial and temporal scales. We will consider scales greater than (about 300 times the size of ions such as or , including solvatation). This scale is chosen large enough for classic electromagnetism theory to apply without ambiguity (although in principle it can apply to scales down to a few nm). We will define as “microscopic”, scales of the order of 1 m, while “macroscopic” scales will be of the order of 50 m or more. We will also consider the typical range of frequencies of electrophysiological signals, up to , for which the quasistatic approximation is valid.

ii.2 Mean-field Maxwell theory

We start from Maxwell equations where we consider the spatial averages of the fields and electric parameters, which will be denoted here by brackets . The spatial average is made over some reference volume, which is invariant. Because of the scale invariance of Maxwell equations (e.g., see Jackson et al.  (1996), chapter 4), the spatial averages of electric field , electric displacement , magnetic induction and magnetic field are linked by the following linear operatorial equations:


where and are the spatial averages of the current density and free charge density, respectively.

These equations allow one to find the general regularities that all models must satisfy. For example, the laws of energy conservation and momentum conservation can be deduced from this set of equations Jackson et al.  (1996). In particular, by using the relation (which is in general true for all vectorial fields of class ), one can deduce the differential law of charge conservation:


However, because the above equations relate the spatial averages of interaction fields (,) with the spatial averages of the two matter fields (), it is necessary to complete them with a specific physical model to apply them to a given biological medium. This specific model must allow measuring spatial averages at a scale which is determined by the measurement method (type of electrode for example). Thus, the measurement system determines a minimal reference volume, which necessarily implies to use a mean-field formalism.

In general, (for all media of class ), the fields , , and are linked by the following general equations:


where and are respectively the absolute magnetic permeability and absolute electric permittivity tensors. Here, we have defined where is the source field resulting from conductance variations. Note that in classic electromagnetism, one defines the electric displacement relative to vacuum permittivity by (in frequency space)333 is the effective electric field and the polarization field is produced by polarization of molecules and cell surface polarization. In general, the relation between these vectors is algebraic in Fourier space, and thus a convolution integral in temporal space., which expresses the fact that the polarization field is proportional to the electric field through electric susceptibility (). It follows that .

Considering Maxwell-Gauss’ law () and the definitions of the interaction fields imply the following relations between charge density and interaction fields:


where represents the average variation of free charge density produced by conductance variations, and is the average free charge density which does not result from membrane conductance variations. Note that the divergence of the field depends on the exact mechanism of conductance variation. If this mechanism does not produce monopoles, then this divergence is zero. Also note that the field is generally assumed to be independent of the field , which is a valid assumption for biological media in general, except if ephaptic interactions must be taken into account.

We also have


where is the electric conductivity and is the mean ionic diffusion tensor 444Because the law of ionic diffusion is given by when the units of are and when we have only one type of ion, we have multiplied the expression of by to yield the electric current density associated to each ionic species. z is the valence of the ions considered, and is the Faraday constant. The choice of the sign is according to the standard convention. Note that if the fundamental charge is taken as that of the proton, then one must multiply by the factor , but if it is that of the electron, then the multiplying factor is .. We define as follows:

where the sums run over the different ionic species. Thus, we can write that the part of current density caused by concentration changes equals 555Note that the spatial average of will have similar values for different ionic species because the diffusion coefficients of the main ions () have similar values for biological tissues in physiological conditions (see Nicholson (1998)). In expression (6), we have separated the current produced by ionic diffusion from the current produced by other physical causes such as Ohm’s effect, polarization, etc. Note that this separation was made here for simplicity, but it is also possible to integrate diffusion effects in the expression of the mean conductivity (see Eq. 42 in Section A.1).

It is important to note that the first term in the righthand side of Eq. 6 is not exclusively due to Ohm’s law (which relates to energy dissipation), but can reduce to it in some cases Bédard and Destexhe (2009); Bédard et al.  (2006b). In general, Eq. 6 gives a time-dependent electric conductivity (or frequency-dependent in frequency space), which is not the case for Ohm’s law in general (see Appendix A.3). Also note that the integrals in Eqs. 4 and 6 can be seen as convolution products relative to time, in which case they take the form of a simple product in frequency space.

Thus, according to this theory, it is sufficient to model the physical and geometrical nature of the extracellular medium by using electromagnetic parameters and diffusion coefficients to simulate the interaction fields when the current sources are known – this is usually called the forward problem. Inversely, we can also deduce the physical characteristics of the sources from the knowledge of the electromagnetic parameters and interaction fields, as well as their spatial and temporal variations – this is known as the inverse problem.

Finally, the integrals in Eqs. 4 and 6 must satisfy the causality principle, according to which the future cannot determine the present state of the system. For example, the value of electric field at time must not influence the value of electric displacement at time t. Thus, the causality principle determines a supplementary constraint on the possible values of tensors , and , which limits the number of possible mathematical models of the extracellular medium. For instance, as detailed below in Section II.6, this principle imposes mathematical relations between the electric parameters and , which are called Kramers-Kronig relations for linear media.

The set of equations above define a mean-field formalism in which Maxwell equations are formulated with spatial averages. In the next sections, we consider different approximations to this formalism.

ii.3 The quasi-static approximation in mean field

The first approximation to the Maxwell equations is the quasi-static approximation, which consists of de-coupling electric and magnetic variables. In general, the time variation of produces an electric field (Lenz-Faraday effect), the electric and magnetic variables are coupled in Maxwell equations (Eq. 2) by the following expression:


It was shown that for biological media and current sources, the Lenz-Faraday effect is negligible Bosetti et al (2008). In such conditions, we can write:


Under this quasi-static approximation, the electric and magnetic variables are de-coupled in Maxwell equations, and the electric field obeys:

This approximation is also called the a priori quasistatic approximation, by opposition to the a posteriori quasistatic approximation, which consists of finding the general solution of Maxwell equations and later de-couple the electric and magnetic variables (see details in Bosetti et al (2008)). Although these two approximations are not strictly equivalent, we will only consider the a priori quasistatic approximation in the remaining of this paper.

According to this approximation, Maxwell equations simplify to the following expressions:


where the current density is linked to the electric field by Eq. 6. Note that, contrary to the static cases (electrostatics and magnetostatics), the fields , , and are here space and time-dependent.

In the following, we consider the complex Fourier transform

where . Note that because of the linearity of the spatial average, we have .

Applying the complex Fourier transform to Eqs.4, 6 and 9 leads to:


The first of these equations is Poisson’s law in mean-field. Although in some cases (electrostatics), the Poisson equation is sufficient to determine the solution of the system, it is not sufficient in the quasistatic case, and the second equation is necessary to close the system. This second equation is the differential law of charge conservation (in the presence of diffusion) and takes into account the time variations of the electric potential

By multiplying the first equation of Eqs. 10 by , and adding the result to the second equation, leads to:


where is the admittance of the extracellular medium. This equation is general and can be used to calculate the extracellular potential in an extracellular medium with arbitrarily complex properties (i.e., when the electric parameters depend on frequency and space). It is a generalization of expressions obtained previously Bédard et al.  (2004); Bédard and Destexhe (2009). The righthand term accounts for ionic diffusion.

Because Maxwell equations are scale invariant, the expression above (Eq. 11) is valid at all scales. Like in any mean-field approach, the spatial scale can be chosen according to the scale of the phenomenon that needs to be modeled, as well as the physical size and distance between electrodes. For example, in the case of CSD of mammalian cerebral cortex, one must consider scales of the order of 50 m to resolve the field produced by each cortical layer.

Note that in the quasistatic approximation, the explicit dependence of the electric field on magnetic permeability completely disappears. However, there can still be an implicit dependence through in nonlinear media, because the electric field does not depend explicitly on magnetic induction anymore.

ii.4 The quasistatic approximation at larger scales

At small scales (m), biological media such as the cerebral cortex are far from homogeneous and isotropic. The electric parameters can display large variations, for instance between fluids and membranes. However, at larger scales (m), such media can be considered as homogeneous and isotropic. In such a case, the tensors and can reduce to scalar quantities. Note that the fact of considering larger scales suppresses the directional dependence of the propagation of currents by a statistical equivalent, without changing the frequency dependence produced by physical phenomena at small scales. The transition from small scales to larger scales gives the same form as Eqs. 11, but with scalar parameters which will have an explicit dependence on space, frequency and the values of the field in general666Note that the space dependence is much smaller than the large variations seen at microscopic scales, for example between fluids and membranes.. The rate of spatial variation of these parameters at scales of the order of m is approximately zero, such that:


Note that there can be a frequency dependence of the current propagation which results from microscopic inhomogeneities of the electric parameters Bédard et al.  (2004), from polarization phenomena Bédard et al.  (2006b) or from ionic diffusion Bédard and Destexhe (2009). This frequency dependence of the current will not disappear when considering larger scales. On the other hand, new frequency dependencies may appear, such as for example the transformation of a frequency-independent conductivity tensor at small scales (m) to a scalar conductivity at large scales (m) will be associated to a frequency dependence of this macroscopic conductivity (for details, see Bédard and Destexhe (2009)). In agreement with this, measurements of the macroscopic conductivity demonstrated strong frequency dependence in different biological tissues Gabriel et al.  (1996).

ii.5 The linear approximation in mean-field

Still within the quasistatic approximation, we now consider the further simplification that the extracellular medium is linear. In a linear medium, the electric parameters are independent of the values of the fields (note that this linearity is different than that of Maxwell equations, which are always linear). In this case, the electric parameters only depend on space and time (or space and frequency).

In this case, the system of equations (10) becomes linear at small scales. This linear approximation is easy to justify for the magnetic field, given the small currents involved (for example, in neocortex) and the gradient of is almost zero (see Hamalainen93 ()). In contrast, the linear approximation is less trivial in the case of the electric field because of the many nonlinearities involved. For example, several ionic conductances are strongly voltage-dependent (such as the Na/K conductances involved in action potentials), which will make the electric parameters of membranes strongly dependent on the electric field. Nevertheless, the total volume of tissue occupied by membranes is small compared to other regions where the linear approximation is valid, so biological tissues can in general be considered as linear. Note that this linearity is evident for low frequencies (10 Hz), but it is less evident for high frequencies (100 Hz), where nonlinear phenomena such as action potentials can have a major contribution.

ii.6 The Kramers-Kronig relations under the linear approximation

As discussed above in Section II.2, the causality principle determines a supplementary constraint on the possible values of tensors , and . In the linear approximation, one can show that, in general (for isotropic media of class ), the linking equation between the electric displacement and electric field takes the following form:


In this case, one can show that the frequency dependence of electric parameters is not arbitrary but is linked by the Kramers-Kronig relations (see Section 82 in Landau1981 ()):


where principal integrals () are used. is the absolute electric permittivity of vacuum and is the static electric conductivity (). Note that these relations can be seen as a direct and inverse transform. The Maxwell-Wagner time represents the characteristic time (or “inertia”) for settling into a stationary regime, and can be strongly frequency dependent Bédard et al.  (2006b). Interestingly, this ratio is mathematically analogous to the time-frequency uncertainty principle in Fourier transforms when the electric conductivity of the extracellular medium becomes very small at zero frequency (see Appendix A.5). Note that there exists no such relation for the spatial variations of electric parameters, which are specific to each medium.

It is important to note that the Kramers-Kronig relations have a strong consequence on the plausibility of purely resistive media. If a medium is purely resistive, then both conductivity and permittivity are constant and independent of frequency. However, if one takes into account a very weak frequency dependence of conductivity (“quasi-resistive” media), such as with , then the Kramers-Kronig relations impose that we necessarily have . In such a case, the permittivity will be strongly frequency dependent, so will be the Maxwell-Wagner time . Thus, a purely resistive extracellular medium is a singularity, and is not likely to be realistic model for complex biological media.

Iii Application of the quasistatic mean-field theory to linear media

In this section, we apply the theory outlined above to media which are linear in the electromagnetic sense, and which are homogeneous and isotropic (at macroscopic scales). This is equivalent to assume that the macroscopic parameters , and are scalars independent of space, such that and . This approximation is certainly valid for relatively large distances (greater than ). We also consider the system under the quasistatic approximation as defined above.

In the following sections, we examine different limit cases. The first case corresponds to the standard model with dipolar sources and a resistive (or quasi-resistive) extracellular medium. A second case will consider the same model, but with additional monopolar sources. The third case will consider ionic diffusion (not present in the two first cases), which also implies monopolar sources. In each case, we will derive the expression to use for CSD analysis.

iii.1 Dipole sources in resistive and quasi-resistive media

We start with the standard model in which the electric conductivity of the extracellular medium is constant (in space and frequency), scale invariant, and isotropic. We also consider that diffusion is negligible. Under these hypotheses, we have (homogeneous and isotropic medium) and , and the general formalism (Eq 11) reduces to:


We can also consider a slightly more realistic model of the extracellular medium by assuming that it is quasi-resistive () instead of resistive () (because the latter represents a singularity as outlined above in Section II.6), isotropic and homogeneous for large scales (m; the medium is allowed to be non-homogeneous for smaller scales of the order of m). Note that in general, the non-homogeneity of conductivity at smaller scales can induce a frequency dependence at larger scales (see Bédard et al.  (2004); Bédard and Destexhe (2009)). Thus, the hypothesis that the macroscopic conductivity is independent of frequency is equivalent to assume that there is no significant variations of impedance at microscopic scales.

In such conditions, Eq. 11 again reduces to:


In some formulations, the standard model does not consider the possibility of microscopic (m) monopolar sources Mitzdorf (1985); Pettersen (2006). This is equivalent to hypothesize that, at every time, each portion of cell membrane has an equal number of positive and negative charges at opposite sides of the membrane, such that it is locally neutral. This hypothesis is also equivalent to state that the extracellular electric field is produced by dipoles (or more complex multipolar arrangements), and that the monopolar component of the field is negligible at scales of m777Note that one cannot say that the monopolar component is rigorously zero, because there is at least a monopolar component in the ion channels themselves, because of ion selectivity.. This implies that the attenuation of the extracellular potential follows a law which varies as (or , … for multipoles of higher order) when , where is the distance to the source. Thus, in the standard model, the electric displacement in frequency space (see Eqs. 4 and 5) is given by:




at large scales (m). Taking the inverse Fourier transform, one obtains:


According to this model, the inverse solution (CSD method) can be obtained assuming that the voltages measured at different extracellular sites are solution of Laplace equation. According to the superposition principle, the extracellular potential can be considered as resulting from a sum of macroscopic dipolar sources for sufficiently large . Note that the value of is determined by Shannon’s sampling theorem, according to which the number of samples (number of electrodes ) must be twice larger than the larger spatial frequency of the field. To evaluate these dipolar sources, one can simply apply the inversion of the matrix linking the measured voltages with the dipolar sources according to the “forward” solution of Laplace equation (see e.g., Pettersen (2006)). Note that this approach is different than the classic CSD method proposed by Mitzdorf Mitzdorf (1985), which is based on a Poisson type equation.

The hypothesis of local neutrality in a homogeneous and isotropic extracellular medium implies that the frequency dependence of the measured signal is only due to the frequency dependence of the source (for example the effect of morphology – see Linden2010 (), the exponential or bi-exponential nature of synaptic conductances, correlations in synaptic activity, action potentials, etc), because Laplace equation does not explicitly depend on frequency. Thus, according to the standard model, there is no filtering due to extracellular space and the power spectrum of the extracellular potential is identical to that of the current sources.

Finally, it is important to note that, in a resistive extracellular medium, if we express the extracellular potential as a function of the dipole moments instead of the current sources, then the power spectral density (PSD) of the electric potential will necessarily have a supplementary frequency dependence of the form compared to that of the current. This is due to the fact that the current is proportional to the temporal derivative of the the dipole moment (defined from the charge distribution). However, the situation is different if the medium is quasi-resistive. In this case, the Kramers-Kronig relations give , and thus the ratio will have little frequency dependence. It follows that the PSD of the extracellular potential will have approximately the same frequency dependence as the current sources in a quasi-resistive medium. This is a striking difference between resistive and quasi-resistive media. As discussed above, the latter is a more realistic situation because any spatial variation of microscopic conductivity will necessarily induce a frequency dependence of the macroscopic conductivity (see Bédard et al.  (2004)).

iii.2 Monopolar sources in resistive and quasiresistive media

In the previous section, we hypothesized that the extracellular medium is locally neutral at microscopic scales (m), and thus, that the sources of the electric potential are dipoles. We now relax this hypothesis, and allow significant electric monopoles to appear in addition to conductance variations, so that the field results from both monopolar and dipolar contributions. Electric monopoles could result from different physical sources, such as the ionic selectivity of synaptic ion channels (similar to a “Maxwell Daemon”), combined with the finite velocity of charge movement Bédard and Destexhe (2008). These factors should create some accumulation of charge in the vicinity of the synapse when synaptic conductances are activated888Note that it is important here to take into account the spatial extent of the synapse, contrary to the standard theory where synapses are considered as point processes.. Note that monopoles are transient by definition, and equivalently, one could consider that the conductance variations determine a non-stationary regime (see Appendix A.2). In this transient regime, Kirchhoff’s “point rule” does not apply (it is based on the law of current conservation ) and would apply only when the system reaches a stationary state. However, Kirchhoff’s “loop rule” is always valid under the quasistatic approximation, because we have , and consequently , which is at the basis of the latter rule.

Contrary to the assumptions of the dipole model, monopolar sources imply that integrating the charge density over a closed surface surrounding each source is non-zero. To include the contribution of monopolar current sources, we have


where the total charge contained in the source. Note that this relation shows that the monopolar component is linked to the current through a temporal derivative, which is a consequence of the charge conservation law. Consequently, the extracellular potential (which is here proportional to the charge) will not have the same power spectrum as the source, and will have an additional component for a resistive extracellular medium. However, similarly to the case of dipolar sources in previous section, the situation is different for a quasi-resistive medium. The Kramers-Kronig relations imply , and the ratio will have very little frequency dependence and the PSD of the extracellular potential will be very similar to that of the sources.

If we consider the same conditions as for the standard model (resistive or quasi-resistive media), we obtain


In such conditions, Eq. 11 becomes:


for resistive and quasi-resistive extracellular media. Thus, in temporal space, we have the same equation for both cases:


However, if we take into account monopolar current sources and the law of charge conservation, then we have in general:


where .

Thus, the model with monopolar current sources has a different structure than the dipole model in Section III.1 because . Local neutrality in a homogeneous and isotropic extracellular medium implies an identical frequency dependence of the current source ) and the extracellular potential. Like the standard model, there is no “filter” in this case. There is a notable difference, however. The law of attenuation with distance varies here in instead of for . If the number of electrodes is large enough to respect Shannon’s sampling theorem, then the current source densities can be simply evaluated by approximating the Laplace equation using finite difference methods, as well as the knowledge of the “forward” solutions of this equation (see Pettersen (2006)). We will see in the next section that these conclusions are different if ionic diffusion is taken into account.

iii.3 Models with ionic diffusion

While the influence of ionic diffusion was neglected in the previous sections, we now consider this case more explicitly without any other hypothesis about the medium. If a selective ion channel opens, the flow of ions may induce accumulation of charges in the region adjacent to the channel if ions diffuse faster than the time needed for passing through the channel (which will generally be the case). The electric field resulting from conductance variations is not selective on the type of ion, such that the positive ions are attracted and negative ions are repulsed if the field is negative (and vice-versa for a positive field). This is contrary to the flow inside the channel because it is selective to only a subset of ionic species. The combination of these factors makes it unavoidable that there will be charge accumulation around open ion channels. In the standard model, this charge accumulation is considered as negligible.

We now evaluate the consequences of this phenomenon on the frequency dependence of the field produced by ionic conductances in the subthreshold regime. If we consider a homogeneous extracellular medium with constant electric parameters (independent of frequency at large scales, m), then we have:


where the parameters , and do not depend on space.

According to those hypotheses, variations of ionic concentrations appear in the vicinity of the open ion channels, and these variations are opposite to the current produced by the electric field resulting from conductance variations. It thus appears that the conditions of current propagation at microscopic scales (m) cannot fulfill the condition of homogeneous ion concentration which is at the basis of Ohm’s differential law (see Appendix A.3). In such conditions, the electric parameters of the extracellular medium have the following form at microscopic scales (m):




This last equation is necessary to keep the consistency between microscopic (m) and macroscopic (m) scales.

According to this model, the current density (at microscopic scales, m) is given by:


where we have (see Eq. 49)


This expression can be deduced by separating the domain into sufficiently small elements such that ion density can be considered as spatially homogeneous, and sufficiently large for Ohm’s law to apply.

According to Boltzmann distribution (see Appendix A.1), we have


By taking into account Eqs. 28 and 29, we obtain


where is an effective diffusion coefficient. Note that the value of is smaller than the mean diffusion coefficient because must be positive. The value of also depends on the values of ionic concentrations because several parameters in Eq. 31 are concentration-dependent and is proportional to temperature because the ionic diffusion coefficient is itself proportional to temperature (see for example the Einstein relation for diffusion).

Applying the differential law of charge conservation, we get


Thus, the charge density produced in the vicinity of the ion channel is solution of a parabolic differential equation similar to the diffusion equation999 Note that the coefficient depends on ion concentrations via , and , and thus could vary greatly according to the activity of the surrounding neurons..

It follows that the charge density obeys:


At microscopic scales (m), we obtain (see Eq. 11):


Here, the proportionality between and can be deduced from the linear (first-order) approximation of Eq. 41 (see Appendix A.1). The second-order approximation would give a cubic term in .

Applying the consistency equation between scales by assuming the statistical independence of the parameters leads to the following equality:




where is the ratio between the reference volumes at macroscopic and microscopic scales (note that to simplify the formalism, we have approximated the macroscopic mean by a discrete summation over microscopic means). The second term of the lefthand side of Eq. 34 becomes zero at macroscopic scales (see consistency equation Eq. 27). Note that the means over parameters and are harmonic means, while the means over matter fields are arithmetic means.

Finally, by applying the inverse Fourier transform, we obtain (for Maxwell-Wagner times much smaller than unity):


Thus, the CSD method in the presence of ionic diffusion takes a form which is very close to the Mitzdorf model Mitzdorf (1985), because we have one source term. However, there are two notable differences: first, the frequency dependence of charge density implies that the extracellular medium will be frequency dependent according to an impedance which varies as (see Appendix A.4). Second, the extracellular potential attenuates with distance according to a Yukawa potential instead of , as in the standard model. In this case, we have


Here, the extracellular potential is proportional to the charge density (see Eqs. 41, 52, and 53 in Appendix A.4) under the linear approximation and for a spherical source. It is interesting to note that the exponential term increases with frequency such that the extracellular medium favors the propagation of low frequencies (low-pass filter), as shown in Fig. 1). This type of attenuation law in Fourier space is also consistent with an exponentially-decaying impedance. If the frequency spectrum is narrow, one can replace by its maximal value , which leads to an attenuation law for the potential as . We can thus write that the electric field is approximately equal to . In a resistive medium, this leads to an electric resistivity given by (note that the current is conserved in a resistive medium, and the field is given by ). Note that this particular distance profile of the potential was calculated by assuming that the medium is homogeneous (see Appendix A.4), which makes it applicable only at short distances from the membrane (of the order of 10 to 50 nm).

[ ————– Figure 1 here ————– ]

Finally, it is important to note that this frequency dependence cannot be removed because it is an effect of the feedback caused by ionic diffusion when ion channels open, and this is inherent to biological tissue. Because the PSD of the extracellular voltage is of the form (see Eq. 53 in Appendix A.4), one can view the effect of ionic diffusion as a “1/f filter”, as found previously Bédard and Destexhe (2009).

iii.4 Comparison between the different models

We now compare the different cases examined here. From the point of view of the differential equations involved, in the “standard” model based on dipoles, as well as with monopoles, the extracellular potential is solution of the Laplace equation, which is elliptic. In the third model with ionic diffusion, the extracellular potential is solution of a Poisson type equation where the source term is proportional to the time derivative of the voltage (under the linear approximation), which gives a parabolic equation. As outlined above, the diffusion model is closer to the monopole model (as diffusion can have monopolar effects) in a resistive extracellular medium, but leads to a fundamentally different mathematical form. The physical reason for this difference is that the ionic diffusion at the interface ion channel/medium increases the inertia of the system as a function of frequency.

At the point of view of the CSD analysis method, different algorithms must be used according to which model of the extracellular medium is assumed. In the two first cases, one must use a “forward” solution because Laplace equation is non invertible. In this case, it is necessary to explicitly include the distance dependence of the extracellular potential , which varies as for dipoles (for distances sufficiently large compared to the size of the dipole) and for the model based on monopoles. Note that if the distance to the sources is not large enough (compared to the typical size of the sources), or if the dipolar moments are very large compared to monopolar moments, then the attenuation will be closer to a linear combination of and .

In the diffusive model, however, the approach is totally different because of the parabolic nature of the equations. In this case, it is enough to apply the Laplace operator to recover the sources. Two strategies are possible. First, one could simply apply Laplace operator on the extracellular voltage to yield estimates of the current source densities. Second, one could use a “forward” model and consider an attenuation law following a Yukawa potential and apply the same procedure as for the other models.

Perhaps the most interesting aspect is that the three different models considered here have a different spectral signature. In the dipole model (or monopole model in a quasi-resistive medium), the PSD of the extracellular potential is identical to that of the sources. The resistive monopole and the resistive dipole model exert a filtering effect of type, whereas the diffusive model is equivalent to a filter. Thus, the frequency characteristics of the signal can serve as a criterion to determine the most appropriate model. For example, if the PSD of the extracellular voltage has structure, this automatically discards a pure monopole model, as well as dipole models in resistive or quasiresistive media, and would suggest diffusive type models.

iii.5 Synthesis and applications to experimental data

In this section, we synthesize the theoretical developments provided here, and suggest a guide of how to apply them to experimental data. The generalization of the CSD method for different cases of current sources and type of extracellular medium, is summarized in Table 1. The table considers monopolar and dipolar current sources, as well as different types of resistive and non-resistive media.

To perform a CSD analysis by allowing non-resistive properties of the extracellular medium, we suggest the following procedure.

  1. Estimate the type of extracellular medium from the power spectral structure of CSD signals. As detailed above (Sections III.1, III.2, III.3, the type of medium (resistive, quasi-resistive, diffusive, etc) and type of current sources (monopolar, dipolar, etc) can be inferred from the power spectral structure of the extracellular potential. This analysis should be done on non-filtered data to set constraints on the possible combinations of sources type of medium.

  2. Identify the correct CSD expression compatible with the type of source/medium inferred from power spectra. Table 1 summarizes the different cases considered here. The expression identified is then used to calculate the current sources from the extracellular potential recordings.

Source Medium Law
1-pole res.
gen. gen.
2-pole res.
gen. gen.
Mitzdorf (2-pole) res.
Table 1: Different generalizations of the CSD method. The table shows the mean-field equations for different types of media, and for monopolar or dipolar sources. The Mitzdorf model is shown apart, because it does not correspond to any of these mean-field scenarios. Abbreviations: res. resistive homogeneous medium, quasi-res. quasi-resistive homogeneous medium, res. + dif. resistive homogeneous medium + ionic diffusion, gen. general (non-homogeneous, with spatial and frequency-dependent variations of electric parameters), 1-pole monopole, 2-pole dipole. Note that the frequency-dependence of the permittivity and conductivity are not independent but are linked by the Kramers-Kronig relations. The quantity is the ratio between the Fourier transform of the extracellular potential and the Fourier transform of each point current-source which produce the field (asymptotic solution, far from the sources). The function in “res+dif” determines a Yukawa type potential (see Fig. 1).

Iv Discussion

In this paper, we have formulated a series of generalizations of the CSD analysis method applicable to extracellular recordings in brain tissue. This generalization is based on a general theory that we derived and which aims at linking the extracellular potential with current source densities in the tissue. We have considered a mean-field version of Maxwell equations by considering the different fields as averages over some reference volume. By varying the size of this volume, one can apply the same theory to different scales. At microscopic scales (m and smaller), the theory must use the microscopic values for electric parameters (for example, the very different resistivities of fluids or membranes). For mesoscopic or macroscopic scales (m and larger), the theory can directly include the “macroscopic” measurements of conductivity and permittivity, as well as their possible frequency dependence if needed. Note that this mean-field approach takes into account the physical and biological properties of the sources, and thus is more general than previous approaches Bédard and Destexhe (2009) which did not consider source densities.

We have examined different limit cases, such as a purely resistive extracellular medium with current sources consisting exclusively of dipoles, in which case the theory recovers the standard model. In this standard model, the mean-field theory shows that the electric potential must be solution of Laplace equation, such that the “classic” CSD approach of Mitzdorf Mitzdorf (1985) does not apply. To inverse the CSD in this model, one must apply the forward solutions of Laplace equation because the associated operator is non-invertible (see Pettersen (2006)). In resistive media, the extracellular potential must have an additional frequency dependence of relatively to that of the current. Interestingly, we found that Laplace equation remains valid for extracellular media which are quasi-resistive (where the electric parameters weakly depend on frequency). In this case, the frequency dependence of the extracellular potential is similar to that of the current. A weak frequency dependence was indeed found in some experimental measurements of resistivity Logo2007 (); Ranck (1963), while other experiments Gabriel et al.  (1996) displayed a much more pronounced frequency dependence. With respect to the attenuation with distance, the standard model predicts an attenuation law as , for both resistive or quasi-resistive media.

We also examined the case of monopolar sources. If such monopolar sources are present in addition to dipolar sources, within resistive or quasi-resistive media, then the CSD equation takes a slightly different form predicting that the potential will attenuate asymptotically with the inverse of distance (), while the standard dipole model predicts a square dependence (). With monopolar sources, the potential in the extracellular medium is also solution of Laplace equation, and thus the inverse algorithm of the CSD method does not apply identically as for dipoles. To find the inverse CSD, one proceeds similarly as the standard model by using the “forward” solution of Laplace equation. However, in this case the sources must be considered as a linear combination of terms varying as (monopoles) and (dipoles) in this forward solution.

As a third model, we examined the case of ionic diffusion within resistive or quasi-resistive media. In this case, the CSD takes a form very close to the “monopolar” CSD discussed above, but we found that charge density is frequency dependent according to a Warburg impedance in (see Appendix A.4 and Section III.3). This result is in agreement with a previous modeling study of extracellular potentials in the presence of ionic diffusion Bédard and Destexhe (2009). Another consequence is that, for spherical symmetry, the attenuation with distance follows a Yukawa potential , which decays faster than the different laws considered above for large enough frequency (see Fig. 1). This particular form is responsible for a low-pass filtering of the extracellular medium. Note that this form is obtained for spherical symmetry, but other forms may be obtained in different geometries.

It is important to note that the CSD theory was originally designed without specific hypotheses about the nature of current sources Plonsey (2007); Nicholson (1975); Pitts (1952), other versions of the CSD theory clearly assumed that current sources are dipoles Mitzdorf (1985); Pettersen (2006). Assuming dipolar sources is equivalent to assume that we have stationary current conditions at all scales. However, we show here (Appendix III.2) that at small scales (synapses), such a stationary current condition is not necessarily met. A first possible source of monopolar effects is the inertia of charge movement along membranes together with ion-channel selectivity. Following the opening of ion channels, the flow of ions will entrain a re-equilibration of the charges adsorbed on both sides of the membrane. While this process is usually considered as instantaneous, together with neglecting ion-channel selectivity, these processes may have important consequences. Indeed, if one takes into account the fact that charges do not move instantaneously and ion-channel selectivity, this will necessarily create transient charge accumulation and monopoles. A similar effect will occur through ionic diffusion and electric field, at the interface between the ion channel and the extracellular medium, because ions diffuse faster than their mean passage time through the channel, which will also create charge accumulation and monopolar effects. Note that when this mechanism produces an external electric field which will contribute to the extracellular field (in addition to transmembrane currents). These effects will contribute to transient monopoles, during which Kirchhoff’s node law will not apply. The fact that ions move considerably slower than electrons in a metal conductor will also participate to deviations from Kirchhoff laws. Whether this transient time is significant, and whether the system could be continuously “outside of equilibrium” due to sustained synaptic activity, should be investigated by future work.

In the diffusion model, one can directly use the Laplace operator to inverse the CSD, contrary to the other models. Taking into account ionic diffusion requires to revise the “forward” approach, because the attenuation law does not follow a law, but rather a Yukawa-type law while the extracellular medium is associated to a Warburg-type impedance. In a previous study, we showed that indeed a Warburg type impedance could account for the transfer function between intracellular and extracellular potentials BedDes2010 () (for frequencies comprised between 3 and 300 Hz). It is also consistent with measurements of conductivity and permittivity Gabriel et al.  (1996) (but see Logo2007 ()). Note that the linear approximation in the diffusion model is not valid for high values of the potential (larger than 50 mV; see A.4), so this model applies well to subthreshold activity, but may need to be revised for action potentials. Similarly, corrections to the CSD given by the “forward” approach (see for example Leski2010 ()) may also need to be reformulated for non-resistive media.

Thus, with respect to the paradox of the CSD method, as described in the introduction, our study suggests that it is naturally solved by taking into account ionic diffusion. This introduces an additional source term in the general equation for the electric potential (see Eq. 11). This additional term gives a Poisson type equation for the potential (instead of Laplace equation), similar to the classic CSD approach. Contrary to the cases with resistive and quasi-resistive media, the classic algorithm of CSD inversion given by Mitzdorf Mitzdorf (1985) is applicable here. Thus, the results obtained with the classic CSD analysis are perfectly consistent with ionic diffusion because diffusion gives a source term which is very close to the phenomenological model of current source density introduced by Pitts and Mitzdorf Mitzdorf (1985); Pitts (1952), but in a manner consistent with Maxwell-Gauss law. So, we conclude that the usual approach for CSD inversion, although paradoxical, should nevertheless give results equivalent to a model with ionic diffusion and consistent with Maxwell-Gauss law.

Finally, the few limit cases considered here are by no means exhaustive. For example, we neglected the Maxwell-Wagner time of the extracellular medium and the microscopic variations of impedance. The theory outlined here is general enough to include these effects if needed, which is another way to solve the paradox. For instance, considering phenomena such as “reactive” extracellular media, which react to the electric field (for example through polarization of cell membranes), can be done by taking into account the Maxwell-Wagner time of the medium (see details in Bédard and Destexhe (2009); Bédard et al.  (2006b)). According to Gabriel et al. Gabriel et al.  (1996), the macroscopic electric permittivity becomes larger while macroscopic conductivity becomes smaller for smaller frequencies, when the electric field is imposed according to a well-defined direction. In these measurements, for frequencies larger than 10 Hz, but may be considerably larger for lower frequencies Bédard and Destexhe (2009), where electric polarization may play an important role. The second term in the lefthand side of Eq. 11 would then not be negligible anymore. Because this term can be considered as an additional source term (see Bédard et al.  (2004)), similar to the case of diffusion, this also solves the paradox described in the introduction.

In conclusion, we have provided here a generalized CSD approach valid for more realistic properties of the extracellular medium, taking into account ionic diffusion or polarization effects, usually neglected in the standard CSD analysis Mitzdorf (1985); Pettersen (2006). We found that including such effects may have deep consequences on the expression to be used for estimating current sources, and thus may also have consequences on the values of current sources estimated from experimental recordings. For example, the potential due to monopolar sources will decay slower than for dipoles, which will necessarily affect the recorded potential at the electrode. Similarly, considering “reactive” aspects of the extracellular medium by including a significant Maxwell-Wagner time leads to a different CSD expression, close to the form derived for ionic diffusion. Future work should apply these expressions to extracellular recordings in brain tissue, with the aim of identifying which of these phenomena are most consistent with experimental data.

Appendix A Appendix

a.1 Impedance for systems with ionic diffusion

In this appendix, we consider ionic diffusion at the interface between ion channels and the extracellular medium, as well as at the interface with the cytoplasm. We use the quasistatic approximation (in the thermodynamic sense), which implies that the net charge density must be solution of a parabolic partial differential equation, as for pure diffusion phenomena. We will next consider system in spherical symmetry, in which case the impedance is equivalent to a Warburg impedance.

a.1.1 Ionic diffusion under the quasistatic approximation in the thermodynamic sense

The current density at microscopic scales obeys the equation:


with (see Eq. 49)


Let us assume that the system is in a quasi-static case in the sense of thermodynamics. As shown by application of Maxwell distribution of velocity distribution and the principle of detailed balance Vas1983 (), we can deduce the Boltzmann distribution for a field which varies infinitely slow. This approximation is valid here because the drift velocity of ions under an electric field is much lower than the absolute velocity of ions (which is of the order of sound velocity). Within this quasistatic approximation, we can apply the Boltzmann distribution to obtain the number of ions per unit volume as a function of time and space:

when we assume that at infinite distance, where is the number of ions per unit volume at an infinite distance from the source (“far distance”) and is the mean absolute charge. is the Boltzmann constant and is the temperature in degrees Kelvin. It follows that the net charge density is related to the value of the electric potential according to:


where is the average net charge density. Note that the sign in the second righthand term comes from the sign of the charge. Also note that this relation implies that the net charge density is zero at an infinite distance, and is linked to the electric potential by a nonlinear relation101010A consequence of this nonlinear relation is that the medium will become nonlinear for high values of the electric potential..

Applying the operator on the net charge density gives:


Note that no such relation would be possible outside of the quasistatic approximation in the thermodynamic sense.

a.2 Non-stationary aspect of the electric field produced by membrane conductance variations

In this appendix, we show that if the transmembrane current is non-zero, this necessarily implies a charge variation at the interior of the compartment according to a non-stationary regime

We define as “interior” the domain delimited by the inner surface of the cell membrane as indicated in Fig. 2. We will show that for different values of stationa