- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Annual Issues
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
Computational and Mathematical Methods in Medicine
Volume 2013 (2013), Article ID 353849, 10 pages
Effective Admittivity of Biological Tissues as a Coefficient of Elliptic PDE
1Department of Computational Science and Engineering, Advanced Science and Technology Center (ASTC), Yonsei University, 50 Yonsei-Ro, 134 Sinchon-dong, Seodaemun-gu, Seoul 120 749, Republic of Korea
2J. Crayton Pruitt Family Department of Biomedical Engineering, University of Florida, Biomedical Sciences Building JG-5, P.O. Box 116131, Gainesville, FL 32611, USA
Received 26 October 2012; Accepted 15 January 2013
Academic Editor: Eung Je Woo
Copyright © 2013 Jin Keun Seo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
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.
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 [1–5]. 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 . 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  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  in 1826 and Faraday in 1827 dealt with the case of a suspension of infinitely conducting spheres in a background. In 1924, Fricke  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 [11–17].
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 . 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č  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 : Hence, the solution of the equation should be contained in the set . 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 . Interfaces play a significant role in the frequency dependence of complex materials, particularly at audio and subaudio frequencies . 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 . Dielectric dispersion in biological tissues can be assumed to depend upon the permittivity (Figure 5) of tissue material with applied electric field frequency . In other words, a significant change in dielectric properties over a frequency range, by convention, is called a dielectric dispersion .
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  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 . 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 , 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 . The -dispersion arises, principally, from interfacial polarization (Maxwell-Wagner effect)  of cell membranes . 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  caused by cell membranes . A large magnitude, low frequency -dispersion was observed by Schwan  in a muscle tissue. This effect is related, in part at least, to the tubular shape of muscle fibers . The theoretical aspects of the low frequency dispersion of colloid particles in electrolyte solution have been studied by Schwarz in 1962 .(iii)-dispersion . The -dispersion in biological tissues is caused by the reorientation of water molecules . This dispersion has been well studied and has found many applications [31–33]. Rajewsky and Schwan  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  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  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  and adds to the -effects , 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 . A second Maxwell-Wagner dispersion  which occurs at frequencies well above those of the main -dispersion  is a characteristic of suspended particles surrounded by a shell and usually of small magnitude .
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
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  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 .
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 . Under the assumption that three quantities , and are very small, Pauly and Schwan  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.
The National Research Foundation of Korea (NRF) Grant funded by the Korea Government (MEST) No. 2011-0028868 (2012R1A2A1A03670512).
- L. A. Geddes and L. E. Baker, “The specific resistance of biological material-a compendium of data for the biomedical engineer and physiologist,” Medical and Biological Engineering, vol. 5, no. 3, pp. 271–293, 1967.
- S. Gabriel, R. W. Lau, and C. Gabriel, “The dielectric properties of biological tissues: II. Measurements in the frequency range 10 Hz to 20 GHz,” Physics in Medicine and Biology, vol. 41, no. 11, pp. 2251–2269, 1996.
- C. Gabriel, S. Gabriel, and E. Corthout, “The dielectric properties of biological tissues: I. Literature survey,” Physics in Medicine and Biology, vol. 41, no. 11, pp. 2231–2249, 1996.
- S. Gabriel, R. W. Lau, and C. Gabriel, “The dielectric properties of biological tissues: III. Parametric models for the dielectric spectrum of tissues,” Physics in Medicine and Biology, vol. 41, no. 11, pp. 2271–2293, 1996.
- S. Grimnes and G. Martinsen, Bioimpedance and Bioelectricity Basics, Academic Press, San Diego, Calif, USA, 2008.
- A. Bensoussan, J. L. Lions, and G. Papanicolaou, Asymptotic Analysis for Periodic Structures, Elsevier, 1978.
- H. Fricke, “A mathematical treatment of the electrical conductivity of colloids and cell suspensions,” Journal of General Physiology, vol. 4, pp. 375–383, 1924.
- O. G. Martinsen, S. Grimnes, and H. P. Schwan, “Interface phenomena and dielectric properties of biological tissue,” Encyclopedia of Surface and Colloid Science, vol. 20, pp. 2643–2653, 2002.
- J. C. Maxwell, Electricity and Magnetism, 1st edition, 1873.
- S. D. Poisson, Mémoires De L'Académie RoyaLe Des Sciences De L'Institut De France, vol. 5, p. 488, 1826.
- S. S. Dukhin, “Dielectric properties of disperse systems,” in Surface and Colloid Science, E. Matijevic, Ed., vol. 3, pp. 83–165, Wiley-Interscience, New York, NY, USA, 1971.
- M. Pavlin, T. Slivnik, and D. Miklavčič, “Effective conductivity of cell suspensions,” IEEE Transactions on Biomedical Engineering, vol. 49, no. 1, pp. 77–80, 2002.
- M. Pavlin, N. Pavšelj, and D. Miklavčič, “Dependence of induced transmembrane potential on cell density, arrangement, and cell position inside a cell system,” IEEE Transactions on Biomedical Engineering, vol. 49, no. 6, pp. 605–612, 2002.
- H. Pauly and H. P. Schwan, “Über die impedanz einer suspension von kugelformigen teilchen mit einer schale,” Zeitschrift für Naturforschung B, vol. 14, pp. 125–131, 1959.
- H. Ammari, H. Kang, and K. Touibi, “Boundary layer techniques for deriving the effective properties of composite materials,” Asymptotic Analysis, vol. 41, no. 2, pp. 119–140, 2005.
- H. Ammari and H. Kang, Polarization and Moment Tensors: With Applications to Inverse Problems and Effective Medium Theory, vol. 162, Springer, 2007.
- H. Ammari, P. Garapon, H. Kang, and H. Lee, “Effective viscosity properties of dilute suspensions of arbitrarily shaped particles,” The Journal of Physical Chemistry, vol. 80, pp. 189–211, 2012.
- W. Graeme, Milton the Theory of Composites, Cambridge University Press.
- M. Pavlin and D. Miklavčič, “Effective conductivity of a suspension of permeabilized cells: A theoretical analysis,” Biophysical Journal, vol. 85, no. 2, pp. 719–729, 2003.
- J. K. Seo and E. J. Woo, Nonlinear Inverse Problems in Imaging, John Wiley & Sons, 2012.
- W. Kuang and S. O. Nelson, “Low-frequency dielectric properties of biological tissues: a review with some new insights,” Transactions of the American Society of Agricultural Engineers, vol. 41, no. 1, pp. 173–184, 1998.
- W. E. Vaughan, “Dielectric relaxation,” Annual Review of Physical Chemistry, vol. 30, pp. 103–124, 1979.
- H. P. Schwan, “Electrical properties of tissue and cell suspensions,” in Advances in Biological and Medical Physics, J. H. Lawrence and C. A. Tobias, Eds., vol. V, pp. 147–209, Academic Press, New York, NY, USA, 1957.
- H. P. Schwan, “Electrical properties of tissues and cell suspensions: mechanisms and models,” in Proceedings of the 16th Annual International Conference of the IEEE, vol. 1, pp. A70–A71, 1994.
- K. R. Foster and H. P. Schwan, “Dielectric properties of tissues and biological materials: a critical review.,” Critical Reviews in Biomedical Engineering, vol. 17, no. 1, pp. 25–104, 1989.
- T. Hao, Electrorheological Fluids: The Non-Aqueous Suspensions, Elsevier, 2005.
- K. S. Cole, Membranes, Ions and Impulses, University of California Press, Berkeley, Calif, USA.
- H. P. Schwan, “Electrical properties of muscle tissue at low frequencies,” Zeitschrift für Naturforschung B, vol. 9, p. 245, 1954.
- G. Falkand and P. Fatt, “Linear electrical properties of striated muscle fibers observed with intracellular electrodes,” Proceedings of the Royal Society of London B, vol. 160, pp. 69–123, 1964.
- G. Schwarz, “A theory of the low-frequency dielectric dispersion of colloidal particles in electrolyte solution,” Journal of Physical Chemistry, vol. 66, no. 12, pp. 2636–2642, 1962.
- J. Thuery, Microwaves: Industrial, Scientific and Medical Applications, Artech House, Boston, Mass, USA, 1991.
- S. O. Nelson, “Dielectric properties of agricultural products measurements and applications,” IEEE Transactions on Electrical Insulation, vol. 26, no. 5, pp. 845–869, 1991.
- A. C. Metaxas and R. J. Meredith, Industrial Microwave Heating, Peter Peregrinus, London, UK, 1983.
- B. Rajewsky and H. P. Schwan, “The dielectric constant and conductivity of blood at ultrahigh frequencies,” Naturwissenschaften, vol. 35, p. 315, 1948.
- H. P. Schwan, “Mechanisms responsible for electrical properties of tissues and cell suspensions,” Medical Progress through Technology, vol. 19, no. 4, pp. 163–165, 1993.
- H. M. Fishman, D. J. M. Poussart, L. E. Moore, and E. Siebenga, “K+ conduction description from the low frequency impedance and admittance of squid axon,” Journal of Membrane Biology, vol. 32, no. 3-4, pp. 255–290, 1977.
- R. D. Stoy, K. R. Foster, and H. P. Schwan, “Dielectric properties of mammalian tissues from 0.1 to 100 MHz; a summary of recent data,” Physics in Medicine and Biology, vol. 27, no. 4, pp. 501–513, 1982.
- K. W. Wagner, “Erklarung der dielectrischen nachwirkungsworgange auf grund maxwellscher vorstellungen,” Archiv für Elektrotechnik, vol. 2, pp. 371–387, 1914.
- H. Fricke, “The complex conductivity of a suspension of stratified particles of spherical or cylindrical form,” Journal of Physical Chemistry, vol. 59, no. 2, pp. 168–170, 1955.
- D. Gakkai, “Japan society for medical electronics and biological engineering,” in Proceedings of 5th International Conference on Electrical Bio-Impedance, p. 261, Business Center for Academic Societies Japan, 1888.