Dataset Information

Boundary Element Fast Multipole Method for Enhanced Modeling of Neurophysiological Recordings.

ABSTRACT:

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.

SUBMITTER: Makarov SN

PROVIDER: S-EPMC7704617 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Objective.In our recent work pertinent to modeling of brain stimulation and neurophysiological recordings, substantial modeling errors in the computed electric field and potential have sometimes been observed for standard multi-compartment head models. The goal of this study is to quantify those errors and, further, eliminate them through an adaptive mesh refinement (AMR) algorithm. The study concentrates on transcranial magnetic stimulation (TMS), transcranial electrical stimulation (TES), and electroencephalography (EEG) forward problems.Approach.We propose, describe, and systematically investigate an AMR method using the boundary element method with fast multipole acceleration (BEM-FMM) as the base numerical solver. The goal is to efficiently allocate additional unknowns to critical areas of the model, where they will best improve solution accuracy. The implemented AMR method's accuracy improvement is measured on head models constructed from 16 Human Connectome Project subjects under problem classes of TES, TMS, and EEG. Errors are computed between three solutions: an initial non-adaptive solution, a solution found after applying AMR with a conservative refinement rate, and a 'silver-standard' solution found by subsequent 4:1 global refinement of the adaptively-refined model.Main results.Excellent agreement is shown between the adaptively-refined and silver-standard solutions for standard head models. AMR is found to be vital for accurate modeling of TES and EEG forward problems for standard models: an increase of less than 25% (on average) in number of mesh elements for these problems, efficiently allocated by AMR, exposes electric field/potential errors exceeding 60% (on average) in the solution for the unrefined models.Significance.This error has especially important implications for TES dosing prediction-where the stimulation strength plays a central role-and for EEG lead fields. Though the specific form of the AMR method described here is implemented for the BEM-FMM, we expect that AMR is applicable and even required for accurate electromagnetic simulations by other numerical modeling packages as well.