A Miniaturized Ultra-Focal Magnetic Stimulator and Its Preliminary Application to the Peripheral Nervous System

Transcranial magnetic stimulation (TMS) is a noninvasive brain stimulation technique used in the clinic to treat several neurological disorders and psychiatric diseases. One of TMS’s significant limitations is its low spatial resolution, which often results in a mismatch between the target area in the brain and the stimulation site on the scalp. To enhance its spatial resolution, we designed and built a complete stimulation system complete with a millimetric-diameter coil and microscopic traces (μCoil). The first tests conducted on healthy volunteers showed that the μCoil stimulation of the radial nerve in the wrist could indeed evoke somatosensory nerve action potentials (SNAPs). In this chapter, we study this nerve stimulation system with electromagnetic and neuron simulators on a neurofunctionalized model from the Virtual Population (ViP v.4) and a μCoil figure-8 geometry. In particular, we study how changes in the μCoil geometry, such as the number of layers, shape, and length of an iron or air core, may help to promote the generation of somatosensory nerve action potentials.

major restriction of TMS in its present form is its limited spatial resolution, which makes the task of focusing the stimulation exclusively to the targeted cortical region very challenging. The spatial resolution of TMS depends on the stimulating coil's geometry and dimensions. A deep investigation of the focality of the electric field induced by 50 TMS coils has been conducted by Deng et al. [20]. In this work, the authors showed that figure-of-eight coils are more focal than circular coils. For example, the Magstim 70 mm figure-8 coil (P/N 9925, 3190) has a 15 cm 2 focality [20,21], while the Magstim 90 mm circular coil (P/N 3192) has a 70 cm 2 focality [20,21]. Also, smaller coils are, in general, more focal than larger coils, with focality that can be as small as 5 cm 2 [20]. For this reason, several attempts have been made to increase the spatial resolution of TMS by reducing the dimension of the stimulating coil [22][23][24]. To date, however, all the coils with reduced dimensions were designed for either invasive [24] or noninvasive animal applications [23]. In fact, as the dimension of the coil shrinks (diameter less than 2 cm), larger electric currents (few hundreds of kiloamperes) are required to produce enough induced electric field in the human cortex (~60-100 V/m). Given the conventional engineering approach which uses a few turns (<30) of wire winding inside an insulation housing, passing such large currents would produce excessive amounts of heat and large magnetic forces that may raise safety concerns or exceed specifications of conventional wires. To overcome such technological barriers, we have introduced flex circuit technology to design a new generation of miniaturized coils, which allows for a high number of turns (100 or greater) in a multilayer structure characterized by a reduced diameter (15-20 mm). Such an approach has been at first used to develop microscopic coils for invasive magnetic stimulation (μMS) capable of activating neuronal circuitry both in vitro [25] and in vivo [26]. Based on these results, a flex circuit technology was adapted to develop the first generation of miniaturized coils for noninvasive magnetic stimulation (μCoil). The concept behind the structure of the μCoil consists of a certain number of parallel copper traces, named layers, deposited on a long flexible substrate wound in N turns around a core. The prototype of μCoil was built using a four-trace Kapton sheet, wound in 123 turns around a copper pin with a 1 mm diameter. Since each layer was 1 mm high, with a distance of 0.5 mm between two adjacent layers, we fabricated a single 8-mm-tall circular solenoid (μCoil), with an outer diameter of 15 mm. A figure-of-eight stimulator was obtained by pairing two single μCoils, thus obtaining a structure with a maximum dimension of 30 mm. To date, these coils are the only figure-of-eight coils capable of performing stimulation of the peripheral nervous system on healthy volunteers. In this chapter, we studied numerically the behavior of six different μCoil geometries to show how each coil geometry affects the peripheral nerve stimulation threshold.

μCoil Modeling
The electromagnetic simulation software Sim4Life (v.5, Zurich MedTech, Zurich) was adopted in all of the simulations. A planar loop coil geometry represented the actual copper trace configuration. The space between two consecutive loops of 1.5 mm was selected to model the trace thickness and spacing. Six different geometries were considered, as shown in Fig. 1. The first geometry studied was a two-layered single circular coil, with each layer wound in 123 turns around a 1 mm air core. Second, two other layers, with the same characteristics as above, were added on top to make a four-layered single circular coil. The third geometry herein studied was the figure-of-eight configuration, which consisted of two four-layered single circular coils paired together and with currents circulating in the opposite phase. Three additional figure-of-eight coil models were built to study how the presence and length of an ideal iron core (i.e., linearity assumption: B ¼ μ r H with μ r ¼ 60) would modify the magnetic and electric field generated by the μCoil. Each solenoid was wound in 100 turns and had an inner and outer diameter of 10 mm and 21 mm, respectively. Three lengths of the iron core were considered: 9 mm, 23 mm, and 41 mm. Of all the geometries, the circular four-layered μCoil and the figure-ofeight μCoil were manufactured and tested. The iron-core μCoils were studied to Fig. 1 Four figure-of-eight μCoil models. μCoil models and placement over the right arm of the neurofunctionalized model Yoon-Sun and the point where the two figure-of-eight μCoils touch was precisely above the nerve trajectory in the wrist area enhance the magnetic field strength of the newly developed μTMS coil for future cortical applications. The electromagnetic simulations were conducted using the Magneto Quasi-Static solver of the software Sim4life (v.5.9, Zurich MedTech, Zurich). This solver implements the quasi-static approximation that allowed decoupling the electric field from the magnetic field. Considering the source current density J, magnetic vector potential A was calculated from Ampere's law and used to compute the B and E field. EM simulations were performed by feeding a sinusoidal current with a 1.9 kHz frequency in each μCoil and a phase of 0 in the first μCoil and of 180 in the second μCoil when considering the figure-of-eight coil. The two coils were placed in a figure-of-eight coil over the arm of the anatomical, neurofunctionalized body model Yoon-Sun (ViP., v4 [27]), in correspondence of the superficial branch of the radial nerve, as shown in Fig. 1. The Yoon-Sun model is characterized by a detailed representation of the peripheral nervous system and the various nerve fibers. This simulation setup was validated by experiments conducted on healthy volunteers [28]. The μCoil driving sinusoidal currents were set to 1 A for comparison among different configurations.

Modeling Peripheral Nerve Stimulation: Titration Analysis
To study the interaction between the induced electric field and the peripheral nervous system, the results of the EM simulations were fed into the dynamic Neuron solver version (v. 5.9) embedded in Sim4Life in a Multiphysics approach. The excitable behavior of the fibers was simulated using the McIntyre-Richardson-Grill model, which is based on a double-cable representation of the axon that allows separating electrical representations for the myelin and underlying internodal axolemma [29,30]. Nerves can be excited by an external potential, V ext, defined as in Eq. 1, where C is the linear path connecting each data point to an arbitrary zero voltage reference [31].
In this work, a titration analysis was conducted for all the simulated μCoil geometries. Titration is the iterative process that allows determining the stimulation threshold of a nerve by stimulating it with pulses of increasing intensity. Figure 2 shows the vector view of the B field generated by each μCoil on the central cross-sectional plane. Inside the coil, the presence of more layers directs the B field along the axis of the coil. The magnetic field lines of forces diverge outside the coil and form closed loops. Doubling the number of layers increased the maximum intensity of B induced on the skin from 30 to 40 mT, thus by a factor of 33% without changing its direction. When forming the eight-shaped coil by pairing two circular μCoils fed in counter phase, the B field lines formed a third loop at the center of the figure-of-eight coil, due to the coupling between the two coils. The use of the iron core modified the direction of the B field that was along the coil axis for the whole length of the core and parallel to the surface of the coil in the space between the two cores. At the interface between iron and air, vector B was deviated to respect the continuity of the normal component to the surface. Moreover, at the level of the nerve (i.e., 2.5 mm from the surface of the skin) B field intensity for the coil with the 41 mm, long core is higher than that generated with the 23 mm long core and with the 9 mm long core, by 15% and 50%, respectively.

Electric Field Induced by the μCoils
The single circular coil that was composed of four layers induced inside the tissues an electric field stronger than the field induced by the two-layered coil: the maximum E field induced at the level of the nerve with 1 A was 0.27 V/m for the two-layered circular coil and 0.40 V/m for the four-layered one (see Fig. 3, Panel a and b). The figure-of-eight coil configuration increased the intensity of the induced electric field on the nerve up to 0.6 V/m (Fig. 3, panel C). Moreover, the direction of the B field (Fig. 2) induced an E field-oriented mostly along the longitudinal axis of the arm and (therefore) along the radial nerve. The presence of an iron core further oriented E along the nerve, as shown by the vector arrows in Fig. 3 d-f. Moreover, the longer the core, the higher was the maximum intensity of the induced E: 2.5 V/m for 41 mm, 2.3 V/m for 23 mm, and 1.6 V/m for 9 mm. Figure 4 shows electric field components tangential (E tan ) and normal (E norm ) to the nerve fibers, for the two circular μCoils and the air core figure-of-eight μCoil. The figure-of-eight μCoil induced a tangential component that was five and seven times higher than that induced by a four-layered and a two-layered single circular μCoil, respectively, and a normal component that was two times lower.

Variation of the Peripheral Nerve Stimulation Threshold
When simulating exposure of the arm to the two-layered single circular coil, the computed stimulation threshold for the radial nerve fibers varied between 259 A and 269 A. Nerve action potentials were generated near the extremities of the coil. Increasing the number of layers from two to four reduced the coil current threshold in each layer/coil to values between 177.56 A and 184.39 A, while the activation site of the neuronal response was not affected. The stimulation threshold was further reduced when simulating with a figure-of-eight coil: the five times higher induced tangential electric field of the figure-of-eight-coil reduced the threshold to values between 58 A and 58.5 A; such a small range of variability may be due to the augmented focality. For all fibers, the activation site of the neuronal response was at the center of the figure-of-eight coil. Onset timing was not affected by the number of layers or by the change of coil shape and remained around an average time of 0.8 ms. A reduction of the onset timing from 0.8 to 0.5 ms, meaning an earlier spike, was obtained when changing from air to iron core. The presence of the core also reduced the stimulation threshold from 58 A to 14 A, 10.1 A, and 9.4 A for the short, medium, and long core, respectively, and further oriented E field along the direction of the nerve.

Discussion and Conclusion
The lack of focality of TMS stimulators is a significant barrier for the future deployment of this type of noninvasive neurostimulator in the clinics or neuroscience research. Several solutions have been proposed during the past years to overcome this limitation. For example, in 2010, Talebinejal and Musallam conducted a numerical study on a multilayer miniaturized coil made up of braided Litz wire [23] and demonstrated the feasibility of realizing a TMS coil reduced in size. Moreover, they showed that adding more layers increased magnetic field intensity while keeping the coil area constant. Following a similar geometry, in Tischler et al. 2011, a mini-coil with an outer diameter of 26.5 mm was developed. It is crucial to notice that all the attempts made to reduce dimensions of the TMS coils were conceived for either invasive [24] or noninvasive [23] animal models. To our knowledge, the newly developed figure-of-eight μCoil, based on flex circuit technology [32], is the first prototype of miniaturized coil manufactured for on-human applications. Experiments conducted on the peripheral nervous system of healthy volunteers showed that the μCoil could elicit somatosensory nerve action potentials (SNAPs) when the coil is placed over the superficial branch of the radial nerve. In this chapter, we conducted a numerical study that aimed to show how different possible μCoil geometries can affect the radial nerve stimulation threshold. Six different topologies of μCoils were modeled, based on the three different μCoils actually manufactured: the single circular four-layered and the figure-of-eight coils were originally designed to be tested on the PNS, and the iron core figure-of-eight coil developed to target the cerebral cortex. Each coil was placed over the radial nerve superficial branch of the neurofunctionalized model Yoon-Sun [27]. A multiphysics approach was herein considered, and electromagnetic simulations were coupled with the dynamic neuronal solution. To study the effect of adding more layers, we first studied a two-layered single circular μCoil and compared it with the four-layered one. The results showed that increasing the number of layers from two to four increased the intensity of the B field by 50%, thus inducing a stronger E field. As a consequence, the stimulation threshold was lower for the four-layered coil. A figure-of-eight coil was obtained by pairing together two four-layered circular coils. Such μCoil was able to generate a stronger magnetic field, thus inducing a stronger electric field on the nerve. Moreover, the figure-of-eight coil oriented the E field along the trajectory of the nerve fibers, leading to a better coupling between the coil and the nerve and a lower stimulation threshold. This coupling is further enhanced by addition of an iron core. Additionally, the results showed that the presence of the core not only increased intensity of the induced electric field, as expected, but also reduced the onset timing of the first spike. Nevertheless, it should be noted that saturation of the iron core was not considered in this study. Thus, computed stimulation threshold was underestimated, due to the overestimation of both B and E field intensities. However, we are planning to model next the nonlinearities of a real core. We confirmed that the figure-of-eight coil was more efficient than the circular coil as it enhanced the focality and the coupling between nerve fibers and E field. Moreover, the results showed that the presence of an iron core further reduced stimulation threshold while increasing depth, making this geometry suitable for investigating the μCoil stimulation of the brain cortex.
Open Access This chapter is licensed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license and indicate if changes were made.
The images or other third party material in this chapter are included in the chapter's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the chapter's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.