Magnetic resonance imaging with an optical atomic magnetometer
Commentary
August 22, 2006
Abstract
We report an approach for the detection of magnetic resonance imaging without superconducting magnets and cryogenics: optical atomic magnetometry. This technique possesses a high sensitivity independent of the strength of the static magnetic field, extending the applicability of magnetic resonance imaging to low magnetic fields and eliminating imaging artifacts associated with high fields. By coupling with a remote-detection scheme, thereby improving the filling factor of the sample, we obtained time-resolved flow images of water with a temporal resolution of 0.1 s and spatial resolutions of 1.6 mm perpendicular to the flow and 4.5 mm along the flow. Potentially inexpensive, compact, and mobile, our technique provides a viable alternative for MRI detection with substantially enhanced sensitivity and time resolution for various situations where traditional MRI is not optimal.
Sign up for PNAS alerts.
Get alerts for new articles, or get an alert when an article is cited.
Magnetic resonance imaging (MRI) is used in fields as diverse as medicine, material science, and geology to noninvasively extract unique information from samples ranging from human tissues to solid silicone rubber to porous rocks (1, 2). Traditional MRI is limited, however, by the use of high magnetic fields. MRI generally uses Faraday induction (which at a given initial level of polarization yields a signal with magnitude proportional to the static magnetic field) to detect the changes in magnetic flux generated by the Larmor precession of the nuclear spins in the sample (3). A superconducting magnet is typically used to provide a strong magnetic field and increase signal; the resulting cost and constrained accessibility limit the applicability of MRI. Furthermore, high fields are not favorable for certain measurements, including those involving heterogeneous samples with large magnetic-susceptibility gradients. Low-field MRI would eliminate the costs and constraints of the superconducting magnet and allow experiments that are currently inhibited by high-field susceptibility issues. Measurements in low fields are typically problematic, however, because traditional Faraday detection suffers from low sensitivity at fields much less than a Tesla.
For improved sensitivity at low fields, superconducting quantum interference devices (4–6) have been used in MRI (7, 8). Unfortunately, these devices require cryogenics, which add significant cost and exclude applications where cryogenics are not readily available, such as in open fields and space. The sensitivity to a dc magnetic field is also rather poor. Furthermore, because of fairly complicated electronics, implementation of multiple-channel detection to enhance detection sensitivity and efficiency is not convenient. A different technique is needed to overcome these limitations. Recent developments in optical atomic magnetometry (9–11), which is based on magneto-optical rotation in a “sensor” atomic vapor, have made it a viable alternative for detecting weak magnetic fields, eliminating the need for cryogenics and complicated electronics (12). Applications of atomic magnetometers have been demonstrated only with simple free induction decay of nuclear magnetization (13, 14). We here report MRI by optical atomic magnetometry.
Another major factor affecting detection sensitivity is the filling factor of the sample. When the nuclear spins are spread out over a large volume, detection becomes less sensitive. For fluidic analytes, the spins can be concentrated in a detection region and the problem can be circumvented by remote detection (15–17): a sample is first prepolarized; then it enters an encoding region where large radio frequency coils that cover the whole sample area store spectroscopic or spatial information as nuclear magnetization along the longitudinal axis; downstream the sample is concentrated and its magnetization is measured with improved filling factor.
By adopting a remote-detection scheme, the sensitive detectors become immune to the encoding magnetic fields and allow measurements to be performed continuously. The possible coupling of various detection methods with remote detection is schematically illustrated in Fig. 1.
Fig. 1.
In this article, we demonstrate the application of optical atomic magnetometry to low-field MRI, with time-resolved images of water flow obtained by using remote detection.
Results and Discussion
A typical flow profile (average of 30 measurements) is shown in Fig. 2. After the encoding π pulse, the signal (the absolute magnetization difference) reaches its maximum at 1.1 s, with a time dispersion of 1.0 s (full-width half-maximum). The dispersion results from time variation in the transfer of water from the encoding channels to the thinner transport tube.
Fig. 2.
Magnetic resonance images obtained by phase encoding are demonstrated in Fig. 3. Fig. 3A shows a schematic of the encoding volume, which is composed of two parallel cylindrical channels. Fig. 3B shows the image of the cross-section (the xy plane) of the encoding volume perpendicular to the flow direction. For imaging in the plane of the flow direction, gradients Gy and Gz are applied. The resulting images are shown in Fig. 3C. The lower parts (outlet side of the encoding volume) of the images represent those portions of the sample closer to the magnetometer and therefore arriving at the detection region earlier, starting at 0.5 s. In the frame at 1.1 s, the intensity distribution shows that the detection volume is filled with the encoded water, consistent with the maximum difference signal in the flow profile of Fig. 2. Then water flows out of the detection region, with the upper parts (inlet side of the encoding volume) lasting until 2.3 s. Therefore the images demonstrate the flow behavior, temporally and spatially resolved, of the encoded water. The temporal resolution is 0.1 s. The spatial resolution of the images is 1.6 mm cross the flow and 4.7 mm along the flow direction, each being 1/5 of the overall length along the respective axis, corresponding to a detection limit of 10 μl.
Fig. 3.
The imaging resolution can be further improved. Currently, it is affected by the long time delay between samples leaving the prepolarization and arriving at the detection region, as the polarization of water relaxes exponentially during this period. The total time delay in the present experiments, 2.5 s, can be shortened to ≈0.1 s by reducing the travel volume, resulting in a detection limit of ≈1 μl of thermally polarized water, with a temporal resolution ≈0.1 s. In addition, the prepolarization field can be increased to 1 T with permanent magnets, bringing the detection limit into the nanoliter regime, corresponding to micrometer resolution in all three dimensions. With further optimization and improvement to the magnetometer, such as adding more detection channels decreasing distance between the rubidium cells and the sample for a better detection filling factor, we expect to gain another order of magnitude in sensitivity.
Our technique has several advantages over conventional MRI. It directly measures the magnetic field generated by the sample, as opposed to the conventional measurement in the transverse plane; because no radio frequency pulses are applied to the encoded sample for detection, signals can be continuously monitored. The time-resolved images illustrate this as we observe the complete filling and emptying process of the detection volume. Low-field detection also opens up avenues of investigations of samples that contain high magnetic-susceptibility gradients, such as the study of porous material with paramagnetic impurities and medical imaging in the presence of metal implants. This technique is also useful for samples that are too big to be introduced into the bore of a superconducting magnet. Furthermore, atomic magnetometers require no cryogenics, making MRI applicable for situations where no cryogenics are available, while dramatically reducing cost, size, and maintenance of the apparatus. Finally, the system can be designed to use an economical miniature vertical-cavity-surface-emitting laser, and the electronics and optics can be consolidated for multiple channel detection. We anticipate that the whole apparatus will become quite compact and deployable as a battery-powered portable device.
A small, highly sensitive MRI scanner will be suitable for remote microfluidic MRI recently developed in our laboratory by using high-field detection (18), making it possible to characterize noninvasively the flow and reactions in microchannels in low magnetic fields. The apparatus can be used as an in-line analytic tool where a large high-field instrument is impractical to perform measurements without disturbing ongoing chemical or biochemical reactions. Our results also show that the spatial information encoded as longitudinal magnetization is preserved during the course of travel and subsequently resolvable by detection of the magnetization in the detection region. This spatial resolution will be useful in low magnetic fields where chemical shifts are not resolvable: chemically different samples located in separate microchannels can be encoded spatially and still give chemical-specific information based solely on physical location. Spatial encoding also allows for the study of flow dynamics of multiple components simultaneously, without the need for chemical differentiation.
In conclusion, we have presented results for low-field MRI using an optical atomic magnetometer. With its outstanding sensitivity in low fields, mobility, and low cost, this technique is ideal for many situations in which conventional MRI is not feasible.
Experimental Procedures
Our recently constructed apparatus, details of which will be provided elsewhere, measures time-dependent nonlinear magneto-optical rotation of the polarization of a frequency-modulated laser beam interacting with a rubidium vapor (Fig. 4). A resonance is generated when the modulation frequency equals twice the Larmor precession frequency gμB/h, where B is the magnetic field, g is the known atomic gyromagnetic ratio, μ is the Bohr magneton, and h is Planck’s constant. A bias magnetic field establishes a baseline resonance frequency and defines the detection axis of the magnetometer. When polarized nuclear spins are introduced, the change in the resonance frequency is proportional to the sample field along the detection axis. The atomic magnetometer uses two rubidium cells that form a gradiometer to cancel environmental common-mode noise. The near-dc (0.1 Hz) sensitivity of this magnetometer is 80 fT/Hz1/2.
Fig. 4.
During the experiments, a continuous stream of water, driven by high-pressure nitrogen (8 atm), first spends ≈15 s in a 0.3-T magnetic field for prepolarization. It then flows through two cylindrical channels in an encoding field of 3.1 mT (B0) (corresponding to a proton Larmor frequency of 132 kHz) directed parallel to the water flow. Finally, the water is delivered through a thin transport tube to the magnetometer where the magnetization is read out. The flow scheme of water is also demonstrated in Fig. 4. The polarization adiabatically follows the local magnetic field, which is the laboratory field between regions.
To obtain the flow profile, a π pulse was applied to invert the spins in the encoding volume, causing a change in magnetic field from the baseline magnetization produced by prepolarization (Fig. 5A).
Fig. 5.
Two-dimensional images of the water flow were obtained by using a phase-encoding pulse sequence (Fig. 5B): a π/2 pulse first rotates the magnetization into the plane perpendicular to the encoding magnetic field; the resulting transverse magnetization processes in the presence of the gradient magnetic fields; a second π/2 pulse stores one component of the phase-encoded magnetization along the longitudinal axis, which can then be remotely detected by the atomic magnetometer. Both the x and y components of the transverse magnetization are measured by performing a phase cycle on the second π/2 pulse. By incrementing the gradient amplitudes along respective spatial encoding axes, the corresponding image is obtained through Fourier transformation. The imaging sequences were repeated 16 times and averaged to create the image.
Acknowledgments
We thank Josef Granwehr, Christian Hilty, and Elad Harel for helpful discussions. This work was supported by the Director of the Office of Science (Materials Sciences Division of the U.S. Department of Energy). D.B. acknowledges support by an Office of Naval Research Multidisciplinary University Research Initiative grant.
References
1
P. T. Callaghan Principles of Nuclear Magnetic Resonance Microscopy (Clarendon, New York, 1991).
2
B. Blümich NMR Imaging of Materials (Oxford Univ. Press, Oxford, 2000).
3
D. I. Hoult, R. E. Richards J. Magn. Reson. 24, 71–85 (1976).
4
R. Kleiner, D. Koelle, F. Ludwig, J. Clarke 92, 1534–1548 (2004).
5
N. Oukhanski, R. Stolz, V. Zakosarenko, H. G. Meyer Physica C 368, 166–170 (2002).
6
R. McDemont, A. H. Trabesinger, M. Mück, E. L. Hahn, A. Pines, J. Clarke Science 295, 2247–2249 (2002).
7
A. Wong-Foy, S. Saxena, A. J. Moulé, H.-M. L. Bitter, J. A. Seeley, R. McDermott, J. Clarke, A. Pines J. Magn. Reson. 157, 235–241 (2002).
8
S. K. Lee, M. Möβle, W. Myers, N. Kelso, A. H. Trabesinger, A. Pines, J. Clarke Magn. Reson. Med. 53, 9–14 (2005).
9
C. Cohen-Tannoudji, J. DuPont-Roc, S. Haroche, F. Lahoë Phys. Rev. Lett. 22, 758–760 (1969).
10
D. Budker, W. Gawlik, D. F. Kimball, S. M. Rochester, V. V. Yashchuk, A. Weis Rev. Mod. Phys. 74, 1153–1202 (2002).
11
I. M. Savukov, S. J. Seltzer, M. V. Romalis, K. L. Sauer Phys. Rev. Lett. 95, 063004 (2005).
12
I. K. Kominis, T. W. Kornack, J. C. Allred, M. V. Romalis Nature 422, 596–599 (2003).
13
V. V. Yashchuk, J. Granwehr, D. F. Kimball, S. M. Rochester, A. H. Trabesinger, J. T. Urban, D. Budker, A. Pines Phys. Rev. Lett. 93, 160801 (2004).
14
I. M. Savukov, M. V. Romalis Phys. Rev. Lett. 94, 123001 (2005).
15
A. J. Moulé, M. M. Spence, S.-I. Han, J. A. Seeley, K. L. Pierce, S. Saxena, A. Pines Proc. Natl. Acad. Sci. USA 100, 9122–9127 (2003).
16
J. A. Seeley, S.-I. Han, A. Pines J. Magn. Reson. 167, 282–290 (2004).
17
J. Granwehr, E. Harel, S.-I. Han, S. Garcia, A. Pines, P. N. Sen, Y. Q. Song Phys. Rev. Lett. 95, 075503 (2005).
18
C. Hilty, E. E. McDonnell, J. Granwehr, K. L. Pierce, S.-I. Han, A. Pines Proc. Natl. Acad. Sci. USA 102, 14960–14963 (2005).
Information & Authors
Information
Published in
Classifications
Copyright
© 2006 by The National Academy of Sciences of the USA.
Submission history
Published online: August 22, 2006
Published in issue: August 22, 2006
Keywords
Acknowledgments
We thank Josef Granwehr, Christian Hilty, and Elad Harel for helpful discussions. This work was supported by the Director of the Office of Science (Materials Sciences Division of the U.S. Department of Energy). D.B. acknowledges support by an Office of Naval Research Multidisciplinary University Research Initiative grant.
Notes
See Commentary on page 12657.
Authors
Competing Interests
Conflict of interest statement: No conflicts declared.
Metrics & Citations
Metrics
Citation statements
Altmetrics
Citations
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.
Cited by
Loading...
View Options
View options
PDF format
Download this article as a PDF file
DOWNLOAD PDFLogin options
Check if you have access through your login credentials or your institution to get full access on this article.
Personal login Institutional LoginRecommend to a librarian
Recommend PNAS to a LibrarianPurchase options
Purchase this article to access the full text.