Software Toolkit for Fast High-Resolution TMS Modeling

Abstract

Computational modeling of Transcranial Magnetic Stimulation (TMS) is a tradeoff between computational speed vs. spatial precision. In this study, we introduce a software toolkit for high-resolution TMS modeling, which may offer the best of both. The toolkit employs the recently developed boundary element fast multipole method (BEM-FMM) with accurate solution computations close to tissue interfaces. It operates within the MATLAB platform and is designed for a broad, medically-oriented computational research community. To enable easy, subject-specific operation, the package is compatible with human head models generated with automated tools such as SimNIBS and FreeSurfer. Both coil design in free space as well as efficacy and focality of TMS for a specific subject could be modeled by the package, including optional user-defined parametric loops. Detailed and widely used coil models, generated in both CAD and wire formats, may include several hundred thousand elementary current elements and observation spaces with approximately 1 M field points; the corresponding computational times are on the order of 1 sec. Detailed head models with approximately 1 M triangular facets and a mesh resolution of 0.6 points per mm2 are simulated in approximately 1.5 min which is arguably the fastest time to date. Further reduction of computational times is foreseen. The toolkit is augmented with a population of 16 ready-to-use head models for performing simulations for computational studies that do not involve MRI data collection. The toolkit is also augmented with a coil geometry generator capable of creating accurate coil models.