- Research
- Open Access

# A compact theory of magnetic nerve stimulation: predicting how to aim

- Charles F Babbs
^{1}Email author

**13**:53

https://doi.org/10.1186/1475-925X-13-53

© Babbs; licensee BioMed Central Ltd. 2014

**Received:**24 December 2013**Accepted:**17 April 2014**Published:**30 April 2014

## Abstract

### Background

A compact theory that predicts quantitatively when and where magnetic neurostimulation will occur is needed as a guide to therapy, ideally providing a single equation that defines the target volume of tissue excited by single or dual coils.

### Methods

A first-principles analysis of magnetic stimulation incorporating a simplified description of electromagnetic fields and a simplified cable theory of the axon yields a mathematical synthesis predicting how to aim.

### Results

Nerve stimulation produced by a single circular coil having one or more closely packed turns occurs in donut shaped volume of tissue beneath the coil. Axons spanning several millimeters are the sites of magnetic stimulation. The sites of maximal transmembrane depolarization in nerve fibers correspond to points where the axons enter or exit this volume of magnetically induced voltage and current. The axonal membrane at one end is depolarized locally during the rising phase of current in the coil. The axonal membrane at the opposite end is depolarized locally during the falling phase of current in the coil. Penetration depths of several centimeters from the skin surface or approximately one to two coil radii are practical. With two coils placed in a figure-of-eight configuration the separate clockwise and counterclockwise currents generate magnetic fields that add, producing maximal stimulation of a spindle shaped volume, centered at a depth of one-third to one-half coil radius from the body surface.

### Conclusions

This condensed synthesis of electromagnetic theory and cable theories of axon physiology provides a partial solution to the targeting problem in peripheral and in transcranial magnetic stimulation.

## Keywords

- Transcranial Magnetic Stimulation
- Magnetic Flux
- Magnetic Stimulation
- Myelinated Axon
- Magnetic Pulse

## Background

Magnetic neurostimulation is a remarkable phenomenon. It is electrodeless, nearly painless, and dry—easily penetrating skin and bone, which usually provide high resistance barriers to current injected through ordinary contact electrodes [1–4]. Stimulating magnetic fields are produced by high intensity electric current, flowing in thick wire coils a few centimeters in diameter that are placed on or near the skin surface. Direct contact with the body surface is not needed. An air gap is perfectly acceptable. Coils placed near the skull can stimulate neurons in the cerebral cortex, a process known as transcranial magnetic stimulation or TMS [5–9]. When coils are placed over the motor strips of cerebral cortex, just above the ears, TMS can be used to test the function of motor tracts extending from head to toe or from head to fingertips [8, 10, 11]. Such testing can be useful in monitoring the functional integrity of the spinal cord during neurosurgery, helping to avoid the dreaded complication of postoperative paralysis [3, 12]. When coils are placed further forward over pre-frontal cortex TMS can be used treat psychic depression [8].

The current-carrying coils used to produce magnetic stimulation have one major drawback: they are hard to aim, producing diffuse stimulation. To reach its full potential magnetic neural stimulation needs to be better targeted. Although mapping of the magnetic fields in three dimensions is well understood [13, 14]; less work has been done to characterize the induced electric fields and eddy currents [2, 4, 15–17] and even less to understand quantitatively the physiological mechanisms by which nerve stimulation happens in this unusual setting [18, 19]. A successful working theory should allow one to predict which neurons will be stimulated at which orientations and at which distances from the axes of the coils [19]. This latter problem requires particular attention to the effects on subthreshold membrane potentials of neurons [20]. Practitioners need to be able to better visualize which nerve fibers in which volumes of space are likely to be excited by a defined magnetic pulse from single or double coils with a defined geometry. What is needed as a guide to therapy is a compact model that predicts quantitatively when and where neurostimulation will occur—ultimately, perhaps, providing a single equation that defines a target volume of tissue.

Magnetic fields are easily mapped and computed [21], and they can be measured routinely in a dry laboratory using small search coils [22], although even the search coil method has been criticized [23]. However, the induced electric fields in volume conductors are not as easily computed or measured. Induced electrical potential gradients depend on the size and shape of any particular current loop or path and the presence of boundaries and barriers to current flow [16]. The induced electric field can be predicted from the negative of the time-derivative of the vector potential field using Finite Element Method packages such as Comsol®, but the computations are not easy for most clinical practitioners. Induced voltage gradients are also challenging to measure in the laboratory owing to electromagnetic interference from the high coil currents needed to produce the magnetic fields [23, 24].

In this context a mathematical modeling approach can be insightful. Many thinkers have been attracted to the complex and interesting problems of predicting the magnetic fields created by TMS coils of various shapes as well as predicting the induced electric fields in various anatomical models. Several quite elaborate models have been described [2, 15, 25] to predict magnetically induced eddy currents. These are based upon complex analytical treatments or upon finite element and finite difference models implemented by computer programs, some of which taking over an hour of computer time to execute [4] and most of which taking several tens of hours of human thinking time to fully comprehend. Once understood, such rigorous and detailed solutions for specific cases are not especially transportable or portable from one patient to the next or one biological application to the next.

For many relatively simple and clinically relevant geometries a compact analytical model could allow a more intuitive understanding of the influence of dominant parameters, based on the algebraic form of the final equation. Once derived, such an equation might well encapsulate much knowledge in compact notation and still be applicable to a large number of specific cases. This approach to understanding magnetic stimulation might be especially useful to clinicians not supported by a department of biomedical engineering that is capable of numerical modeling of many individual patients.

Thus the field of magnetic neural stimulation could benefit from a relatively straightforward, working theory that provides a way to visualize the three dimensional distribution of neurons that will be stimulated by a specified external coil—in other words, a pattern that tells the user how to aim the coil to stimulate a particular set of targets, including peripheral nerves, brain cortical regions, or even deeper brain structures, including if possible the directions of the nerve fibers in space that are likely to be activated. Such a simple, working scheme for visualization is not out of reach. The general features of magnetically induced voltage gradients, sufficient for modeling the physiology of nerve stimulation, can be calculated analytically because of the circular symmetry of the problem in a relatively uniform volume conductor like brain tissue. Simple boundary conditions, such as a flat or gently curved insulating surface and unlimited depth and width can be utilized [20]. Such conditions apply to brain covered by layers of skin, skull, and cerebrospinal fluid, or peripheral nerves covered by flat skin and subcutaneous tissue. Here the mathematics of electromagnetism and the mathematics of nerve impulse generation can be combined to create an analytical model of magnetic stimulation of brain and nerves that is sufficiently accurate for medical and biological applications, yet is mathematically compact, three dimensional, and easy to visualize.

This paper presents a fresh, clean-sheet analysis along these lines to clarify the underlying physics and to answer the following questions. How do magnetic fields stimulate nerve cells? What parts of the neuron are depolarized sufficiently by magnetic stimulation to generate an action potential? Where are the zones of stimulation located in the field beneath the exciting coils? How deep can stimulation be achieved for a given coil geometry? How does the spatial orientation of the nerve fiber tracts influence the ability of nerve fibers to be stimulated? How might the use of figure of eight coil configurations better target particular brain regions for focal stimulation?

The approach followed here is to begin with first principles (1) to describe the physics by which electric currents and voltages are produced in tissue by magnetic fields and then (2) to explore the mechanisms by which action potentials are subsequently initiated in neurons. The goal is to predict the approximate three dimensional patters of magnetically induced voltage and current and also to predict in terms of easily measured variables the target volume of space in which neurons are brought from their resting membrane potential to the threshold membrane potential for initiation of action potentials.

## Theory and results

### Part one: quantitative estimates of magnetically induced electric fields

#### Electromagnetism and electromagnetic induction

**B**, is the magnetic flux density in units of Tesla. For a circular coil, the magnetic or B-field induced by the coil has radial symmetry. In general, off-axis magnetic field vectors,

**B**, have both axial and radial components, B

_{z}and B

_{r}. The axial component, B

_{z}, is most important for magnetic neural stimulation. Along the central z-axis at x = 0, y = 0 the outward radial components, B

_{r}, created by current, I, flowing in paired short segments along the circumference on opposite sides of the coil cancel exactly, and so B

_{r}= 0. However, the z-axis components, B

_{z}, created by current, I, flowing in various short segments of the circumference of the coil add or reinforce one another to induce a well characterized [14] z-axis magnetic field.

_{0}is the magnetic permeability of free space in units of Henries/meter. This and other variables and mathematical symbols are defined for easy reference in Table 1, “Nomenclature”. If the coil has N closely packed turns, the result is multiplied by N. It is helpful to express tissue penetration distance, z, relative to the radius of the coil. This normalized distance is $\widehat{\mathrm{z}}=\mathrm{z}/\mathrm{R}$. Then the z-axis magnetic field is given by

**Nomenclature**

Variable | Units | Description |
---|---|---|

a, b | Regression constants describing radii at which axial field strength falls to zero | |

B | Tesla | Magnetic field strength |

B | Tesla | Maximal z-axis component of magnetic field |

C | Farads/cm | Specific membrane capacitance of nerve cells |

$\overrightarrow{\mathit{E}}$ | V/m | Vector of magnetically induced electric field |

E | V/m | Signed scalar magnitude of induced electric field around a circular path in homogeneous models |

| V/m | Scalar component of the voltage gradient along the path of an axon |

| Volts | Change in transmembrane potential of an axon |

I | Amps | Coil current |

I | Amps | Maximum coil current in time |

L | m | Length of an axon segment |

μ | Henries/m | Magnetic permeability of free space |

n | Number of coil turns or fold-increase in coil current | |

Φ | Tesla⋅m | Magnetic flux |

R | m | Coil radius |

r | m | Radial distance from z-axis in space |

$\widehat{\mathit{r}}$ | Normalized radial distance, r/R | |

R | Ohms | Resistance of axoplasm |

ρ | Ohm-m | Resistivity of intracellular fluid |

S | m | Surface area for magnetic flux |

s | m | Span or width of cell |

t | sec | Time |

θ | Radians | Angle between induced electric field and an axon segment |

| | Volts | Absolute value of voltage appearing along the length, L, of axon |

x, y, z | m | Spatial coordinates |

$\widehat{\mathit{x}},\phantom{\rule{0.5em}{0ex}}\widehat{\mathit{y}},\phantom{\rule{0.5em}{0ex}}\widehat{\mathit{z}}$ | Spatial coordinates normalized by coil radius, R | |

x | m | Radial distance from z-axis at which axial component of magnetic field becomes zero |

_{z}(z, r), diminishes with radial distance, r, away from the z-axis in a way characterized by much more complex equations [13], but easily appreciated graphically. In general, except when z < < R, the magnitude of B

_{z}(z, r) as a function of r is roughly Gaussian or dome shaped, with maximal values near r = 0, diminishing to near zero values at larger radii (Figure 2) [21, 22]. Very near the plane of the coil itself (z < < R) the field is more flat [21, 22] and less dome shaped. This near field region is typically outside the body during magnetic neural stimulation and its extra complexities can be ignored for present purposes.

At any particular depth, z, the magnetic field has axial and radial components, B_{z}(z, r) and B_{r}(z, r), sketched approximately in Figure 2 to show relative shapes and strengths. Of note are the radii r_{0}, at which the axial field strength falls to zero. The values of r_{0} increase gradually as a function of depth, z. As distance, z, from the plane of the coil increases, the z-directed components of the magnetic field spread out and diminish in intensity.

#### Induced electric fields in homogeneous underlying tissues

**B**, normal or perpendicular to the plane of dS, that is

**E**, around the loop is

**E**around the loop being defined by the right hand rule. A more modern and more general statement of Faraday’s law is the integral form. If we define the electric field, $\overrightarrow{\mathrm{E}}$ (volts/meter), along a wire or conducting loop for any shape of closed path consisting of multiple steps, $\mathrm{d}\overrightarrow{\mathrm{L}}$, in a time-varying magnetic field, then the result of Faraday’s law is

where the integral is taken over the area bordered by the closed path [14]. In Figure 1 when the current, I, in the coil changes rapidly, the induced magnetic field components B_{z} are normal to the x-y plane and correspond to B_{n} in Equation (2b).

Measuring EMF induced by a rapidly changing magnetic field using a search coil in a dry laboratory is a straightforward task. One simply fashions a small loop of insulated wire about 1 cm in diameter and measures the voltage appearing between the ends of the wire. By tilting the plane of the coil one can find the direction of the field. In a volume conductor, however, the induced currents are constrained by insulated boundaries, and some further insight is required.

#### Radial components of the B-field can be ignored

_{r}of the magnetic field are greatly reduced by destructive interference. As a practical matter, if the plane of the coil is parallel or tangent to the body surface, the radial components B

_{r}(z, r) can be ignored in neural stimulation [26].

#### Estimating z-axis components of the B field

_{z}(z, r), because of the radial symmetry of the problem, the differential surfaces dS in Equation (2b) are nested rings of diameters, 2r, circumference 2πr, and width, dr. The magnetic flux at any depth, z, is therefore,

_{z}(z, r), describing axial field B

_{z}as a function of radius r is complex and contains elliptic integrals that cannot be simply computed [13, 14]. However except for cases of z < < R, the shape of B

_{z}(r) in three dimensions is generally a dome or bell shaped surface. In this region of biological interest, the strength of B

_{z}(r) falls from a maximal value at r = 0 to zero at r = r

_{0}[22]. The value of r

_{0}gradually increases with depth as the strength of B

_{z}(r) weakens (Figure 4). Here the extent of the dome shaped region of positive z-directed magnetic field is computed from the exact formula [27] and plotted as a function of axial depth beneath a simple circular coil. The dashed line is the linear regression ${\widehat{\mathrm{r}}}_{0}\left(\mathrm{z}\right)=1.0\phantom{\rule{0.12em}{0ex}}\left(\widehat{\mathrm{z}}+0.87\right),$ where the distances, r and z, embellished with the “hat” symbols, are normalized by the coil radius, R, so that $\widehat{\mathrm{r}}=\mathrm{r}/\mathrm{R}$ and $\widehat{\mathrm{z}}=\mathrm{z}/\mathrm{R}.$

Approximating the dome shaped function B_{z}(z, r) as the parabola ${\mathrm{B}}_{\mathrm{z}}\left(\mathrm{z},0\right)\left(1-\frac{{\mathrm{r}}^{2}}{{\mathrm{r}}_{0}^{2}\left(\mathrm{z}\right)}\right)$, in keeping with more detailed and exact computations [13, 22], one can take advantage of the circular symmetry of the problem to find the total magnetic flux through circles of any radius, 0 ≤ r ≤ r_{0}, centered along the z-axis and perpendicular to it. This approximation greatly simplifies the calculation of magnetic and induced electric fields and is sufficiently accurate for the purpose of understanding the biological response to magnetic stimulation. If desired, the parabolic representation of B_{z}(z, r) can be expanded to a higher order polynomial in r for any desired degree of accuracy. The following approach to integration will still be valid.

_{0}is therefore

where B_{z}(z,0) is given by Equation (1a). For radii r > r_{0} one can estimate the magnetic flux through a loop of electrical conductor with radius r as *Φ*(*z*, *r*) = *Φ*(*z*, *r*
_{0}), ignoring low amplitude reversed flux at outer radii, which are outside the region of interest in magnetic stimulation. (Alternatively, a higher order polynomial can be used.) In this way the magnetic flux through a circular domain of the volume conductor with radius r can be estimated from the peak on-axis value in Equation (1), multiplied by the area of the circle and the attenuation factor in parentheses in Equation (5b).

_{0}, assume that dΦ(z, r) ≈ 0 for r > r

_{0}, so that Φ(z, r) = Φ(z, r

_{0}). Then using Equation (5) it is straightforward to show that

For example, one can compute the magnitude of induced EMF caused by a typical magnetic pulse, for which in practical units the maximal current output I_{max} ≈ 10,000 amperes, Δt ≈ 100 microseconds, giving dI/dt ≈ 10^{8} A/sec [4], with the magnetic permeability constant, μ_{0} = 4π × 10^{−7} V · sec /(A · m). Current in the coil is typically produced by repeated capacitor discharges, which are damped by the inductance of the coil and opposed also by the small coil resistance. In such a discharge circuit the current rises quickly to a maximal level and then decays quickly [28] the rise time Δt_{1} is on the order of 100 microseconds and the fall time Δt_{2} is on the order of 200 microseconds [28], leading to reasonable estimates of (dI/dt).

Figure 5 shows the predicted strength of the magnetic field in Volts/meter given by Equations (7c) and (7d). The induced voltage gradients and eddy currents flow parallel to the surface in circles concentric with the coil, and opposite in direction to the direction of current increase in the coil, as expected from Faraday’s Law. In this simple system there are no insulating boundaries that obstruct current flow and no accumulation of surface charge at insulating boundaries [4]. The resistivity of the tissue of the volume conductor in this model is also homogeneous. The results are quite similar to those reported by Roth and Basser [18] and by Tofts [16] using much more elaborate methods.

The direction of E depends on the direction of current flow in the coil and also on whether the current is increasing or decreasing. If the B-field is rising to its maximal value the induced EMF will be negative in the sense of the right hand rule and Equation (6). When the B-field is falling from its maximal value the direction of the induced EMF will be reversed. The results in Figure 5 are for a single coil with a single turn. For multi-turn coils the current derivative term (dI/dt) is multiplied by the number of turns. As explained subsequently, the same approach can be applied to double figure of eight coils by doubling the induced field in the region of overlap.

### Part two: interaction of neurons with induced electric fields

#### Neuron anatomy and physiology

Anatomically typical nerve cells, or neurons, are most unusually shaped, having a cell body at one end that is perhaps 25 micrometers in diameter and a very long, thin arm or projection of membrane covered cytoplasm known as the axon that is roughly one micron in diameter and can span distances of one thousand to one million micrometers in some cases [29]. Excitation of a nerve cell means the initiation of a self-propagating wave of depolarization known as an action potential, usually starting at one end near the cell body and continuing at speeds of several meters per second down the length of the axon. This special anatomy and physiology allows neurons to send signals at high speed over macroscopic distances from one part of the body to another [30].

Normal neurons are polarized in the resting state by a charge difference across the outer membrane of the cell of (inside negative, outside positive) of approximately −85 mV [30]. In order for an action potential to be generated the transmembrane potential at a particular site on the surface of the cell must be brought from the resting level of −85 mV to a threshold level near −55 mV, at which voltage sensitive sodium channels in the membrane open to allow depolarizing current to flow in the form of charged Na^{+} ions. This current eliminates the transmembrane charge difference locally and triggers depolarization of adjacent membrane, sending a propagated signal along the length of the axon [31].

#### Cell bodies are not stimulated directly

_{ICF}C/2, the product of the intracellular fluid resistance and net capacitance, or

_{m}, per unit area of cell membrane. For ρ = 200 ohm-cm [33] and C

_{m}= 1 microfarad per square centimeter [30] and s = 25 × 10

^{−4}cm [29] the time constant, τ, is here about 0.25 microseconds. Magnetically induced pulses last about 100 microseconds [28] so there is plenty of time to fully charge the cell body after a magnetically induced E-field appears across it.

The full magnitude of the potential appearing across the cell, however, is not significant physiologically. In the case of Figure 5, taking the maximal gradient as 25 V/m, the total voltage developed across the cell is only 25 μV/μm × 25 μm = 625 μV or 0.6 mV, a far cry from the 30 or so mV needed to depolarize one side of the cell from the resting membrane potential to the threshold potential. This same reasoning also applies to larger cell bodies, like pyramidal cells in the brain. In the same vein it is obvious that electric fields that appear at right angles to axons in nerve fiber tracts cannot induce a threshold membrane potential change, since the diameters of axons are on the order of only 1 micrometer. This conclusion agrees with the prior work of Nowak and coworkers [34, 35].

#### Axons spanning several millimeters are the sites of magnetic stimulation

_{a}, of the axoplasm in the segment, as shown in Figure 7. During sub-threshold stimulation (before voltage sensitive ion channels open) this electrical arrangement is very similar to the cell model of Figure 6, with a very short, sub-microsecond time constant for charging. This means that all segments of small length, ΔL, will come up to the full gradient voltage E

_{L}ΔL very quickly. When the membrane capacitances along the cable model became fully charged, then the steady-state potential difference between the two widely separated points along the cable model, P

_{1}and P

_{2}, would be

where *E*
_{
L
} is the scalar component of the electric field along the path of the axon, and points P_{1} and P_{2} are points where the axon enters and exits the induced voltage field. In this scenario ± 0.5 ΔV_{L} will appear very quickly across the axonal membrane segments at P_{1} and P_{2}.

_{n}, between nodes of Ranvier in the brain is about 300 micrometers [36], and the axon diameter is about 0.5 micrometers, the gap, s, is about 0.3 micrometers. In turn, the time constant for charging segments of a myelinated axon is, in keeping with Equation (10),

or about 2.4 microseconds. This means that after about 5 microseconds the membrane sections are almost fully charged and the whole voltage gradient appears between the ends of a myelinated axon segment. Since the rise or fall times, Δt, for magnetic pulses are approximately 100 microseconds long, there is plenty of time to charge the membranes of myelinated axons, as well as non-myelinated ones, to the full gradient potential, ΔV_{L}, created by magnetic induction.

_{L}, consider a nerve axon having length, L, that passes through the torus of magnetically induced current and voltage in three dimensions. Let the midpoint of length, L, be located at coordinates z, r (Figure 9), where the induced voltage gradient is $\overrightarrow{\mathrm{E}}\left(\mathrm{z},\mathrm{r}\right)$. The axon makes angle θ in three dimensions with the induced electric field vector $\overrightarrow{\mathrm{E}},$ which itself is perpendicular to r.

_{1}and P

_{2}will be either hyperpolarized or depolarized by amount.

_{0}defined as in Figure 2, gives a simple analytical expression to describe the effect of magnetic stimulation on the transmembrane potentials of affected neurons in terms of the rate of change in coil current, dI/dt, the distance from the plane of the coil in units of coil radius ($\widehat{\mathrm{z}}=\mathrm{z}/\mathrm{R}$), the normalized radial distance $\widehat{\mathrm{r}}=\mathrm{r}/\mathrm{R}$ from the z-axis through the center of the coil, the length, L, of the axon segment in the induced electric field, and the angle, θ, of the axon segment with respect to the field.

^{8}A/sec. Distances are normalized to the radius, R, of the coil. The projection distance, L cos θ, is 1 cm. The values of the contours can therefore be interpreted as millivolts of membrane potential change per centimeter of projection distance along the axis of the induced field. The plotted values are also proportional to the number of turns of the coil and to the peak coil current. Axons travelling through the indicated regions of space more or less parallel with the induced electric field (into and out of the page) will be stimulated to the extent shown in Figure 11. In three dimensions the stimulation zones for a simple circular coil are toroidal or donut shaped.

#### Figure of eight coils

_{z}, are opposite in sign and may interact as sketched in Figure 13. The forgoing analysis of individual circular coils can be applied with modification to the figure of eight coils, here in a semi-quantitative way. B

_{z}-fields add algebraically in the vicinity of the z-axis, but remain largely undisturbed elsewhere by the presence of the other coil. One can imagine the domes of instantaneous magnetic field strength as two bowls placed side-by-side, one inverted, with the complete rims of both bowls in the same x-y plane (Figure 13, bottom).

Helpful reinforcing effects are possible in simple cases without severely restricting lateral boundaries. In a semi-infinite space a double torus of induced voltage and current is created. As the absolute intensities of the paired B_{z}-fields rise and fall in time there are counter rotating induced currents that can combine and reinforce the voltage gradient near the common z-axis. As coil current rises and falls the heights of the bowls rise and fall. If current reverses direction, the bowls flip vertically. Parallel components of induced current near the z-axis add to reinforce each other. At a depth of one third to one half coil radius (see below) a sweet spot exists where induced voltage gradients are nearly double what they would have been at the same z-level for a single coil [37]. In the far field the overlapping B-fields near the z-axis cancel, but only in the zone of overlap. Elsewhere around the circumferences of the individual loops of the figure of eight, circular symmetry prevails and eddy currents are still induced, as they would be for a single coil.

_{max}, is adjusted so that the single coil stimulus never exceeds threshold, but twice that value does. The focal region will be roughly spindle shaped and correspond to the overlapping region for two toroidal volumes of high induced EMF from either coil. This scheme provides a practical solution to the problem of focusing magnetic neurostimulation and agrees with the prior results of Ravazzini [26]. This concept can be extended to 4 coil fields but vector addition and cancellation limit the boost in net induced EMF from two to four coils to a factor of no more than the . The reinforcing double toroidal patterns of induced current flow predicted by the present compact analysis are quite similar to those found by Wagner et al. [25] using much more elaborate computational methods. The maps of induced current predicted here are also similar to those described by Pascual-Leone [5]. Figure of eight coils produce neurostimulation that is roughly twice as strong as that produced by single coils and confined to a much smaller, targeted volume of tissue.

## Discussion

Axons running in many directions, but not directly parallel to the z-axis of an electromagnetic coil carrying rapidly changing current can be depolarized or hyperpolarized by induced electric fields at depths up to one coil diameter from the plane of the coil. Both myelinated and non-myelinated axons can be stimulated. The sites of greatest membrane depolarization are located toward the ends of the axon segments traversing the region of induced voltage and current. The target volume for magnetic stimulation takes the shape of a toroid centered along the z-axis with major radius similar to that of the current carrying coil and its minor radius varying in size, depending on the coil current. When figure of eight coils are used there is the opportunity to create spindle shaped target volumes of stimulation, centered on the z-axis at a depth of about one third to one half coil radius.

In the forgoing analytical treatment many simplifying assumptions have been introduced to make the mathematics more tractable and useful, realizing than in any particular biological experiment there is some lack of precision in specifying geometric distances and physiological conditions. Hence, only the most dominant terms in the mathematics need to be considered in a biologically satisfactory solution to the targeting problem. For example, the position of a coil in a clinical stimulation protocol can only be specified within about one millimeter. The thickness of a practical coil is on the order of several millimeters, perhaps better characterized as a bundle of one dimensional wire loops, blurring the idealized fields calculated for any single loop. The positions of underlying nerves and brain structures vary from subject to subject and typically cannot be known exactly in a particular treatment situation or experiment. Similarly, the intracellular and extracellular sodium and potassium ion concentrations that determine resting and threshold membrane potentials can vary among subjects, the intracellular values being especially hard to measure. Additionally there are motion artifacts, including subtle movement in time and space with breathing, fidgeting, as well as any muscle movement caused by nerve stimulation itself.

Accordingly, the positions of biologic structures with respect to magnetic and electric fields can only be specified within a millimeter or two at best. For such reasons over-precise prediction is a fool’s errand. What is more useful is the ability to predict within, say, 10 percent, the locations and directions of axons likely to be stimulated by a given practical apparatus, and especially the ability to visualize in three dimensions the locations and orientations of nerve axons likely to be stimulated by a given coil and current intensity. Equations (14a) and (14b) provide a compact, closed form solution for the effects of magnetic neurostimulation in terms of the dominant variables: the slew rate of coil current, the axial and radial distances of the target axon from the center of the coil, the length of the axon within the induced electric field, and the angle of the axon in space. This kind of predictive capability has not been available heretofore.

The theory presented here is validated by experimental observations. When using a single circular coil to stimulate the median nerve in the wrist, Maccabee et al. [37] found the most sensitive position for stimulation the nerve was located beneath the middle circumference of the coil in the position of maximal toroidal current. However when figure of eight coils were used, the most sensitive position for stimulation of the nerve was located at the midpoint between the two loops of the coil. Rudiak and Marg [38] used a clever approach to estimate the effective depth of magnetic brain stimulation in human subjects. For figure of eight coils with 10 cm diameter (5 cm radius) loops, the focal depth of stimulation was between 1.8 and 2.1 cm, or between 36% and 42% of coil radius, very close to that predicted by the present analysis (Figure 14).

## Conclusion

Magnetic neural stimulation is a remarkable and subtle phenomenon, able to penetrate highly resistive skin and bone with ease and stimulate underlying nerve fibers. This paper presents a condensed synthesis of electromagnetic theory and cable theories of axon physiology to better inform further development and clinical practice of magnetic neurostimulation, including transcranial magnetic stimulation. With knowledge of the critical variables and a little imagination it is possible to visualize in three dimensions the best way to arrange and orient surface coils to achieve stimulation within a defined volume of underlying tissue.

## Declarations

## Authors’ Affiliations

## References

- Alfonsetti M, Clementi V, Iotti S, Placidi G, Lodi R, Barbiroli B, Sotgiu A, Alecci M:
**Versatile coil design and positioning of transverse-field RF surface coils for clinical 1.5-T MRI applications.***MAGMA*2005,**18:**69–75. 10.1007/s10334-004-0090-4View ArticleGoogle Scholar - Esselle KP, Stuchly MA:
**Neural stimulation with magnetic fields: analysis of induced electric fields.***IEEE Trans Bio-Med Eng*1992,**39:**693–700. 10.1109/10.142644View ArticleGoogle Scholar - Kobayashi M, Pascual-Leone A:
**Transcranial magnetic stimulation in neurology.***Lancet Neurol*2003,**2:**145–156. 10.1016/S1474-4422(03)00321-1View ArticleGoogle Scholar - Roth BJ, Saypol JM, Hallett M, Cohen LG:
**A theoretical calculation of the electric field induced in the cortex during magnetic stimulation.***Electroencephalogr Clin Neurophysiol*1991,**81:**47–56. 10.1016/0168-5597(91)90103-5View ArticleGoogle Scholar - Pascual-Leone A, Bartres-Faz D, Keenan JP:
**Transcranial magnetic stimulation: studying the brain-behaviour relationship by induction of 'virtual lesions'.***Philos Trans R Soc Lond B Biol Sci*1999,**354:**1229–1238. 10.1098/rstb.1999.0476View ArticleGoogle Scholar - Pascual-Leone A, Nguyet D, Cohen LG, Brasil-Neto JP, Cammarota A, Hallett M:
**Modulation of muscle responses evoked by transcranial magnetic stimulation during the acquisition of new fine motor skills.***J Neurophysiol*1995,**74:**1037–1045.Google Scholar - Roth Y, Zangen A, Hallett M:
**A coil design for transcranial magnetic stimulation of deep brain regions.***J Clin Neurophysiol*2002,**19:**361–370. 10.1097/00004691-200208000-00008View ArticleGoogle Scholar - Terao Y, Ugawa Y:
**Basic mechanisms of TMS.***J Clin Neurophysiol*2002,**19:**322–343. 10.1097/00004691-200208000-00006View ArticleGoogle Scholar - Zangen A, Roth Y, Voller B, Hallett M:
**Transcranial magnetic stimulation of deep brain regions: evidence for efficacy of the H-coil.***Clin Neurophysiol*2005,**116:**775–779. 10.1016/j.clinph.2004.11.008View ArticleGoogle Scholar - Stuchly MA, Esselle KP:
**Factors affecting neural stimulation with magnetic fields.***Bioelectromagnetics*1992,**13**(Suppl 1):191–204.View ArticleGoogle Scholar - Salvador R, Silva S, Basser PJ, Miranda PC:
**Determining which mechanisms lead to activation in the motor cortex: a modeling study of transcranial magnetic stimulation using realistic stimulus waveforms and sulcal geometry.***Clin Neurophysiol*2011,**122:**748–758. 10.1016/j.clinph.2010.09.022View ArticleGoogle Scholar - Takahashi S, Vajkoczy P, Picht T:
**Navigated transcranial magnetic stimulation for mapping the motor cortex in patients with rolandic brain tumors.***Neurosurg Focus*2013,**34:**E3.View ArticleGoogle Scholar - Schill RA:
**General relation for the vector magnetic field of a circular current loop: a closer look.***IEEE Trans Magnetics*2003,**39:**961–966. 10.1109/TMAG.2003.808597View ArticleGoogle Scholar - Kip AF:
*Fundamentals of Electricity and Magnetism*. 2nd edition. New York: McGraw-Hill; 1969.Google Scholar - Salinas FS, Lancaster JL, Fox PT:
**3D modeling of the total electric field induced by transcranial magnetic stimulation using the boundary element method.***Phys Med Biol*2009,**54:**3631–3647. 10.1088/0031-9155/54/12/002View ArticleGoogle Scholar - Tofts PS:
**The distribution of induced currents in magnetic stimulation of the nervous system.***Phys Med Biol*1990,**35:**1119–1128. 10.1088/0031-9155/35/8/008View ArticleGoogle Scholar - Mouchawar A, Nyenhuis J, Bourland JD, Geddes LA, Schaefer DJ, Riehl MR:
**Magnetic stimulation of excitable tissue: calculation of induced eddy-currents with a three-dimensional finite-element model.***IEEE Trans Magnetics*1993,**29:**3355–3357. 10.1109/20.281174View ArticleGoogle Scholar - Roth BJ, Basser PJ:
**A model of the stimulation of a nerve fiber by electromagnetic induction.***IEEE Trans Biomed Eng*1990,**37:**588–597. 10.1109/10.55662View ArticleGoogle Scholar - Amassian VE, Eberle L, Maccabee PJ, Cracco RQ:
**Modelling magnetic coil excitation of human cerebral cortex with a peripheral nerve immersed in a brain-shaped volume conductor: the significance of fiber bending in excitation.***Electroencephalogr Clin Neurophysiol*1992,**85:**291–301. 10.1016/0168-5597(92)90105-KView ArticleGoogle Scholar - Rotem A, Moses E:
**Magnetic stimulation of one-dimensional neuronal cultures.***Biophys J*2008,**94:**5065–5078. 10.1529/biophysj.107.125708View ArticleGoogle Scholar - Alfonsetti M, Sotgiu A, Alecci M:
**A theoretical and experimental study on transverse field radio frequency surface coils.***Measurement*2010,**43:**1503–1515. 10.1016/j.measurement.2010.08.017View ArticleGoogle Scholar - Cohen LG, Roth BJ, Nilsson J, Dang N, Panizza M, Bandinelli S, Friauf W, Hallett M:
**Effects of coil design on delivery of focal magnetic stimulation. Technical considerations.***Electroencephalogr Clin Neurophysiol*1990,**75:**350–357. 10.1016/0013-4694(90)90113-XView ArticleGoogle Scholar - Tofts PS, Branston NM:
**The measurement of electric field, and the influence of surface charge, in magnetic stimulation.***Electroencephalogr Clin Neurophysiol*1991,**81:**238–239. 10.1016/0168-5597(91)90077-BView ArticleGoogle Scholar - Basham E, Zhi Y, Wentai L:
**Circuit and coil design for in-vitro magnetic neural stimulation systems.***IEEE Trans Biomed Circuits Syst*2009,**3:**321–331.View ArticleGoogle Scholar - Wagner TA, Zahn M, Grodzinsky AJ, Pascual-Leone A:
**Three-dimensional head model simulation of transcranial magnetic stimulation.***IEEE Trans Biomed Circuits Syst*2004,**51:**1586–1598. 10.1109/TCSI.2004.832793View ArticleGoogle Scholar - Ravazzani P, Ruohonen J, Grandori F, Tognola G:
**Magnetic stimulation of the nervous system: induced electric field in unbounded, semi-infinite, spherical, and cylindrical media.***Ann Biomed Eng*1996,**24:**606–616. 10.1007/BF02684229View ArticleGoogle Scholar - Calculator for Off-Axis Fields Due to a Current [http://www.netdenizen.com/emagnet/offaxis/iloopcalculator.htm]
- Davey K, Epstein CM:
**Magnetic stimulation coil and circuit design.***IEEE Trans Biomed Circuits Syst*2000,**47:**1493–1499.Google Scholar - Di Fiori MSH:
*Atlas of Human Histology*. 4th edition. Philadelphia: Lea & Febiger; 1981.Google Scholar - Boron WF, Boulpaep EL:
*Medical Physiology*. 2nd edition. Philadelphia: Elsevier; 2005.Google Scholar - Plonsey R, Barr RC:
**Electric field stimulation of excitable tissue.***IEEE Eng Med Biol Mag*1998,**17:**130–137. 10.1109/51.715497View ArticleGoogle Scholar - Brophy JJ:
*Basic Electronics for Scientists*. New York: McGraw-Hill; 1966.Google Scholar - Geddes LA, Baker LE:
*Prinicples of Applied Biomedical Instrumentation*. 2nd edition. 1975.Google Scholar - Nowak LG, Bullier J:
**Axons, but not cell bodies, are activated by electrical stimulation in cortical gray matter. II. Evidence from selective inactivation of cell bodies and axon initial segments.***Exp Brain Res*1998,**118:**489–500. 10.1007/s002210050305View ArticleGoogle Scholar - Nowak LG, Bullier J:
**Axons, but not cell bodies, are activated by electrical stimulation in cortical gray matter. I. Evidence from chronaxie measurements.***Exp Brain Res*1998,**118:**477–488. 10.1007/s002210050304View ArticleGoogle Scholar - Lasiene J, Matsui A, Sawa Y, Wong F, Horner PJ:
**Age-related myelin dynamics revealed by increased oligodendrogenesis and short internodes.***Aging cell*2009,**8:**201–213. 10.1111/j.1474-9726.2009.00462.xView ArticleGoogle Scholar - Maccabee PJ, Eberle L, Amassian VE, Cracco RQ, Rudell A, Jayachandra M:
**Spatial distribution of the electric field induced in volume by round and figure '8' magnetic coils: relevance to activation of sensory nerve fibers.***Electroencephalogr Clin Neurophysiol*1990,**76:**131–141. 10.1016/0013-4694(90)90211-2View ArticleGoogle Scholar - Rudiak D, Marg E:
**Finding the depth of magnetic brain stimulation: a re-evaluation.***Electroencephalogr Clin Neurophysiol*1994,**93:**358–371. 10.1016/0168-5597(94)90124-4View ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.