Fast and accurate computational E-field dosimetry for group-level transcranial magnetic stimulation targeting

Transcranial magnetic stimulation (TMS) is used to study brain function and treat mental health disorders. During TMS, a coil placed on the scalp induces an E-field in the brain that modulates its activity. TMS is known to stimulate regions that are exposed to a large E-field. Clinical TMS protocols prescribe a coil placement based on scalp landmarks. There are inter-individual variations in brain anatomy that result in variations in the TMS-induced E-field at the targeted region and its outcome. These variations across individuals could in principle be minimized by developing a large database of head subjects and determining scalp landmarks that maximize E-field at the targeted brain region while minimizing its variation using computational methods. However, this approach requires repeated execution of a computational method to determine the E-field induced in the brain for a large number of subjects and coil placements. We developed a probabilistic matrix decomposition-based approach for rapidly evaluating the E-field induced during TMS for a large number of coil placements due to a pre-defined coil model. Our approach can determine the E-field induced in over 1 Million coil placements in 9.5 h, in contrast, to over 5 years using a brute-force approach. After the initial set-up stage, the E-field can be predicted over the whole brain within 2-3 ms and to 2% accuracy. We tested our approach in over 200 subjects and achieved an error of <2% in most and <3.5% in all subjects. We will present several examples of bench-marking analysis for our tool in terms of accuracy and speed. Furthermore, we will show the methods' applicability for group-level optimization of coil placement for illustration purposes only. The software implementation link is provided in the appendix.