Kramers-Kronig Relations and the Properties of Conductivity and Permittivity in Heterogeneous Media

Show more

1. Introduction

The theory of electromagnetism can be applied to complex media like inhomogeneous materials or biological tissue. A first approach to such media is to explicitly consider their microscopic structure, and the associated variations of electric conductivity or permittivity, but this approach requires a detailed mapping of these electric parameters and includes this in detailed simulations. Another approach is to use a mean-field electromagnetic theory by considering scales larger than the typical scales of inhomogeneities in the medium. In this case, the mean-field theory relates to macroscopic measurements of conductivity and permittivity. For example, in neural tissue, macroscopic measurements of these quantities were done in a number of studies [1] [2] [3] [4] who measured the frequency dependence of electric parameters in different conditions (reviewed in ref. [5] ). Theoretical work showed that this frequency dependence can be accounted by physical phenomena such as ionic diffusion or cell polarization [6] [7] . In the present paper, we investigate the relations between the apparent macroscopic conductivity and permittivity^{1} and show that they cannot take arbitrary values but are strongly constrained to be consistent with Maxwell’s theory of electromagnetism.

An important aspect of the equations of electromagnetism when applied to media with a given electric permittivity $\epsilon \left(r\mathrm{,}t\right)$ , is that the real ${\epsilon}^{\prime}$ and imaginary ${\epsilon}^{\u2033}$ parts of the complex Fourier transform of $\epsilon \left(r\mathrm{,}\omega \right)$ cannot take independent values, but they are linked together by a set of mathematical relations called Kramers-Kronig relations [8] [9] [10] . These relations are a direct consequence of the principle of causality, namely that the future cannot influence the past [8] , and were shown to be equivalent to the Hilbert transform in a perfect dielectric [11] [12] .

Here, we will re-examine these Kramers-Kronig relations and their equivalence to Hilbert transforms in heterogeneous media. The goal is to determine the constraints that the frequency dependence of electric conductivity and permittivity must satisfy. We will also consider the application of this formalism to a few concrete examples.

2. Theory

We begin by setting the framework of the present study, by outlining the electromagnetic theory used here, and how to apply it to heterogeneous media. We next consider the Kramers-Kronig relations within this framework.

2.1. General Framework

^{1}Here, “apparent” refers to the values that can be measured experimentally.

^{2}Note that there also exists a magnetic quasi-static approximation, where electromagnetic induction is not negligible, but the displacement current is rather neglected. The latter approximation applies when
$\Vert E\Vert \ll \mathrm{}c\Vert B\Vert $ .

To derive a formalism applicable to heterogeneous media, we consider linear media within the electric quasi-static approximation in mean-field^{2} [13] [14] , which corresponds to a physical situation where electromagnetic induction can be neglected. This is the case for neural tissue, where we have an excellent approximation of the electric field if we assume

$\{\begin{array}{l}\nabla \cdot D={\rho}^{\text{free}}\\ \nabla \times E=0\end{array}$ (1)

where the fields $E\mathrm{,}D\mathrm{,}{\rho}^{\text{free}}$ are mean-fields over a base volume $\mathcal{V}$ [7] . By definition, we have $X={\langle X\rangle |}_{\mathcal{V}}$ .

In this approximation of the Maxwell-Heaviside equations, the electric field and magnetic induction are such that $\Vert E\Vert \gg \mathrm{}c\Vert B\Vert $ where c is the velocity of electromagnetic waves [13] . Thus, in this approximation, we can calculate $E$ independently of $B$ and $H$ if we know the linking equation between $D$ and $E$ .

On the other hand, the link between the electric displacement field and magnetic field always hold because we have

$\{\begin{array}{l}\nabla \cdot B=0\\ \nabla \times B={\mu}_{o}\left[{j}^{\text{free}}+\frac{\partial D}{\partial t}\right]\end{array}$ (2)

where ${\mu}_{o}$ is the magnetic permeability of vacuum. We do not consider here the situation where this permeability would be different from vacuum, which normally should be a good approximation because of the absence of large amounts of ferromagnetic, paramagnetic or diamagnetic materials in neural tissue.

Moreover, if we restrict to isotropic media (when the inhomogeneity of the medium similarly affects all directions), we have:

$\{\begin{array}{l}D\left(r,t\right)={\displaystyle {\int}_{-\infty}^{+\infty}}\epsilon \left(r,\tau \right)E\left(r,t-t\right)\text{d}\tau \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{(a)}\\ {j}^{\text{free}}\left(r,t\right)={\displaystyle {\int}_{-\infty}^{+\infty}}{\sigma}_{e}\left(r,\tau \right)E\left(r,t-\tau \right)\text{d}\tau \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{(b)}\end{array}$ (3)

where $\epsilon \left(r\mathrm{,}t\right)$ and ${\sigma}_{e}\left(r,t\right)$ are real-valued functions.

This formalism applies to neural tissue, which can be considered as a heterogeneous isotropic medium within a mean-field context, when the base volume is sufficiently large (>1 mm^{3}; see ref. [7] ). The linking equations imply that the values of
$D$ and
${j}^{\text{free}}$ at a given time depend on the past values of the electric field in general, which can be seen as a kind of “memory”. The only way to avoid such a memory is to assume
${\sigma}_{e}\left(r,t\right)={\sigma}_{c}\left(r\right)\mathrm{}\delta \left(t\right)$ and
$\epsilon \left(r,t\right)={\epsilon}_{c}\left(r\right)\mathrm{}\delta \left(t\right)$ where
${\sigma}_{c}$ and
${\epsilon}_{c}$ are not time-dependent. Note that such a memory is equivalent to having a frequency dependence of the electric permittivity and/or conductivity, when we formulate the problem in Fourier frequency space. In this case, the relations (3) imply
$D\left(r\mathrm{,}\omega \right)=\epsilon \left(r\mathrm{,}\omega \right)E\left(r\mathrm{,}\omega \right)$ and
$j\left(r\mathrm{,}\omega \right)={\sigma}_{e}\left(r\mathrm{,}\omega \right)E\left(r\mathrm{,}\omega \right)$ .

Finally, recalling that $\nabla \cdot D$ allows us to calculate the free charge density in a region $\mathcal{D}$ , we have:

${Q}^{\text{free}}\left(t\right)={\displaystyle \underset{\mathcal{D}}{\iiint}}\text{\hspace{0.05em}}\nabla \cdot D\left(r,t\right)\text{d}x\text{d}y\text{d}z$ . (4)

The electric permittivity $\epsilon $ measures the amount of free charges in a given region. The higher the density, the larger the permittivity.

2.2. Electric Conductivity and Permittivity in a Mean-Field Model of Isotropic Media

In the following, we will re-examine the Kramers-Kronig relations in the case of a heterogeneous and conductive medium, where the electric field is time dependent. In a first step, we show that the Kramers-Kronig relations are equivalents to the Hilbert transform if we apply them to $\omega \epsilon \left(r\mathrm{,}\omega \right)$ instead of $\epsilon \left(r\mathrm{,}\omega \right)$ . In a second step, we show that we have the same relations for the conductivity ${\sigma}_{e}$ .

2.2.1. General Expression for the Absolute Electric Permittivity

According to relations (3), the electric displacement field $D$ is linked to the resulting electric field $E$ by a convolution integral:

$D\left(r,t\right)={\displaystyle {\int}_{-\infty}^{+\infty}}\epsilon \left(r,\tau \right)\mathrm{}E\left(x,t-\tau \right)\text{d}\tau ={\displaystyle {\int}_{-\infty}^{+\infty}}\epsilon \left(x,t-\tau \right)E\left(r,\tau \right)\text{d}\tau $ . (5)

However, we must also assume that $\epsilon \left(x\mathrm{,}t\right)=0$ when $\tau <0$ because the future cannot influence the past (causality principle). Thus, we can write

$D\left(r,t\right)={\displaystyle {\int}_{{0}^{-}}^{+\infty}}\epsilon \left(r,\tau \right)E\left(r,t-\tau \right)\text{d}\tau $ (6)

which is always valid in a plausible physical system.

Moreover, in general, we can write

$\epsilon \left(r,t\right)={\epsilon}_{o}\left[\delta \left(t\right)+\chi \left(r,t\right)\right]$ (7)

which is always valid because $\chi $ is an arbitrary function (or distribution). The parameter ${\epsilon}_{o}=8.875\times {10}^{-12}\text{\hspace{0.17em}}\text{F}/\text{m}$ is the vacuum permittivity and $\chi $ is the electric susceptibility expressed in temporal space. $\chi $ measures the amount of free charges induced by applying an electric field in the medium. Note that this phenomenon will be necessarily present in a heterogeneous medium such as neural tissue because we have charge accumulation in membranes when applied to an electric field. Note that the phenomenon of charge induction can change the law of attenuation with distance, as shown before [15] .

Thus, we can write

$\begin{array}{c}D\left(r,t\right)={\epsilon}_{o}{\displaystyle {\int}_{0}^{+\infty}}\left[\delta \left(t\right)+\chi \left(r,\tau \right)\right]E\left(r\mathrm{,}t-\tau \right)\text{d}\tau \\ ={\epsilon}_{o}E+{\displaystyle {\int}_{0}^{+\infty}}{\epsilon}_{o}\chi \left(r,\tau \right)E\left(r\mathrm{,}t-\tau \right)\text{d}\tau \mathrm{.}\end{array}$ (8)

It follows that, if the Fourier transform relative to
$\omega =2\text{\pi}\nu $ exists^{3}, then we have

$D\left(r\mathrm{,}\omega \right)={\epsilon}_{o}\left[1+\chi \left(x\mathrm{,}\omega \right)\right]E\left(r\mathrm{,}\omega \right)=\epsilon \left(r\mathrm{,}\omega \right)E\left(r\mathrm{,}\omega \right)$ (9)

^{3}Note that the Fourier transform may not necessarily exist for all possible mathematical models. However, in the case of heterogeneous media like neural tissue, we can safely assume that the Fourier transform of electric field and parameters associated to the medium exist because they have a finite duration. Thus, the time integral does not go to infinity, which means it will necessarily converge because the induced electric charge and susceptibility are bounded.

where

$\epsilon \left(r\mathrm{,}\omega \right)={\epsilon}_{o}\left[1+{\displaystyle {\int}_{0}^{+\infty}}\chi \left(x\mathrm{,}\tau \right){\text{e}}^{-i\omega \tau}\text{d}\tau \right]$ . (10)

We see that the electric permittivity in Fourier frequency space is in general a complex function, and its variation with respect to that of vacuum is such that

$\stackrel{\xaf}{\Delta}\epsilon =\epsilon -{\epsilon}_{o}={\displaystyle {\int}_{0}^{+\infty}}{\epsilon}_{o}\chi \left(x\mathrm{,}\tau \right){\text{e}}^{-i\omega \tau}\text{d}\tau \mathrm{.}$ (11)

The latter equation implies that we have $\stackrel{\xaf}{\Delta}\epsilon =0$ when $\omega \to \infty $ . Thus, the permittivity tends to that of vacuum when the frequency tends to infinity, which is in good agreement with experimental measurements. The real part of $\stackrel{\xaf}{\Delta}\epsilon \left(\omega \right)$ is an even function relative to the frequency $\omega $ and its imaginary part is an odd function because ${\epsilon}_{o}\chi \left(t\right)$ is a real function.

2.2.2. Analytic Continuation of we and the Kramers-Kronig or Hilbert Transform

We have shown that the imaginary part of $\stackrel{\xaf}{\Delta}\epsilon \left(r\mathrm{,}\omega \right)$ in Fourier frequency space must be odd and that the real part must be even if we want the system to be physically plausible. Moreover, it was shown that the real and imaginary parts of permittivity are linked by the Kramers-Kronig relations. However, we will develop this transform for $\omega \stackrel{\xaf}{\Delta}\epsilon \left(r\mathrm{,}\omega \right)$ instead of $\stackrel{\xaf}{\Delta}\epsilon \left(r\mathrm{,}\omega \right)$ to make the formalism more uniform between conductivity and permittivity.

According to expression (11), we have

$\omega \stackrel{\xaf}{\Delta}\epsilon ={\displaystyle {\int}_{0}^{+\infty}}\omega {\epsilon}_{o}\chi \left(x\mathrm{,}\tau \right){\text{e}}^{-i\omega \tau}\text{d}\tau $ . (12)

If we analytically continue the real frequency $\omega $ over the complex plane by setting $\omega ={\omega}^{\prime}+i{\omega}^{\u2033}$ then the last integral becomes

${\int}_{0}^{+\infty}}{\epsilon}_{o}\chi \left(r,\tau \right)\omega {\text{e}}^{-i\omega \tau}\text{d}\tau =\underset{\delta \to 0}{\mathrm{lim}}{\displaystyle {\int}_{\delta}^{1/\delta}}{\epsilon}_{o}\chi \left(r,\tau \right)\left[{\omega}^{\prime}+i\left({\omega}^{\u2033}-\delta \right)\right]{\text{e}}^{-i{\omega}^{\prime}\tau}{\text{e}}^{\left({\omega}^{\u2033}-\delta \right)\tau}\text{d}\tau $ (13)

where $\delta >0$ . The last integral converges when ${\omega}^{\u2033}\le 0$ and tends to zero when ${\omega}^{\u2033}\to -\infty $ . It diverges when ${\omega}^{\u2033}>0$ . Note that it is because we have applied the causality principle that the integral can converge for ${\omega}^{\u2033}<0$ ; omitting causality would prevent convergence because one would need to integrate between $-\infty $ and $+\infty $ (see Equation (5)).

Thus, the causality principle has the consequence that the function
$\omega \stackrel{\xaf}{\Delta}\epsilon $ is holomorphic in the complex inferior half-plane, as well as on the real axis^{4}. If we call this Region
$\mathcal{P}$ and apply the Cauchy integral of a holomorphic function, which gives

$\omega \stackrel{\xaf}{\Delta}\epsilon \left(r\mathrm{,}\omega \right)=\frac{1}{2\text{\pi}i}{\displaystyle \oint}\frac{{\omega}_{a}\stackrel{\xaf}{\Delta}\epsilon \left(r\mathrm{,}{\omega}_{a}\right)}{{\omega}_{a}-\omega}\text{d}{\omega}_{a}$ (14)

^{4}Note that we take into account that the permittivity may vary in time, and its Fourier transform reflects these time variations. But one could also consider a specific frequency profile of permittivity, and in that case, the temporal variations are given by its inverse Fourier transform (see ref. [8] [11] ). In the latter case, one must consider the holomorphic transformation in the superior half-plane.

if the path is closed within Region $\mathcal{P}$ , and if all points in this path are all situated at a finite distance from the origin. Note that, under this constraint, the function in the integral has only one singular point at ${\omega}_{a}=\omega $ because $\omega \stackrel{\xaf}{\Delta}\epsilon $ is a holomorphic function inside $\mathcal{P}$ . This is not the case for $\stackrel{\xaf}{\Delta}\epsilon $ in a conductive medium.

An interesting choice of path is a clockwise along the real axis, going around the singularity ${\omega}_{a}=\omega $ along a half-circle of radius $\delta $ and a large half-circle of radius $R=1/\delta $ (see Figure 1). In this case, there is only the principal value of the integral on the real axis, and we can write

$\omega \stackrel{\xaf}{\Delta}\epsilon \left(r\mathrm{,}\omega \right)=-\frac{1}{\text{\pi}i}{{{\displaystyle \overline{)\int}}}^{\text{}}}_{-\infty}^{+\infty}\frac{{\omega}_{a}\stackrel{\xaf}{\Delta}\epsilon \left(r\mathrm{,}\omega \right)}{{\omega}_{a}-\omega}\text{d}{\omega}_{a}$ . (15)

For ${\omega}_{a}\in \mathbb{R}$ and $\left|{\omega}_{a}\right|\to \infty $ , we have $\left|{\omega}_{a}\stackrel{\xaf}{\Delta}\epsilon \left(r\mathrm{,}\omega \right)\right|\to 0$ . Note that this condition is equivalent to postulate that $\left|\stackrel{\xaf}{\Delta}\epsilon \left(r\mathrm{,}{\omega}_{a}\right)\right|$ tends to 0 faster than $1/{\omega}_{a}$ when frequency tends to infinity. This insures that the integral in Equation (15) converges. This is obtained because for $R\to \infty $ ( $\delta \to 0$ ), the integral over the large circle is zero (see Equation 10) while the integral over the small circle is equal to $-\omega \stackrel{\xaf}{\Delta}\epsilon \left(r\mathrm{,}\omega \right)/2$ when $\delta \to 0$ (see Figure 1). We note the principal part of the integral by ${{\displaystyle \overline{)\int}}}^{\text{}}$ .

We can separate the real and imaginary parts of $\omega \stackrel{\xaf}{\Delta}\epsilon \left(r\mathrm{,}\omega \right)$ , leading to

$\{\begin{array}{l}\omega \stackrel{\xaf}{\Delta}{\epsilon}^{\u2033}\left(r\mathrm{,}\omega \right)=\mathcal{H}\left(\omega \stackrel{\xaf}{\Delta}{\epsilon}^{\prime}\right)=\frac{1}{\text{\pi}}{{{\displaystyle \overline{)\int}}}^{\text{}}}_{-\infty}^{+\infty}\frac{{\omega}_{a}\stackrel{\xaf}{\Delta}{\epsilon}^{\prime}\left(r,{\omega}_{a}\right)}{{\omega}_{a}-\omega}\text{d}{\omega}_{a}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{(}\text{a}\mathrm{)}\\ \omega \stackrel{\xaf}{\Delta}{\epsilon}^{\prime}\left(r\mathrm{,}\omega \right)=-\mathcal{H}\left(\omega \stackrel{\xaf}{\Delta}{\epsilon}^{\u2033}\right)=-\frac{1}{\text{\pi}}{{{\displaystyle \overline{)\int}}}^{\text{}}}_{-\infty}^{+\infty}\frac{{\omega}_{a}\stackrel{\xaf}{\Delta}{\epsilon}^{\u2033}\left(r,{\omega}_{a}\right)}{{\omega}_{a}-\omega}\text{d}{\omega}_{a}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{(}\text{b}\mathrm{)}\end{array}$ (16)

Note that Equations (16a) and (16b) are of opposite sign as the Kramers-Kronig relations as presented by Landau & Lifshitz [8] and Forster & Schwan [10] . The reason for this difference is that our analytic continuation of the electric parameters are holomorphic in the negative half-plane, instead of the positive half-plane [11] ^{5}.

^{5}The starting point of Landau & Lifshitz [8] is to assume that the permittivity depends on frequency a priori, which defines the temporal evolution of the permittivity as the direct Fourier transform of the permittivity defined in frequency space. In our approach, we consider that the permittivity is first defined in time, while the permittivity in frequency space is given by its direct Fourier transform. These two definitions are mathematically equivalent but with a change of sign. This is why in our formalism, we apply the inverse Fourier transform of the permittivity in frequency space to obtain the temporal evolution of the permittivity. The two approaches are equivalent, but our approach is more appropriate in the case there are time variations of electric parameters (such as the opening and closing of ion channels).

Figure 1. Illustration of the integration path in Region $\mathbb{P}$ . The singular point ${\omega}_{a}=\omega $ is on the real segment inside the half circle indicated. If $r\to 0$ and $R\to \infty $ , then applying the Cauchy integral gives the Hilbert transform.

Consequently, the real and imaginary parts of
$\omega \stackrel{\xaf}{\Delta}\epsilon \left(r\mathrm{,}\omega \right)$ are linked by the Hilbert transform. To calculate the real part from the knowledge of the imaginary part, one can apply the inverse Hilbert transform, while the direct transform is used to calculate the imaginary part from the real part. Note that this transform is totally equivalent to the Kramers-Kronig relations^{6}.

Finally, we can write

$\omega \stackrel{\xaf}{\Delta}\epsilon =\omega \stackrel{\xaf}{\Delta}{\epsilon}^{\prime}+i\mathcal{H}\left[\omega \stackrel{\xaf}{\Delta}{\epsilon}^{\prime}\right]=-\mathcal{H}\left[\omega \stackrel{\xaf}{\Delta}{\epsilon}^{\u2033}\right]+i\omega \stackrel{\xaf}{\Delta}{\epsilon}^{\u2033}$ . (17)

We can see that the phase $\varphi $ of $\omega \stackrel{\xaf}{\Delta}\epsilon $ is such that we have:

$\mathrm{tan}\left(\varphi \right)=\frac{\mathcal{H}\left[\omega \stackrel{\xaf}{\Delta}{\epsilon}^{\prime}\right]}{\omega \stackrel{\xaf}{\Delta}{\epsilon}^{\prime}}=-\frac{\omega \stackrel{\xaf}{\Delta}{\epsilon}^{\u2033}}{\mathcal{H}\left[\omega \stackrel{\xaf}{\Delta}{\epsilon}^{\u2033}\right]}$ .

Thus, the phase is completely determined when either the real or imaginary part of $\omega \stackrel{\xaf}{\Delta}\epsilon $ is known.

2.3. Electric Conductivity within the Quasi-Static Approximation

We now consider the temporal and frequency dependence of the electric conductivity in linear media within the electric quasi-static approximation. In the first section below, we review the constraints that must be satisfied in Fourier frequency space on the electric conductivity to simulate a physically plausible system. We show that the electric conductivity also obeys to a Hilbert transform, similar to that shown above for $\omega \stackrel{\xaf}{\Delta}\epsilon $ .

2.3.1. Electric Conductivity and Free-Charge Current in Linear Electromagnetism

It is well known that the most general relation between the free-charge current density ${j}^{\text{free}}\left(r\mathrm{,}t\right)$ and the electric field $E\left(r\mathrm{,}t\right)$ in linear electromagnetism is given by the convolution integral:

${j}^{\text{free}}\left(r\mathrm{,}t\right)={\displaystyle {\int}_{-\infty}^{+\infty}}{\sigma}_{e}\left(r,t-\tau \right)E\left(r,\tau \right)\text{d}\tau ={\displaystyle {\int}_{-\infty}^{+\infty}}{\sigma}_{e}\left(x,\tau \right)E\left(x,t-\tau \right)\text{d}\tau $ (18)

where
${\sigma}_{e}\ge 0$ ^{7} Expression (1) can be written as:

${j}^{\text{free}}\left(r\mathrm{,}\omega \right)={\sigma}_{e}\left(r\mathrm{,}\omega \right)E\left(r\mathrm{,}\omega \right)$ (19)

^{6}One can recover Kramers-Kronig by considering the parity of the real and imaginary parts.

^{7}This function is necessarily positive or zero, because experiments show that the current density is always in the same direction as the electric field. The terms under the integral are real functions, and the function
${\sigma}_{e}\left(t\right)$ is in [S/ms].

in Fourier frequency space. Here, the convolution in temporal space corresponds to a simple product in frequency space. Note that the function ${\sigma}_{e}\left(r\mathrm{,}\omega \right)$ depends on the value of the electric field when the relation between the current density and eclectic field is nonlinear, but this dependence vanishes if the system is linear.

Finally, Ohm’s law corresponds to the simplest model expressed by Equations (17), (18). In this particular case, we have
${\sigma}_{e}\left(x,t-\tau \right)=\sigma \left(x\right)\delta \left(t-\tau \right)$ where
$\sigma $ does not depend on time.
$\delta $ is the Dirac distribution. This law corresponds to an idealized physical system without memory (see Section 2), were the work produced by the electric field on the free charges dissipates almost instantaneously^{8} in the system. Note we also have in this case
${j}^{\text{free}}\left(r\mathrm{,}t\right)=\sigma \left(r\right)E\left(r\mathrm{,}t\right)$ in temporal space, and we have a similar relation
${j}^{\text{free}}\left(r\mathrm{,}\omega \right)=\sigma \left(r\right)E\left(r\mathrm{,}\omega \right)$ in frequency space. Thus, we have the same algebraic relation between
${j}^{\text{free}}$ and
$E$ in both spaces.

2.3.2. Constraints Imposed by the Causality Principle

To be physically plausible, a system must obey the causality principle. This principle determines a constraint on the relation between the free-charge current density ${j}^{\text{free}}$ and the electric field $E$ (within the linear electromagnetic theory). According to this principle, the future cannot influence the past, and thus we can write that, at a given time t, ${j}^{\text{free}}$ and $E$ are related as

${j}^{\text{free}}\left(r\mathrm{,}t\right)={\displaystyle {\int}_{0}^{+\infty}}{\sigma}_{e}\left(r,\tau \right)E\left(r,t-\tau \right)\text{d}\tau $ (20)

because the values of the electric field at times greater than t cannot influence the current density at time t. Note that this is equivalent to assume that ${\sigma}_{e}\left(r\mathrm{,}\tau \right)=0$ when $\tau <0$ in Equation (19). This constraint is general and must be included in all mathematical models of conductivity to be physically plausible.

2.3.3. Electrical Conductivity in a Heterogeneous Medium

^{8}“Almost instantaneously” means that the dissipation of the energy brought by the electric field is produced at
${0}^{+}=0+\left|dt\right|$ .

We now consider the case of a heterogeneous medium composed of different cells and various processes immersed in a conductive fluid, such that the distance between different elements ${\delta}_{c}$ is always greater than zero. We also assume that the electric conductivity of this medium tends asymptotically to ${\sigma}_{\infty}\left(r,t\right)$ when the frequency $\nu $ tends to infinity. For high frequencies, there is a portion of free charges which does not meet any process, because their mean displacement becomes smaller than ${\delta}_{c}$ . However, for sufficiently low frequencies, the presence of cells will impact all free charges, and the conductivity will be affected and will be different as that of high frequencies. Thus, the conductivity of a heterogeneous medium will necessarily be frequency dependent within some frequency range.

This intuitive explanation can be formulated more quantitatively. The electric conductivity can be expressed as

${\sigma}_{e}\left(r,t\right)={\sigma}_{\infty}\left(r,t\right)\left[1+{\stackrel{\xaf}{\Delta}}_{e}\left(r,t\right)\right]$ (21)

where ${\stackrel{\xaf}{\Delta}}_{e}\left(r\mathrm{,}t\right)\ne 0$ . It follows that Expression (5) can be written as

${j}^{\text{free}}\left(r\mathrm{,}t\right)=\mathrm{}{\displaystyle {\int}_{0}^{+\infty}}{\sigma}_{\infty}\left(x,\tau \right)\left[1+{\stackrel{\xaf}{\Delta}}_{e}\left(x,\tau \right)\right]E\left(x,t-\tau \right)\text{d}\tau $ . (22)

Taking the Fourier transform, we obtain

${j}^{\text{free}}\left(r\mathrm{,}\omega \right)=\left[{\sigma}_{\infty}\left(r\mathrm{,}\omega \right)+{\displaystyle {\int}_{0}^{+\infty}}{\sigma}_{\infty}\left(x,\tau \right){\stackrel{\xaf}{\Delta}}_{e}\left(x,\tau \right){\text{e}}^{-i\omega \tau}\text{d}\tau \right]E\left(r\mathrm{,}\omega \right)$ . (23)

Thus, we can write

${\sigma}_{e}\left(r\mathrm{,}\omega \right)={\sigma}_{\infty}\left(r\mathrm{,}\omega \right)+{\displaystyle {\int}_{0}^{+\infty}}{\sigma}_{\infty}\left(r,\tau \right){\stackrel{\xaf}{\Delta}}_{e}\left(r,\tau \right)\mathrm{}{\text{e}}^{-i\omega \tau}\text{d}\tau $ (24)

where the function ${\sigma}_{\infty}\left(r,\tau \right){\stackrel{\xaf}{\Delta}}_{e}\left(r,\tau \right)$ is such that the integral in Expression (24) tends to zero when $\omega \to \infty $ . Thus, ${\sigma}_{e}\left(r\mathrm{,}\omega \to \infty \right)$ tends to ${\sigma}_{\infty}\left(r\mathrm{,}\omega \right)$ . Note that the integral in the right-hand side converges because it is the Fourier transform of $\stackrel{\xaf}{\Delta}{\sigma}_{e}={\sigma}_{e}-{\sigma}_{\infty}$ .

2.3.4. Similar Relations between $\omega \stackrel{\xaf}{\Delta}\epsilon $ and $\stackrel{\xaf}{\Delta}{\sigma}_{e}$

We know that the imaginary part of $\stackrel{\xaf}{\Delta}{\sigma}_{e}$ in frequency space is an odd function, while the real part is even, if the model is physically plausible because the time-dependent conductivity is a real function. However, these parts could still be independent of each-other. We now show that if we apply the causality principle, the imaginary part of $\stackrel{\xaf}{\Delta}{\sigma}_{e}$ (in frequency space) is completely determined by its real part, via a Hilbert transform, exactly like $\omega \stackrel{\xaf}{\Delta}\epsilon $ .

If we analytically continue the frequency $\omega $ over the complex plane by taking $\omega ={\omega}^{\prime}+i{\omega}^{\u2033}$ , then the integral in Expression (23) becomes

${\int}_{0}^{+\infty}}{\sigma}_{\infty}\left(r,\tau \right){\stackrel{\xaf}{\Delta}}_{e}\left(r,\tau \right){\text{e}}^{-i\omega \tau}\text{d}\tau =\underset{\delta \to 0}{lim}{\displaystyle {\int}_{\delta}^{1/\delta}}{\sigma}_{\infty}\left(r,\tau \right){\stackrel{\xaf}{\Delta}}_{e}\left(r,\tau \right)\mathrm{}{\text{e}}^{-i{\omega}^{\prime}\tau}{\text{e}}^{\left({\omega}^{\u2033}-\delta \right)\tau}\text{d}\tau $ (25)

where $\delta >0$ . This integral converges when ${\omega}^{\u2033}<0$ and tends to 0 for ${\omega}^{\u2033}\to -\infty $ . On the other hand, it diverges when ${\omega}^{\u2033}>0$ . Note that the principle of causality allows the convergence of this integral when the imaginary part of $\omega $ is smaller than zero. This would not be the case if the causality principle is not used, because one would need to integrate between $-\infty $ and $+\infty $ .

Thus, the situation is completely analogous to the analytic continuation of $\omega \stackrel{\xaf}{\Delta}\epsilon $ , and we can write:

$\{\begin{array}{l}\stackrel{\xaf}{\Delta}{{\sigma}^{\u2033}}_{e}=\mathcal{H}\left(\stackrel{\xaf}{\Delta}{{\sigma}^{\prime}}_{e}\right)=\frac{1}{\text{\pi}}{{{\displaystyle \overline{)\int}}}^{\text{}}}_{-\infty}^{+\infty}\frac{\stackrel{\xaf}{\Delta}{{\sigma}^{\prime}}_{e}\left(r,{\omega}_{a}\right)}{{\omega}_{a}-\omega}\text{d}{\omega}_{a}\mathrm{}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}(\text{a})\\ \stackrel{\xaf}{\Delta}{{\sigma}^{\prime}}_{e}=-\mathcal{H}\left(\stackrel{\xaf}{\Delta}{{\sigma}^{\u2033}}_{e}\right)=-\frac{1}{\text{\pi}}{{{\displaystyle \overline{)\int}}}^{\text{}}}_{-\infty}^{+\infty}\frac{\stackrel{\xaf}{\Delta}{{\sigma}^{\u2033}}_{e}\left(r,{\omega}_{a}\right)}{{\omega}_{a}-\omega}\text{d}{\omega}_{a}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.17em}}(\text{b})\mathrm{}\end{array}$ (26)

${\sigma}_{e}\left(r\mathrm{,}\omega \right)={{\sigma}^{\prime}}_{e}\left(r\mathrm{,}\omega \right)+i{{\sigma}^{\u2033}}_{e}\left(r\mathrm{,}\omega \right)$ and ${\sigma}_{\infty}={{\sigma}^{\prime}}_{\infty}+i{{\sigma}^{\u2033}}_{\infty}$ . ${{\sigma}^{\prime}}_{e}\left(r\mathrm{,}\omega \right)$ and ${{\sigma}^{\prime}}_{\infty}$ are respectively the real part of ${\sigma}_{e}$ and ${\sigma}_{\infty}$ , while ${{\sigma}^{\u2033}}_{e}\left(r\mathrm{,}\omega \right)$ and ${{\sigma}^{\u2033}}_{\infty}$ are respectively their imaginary part. Consequently, if we know the conductivity at very high frequencies and the real part of its frequency dependence, then we can determine the corresponding imaginary part using the Hilbert transform, and vice et versa. Finally, note that we apply here the Hilbert transform to the variation of conductivity relative to ${\sigma}_{e}\left(r\mathrm{,}\infty \right)$ to make sure that the integral converges when the frequency tends to infinity.

2.4. Apparent Conductivity and Permittivity between Two Isopotential Surfaces

Within the electric quasi-static approximation, the electric field can be expressed as the gradient of the potential,
$E=-\nabla V$ . We also know that the free-charge current is not conserved in general in a heterogeneous medium, for example because charge accumulation can occur^{9} The density of the generalized current

is such that $\nabla \cdot {j}^{g}=0$ , where ${j}^{g}={j}^{\text{free}}+\frac{\partial D}{\partial t}$ . Note that this does not represent

a stationary law, because
${j}^{g}$ is time dependent, but it is rather a conservation law [16] . The generalized current entering a given domain is always equal to the generalized current exiting that domain, even if charge accumulation occurs^{10}.

Thus, according to these laws, the generalized current density between two close-by equipotential surfaces is given by:

${j}^{g}\left(r\mathrm{,}\omega \right)=-\gamma \left(r\mathrm{,}\omega \right)\nabla V\left(r\mathrm{,}\omega \right)=-\left({\sigma}_{e}+i\omega \epsilon \right)\nabla V\left(r\mathrm{,}\omega \right)$ (27)

where $\gamma $ is the admittance of the medium. ${\sigma}_{e}$ is the link between ${j}^{\text{free}}$ and $E$ in frequency space, while $\epsilon $ is the link between $E$ and $D$ . The latter link allows one to calculate the induced charge in a given region of the medium by a given electric field, expressed in Fourier frequency space. Note that these links are all complex numbers in general because there can be a non-zero phase, similar to a capacitance.

If we now assume the following equalities:

$\gamma \left(r,\infty \right)=\sigma \left(r,\infty \right)+i\omega \epsilon \left(r,\infty \right)\mathrm{}et\mathrm{}\stackrel{\xaf}{\Delta}\gamma =\gamma \left(r\mathrm{,}\omega \right)-{\gamma}_{\infty}$ . (28)

Applying Equations (16) and (25), we can write:

$\stackrel{\xaf}{\Delta}\gamma =\stackrel{\xaf}{\Delta}{\gamma}^{\prime}+i\stackrel{\xaf}{\Delta}{\gamma}^{\u2033}=\left[\stackrel{\xaf}{\Delta}{{\sigma}^{\prime}}_{e}-\mathcal{H}\left(\omega \stackrel{\xaf}{\Delta}{\epsilon}^{\prime}\right)\right]+i\mathrm{}\left[\omega \stackrel{\xaf}{\Delta}{\epsilon}^{\prime}+\mathcal{H}\left(\stackrel{\xaf}{\Delta}{{\sigma}^{\prime}}_{e}\right)\right]$ . (29)

It follows that

$\{\begin{array}{l}\stackrel{\xaf}{\Delta}{\gamma}^{\u2033}=\mathcal{H}\left(\stackrel{\xaf}{\Delta}{\gamma}^{\prime}\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}(\text{a})\\ \stackrel{\xaf}{\Delta}{\gamma}^{\prime}=-\mathcal{H}\left(\stackrel{\xaf}{\Delta}{\gamma}^{\u2033}\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}(\text{b})\end{array}$ (30)

because the Hilbert transform obeys $-{\mathcal{H}}^{2}\left(F\right)=F$ . We can also write

$\stackrel{\xaf}{\Delta}\gamma =\stackrel{\xaf}{\Delta}{\gamma}^{\prime}+i\mathcal{H}\left(\stackrel{\xaf}{\Delta}{\gamma}^{\prime}\right)$ (31)

^{9}In the Debye layers of ions around a membrane, for example.

^{10}Note that this is only true at chemical steady-state, when there is no current created by chemical reactions (see discussion in ref. [17] ).

such that the phase $\varphi $ obeys

$\mathrm{tan}\left(\varphi \right)=\frac{\mathcal{H}\left(\stackrel{\xaf}{\Delta}{\gamma}^{\prime}\right)}{\stackrel{\xaf}{\Delta}{\gamma}^{\prime}}$ . (32)

For a sufficiently large base volume, we can always assume that the conductivity and permittivity do not depend on position. In this particular case, we can write

$\{\begin{array}{l}{I}^{g}={\displaystyle \underset{{\mathcal{S}}_{eq}}{\iint}}j\cdot \stackrel{^}{n}\text{d}S={\langle \Vert j\Vert \rangle}_{\mathcal{S}}S\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}(\text{a})\\ \Delta V={\displaystyle {\int}_{a}^{b}}\nabla V\cdot \stackrel{^}{n}\text{d}s=\frac{{\langle \Vert j\Vert \rangle}_{l}}{\gamma \left(\omega \right)}l\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}(\text{b})\end{array}$ (33)

where the current ${I}^{g}$ is conserved. ${\langle \Vert j\Vert \rangle}_{\mathcal{S}}$ is the mean current density over an equipotential surface and ${\langle \Vert j\Vert \rangle}_{l}$ is the mean current density along a current line between two equipotential surfaces. These two means are different in general, but if we take the curve that goes through the mean current density ${\langle \Vert j\Vert \rangle}_{\mathcal{S}}$ over each equipotential surface, then we can say that this curve is a current line. Consequently, there exists an equipotential surface and a current line such that ${\langle \Vert j\Vert \rangle}_{\mathcal{S}}={\langle \Vert j\Vert \rangle}_{l}$ . It follows that, if ${A}_{\mathcal{C}}$ is the area of this surface, and ${l}_{\mathcal{C}}$ is the distance between two isopotential surfaces, then we can write

$\frac{{I}^{g}}{\Delta V}=\gamma \left(\omega \right)\mathrm{}\frac{{A}_{\mathcal{C}}}{{l}_{\mathcal{C}}}$ . (34)

This form is identical to that of a plane capacitor which possesses a leak current, because the real part is non-zero, or a conductance with a non-negligible capacitive effect, because the imaginary part is also non-zero.

If we define the apparent conductivity as ${\sigma}_{\mathcal{A}}$ as the real part of $\gamma $

${\sigma}_{\mathcal{A}}={\gamma}^{\prime}={{\sigma}^{\prime}}_{e}-\omega {\epsilon}^{\u2033}$ (35)

and the apparent permittivity ${\epsilon}_{\mathcal{A}}$ times the angular frequency $\omega $ as the imaginary part of $\gamma $

$\omega {\epsilon}_{\mathcal{A}}={\gamma}^{\u2033}=\omega {\epsilon}^{\prime}+{{\sigma}^{\u2033}}_{e}\mathrm{},$ (36)

we can then write where these apparent parameters are real and linked by a Hilbert transform. We have $\omega \stackrel{\xaf}{\Delta}{\epsilon}_{\mathcal{A}}=\mathcal{H}\left(\stackrel{\xaf}{\Delta}{\sigma}_{\mathcal{A}}\right)$ and $\stackrel{\xaf}{\Delta}{\sigma}_{\mathcal{A}}=-\mathcal{H}\left(\omega \stackrel{\xaf}{\Delta}{\epsilon}_{\mathcal{A}}\right)$ .

Thus, the knowledge of one parameter is sufficient to deduce the other. For a given frequency, it is always possible to simulate the current-voltage relation between two isopotential surfaces as a plane capacitor (see Figure 2). The leak current of this capacitor is determined by the Hilbert transform of its admittance, or by a conductance possessing a capacitive effect determined by Hilbert transform of the admittance. Finally, ${\sigma}_{\mathcal{A}}$ and ${\epsilon}_{\mathcal{A}}$ are an even function (see Equations (35) and (36)) relative to $\omega $ .

Figure 2. Capacitance with a leak current ${I}_{f}$ . In Fourier frequency space, the voltage difference V and the current density j obey $j=\left[{\sigma}_{e}+i\omega \epsilon \right]\mathrm{}\frac{V}{d}$ when the electric parameters do not depend on frequency. d is the distance between the arms of the capacitor. The density of the leak current is given by ${j}_{f}={\sigma}_{e}V$ that of the capacitive current is ${j}_{c}=i\omega \epsilon V$ . We also have ${I}_{f}={j}_{f}A$ and ${I}_{c}={j}_{c}A$ where A is the area of the arms.

It is important to note that, for the apparent parameters, as defined in Equations (35) and (36), we have the following properties. First, the apparent conductivity is equal to the electric conductivity when $\omega =0$ , or when the imaginary part of the electric permittivity is zero. Second, the apparent permittivity becomes infinitely large for $\omega =0$ when the imaginary part of the electric conductivity is different from zero. Third, if the imaginary parts of both electric conductivity and permittivity are zero, then the apparent parameters are identical to electric parameters. Finally, if the physical effects associated to the electric permittivity (density of induced charges) are negligible compared to that of electric conductivity, then the apparent conductivity is approximately given by the real part of electric conductivity, and the apparent permittivity is approximately to the imaginary part of electric conductivity divided by $\omega $ .

3. Applications

In this section, we consider models of different known media to illustrate the consequences of the relations outlined in the theoretical part. All models considered are built within the quasistatic approximation of the linear electromagnetic theory of Maxwell-Heaviside. In other words, we do not consider physical phenomena associated to electromagnetic induction (
$\nabla \times E=0$ ). In such conditions, we can apply the Kramers-Kronig relations (or Hilbert transform) over the approximated electric parameters^{11}.

3.1. Constraints on the Conductivity Given by the Kramers-Kronig Relations

3.1.1. First Example

As a first example, we consider an isotropic medium where the base volume of the mean-field does not depend on position. We suppose that the medium is such that ${\sigma}_{e}$ depends on frequency with a null imaginary part, and that $\underset{\omega \to \infty}{\mathrm{lim}}{\sigma}_{e}=k$ is real and frequency independent. This model would correspond to a heterogeneous and isotropic medium where the conductivity depends on frequency $\nu $ according to a law ${\sigma}_{e}\left(\omega =2\text{\pi}\nu \right)={{\sigma}^{\prime}}_{e}\left(\omega \right)+i{{\sigma}^{\u2033}}_{e}\left(\omega \right)={{\sigma}^{\prime}}_{e}\left(\omega \right)$ .

In this case, according to Equation (26b), we must have

${{\sigma}^{\prime}}_{e}\left(\omega \right)=\sigma \left(\infty \right)=k$ (37)

^{11}By “approximated” parameters, we mean the electric parameters within the quasistatic electric approximation. They constitute an excellent approximation for low frequencies, <~10 kHz.

because $\stackrel{\xaf}{\Delta}{{\sigma}^{\prime}}_{e}=0$ when $\stackrel{\xaf}{\Delta}{{\sigma}^{\u2033}}_{e}=0$ . Therefore, according to Kramers-Kronig relations, ${{\sigma}^{\prime}}_{e}$ would not depend on frequency, which is contradictory with the initial hypothesis.

Thus, a model of this kind is not acceptable, because we have a contradiction with the causality principle because the Kramers-Kronig relations are a direct consequence of this principle in a linear system.

For example, it is impossible to have a conductivity law of the form ${\sigma}_{e}={\omega}^{x}{\text{e}}^{-\omega /c}+k$ where x is arbitrary. One must add an imaginary part, which implies that the phase is non-zero. Thus, a model with a frequency-dependent conductivity but a null phase is physically impossible in a model where conductivity becomes resistive or if the phase becomes zero at high frequencies. Note that, in general, the electric conductivity is not necessarily real for large frequencies. It was shown before that the linear approximation of a physical system with ionic diffusion gives an electric conductance of the form $\sqrt{\omega}\left(a+ib\right)$ [7] . Thus, the asymptotic behavior of the electric conductivity at high frequencies is very different from the electric permittivity, because the latter tends to the permittivity of vacuum, which is real.

For example, let us consider a model of a medium in which the electric conductivity is frequency-dependent, but the permittivity is constant. According to above, such a medium is physically impossible, because it would violate the principle of causality. It is important to note that the standard model of a resistive extracellular medium and the capacitive effects are neglected.

This model is asymptotically resistive at high frequencies ( $\underset{\omega \to \infty}{\mathrm{lim}}{\sigma}_{e}=k$ ), and thus, we can conclude that, if the conductivity depends on frequency, one must necessarily have both real and imaginary parts different from zero for the model to be physically plausible.

3.1.2. Second Example

In this section, we suppose that the apparent macroscopic permittivity ${\epsilon}_{\mathcal{A}}$ between two equipotential surfaces is given by

${\epsilon}_{\mathcal{A}}={\epsilon}_{o}+\frac{\kappa}{{\left|\omega \right|}^{2+x}}$ (38)

where $x\mathrm{,}\kappa \in \mathbb{R}$ and $x\ge -1$ , and $\kappa \ge 0$ and frequency-independent. Note that this model respects the parity of ${\epsilon}_{\mathcal{A}}$ (see Equation (36)) because we have ${\epsilon}_{\mathcal{A}}\left(+\omega \right)={\epsilon}_{\mathcal{A}}\left(-\omega \right)$ . This simple model corresponds to a physical situation where the apparent permittivity diminishes for increasing frequency, and approaches that of vacuum for very high frequencies. This situation is often encountered in heterogeneous media such as biological tissues (see for example, the experimental measurements of Gabriel et al. [1] ).

Note that if the apparent permittivity ${\epsilon}_{\mathcal{A}}$ tends to that of vacuum, then the electric conductivity ${\sigma}_{e}$ must necessarily have a zero phase when frequency tends to infinity (see Equation 36) because the real part of the electric permittivity $\epsilon $ tends to that of vacuum.

By applying Equation (30b) and assuming that $y=\frac{{\omega}_{a}}{w}$ , we obtain:

$\stackrel{\xaf}{\Delta}{\sigma}_{\mathcal{A}}\left(\omega \right)=-\mathcal{H}\left(\stackrel{\xaf}{\Delta}\omega {\epsilon}_{\mathcal{A}}\right)=-\frac{2\kappa}{\text{\pi}{\omega}^{1+x}}{{{\displaystyle \overline{)\int}}}^{\text{}}}_{0}^{+\infty}\frac{1}{\left({y}^{2}-1\right){y}^{x}}\text{d}y=-\frac{{\kappa}_{1}}{{\omega}^{1+x}}$ (39)

where we have

${\kappa}_{1}=\frac{2\kappa}{\text{\pi}}{{{\displaystyle \overline{)\int}}}^{\text{}}}_{0}^{+\infty}\frac{1}{\left({y}^{2}-1\right){y}^{x}}\text{d}y$ . (40)

Consequently, the phase of the admittance variation is such that

$\mathrm{tan}\left(\varphi \right)=\frac{\stackrel{\xaf}{\Delta}\left(\omega {\epsilon}_{\mathcal{A}}\right)}{\stackrel{\xaf}{\Delta}\left({\sigma}_{\mathcal{A}}\right)}=-\frac{\kappa}{{\kappa}_{1}}$ (41)

when $\omega >0$ and $+\frac{\kappa}{{\kappa}_{1}}$ for $\omega <0$ . Note that ${\kappa}_{1}$ depends on the value of the

exponent x. Also note that the phase of the admittance variation $\Delta \gamma $ does not depend on frequency in this type of model.

We now calculate the phase for different values of the exponent x. The principal part of the integral in expression (40) is explicitly given by the following expression:

${\kappa}_{1}=\underset{\delta \to 0}{\mathrm{lim}}\left[{\displaystyle {\int}_{\delta}^{1-\delta}}\frac{1}{\left({y}^{2}-1\right){y}^{x}}\text{d}y+{\displaystyle {\int}_{1+\delta}^{1/\delta}}\frac{1}{\left({y}^{2}-1\right){y}^{x}}\text{d}y\right]$ (42)

where $\delta >0$ .

Case x = 0.

For $x=0$ , Equation (42) becomes

${\kappa}_{1}=\frac{2\kappa}{\text{\pi}}\mathrm{}\underset{\delta \to 0}{\mathrm{lim}}\mathrm{ln}\left[\frac{-\delta}{2-\delta}\cdot \frac{\delta +1}{\delta -1}\cdot \frac{1-\delta}{1+\delta}\cdot \frac{2+\delta}{\delta}\right]=0$ . (43)

It follows that the phase of the admittance variation is given by $\varphi =-\frac{\text{\pi}}{2}$ for ( $\omega >0$ ), and $\varphi =\frac{\text{\pi}}{2}$ for ( $\omega <0$ ).

Moreover, we see that the variations of apparent conductivity with frequency are such that $\stackrel{\xaf}{\Delta}{\sigma}_{\mathcal{A}}=0$ . Thus, in this model, the apparent conductivity is equal to the real part of the electric conductivity because, for infinite frequencies, the imaginary part of electric permittivity ${\epsilon}^{\u2033}$ tends to zero and ${\epsilon}^{\prime}$ tends to that of

vacuum. Thus, we have ${\epsilon}_{\mathcal{A}}={\epsilon}_{o}+\frac{\kappa}{{\left|\omega \right|}^{2}}$ and ${\sigma}_{\mathcal{A}}={{\sigma}^{\prime}}_{e}\left(\omega \to \infty \right)$ . Consequently, if

we assume that the medium is resistive for very high frequencies, then we see that the phase is necessarily frequency-dependent because the phase of the

admittance $\gamma $ is given by $\mathrm{tan}\left(\varphi \right)=\frac{\omega {\epsilon}_{0}+\frac{\kappa}{\omega}}{{{\sigma}^{\prime}}_{e}}$ for $\omega >0$ .

Case $x=\pm 1/2$ .

In this case, the primitives of the integrals in expression (42) are:

$\frac{1}{2}\mathrm{ln}\mathrm{}\left[\frac{1-\sqrt{y}}{1+\sqrt{y}}\right]\mp {\text{tg}}^{-1}\sqrt{y}$ . (44)

It follows that

${\kappa}_{1}=\frac{2\kappa}{\text{\pi}}\mathrm{}\underset{\delta \to 0}{\mathrm{lim}}\mathrm{ln}\mathrm{}\left[\frac{1-\sqrt{1-\delta}}{1+\sqrt{1+\delta}}\cdot \frac{1+\sqrt{\delta}}{1-\sqrt{\delta}}\cdot \frac{1-\sqrt{1/\delta}}{1+\sqrt{1/\delta}}\cdot \frac{1+\sqrt{1+\delta}}{1-\sqrt{1+\delta}}\right]+t=t$ (45)

where

$t=\mp \frac{2\kappa}{\text{\pi}}\mathrm{}\underset{\delta \to 0}{\mathrm{lim}}\left[{\text{tg}}^{-1}\sqrt{1-\delta}-{\text{tg}}^{-1}\sqrt{\delta}+{\text{tg}}^{-1}\sqrt{1/\delta}-{\text{tg}}^{-1}\sqrt{1+\delta}\right]=\mp \kappa $ .

Thus, the phase (Equation (41)) of the admittance variation is given by

$\varphi =\pm \frac{\text{\pi}}{4}$ , because we have

$\mathrm{tan}\left(\varphi \right)=\frac{\stackrel{\xaf}{\Delta}\left(\omega {\epsilon}_{\mathcal{A}}\right)}{\stackrel{\xaf}{\Delta}\left({\sigma}_{\mathcal{A}}\right)}=\pm 1$ . (46)

Consequently, the admittance variation $\stackrel{\xaf}{\Delta}{\gamma}_{\mathcal{A}}$ is such that $\stackrel{\xaf}{\Delta}{\gamma}_{\mathcal{A}}=\stackrel{\xaf}{\Delta}\left({\sigma}_{\mathcal{A}}\right)\left(1-i\right)$ and $\stackrel{\xaf}{\Delta}{\gamma}_{\mathcal{A}}=\stackrel{\xaf}{\Delta}\left({\sigma}_{\mathcal{A}}\right)\left(1+i\right)$ when we have, respectively, ${\epsilon}_{\mathcal{A}}={\epsilon}_{o}+\frac{\kappa}{{\left|\omega \right|}^{3/2}}$ and ${\epsilon}_{\mathcal{A}}={\epsilon}_{o}+\frac{\kappa}{{\left|\omega \right|}^{5/2}}$ .

This result is interesting for the following reason: we note that if we consider the case where the exponent equals −3/2 for $\stackrel{\xaf}{\Delta}{\epsilon}_{\mathcal{A}}$ , then we have $\stackrel{\xaf}{\Delta}{\sigma}_{\mathcal{A}}<0$ . This corresponds to a situation where conductivity increases with frequency. On the other hand, if the exponent equals −5/2, then the apparent conductivity decreases with frequency. The first case seems in agreement with experimental measurements (see [1] ).

It follows that, for the first case (exponent of −3/2), we can write:

${\gamma}_{\mathcal{A}}=\stackrel{\xaf}{\Delta}{\gamma}_{\mathcal{A}}+{\gamma}_{\mathcal{A}}\left(\omega =\infty \right)=\stackrel{\xaf}{\Delta}\left({\sigma}_{\mathcal{A}}\right)\left(1-i\right)+{\gamma}_{\mathcal{A}}\left(\omega =\infty \right)$ . (47)

We see that if the apparent admittance tends asymptotically to a constant C (independent of frequency) for high frequencies, then we have $C\approx \stackrel{\xaf}{\Delta}\left({\sigma}_{\mathcal{A}}\right)$ , so there will be a non-negligible phase at low frequencies. We recall that, by

definition, we have $\underset{\omega \to \infty}{\mathrm{lim}}\stackrel{\xaf}{\Delta}\left({\sigma}_{\mathcal{A}}\right)=0$ (see Equation (28)).

We also see that if the product $\omega \epsilon $ is negligible compared to electric conductivity, and if the latter tends to a Warburg impedance for high frequencies (as in ref. [18] ), then, according to the definition of apparent parameters (Equations (35) and (36)), we can write:

${\gamma}_{\mathcal{A}}={\gamma}_{e}=\left(\stackrel{\xaf}{\Delta}\left({\sigma}_{\mathcal{A}}\right)+k\right)\sqrt{\omega}\left(1\pm i\right)$ (48)

where ${\gamma}_{e}$ is the electric admittance.

Case x arbitrary and greater than −1.

For an arbitrary value of x, we can write

$\mathrm{tan}\left(\varphi \right)=\frac{\stackrel{\xaf}{\Delta}\left(\omega {\epsilon}_{\mathcal{A}}\right)}{\stackrel{\xaf}{\Delta}\left({\sigma}_{\mathcal{A}}\right)}=cst$ (49)

such that we can write

${\gamma}_{\mathcal{A}}=\stackrel{\xaf}{\Delta}{\gamma}_{\mathcal{A}}+{\gamma}_{\mathcal{A}}\left(\omega =\infty \right)=\stackrel{\xaf}{\Delta}\left({\sigma}_{\mathcal{A}}\right)\left(1+\mathrm{tan}\left(\varphi \right)i\right)+{\gamma}_{\mathcal{A}}\left(\omega =\infty \right)$ (50)

where the phase ${\gamma}_{\mathcal{A}}$ is frequency dependent because $\stackrel{\xaf}{\Delta}\left({\sigma}_{\mathcal{A}}\right)$ depends on frequency.

Finally, we see that if the electric permittivity varies according to the following law:

${\epsilon}_{\mathcal{A}}={\epsilon}_{o}+{\displaystyle \underset{j=1}{\overset{N}{\sum}}}\frac{{\kappa}_{j}}{{\left|\omega \right|}^{2+{x}_{j}}}$ (51)

where ${x}_{j}>-1$ , then we can calculate the phase of the apparent electric admittance, by fitting the expression (51) to the measured electric permittivity. Note that in these examples, this phase is non-negligible when the electric or apparent parameters (admittance, conductivity, permittivity) are frequency dependent.

4. Discussion

In this paper, we have re-examined electromagnetism theory in heterogeneous media, in particular focusing on the Kramers-Kronig relations, following the Landau & Lifshitz formalism [8] . Our main finding is that, similar to the well-known Kramers-Kronig relations linking the real and imaginary parts of the electric permittivity [8] , one can also derive similar relations for electric conductivity. Thus, in heterogeneous media, the two electric parameters obey symmetric dependencies in Fourier frequency space. This finding is general, and also applies to a homogeneous medium (which is a particular case of heterogeneous media). We discuss below the significance of these results.

As a first example, we considered the model of a medium where electric conductivity is assumed to be frequency-dependent, but with a constant permittivity. We showed that, according to Kramers-Kronig relations, this model is physically impossible, as it would violate the principle of causality. This model is also asymptotically resistive at high frequencies ( $\underset{\omega \to \infty}{\mathrm{lim}}{\sigma}_{e}=k$ ), and thus, we can conclude that, if the conductivity depends on frequency, one must necessarily have both real and imaginary parts different from zero for such a model to be physically plausible. It is important to note that the NEURON simulator [19] used for realistic neuronal simulations cannot deal with complex electric parameters (non-negligible phase), so cannot be used to simulate physically-plausible non-resistive situations.

In a second example, we considered experimental measurements suggesting that the extracellular medium around neurons is non-resistive [3] [5] . We showed here that these measurements are consistent with Kramers-Kronig relations, and the principle of causality. However, other measurements suggesting resistive media [2] [4] , are also consistent with Kramers-Kronig relations. All these experimental measurements are thus physically plausible and self-consistent. On the modeling point of view, the models of non-resistive media [3] [6] [7] [16] , as well as those of resistive media [4] , are all consistent with Kramers-Kronig relations as well. Note that ionic diffusion was proposed as a mechanism to explain the non-resistive measurements [3] [5] [6] , and according to this mechanism, the phase of the apparent admittance should tend to $\text{\pi}/4$ at high frequencies (Warburg impedance). This value seems in agreement with phase measurements in cerebral cortex [3] [20] and retina [18] , and is also consistent with Kramers-Kronig relations.

A further interesting property, developed Section 2.4, is that the apparent permittivity is far from negligible at very low frequencies when the imaginary part of the electric conductivity is non zero, even if it is very small. This property must be related to the difficulty of measuring the impedance (or admittance) at low frequencies [1] [2] [3] [20] . This difficulty could be the sign that the apparent permittivity increases rapidly at low frequencies (lower than ~10 Hz). This phenomenon should not appear in a resistive medium.

5. Conclusion

We conclude that, given the constraints imposed by Kramers-Kronig relations, the previous experimental measurements seem internally consistent for frequencies larger than 10 Hz, either for resistive [2] [4] or non-resistive media [1] [3] . Further experiments are needed to distinguish between these two alternatives, as reviewed in detail previously [5] .

Acknowledgements

Research funded by the CNRS, the Paris-Saclay excellence network (IDEX) and the European Community (Human Brain Project, H2020-720270).

References

[1] Gabriel, S., Lau, R.W. and Gabriel, C. (1996) The Dielectric Properties of Biological Tissues: II. Measurements in the Frequency Range 10 Hz to 20 GHz. Physics in Medicine & Biology, 41, 2251-2269.

https://doi.org/10.1088/0031-9155/41/11/002

[2] Logothetis, N.K., Kayser, C. and Oeltermann, A. (2007) In Vivo Measurement of Cortical Impedance Spectrum in Monkeys: Implications for Signal Propagation. Neuron, 55, 809-823.

https://doi.org/10.1016/j.neuron.2007.07.027

[3] Gomes, J.M., Bédard, C., Valtcheva, S., Nelson, M., Khokhlova, V., Pouget, P., Venance, L., Bal, T. and Destexhe, A. (2016) Intracellular Impedance Measurements Reveal Non-Ohmic Properties of the Extracellular Medium around Neurons. Biophysical Journal, 110, 234-246.

https://doi.org/10.1016/j.bpj.2015.11.019

[4] Miceli, S., Ness, T.V., Einevoll, G.T. and Schubert, D. (2017) Impedance S pectrum in Cortical Tissue: Implications for Propagation of LFP Signals on the Microscopic Level. eNeuro, 4, e0291.

[5] Bédard, C., Gomes, J.-M., Bal, T. and Destexhe, A. (2017) A Framework to Reconcile Frequency Scaling Measurements, from intracellular Recordings, Local-Field Potentials, up to EEG and MEG Signals. Journal of Integrative Neuroscience, 16, 3-18.

https://doi.org/10.3233/JIN-160001

[6] Bédard, C. and Destexhe, A. (2009) Macroscopic Models of Local Field Potentials and the Apparent 1/f Noise in Brain Activity. Biophysical Journal, 96, 2589-2603.

https://doi.org/10.1016/j.bpj.2008.12.3951

[7] Bédard, C. and Destexhe, A. (2011) A Generalized Theory for Current-Source Density Analysis in Brain Tissue. Physical Review E, 84, 041909.

https://doi.org/10.1103/PhysRevE.84.041909

[8] Landau, L.D. and Lifshitz, E.M. (1981) Electrodynamics of Continuous Media. Pergamon Press, Moscow, Russia.

[9] Kronig, R.D.L. (1926) On the Theory of Dispersion of X-Rays. Journal of the Optical Society of America, 12, 547.

https://doi.org/10.1364/JOSA.12.000547

[10] Foster, K.R. and Schwan, H.P. (1989) Dielectric Properties of Tissues and Biological Materials: A Critical Review. Crit. Reviews Biomed. Engineering, 17, 25-104.

[11] Appel, W. (2007) Mathematics for Physics and Physisics. Princeton University Press, Princeton, NJ.

[12] Schönleber, M., Klotz, D. and Ivers-Tiffée, E. (2014) A Method for Improving the Robustness of linear Kramers-Kronig Validity Tests. Electrochimica Acta, 131, 20-27.

https://doi.org/10.1016/j.electacta.2014.01.034

[13] Le Bellac, M. and Levy-Leblond, J.-M. (1973) Galilean Electromagnetism. Nuevo Cimento, 14B, N.2

[14] Rousseaux, G. (2005) Lorenz or Coulomb in Galilean Electromagnetism? Europhysics Letters, 71, 15-20.

https://doi.org/10.1209/epl/i2005-10059-5

[15] Bedard, C. and Destexhe, A. (2012) Modeling Local Field Potentials and Their Interaction with the Extracellular Medium. In: Brette R. and Destexhe A., Eds., Handbook of Neural Activity Measurement, Cambridge University Press, Cambridge, 136-191.

[16] Bédard, C. and Destexhe, A. (2013) Generalized Cable Theory for Neurons in Complex and Heterogeneous Media. Physical Review E, 88, 022709.

https://doi.org/10.1103/PhysRevE.88.022709

[17] Bédard, C. and Destexhe, A. (2014) Generalized Cable Formalism to Calculate the Magnetic Field of Single Neurons and Neuronal Populations. Physical Review E, 90, 042723.

https://doi.org/10.1103/PhysRevE.90.042723

[18] Pham, P., Roux, S., Matonti, F., Dupont, F., Agache, V. and Chavane, F. (2013) Post-Implantation Impedance Spectroscopy of Subretinal Microelectrode Arrays, OCT Imaging and Numerical Simulation: Towards a More Precise Neuroprosthesis Monitoring Tool. Journal of Neural Engineering, 10, 046002.

https://doi.org/10.1088/1741-2560/10/4/046002

[19] Hines, M.L. and Carnevale, N.T. (1997) The NEURON Simulation Environment. Neural Computation, 9, 1179-1209.

https://doi.org/10.1162/neco.1997.9.6.1179

[20] Bédard, C. and Destexhe, A. (2017) Is the Extracellular Impedance High and Non-Resistive in Cerebral Cortex? Biophysical Journal, 113, 1639-1642.

https://doi.org/10.1016/j.bpj.2017.08.021