Measurement of the Permanent Electric Dipole Moment of the Neutron.

We present the result of an experiment to measure the electric dipole moment (EDM) of the neutron at the Paul Scherrer Institute using Ramsey's method of separated oscillating magnetic fields with ultracold neutrons. Our measurement stands in the long history of EDM experiments probing physics violating time-reversal invariance. The salient features of this experiment were the use of a ^{199}Hg comagnetometer and an array of optically pumped cesium vapor magnetometers to cancel and correct for magnetic-field changes. The statistical analysis was performed on blinded datasets by two separate groups, while the estimation of systematic effects profited from an unprecedented knowledge of the magnetic field. The measured value of the neutron EDM is d_{n}=(0.0±1.1_{stat}±0.2_{sys})×10^{-26} e.cm.

New Limit on the Permanent Electric Dipole Moment of

We report results of a new technique to measure the electric dipole moment of ^{129}Xe with ^{3}He comagnetometry. Both species are polarized using spin-exchange optical pumping, transferred to a measurement cell, and transported into a magnetically shielded room, where SQUID magnetometers detect free precession in applied electric and magnetic fields. The result from a one week measurement campaign in 2017 and a 2.5 week campaign in 2018, combined with detailed study of systematic effects, is d_{A}(^{129}Xe)=(1.4±6.6_{stat}±2.0_{syst})×10^{-28} e cm. This corresponds to an upper limit of |d_{A}(^{129}Xe)|<1.4×10^{-27} e cm (95% C.L.), a factor of 5 more sensitive than the limit set in 2001.

Development of a vector-tensor system to measure the absolute magnetic flux density and its gradient in magnetically shielded rooms.

Several experiments in fundamental physics demand an environment of very low, homogeneous, and stable magnetic fields. For the magnetic characterization of such environments, we present a portable SQUID system that measures the absolute magnetic flux density vector and the gradient tensor. This vector-tensor system contains 13 integrated low-critical temperature (LTc) superconducting quantum interference devices (SQUIDs) inside a small cylindrical liquid helium Dewar with a height of 31 cm and 37 cm in diameter. The achievable resolution depends on the flux density of the field under investigation and its temporal drift. Inside a seven-layer mu-metal shield, an accuracy better than ±23 pT for the components of the static magnetic field vector and ±2 pT/cm for each of the nine components of the gradient tensor is reached by using the shifting method.

Ultra-low-noise EEG/MEG systems enable bimodal non-invasive detection of spike-like human somatosensory evoked responses at 1 kHz.

Non-invasive EEG detection of very high frequency somatosensory evoked potentials featuring frequencies up to and above 1 kHz has been recently reported. Here, we establish the detectability of such components by combined low-noise EEG/MEG. We recorded SEP/SEF simultaneously using median nerve stimulation in five healthy human subjects inside an electromagnetically shielded room, combining a low-noise EEG custom-made amplifier (4.7 nV/√Hz) and a custom-made single-channel low-noise MEG (0.5 fT/√Hz @ 1 kHz). Both, low-noise EEG and MEG revealed three spectrally distinct and temporally overlapping evoked components: N20 (<100 Hz), sigma-burst (450-750 Hz), and kappa-burst (850-1200 Hz). The two recording modalities showed similar relative scaling of signal amplitude in all three frequencies domains (EEG [10 nV] ≅ MEG [1 fT]). Pronounced waveform (peak-by-peak) overlap of EEG and MEG signals is observed in the sigma band, whereas in the kappa band overlap was only partial. A decreasing signal-to-noise ratio (SNR; calculated for n = 12.000 averages) from sigma to kappa components characterizes both, electric and magnetic field recordings: Sigma-band SNR was 12.9  ±  5.5/19.8  ±  12.6 for EEG/MEG, and kappa-band SNR at 3.77  ±  0.8/4.5  ±  2.9. High-frequency performance of a tailor-made MEG matches closely with simultaneously recorded low-noise EEG for the non-invasive detection of somatosensory evoked activity at and above 1 kHz. Thus, future multi-channel dual-mode low-noise technology could offer complementary views for source reconstruction of the neural generators underlying such high-frequency responses, and render neural high-frequency processes related to multi-unit spike discharges accessible in non-invasive recordings.

A magnetically shielded room with ultra low residual field and gradient.

A versatile and portable magnetically shielded room with a field of (700 ± 200) pT within a central volume of 1 m × 1 m × 1 m and a field gradient less than 300 pT/m, achieved without any external field stabilization or compensation, is described. This performance represents more than a hundredfold improvement of the state of the art for a two-layer magnetic shield and provides an environment suitable for a next generation of precision experiments in fundamental physics at low energies; in particular, searches for electric dipole moments of fundamental systems and tests of Lorentz-invariance based on spin-precession experiments. Studies of the residual fields and their sources enable improved design of future ultra-low gradient environments and experimental apparatus. This has implications for developments of magnetometry beyond the femto-Tesla scale in, for example, biomagnetism, geosciences, and security applications and in general low-field nuclear magnetic resonance (NMR) measurements.

Distinction between added-energy and phase-resetting mechanisms in non-invasively detected somatosensory evoked responses.

Non-invasively recorded averaged event-related potentials (ERP) represent a convenient opportunity to investigate human brain perceptive and cognitive processes. Nevertheless, generative ERP mechanisms are still debated. Two previous approaches have been contested in the past: the added-energy model in which the response raises independently from the ongoing background activity, and the phase-reset model, based on stimulus-driven synchronization of oscillatory ongoing activity. Many criteria for the distinction of these two models have been proposed, but there is no definitive methodology to disentangle them, owing also to the limited information at the single trial level. Here, we propose a new approach combining low-noise EEG technology and multivariate decomposition techniques. We present theoretical analyses based on simulated data and identify in high-frequency somatosensory evoked responses an optimal target for the distinction between the two mechanisms.

Constraints on spin-dependent short-range interaction between nucleons.

We search for a spin-dependent P- and T-violating nucleon-nucleon interaction mediated by light pseudoscalar bosons such as axions or axionlike particles. We employ an ultrasensitive low-field magnetometer based on the detection of free precession of colocated 3He and 129Xe nuclear spins using SQUIDs as low-noise magnetic flux detectors. The precession frequency shift in the presence of an unpolarized mass was measured to determine the coupling of pseudoscalar particles to the spin of the bound neutron. For boson masses between 2 and 500 µeV (force ranges between 3×1(-4) m and 10(-1) m) we improved the laboratory upper bounds by up to 4 orders of magnitude.

Towards non-invasive multi-unit spike recordings: mapping 1kHz EEG signals over human somatosensory cortex.

OBJECTIVE: Scalp-derived human somatosensory evoked potentials (SEPs) contain high-frequency oscillations (600 Hz; 'sigma-burst') reflecting concomitant bursts of spike responses in primary somatosensory cortex that repeat regularly at 600 Hz. Notably, recent human intracranial SEP have revealed also 1 kHz responses ('kappa-burst'), possibly reflecting non-rhythmic spiking summed over multiple cells (MUA: multi-unit activity). However, the non-invasive detection of EEG signals at 1 kHz typical for spikes has always been limited by noise contributions from both, amplifier and body/electrode interface. Accordingly, we developed a low-noise recording set-up optimised to map non-invasively 1 kHz SEP components. METHODS: SEP were recorded upon 4 Hz left median nerve stimulation in 6 healthy human subjects. Scalp potentials were acquired inside an electrically and magnetically shielded room using low-noise custom-made amplifiers. Furthermore, in order to reduce thermal Johnson noise contributions from the sensor/skin interface, electrode impedances were adjusted to ≤ 1 kΩ. Responses averaged after repeated presentation of the stimulus (n=4000 trials) were evaluated by spatio-temporal pattern analyses in complementary spectral bands. RESULTS: Three distinct spectral components were identified: N20 (<100 Hz), sigma-burst (450-750 Hz), and kappa-burst (850-1200 Hz). The two high-frequency bursts (sigma, kappa) exhibited distinct and partially independent spatiotemporal evolutions, indicating subcortical as well as several cortical generators. CONCLUSIONS: Using a dedicated low-noise set-up, human SEP 'kappa-bursts' at 1 kHz can be non-invasively detected and their scalp distribution be mapped. Their topographies indicate a set of subcortical/cortical generators, at least partially distinct from the topography of the 600 Hz sigma-bursts described previously. SIGNIFICANCE: The non-invasive detection and surface mapping of 1 kHz EEG signals presented here provides an essential step towards non-invasive monitoring of multi-unit spike activity.

Extension of non-invasive EEG into the kHz range for evoked thalamocortical activity by means of very low noise amplifiers.

Ultrafast electroencephalographic signals, having frequencies above 500 Hz, can be observed in somatosensory evoked potential measurements. Usually, these recordings have a poor signal-to-noise ratio (SNR) because weak signals are overlaid by intrinsic noise of much higher amplitude like that generated by biological sources and the amplifier. As an example, recordings at the scalp taken during electrical stimulation of the median nerve show a 600 Hz burst with submicro-volt amplitudes which can be extracted from noise by the use of massive averaging and digital signal processing only. We have investigated this signal by means of a very low noise amplifier made in-house (minimal voltage noise 2.7 nV Hz(-1/2), FET inputs). We examined how the SNR of the data is altered by the bandwidth and the use of amplifiers with different intrinsic amplifier noise levels of 12 and 4.8 nV Hz(-1/2), respectively. By analyzing different frequency contributions of the signal, we found an extremely weak 1 kHz component superimposed onto the well-known 600 Hz burst. Previously such high-frequency electroencephalogram responses around 1 kHz have only been observed by deep brain electrodes implanted for tremor therapy of Parkinson patients. For the non-invasive measurement of such signals, we recommend that amplifier noise should not exceed 4 nV Hz(-1/2).

The natural line width of low field nuclear magnetic resonance spectra.

The authors suggest a procedure for the determination of the natural nuclear magnetic resonance line width Delta nu of liquids using an air coil system at flux densities from 25 microT to 150 microT. Even if the line broadening caused by instrumental field inhomogeneity is much higher than Delta nu, Delta nu can be found by extrapolating the measured line width's field dependency. For pure water this procedure yielded Delta nu=0.125 Hz+/-0.005 Hz. This is shown to be consistent with the smallest line width found for this sample below 25 microT using a superconducting quantum interference device-based spectrometer.

DC-magnetoencephalography and time-resolved near-infrared spectroscopy combined to study neuronal and vascular brain responses.

The temporal relation between vascular and neuronal responses of the brain to external stimuli is not precisely known. For a better understanding of the neuro-vascular coupling changes in cerebral blood volume and oxygenation have to be measured simultaneously with neuronal currents. With this motivation modulation dc-magnetoencephalography was combined with multi-channel time-resolved near-infrared spectroscopy to simultaneously monitor neuronal and vascular parameters on a scale of seconds. Here, the technique is described, how magnetic and optical signals can be measured simultaneously. In a simple motor activation paradigm (alternating 30 s of finger movement with 30 s of rest for 40 min) both signals were recorded non-invasively over the motor cortex of eight subjects. The off-line averaged signals from both modalities showed distinct stimulation related changes. By plotting changes in oxy- or deoxyhaemoglobin as a function of magnetic field a characteristic trajectory was created, which was similar to a hysteresis loop. A parametric analysis allowed quantitative results regarding the timing of coupling: the vascular signal increased significantly slower than the neuronal signal.

Latency analysis of single auditory evoked M100 responses by spatio-temporal filtering.

Appropriate spatial filtering followed by temporal filtering is well suited for the single-trial analysis of multi-channel magnetoencephalogram or electroencephalogram recordings. This is demonstrated by the results of a single-trial latency analysis obtained for auditory evoked M100 responses from nine subjects using two different stimulation frequencies. Spatial filters were derived automatically from the data via noise-adjusted principle component analysis, and single-trial latencies were estimated from the signal phase after complex bandpass filtering. For each of the two stimulation frequencies, estimated single-trial latencies were consistent with results obtained from a standard approach using averaged evoked responses. The quality of the estimated single-trial latencies was additionally assessed by their ability to separate between the two different stimulation frequencies. As a result, more than 80% of the single trials can be classified correctly by their estimated latencies.

Demagnetization of magnetically shielded rooms.

Magnetically shielded rooms for specific high resolution physiological measurements exploiting the magnetic field, e.g., of the brain (dc-magnetoencephalography), low-field NMR, or magnetic marker monitoring, need to be reproducibly demagnetized to achieve reliable measurement conditions. We propose a theoretical, experimental, and instrumental base whereupon the parameters which affect the quality of the demagnetization process are described and how they have to be handled. It is demonstrated how conventional demagnetization equipment could be improved to achieve reproducible conditions. The interrelations between the residual field and the variability at the end of the demagnetization process are explained on the basis of the physics of ferromagnetism and our theoretical predictions are evaluated experimentally.

Cross-correlation analysis of the correspondence between magnetoencephalographic and near-infrared cortical signals.

OBJECTIVES: The study of neurovascular coupling greatly benefits from combined measurements of neuronal and vascular signals. Two-step signal processing is developed to extract parameters describing the coupling. METHODS: Using a magnetometer in an extremely well shielded room a broadband magnetoencephalogram was simultaneously measured with time-resolved near-infrared spectroscopy during a motor activity paradigm. The raw MEG and NIRS data were denoised separately using independent component analysis. RESULTS: After averaging the resulting signals showed motor activity-related changes. The temporal correspondence between MEG and NIRS was assessed plotting a combined trajectory and calculating a cross-correlation. Compared to the MEG signal, at movement onset the NIRS signal showed an onset delay in the range of seconds. CONCLUSIONS: Multi-variate signal pre-processing followed by temporal delay estimates demonstrated the extraction of neurovascular coupling parameters.

A template-free approach for determining the latency of single events of auditory evoked M100.

The phase of the complex output of a narrow band Gaussian filter is taken to define the latency of the auditory evoked response M100 recorded by magnetoencephalography. It is demonstrated that this definition is consistent with the conventional peak latency. Moreover, it provides a tool for reducing the number of averages needed for a reliable estimation of the latency. Single-event latencies obtained by this procedure can be used to improve the signal quality of the conventional average by latency adjusted averaging.

Systematic latency variation of the auditory evoked M100: from average to single-trial data.

Standard analyses of neurophysiologically evoked response data rely on signal averaging across many epochs associated with specific events. The amplitudes and latencies of these averaged events are subsequently interpreted in the context of the given perceptual, motor, or cognitive tasks. Can such critical timing properties of event-related responses be recovered from single-trial data? Here, we make use of the M100 latency paradigm used in previous magnetoencephalography (MEG) research to evaluate a novel single-trial analysis approach. Specifically, the latency of the auditory evoked M100 varies systematically with stimulus frequency over a well-defined time range (lower frequencies, e.g., 125 Hz, yield up to 25 ms longer latencies than higher frequencies, e.g., 1000 Hz). Here, we show that the complex filtering approach to single-trial analysis recovers this key characteristic of the M100 response, as well as some other important response properties relating to lateralization. The results illustrate (i) the utility of the complex filtering method and (ii) the potential of the M100 latency to be used for stimulus encoding, since the relevant variation can be observed in single trials.

Technique for the direct measurement of DC-like magnetic biosignals demonstrated by the cold reflex of the abdomen.

Very low frequency dc-like signals, such as the cold reflex, could only be measured up to now by moving the subject repeatedly, up to the magnetic detector. PTB's novel magnetically shielded room BMSR 2, together with a low noise 16 channel SQUID magnetometer, allow the recording of dc-like signals without moving the subject; these are direct measurements. The total observed magnetic drifts are limited by 1/f-noise and external disturbances to a value below 6 pT/h. The measurement is continuous in time, therefore provides frequency resolution from dc to several kHz. This allows us to also observe the changing pattern between two different static magnetic states. As an example, the measurement of the cold reflex of the abdomen is shown and discussed. Not only the expected cold reflex, but other periodic and spontaneous signals from the human body can be seen with this method.

Spatial distribution of cardiac magnetic vector fields acquired from 3120 SQUID positions.

An extended measurement of the magnetic vector field of the human heart is presented. It is acquired by sequential recordings, shifting a 16 SQUID vector magnetometer across 195 positions over a healthy subject's thorax. The magnetocardiographic (MCG) signals were synchronized using a simultaneously measured ECG channel. The registration of the field extends over a volume of 1000 mm x 600 mm x 420 mm sampled at 3120 SQUID positions. We present diagrams of the vector amplitude of selected points in 6 planes at increasing distances from the frontal thorax. Each plane contains 76 vector points. Additionally, we measured the vector field at 126 points lateral to the chest. At the edge points of the measurement volume, the absolute value of the magnetic vector signal amplitude exceeds 0.3 pT in all measurement points. The dataset provides an excellent base to study dedicated MCG detection or rejection methods. Examples where rejection of the heart signal is necessary are magnetoencephalography, magnetoneurography and fetal MCG. The knowledge of the spatio-temporal distribution of the magnetic vector field of the heart supports the development and comparison of multi-SQUID systems and will be used to create new MCG interpretation and representation algorithms.

Separation of fetal and maternal magnetocardiographic signals in twin pregnancy using independent component analysis (ICA).

The identification of fetal and maternal signals in magnetocardiograms (MCG) is central to data preprocessing and a prerequisite for data analysis and assessment. This is usually done by creating a template of the signal to be identified and marking data segments correlating to this template before averaging. This procedure is not only cumbersome, but may also lead to problems when there are several overlapping signals of interest such as in MCG recording in single or, more so, in twin pregnancy. Independent component analysis (ICA), which uses higher order statistics to decompose the signal into statistical independent components, has already been used in single pregnancies to distinguish between maternal and fetal signals. We applied the ICA algorithm TDSEP to 9 data sets of twin pregnancies acquired between the 28th and 38th week of pregnancy. Resulting ICA components can be used for further data analysis, e.g., for finding robust triggers or estimating the heart rate and its variability of the twins. The results showed that the maternal and fetal components can be separated from each other as well as from other sources of noise and artifacts. Differences between averaged ICA time curves and averaged raw data are not significant. Limitations include a concurrence of heart rates and changes in signal morphology due to gross movement. Nonetheless, ICA offers a fast and efficient approach for the preprocessing of MCGs with multiple signals of interest.

Discrimination of multiple sources using a SQUID vector magnetometer.

For many biomagnetic applications the discrimination between simultaneously active sources is required. To evaluate the performance of a given SQUID system in this respect, the angle between the signal vectors of different sources is used. If the angle reaches large values, discrimination between the multiple sources is possible. We tested this approach with the first module of a new vector magnetometer system consisting of 19 identical modules. Two examples of measurements illustrate the differentiation of multiple sources, i.e. the fetal and the mother's heart signal, and alpha rhythm and heart signal in MEG recordings. This first module of a vector magnetometer system containing 16 SQUIDs is operated at PTB in the new Berlin Magnetically Shielded Room (BMSR 2) The spatial configuration of the 16 integrated SQUID magnetometers of the module is such that all three vector components of the magnetic field can be calculated in three measurement planes at 1.5 cm, 5 cm, and 10.5 cm above the Dewar bottom, respectively. The SQUID magnetometer channels have a typical white noise level of less than 2.3 fT/square root of Hz1/2 at 1 kHz.