# Helium optically pumped magnetometers can detect epileptic abnormalities as well as squids as shown by intracerebral recordings

Jean-Michel Badier, Denis Schwartz, Christian-George Bénar, Khoubeib

Kanzari, Sébastien Daligault, Rudy Romain, Sergey Mitryukovskiy, William

Fourcault, Vincent Josselin, Matthieu Le Prado, et al.

► To cite this version:

Jean-Michel Badier, Denis Schwartz, Christian-George Bénar, Khoubeib Kanzari, Sébastien Daligault, et al. Helium optically pumped magnetometers can detect epileptic abnormalities as well as squids as shown by intracerebral recordings. eNeuro, 2023, 10 (12), pp.ENEURO.0222 - 23.2023. 10.1523/eneuro.0222-23.2023. cea-04756789

HAL Id: cea-04756789

https://cea.hal.science/cea-04756789v1

Submitted on 28 Oct 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Helium Optically Pumped Magnetometers Can Detect Epileptic Abnormalities as Well as SQUIDs as Shown by Intracerebral Recordings

Jean-Michel Badier,¹ Denis Schwartz,² Christian-George Bénar,¹ Khoubeib Kanzari,¹

Sébastien Daligault,² Rudy Romain,³,⁴ Sergey Mitryukovskiy,³,⁴ William Fourcault,³ Vincent Josselin,³

Matthieu Le Prado,³,⁴ Julien Jung,⁵ Augustin Palacios-Laloy,³,⁴ Carron Romain,¹,⁶ Fabrice Bartolomei,¹,⁷

Etienne Labyt,³,⁴ and Francesca Bonini¹,⁷

https://doi.org/10.1523/ENEURO.0222-23.2023

¹Institut de Neurosciences des Systèmes, Institut National de la Santé et de la Recherche Médicale, Aix Marseille Université, Marseille 13005, France, ²MEG Département, CERMEP-Imagerie du Vivant, Lyon 69003, France, ³CEA-LETI, MINATEC, Université Grenoble Alpes, Grenoble 38054, France, ⁴MAG4Health, Grenoble 38000, France, ⁵Centre de Recherche en Neurosciences de Lyon, Unité Mixte de Recherche S1028, Centre National de la Recherche Scientifique, Hospices Civils de Lyon, Institut National de la Santé et de la Recherche Médicale, Université Lyon 1, Lyon 69002, France, ⁶Department of Functional and Stereotactic Neurosurgery, Hôpital de la Timone, Assistance Publique-Hôpitaux de Marseille, Marseille 3005, France, and ⁷Department of Epileptology and Cerebral Rythmology, Hôpital de la Timone, Assistance Publique-Hôpitaux de Marseille, Marseille 3005, France

## Abstract

Magnetoencephalography based on superconducting quantum interference devices (SQUIDs) has been shown to improve the diagnosis and surgical treatment decision for presurgical evaluation of drug-resistant epilepsy. Still, its use remains limited because of several constraints such as cost, fixed helmet size, and the obligation of immobility. A new generation of sensors, optically pumped magnetometers (OPMs), could overcome these limitations. In this study, we validate the ability of helium-based OPM ($^4$He-OPM) sensors to record epileptic brain activity thanks to simultaneous recordings with intracerebral EEG [stereotactic EEG (SEEG)]. We recorded simultaneous SQUIDs-SEEG and 4He-OPM-SEEG signals in one patient during two sessions. We show that epileptic activities on intracerebral EEG can be recorded by OPMs with a better signal-to noise ratio than classical SQUIDs. The OPM sensors open new venues for the widespread application of magnetoencephalography in the management of epilepsy and other neurologic diseases and fundamental neuroscience.

**Key words:** epilepsy; interictal epileptic discharges; magnetoencephalography; optically pumped magnetometers; simultaneous recording; stereotactic-EEG

## Significance Statement

We performed a simultaneous recording of $^4$He-OPM and intracerebral EEG and validate for the first time OPM results with signals recorded directly within the brain. We demonstrate that epileptic abnormalities seen on intracerebral electrodes are detected by OPMs with a better signal-to noise ratio than classical magnetoencephalography. This represents a significant step toward the validation of OPM-based recordings for epilepsy diagnosis and for clinical and neuroscience research.

Received June 27, 2023; accepted September 26, 2023; First published November 6, 2023.

M.L.-P., A.P.-L., and E.L. hold founding equity in Mag4Health, which is developing and commercializing MEG systems based on He-OPM technology. R.R. and S.M. are employees of Mag4Health. Mag4Health provided technical support for data acquisition. At the time the recordings were made, R.R., S.S., M.L.-P., A.P.-L., and E.L. were employees of CEA LETI.

Author contributions: J.-M.B., E.L., and F.Bonini designed research; J.-M.B., K.K., R.R., S.M., M.L., A.P.-L., C.R., and E.L. performed research; S.D., W.F., V.J., and J.J. contributed unpublished reagents/analytic tools; J.-M.B. and F.Bonini analyzed data; J.-M.B., D.S., C.B., F.Bartolomei, E.L., and F.Bonini wrote the paper.
## Introduction

Magnetoencephalography (MEG) is a noninvasive electro-
physiological recording technology that source reconstruc-
tion methods are applied to to localize brain signals with
millisecond resolution. The analysis of magnetic brain activity
allows for a deeper understanding of the neural substrates of
brain pathologies. In the management of drug-resistant epi-
lepsy, MEG contributes to the success of surgical treatment
through noninvasive localization of interictal epileptic dis-
charges (IEDs) and can also guide the implantation of depth
electrodes when invasive recordings are required (Fischer
et al., 2005; Murakami et al., 2016; Alkawadri et al., 2018;
Rampp et al., 2019; Corona et al., 2023).

Still, the spread of MEG is strongly limited by the con-
straints imposed by the use of superconducting quantum
interference devices (SQUIDs; Cohen, 1972; Hämäläinen
et al., 1993), requiring cooling at a very low temperature
(4.2 K). Thus, the sensors are isolated and enclosed in a
fixed array inside a rigid dewar, at least 3 cm from the
brain. The critical consequences are a substantial reduc-
tion in the magnetic signal amplitude, inhomogeneous
coverage in terms of brain source-sensor distance, and
the need for total immobility during signal acquisition as
well as between recording runs.

The emergence of a new generation of MEG sensors,
the optically pumped magnetometers (OPMs; Budker and
Romalis, 2007), could overcome these limitations. OPMs
are quantum sensors that exploit the interaction between
an atomic gas and laser light to obtain very precise meas-
urements of tiny magnetic fields. OPMs do not require ex-
treme cooling and can be placed near the scalp in a
wearable system thus allowing the subject's movement
(Johnson et al., 2013; Boto et al., 2017; livanainen et al.,
2017, Boto et al., 2018, 2022; Seymour et al., 2021). The
development and improvement of this new technology
has taken a quantum leap in recent years, with evolution
in miniaturization and sensitivity (for review, see Tierney et
al., 2019; Brookes et al., 2022; Labyt et al., 2022). The first
commercially available OPMs, based on Alkali atoms
(Budker and Romalis, 2007), have a sensitivity of 20 fem-
totesla/rootHz (fT/rtHz) in triaxis mode compared with 5
fT/rtHz for SQUID sensors. However, their placement
close to the scalp allows a threefold to eightfold increase
in signal power (Budker and Romalis, 2007; livanainen et
al., 2017; Johnson et al., 2013), with a significant neuro-
magnetic signal enhancement compared with SQUIDs,
namely for superficial sources (Boto et al., 2016, 2017;
livanainen et al., 2017).

This study was supported by Agence Nationale de la Recherche Grants ANR-22-CE19-0012-01, Grant ANR-11-INBS-0006 to France Life Imaging network, and ANR-16-CONV-0002 to Institute of Language, Communication, and the Brain.

Acknowledgment: We thank the nurses involved in the preparation and
recording.

Correspondence should be addressed to Francesca Bonini at francesca.bonini@univ-amu.fr.

https://doi.org/10.1523/ENEURO.0222-23.2023

Copyright © 2023 Badier et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

New OPMs based on helium atoms (⁴He-OPMs) have
been developed (Beato et al., 2018), which have a large dy-
namic range (up to 250 nanotesla, nT) and a large frequency
bandwidth (up to 2 kHz) with negligible heat dissipation
(10 mW per sensor). Their properties allow recording in
a standard magnetic shielding room, facilitating data
acquisition.

A few studies have shown that alkali OPMs can record IEDs comparable with those observed with SQUID-MEG (Feys et al., 2022) or previously obtained EEG (Vivekananda et al., 2020). However, to compare OPM to SQUID-MEG, it is necessary to ensure that both modalities are recording equivalent IEDs. By using intracerebral stereotactic EEG (SEEG) as a ground truth, simultaneously recorded with OPMs and SQUIDs (Badier et al., 2017), we aimed to demonstrate that ⁴He-OPMs can perform as well as SQUID-MEG but at an expected lower cost and with greater convenience. These unique simultaneous recordings allowed selecting similar IEDs to compare SQUID-MEG and ⁴He-OPM and thus evaluate the signal-to-noise ratio (SNR) of both techniques.

## Materials and Methods

The simultaneous recordings were performed inside a
two-layer Mu-metal magnetic shielded room (MSR) on
a patient with drug-resistant focal epilepsy undergoing
SEEG. Fourteen intracerebral electrodes with a total of
119 recording contacts were implanted mainly in the left
temporal structures (Fig. 1d). These electrodes, and par-
ticularly those exploring the anterior medial and lateral
temporal regions, disclosed abundant interictal epileptic
abnormalities. A first SQUID-MEG/SEEG recording ses-
sion (Fig. 1a) of 20 min at rest was followed by a compar-
able ⁴He-OPM/SEEG session. The four ⁴He-OPM sensors
(2 × 2 × 5 cm) fixed on a rigid helmet (Fig. 1b) were placed
over the left central and temporal regions, in contact with
the bandage covering the scalp (Fig. 1c, red dots). The
four SQUID sensors closest to the four ⁴He-OPMs sen-
sors were selected for analysis and for comparison with
the magnetic activity recorded by the OPM sensors (Fig.
1c, green dots). These four SQUIDs were between 2.72
and 3.23 cm away from the OPM sensors.

## Patient

The patient was a 33-year-old woman with intractable
temporal epilepsy associated with a cavernoma of the left
parahippocampal gyrus. An SEEG was performed to de-
fine the epileptogenic zone to be removed to study its re-
lationship with the functional areas of language and to
evaluate the functionality of the hippocampus and the
possible functional risk in terms of memory in case of epi-
lepsy surgery. Most electrodes (12/14) were aimed at an
extensive exploration of the left temporal structures, as
well as the anterior insula and the orbitofrontal cortex; two
electrodes were implanted in the right anterior temporal
region. At the end of the SEEG, the epileptogenic zone
network could be defined as involving the left mesial tem-
poral structures including the amygdala, hippocampus,
rhinal cortex, left parahippocampal cortex and collateral
sulcus (the latter being posterior to the cavernoma), and the
left temporal pole. The simultaneous recording session was
<figure><img src="image_2.png" /><figcaption>**Figure 1.** SQUID-MEG, ⁴He-OPM-MEG, and SEEG recording setup. **a**, Simultaneous SQUID-MEG/SEEG. The classic cryogenic MEG system, measuring 120 × 100 cm and weighing ~300 kg with 265 SQUID sensors in a fixed array requiring subject immobility during data acquisition. **b**, Simultaneous ⁴He-OPM-MEG/SEEG: recording configuration composed of four sensors integrated into a wearable helmet placed on the scalp and in contact with the bandage covering the SEEG electrode inserts (the cables and connectors outside the helmet are visible). Insert, A photo of the ⁴He-OPM sensor. **c**, Three-dimensional reconstruction of the patient's head from MRI, with ⁴He-OPM (red, named Cant, Cpost, Tant, Tpost), SQUID sensors (blue), and SEEG electrodes entry points (orange). The four SQUID sensors (named A95, A68, A179 and A157) closest to the OPMs are in green. Note the distance between the SQUIDs and the scalp (at least 3 cm), whereas the ⁴He-OPMs are in contact with it. **d**, SEEG implantation consists of 14 intracerebral electrodes, 2 in the right hemisphere (not shown) and 12 in the left hemisphere exploring the whole left temporal structures. A', Amygdala; TB', rhinal cortex; C', posterior hippocampus; GPH', parahippocampal gyrus; T', anterior insula/lateral T1; H', thalamus/Heschl gyrus; la', anterior insula/F2; FCA', lingual gyrus; GC', posterior cingulate/T1; OR', orbitofrontal cortex/middle frontal sulcus.</figcaption></figure>

conducted at the end of the long-term video SEEG record-
ing, at J11 from electrode implantation once all clinical data
were acquired. Ethical approval was obtained at the Comité
de Protection des Personnes Sud Méditerranée I under ID
RCB 2020-A01830-39, and the patient gave informed
consent.

SEEG recordings

SEEG exploration was performed using intracerebral
multiple-contacts electrodes placed intracranially with ro-
botic assistance (ROSA, Zimmer Biomet) and intraope-
rative Mobius Airo CT-scan (Stryker) verification. Small
insertion screws (catalog #2023-VG-C-10 or #2023-VG-
C-10-15, Alcis) were used to minimize the bulk of the elec-
trodes around the skull and to allow for simultaneous re-
cordings. The electrodes had a diameter of 0.8 mm and
contained 10–15 contacts. Each contact consisted of 2-

mm-long platinum iridium and was separated from each
other by 1.5 mm of insulated material (Alcis). To accu-
rately define the anatomic position of each SEEG contact
along the electrode trajectory, a CT-scan/MRI data fusion
was performed using the in-house software GARDEL (a
Graphical User Interface for Automatic Registration and
Depth Electrodes Localization). This MATLAB-based tool
is able to coregister the MRI to the CT scan and automati-
cally segment and localize depth electrodes contacts by
image processing. The signals were formatted in a bipolar
configuration keeping only noncontiguous bipolar channels.
Recording was made with a Brain Products BrainAmp DC
amplifier. SEEG data have been sampled at 2500 Hz.

SQUID-MEG recordings

The first recording was made with a SQUID MEG system
simultaneously with the SEEG (see above) while the patient
was in a supine position. MEG signals were acquired on a 4D Neuroimaging 3600 whole-head system at a sampling rate of 2034.51 Hz with a total of 248 magnetometers. Additionally, three magnetometers and nine gradiometers were used for noise compensation. Electrocardiographic and the electrooculographic activity was recorded on bipo- lar EEG channels.

⁴He-OPM recordings

The second recording was made with a prototype of five ⁴He-OPM sensors simultaneously with SEEG, also while the patient was in a supine position. The ⁴He-OPMs are sensors that measure the brain magnetic field along the three axes with continuous self-compensation of the magnetic field on all axes. The magnetic field measure- ment relies on a measure of the variation of the light ab- sorption caused by the deviation of the electronic spin of ⁴He atoms from the state originally set by laser pumping. Technical information and physical principles used in our sensors can be found in a previous publication (Fourcault et al., 2021) and are summarized below. The OPM used in this study is based on parametric resonance of helium-4 metastable atoms at near zero magnetic field (DuPont- Roc, 1971; Beato et al., 2018). The cell containing the ⁴He gas is a cylinder 1 cm in diameter and 1 cm in height. This cell, placed at the bottom of the sensor, is surrounded by small three-axis Helmholtz coils that are used to apply both the radiofrequency (RF) fields and the compensation fields (see below). A high-frequency (HF) discharge (ex- cited between 10 and 20 MHz and consuming ~10 mW power) excites the ⁴He atoms from their ground state to the metastable triplet state, which has three Zeeman sub- levels. Selective optical pumping (with a linearly polarized beam tuned on the D0 line at 1083 nm) was performed to prepare a macroscopic magnetic moment on the gas, which evolves in the magnetic field created by the brain. In our OPMs, to derive a vector measurement of the three components of the magnetic field, two RF fields are ap- plied to the gas, B<sub>Ω</sub>cosΩt along one tangential x-axis and B<sub>ω</sub>cosωt along the radial y-axis. Both are orthogonal to the polarization of the pump laser beam. Using this scheme, first introduced by DuPont-Roc (1971), three resonance signals were detected on the transmitted pump light at Ω, ω, and ω ± Ω. At first order, the amplitude of each reso- nance is respectively proportional to one of the three com- ponents of the magnetic field to be measured (Bx, By, and Bz, respectively).

Each sensor is operated in a closed-loop mode on the three axes. This consists in continuously cancelling the three components of the magnetic field by applying an op- posite field with the three-axis Helmholtz coils. The value of each magnetic field component is deduced from the cur- rent injected in the compensation coil. In this way, the sen- sor becomes self-calibrated, that is, its output can only be affected by variations of the transfer function between the current and magnetic field set by the coil geometry and not by other operating parameters (light intensity, HF power, etc.). This closed-loop mode suppresses the cross-axis ef- fects (Cohen-Tannoudji et al., 1970a,b; DuPont-Roc, 1970, 1971) by which the measurement of one axis becomes

dependent on the field along another axis. This phenomen- on has been recently referred as the Cross-Axis Projection Error (CAPE; Borna et al., 2022), yielding both phase errors and a tilt of the sensing axis. Previous studies with alkali OPM operated in an open-loop mode have characterized an axis tilt of 3.3°/nT at low frequencies (Borna et al., 2022), and offsets change as small as ±1.5 nT resulted in gain er- rors of ~4% (Boto et al., 2018). Another study, focusing on an active shielding system, reported that the calibration er- rors of the alkaline OPM were three to four times higher when the ambient field was not compensated (Iivanainen et al., 2019). Thereby, ⁴He OPM is the first sensor, to our knowledge, to provide a measurement of the magnetic field components in a closed-loop mode along the three axes, guaranteeing the reliability of the measurement and avoiding any CAPE. Another important advantage of closed- loop operation is the possibility of broadening the dynamic range well above the atomic line width. A dynamic range of ±250 nT is currently achieved for our ⁴He OPMs. The sensi- tivity of our magnetometer operating in the closed-loop triax- ial mode is better than 50 ft/Hz¹/² on two of the three axes (the radial and one tangential) with a bandwidth going from DC to 2 kHz. The sensitivity of the third axis was 200 ft/Hz¹/². This axis was not used in further data processing steps in the present study.

However, if the closed-loop mode avoids CAPE, it has some unwanted consequences because of the cross talks that unavoidably exist between the sensors within the OPM array. This problem can be solved by appropriate postprocessing as far as the cross talks are appropriately characterized (Fourcault et al., 2021, provides a detailed de- scription). The measured cross talk matrix for an array of four ⁴He OPM sensors with only 2 mm spacing, which corre- sponds to an extremely unfavorable situation compared with a real OPM MEG setup, revealed low cross talk (<10%) and showed a good agreement with the estimated matrix from the Biot-Savart calculations. Knowing this cross talk matrix, minor cross-talk-related errors are corrected in the measurement by adequate postprocessing. OPM data have been sampled at 11161 Hz and downsampled to 4 kHz.

OPM, SQUID, and SEEG spatial coregistration

Spatial coregistration was performed for both sessions using a 3D Polhemus digitizer, based on three fiducial markers (nasion and left and right preauricular points). The quality of the coregistration was checked using the digitization of the facial mask. Following the standard proce- dure for 4D neuroimaging, the position of nasion and both tragi were digitized to allow the construction of the patient frame. The positions of five coils located on the subject’s head were also digitized. Activations of these five coils be- fore and after recording allowed to determine the location of the MEG sensors within the patient frame. For the ⁴He OPM-MEG/SEEG session, the same three fiducial markers were digitized in the same reference frame. In this case, the sensors were linked to the subject’s head by the ⁴He OPM headset. This headset was positioned on the subject’s head before insertion of the ⁴He OPM sensors, and the corre- sponding slots were digitized giving the location of each ⁴He OPM sensor within the same patient reference frame.
Finally, both types of sensors were localized in the same reference frame, that is, the subject's head. It was then possible to select the SQUID sensors closest to each ⁴He OPM based on Euclidean distance. A preliminary SQUID MEG recording, performed several months before in the context of a noninvasive presurgical evaluation, allowed us to appreciate the typical topography of IEDs of this patient and therefore to choose the placement of OPMs. The four OPM sensors were placed according to this topography, with the constraint of the limited number of locations on the helmet. The signal recorded from OPMs therefore does not necessarily correspond to the maximum amplitude that can be recorded on this patient. The signal obtained from each OPM is compared with the one recorded from the closest SQUID sensor. This one, given the geometry of the array of sensors, lies on a straight line starting from a point near the center of the head and passing through the OPM sensor. This makes it the best choice for comparing the two signals.

## Data processing

For noise correction, a noise compensation was performed for SQUIDs thanks to recordings from reference sensors (magnetometers and gradiometers). Noise compensation is obtained by subtracting the contribution of the noise measured by the references for each sensor. As there are not enough sensors to perform this, there is no such compensation for the OPM. In that case noise cancellation is only performed by filtering (bandpass filter 2–70 Hz and a notch filter at 50 Hz). For downsampling and temporal registration, because the data were acquired with different devices, postprocessing was necessary to temporally register the simultaneous recordings with the same sampling rate and the same number of samples. To do so, randomly distributed triggers were sent to both kinds of MEG sensors (SQUID or OPM) and SEEG. Data were then processed by an in-house MATLAB function to match the sequence of triggers and to resample the data (from SEEG to an MEG time frame). The resampling procedure did not result in differences in delays larger than one sample. The end results for each session (SQUID-SEEG and OPM-SEEG) were two files that contained the same number of samples with synchronized SEEG and MEG data. For filtering, all data were bandpass filtered at 2–70 Hz; a notch filter at 50 Hz was added. For SNR computation, to investigate the amplitude of events, the SNR was computed as the following ratio: the max amplitude of the event divided by the SD of a baseline (2 s of signal before the event for the single event presented in Fig. 2 and 500 ms for the averages in Fig. 3 and Fig. 4).

For automatic extraction of events of interest, we had to compare data recorded in the two separate sessions. To do so, we used a procedure to identify identical events across recordings. As a first step, an expert neurologist (F.Bonini) manually selected four kinds of events of interest from the SEEG recording of the SQUID session. The selection was done using a montage of one bipolar SEEG derivation per explored brain area. These reference events were chosen to be representative of IEDs involving lateral temporal structures (spike I) and IEDs involving both medial and lateral temporal structures (spike II). We

used the function findsignal in MATLAB (version 202b) to find similar intracerebral events for both the SQUID session and the OPM session. Three to 10 events were kept for averaging. From this set of events, we chose a representative occurrence of spike II to illustrate a single event (Fig. 1).

## Results

As a main result, the analysis of the simultaneous recordings globally shows that the activities recorded by the intracerebral electrodes placed in the left temporal pole are detected on the scalp surface by both systems, conventional SQUID-MEG and OPM-MEG (Figs. 2, 3, 4). Notably ⁴He-OPM-MEG recordings showed an SNR better than that of SQUID-MEG.

Single interictal events (epileptic spikes), are visible on both OPM and SQUID recordings. Figure 2 presents single epileptic spikes with the corresponding SEEG traces. Using bandpass filtering (2–70 Hz) and a notch filter (50 Hz) with no further processing, the ⁴He-OPM-MEG can clearly record, with a high SNR, the IEDs identified with the SEEG electrodes. A single epileptic spike arising from the left temporal pole and the adjacent anterior third temporal gyrus (Fig. 2c, intracerebral traces) is distinctly identified on the scalp by the ⁴He-OPM-MEG (SNR = 7, peak-to-peak amplitude = 5 pT) (Fig. 2d). In comparison, the four closest SQUID-MEG sensors, using advanced denoising based on reference sensors, detect an equivalent intracerebral spike (Fig. 2a, visible on simultaneous SEEG) but with a lower SNR of five (Fig. 2b). This can be explained by the reduced distance from the neuronal sources of ⁴He-OPM-MEG, yielding increased signal power compared with SQUID-MEG (Labyt et al., 2019), even if this should be weighed against the fact that the SEEG has a slightly higher SNR of the same order of magnitude. Furthermore, it is interesting to note that different time courses and polarities are clearly identifiable, depending on the location of the channel and the orientation of the magnetic field (Fig. 2d). The ⁴He-OPM sensors are natively sensitive to two orthogonal orientations of the magnetic field so that two different signals corresponding to the tangential and radial components (Fig. 2d, red lines and blue lines, respectively) are output for each sensor.

These results observed on single spikes are even more evident when comparing the average of ⁴He-OPM-MEG and SQUID-MEG signals (Fig. 3). To compare a sufficient number of equivalent epileptic events recorded by OPMs and SQUID, we manually identified, using a selective SEEG electrodes montage, four different types of epileptic spikes recorded in the two separate simultaneous sessions (see above, Materials and Methods; Extended Data Figs. 1-1, 2-1). These IEDs have been subsequently automatically extracted from the entire time series and averaged by type. Two types of IEDs are illustrated, an average spike arising from the left temporal pole only (Fig. 3a,c, red line) and an average spike arising from the temporal pole, the third anterior temporal gyrus (minimally), and the anterior hippocampus (Fig. 4a,c, red, orange, and blue lines, respectively). In both cases, the simultaneous SQUID-MEG averaged data reveal a small deflection (Figs. 3b, 4b, respectively, −6 pT
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Individual spike results. **a**, **b**, The signal collected during the SQUID-MEG/SEEG simultaneous session. Bipolar SEEG data (**a**). The name of the SEEG electrodes and the recording contact label are shown. The labels range from 1 (deeper location) to 11 (more superficial location). The spike (SNR = 14.9) clearly involves both deep and more superficial structures (amygdala, anterior and posterior hippocampus, third anterior temporal gyrus, and temporal pole). Simultaneous SQUID-MEG data collected on the four sensors closest to the ⁴He-OPM channels (**b**); an interictal epileptic spike appears at ~82.5 s, with a peak-to-peak amplitude of 1.1 pT. **c**, **d**, The signal collected during the ⁴He-OPM-MEG/SEEG simultaneous session. Bipolar SEEG data (**c**). Note the similarity between two intracerebral spikes disclosing the same anatomic location and time course (SNR = 24.8). Simultaneous ⁴He-OPM-MEG data collected on four sensors (**d**). t, Tangential magnetic field (red lines); r, radial magnetic field (blue lines). A spike appears at 342.5 s with a peak-to-peak amplitude of 2,5 pT. The vertical scale is identical to that of the SQUID data in **a**, **b**.</figcaption></figure>

and 3.7 pT maximum value) in correspondence with the averaged intracerebral spikes (Figs. 3a, 4a). In contrast, the ⁴He-OPM-MEG averaged data disclose clear averaged spikes (Figs. 3d, 4d) occurring along with the two averaged intracerebral spikes (Figs. 3c, 4c). Both ⁴He-OPM-MEG spikes have a higher amplitude (spike I = -9.5 pT; spike II = -10 pT) and a higher SNR (spike I, SNR, *r* = 6.4; spike II, SNR, *r* = 16.7) than those detected in the SQUID-MEG data. Regarding averaged spike II, it is interesting to note that the time course of the OPM signal appears to be correlated with the decay of the temporopolar spike, as recorded by temporal pole (TP') electrodes 3-4, whereas the hippocampal

activity recorded by anterior hippocampus (B') electrodes 2-3 is not detected by MEG sensors, SQUIDs nor OPMs. On the other hand, variable time courses, polarities, and amplitudes between radial and tangential measurement axes can be observed on the OPM signals, particularly on spike II (Fig. 4d). The radial components of the two anterior sensors show a negative deflection, whereas for the corresponding tangential components, the deflection is positive and slightly delayed. This suggests that ⁴He-OPM-MEG provides more information than SQUID-MEG about the spatiotemporal organization of IEDs across the cerebral cortex, thanks to their tangential component.
<figure><img src="image_4.png" /></figure>

**Figure 3.** Averaged Type I spikes (n = 10). **a**, **b**, The averaged signal collected during the SQUID-MEG/SEEG simultaneous session. **c**, **d**, The averaged signal collected during the ⁴He-OPM-MEG/SEEG simultaneous session. The names of the SEEG electrodes and recording contacts are shown. Contacts range from 1 (deeper location) to 9 (most surface location). Bipolar averaged SEEG data for the SQUID session (**a**). The averaged spike involves only the bipolar TP'3-TP'4 recording. Simultaneous SQUID-MEG averaged data on the four sensors closest to the ⁴He-OPM sensors (**b**); a small interictal spike appears ~25 ms with a peak-to-peak amplitude of 8 pT. Bipolar averaged SEEG data for the OPM session (**c**). Note the similarity between the intracerebral spikes of the two sessions, disclosing the same anatomic location and time course. ⁴He-OPM-MEG averaged data collected on four sensors, T<sub>pos</sub>, C<sub>ant</sub>, C<sub>post</sub>, and T<sub>ant</sub> (**d**). pos, Posterior; ant, anterior. The subscript indicates the orientation. T, Tangential magnetic field (dashed lines); R, radial magnetic field (solid lines). A spike appears at 35 ms with a maximum peak-to-peak amplitude of 15 pT. The vertical scale is identical to that of the SQUID data in **a**, **b**. Extended Data Figures 1-1 and 2-1 show more details on spikes involving medial temporal structures.

Discussion

In this study we report the results of two sets of simulta-
neous intracerebral recordings of IEDs, one with ⁴He-
OPM-MEG and the other one with SQUID-MEG. Using

SEEG as reference and ground truth, we correlated the
IEDs recorded by intracerebral electrodes with both ⁴He-
OPM and SQUID sensors. We obtained the first direct
validation of the ability of ⁴He-OPM sensors to record

<figure><img src="image_5.png" /></figure>

**Figure 4.** Averaged Type II spikes (n = 3). Arrangement of figure is similar to that of Figure 3. **a**, Bipolar averaged SEEG data for the SQUID session. The spikes involve a larger network encompassing TP', TB' and B' electrodes (medial and lateral temporal pole, anterior hippocampus, third anterior temporal gyrus). TB', Rhinal cortex. **b**, Simultaneous SQUID-MEG averaged data on the four sensors closest to the ⁴He-OPM sensors; a very faint spike appears at ~75 ms with a peak-to-peak amplitude of 4 pT. **c**, Bipolar averaged SEEG data for the OPM session. Note the similarity between the intracerebral spikes of the two sessions, disclosing the same anatomic location and time course. **d**, ⁴He-OPM-MEG averaged data collected on four sensors, T<sub>pos</sub>, C<sub>ant</sub>, C<sub>post</sub>, T<sub>ant</sub>. pos, Posterior; ant, anterior. The subscript indicates the orientation. T, Tangential magnetic field (dashed lines); R, radial magnetic field (solid lines). A spike appears at 75 ms with a peak-to-peak amplitude of 9.3 pT maximum. The vertical scale is identical to that of the SQUID data in **a**, **b**. Extended Data Figures 1-1 and 2-1 show more details on spikes involving medial temporal structures.
epileptic activities, and we demonstrate that the new ⁴He-OPM system performs better in recording IEDs than the SQUID-MEG system, as evidenced by its higher SNR. Notably, these recordings have been achieved in a regular clinical environment without advanced noise correction. Because of their native 3D measurement of the magnetic field, OPM signals disclose variations in time courses, polarities, and amplitudes between the radial and tangential components of the recorded activity.

Simultaneous acquisition of MEG and intracerebral EEG is a technical feat (Dalal et al., 2009; Santiuste et al., 2008; Dubarry et al., 2014) but can now be performed without major difficulties (Badier et al., 2017). It has several key advantages in comparison with separate acquisition. It allows capturing the same activity at the surface and in depth, avoiding potential differences in brain state and medication. This is particularly important for IEDs that are spontaneous events that can widely vary in extent from one event to the other (Badier and Chauvel, 1995). Simultaneity allows correlating signals across events (Dubarry et al., 2014), and in the current study, allowed finding similar events to be compared in SQUID and OPM sessions based on SEEG topography. With this reliable comparison, we can establish that ⁴He-OPMs are at least as capable of recording epileptic activity as SQUIDs.

Currently, alkali-based OPM are mainly used for MEG recordings in healthy volunteers (for review, see Brookes et al., 2022), and few studies in epileptic patients have been performed (Tierney et al., 2021; Vivekananda et al., 2020; Feys et al., 2022), with one seizure recording recently reported (Feys et al., 2023). Present results strengthen the perspective that OPMs are an accurate and valid alternative to SQUIDs.

Nonetheless, although very promising, clinical adoption of OPMs remains challenging because of some limitations of the current technology. Alkali OPMs have a limited bandwidth (1–100 Hz) and a small dynamic range (5 nT), requiring higher attenuation than a conventional shielded room for SQUID-MEG, demanding a complementary system of magnetic shielding coils to compensate for the remaining magnetic field and to reduce the cross talk. Another challenge posed by the heat dissipation by sensors must be solved, possibly with a helmet design including insulation or an air-flowing system. However, the greater the number of sensors used, the greater the amount of heat dissipated and the risk of the system becoming uncomfortable, probably requiring a more complex active cooling system. Despite these constraints and because of technological advances to overcome them, alkali-based OPMs have given numerous proofs of their good sensitivity to biomagnetic measurement (Labyt et al., 2019; Tierney et al., 2019; Brookes et al., 2022). As evidence of the rapid evolution of the technology, 90-channel OPM systems offering triaxial magnetic field detection have been shown to improve cortical coverage successfully (Boto et al., 2022; Rea et al., 2022). The additional information offered by a vectorial measurement of the magnetic field, also reported in the present study, will be of interest to better characterize the spatiotemporal dynamics of epileptic activity and to interpret clinical data (livanainen et al., 2017; Zahran et al., 2022). In our study,

the delay of the tangential component (see averaged spike II in Figure 4d) is potentially informative on the propagation of the interictal activity, as shown by comparing radial SQUID-MEG and EEG (Merlet et al., 1997).

Helium-based OPMs could overcome some of the limitations of alkali OPM as they can be placed nearest the scalp without any discomfort for the patient because they do not require thermal insulation. ⁴He-OPM are currently the only sensors to be self-compensated on their three measurement axes in a closed-loop operating mode. This ensures the highly reliable measurement of the cerebral magnetic field, with no cross-axis projection errors (Fourcault et al., 2021), and improves the stability of the scale factor of the sensors over time. This also allows for an increased dynamic range (up to 250 nT), larger than the one corresponding to the helium resonance line width, which eliminates the constraint of a strict field nulling system. In this study, ⁴He-OPMs have been used in a standard two-layer MSR in the hospital environment, known to be particularly magnetically noisy. On the other hand, in their current release, ⁴He-OPMs have a worse sensitivity (40 fT/rtHz; Fourcault et al., 2021) compared with alkali OPMs (20 fT/rtHz in tri axis mode; Boto et al., 2022). However, in a recent study in a large group of healthy subjects, the ⁴He-OPMs showed very similar results to the classical SQUID-MEG system because of their shorter distance to the brain (Gutteling et al., 2023). In our study, we observed that ⁴He-OPM can record IEDs with a better SNR compared with SQUIDs. This result, together with the recent data from 18 volunteers, shows that ⁴He-OPMs are able to deliver high-quality brain recordings (Gutteling et al., 2023).

Some limitation of this work comes from the limited number of OPM sensors. Nevertheless, the full coverage of the SQUID system allowed selecting of the colocalized SQUID sensors and ensuring a reasonable comparison. We found that the amplitude of the OPM signals was higher than that of SQUID signals, as expected from both simulations and real data (Boto et al., 2016, 2017; livanainen et al., 2017). Still, the MEG IEDs, measured with ⁴He-OPMs and SQUIDs, arose from neocortical structures. We could not show activity from deep structures, such as the hippocampus, on either type of sensor. Sensitivity to deep sources remains a challenge. It will be necessary to use a larger number of sensors to verify whether source separation techniques such as independent component analysis will allow to extract activity from deep sources as previously demonstrated (Pizzo et al., 2019).

With all its advantages, OPM technology could extend the use of MEG to many clinical and research applications. Currently, only very few clinical centers have an MEG facility. By offering more affordability, higher signal sensitivity and bandwidth, and portability without additional equipment, the OPM technology paves the way for the democratization of this unique noninvasive method of high-resolution brain exploration, with a potentially powerful impact on both clinical practice and neuroscience.

## References

* Alkawadri R, Burgess RC, Kakisaka Y, Mosher JC, Alexopoulos AV (2018) Assessment of the utility of ictal magnetoencephalography
in the localization of the epileptic seizure onset zone. JAMA Neurol 75:1264-1272.

Badier JM, Chauvel P (1995) Spatio-temporal characteristics of paroxysmal interictal events in human temporal lobe epilepsy. J Physiol Paris 89:255-264.

Badier JM, Dubarry AS, Gavaret M, Chen S, Trébuchon AS, Marquis P, Régis J, Bartolomei F, Bénar CG, Carron R (2017) Technical solutions for simultaneous MEG and SEEG recordings: towards routine clinical use. Physiol Meas 38:N118-N127.

Beato F, Belorizky E, Labyt E, Le Prado M, Palacios-Laloy A (2018) Theory of a He 4 parametric-resonance magnetometer based on atomic alignment. Phys Rev A 98:053431.

Borna A, Iivanainen J, Carter TR, McKay J, Taulu S, Stephen J, Schwindt PDD (2022) Cross-Axis projection error in optically pumped magnetometers and its implication for magnetoencephalography systems. Neuroimage 247:118818.

Boto E, Bowtell R, Krüger P, Fromhold TM, Morris PG, Meyer SS, Barnes GR, Brookes MJ (2016) On the potential of a new generation of magnetometers for MEG: a beamformer simulation study. PLoS One 11:e0157655-24.

Boto E, Meyer SS, Shah V, Alem O, Knappe S, Kruger P, Fromhold TM, Lim M, Glover PM, Morris PG, Bowtell R, Barnes GR, Brookes MJ (2017) A new generation of magnetoencephalography: room temperature measurements using optically-pumped magnetometers. Neuroimage 149:404-414.

Boto E, Holmes N, Leggett J, Roberts G, Shah V, Meyer SS, Muñoz LD, Mullinger KJ, Tierney TM, Bestmann S, Barnes GR, Bowtell R, Brookes MJ (2018) Moving magnetoencephalography towards real-world applications with a wearable system. Nature 555:657-661.

Boto E, Shah V, Hill RM, Rhodes N, Osborne J, Doyle C, Holmes N, Rea M, Leggett J, Bowtell R, Brookes MJ (2022) Triaxial detection of the neuromagnetic field using optically-pumped magnetometry: feasibility and application in children. Neuroimage 252:119027.

Brookes MJ, Leggett J, Rea M, Hill RM, Holmes N, Boto E, Bowtell R (2022) Magnetoencephalography with optically pumped magnetometers (OPM-MEG): the next generation of functional neuroimaging. Trends Neurosci 45:621-634.

Budker D, Romalis M (2007) Optical magnetometry. Nature Phys 3:227-234.

Cohen D (1972) Magnetoencephalography: detection of the brain's electrical activity with a superconducting magnetometer. Science 175:664-666.

Cohen-Tannoudji C, Dupont-Roc J, Haroche S, Laloë F (1970a) Diverses résonances de croisement de niveaux sur des atomes pompés optiquement en champ nul I. Théorie. Rev Phys Appl (Paris) 5:95-101.

Cohen-Tannoudji C, Dupont-Roc J, Haroche S, Laloë F (1970b) Diverses résonances de croisement de niveaux sur des atomes pompés optiquement en champ nul II. Applications à la mesure de champs faibles. Rev Phys Appl (Paris) 5:102-108.

Corona L, Familia E, Perry MS, Madsen JR, Bolton J, Stone SSD, Stufflebeam SM, Pearl PL, Papadelis C (2023) Non-invasive mapping of epileptogenic networks predicts surgical outcome. Brain 146:1916-1931.

Dalal SS, Baillet S, Adam C, Ducorps A, Schwartz D, Jerbi K, Bertrand O, Gamero L, Martinerie J, Lachaux JP (2009) Simultaneous MEG and intracranial EEG recordings during attentive reading. Neuroimage 45:1289-1304.

Dubarry A-S, Badier J-M, Trébuchon-Da Fonseca A, Gavaret M, Carron R, Bartolomei F, Liégeois-Chauvel C, Régis J, Chauvel P, Alario F-X, Bénar C-G (2014) Simultaneous recording of MEG, EEG and intracerebral EEG during visual stimulation: from feasibility to single-trial analysis. Neuroimage 99:548-558.

Dupont-Roc J (1970) Détermination par des méthodes optiques des trois composantes d'un champ magnétique très faible. Rev Phys Appl (Paris) 5:853-864.

Dupont-Roc J (1971) Etude théorique de diverses résonances observables en champ nul sur des atomes 'habillés' par des photons de radiofréquence. J Phys France 32:135-144.

Feys O, Corvilain P, Aeby A, Sculier C, Holmes N, Brookes M, Goldman S, Wens V, De Tiège X (2022) On-scalp optically pumped magnetometers versus cryogenic magnetoencephalography for diagnostic evaluation of epilepsy in school-aged children. Radiology 304:429-434.

Feys O, Corvilain P, Van Hecke A, Sculier C, Rikir E, Legros B, Gaspard N, Leurquin-Sterk G, Holmes N, Brookes M, Goldman S, Wens V, De Tiège X (2023) Recording of ictal epileptic activity using on-scalp magnetoencephalography. Ann Neurol 93:419-421.

Fischer MJM, Scheler G, Stefan H (2005) Utilization of magnetoencephalography results to obtain favourable outcomes in epilepsy surgery. Brain 128:153-157.

Fourcault W, Romain R, Le Gal G, Bertrand F, Josselin V, Le Prado M, Labyt E, Palacios-Laloy A (2021) Helium-4 magnetometers for room-temperature biomedical imaging: toward collective operation and photon-noise limited sensitivity. Opt Express 29:14467-14475.

Gutteling TP, Bonnefond M, Clausner T, Daligaitz S, Romain R, Mitryukovskiy S, Fourcault W, Josselin V, Le Prado M, Palacios-Laloy A, Labyt E, Jung J, Schwartz D (2023) A new generation of OPM for high dynamic and large bandwidth MEG: the ⁴He OPMs first applications in healthy volunteers. Sensors (Basel) 23:2801.

Hämäläinen M, Hari R, Ilmoniemi RJ, Knuutila J, Lounasmaa OV (1993) Magnetoencephalography theory, instrumentation, and applications to noninvasive studies of the working human brain. Rev Mod Phys 65:413-497.

Iivanainen J, Stenroos M, Parkkonen L (2017) Measuring MEG closer to the brain: performance of on-scalp sensor arrays. Neuroimage 147:542-553.

Iivanainen J, Zetter R, Grön M, Hakkarainen K, Parkkonen L (2019) On-scalp MEG system utilizing an actively shielded array of optically-pumped magnetometers. Neuroimage 194:244-258.

Johnson CN, Schwindt PDD, Weisend M (2013) Multi-sensor magnetoencephalography with atomic magnetometers. Phys Med Biol 58:6065-6077.

Labyt E, Corsi M-C, Fourcault W, Palacios Laloy A, Bertrand F, Lenouvel F, Cauffet G, Le Prado M, Berger F, Morales S (2019) Magnetoencephalography with optically pumped ⁴He magnetometers at ambient temperature. IEEE Trans Med Imaging 38:90-98.

Labyt E, Sander T, Wakai R (2022) Flexible high performance magnetic field sensors: on-scalp magnetoencephalography and other applications. Cham, Switzerland: Springer.

Merlet I, Paetau R, García-Larrea L, Uutela K, Granström ML, Mauguière F (1997) Apparent asynchrony between interictal electric and magnetic spikes. Neuroreport 8:1071-1076.

Murakami H, Wang ZI, Marashly A, Krishnan B, Prayson RA, Kakisaka Y, Mosher JC, Bulacio J, Gonzalez-Martinez JA, Bingaman WE, Najm IM, Burgess RC, Alexopoulos AV (2016) Correlating magnetoencephalography to stereo-electroencephalography in patients undergoing epilepsy surgery. Brain 139:2935-2947.

Pizzo F, Roehri N, Medina Villalon S, Trébuchon A, Chen S, Lagarde S, Carron R, Gavaret M, Giusiano B, McGonigal A, Bartolomei F, Badier JM, Bénar CG (2019) Deep brain activities can be detected with magnetoencephalography. Nat Commun 10:971.

Rampp S, Stefan H, Wu X, Kaltenhäuser M, Maess B, Schmitt FC, Wolters CH, Hamer H, Kasper BS, Schwab S, Doerfler A, Blü I, Rö Ssler K, Buchfelder M (2019) Magnetoencephalography for epileptic focus localization in a series of 1000 cases. Brain 142:3059-3071.

Fea M, Boto E, Holmes N, Hill R, Osborne J, Rhodes N, Leggett J, Rier L, Bowtell R, Shah V, Brookes MJ (2022) A 90-channel triaxial magnetoencephalography system using optically pumped magnetometers. Ann N Y Acad Sci 1517:107-124.

Santiuste M, Nowak R, Russi A, Tarancon T, Oliver B, Ayats E, Scheler G, Graetz G (2008) Simultaneous magnetoencephalography and intracranial EEG registration: technical and clinical aspects. J Clin Neurophysiol Off Neurophysiol 25:331-339.
Seymour RA, Alexander N, Mellor S, O'Neill GC, Tierney TM, Barnes GR, Maguire EA (2021) Using OPMs to measure neural activity in standing, mobile participants. Neuroimage 244:118604.

Tierney TM, Holmes N, Mellor S, López JD, Roberts G, Hill RM, Boto E, Leggett J, Shah V, Brookes MJ, Bowtell R, Barnes GR (2019) Optically pumped magnetometers: from quantum origins to multi-channel magnetoencephalography. Neuroimage 199:598-608.

Tierney TM, Levy A, Barry DN, Meyer SS, Shigihara Y, Everatt M, Mellor S, Lopez JD, Bestmann S, Holmes N, Roberts G, Hill RM, Boto E, Leggett J, Shah V, Brookes MJ, Bowtell R, Maguire EA,

Barnes GR (2021) Mouth magnetoencephalography: a unique perspective on the human hippocampus. Neuroimage 225:117443.

Vivekananda U, Mellor S, Tierney TM, Holmes N, Boto E, Leggett J, Roberts G, Hill RM, Litvak V, Brookes MJ, Bowtell R, Barnes GR, Walker MC (2020) Optically pumped magnetoencephalography in epilepsy. Ann Clin Transl Neurol 7:397-401.

Zahran S, Mahmoudzadeh M, Wallois F, Betrouni N, Derambure P, Le Prado M, Palacios-Laloy A, Labyt E (2022) Performance analysis of optically pumped ⁴He magnetometers vs. conventional SQUIDs: from adult to infant head models. Sensors (Basel) 22:3093.