Boundary element fast multipole method for modeling electrical brain stimulation with voltage and current electrodes.

Objective. To formulate, validate, and apply an alternative to the finite element method (FEM) high-resolution modeling technique for electrical brain stimulation-the boundary element fast multipole method (BEM-FMM). To include practical electrode models for both surface and embedded electrodes.Approach. Integral equations of the boundary element method in terms of surface charge density are combined with a general-purpose fast multipole method and are expanded for voltage, shunt, current, and floating electrodes. The solution of coupled and properly weighted/preconditioned integral equations is accompanied by enforcing global conservation laws: charge conservation law and Kirchhoff's current law.Main results.A sub-percent accuracy is reported as compared to the analytical solutions and simple validation geometries. Comparison to FEM considering realistic head models resulted in relative differences of the electric field magnitude in the range of 3%-6% or less. Quantities that contain higher order spatial derivatives, such as the activating function, are determined with a higher accuracy and a faster speed as compared to the FEM. The method can be easily combined with existing head modeling pipelines such as headreco or mri2mesh.Significance.The BEM-FMM does not rely on a volumetric mesh and is therefore particularly suitable for modeling some mesoscale problems with submillimeter (and possibly finer) resolution with high accuracy at moderate computational cost. Utilizing Helmholtz reciprocity principle makes it possible to expand the method to a solution of EEG forward problems with a very large number of cortical dipoles.

Vertical open-bore MRI scanners generate significantly less radiofrequency heating around implanted leads: A study of deep brain stimulation implants in 1.2T OASIS scanners versus 1.5T horizontal systems.

PURPOSE: Patients with active implants such as deep brain stimulation (DBS) devices are often denied access to MRI due to safety concerns associated with the radiofrequency (RF) heating of their electrodes. The majority of studies on RF heating of conductive implants have been performed in horizontal close-bore MRI scanners. Vertical MRI scanners which have a 90° rotated transmit coil generate fundamentally different electric and magnetic field distributions, yet very little is known about RF heating of implants in this class of scanners. We performed numerical simulations as well as phantom experiments to compare RF heating of DBS implants in a 1.2T vertical scanner (OASIS, Hitachi) compared to a 1.5T horizontal scanner (Aera, Siemens). METHODS: Simulations were performed on 90 lead models created from post-operative CT images of patients with DBS implants. Experiments were performed with wires and commercial DBS devices implanted in an anthropomorphic phantom. RESULTS: We found significant reduction of 0.1 g-averaged specific absorption rate (30-fold, P < 1 × 10-5 ) and RF heating (9-fold, P < .026) in the 1.2T vertical scanner compared to the 1.5T conventional scanner. CONCLUSION: Vertical MRI scanners appear to generate lower RF heating around DBS leads, providing potentially heightened safety or the flexibility to use sequences with higher power levels than on conventional systems.