Abstract

The electrical properties of biological tissues can be described by a complex tensor comprising a simple expression of the effective admittivity. The effective admittivities of biological tissues depend on scale, applied frequency, proportions of extra- and intracellular fluids, and membrane structures. The effective admittivity spectra of biological tissue can be used as a means of characterizing tissue structural information relating to the biological cell suspensions, and therefore measuring the frequency-dependent effective conductivity is important for understanding tissue’s physiological conditions and structure. Although the concept of effective admittivity has been used widely, it seems that its precise definition has been overlooked. We consider how we can determine the effective admittivity for a cube-shaped object with several different biologically relevant compositions. These precise definitions of effective admittivity may suggest the ways of measuring it from boundary current and voltage data. As in the homogenization theory, the effective admittivity can be computed from pointwise admittivity by solving Maxwell equations. We compute the effective admittivity of simple models as a function of frequency to obtain Maxwell-Wagner interface effects and Debye relaxation starting from mathematical formulations. Finally, layer potentials are used to obtain the Maxwell-Wagner-Fricke expression for a dilute suspension of ellipses and membrane-covered spheres.

1. Introduction

The human body can be regarded as a complex electrical conductor comprising many tissues that have distinct electrical properties. Measurements of the electrical properties of biological tissues have shown that effective conductivity () and permittivity () values of biological tissues in the frequency range from a few Hz to MHz are influenced by physiological and pathological conditions [15]. The effective admittivity of a biological tissue under the influence of a time-harmonic electric field at an angular frequency is determined by its ion concentrations in extra- and intracellular fluids, cellular structure and density, molecular compositions, membrane characteristics, and other factors. Cell membranes contribute to capacitance; the intracellular fluid gives rise in an intracellular resistance; the extracellular fluid contributes to effective resistance. As a result, biological tissues show a variable response over the frequency range from a few Hz to MHz. For most biological tissues, at low frequencies below 10 kHz, whereas the term is not negligible beyond 10 kHz due to the abundant membranous structures in organisms.

The effective admittivity can be regarded as a function of tissue composition and the applied angular frequency . Assume that a biological subject under consideration is a mixture of homogeneous tissue at macroscopic length scale and has a constant effective admittivity in a particular cubic sample voxel. The can be viewed as the effective tensor according to the well-known concept of homogenization when the admittivity is periodic [6]. The effective admittivity as a function of and the voxel can be determined by Ohm's law: for a time-harmonic electric field and the corresponding current density at angular frequency . Here, is a symmetric matrix and is the position within the voxel. If the quantity for any pair of electric field and current density , then is scalar and the subject is isotropic. Otherwise, a subject is said to be anisotropic. Depending on the measurement scale used, anisotropy may or may not be detected. There have been numerous studies and models formulated for admittivity spectra of biological tissue as a means of characterizing tissue structural information relating to biological cell suspensions [7, 8]. In 1873, Maxwell [9] derived an expression of for the special case of a strongly dilute suspension of spherical particles and . Wagner extended the expression to a general . Poisson [10] in 1826 and Faraday in 1827 dealt with the case of a suspension of infinitely conducting spheres in a background. In 1924, Fricke [7] provided an expression for by considering the capacity due to a polarization at the interphases or the presence in the interphases of thin poorly conducting membranes. There have been also studies on an effective conductivity of cell suspensions, both analytically and numerically [1117].

The concept of effective admittivity has been used widely, but it seems that its precise definition has been overlooked. How can we determine the effective admittivity of a given cubic region? In this paper we give precise definitions of effective admittivity to provide a way to measure it from boundary current and voltage data. As in the homogenization theory, the effective admittivity can be computed from pointwise admittivity by solving Maxwell equations. We compute the effective admittivity of simple models as a function of frequency to observe the Maxwell-Wagner interface effect and Debye relaxation using a mathematical point of view. Single layer and double layer potentials are used to produce the Maxwell-Wagner-Fricke expression for a dilute suspension of ellipses and a membrane-covered conductor, respectively. We also note that Maxwell equations make both microscopic and macroscopic senses.

2. Effective Admittivity Spectra of Biological Tissues

The concept of admittivity contains four key definitions: pointwise admittivity, effective admittivity, apparent admittivity, and equivalent admittivity. (i)Pointwise admittivity refers to electrical properties at microscopic scale. (ii)Effective admittivity is defined at macroscopic scale. Homogenization methods have been used to compute effective property of a periodic heterogeneous subject from its pointwise structure [18]. It is used to describe the linear relationship between the ensemble mean current density and the ensemble mean electrical field. Effective admittivity depends only on the electrical properties of the sample. (iii)Apparent admittivity is defined as the admittivity of electrically homogeneous and isotropic medium that could yield the potential measured on the heterogeneous subject using the same applied current and arrangement of the electrodes. (iv)Two expressions that have the same effective admittivity are called equivalent admittivity. Pavlin and Miklavčič [19] use a simpler equivalent conductivity of a single cell for the purpose of computation of effective conductivity of a suspension of permeabilized cells.Let be a three-dimensional domain with a pointwise admittivity of , where the conductivity and the permittivity values are assumed to depend only on position , and both are isotropic. Then, the domain can be viewed as a union of many voxels , and the effective properties mainly depend on the choice of voxels. With a given voxels, we can define the effective admittivity that is a constant on each voxel . The effective admittivity is a tensor-valued function of the voxel and the angular frequency such that The must be the best approximation of the average of the pointwise admittivity over the voxel in the sense that The physically meaningful solution must have a finite energy [20]: Hence, the solution of the equation should be contained in the set [20]. Here, is the standard Sobolev space equipped with norm .

However, there is no such tensor satisfying (3) exactly. Hence, we may take an appropriate satisfying (3) approximately, and the choice of may differ for a biological sample. To clearly define the effective admittivity, we need to select suitable potentials satisfying in .

2.1. Definition of Effective Admittivity for a Cubic Sample

Let us consider a rectangular-shaped tissue sample (Figure 1) occupied in the unit cube with its three pairs of facing surfaces (Figure 2): Assume that the admittivity distribution of the sample at frequency is given by , where the conductivity and the permittivity values are scalar-valued functions depending only on position . Hence, is isotropic on a microscopic scale. If we apply a current of through the pair of electrodes attached on and , then the resulting time-harmonic potential satisfies the following equation from a suitable arrangement of Maxwell equations (at frequencies below about 100 kHz): where is the unit outward normal vector on . For each , we denote the voltage difference by

Lemma 1 (reciprocity). For , one has

Proof. From the boundary conditions of and divergence theorem, we have Hence, the symmetry (8) follows from the reciprocity relation
If (homogeneous), then and and must be . If the effective admittivity is a diagonal matrix satisfying then it must be , and due to the following theorem.

Theorem 2. If is the solution of (6), then one has

Proof. We will only  prove  . From the definition of , , and the divergence theorem, we have Since and , This completes the proof of (12).

Now, we are ready to define the effective admittivity tensor .

Definition 3. For a given unit cubic and each , let be the potential difference given in (7). Then the effective admittivity tensor is defined by

The proposed definition may not have coordinate invariance due to its limitation of the tensor expression. For a proper invariance, we need to compute all the tensors (15) by rotating the coordinate system. We may define the effective admittivity tensor as the best fit of the minimization problem described in (3).

Next, we study how the distribution of is related to the frequency-dependent behavior of .

2.2. One Dimensional Sample

We begin by considering a special sample (Figure 3) with depending only on the -variable: where , and are constants and . For this sample, the potential in (6) depends only on -variable, and Since is a constant, Writing , we have Hence, it follows from the definition of (15) of that This means that is the harmonic average of the admittivity that can be expressed as From this, we have where Writing and , we have Here, is referred to as a relaxation time, since its value controls polarization time [8, 21]. It is remarkable to observe that the relaxation time may be obtained by solving the elliptic PDE (17).

Using (23), the average current density generated inside the dielectric due to the average electric field is given by which can be expressed as where and Here, and are referred to as the dielectric constant and loss factor of the dielectric material, respectively. The average current density can also be written as

In biological materials is produced by the ionic conduction and is produced by dielectric relaxation.

The dielectric response of biological tissues is always frequency dependent, and the electric charge movement inside the material in response to an externally applied electric field is controlled by the dielectric properties of the material. The free charge movement inside a material affected by an external field is controlled by its conductivity . Biological tissues display extremely high dielectric constants at low frequencies, and as the excitation frequency is increased, the dielectric constants of the tissues fall off in more or less distinct steps [8]. Interfaces play a significant role in the frequency dependence of complex materials, particularly at audio and subaudio frequencies [8]. The frequency response of biological tissue admittivity is highly influenced by the dielectric polarization, dielectric relaxation, and dielectric dispersion.

Electric polarization (Figure 4) may be defined as the electric-field-induced disturbance (shift from average equilibrium positions) of the charge distribution in a region [8]. Dielectric dispersion in biological tissues can be assumed to depend upon the permittivity (Figure 5) of tissue material with applied electric field frequency [8]. In other words, a significant change in dielectric properties over a frequency range, by convention, is called a dielectric dispersion [21].

As there is always a lag between the changes in an applied electric field and changes in polarization, the permittivity of the biological tissues is a complex-valued function of the frequency of the applied electric field. The term dielectric relaxation [22] in a biological tissue connotes the delay or lag in its response to create the dielectric polarization following the application of electric field across the tissue sample. In other words, the dielectric relaxation of a tissue can be defined as the lag (momentary delay) in the dielectric constant which is usually caused by the delay in molecular polarization with respect to a change in applied electric field. According to the previous simple computations of (23) in the 1D model (17), the central frequency of the dispersion is .

Schwan [23, 24] studied the properties of biological tissue and cell suspensions over a large frequency range and observed that the dielectric properties of biological tissues are characterized by three major dispersions, -dispersion [23, 24], -dispersion [23, 24], and -dispersion [23, 24] occurring at low frequency, radio frequency, and microwave frequency, respectively. We consider each of these dispersions below.(i)-dispersion : The -dispersion is associated with tissue interfaces such as membranes [23]. Below about 10 kHz, the dielectric studies of biological or any other electrolyte systems become very complex and difficult to characterize. Foster and Schwan, 1989 [25], reported that -dispersion is believed to be associated with a counterion layer (electrical double layer) polarization in tissues. (ii)-dispersion : In biological tissues, the -dispersion is caused by the polarization of cellular membranes and polarization of protein and other organic macromolecules [23]. The -dispersion arises, principally, from interfacial polarization (Maxwell-Wagner effect) [26] of cell membranes [21]. In the frequency range , the dielectric behavior of the tissues is dominated by the heterogeneous composition and ionic activities inside the biological tissue. These effects are principally responsible for the -dispersion. The radio frequency dispersion or -dispersion has been recognized as a Maxwell-Wagner relaxation [26] caused by cell membranes [27]. A large magnitude, low frequency -dispersion was observed by Schwan [28] in a muscle tissue. This effect is related, in part at least, to the tubular shape of muscle fibers [29]. The theoretical aspects of the low frequency dispersion of colloid particles in electrolyte solution have been studied by Schwarz in 1962 [30].(iii)-dispersion . The -dispersion in biological tissues is caused by the reorientation of water molecules [18]. This dispersion has been well studied and has found many applications [3133]. Rajewsky and Schwan [34] noted the -dispersion at microwave frequencies which is understood to be caused by abundant tissue water. Schwan conducted the extensive studies on the electrical properties of biological cell suspensions in 1993 [35] over a broad frequency range extending from less than 1 Hz to many GHz and summarized the mechanisms which contribute to the total frequency response. He studied the mechanisms responsible for electrical properties of tissues and cell suspensions, and he observed that the frequency changes of these properties obey causality, that is, the Kramers-Kronig relationships [35] which relate changes of dielectric constants to conductivity changes. A number of mechanisms which reflect the various compartments of the biological materials were identified such as membranes and their properties, biological macromolecules, and fluid compartments inside and outside membranes.

Membrane relaxation is anticipated from the Hodgkin-Huxley membrane model [36] and adds to the -effects [8], and hence a number of -effects of small magnitude occur at the tail of the -dispersion caused by proteins, proteinbound water (called -dispersion), and cell organelles such as mitochondria [37]. A second Maxwell-Wagner dispersion [26] which occurs at frequencies well above those of the main -dispersion [8] is a characteristic of suspended particles surrounded by a shell and usually of small magnitude [14].

Figure 6(a) shows Cole-Cole plot explaining . Figure 6(b) shows versus line describing . Figure 6(c) shows versus line.

Remark 4. In the case when is sufficiently small (dilute suspension) so that , (21) can be expressed as because . Neglecting in (30), we get In three dimensional heterogeneous medium, this type of dilute suspension model with neglecting had been used in computation of the effective admittivity .

Next, we will investigate the effective admittivity for dilute suspensions of membrane of materials. We will express potential of models comprising suspension of arbitrary-shaped membrane, using double layer potential technique.

2.3. Dilute Single Suspension of Ellipses in a Cube

Maxwell [9] and Wagner [38] analyzed expressions for the effective admittivity of a strongly diluted suspension of spheres [7].

Let   be a cube, and let be an ellipsoid with . As in the previous section, let the admittivity distribution (Figure 7) be given by If is a potential satisfying in , then it can be expressed as a sum of harmonic function in and a single layer potential: where is determined by Assuming that the volume fraction is small, Fricke obtained the following approximation: where

When ( is sphere) and the volume fraction is small, Maxwell-Wagner [26] formula for is given by With the aid of , we can derive the Debye dispersion function for as follows:

Next, we will investigate the effective admittivity for dilute suspensions of membranes.

2.4. Dilute Single Suspension of Membrane

Finally, consider the case of a dilute single suspension of membranes (Figure 8) to investigate the role of the thin insulating membrane influencing the frequency-dependent behavior of the effective admittivity. In the model comprising a suspended membrane, there exists a thin membrane of a thickness , as shown in Figure 6, within our target voxel. Assume that the admittivity distribution changes abruptly across the membrane

In the case of a dilute suspension of a single membrane, any potential satisfying can be expressed as where is the surface of the membrane and The density is determined by membrane structure and the refraction index: Here, we recall on the surface .

In 1955, Fricke [39] studied the equivalent admittivity for the case of a spherical membrane suspension using its pointwise admittivity distribution (Figure 7):

Fricke's expression for the equivalent admittivity for the internal domain including the membrane and is Then the effective admittivity can be computed by substituting this equivalent admittivity in Maxwell-Wagner-Fricke formula [40]. Under the assumption that three quantities , and are very small, Pauly and Schwan [14] obtained where

Biological tissues are comprised of cells and an extracellular matrix of macromolecules and extracellular fluid. Cells have different shapes and sizes in the order of m. The cell has a membrane whose thickness is in the order of several nm. Enclosed within the cell membrane is the cytoplasm containing the nucleus, organelles, and intracellular fluid. Although they are very thin and resistive, cell membranes play an important role in determining admittivity values of tissues. They usually account for large susceptivity values, even though there is little quantitative analysis on how membrane properties affect them.

Membrane suspensions in homogeneous media can be viewed as an equivalent biological tissue admittivity. Analysis of the simple membrane structures should enable us to interpret admittivity images from phantom experiments and numerical simulations.

3. Discussion and Conclusions

Tomographic imaging of the admittivity distributions inside biological subjects such as the human body has been an active research goal in electrical impedance tomography (EIT). Can EIT technique distinguish between cucumber and carrot or lung and liver? Measuring the frequency-dependent behavior of effective admittivity increases distinguishability and has a potential of expanding clinical applications. The effective conductivity of biological tissue is associated with the forward problem of an elliptic PDE, the Laplace equation, with a complex material parameters and thin insulating membranes, and there is little understanding about this kind of forward problem. Since traditional electrical impedance tomography can only provide change in effective admittivity on a coarse scale due to its low spatial resolution, we will have a better understanding of EIT images if we can link microscopic cell structures to macroscopic (or effective) tissue admittivity images obtained from multifrequency EIT systems.

We defined the effective admittivity of a cubic voxel from its pointwise admittivity and the Maxwell equations. Because of the reciprocity property, the effective admittivity must be symmetric. Direct computation of effective admittivity using simple models produced the Maxwell-Wagner interface effect and Debye relaxation. The complex permittivity, dielectric constant, and loss factors were derived and modeled as the functions of frequencies, and the dielectric polarization, dispersions, and the relaxation time were studied in a mathematical framework. The single and double layer potentials were used to present a mathematical expression for the Maxwell-Wagner-Fricke expression for the subject containing various geometries of cells or membranes while previous work of Maxwell and Wagner only handles ellipses in a cube.

Conflict of Interests

All the authors declare that they have no conflict of interests.

Acknowledgment

The National Research Foundation of Korea (NRF) Grant funded by the Korea Government (MEST) No. 2011-0028868 (2012R1A2A1A03670512).