Numerical simulation of electromagnetic, thermal, and mechanical responses of the human body to different stimuli in magnetic resonance imaging safety, antenna research, electromagnetic tomography, and electromagnetic stimulation is currently limited by the availability of anatomically adequate and numerically efficient cross-platform computational models or “virtual humans.” The objective of this study is to provide a comprehensive review of modern human models and body region models available in the field and their important features.

Objective: We develop a new accurate version of the boundary element fast multipole method for transcranial magnetic stimulation (TMS) related problems. This method is based on the surface-charge formulation and is using the highly efficient fast multipole accelerator along with analytical computations of neighbor surface integrals. Results: The method accuracy is demonstrated by comparison with the proven commercial finite-element method (FEM) software ANSYS Maxwell 18.2 2017 operating on unstructured grids and with adaptive mesh refinement. Five realistic high-definition head models from the Population Head Repository (IT'IS Foundation, Switzerland) have been acquired and augmented with a commercial TMS coil model (MRi-B91, MagVenture, Denmark). For each head model, simulations with our method and simulations with the FEM software ANSYS Maxwell 18.2 2017 have been performed. These simulations have been compared with each other and an excellent agreement was established in every case. Significance: At the same time, our new method runs approximately 500 times faster than the ANSYS FEM, finishes in about 200 s on a standard server, and naturally provides a submillimeter field resolution, which is justified using mesh refinement. Conclusions: Our method can be applied to modeling of brain stimulation and recording technologies such as TMS and magnetoencephalography, and has the potential to become a real-time high-resolution simulation tool.

Objective: To present and disseminate our transcranial magnetic stimulation (TMS) modeling software toolkit, including several new algorithmic developments, and to apply this software to realistic TMS modeling scenarios given a high-resolution model of the human head including cortical geometry and an accurate coil model. Approach: The recently developed charge-based boundary element fast multipole method (BEM-FMM) is employed as an alternative to the 1st order finite element method (FEM) most commonly used today. The BEM-FMM approach provides high accuracy and unconstrained numerical field resolution close to and across cortical interfaces. Here, the previously proposed BEM-FMM algorithm has been improved in several novel ways. Main results.: The improvements resulted in a threefold increase in computational speed while maintaining the same solution accuracy. The computational code based on the MATLAB ^® platform is made available to all interested researchers, along with a coil model repository and examples to create custom coils, head model repository, and supporting documentation. The presented software toolkit may be useful for post-hoc analyses of navigated TMS data using high-resolution subject-specific head models as well as accurate and fast modeling for the purposes of TMS coil/hardware development. Significance.: TMS is currently the only non-invasive neurostimulation modality that enables painless and safe supra-threshold stimulation by employing electromagnetic induction to efficiently penetrate the skull. Accurate, fast, and high resolution modeling of the electric fields may significantly improve individualized targeting and dosing of TMS and therefore enhance the efficiency of existing clinical protocols as well as help establish new application domains.

Objective . A study pertinent to the numerical modeling of cortical neurostimulation is conducted in an effort to compare the performance of the finite element method (FEM) and an original formulation of the boundary element fast multipole method (BEM-FMM) at matched computational performance metrics. Approach . We consider two problems: (i) a canonic multi-sphere geometry and an external magnetic-dipole excitation where the analytical solution is available and; (ii) a problem with realistic head models excited by a realistic coil geometry. In the first case, the FEM algorithm tested is a fast open-source getDP solver running within the SimNIBS 2.1.1 environment. In the second case, a high-end commercial FEM software package ANSYS Maxwell 3D is used. The BEM-FMM method runs in the MATLAB ^® 2018a environment. Main results . In the first case, we observe that the BEM-FMM algorithm gives a smaller solution error for all mesh resolutions and runs significantly faster for high-resolution meshes when the number of triangular facets exceeds approximately 0.25 M. We present other relevant simulation results such as volumetric mesh generation times for the FEM, time necessary to compute the potential integrals for the BEM-FMM, and solution performance metrics for different hardware/operating system combinations. In the second case, we observe an excellent agreement for electric field distribution across different cranium compartments and, at the same time, a speed improvement of three orders of magnitude when the BEM-FMM algorithm used. Significance . This study may provide a justification for anticipated use of the BEM-FMM algorithm for high-resolution realistic transcranial magnetic stimulation scenarios.

When modeling transcranial electrical stimulation (TES) and transcranial magnetic stimulation (TMS) in the brain, the meninges - dura, arachnoid, and pia mater - are often neglected due to high computational costs. We investigate the impact of the meningeal layers on the cortical electric field in TES and TMS while considering the headreco segmentation as the base model. We use T1/T2 MRI data from 16 subjects and apply the boundary element fast multipole method with adaptive mesh refinement, which enables us to accurately solve this problem and establish method convergence at reasonable computational cost. We compare electric fields in the presence and absence of various meninges for two brain areas (M1_HAND and DLPFC) and for several distinct TES and TMS setups. Maximum electric fields in the cortex for focal TES consistently increase by approximately 30% on average when the meninges are present in the CSF volume. Their effect on the maximum field can be emulated by reducing the CSF conductivity from 1.65 S/m to approximately 0.85 S/m. In stark contrast to that, the TMS electric fields in the cortex are only weakly affected by the meningeal layers and slightly (∼6%) decrease on average when the meninges are included. Our results quantify the influence of the meninges on the cortical TES and TMS electric fields. Both focal TES and TMS results are very consistent. The focal TES results are also in a good agreement with a prior relevant study. The solver and the mesh generator for the meningeal layers (compatible with SimNIBS) are available online.

Objective: A new numerical modeling approach is proposed which provides forward-problem solutions for both noninvasive recordings (EEG/MEG) and higher-resolution intracranial recordings (iEEG). Methods: The algorithm is our recently developed boundary element fast multipole method or BEM-FMM. It is based on the integration of the boundary element formulation in terms of surface charge density and the fast multipole method originating from its inventors. The algorithm still possesses the major advantage of the conventional BEM - high speed - but is simultaneously capable of processing a very large number of surface-based unknowns. As a result, an unprecedented spatial resolution could be achieved, which enables multiscale modeling. Results: For non-invasive EEG/MEG, we are able to accurately solve the forward problem with approximately 1 mm anatomical resolution in the cortex within 1-2 min given several thousand cortical dipoles. Targeting high-resolution iEEG, we are able to compute, for the first time, an integrated electromagnetic response for an ensemble (2,450) of tightly packed realistic pyramidal neocortical neurons in a full-head model with 0.6 mm anatomical cortical resolution. The neuronal arbor is comprised of 5.9 M elementary 1.2 μm long dipoles. On a standard server, the computations require about 5 min. Conclusion: Our results indicate that the BEM-FMM approach may be well suited to support numerical multiscale modeling pertinent to modern high-resolution and submillimeter iEEG. Significance: Based on the speed and ease of implementation, this new algorithm represents a method that will greatly facilitate simulations at multi-scale across a variety of applications.