0408

Fast PNS characterization of MRI gradient coils using a Huygens’ PNS model: Application to multiple patient positions and orientations
Mathias Davids1,2,3, Bastien Guerin1,2, and Lawrence L Wald1,2,4
1Martinos Center for Biomedical Imaging, Boston, MA, United States, 2Harvard Medical School, Boston, MA, United States, 3Computer Assisted Clinical Medicine, Mannheim, Germany, 4Harvard-MIT, Division of Health Sciences and Technology, Cambridge, MA, United States

Synopsis

We apply a fast PNS model evaluated on a Huygens’ surface to characterize PNS performance as a function of patient position and pose in previously described PNS optimized gradients. The coils were initially PNS-optimized for a female body model in a standard brain imaging position. The Huygens’ approach allows us to assess other body positions and demonstrate the PNS benefits/penalties associated with other imaging applications and in both a female and male body model by dramatically speeding up the PNS characterization (to a couple of seconds per body position/orientation). The findings support a broad benefit from the PNS optimized windings.

Target audience

MRI gradient designers and safety researchers

Purpose

Peripheral Nerve Stimulation (PNS) significantly limits the usability of the latest generation gradient coils in humans [1-3] where rapid gradient slewing induces strong E-fields that can excite nerves. We recently applied PNS optimization to whole-body and head-only gradients [4-5] using a single body model (adult female) and single body position (head-first supine brain imaging). This approach yielded up to 2-fold PNS improvements at the cost of moderate inductance and linearity penalties. It was not clear, however, that this improvement extends to non-head imaging applications. Here we apply a Huygens’ surface approach [6] to enable systematic characterization of PNS performance in a wide range of body positions/orientations and an additional body model in a matter of seconds.

Methods

Huygens P-matrix approach: The PNS analysis was performed using a linear PNS matrix (P-matrix) for a current basis function set defined on a Huygens’ surface tightly enclosing one of our body models [7,8]. The PNS responses are precomputed for all basis functions on the Huygens’ surface (small current loops, Fig. 1A). For each current basis, we simulate the induced E-fields and extract the linear PNS oracle values [9] (reciprocal PNS thresholds) along all nerves (Fig. 1A/B). After precomputation and assembly of the full Huygens’ P-matrix (Fig. 1C) which takes a couple of days per body model, the Huygens’ P-matrix can be mapped to a new discrete coil winding or coil geometry and patient position/pose, yielding a coil-specific P-matrix. This mapping only takes seconds, as it relies on fast Biot-Savart calculations of the magnetic fields and solving a small linear system of equations.
Studied coil solutions and body positions: We analyze previously described Y-axis body gradients designed with 5% field non-linearity [5]. We focus our analysis on one conventionally designed coil without PNS optimization (coil B1), and one coil that was optimized with a PNS constraint from the female model in a head-imaging position (coil B2). We evaluate these coils in both female and male models in both head-first and feet-first body orientations for a variety of body z-positions ranging between −40 cm to +120 cm in 5 cm steps (z = 0 cm represents head at isocenter).

Results

Figure 2 (top) shows the tradeoff L-curve between the maximum PNS oracle (reciprocal PNS threshold) and coil inductance for body Y-axis gradients optimized for the female model for head imaging (head-first supine). Every point on the L-curve is a winding solution, including the two coils B1 (without PNS optimization) and B2 (with PNS optimization). Figure 2 (bottom) extends the PNS analysis to show how the coils on the L-curve perform for other patient positions. The additional dimension corresponds to a new z-position of the female model in the coils. Optimization of the coils for head imaging (head-first supine) led to 35% PNS oracle improvements (55% higher thresholds) for that body position, but also improved PNS characteristics for cardiac imaging (25% PNS oracle reduction, 35% increased PNS threshold).
Figure 3 shows a more detailed PNS threshold characterization of the unoptimized coil B1 (solid curve) and the optimized coil B2 (dashed curve) as a function of position and orientation (including feet-first). Regions where the solid line lies above the dashed line indicate use-cases that retain some benefit from the PNS optimized design B2. The color corresponds to the origin of the nerve activation in the body (blue: head/chest; red: abdomen/pelvis).
The optimized coil B2 achieves higher PNS thresholds for head/cardiac imaging and lower thresholds for abdominal/pelvic/knee imaging. Note that the origin of PNS changes multiple times between head/chest (blue ROI) and abdominal/pelvic activation (red ROI). The PNS performance for abdominal/pelvic/knee imaging is improved by switching to the more common feet-first orientation (Fig. 3, bottom), at the cost of PNS degradation for head/cardiac imaging (which would not be done feet-first). Figure 4 shows a similar analysis for the male body model. Here, using coil B2 with the head-first orientation showed PNS improvements for head/cardiac imaging (+90% and +5%, resp.), but degraded PNS performance for head-first abdominal/pelvic/knee imaging (thresholds lowered by 30%, 10%, and 14%). These results are qualitatively similar to those found in the female model, although PNS thresholds varied more strongly with body position for the male model.
Figure 5 shows maps of PNS hot-spots (sites of strongest activation) in the female model for both coils B1 and B2 at common clinically used body locations and orientations. The distribution of activation patterns was similar in many cases, with typical activation sites in the shoulders, arms, intercostal nerves, neck, facial and pelvic area nerves. In all cases (except for the head-first abdominal imaging case), the optimized coil B2 improved PNS performance (increased PNS thresholds) compared to the unoptimized coil B1, even though this coil was only optimized for head-first supine head imaging.

Conclusion

We applied the recently described Huygens’ P-matrix approach to characterize whole-body gradients for varying scan positions and body orientations. The data demonstrates that gradient coils optimized for a specific scan position also improve PNS for other imaging applications, underscoring a robust ability of PNS optimized gradient coil windings to provide PNS improvement in a wide range of imaging applications.

Acknowledgements

The authors would like to acknowledge the help of past and present members of the gradient coil group at Siemens Healthineers, including Peter Dietz, Gudrun Ruyters, Axel vom Endt, Ralph Kimmlingen, Eva Eberlein, and Franz Hebrank. Research reported in this publication was supported by the National Institute of Biomedical Imaging and Bioengineering, and the National Institute for Mental Health of the National Institutes of Health under award numbers R24MH106053, U01EB026996, U01EB025162, U01EB025121, R01EB028250. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

References

[1] Setsompop et al., “Pushing the limits of in vivo diffusion MRI for the human connectome project”, NeuroImage, vol. 80, pp. 220 – 233, 2013

[2] McNab et al., “The human connectome project and beyond: Initial applications of 300 mT/m gradients”, NeuroImage, vol. 80, pp. 234 – 245, 2013

[3] Tan et al., “Peripheral nerve stimulation limits of a high amplitude and slew rate magnetic field gradient coil for neuroimaging”, Magnetic resonance in medicine, vol. 83, no. 1, pp. 352–366, 2020.

[4] Davids et al., “Peripheral Nerve Stimulation (PNS) constrained gradient coil design within a Boundary Element Method Stream Function (BEM-SF) optimization”, Proceedings of the 27th Annual Meeting of ISMRM, Montreal, Canada, 2019

[5] Davids et al., “Optimization of MRI Gradient Coils with Explicit Peripheral Nerve Stimulation Constraints”, IEEE Transactions on Medical Imaging (in print), 2020

[6] Davids et al., “Assembly of a PNS predicting "P-matrix" on a Huygens' surface for rapid PNS assessment of 2D or 3D gradient coil windings", Proceedings of the 28th Annual Meeting of ISMRM, 2020

[7] Davids et al. “Predicting magnetostimulation thresholds in the peripheral nervous system using realistic body models”, Sci Rep, vol. 7, no. 1, p. 5316, 2017.

[8] Davids et al. "Prediction of peripheral nerve stimulation thresholds of MRI gradient coils using coupled electromagnetic and neurodynamic simulations", Magnetic resonance in medicine vol. 81, no. 1 p. 686-701, 2019

[9] Davids et al. “Optimizing selective stimulation of peripheral nerves with arrays of coils or surface electrodes using a linear peripheral nerve stimulation metric”, J Neural Eng, vol. 17, no. 1, p. 016029, 2020.


Figures

Figure 1: A: Definition of the female model Huygens’ surface populated with 2497 magnetic basis functions. B: E-field induced by single basis near the heart. C: After performing E-field simulations and extracting PNS oracle values for all basis functions, we assemble the Huygens’ P-matrix linking all 2497 bases (columns) to PNS responses of all nerves (rows). D: PNS oracle for single 0.1 mm segment of the cauda equina mapped onto the Huygens’ surface (highlighted row in P-matrix). E: PNS oracle of all nerve segments for a single Huygens’ basis function (highlighted column in P-matrix).

Figure 2: Top: L-curve for optimized Y-axis gradients showing tradeoff between PNS oracle (reciprocal PNS threshold) and coil inductance. The coils were optimized for the female model for head imaging (head-first supine). Every circle corresponds to a coil solution, including the coils marked B1 and B2. Bottom: In order to assess the PNS performance of solutions B1 and B2 for non-head imaging, the L-curve is extended to other body positions within the coils, resulting in the shown surface. The PNS thresholds for coils B1 and B2 as a function of body position are shown by dashed lines.

Figure 3: Analysis of coils B1 and B2 (two dashed curves in Fig. 2, bottom), shown as PNS thresholds as a function of body position of the female model. This describes how the 55% PNS improvement of coil B2 (which was optimized for head at iso-center) changes for other body positions and poses. The color refers to the location of the activation. Regions where the solid line (coil B2 thresholds) lies above the dashed line (non-PNS constrained coil B1) indicate use-cases that retain some benefit from the optimized coil even though optimization only considered a head-at-isocenter position.

Figure 4: Results from the male model for the analysis shown in Fig. 3. Use-cases where the optimized coil B2 coil yields higher PNS thresholds than coil B1 are similar to the results in the female model: PNS thresholds increase for head-first supine head and cardiac imaging (top graph) and for feet-first supine abdominal, pelvic, and knee imaging (bottom graph).

Figure 5: PNS hot-spots in the female model expected for different body positions in the coils B1 (top) and B2 (bottom). Colored spheres show hot-spots with largest PNS oracle (smallest PNS threshold). For abdominal imaging, we show both head-first and feet-first supine poses; all other scan positions use either head-first (head/cardiac imaging) or feet-first pose (pelvic/knee imaging). In all cases, meant to correspond to conventional clinical patient positions, the optimized coil retains some value in raising PNS thresholds, except for the head-first supine abdominal imaging.

Proc. Intl. Soc. Mag. Reson. Med. 29 (2021)
0408