The multipole approach for EEG forward modeling using the finite element method

For accurate EEG forward solutions, it is necessary to apply numerical methods that allow to take into account the realistic geometry of the subject's head. A commonly used method to solve this task is the finite element method (FEM). Different approaches have been developed to obtain EEG forward solutions for dipolar sources with the FEM. The St. Venant approach is frequently applied, since its high numerical accuracy and stability as well as its computational efficiency was demonstrated in multiple comparison studies. In this manuscript, we propose a variation of the St. Venant approach, the multipole approach, to improve the numerical accuracy of the St. Venant approach even further and to allow for the simulation of additional source scenarios, such as quadrupolar sources. Exploiting the multipole expansion of electric fields, we demonstrate that the newly proposed multipole approach achieves even higher numerical accuracies than the St. Venant approach in both multi-layer sphere and realistic head models. Additionally, we exemplarily show that the multipole approach allows to not only simulate dipolar but also quadrupolar sources.