The octapolic ellipsoidal term in magnetoencephalography

The forward problem of magnetoencephalography (MEG) in ellipsoidal geometry has been studied by Dassios and Kariotou ["Magnetoencephalography in ellipsoidal geometry," J. Math. Phys. 44, 220 (2003)] using the theory of ellipsoidal harmonics. In fact, the analytic solution of the quadrupolic term for the magnetic induction field has been calculated in the case of a dipolar neuronal current. Nevertheless, since the quadrupolic term is only the leading nonvanishing term in the multipole expansion of the magnetic field, it contains not enough information for the construction of an effective algorithm to solve the inverse MEG problem, i.e., to recover the position and the orientation of a dipole from measurements of the magnetic field outside the head. For this task, the next multipole of the magnetic field is also needed. The present work provides exactly this octapolic contribution of the dipolar current to the expansion of the magnetic induction field. The octapolic term is expressed in terms of the ellipsoidal harmonics of the third degree, and therefore it provides the highest order terms that can be expressed in closed form using long but reasonable analytic and algebraic manipulations. In principle, the knowledge of the quadrupolic and the octapolic terms is enough to solve the inverse problem of identifying a dipole inside an ellipsoid. Nevertheless, a simple inversion algorithm for this problem is not yet known