Rapid computation of TMS-induced E-fields using a dipole-based magnetic stimulation profile approach

• Published in: NeuroImage (Orlando, Fla.) Vol. 237; p. 118097
• Main Authors: Daneshzand, Mohammad, Makarov, Sergey N., de Lara, Lucia I. Navarro, Guerin, Bastien, McNab, Jennifer, Rosen, Bruce R., Hämäläinen, Matti S., Raij, Tommi, Nummenmaa, Aapo
• Format: Journal Article
• Language: English
• Published: United States Elsevier Inc 15.08.2021; Elsevier Limited; Elsevier
• Subjects: Brain research; Conductivity; Dipole basis functions; Dosage; Dosing; Electric fields; Magnetic fields; Magnetic resonance imaging; Magnetic stimulation profile; Neuronavigation; Spatial distribution; Targeting; Tissues; Transcranial magnetic stimulation; Dipole basis functions; Dosing; Magnetic stimulation profile; Transcranial magnetic stimulation; Targeting; Neuronavigation
• ISSN: 1053-8119; 1095-9572; 1095-9572
• DOI: 10.1016/j.neuroimage.2021.118097

TMS neuronavigation with on-line display of the induced electric field (E-field) has the potential to improve quantitative targeting and dosing of stimulation, but present commercially available solutions are limited by simplified approximations. Developing a near real-time method for accurate approximation of TMS induced E-fields with subject-specific high-resolution surface-based head models that can be utilized for TMS navigation. Magnetic dipoles are placed on a closed surface enclosing an MRI-based head model of the subject to define a set of basis functions for the incident and total E-fields that define the subject's Magnetic Stimulation Profile (MSP). The near real-time speed is achieved by recognizing that the total E-field of the coil only depends on the incident E-field and the conductivity boundary geometry. The total E-field for any coil position can be obtained by matching the incident field of the stationary dipole basis set with the incident E-field of the moving coil and applying the same basis coefficients to the total E-field basis functions. Comparison of the MSP-based approximation with an established TMS solver shows great agreement in the E-field amplitude (relative maximum error around 5%) and the spatial distribution patterns (correlation >98%). Computation of the E-field took ~100 ms on a cortical surface mesh with 120k facets. The numerical accuracy and speed of the MSP approximation method make it well suited for a wide range of computational tasks including interactive planning, targeting, dosing, and visualization of the intracranial E-fields for near real-time guidance of coil positioning.