1
|
Makaroff SN, Qi Z, Rachh M, Wartman WA, Weise K, Noetscher GM, Daneshzand M, Deng ZD, Greengard L, Nummenmaa AR. A fast direct solver for surface-based whole-head modeling of transcranial magnetic stimulation. Sci Rep 2023; 13:18657. [PMID: 37907689 PMCID: PMC10618282 DOI: 10.1038/s41598-023-45602-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2023] [Accepted: 10/21/2023] [Indexed: 11/02/2023] Open
Abstract
When modeling transcranial magnetic stimulation (TMS) in the brain, a fast and accurate electric field solver can support interactive neuronavigation tasks as well as comprehensive biophysical modeling. We formulate, test, and disseminate a direct (i.e., non-iterative) TMS solver that can accurately determine global TMS fields for any coil type everywhere in a high-resolution MRI-based surface model with ~ 200,000 or more arbitrarily selected observation points within approximately 5 s, with the solution time itself of 3 s. The solver is based on the boundary element fast multipole method (BEM-FMM), which incorporates the latest mathematical advancement in the theory of fast multipole methods-an FMM-based LU decomposition. This decomposition is specific to the head model and needs to be computed only once per subject. Moreover, the solver offers unlimited spatial numerical resolution. Despite the fast execution times, the present direct solution is numerically accurate for the default model resolution. In contrast, the widely used brain modeling software SimNIBS employs a first-order finite element method that necessitates additional mesh refinement, resulting in increased computational cost. However, excellent agreement between the two methods is observed for various practical test cases following mesh refinement, including a biophysical modeling task. The method can be readily applied to a wide range of TMS analyses involving multiple coil positions and orientations, including image-guided neuronavigation. It can even accommodate continuous variations in coil geometry, such as flexible H-type TMS coils. The FMM-LU direct solver is freely available to academic users.
Collapse
Affiliation(s)
- S N Makaroff
- Electrical and Computer Engineering Department, Worcester Polytechnic Institute, Worcester, MA, 01609, USA
- Athinoula A. Martinos Ctr. for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, 02129, USA
| | - Z Qi
- Electrical and Computer Engineering Department, Worcester Polytechnic Institute, Worcester, MA, 01609, USA.
| | - M Rachh
- Center for Computational Mathematics, Flatiron Institute, New York, NY, 10010, USA
| | - W A Wartman
- Electrical and Computer Engineering Department, Worcester Polytechnic Institute, Worcester, MA, 01609, USA
| | - K Weise
- Max Planck Institute for Human Cognitive and Brain Sciences, Stephanstr. 1a, 04103, Leipzig, Germany
- Advanced Electromagnetics Group, Technische Universität Ilmenau, Helmholtzplatz 2, 98693, Ilmenau, Germany
| | - G M Noetscher
- Electrical and Computer Engineering Department, Worcester Polytechnic Institute, Worcester, MA, 01609, USA
| | - M Daneshzand
- Athinoula A. Martinos Ctr. for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, 02129, USA
| | - Zhi-De Deng
- Computational Neurostimulation Research Program, Noninvasive Neuromodulation Unit, Experimental Therapeutics and Pathophysiology Branch, National Institute of Mental Health, NIH 10 Center Drive, Bethesda, MD, 20892, USA
| | - L Greengard
- Center for Computational Mathematics, Flatiron Institute, New York, NY, 10010, USA
- Courant Institute of Mathematical Sciences, 251 Mercer Street, New York, NY, 10012, USA
| | - A R Nummenmaa
- Athinoula A. Martinos Ctr. for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, 02129, USA
| |
Collapse
|
2
|
Makaroff SN, Qi Z, Rachh M, Wartman WA, Weise K, Noetscher GM, Daneshzand M, Deng ZD, Greengard L, Nummenmaa AR. A fast direct solver for surface-based whole-head modeling of transcranial magnetic stimulation. Res Sq 2023:rs.3.rs-3079433. [PMID: 37503106 PMCID: PMC10371170 DOI: 10.21203/rs.3.rs-3079433/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/29/2023]
Abstract
Background When modeling transcranial magnetic stimulation (TMS) in the brain, a fast and accurate electric field solver can support interactive neuronavigation tasks as well as comprehensive biophysical modeling. Objective We formulate, test, and disseminate a direct (i.e., non-iterative) TMS solver that can accurately determine global TMS fields for any coil type everywhere in a high-resolution MRI-based surface model with ~200,000 or more arbitrarily selected observation points within approximately 5 sec, with the solution time itself of 3 sec. Method The solver is based on the boundary element fast multipole method (BEM-FMM), which incorporates the latest mathematical advancement in the theory of fast multipole methods - an FMM-based LU decomposition. This decomposition is specific to the head model and needs to be computed only once per subject. Moreover, the solver offers unlimited spatial numerical resolution. Results Despite the fast execution times, the present direct solution is numerically accurate for the default model resolution. In contrast, the widely used brain modeling software SimNIBS employs a first-order finite element method that necessitates additional mesh refinement, resulting in increased computational cost. However, excellent agreement between the two methods is observed for various practical test cases following mesh refinement, including a biophysical modeling task. Conclusion The method can be readily applied to a wide range of TMS analyses involving multiple coil positions and orientations, including image-guided neuronavigation. It can even accommodate continuous variations in coil geometry, such as flexible H-type TMS coils. The FMM-LU direct solver is freely available to academic users.
Collapse
Affiliation(s)
- S N Makaroff
- Electrical and Computer Engineering Department, Worcester Polytechnic Institute, Worcester, MA 01609 USA
- Athinoula A. Martinos Ctr. for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA 02129 USA
| | - Z Qi
- Electrical and Computer Engineering Department, Worcester Polytechnic Institute, Worcester, MA 01609 USA
| | - M Rachh
- Center for Computational Mathematics, Flatiron Institute, New York, NY 10010 USA
| | - W A Wartman
- Electrical and Computer Engineering Department, Worcester Polytechnic Institute, Worcester, MA 01609 USA
| | - K Weise
- Max Planck Institute for Human Cognitive and Brain Sciences, Stephanstr. 1a, 04103, Leipzig Germany
- Technische Universität Ilmenau, Advanced Electromagnetics Group, Helmholtzplatz 2, 98693 Ilmenau Germany
| | - G M Noetscher
- Electrical and Computer Engineering Department, Worcester Polytechnic Institute, Worcester, MA 01609 USA
| | - M Daneshzand
- Athinoula A. Martinos Ctr. for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA 02129 USA
| | - Zhi-De Deng
- Computational Neurostimulation Research Program, Noninvasive Neuromodulation Unit, Experimental Therapeutics & Pathophysiology Branch, National Institute of Mental Health, NIH 10 Center Drive, Bethesda, MD 20892 USA
| | - L Greengard
- Center for Computational Mathematics, Flatiron Institute, New York, NY 10010 USA
- Courant Institute of Mathematical Sciences, 251 Mercer Street, New York, NY 10012 USA
| | - A R Nummenmaa
- Athinoula A. Martinos Ctr. for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA 02129 USA
| |
Collapse
|