Frontiers in Human Neuroscience (Aug 2023)

CutFEM forward modeling for EEG source analysis

  • Tim Erdbrügger,
  • Tim Erdbrügger,
  • Andreas Westhoff,
  • Malte Höltershinken,
  • Malte Höltershinken,
  • Jan-Ole Radecke,
  • Jan-Ole Radecke,
  • Yvonne Buschermöhle,
  • Yvonne Buschermöhle,
  • Alena Buyx,
  • Fabrice Wallois,
  • Sampsa Pursiainen,
  • Joachim Gross,
  • Joachim Gross,
  • Rebekka Lencer,
  • Rebekka Lencer,
  • Rebekka Lencer,
  • Rebekka Lencer,
  • Christian Engwer,
  • Carsten Wolters,
  • Carsten Wolters

DOI
https://doi.org/10.3389/fnhum.2023.1216758
Journal volume & issue
Vol. 17

Abstract

Read online

IntroductionSource analysis of Electroencephalography (EEG) data requires the computation of the scalp potential induced by current sources in the brain. This so-called EEG forward problem is based on an accurate estimation of the volume conduction effects in the human head, represented by a partial differential equation which can be solved using the finite element method (FEM). FEM offers flexibility when modeling anisotropic tissue conductivities but requires a volumetric discretization, a mesh, of the head domain. Structured hexahedral meshes are easy to create in an automatic fashion, while tetrahedral meshes are better suited to model curved geometries. Tetrahedral meshes, thus, offer better accuracy but are more difficult to create.MethodsWe introduce CutFEM for EEG forward simulations to integrate the strengths of hexahedra and tetrahedra. It belongs to the family of unfitted finite element methods, decoupling mesh and geometry representation. Following a description of the method, we will employ CutFEM in both controlled spherical scenarios and the reconstruction of somatosensory-evoked potentials.ResultsCutFEM outperforms competing FEM approaches with regard to numerical accuracy, memory consumption, and computational speed while being able to mesh arbitrarily touching compartments.DiscussionCutFEM balances numerical accuracy, computational efficiency, and a smooth approximation of complex geometries that has previously not been available in FEM-based EEG forward modeling.

Keywords