0580

Increasing Peripheral Nerve Stimulation Thresholds Using Gradient Array Coils
Reza Babaloo1,2, Manouchehr Takrimi2, and Ergin Atalar1,2
1Department of Electrical and Electronics Engineering, Bilkent University, Ankara, Turkey, 2National Magnetic Resonance Research Center (UMRAM), Bilkent University, Ankara, Turkey

Synopsis

Recent high-performance gradient coil designs are primarily limited by Peripheral Nerve Stimulation thresholds rather than hardware. Electrical fields induced by gradient switching may stimulate nerves within the human body if they exceed predefined thresholds. Gradient array coils can reduce the maximum electric field by optimizing feeding currents while generating the desired magnetic field profile within a specified region of interest. Another degree of freedom provided by array configuration is the ability to generate a flexible shaped region of interest, which allows for electrical fields manipulation to increase peripheral nerve stimulation thresholds.

Introduction

The use of high-performance gradient coils in MRI has become significantly limited due to Peripheral Nerve Stimulation (PNS). The rapid switching of strong gradients induces electric fields (E-fields) within the human body results in stimulating nerves. Head-insert gradient coils with smaller linear Region of Interest (ROI) were developed to increase PNS thresholds1-3. Despite these efforts, recent gradient design algorithms incorporate constraints on the E-fields leading to PNS-optimized gradient coils4,5; However, further optimization cannot be performed after coil fabrication. On the other hand, gradient array coil consists of multiple coils6-9, which can be independently driven by separate power amplifiers, are capable of generating tunable field profiles (B- and E-field) within a flexible ROI by optimizing the feeding waveforms. This work shows that E-field produced by a gradient coil can be significantly reduced when linearity is needed in a small ROI, demonstrating the power of a gradient array coil.

Methods

A z-gradient array coil was used to demonstrate the substantial reduction of the electric field. The coil comprises primary and shield coils, each with 12 pairs of wire bundles (10 wires per bundle) and diameters/heights of 24/30 and 30/35 cm, respectively. This coil generates gradient strength of 122 mT/m within the spherical ROI of 140 mm diameter, where the feeding currents RMS is 78.2A (linearity error less than 10%). By optimizing the feeding currents, the coil can also produce the same gradient strength within a customizable disk-shaped ROI. The disk position is adjustable, and the required spherical ROI for imaging can be covered by shifting the disk to the left and right of the isocenter. A conventional symmetric z-gradient coil (of the exact dimensions as the array coil) was designed with Sim4Life (Zurich Med Tech) to provide similar performance. First, we compare the electric fields induced by these two coils while keeping the same B-field profiles. E-fields of the disk-shaped ROI were then simulated in two positions: (1) at isocenter, (2) at z=-60mm. Simulations were performed using the low-frequency Magneto Quasi-Static solvers available in Sim4Life. Sinusoidal currents with dedicated peak amplitudes were applied through the coils at the frequency of 1 kHz. The magnetic fields were extracted from the vector potentials calculated using the Biot-Savart equation within the coils. The solver employs tricubic interpolation to compute the E-fields. We report maximum intensity projection of E-fields within the volume of an arm model (Yoon-sun Arm, Sim4Life). Winding patterns, the position of the arm model, and B-fields (BZ) are shown in Figure 1.

Results

Magnetic field profile simulations for symmetric conventional and array gradient coils are shown in Figures 1C and 1D, respectively. In a spherical ROI (140mm diameter), both of these coils generate nearly 120 mT/m gradient strength, but the maximum magnetic field is lower in an array configuration. Since the arm model is not homogeneous, the maximum E-field may occur within the volume rather than on the surface. Therefore, we demonstrate the maximum intensity projection of E-fields. Figure 2 shows the E-fields of conventional and array gradient coils. Different patterns of array winding result in significantly reduced E-field magnitudes compared to the conventional one (for the same gradient B-field), allowing for higher PNS thresholds. Alternatively, at a constant maximum E-field for both coils, it is possible to achieve higher gradient strengths and slew rates with the array coil without hitting nerves. Figure 3 depicts the B-fields and maximum E-field projections for a disk-shaped ROI at isocenter and z=-60mm. Because array coils can generate dynamic gradient fields by optimizing feeding currents, ROI dimensions can be changed after the coil is manufactured. As shown in Figure 3, shrinking the ROI (gradient strength remains at 120mT/m) significantly reduces the maximum E-field, even when compared to the array coil's functionality in spherical ROI. Therefore, much higher gradients can be generated in a small ROI without exceeding PNS thresholds, promising for applications like diffusion-weighted imaging. Moving the disk to different positions alters the E-field maps and may result in some hotspots outsides of ROI, as they are visible in Figure 3B; however, this can be resolved by considering the E-field thresholds during the array coil design.

Discussion

In this work, we compared E-field simulations of conventional symmetric and array gradient coils. The results indicate that the maximum E-field induced on the surface and inside of the arm model in the array configuration is much less than the maximum E-field induced by conventional coil, while the generated B-field profile and gradient strength are the same within the desired ROI. As a result, gradient array coils can be used with higher performances without surpassing PNS thresholds. The array design's ability to produce flexible magnetic fields also allows for different ROI sizes and shapes, resulting in further reduction in maximum E-field. Here, the array coil design and feeding currents optimization were performed subject to the gradient strength, ROI, and linearity error; however, adding the E-fields as a constraint to the optimization process will result in a PNS-optimized design.

Acknowledgements

The authors acknowledge “Sim4Life by ZMT, www.zurichmeditech.com” for providing an Academic License.

References

1. Zhang B, Yen YF, Chronik BA, McKinnon GC, Schaefer DJ, Rutt BK. Peripheral nerve stimulation properties of head and body gradient coils of various sizes. Magn Reson Med. 2003;50:50-58.

2. Tang F, Liu F, Freschi F, et al. An improved asymmetric gradient coil design for high-resolution MRI head imaging. Phys Med Biol. 2016;61:8875-8889.

3. Roemer PB, Wade T, Alejski A, McKenzie CA, Rutt BK. Electric field calculation and peripheral nerve stimulation prediction for head and body gradient coils. Magn Reson Med. 2021;86:2301-2315.

4. Davids M, Guerin B, Klein V, Wald LL. Optimization of MRI gradient coils with explicit peripheral nerve stimulation constraints. IEEE Trans Med Imaging. 2021;40:129-142.

5. Roemer PB, Rutt BK. Minimum electric-field gradient coil design: theoretical limits and practical guidelines. Magn Reson Med. 2021;86:569-580.

6. Ertan, K., Taraghinia, S., Sadeghi, A. and Atalar, E. (2018), A z‐gradient array for simultaneous multi-slice excitation with a single‐band RF pulse. Magn. Reson. Med, 80: 400-412.

7. Ertan, K, Taraghinia, S, Atalar, E. Driving mutually coupled gradient array coils in magnetic resonance imaging. Magn Reson Med. 2019; 82: 1187– 1198.

8. Babaloo, R., Taraghinia, S. Acikel, V., Takrimi, M., and Atalar, E. Digital Feedback Design for Mutual Coupling Compensation in Gradient Array System. Proceedings of the 28th Annual Meeting of ISMRM, 2020.

9. Takrimi, M. and Atalar, E. A Programmable Set of Z-Gradient Array and Active-Shield Array for Magnetic Resonance Imaging. Proceedings of the 28th Annual Meeting of ISMRM, 2020.

Figures

Figure 1. Winding patterns and the arm model positioning the hand at the isocenter are shown in A and B for conventional and array coils, respectively. C and D depict magnetic field profiles (Bz) on the xz-plane, which are very similar for both coils (120 mT/m gradient strength) within ROI indicated by a dashed red circle.

Figure 2. Maximum intensity projection of E-field magnitude (on xz and xy planes) are shown for conventional and array coils in the left and right columns, respectively. The color scale is the magnitude of the E-field in units of V/m. The array coil E-field shows a 14% reduction compared to the conventional one. Because the arm model is inhomogeneous, the maximum E-field occurs within the volume rather than on the surface.

Figure 3. B-field and E-field of array coil when optimized for a disk-shaped ROI at the isocenter (A) and z=-60mm (B). E-field maps show a significant reduction in maximum E-field magnitude for both optimizations, except for some points within the shifted ROI (B). The voxels size in the simulation process can also cause these hotspots, which can be eliminated by reducing voxels size.

Proc. Intl. Soc. Mag. Reson. Med. 30 (2022)
0580
DOI: https://doi.org/10.58530/2022/0580