Magnetic resonance imaging with an optical atomic magnetometer

August 22, 2006
103 (34) 12668-12671
Commentary
MRI without the magnet
Daniel Raftery

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.
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 (46) 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 (911), 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 (1517): 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.
Various prepolarization and detection methods for MRI in combination with remote detection. To overcome the limited nuclear polarization in low magnetic fields, prepolarization methods can be applied, including spin-exchange optical pumping by lasers, polarization by cryogenics, and thermal polarization by permanent magnets. In the encoding stage, excitation pulses and gradient pulses are applied to store spatial information as nuclear spin magnetization. For detection, the methods of optical atomic magnetometer, radio frequency (RF) coil, and superconducting quantum interference device (SQUID) are shown.
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.
Flow profile of water measured by magnetization inversion by a single π pulse.
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.
Magnetic resonance images. (A) The encoding volume. The two channels are 3.2 mm in diameter and 25 mm long each, with a center-to-center spacing of 5.1 mm. (B) Image of the cross-section of the encoding volume perpendicular to the flow (xy plane) at t = 1.1 s. (C) Time-resolved images in the yz plane. Measurements were obtained with a time interval of 0.1 s. All of the images are color-mapped at the same scale, as indicated below the images. The total experimental time for these flow images is 12 h, which is dominated by the waiting time between measurements to allow the sample from the previous measurement cycle to clear the system. The overall time will be reduced to minutes with shorter travel distances.
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.
Schematic of the experimental setup. The magnetometer consists of two rubidium (87Rb) vapor cells forming a first-order gradiometer. The cells are cubic, with a 1-cm side length, and maintained at 43°C. The beam from a single laser (whose wavelength is resonant with the rubidium D1 transition) is split equally into two, one for each cell. For each arm of the gradiometer, polarizing and analyzing prisms are oriented at 45° to each other to detect optical rotation occurring in the vapor cell. The intrinsic resonance linewidth is ≈5 Hz. A piercing solenoid provides a 0.05-mT leading field (Bl). The geometry is such that the rubidium atoms are not subject to this field. A bias field of 70 nT (Bb) gives a resonance frequency of ≈1,000 Hz for the modulation of the laser when no sample is introduced. The magnetometer measures the magnetic field change along the direction of the bias field. The magnetized sample in the detection region produces magnetic fields of opposite direction in the two cells. The frequency of magneto-optical resonance on one arm of the gradiometer is fed back to the laser modulation to maintain this arm on resonance. Thus the optical rotation in the other cell represents the difference field between the two cells, equal to twice the magnetic field produced by the sample. N, north pole; S, south pole; PD, photo diode; PP, polarization prism; MS, magnetic shield; BS, beam splitter.
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.
Pulse sequences used in the experiments. (A) A π pulse used to invert the spins for measuring the flow profile. The pulse duration is 116 μs. The detection timing of the magnetometer is relative to the π pulse. (B) Phase-encoding pulse sequence for magnetic resonance imaging. The pulse duration is 58 μs for the π/2 pulses and 1.2 ms for the phase-encoding gradient pulses (GPE). The gradient step sizes were 0.3, 0.22, and 0.1 G/cm, for axes x, y, and z (flow direction), respectively. The number of steps are 10, 14, and 10 for x, y, and z, respectively.
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

Go to Proceedings of the National Academy of Sciences
Proceedings of the National Academy of Sciences
Vol. 103 | No. 34
August 22, 2006
PubMed: 16885210

Classifications

Submission history

Published online: August 22, 2006
Published in issue: August 22, 2006

Keywords

  1. low field
  2. remote detection

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

Affiliations

Shoujun Xu
Materials Sciences Division,
Departments of Chemistry and
Valeriy V. Yashchuk
Advanced Light Source, Lawrence Berkeley National Laboratory, Berkeley, CA 94720; and
Marcus H. Donaldson
Materials Sciences Division,
Departments of Chemistry and
Simon M. Rochester
Physics, University of California, Berkeley, CA 94720
Dmitry Budker
Nuclear Sciences Division, and
Physics, University of California, Berkeley, CA 94720
Alexander Pines [email protected]
Materials Sciences Division,
Departments of Chemistry and

Notes

To whom correspondence should be addressed. E-mail: [email protected]
Contributed by Alexander Pines, June 28, 2006
Author contributions: S.X., V.V.Y., S.M.R., D.B., and A.P. designed research; and S.X. and M.H.D. performed research.

Competing Interests

Conflict of interest statement: No conflicts declared.

Metrics & Citations

Metrics

Note: The article usage is presented with a three- to four-day delay and will update daily once available. Due to ths delay, usage data will not appear immediately following publication. Citation information is sourced from Crossref Cited-by service.


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 PDF

    Login options

    Check if you have access through your login credentials or your institution to get full access on this article.

    Personal login Institutional Login

    Recommend to a librarian

    Recommend PNAS to a Librarian

    Purchase options

    Purchase this article to access the full text.

    Single Article Purchase

    Magnetic resonance imaging with an optical atomic magnetometer
    Proceedings of the National Academy of Sciences
    • Vol. 103
    • No. 34
    • pp. 12655-12951

    Media

    Figures

    Tables

    Other

    Share

    Share

    Share article link

    Share on social media