Calibration and Localization of Optically Pumped Magnetometers Using Electromagnetic Coils.

In this paper, we propose a method to estimate the position, orientation, and gain of a magnetic field sensor using a set of (large) electromagnetic coils. We apply the method for calibrating an array of optically pumped magnetometers (OPMs) for magnetoencephalography (MEG). We first measure the magnetic fields of the coils at multiple known positions using a well-calibrated triaxial magnetometer, and model these discreetly sampled fields using vector spherical harmonics (VSH) functions. We then localize and calibrate an OPM by minimizing the sum of squared errors between the model signals and the OPM responses to the coil fields. We show that by using homogeneous and first-order gradient fields, the OPM sensor parameters (gain, position, and orientation) can be obtained from a set of linear equations with pseudo-inverses of two matrices. The currents that should be applied to the coils for approximating these low-order field components can be determined based on the VSH models. Computationally simple initial estimates of the OPM sensor parameters follow. As a first test of the method, we placed a fluxgate magnetometer at multiple positions and estimated the RMS position, orientation, and gain errors of the method to be 1.0 mm, 0.2°, and 0.8%, respectively. Lastly, we calibrated a 48-channel OPM array. The accuracy of the OPM calibration was tested by using the OPM array to localize magnetic dipoles in a phantom, which resulted in an average dipole position error of 3.3 mm. The results demonstrate the feasibility of using electromagnetic coils to calibrate and localize OPMs for MEG.

Spatial sampling of MEG and EEG based on generalized spatial-frequency analysis and optimal design.

In this paper, we analyze spatial sampling of electro- (EEG) and magnetoencephalography (MEG), where the electric or magnetic field is typically sampled on a curved surface such as the scalp. By simulating fields originating from a representative adult-male head, we study the spatial-frequency content in EEG as well as in on- and off-scalp MEG. This analysis suggests that on-scalp MEG, off-scalp MEG and EEG can benefit from up to 280, 90 and 110 spatial samples, respectively. In addition, we suggest a new approach to obtain sensor locations that are optimal with respect to prior assumptions. The approach also allows to control, e.g., the uniformity of the sensor locations. Based on our simulations, we argue that for a low number of spatial samples, model-informed non-uniform sampling can be beneficial. For a large number of samples, uniform sampling grids yield nearly the same total information as the model-informed grids.

Potential of on-scalp MEG: Robust detection of human visual gamma-band responses.

Electrophysiological signals recorded intracranially show rich frequency content spanning from near-DC to hundreds of hertz. Noninvasive electromagnetic signals measured with electroencephalography (EEG) or magnetoencephalography (MEG) typically contain less signal power in high frequencies than invasive recordings. Particularly, noninvasive detection of gamma-band activity (>30 Hz) is challenging since coherently active source areas are small at such frequencies and the available imaging methods have limited spatial resolution. Compared to EEG and conventional SQUID-based MEG, on-scalp MEG should provide substantially improved spatial resolution, making it an attractive method for detecting gamma-band activity. Using an on-scalp array comprised of eight optically pumped magnetometers (OPMs) and a conventional whole-head SQUID array, we measured responses to a dynamic visual stimulus known to elicit strong gamma-band responses. OPMs had substantially higher signal power than SQUIDs, and had a slightly larger relative gamma-power increase over the baseline. With only eight OPMs, we could obtain gamma-activity source estimates comparable to those of SQUIDs at the group level. Our results show the feasibility of OPMs to measure gamma-band activity. To further facilitate the noninvasive detection of gamma-band activity, the on-scalp OPM arrays should be optimized with respect to sensor noise, the number of sensors and intersensor spacing.

On-scalp MEG system utilizing an actively shielded array of optically-pumped magnetometers.

The spatial resolution of magnetoencephalography (MEG) can be increased from that of conventional SQUID-based systems by employing on-scalp sensor arrays of e.g. optically-pumped magnetometers (OPMs). However, OPMs reach sufficient sensitivity for neuromagnetic measurements only when operated in a very low absolute magnetic field of few nanoteslas or less, usually not reached in a typical magnetically shielded room constructed for SQUID-based MEG. Moreover, field drifts affect the calibration of OPMs. Static and dynamic suppression of interfering fields is thus necessary for good-quality neuromagnetic measurements with OPMs. Here, we describe an on-scalp MEG system that utilizes OPMs and external compensation coils that provide static and dynamic shielding against ambient fields. In a conventional two-layer magnetically shielded room, our coil system reduced the maximum remanent DC-field component within an 8-channel OPM array from 70 to less than 1 nT, enabling the sensors to operate in the sensitive spin exchange relaxation-free regime. When compensating field drifts below 4 Hz, a low-frequency shielding factor of 22 dB was achieved, which reduced the peak-to-peak drift from 1.3 to 0.4 nT and thereby the standard deviation of the sensor calibration from 1.7% to 0.5%. Without band-limiting the field that was compensated, a low-frequency shielding factor of 43 dB was achieved. We validated the system by measuring brain responses to electric stimulation of the median nerve. With dynamic shielding and digital interference suppression methods, single-trial somatosensory evoked responses could be detected. Our results advance the deployment of OPM-based on-scalp MEG in lighter magnetic shields.

Adaptive neural network classifier for decoding MEG signals.

We introduce two Convolutional Neural Network (CNN) classifiers optimized for inferring brain states from magnetoencephalographic (MEG) measurements. Network design follows a generative model of the electromagnetic (EEG and MEG) brain signals allowing explorative analysis of neural sources informing classification. The proposed networks outperform traditional classifiers as well as more complex neural networks when decoding evoked and induced responses to different stimuli across subjects. Importantly, these models can successfully generalize to new subjects in real-time classification enabling more efficient brain-computer interfaces (BCI).

Requirements for Coregistration Accuracy in On-Scalp MEG.

Recent advances in magnetic sensing has made on-scalp magnetoencephalography (MEG) possible. In particular, optically-pumped magnetometers (OPMs) have reached sensitivity levels that enable their use in MEG. In contrast to the SQUID sensors used in current MEG systems, OPMs do not require cryogenic cooling and can thus be placed within millimetres from the head, enabling the construction of sensor arrays that conform to the shape of an individual's head. To properly estimate the location of neural sources within the brain, one must accurately know the position and orientation of sensors in relation to the head. With the adaptable on-scalp MEG sensor arrays, this coregistration becomes more challenging than in current SQUID-based MEG systems that use rigid sensor arrays. Here, we used simulations to quantify how accurately one needs to know the position and orientation of sensors in an on-scalp MEG system. The effects that different types of localisation errors have on forward modelling and source estimates obtained by minimum-norm estimation, dipole fitting, and beamforming are detailed. We found that sensor position errors generally have a larger effect than orientation errors and that these errors affect the localisation accuracy of superficial sources the most. To obtain similar or higher accuracy than with current SQUID-based MEG systems, RMS sensor position and orientation errors should be [Formula: see text] and [Formula: see text], respectively.

Neural Substrate for Metacognitive Accuracy of Tactile Working Memory.

The human prefrontal cortex (PFC) has been shown to be important for metacognition, the capacity to monitor and control one's own cognitive processes. Here we dissected the neural architecture of somatosensory metacognition using navigated single-pulse transcranial magnetic stimulation (TMS) to modulate tactile working memory (WM) processing. We asked subjects to perform tactile WM tasks and to give a confidence rating for their performance after each trial. We circumvented the challenge of interindividual variability in functional brain anatomy by applying TMS to two PFC areas that, according to tractography, were neurally connected with the primary somatosensory cortex (S1): one area in the superior frontal gyrus (SFG), another in the middle frontal gyrus (MFG). These two PFC locations and a control cortical area were stimulated during both spatial and temporal tactile WM tasks. We found that tractography-guided TMS of the SFG area selectively enhanced metacognitive accuracy of tactile temporal, but not spatial WM. Stimulation of the MFG area that was also neurally connected with the S1 had no such effect on metacognitive accuracy of either the temporal or spatial tactile WM. Our findings provide causal evidence that the PFC contains distinct neuroanatomical substrates for introspective accuracy of tactile WM.

Optical Co-registration of MRI and On-scalp MEG.

To estimate the neural generators of magnetoencephalographic (MEG) signals, MEG data have to be co-registered with an anatomical image, typically an MR image. Optically-pumped magnetometers (OPMs) enable the construction of on-scalp MEG systems providing higher sensitivity and spatial resolution than conventional SQUID-based MEG systems. We present a co-registration method that can be applied to on-scalp MEG systems, regardless of the number of sensors. We apply a structured-light scanner to create a surface mesh of the subject's head and the sensor array, which we fit to the MR image. We quantified the reproducibility of the mesh and localised current dipoles with a phantom. Additionally, we measured somatosensory evoked fields (SEFs) to median nerve stimulation and compared the dipole positions between on-scalp and SQUID-based systems. The scanner reproduced the head surface with <1 mm error. Phantom dipoles were localised with 2.1 mm mean error. SEF dipoles corresponding to the P35m response for OPMs were well localised to the somatosensory cortex, while SQUID dipoles for two subjects were erroneously localised to the motor cortex. The developed co-registration method is inexpensive, fast and can easily be applied to on-scalp MEG. It is more convenient than traditional co-registration methods while also being more accurate.