MEG forward problem formulation using equivalent surface current densities

We present a formulation for the magnetoencephalography (MEG) forward problem with a layered head model. Traditionally the magnetic field is computed based on the electric potential on the interfaces between the layers. We propose to express the effect of the volumetric currents in terms of an equivalent surface current density on each interface, and obtain the magnetic field based on them. The boundary elements method is used to compute the equivalent current density and the magnetic field for a realistic head geometry. We present numerical results showing that the MEG forward problem is solved correctly with this formulation, and compare it with the performance of the traditional formulation. We conclude that the traditional formulation generally performs better, but still the new formulation is useful in certain situations.