Optimal design of on-scalp electromagnetic sensor arrays for brain source localisation

A fully integrated whole-head helium OPM MEG: a performance assessment compared to cryogenic MEG

Magnetoencephalography (MEG) is a neuroimaging technique that measures neuronal activity at a millisecond scale. A few years ago, a new generation of MEG sensors emerged: optically pumped magnetometers (OPMs). The most common OPMs use alkali atoms as the sensing element. These alkali OPM sensors must be heated to approximately 150°C, in contrast to classical MEG sensors [superconducting quantum interference device MEG], which need to be cooled down to -269°C. This article focuses on a new kind of OPM that uses Helium-4 gas as the sensing element, which solves some disadvantages of alkali OPMs. 4He-OPM sensors operate at room temperature, with negligible heat dissipation (10 mW) and thus do not need thermal insulation. They also offer a large dynamic range (±200 nT) and frequency bandwidth (2,000 Hz). The main goal of this study is to characterize the performance of a whole-head MEG system based on 4He OPM sensors (4He OPM MEG). We first simulated different sensor configurations with three different numbers of channels and three different head sizes, from child to adult, in order to assess the signal-to-noise ratio and the source reconstruction accuracy. Experimental testing was also performed using a phantom to simulate brain magnetic activity. The simulation and experiments show equivalent detection capability and localization accuracy on both MEG systems. These results illustrate the benefit of 4He OPM sensors that operate at room temperature and are positioned closer to the scalp.

Magnetically shielded high-resolution visual stimulation for OPM-MEG applications

Many magnetoencephalography (MEG) experiments require visual stimulation (VS) inside a magnetically shielded room (MSR). For conventional MEG utilizing superconducting quantum interference devices (SQUIDs), the participant's head must stay within the semi-spherical surface of a cryogenic storage Dewar. This design allows to have many SQUID sensors as close as possible to the head in order to achieve good signal quality. Because Dewars have very restricted mobility, VS is usually realized using a projector outside of the MSR, some optical elements and a back-projection screen in the line of sight of the participant.Recently, the feasibility of MEG using optically pumped magnetometers (OPMs) was demonstrated. These sensors can be attached directly to the head because they operate near room temperature. OPM-MEG therefore offers more experimental freedom including different postures, movements or hyperscanning, creating the need for a more flexible kind of VS setup.In this paper, we present a compact, high-resolution VS setup which is enclosed by a portable magnetic shield with an opening for the projection. The VS setup is based on a single-board computer which acts as experiment control device to create visual stimuli, process inputs, log participant activity and set off trigger signals. This setup supports the new possibilities of OPM-MEG and can be easily installed into any MSR. We investigate if the shielded VS inside the MSR generates distortion signals above the noise floor of the OPMs. We also show that visual cortex activity can be evoked with our setup and recorded with a custom-made OPM-MEG cap. By applying two well-established visual stimulation paradigms, we demonstrate the ability of our setup to elicit brain activity in different frequency ranges.

Optimizing magnetometers arrays and analysis pipelines for multivariate pattern analysis

BACKGROUND

Multivariate pattern analysis (MVPA) has proven an excellent tool in cognitive neuroscience. It also holds a strong promise when applied to optically-pumped magnetometer-based magnetoencephalography.

NEW METHOD

To optimize OPM-MEG systems for MVPA experiments this study examines data from a conventional MEG magnetometer array, focusing on appropriate noise reduction techniques for magnetometers. We determined the least required number of sensors needed for robust MVPA for image categorization experiments.

RESULTS

We found that the use of signal space separation (SSS) without a proper regularization significantly lowered the classification accuracy considering a sub-array of 102 magnetometers or a sub-array of 204 gradiometers. We also found that classification accuracy did not improve when going beyond 30 sensors irrespective of whether SSS has been applied.

COMPARISON WITH EXISTING METHODS

The power spectra of data filtered with SSS has a substantially higher noise floor that data cleaned with SSP or HFC. Consequently, MVPA decoding results obtained from the SSS-filtered data are significantly lower compared to all other methods employed.

CONCLUSIONS

When designing MEG system based on SQUID magnetometers optimized for multivariate analysis for image categorization experiments, about 30 magnetometers are sufficient. We advise against applying SSS filters without a proper regularization to data from MEG and OPM systems prior to performing MVPA as this method, albeit reducing low-frequency external noise contributions, also introduces an increase in broadband noise. We recommend employing noise reduction techniques that either decrease or maintain the noise floor of the data like signal-space projection, homogeneous field correction and gradient noise reduction.

Applications of OPM-MEG for translational neuroscience: a perspective

Magnetoencephalography (MEG) allows the non-invasive measurement of brain activity at millisecond precision combined with localization of the underlying generators. So far, MEG-systems consisted of superconducting quantum interference devices (SQUIDS), which suffer from several limitations. Recent technological advances, however, have enabled the development of novel MEG-systems based on optically pumped magnetometers (OPMs), offering several advantages over conventional SQUID-MEG systems. Considering potential improvements in the measurement of neuronal signals as well as reduced operating costs, the application of OPM-MEG systems for clinical neuroscience and diagnostic settings is highly promising. Here we provide an overview of the current state-of-the art of OPM-MEG and its unique potential for translational neuroscience. First, we discuss the technological features of OPMs and benchmark OPM-MEG against SQUID-MEG and electroencephalography (EEG), followed by a summary of pioneering studies of OPMs in healthy populations. Key applications of OPM-MEG for the investigation of psychiatric and neurological conditions are then reviewed. Specifically, we suggest novel applications of OPM-MEG for the identification of biomarkers and circuit deficits in schizophrenia, dementias, movement disorders, epilepsy, and neurodevelopmental syndromes (autism spectrum disorder and attention deficit hyperactivity disorder). Finally, we give an outlook of OPM-MEG for translational neuroscience with a focus on remaining methodological and technical challenges.

Single-trial classification of evoked responses to auditory tones using OPM- and SQUID-MEG

Objective.Optically pumped magnetometers (OPMs) are emerging as a near-room-temperature alternative to superconducting quantum interference devices (SQUIDs) for magnetoencephalography (MEG). In contrast to SQUIDs, OPMs can be placed in a close proximity to subject's scalp potentially increasing the signal-to-noise ratio and spatial resolution of MEG. However, experimental demonstrations of these suggested benefits are still scarce. Here, to compare a 24-channel OPM-MEG system to a commercial whole-head SQUID system in a data-driven way, we quantified their performance in classifying single-trial evoked responses.Approach.We measured evoked responses to three auditory tones in six participants using both OPM- and SQUID-MEG systems. We performed pairwise temporal classification of the single-trial responses with linear discriminant analysis as well as multiclass classification with both EEGNet convolutional neural network and xDAWN decoding.Main results.OPMs provided higher classification accuracies than SQUIDs having a similar coverage of the left hemisphere of the participant. However, the SQUID sensors covering the whole helmet had classification scores larger than those of OPMs for two of the tone pairs, demonstrating the benefits of a whole-head measurement.Significance.The results demonstrate that the current OPM-MEG system provides high-quality data about the brain with room for improvement for high bandwidth non-invasive brain-computer interfacing.

A minimum assumption approach to MEG sensor array design

Objective.Our objective is to formulate the problem of the magnetoencephalographic (MEG) sensor array design as a well-posed engineering problem of accurately measuring the neuronal magnetic fields. This is in contrast to the traditional approach that formulates the sensor array design problem in terms of neurobiological interpretability the sensor array measurements.Approach.We use the vector spherical harmonics (VSH) formalism to define a figure-of-merit for an MEG sensor array. We start with an observation that, under certain reasonable assumptions, any array ofmperfectly noiseless sensors will attain exactly the same performance, regardless of the sensors' locations and orientations (with the exception of a negligible set of singularly bad sensor configurations). We proceed to the conclusion that under the aforementioned assumptions, the only difference between different array configurations is the effect of (sensor) noise on their performance. We then propose a figure-of-merit that quantifies, with a single number, how much the sensor array in question amplifies the sensor noise.Main results.We derive a formula for intuitively meaningful, yet mathematically rigorous figure-of-merit that summarizes how desirable a particular sensor array design is. We demonstrate that this figure-of-merit is well-behaved enough to be used as a cost function for a general-purpose nonlinear optimization methods such as simulated annealing. We also show that sensor array configurations obtained by such optimizations exhibit properties that are typically expected of 'high-quality' MEG sensor arrays, e.g. high channel information capacity.Significance.Our work paves the way toward designing better MEG sensor arrays by isolating the engineering problem of measuring the neuromagnetic fields out of the bigger problem of studying brain function through neuromagnetic measurements.

Sensor array design of optically pumped magnetometers for accurately estimating source currents

An optically pumped magnetometer (OPM) is a new generation of magnetoencephalography (MEG) devices that is small, light, and works at room temperature. Due to these characteristics, OPMs enable flexible and wearable MEG systems. On the other hand, if we have a limited number of OPM sensors, we need to carefully design their sensor arrays depending on our purposes and regions of interests (ROIs). In this study, we propose a method that designs OPM sensor arrays for accurately estimating the cortical currents at the ROIs. Based on the resolution matrix of minimum norm estimate (MNE), our method sequentially determines the position of each sensor to optimize its inverse filter pointing to the ROIs and suppressing the signal leakage from the other areas. We call this method the Sensor array Optimization based on Resolution Matrix (SORM). We conducted simple and realistic simulation tests to evaluate its characteristics and efficacy for real OPM-MEG data. SORM designed the sensor arrays so that their leadfield matrices had high effective ranks as well as high sensitivities to ROIs. Although SORM is based on MNE, the sensor arrays designed by SORM were effective not only when we estimated the cortical currents by MNE but also when we did so by other methods. With real OPM-MEG data we confirmed its validity for real data. These analyses suggest that SORM is especially useful when we want to accurately estimate ROIs' activities with a limited number of OPM sensors, such as brain-machine interfaces and diagnosing brain diseases.

Yttrium-Iron Garnet Magnetometer in MEG: Advance towards Multi-Channel Arrays

Recently, a new kind of sensor applicable in magnetoencephalography (MEG) has been presented: a solid-state yttrium-iron garnet magnetometer (YIGM). The feasibility of yttrium-iron garnet magnetometers (YIGMs) was demonstrated in an alpha-rhythm registration experiment. In this paper, we propose the analysis of lead-field matrices for different possible multi-channel on-scalp sensor layouts using YIGMs with respect to information theory. Real noise levels of the new sensor were used to compute signal-to-noise ratio (SNR) and total information capacity (TiC), and compared with corresponding metrics that can be obtained with well-established MEG systems based on superconducting quantum interference devices (SQUIDs) and optically pumped magnetometers (OPMs). The results showed that due to YIGMs' proximity to the subject's scalp, they outperform SQUIDs and OPMs at their respective noise levels in terms of SNR and TiC. However, the current noise levels of YIGM sensors are unfortunately insufficient for constructing a multichannel YIG-MEG system. This simulation study provides insight into the direction for further development of YIGM sensors to create a multi-channel MEG system, namely, by decreasing the noise levels of sensors.

Exploring the limits of MEG spatial resolution with multipolar expansions

The advent of scalp magnetoencephalography (MEG) based on optically pumped magnetometers (OPMs) may represent a step change in the field of human electrophysiology. Compared to cryogenic MEG based on superconducting quantum interference devices (SQUIDs, placed 2-4 cm above scalp), scalp MEG promises significantly higher spatial resolution imaging but it also comes with numerous challenges regarding how to optimally design OPM arrays. In this context, we sought to provide a systematic description of MEG spatial resolution as a function of the number of sensors (allowing comparison of low- vs. high-density MEG), sensor-to-brain distance (cryogenic SQUIDs vs. scalp OPM), sensor type (magnetometers vs. gradiometers; single- vs. multi-component sensors), and signal-to-noise ratio. To that aim, we present an analytical theory based on MEG multipolar expansions that enables, once supplemented with experimental input and simulations, quantitative assessment of the limits of MEG spatial resolution in terms of two qualitatively distinct regimes. In the regime of asymptotically high-density MEG, we provide a mathematically rigorous description of how magnetic field smoothness constraints spatial resolution to a slow, logarithmic divergence. In the opposite regime of low-density MEG, it is sensor density that constraints spatial resolution to a faster increase following a square-root law. The transition between these two regimes controls how MEG spatial resolution saturates as sensors approach sources of neural activity. This two-regime model of MEG spatial resolution integrates known observations (e.g., the difficulty of improving spatial resolution by increasing sensor density, the gain brought by moving sensors on scalp, or the usefulness of multi-component sensors) and gathers them under a unifying theoretical framework that highlights the underlying physics and reveals properties inaccessible to simulations. We propose that this framework may find useful applications to benchmark the design of future OPM-based scalp MEG systems.

Non-invasive measurements of ictal and interictal epileptiform activity using optically pumped magnetometers

Magneto- and electroencephalography (MEG/EEG) are important techniques for the diagnosis and pre-surgical evaluation of epilepsy. Yet, in current cryogen-based MEG systems the sensors are offset from the scalp, which limits the signal-to-noise ratio (SNR) and thereby the sensitivity to activity from deep structures such as the hippocampus. This effect is amplified in children, for whom adult-sized fixed-helmet systems are typically too big. Moreover, ictal recordings with fixed-helmet systems are problematic because of limited movement tolerance and/or logistical considerations. Optically Pumped Magnetometers (OPMs) can be placed directly on the scalp, thereby improving SNR and enabling recordings during seizures. We aimed to demonstrate the performance of OPMs in a clinical population. Seven patients with challenging cases of epilepsy underwent MEG recordings using a 12-channel OPM-system and a 306-channel cryogen-based whole-head system: three adults with known deep or weak (low SNR) sources of interictal epileptiform discharges (IEDs), along with three children with focal epilepsy and one adult with frequent seizures. The consistency of the recorded IEDs across the two systems was assessed. In one patient the OPMs detected IEDs that were not found with the SQUID-system, and in two patients no IEDs were found with either system. For the other patients the OPM data were remarkably consistent with the data from the cryogenic system, noting that these were recorded in different sessions, with comparable SNRs and IED-yields overall. Importantly, the wearability of OPMs enabled the recording of seizure activity in a patient with hyperkinetic movements during the seizure. The observed ictal onset and semiology were in agreement with previous video- and stereo-EEG recordings. The relatively affordable technology, in combination with reduced running and maintenance costs, means that OPM-based MEG could be used more widely than current MEG systems, and may become an affordable alternative to scalp EEG, with the potential benefits of increased spatial accuracy, reduced sensitivity to volume conduction/field spread, and increased sensitivity to deep sources. Wearable MEG thus provides an unprecedented opportunity for epilepsy, and given its patient-friendliness, we envisage that it will not only be used for presurgical evaluation of epilepsy patients, but also for diagnosis after a first seizure.

Optimising the sensing volume of OPM sensors for MEG source reconstruction

Magnetoencephalography (MEG) based on optically pumped magnetometers (OPMs) has been hailed as the future of electrophysiological recordings from the human brain. In this work, we investigate how the dimensions of the sensing volume (the vapour cell) affect the performance of both a single OPM-MEG sensor and a multi-sensor OPM-MEG system. We consider a realistic noise model that accounts for background brain activity and residual noise. By using source reconstruction metrics such as localization accuracy and time-course reconstruction accuracy, we demonstrate that the best overall sensitivity and reconstruction accuracy are achieved with cells that are significantly longer and wider that those of the majority of current commercial OPM sensors. Our work provides useful tools to optimise the cell dimensions of OPM sensors in a wide range of environments.

Spherical harmonic based noise rejection and neuronal sampling with multi-axis OPMs

In this study we explore the interference rejection and spatial sampling properties of multi-axis Optically Pumped Magnetometer (OPM) data. We use both vector spherical harmonics and eigenspectra to quantify how well an array can separate neuronal signal from environmental interference while adequately sampling the entire cortex. We found that triaxial OPMs have superb noise rejection properties allowing for very high orders of interference (L=6) to be accounted for while minimally affecting the neural space (2dB attenuation for a 60-sensor triaxial system). We show that at least 11th order (143 spatial degrees of freedom) irregular solid harmonics or 95 eigenvectors of the lead field are needed to model the neural space for OPM data (regardless of number of axes measured). This can be adequately sampled with 75-100 equidistant triaxial sensors (225-300 channels) or 200 equidistant radial channels. In other words, ordering the same number of channels in triaxial (rather than purely radial) configuration may give significant advantages not only in terms of external noise rejection but also by minimizing cost, weight and cross-talk.

A Comparative Study on the Influence of Undersampling and Oversampling Techniques for the Classification of Physical Activities Using an Imbalanced Accelerometer Dataset

Accelerometer data collected from wearable devices have recently been used to monitor physical activities (PAs) in daily life. While the intensity of PAs can be distinguished with a cut-off approach, it is important to discriminate different behaviors with similar accelerometry patterns to estimate energy expenditure. We aim to overcome the data imbalance problem that negatively affects machine learning-based PA classification by extracting well-defined features and applying undersampling and oversampling methods. We extracted various temporal, spectral, and nonlinear features from wrist-, hip-, and ankle-worn accelerometer data. Then, the influences of undersampilng and oversampling were compared using various ML and DL approaches. Among various ML and DL models, ensemble methods including random forest (RF) and adaptive boosting (AdaBoost) exhibited great performance in differentiating sedentary behavior (driving) and three walking types (walking on level ground, ascending stairs, and descending stairs) even in a cross-subject paradigm. The undersampling approach, which has a low computational cost, exhibited classification results unbiased to the majority class. In addition, we found that RF could automatically select relevant features for PA classification depending on the sensor location by examining the importance of each node in multiple decision trees (DTs). This study proposes that ensemble learning using well-defined feature sets combined with the undersampling approach is robust for imbalanced datasets in PA classification. This approach will be useful for PA classification in the free-living situation, where data imbalance problems between classes are common.

On-scalp magnetocorticography with optically pumped magnetometers: Simulated performance in resolving simultaneous sources

Currently, the gold standard for high-resolution mapping of cortical electrophysiological activity is invasive electrocorticography (ECoG), a procedure that carries with it the risk of serious morbidity and mortality. Due to these risks, the use of ECoG is largely limited to pre-surgical mapping in intractable epilepsy. Nevertheless, many seminal studies in neuroscience have utilized ECoG to explore domains such as visual perception, attention, auditory processing, and sensorimotor behavior. Studies such as these, occurring in patients with epilepsy rather than healthy controls, may lack generalizability, and are limited by the placement of the electrode arrays over the presumed seizure focus. This manuscript explores the use of optically pumped magnetometers (OPMs) to create a non-invasive alternative to ECoG, which we refer to as magnetocorticography. Because prior ECoG studies reveal that most cognitive processes are driven by multiple, simultaneous independent neuronal assemblies, we characterize the ability of a theoretical 56-channel dense OPM array to resolve simultaneous independent sources, and compare it to currently available SQUID devices, as well as OPM arrays with inter-sensor spacings more typical of other systems in development. Our evaluation of this theoretical system assesses many potential sources of error, including errors of sensor calibration and position. In addition, we investigate the influence of geometrical and anatomical factors on array performance. Our simulations reveal the potential of high-density, on-scalp OPM MEG devices to localize electrophysiological brain responses at unprecedented resolution for a non-invasive device.

Simulation Study of Different OPM-MEG Measurement Components

Magnetoencephalography (MEG) is a neuroimaging technique that measures the magnetic fields of the brain outside of the head. In the past, the most suitable magnetometer for MEG was the superconducting quantum interference device (SQUID), but in recent years, a new type has also been used, the optically pumped magnetometer (OPM). OPMs can be configured to measure multiple directions of magnetic field simultaneously. This work explored whether combining multiple directions of the magnetic field lowers the source localization error of brain sources under various conditions of noise. We simulated dipolar-like sources for multiple configurations of both SQUID- and OPM-MEG systems. To test the performance of a given layout, we calculated the average signal-to-noise ratio and the root mean square of the simulated magnetic field; furthermore, we evaluated the performance of the dipole fit. The results showed that the field direction normal to the scalp yields a higher signal-to-noise ratio and that ambient noise has a much lower impact on its localization error; therefore, this is the optimal choice for source localization when only one direction of magnetic field can be measured. For a low number of OPMs, combining multiple field directions greatly improves the source localization results. Lastly, we showed that MEG sensors that can be placed closer to the brain are more suitable for localizing deeper sources.

Transforming and comparing data between standard SQUID and OPM-MEG systems

Optically pumped magnetometers (OPMs) have recently become so sensitive that they are suitable for use in magnetoencephalography (MEG). These sensors solve operational problems of the current standard MEG, where superconducting quantum interference device (SQUID) gradiometers and magnetometers are being used. The main advantage of OPMs is that they do not require cryogenics for cooling. Therefore, they can be placed closer to the scalp and are much easier to use. Here, we measured auditory evoked fields (AEFs) with both SQUID- and OPM-based MEG systems for a group of subjects to better understand the usage of a limited sensor count OPM-MEG. We present a theoretical framework that transforms the within subject data and equivalent simulation data from one MEG system to the other. This approach works on the principle of solving the inverse problem with one system, and then using the forward model to calculate the magnetic fields expected for the other system. For the source reconstruction, we used a minimum norm estimate (MNE) of the current distribution. Two different volume conductor models were compared: the homogeneous conducting sphere and the three-shell model of the head. The transformation results are characterized by a relative error and cross-correlation between the measured and the estimated magnetic field maps of the AEFs. The results for both models are encouraging. Since some commercial OPMs measure multiple components of the magnetic field simultaneously, we additionally analyzed the effect of tangential field components. Overall, our dual-axis OPM-MEG with 15 sensors yields similar information to a 62-channel SQUID-MEG with its field of view restricted to the right hemisphere.

Imaging somatosensory cortex responses measured by OPM-MEG: Variational free energy-based spatial smoothing estimation approach

In recent years, optically pumped magnetometer (OPM)-based magnetoencephalography (MEG) has shown potential for analyzing brain activity. It has a flexible sensor configuration and comparable sensitivity to conventional SQUID-MEG. We constructed a 32-channel OPM-MEG system and used it to measure cortical responses to median and ulnar nerve stimulations. Traditional magnetic source imaging methods tend to blur the spatial extent of sources. Accurate estimation of the spatial size of the source is important for studying the organization of brain somatotopy and for pre-surgical functional mapping. We proposed a new method called variational free energy-based spatial smoothing estimation (FESSE) to enhance the accuracy of mapping somatosensory cortex responses. A series of computer simulations based on the OPM-MEG showed better performance than the three types of competing methods under different levels of signal-to-noise ratios, source patch sizes, and co-registration errors. FESSE was then applied to the source imaging of the OPM-MEG experimental data.

A Proposed Brain-, Spine-, and Mental- Health Screening Methodology (NEUROSCREEN) for Healthcare Systems: Position of the Society for Brain Mapping and Therapeutics

The COVID-19 pandemic has accelerated neurological, mental health disorders, and neurocognitive issues. However, there is a lack of inexpensive and efficient brain evaluation and screening systems. As a result, a considerable fraction of patients with neurocognitive or psychobehavioral predicaments either do not get timely diagnosed or fail to receive personalized treatment plans. This is especially true in the elderly populations, wherein only 16% of seniors say they receive regular cognitive evaluations. Therefore, there is a great need for development of an optimized clinical brain screening workflow methodology like what is already in existence for prostate and breast exams. Such a methodology should be designed to facilitate objective early detection and cost-effective treatment of such disorders. In this paper we have reviewed the existing clinical protocols, recent technological advances and suggested reliable clinical workflows for brain screening. Such protocols range from questionnaires and smartphone apps to multi-modality brain mapping and advanced imaging where applicable. To that end, the Society for Brain Mapping and Therapeutics (SBMT) proposes the Brain, Spine and Mental Health Screening (NEUROSCREEN) as a multi-faceted approach. Beside other assessment tools, NEUROSCREEN employs smartphone guided cognitive assessments and quantitative electroencephalography (qEEG) as well as potential genetic testing for cognitive decline risk as inexpensive and effective screening tools to facilitate objective diagnosis, monitor disease progression, and guide personalized treatment interventions. Operationalizing NEUROSCREEN is expected to result in reduced healthcare costs and improving quality of life at national and later, global scales.