#### Abstract

*Background*. Action currents in the heart produce a magnetic field, which could provide a way to detect the propagation of electrical activity through cardiac tissue using magnetic resonance imaging. However, the magnetic field produced by current in the heart is small. The key question addressed in this study is are cardiac biomagnetic fields large enough to be detectable by MRI?* Results*. A spherical model is used to calculate the magnetic field inside the heart, which has a magnitude of about 14 nT. This field implies a phase shift in the MRI signal of about 0.2°.* Conclusion*. Phase shifts associated with cardiac action currents will be difficult to detect using current MRI technology but may be possible if motion artifacts and other physiological noise can be suppressed.

#### 1. Introduction

Action currents in the heart produce a magnetic field, which could create an artifact in a magnetic resonance image (MRI). Viewed from another perspective, in theory MRI could be used to image action currents by detecting their magnetic field. Indeed, if one could image action currents using MRI, it would be a fundamental advance that might make magnetocardiography and even much of electrocardiography obsolete. However, the magnetic field produced by current in the heart is small. The key question addressed in this study is are biomagnetic fields large enough to be detectable by MRI? The main goal of this paper is to compute how large this magnetic field is and explore the possibility of detecting the magnetic field by magnetic resonance imaging.

Over the last ten years, much work has been done to detect neuronal activation by measuring the magnetic field produced by biocurrents using MRI [1–8]. In general, these magnetic fields are so tiny that they are difficult to measure [3]. However, the heart is the source of the largest biocurrents in the body, so it may be easier to detect cardiac action currents using MRI compared to neural action currents [9]. In order to interpret and predict the MRI signal produced by cardiac action currents, we must know the magnetic field distribution within the heart.

In general, the best starting place for calculating the magnetic field is a simple model that clearly illustrates mechanisms. In this study, we examine a “spherical heart” model [10] with a simple transmembrane potential distribution and calculate the resulting action currents and magnetic field. A simple model allows us to obtain analytical solutions that provide insight into the qualitative behavior and allows us to estimate the quantitative magnitude of the magnetic field.

There are two sources that produce a magnetic field in cardiac tissue: one is the intracellular current in the tissue with the “return” loops through an adjacent volume conductor, and the other is the anisotropy of the tissue [11]. In this study we account for the first of these mechanisms: the heart is surrounded by a volume conductor that provides a path for return currents and results in a magnetic field.

#### 2. Methods

We consider a spherical shell of cardiac tissue surrounding a blood cavity and surrounded by an unbounded conducting bath [10] (Figure 1). We use spherical coordinates (, , and ) to specify position. The inner and outer radii of the shell are and , respectively. The blood (subscript “”) conductivity is and the bath (subscript “”) conductivity is .

We assume the transmembrane potential has azimuthally symmetry. Therefore, it does not depend on the angle around the heart. With an analogy to earth, the transmembrane potential only depends on latitude but not longitude. The fiber geometry in the heart is complex, and we must simplify it significantly in order to obtain an analytical solution. Here the fiber axis is parallel to the meridians, that is, along lines of constant longitude (lie along direction) [10]. With this fiber model, we can calculate analytically the most significant information that we want, that is, the magnitude of the magnetic field .

We use the bidomain model to represent the electrical properties of cardiac tissue, which accounts for anisotropy both inside and outside the myocardial cells [12]. The bidomain model is a continuum model, in the sense that it represents the electrical properties of the tissue averaged over many cells [13]. The electrical conductivities in cardiac tissue (, , , and ) are different in and directions (perpendicular and parallel to the fibers) and in the intracellular (subscript “”) and extracellular (subscript “”) spaces.

We assume the tissue is quasistatic: the governing equations do not depend on time. The transmembrane potential will depend on time as the action potential propagates through the tissue, but, given the transmembrane potential at any moment, the other potentials, currents, and magnetic field can be calculated at that instant quasistatically.

##### 2.1. Calculation of the Electrical Potential

In the cardiac shell, , there are two potential fields, and , which are functions of and . The transmembrane potential is the difference of the intracellular and extracellular potentials:The fundamental equation governing the current density is the equation of continuity Since we only require the sum of and be divergenceless, and individually need not have vanishing divergence. Thus, current can flow from one domain to another.

Using Ohm’s law and the relationship between the potential and the electrical field , we find that , , , and , with and . If we use the above four current density expressions, then (2)–(4) become In order to simplify (6), we do a linear transformation from and to and [10, 14]:with the inverse transformation given byIf we substitute (9) and (10) into the expressions in (6), we get If we make the additional assumption that is not a function of , the second term in the square brackets of (11) does not contribute to (6), and we obtainwhere

##### 2.2. Boundary Conditions

The behavior at the interface between cardiac tissue and an adjacent volume conductor is determined by boundary conditions [15, 16]. At the surfaces of the shell, and , we require that the normal component of the current density is continuous:and the extracellular potential is continuousThe boundary conditions in (15)–(17) give at For the boundary , we replace by in (18).

##### 2.3. Solution of the Equations for the Potential

If we expand in terms of Legendre polynomials [17], , where , then the general solution to (13) is [10]whereThe general solutions to Laplace’s equation in the blood and bath areThe boundary conditions determine the constants , , , and for each . We expand in terms of Legendre polynomialswhereSolving the four linear equations from the boundary conditions, we obtain the constants , , , and for each :where

##### 2.4. Calculation of the Magnetic Field

Using Ohm’s law, we obtain the current density distributionAmpere’s law in its integral form iswhere is the current enclosed by the Amperian loop. We choose this loop to be the edge of a right circular cone with the apex at the center of the heart sphere and the integral surface is the surface of the right circular cone: The magnetic field is in direction, tangential to this loop. By varying and of this cone, we can get the magnetic field in all the spacewhere In order to compute and thereby the magnitude of** B**, we need to substitute the potential ((19), (21), and (22)) into the expression of current density (30), and then replace in (34) by these expressions of current density in each region , , and , respectively. We then sum up the integrals in a piecewise manner for the three spaces in blood, bath, and tissue.

In the blood , (note: ).

In the bath , Since the total current integrated over all the space is zero, we deduceIn the tissue The integral interval of (34) is split into two pieces: is and is . The total enclosed current is then Using (27) and (28), we getAn advantage of this formulation is that by varying and the magnetic field is determined everywhere including within the heart wall.

#### 3. Results

To perform calculations of the magnetic field, we need to consider a specific transmembrane potential distribution and values for the model parameters. We set the transmembrane potential to be = = 20 mV for *θ* < 89° and = = −80 mV for *θ* > 91°. For the interval 89° < *θ* < 91°, decreases linearly from 20 to −80 mV. The heart size is = 30 and = 40 mm [10]. Furthermore, we set the blood and the bath conductivities to = = 1 S/m. The electrical conductivities of cardiac tissue are anisotropic: = 0.02, = 0.2, = 0.08, and = 0.2 S/m [18].

Using the given transmembrane potential, we need to calculate its expansion in terms of Legendre polynomials.

##### 3.1. The Calculation of Transmembrane Potential Expansion

To simplify the calculation, we change variables by letting . We let , , , and , and then (24) becomes For these integral intervals, ① refers to , ② refers to , and ③ refers to . ① isUsing the equation of a line, , we get . Therefore, ② isWe need to obtain and :

③ isThe final result is the sum of integrals of the three segments. Figure 2 is a comparison of the calculated approximate transmembrane potential with the original transmembrane potential, showing excellent agreement.

**(a)**

**(b)**

Figure 3 contains plots of the extracellular potential, calculated using (10), (16), and (17) with (19), (21), and (22). The potential falls off rapidly with distance into the bath or blood but is relatively large inside the cardiac tissue.

##### 3.2. The Upper Bound of l of the Approximated Sums

The magnetic field is calculated in terms of Legendre polynomials. The infinite sums given above ( from zero to infinity) must be approximated as finite sums. Figure 4 shows the change of the calculated peak value of the magnetic field () as the total number of terms, , in the sum increases. changes by less than 3% for > 200, and therefore we use an upper bound = 200, which is enough to determine the magnitude of accurately and is still easy to compute. When we choose* l* = 200, the point with the maximum magnitude, , is about 14 nT.

Figure 5 displays the -component of the magnetic field over a cross section of the heart. The field is largest near the inner and outer surfaces of the tissue and is relatively small throughout much of the heart wall. It has a peak magnitude, , of about 14 nT. Figure 6 shows the total current density (intracellular plus extracellular).

#### 4. Discussion

The calculation of the magnetic field throughout the heart indicates that the peak magnetic field has a magnitude of about 14 nT. This is a similar order of magnitude as calculated by Roth [19] for a cylindrical strand of cardiac tissue, such as a papillary muscle, and is several times larger than the magnetic field at the apex of the heart (2 nT) [20, 21], in planar tissue samples (1 nT) [22], or outside the body during measurements of the magnetocardiogram (0.05 nT) [23]. The magnetic field is much larger at the tissue surface than just a few millimeters away from the surface, either in the surrounding bath or deeper in the tissue, which is again similar to the results found in a tissue strand [19]. This behavior arises because deep in the tissue the intracellular and extracellular currents are in the opposite directions with almost the same magnitude, and therefore they tend to cancel each other.

The calculation is based on a piecewise linear model for the transmembrane potential. The upstroke is distributed over 2° in , which is 0.035 radians. In the middle of the heart wall, the radius is 35 mm, implying that the upstroke occurs over a distance of 1.22 mm. This is a slightly slower upstroke than is typically observed in cardiac tissue. The rapid rise of the action potential occurs in about 1 ms, and the conduction velocity is about 0.5 m/s parallel to the fibers, implying an upstroke distributed over about 0.5 mm. We performed a calculation with the upstroke over 1° and found the peak magnitude of the magnetic field is about 16 nT, a 13% increase compared to 2° case.

The cardiac tissue is assumed to have unequal anisotropy ratios [18, 24]. However, this condition is not essential for producing the magnetic field. We repeated the calculation for equal anisotropy ratios, = 0.032, = 0.2, = 0.032, and = 0.2 S/m [18, 24], and found that changed only slightly, to 12 nT.

Our calculation is based on a simplified model. The heart is spherical, the fiber geometry and wave front are idealized, and the transmembrane potential is uniform across the heart wall. These assumptions are not physiologically accurate but are necessary in order to calculate the biomagnetic field analytically. More realistic calculations would require numerical methods.

##### 4.1. How Action Currents Affect MRI

The groundbreaking work of Joy and his colleagues demonstrated that MRI can be used to measure current by detecting its resulting magnetic field [25–27]. However, much of this work has examined applied current, such as during defibrillation [28]. Detecting biocurrents is a greater challenge, because they are so small.

In living neural or cardiac tissues there exist various physiological processes such as blood flow, diffusion, and chemical exchange. All these physiological processes alter the magnetic field surrounding the atomic nuclei in ways that affect the magnetic resonance behavior of the nuclei [29]. Therefore, special MRI acquisition methods were designed or need to be designed accordingly to observe different physiological functions.

Functional magnetic resonance imaging (fMRI) makes use of the Blood Oxygenation Level Dependent (BOLD) signal changes [30]. Briefly, deoxyhemoglobin in blood is paramagnetic which will distort the external magnetic field of nearby hydrogen nuclei spins and further affect the MR image. When part of brain is active, the vascular system provides more blood and thereby reduces the fraction of blood that is converted from oxyhemoglobin into deoxyhemoglobin, and the magnetic distortion decreases. BOLD fMRI provides a method to measure oxygenation perfusion and thereby to indirectly estimate brain activity. However the spatial and temporal distribution of perfusion may differ from the distribution of electrical activity. Therefore, researchers would like to measure electrical activity directly using MRI.

Can the magnetic field produced by bioelectric currents be used as the gradient field for MRI? Theoretically, it is possible. The question is if the resulting phase shift is large enough to detect by current MRI systems. The phase shift of the MRI signal is given approximately by [5]where *γ* = 267 × 10^{6} radians/(s T) is the gyromagnetic ratio for protons and is the time over which the phase accumulates (the duration of the magnetic field for brief biomagnetic signals or a duration related to the MRI pulse sequence such as the echo time for long duration biomagnetic signals).

The phase shift determines if the electrical activity can be detected by the MRI system. In the brain, although the phase shift produced by a single neuron firing is not detectable, neural firings are not isolated events; in fact they are correlated with each other. There are two ways to increase the MR signal so that it is detectable. One is doing spatial and temporal integrals of neuronal firings over an MRI voxel according to the morphologies and the physiology of the neuronal events; the other one is altering the acquisition method such as using special pulse sequences to increase phase shift. Many researchers have been doing this both experimentally and theoretically [1–8, 31–34]. The MRI signals they measure are very small, and it is still unclear if they are detectable, whether the results are “fantasy, possibility, or reality” [3]. Recently, Sundaram et al. [5] performed neural current imaging in an intact turtle cerebellum using MRI. They were able to detect biomagnetic fields of strength of up to one nanotesla and phase shifts of about a third of a degree. These measurements required signal averaging with thousands of repetitions, imaged long duration metabotropic-induced neural signals lasting nearly a second, and had no interference caused by the vascular response (BOLD signal), tissue motion related to pulsatile blood flow, or other sources of physiological noise. Under these idealized conditions, phase shifts on the order of a third of a degree are detectable.

##### 4.2. Phase Shifts by Action Currents in the Heart

The brain is not the source of the largest biocurrents in the body; the heart is, so it may be easier to detect cardiac action currents using MRI compared to neural action currents. Deep in the heart, the magnetic field is small since the intracellular and the extracellular currents have almost the same value but are in opposite directions. However near the surface of the heart the extracellular current leaves the heart with the “return” loops through an adjacent volume conductor which produces a significant magnetic field [11].

From this analysis, we have achieved our original goal of computing the maximum value of the magnetic field produced by the action currents near the surface of the heart, 14 nT, which is significantly larger than any magnetic field produced by neuronal activation. For a firing time of 1 ms, this magnetic field can produce a maximum phase shift of 0.21° by (47). This phase shift is similar to that measured by Sundaram et al. [5], implying that the cardiac signal may be just measurable using MRI. However, these signals may be smaller than artifacts caused by heart motion and blood flow effects in a beating heart. Therefore, in vivo measurements of cardiac biocurrents in humans remain a challenge. Perhaps more promising would be to search for such signals in a Langendorff-perfused heart with an electromechanical uncoupler to eliminate contractions. In such an in vitro system—somewhat analogous to the in vitro turtle cerebellum model used by Sundaram et al. [5] to study signals from the brain—biocurrents in the heart may be measureable. Our conclusion is that the magnetic field in the heart may be detectable under idealized conditions but that using MRI to measure action currents in the humans will be difficult.

##### 4.3. The Potential of Detecting Cardiac Arrhythmias Using MRI

“Arrhythmia” is a term describing conditions in which the heart is beating too fast (tachycardia), too slow (bradycardia), or irregularly. Abnormal electrical activity exists in the heart in each of these conditions.

The usual way to detect a cardiac arrhythmia is using an electrocardiogram (ECG). Magnetic fields are also sometimes measured (the magnetocardiogram [23], or MCG). In MCG measurements, the biomagnetic field is measured outside the body and a difficult inverse problem must be solved to obtain the current sources that produce the magnetic field [35]. Simplifying assumptions are required to make the inverse problem solvable and stable, but these assumptions make the resulting current source distribution less precise and detailed. MRI would provide a way to measure the magnetic field inside the heart directly, avoiding the inverse problem of deducing the field inside the heart from measurements far outside the heart. Much research has been done that attempts to measure the magnetic field of a nerve using MRI. However, the heart, which produces a much larger magnetic field than nerves, may be a better place to start when attempting to use MRI for measurement of biocurrents. Of course, many challenges exist, not the least of which is the motion of the beating heart and the differences between the paramagnetic susceptibility of the blood within the heart, the cardiac tissue, and the tissues outside the heart [36, 37]. However, if successful, obtaining biocurrents from an MRI of the heart would provide a closer look at the function of the heart’s electrical conduction system and should help to diagnose a large number of cardiac arrhythmias.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This research was supported by the National Institutes of Health, Grant R01EB008421.