Juhani Dabek

Hybrid ultra-low-field MRI and magnetoencephalography system based on a commercial whole-head neuromagnetometer

Ultra‐low‐field MRI uses microtesla fields for signal encoding and sensitive superconducting quantum interference devices for signal detection. Similarly, modern magnetoencephalography (MEG) systems use arrays comprising hundreds of superconducting quantum interference device channels to measure the magnetic field generated by neuronal activity. In this article, hybrid MEG‐MRI instrumentation based on a commercial whole‐head MEG device is described. The combination of ultra‐low‐field MRI and MEG in a single device is expected to significantly reduce coregistration errors between the two modalities, to simplify MEG analysis, and to improve MEG localization accuracy. The sensor solutions, MRI coils (including a superconducting polarizing coil), an optimized pulse sequence, and a reconstruction method suitable for hybrid MEG‐MRI measurements are described. The performance of the device is demonstrated by presenting ultra‐low‐field‐MR images and MEG recordings that are compared with data obtained with a 3T scanner and a commercial MEG device. Magn Reson Med, 2013.

Avoiding eddy-current problems in ultra-low-field MRI with self-shielded polarizing coils.

In ultra-low-field magnetic resonance imaging (ULF MRI), superconductive sensors are used to detect MRI signals typically in fields on the order of 10-100 μT. Despite the highly sensitive detectors, it is necessary to prepolarize the sample in a stronger magnetic field on the order of 10-100 mT, which has to be switched off rapidly in a few milliseconds before signal acquisition. In addition, external magnetic interference is commonly reduced by situating the ULF-MRI system inside a magnetically shielded room (MSR). With typical dipolar polarizing coil designs, the stray field induces strong eddy currents in the conductive layers of the MSR. These eddy currents cause significant secondary magnetic fields that may distort the spin dynamics of the sample, exceed the dynamic range of the sensors, and prevent simultaneous magnetoencephalography and MRI acquisitions. In this paper, we describe a method to design self-shielded polarizing coils for ULF MRI. The experimental results show that with a simple self-shielded polarizing coil, the magnetic fields caused by the eddy currents are largely reduced. With the presented shielding technique, ULF-MRI devices can utilize stronger and spatially broader polarizing fields than achievable with unshielded polarizing coils.

All-planar SQUIDs and pickup coils for combined MEG and MRI

Flux trapping and random flux movement are common problems in superconducting thin-film devices. Ultrasensitive magnetic field sensors based on superconducting quantum interference devices (SQUIDs) coupled to large pickup coils are especially vulnerable to strong external fields. The issue has become particularly relevant with the introduction of SQUID-based ultra-low-field (ULF) nuclear magnetic resonance (NMR) and magnetic resonance imaging (MRI) techniques. In this paper, we study the constraints of thin-film-based magnetometers and gradiometers as exposed to magnetic field sequences of ULF MRI. In particular, we address issues such as response recovery, transient noise, magnetization and behaviour under shielded room conditions after prepolarization. As a result, we demonstrate sensors that are suitable for a combined multi-channel magnetoencephalography (MEG) and MRI imaging system.

The Spatial and Temporal Distortion of Magnetic Fields Applied Inside a Magnetically Shielded Room

In the context of biomagnetism, a magnetically shielded room (MSR) is designed for shielding against external magnetic fields. Recently, several applications, such as combined structural magnetic resonance imaging (MRI) and functional magnetoencephalography (MEG), have emerged that require applying relatively strong magnetic fields inside the MSR. These magnetic fields induce eddy currents and magnetize the MSR walls that are made of materials with high permeability and conductivity. These eddy currents and magnetization generate secondary magnetic fields inside the room that disturb, e.g., combined MEG-MRI by affecting sample spins and by exceeding the available dynamic range of the magnetic field sensors. In this work, static and dynamic magnetic fields applied inside an MSR are studied both theoretically and experimentally. Using a spherical model, analytical expressions are derived for the amplitudes and time constants of the various secondary magnetic field modes. These predictions are validated by comparison with experimental measurements in a rectangular MSR. The results of this study facilitate, e.g., the design of coils compatible with an MSR; a self-shielded coil is presented that decreases the secondary magnetic fields to a small fraction.

Temperature dependence of relaxation times and temperature mapping in ultra-low-field MRI

Ultra-low-field MRI is an emerging technology that allows MRI and NMR measurements in microtesla-range fields. In this work, the possibilities of relaxation-based temperature measurements with ultra-low-field MRI were investigated by measuring T1 and T2 relaxation times of agarose gel at 50 μT-52 mT and at temperatures 5-45°C. Measurements with a 3T scanner were made for comparison. The Bloembergen-Purcell-Pound relaxation theory was combined with a two-state model to explain the field-strength and temperature dependence of the data. The results show that the temperature dependencies of agarose gel T1 and T2 in the microtesla range differ drastically from those at 3T; the effect of temperature on T1 is reversed at approximately 5 mT. The obtained results were used to reconstruct temperature maps from ultra-low-field scans. These time-dependent temperature maps measured from an agarose gel phantom at 50 μT reproduced the temperature gradient with good contrast.

Suppressing Multi-Channel Ultra-Low-Field MRI Measurement Noise Using Data Consistency and Image Sparsity

Ultra-low-field (ULF) MRI (B 0 = 10–100 µT) typically suffers from a low signal-to-noise ratio (SNR). While SNR can be improved by pre-polarization and signal detection using highly sensitive superconducting quantum interference device (SQUID) sensors, we propose to use the inter-dependency of the k-space data from highly parallel detection with up to tens of sensors readily available in the ULF MRI in order to suppress the noise. Furthermore, the prior information that an image can be sparsely represented can be integrated with this data consistency constraint to further improve the SNR. Simulations and experimental data using 47 SQUID sensors demonstrate the effectiveness of this data consistency constraint and sparsity prior in ULF-MRI reconstruction.

SQUID-sensor-based ultra-low-field MRI calibration with phantom images: towards quantitative imaging.

In ultra-low-field magnetic resonance imaging (ULF MRI), measured resonance signals oscillate at Larmor frequencies around 1 kHz compared to even above 100 MHz in high-field MRI. Thus, detection by induction coils in ULF MRI is not feasible, whereas superconducting quantum interference device (SQUID) sensors can measure these femtotesla-level signals. The signal-to-noise ratio is enhanced by prepolarization in a field that is typically 100-1000 times higher than the field during acquisition. Based on both measurements and simulations, a procedure for calibrating a SQUID-sensor-based MRI system with MR images is presented in this article. Magnetoencephalography (MEG) can be integrated with ULF MRI, and may also benefit from such a calibration procedure. Conventionally, electromagnet probe signals have been used for the SQUID-sensor calibration in MEG; the presented ULF-MRI-based approach using an imaging phantom could replace this procedure in hybrid MEG-MRI or ULF MRI alone. The necessary theory is provided here with experimental verification. The calibration procedure opens the possibility of performing quantitative ULF MRI without sample-specific reference scans.

Low field MRI with magnetoresistive mixed sensors

Low-field magnetic resonance imaging (MRI) requires very sensitive magnetometers for detection of weak magnetic signals in femtoTesla range. Magnetoresistive mixed sensors could be used as the sensitive magnetometers [1]. These sensors are the combination of a Giant Magnetoresistive (GMR) field sensor and a superconducting loop used as a flux-to-field transformer. A MRI set-up based on mixed sensor has been mounted. This system operates in static field up to 10mT without shielded room. We will present the LF-MRI setup with the implementation of mixed sensors and the first MRI images obtained at low field without prepolarization.

Efficient Concomitant and Remanence Field Artifact Reduction in Ultra-Low-Field MRI Using a Frequency-Space Formulation

For ultra‐low‐field MRI, the spatial‐encoding magnetic fields generated by gradient coils can have strong concomitant fields leading to prominent image distortion. Additionally, using superconducting magnet to pre‐polarize magnetization can improve the signal‐to‐noise ratio of ultra‐low‐field MRI. Yet the spatially inhomogeneous remanence field due to the permanently trapped flux inside a superconducting pre‐polarizing coil modulates magnetization and causes further image distortion.

Improved determination of FID signal parameters in low-field NMR

In this work, novel methods are suggested for assessing signal parameters of the free induction decay (FID) in nuclear magnetic resonance (NMR) experiments. The FID signal was recorded in a microtesla field and analysed to determine its relaxation time, amplitude, Larmor frequency and phase. The challenge was posed by the narrow line width, whose related effects were investigated through simulations, also. The developed methods give a new view on FID signal estimation in microtesla as well as lower and higher fields. It is shown that the transverse relaxation time of a sample can be accurately determined in the frequency domain by other means than the Lorentz peak half width. Also, with some realistic approximations, a simple functional form for the power spectrum Lorentz peak shape is proposed. As shown in this work, the inspection of the power spectrum instead of the absorption and dispersion Lorentzians is advantageous in the sense that the waveform is independent of the FID phase. The automatic and efficient methods presented in this work incorporate an integral exponential fit, the fit of the power spectrum Lorentz peak and two ways to determine the FID phase. When there are sufficiently many data points in the Lorentz peak, the power spectrum Lorentz peak shape fit provides a quick, simple and accurate way of determining the amplitude, relaxation time and Larmor frequency of the FID. In the measurements of this work, however, the narrow line width led to establishing a more applicable method which is based on the exponential decay of the Lorentz peak with a temporally moving power spectrum window.