New Research In
Physical Sciences
Social Sciences
Featured Portals
Articles by Topic
Biological Sciences
Featured Portals
Articles by Topic
- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology
Nanoscale magnetic resonance imaging
Communicated by Stuart S. P. Parkin, IBM Almaden Research Center, San Jose, CA, December 1, 2008 (received for review August 21, 2008)

Abstract
We have combined ultrasensitive magnetic resonance force microscopy (MRFM) with 3D image reconstruction to achieve magnetic resonance imaging (MRI) with resolution <10 nm. The image reconstruction converts measured magnetic force data into a 3D map of nuclear spin density, taking advantage of the unique characteristics of the “resonant slice” that is projected outward from a nanoscale magnetic tip. The basic principles are demonstrated by imaging the 1H spin density within individual tobacco mosaic virus particles sitting on a nanometer-thick layer of adsorbed hydrocarbons. This result, which represents a 100 million-fold improvement in volume resolution over conventional MRI, demonstrates the potential of MRFM as a tool for 3D, elementally selective imaging on the nanometer scale.
Magnetic resonance imaging (MRI) is well-known in medicine and in the neurosciences as a powerful tool for acquiring 3D morphological and functional information with resolution in the millimeter-to-submillimeter range (1, 2). Unfortunately, despite considerable effort, attempts to push the spatial resolution of conventional MRI into the realm of high-resolution microscopy have been stymied by fundamental limitations, especially detection sensitivity (3, 4). Consequently, the highest resolution MRI microscopes today remain limited to voxel volumes >40 μm3 (5–8). The central issue is that MRI is based on the manipulation and detection of nuclear magnetism, and nuclear magnetism is a relatively weak physical effect. It appears that conventional coil-based inductive detection techniques simply cannot provide adequate signal-to-noise ratio for detecting voxel volumes below the micrometer size. This sensitivity constraint is unfortunate because MRI has much to offer the world of microscopy with its unique contrast modalities, its elemental selectivity, and its avoidance of radiation damage.
Despite the many challenges, there is strong motivation to extend MRI to finer resolution, especially if the nanometer scale can be reached. At the nanometer scale, one might hope to directly and nondestructively image the 3D structure of individual macromolecules and molecular complexes (9). Such a powerful molecular imaging capability could be of particular interest to structural biologists trying to unravel the structure and interactions of proteins, especially for those proteins that cannot be crystallized for X-ray analysis, or are too large for conventional NMR spectroscopy. Nanoscale MRI, with its capacity for true 3D, subsurface imaging, its potential for generating contrast by selective isotopic labeling and its nondestructive nature, would be a welcome complement to the characteristics of electron microscopy. The key to pushing MRI to the nanoscale is detection sensitivity.
Recently, a significant breakthrough in magnetic resonance detection sensitivity has been achieved by using magnetic resonance force microscopy (MRFM) (9–13), resulting in single spin detection for electrons (14) and substantial progress in nuclear spin detection (15–24). Despite the great progress in nuclear spin MRFM, only one previous nanoscale imaging experiment has been demonstrated, and it was limited to 90-nm resolution in 2 dimensions for 19F nuclei in an inorganic test sample (25). Here, we report that MRFM can perform 3D MRI of 1H nuclear spins (protons) in a biological specimen [tobacco mosaic virus (TMV) particles] with a spatial resolution down to 4 nm. This capability is enabled by several key technical advances, including the generation of magnetic field gradients as high as 4 million Tesla (T) per meter, detailed understanding of the MRFM point-spread function, and application of an image-reconstruction technique capable of converting magnetic force measurements into a 3D map of proton density.
Principles
MRFM is based on mechanical measurement of ultrasmall (attonewton) magnetic forces between nuclear spins in a sample and a nearby magnetic tip. Basic elements of our MRFM apparatus are shown in Fig. 1. The test sample consists of individual TMV particles that are deposited onto the flat end of an ultrasensitive silicon cantilever. The end of the cantilever is positioned close to a 200-nm-diameter magnetic tip that produces a strong and very inhomogeneous magnetic field. The magnetic tip sits on a copper “microwire” that serves to efficiently generate a radiofrequency (rf) magnetic field that excites NMR (26). Frequency modulation of the rf field induces periodic inversions of the 1H spins in the sample, resulting in a periodic force that drives the mechanical resonance of the cantilever. Monitoring the cantilever oscillation amplitude while mechanically scanning the magnetic tip with respect to the sample in 3 dimensions provides data that allow the reconstruction of the 1H density. The imaging is performed in vacuum and at low temperature (T = 300 mK).
Configuration of MRFM apparatus. (A) Tobacco mosaic virus particles, attached to the end of an ultrasensitive silicon cantilever, are positioned close to a magnetic tip. A rf current irf passing through a copper microwire generates an alternating magnetic field Brf that induces magnetic resonance in the 1H spins of the virus particles. The resonant slice represents those points in space where the field from the magnetic tip (plus an external field) matches the condition for magnetic resonance. Three-dimensional scanning of the tip with respect to the cantilever, followed by image reconstruction is used to generate a 3D image of the spin density in the virus sample. (B) Scanning electron micrograph of the end of the cantilever. Individual tobacco mosaic virus particles are visible as long, dark rods on the 0.8-μm × 1.3-μm-sized sample platform. (C) Detail of the magnetic tip.
The TMV particles are deposited onto the cantilever in solution and then air dried. (See supporting information (SI) Appendix for preparation details.) As shown in Fig. 1B, the sample consists of both whole virus and smaller fragments. The TMV particles, which have a rod-like geometry with diameter of 18 nm and lengths up to 300 nm (27, 28), were chosen as test objects because they are physically robust and have a size suitable for evaluating our imaging resolution. They also serve to demonstrate that MRFM is capable of imaging native biological specimens. Approximately 95% of the virus mass consists of protein, resulting in a 1H density estimated to be ρ = 4 × 1028 spins per m3. In the future, rapid freezing techniques, such as used in cryoelectron microscopy, could be used to better preserve the structural integrity of fully hydrated biological samples (29, 30).
NMR will only occur if the 1H spins in the sample are at the correct field for satisfying the Larmor resonance condition: B0(r) = ω0/γ ≡ Bres, where ω0 is the rf field frequency, and γ = 2π × 42.57 MHz/T is the proton gyromagnetic ratio. The field B0(r) ≡ ∣Bextz^ + Btip(r)∣ is supplied by the combination of the field from an external superconducting magnet, Bext, and the field from the magnetic tip, Btip(r), where r is the position with respect to the tip apex. For an rf center frequency of ω0 = 2π × 114.8 MHz, the resonance condition is met for B0(r) = 2.697 T. Because the field from the magnetic tip is a strong function of position, the resonance is confined to a thin, approximately hemispherical “resonant slice” that extends outward from the tip (Figs. 1A and 2). The field gradient at the resonant slice can exceed 4 × 106 T/m at a distance of 25 nm from the tip, resulting in a slice thickness that is as thin as a few nanometers. The rf field is frequency modulated with a peak deviation of Δωrf,peak = 2π × 600 kHz in order to drive adiabatic inversions of the protons. The periodicity of the spin inversion is chosen to match the mechanical resonance of the cantilever (≈2.9 kHz). In the presence of the field gradient from the magnetic tip, the spin inversions generate a small oscillating force, typically on the order of 10 aN-rms, that excites a slight (subangstrom) vibration of the cantilever. The vibration is detected by a fiber-optic interferometer and lock-in amplifier.
The spin signal originates from the naturally occurring
Three-dimensional imaging of the sample requires 2 steps: data collection and image reconstruction. First, the spin signal is measured as the magnetic tip is mechanically scanned with respect to the sample in a 3D raster pattern, yielding a map of the spin signal as a function of tip position. Because of the extended geometry of the resonant slice, however, a spatial scan does not directly produce a map of the proton distribution in the sample. Instead, each data point in the scan contains spin signal contributions from a variety of depths and lateral positions. Specifically, the map ξ(rs) of the signal as a function of tip scan position rs is related to the proton distribution ρ(r) by the convolution integral
where K(r) is the 3D point spread function (PSF) associated with the resonant slice. K(r) is defined as the mean spin signal generated by a randomly polarized spin in the sample at a position r with respect to the magnetic tip.
The amplitude of the MRFM point-spread function is set by 2 main factors: one that determines the effectiveness of the spin inversions and confines the response to the near vicinity of the resonant slice and one that reflects the strength of the lateral field gradient at the position of the spin. As discussed in the SI Appendix, we find that the PSF can be approximated by
Here, G(r) ≡ ∂Btip,z/∂x is the lateral field gradient from the tip, ΔB0(r) ≡ B0(r) − ω0/γ is the off-resonance condition, and Δωrf,peak/γ ≈ 14 mT is the peak FM modulation converted to magnetic field units. The proportionality constant A depends on details of the experiment, such as the correlation time of the spin inversions and the bandwidth of the detection. G(r) and B0(r) are key components of K(r), and both require detailed knowledge of the field produced by the magnetic tip. As discussed in SI Appendix, we calculate Btip(r) using a magnetostatic model of the tip and then tune the parameters of the model (for example, tip magnetization and geometry) to be consistent with the measured scan data. Based on Eq. 2 and our best estimates of the magnetic tip parameters, we obtain the PSF shown in Fig. 2. At a distance of 24 nm, where the peak gradient is G = 4.2 mT/nm, we find the shell thickness (full width at half maximum) to be as thin as
Details of the resonant slice and associated point spread function (PSF). (A) Three-dimensional representation corresponding to the conditions Bres = ω0/γ = 2.697 T and Bext = 2.482 T. The center of the tip apex is assumed to be at the origin of the coordinate system. The resonant slice is the hemispherical “shell” outlined in red, representing the region of space for which B0(r) lies within Bres ± Δωrf,peak/γ (here 2.697 ± 0.014 T). Regions to the left and right of the tip (shaded red) contribute most to the signal because this is where the lateral gradient G(r) is largest. (B) Cross-sections of the point spread function at the 4 tip-sample spacings used in the imaging experiment. The PSF was calculated by using Eq. 2, assuming A = 1. The color scale reflects the force variance per spin (zN = zeptonewton = 10−21 N). The size of the tip apex (ra = 100 nm) is indicated by a dotted circle. At z = 24 nm, the PSF lobe thickness reaches a minimum of ≈4.8 nm (FWHM). The small left–right asymmetry is due to a slight (1.7°) tilt of the sample plane with respect to the magnetic tip.
To recover the real-space proton distribution ρ(r) from the 3D scan data ξ(rs), the effect of the PSF must be deconvolved. We use an iterative Landweber algorithm that starts with an initial estimate for the spin density of the object, ρ0(r) and then improves the estimate successively by using the following steps (34, 35):
where ρn(r) is the reconstructed spin density after n iterations, δξn(rs) is the difference between the measured and predicted spin signal maps (the error to be minimized), and α(r) controls the rate of convergence. Because spin density should always be ≥ 0, we enforce this condition by setting any negative values of ρn(r) to zero after each iteration step. The iterations typically proceed for a few thousand steps until the residual error becomes comparable with the measurement noise. In the future, the implementation of more sophisticated image reconstruction algorithms may be advantageous (36).
Results
Fig. 3A shows 3D scan data organized as x–y scans taken at 4 different tip-sample spacings: d = 24, 37, 50, and 62 nm. Each of the approximately 50-nm × 50-nm scan areas contains 60 × 32 data points acquired with an acquisition time of 1 min per point. Peak signal strength ranged from ≈600 aN2 for the closest scan plane to ≈90 aN2 at the most distant plane. The uncertainty of the spin signal estimate (the standard error of the measurement) ranges from ≈40 aN2 for the d = 24 nm scan plane to ≈15 aN2 for d = 62 nm. Thus, at the positions of maximum signal, the signal-to-noise ratio (SNR) varied from ≈15 for the closest scan plane to ≈6 at the most distant plane.
Spin signal scan data and resulting 3D reconstruction of the hydrogen (proton) density distribution. (A) Raw scan data presented as x–y scans of the spin signal at 4 different tip-sample spacings. Pixel spacing is 8.3 nm × 16.6 nm in x × y, respectively. Each data point represents the spin signal variance obtained during a 1-min integration. Bext = 2.482 T. (B) A more finely sampled line scan showing 4-nm lateral resolution. The scanned region is indicated by the dashed line in A. Bext = 2.432 T. (C) Reconstructed 3D 1H spin density. Black represents very low or zero density of hydrogen, whereas white is high hydrogen density. The image is the result of the Landweber reconstruction, followed by a 5-nm smoothing filter. (D) Horizontal slice of C, showing several TMV fragments. (E) Scanning electron micrograph of the same region. (F) Cross-section showing 2 TMV particles on top of a hydrogen-rich background layer adsorbed on the Au surface. (G) Reconstruction is improved if this background layer is included as a priori information by assuming a thin, uniform plane of 1H density as the starting point of the reconstruction.
The correlation time τm of the spin signal was estimated from the signal bandwidth and found to be quite short, on the order of 20 ms. This is significantly shorter than has been seen in some inorganic test samples (26), possibly because of a naturally short rotating frame relaxation time T1ρ.
An indication of the spatial resolution inherent in the raw data can be seen in the line scan shown in Fig. 3B. The line scan, taken with d = 24 nm tip-sample spacing, shows a background spin signal of ≈300 aN2, with signal peaks up to ≈600 aN2. The leading edges of the peaks are sharp and indicate that the resolution in the x direction (before deconvolution) is on the order of 4 nm, which is approximately the thickness of the resonant slice (Fig. 2B).
Perhaps the most striking aspects of Fig. 3A are the apparent “double” images in each of the x–y scans. This feature doubling effect is due to the 2-lobed character of the PSF (Fig. 2B). As expected, the feature pairs are seen to move apart as the tip-sample spacing decreases.
The signal information from the 4 scan planes (together with a fifth “zero” plane at d = 75 nm, where no signal was found to be present) constitutes a 3D dataset that can be deconvolved to find the real-space image of the spin density. The result of the deconvolution is displayed in Fig. 3C, which shows the calculated spin density of the sample as a series of x–y planes at various z positions. A representative x–y plane is highlighted in Fig. 3D and clearly shows an assembly of a few individual virus particles (mostly fragments). Given that the raw MRFM data are spatially undersampled and have only modest SNR, the quality of the reconstruction is remarkable. The observation of significant improvement in image SNR after reconstruction is expected because most spins contribute force signal to more than one position in the scan, and the cumulative effect benefits the SNR of the reconstruction. The resolution appears to be in the 4- to 10-nm range, depending on the direction, with the x direction having the highest resolution. This resolution anisotropy is expected because of the directional dependence of the PSF, which reflects the fact that the cantilever responds only to the x component of magnetic force.
The fidelity of the MRFM reconstruction is confirmed by comparing the results in Fig. 3D to the scanning electron microscope (SEM) image of the same sample region (Fig. 3E). Excellent agreement is found even down to small details. Note that the origin of contrast in the 2 images is qualitatively very different. The MRFM reconstruction is elementally specific and shows the 3D distribution of hydrogen in the sample. Contrast in the SEM image is mainly due to the virus blocking the secondary electrons emitted from the underlying gold-coated cantilever surface.
The depth resolution that is enabled by the 3D image reconstruction is illustrated in Fig. 3F, which shows a vertical cross-section that extends through 2 virus fragments. Interestingly, the cross-sectional cut reveals that the virus particles sit on top of a thin layer of proton-containing material that is not visible by SEM. The presence of this background layer is not surprising because it is well known that surfaces that have been exposed to ordinary laboratory air become coated with a thin layer of adsorbed hydrocarbons and water. It is this proton-rich layer that is responsible for the 300 aN2 baseline signal in Fig. 3B. Taking this baseline signal level, combined with knowledge of the point-spread function and the typical spin density of hydrocarbons, one can estimate the layer thickness. Assuming a hydrocarbon spin density of ρ = 7 × 1028 spins per m3, the layer thickness is on the order of a nanometer.
The reconstructions shown in Fig. 3 C, D, and F use no a priori knowledge of the imaged object, such as the presence of the adsorbate layer on the sample. We find that the inclusion of this information into our initial guess of the spin density ρ0 substantially improves the result of the reconstruction. For Fig. 3G, we initialize ρ0 with a uniform layer of material positioned at the level of the cantilever surface, just below the virus particle. As seen in the figure, this modification allows for a much better separation of the background layer and more clearly resolves the TMV particles. Such constraints on ρ0 (or ρ) offer a general way to systematically include known structural information about the object.
A second MRFM image reconstruction from a neighboring sample area is shown in Fig. 4. Again, a clear 3D view of the hydrogen distribution is seen, including good spatial discrimination between the virus particles and the underlying adsorbed layer.
Imaging results for a second sample region. (A) Three-dimensional reconstruction of 1H spin density for virus particles sitting on adsorbed layer of hydrocarbons. (B) Representative horizontal slice from the 3D reconstruction showing the distribution of hydrogen in the plane located 13 nm above the hydrocarbon layer. Several virus particles are evident. (C) Corresponding scanning electron micrograph.
Discussion
We believe the present work represents substantial progress in developing the capability to probe the 3D chemical/elemental composition of nanostructures. With further progress in resolution and sample preparation, force-detected MRI techniques could have significant impact on the imaging of nanoscale biological structures, even down to the scale of individual molecules. Achieving resolution <1 nm seems realistic because the current apparatus operates almost a factor of 10 away from the best demonstrated force sensitivities (37) and field gradients (38).
Even with a resolution >1 nm, MRFM may allow the basic structure of large molecular assemblies to be elucidated. One can imagine enhancing MRFM image contrast beyond the basic spin-density information by using techniques similar to those developed for clinical MRI and NMR spectroscopy. Such contrast may include selective isotopic labeling (for example, substituting 1H with 2H), selective imaging of different chemical species (like 13C, 15N, or 31P), relaxation-weighted imaging, and spectroscopic imaging that reflects the local chemical environment (20, 39). Some techniques, such as cross-polarization and depolarization between different nuclear spin species, have already been demonstrated for MRFM on the micrometer scale (21, 22). At the nanometer scale, the ability to target and locate specific proteins although selective labeling, for example, could allow direct 3D imaging of the organization and structure of multicomponent macromolecular complexes. Such a capability would be complementary to current techniques, such as cryoelectron microscopy, and could develop into a powerful tool for structural biology.
Acknowledgments
We thank C. Micheel, G. Wallraff, J. Cha, and S. Swanson for help with sample preparation; J. Frommer for AFM support; M. Hart and M. Farinelli for fabrication of the microwire and magnetic tip; M. Ting and A. Hero for discussions on image reconstruction; and B. Chui for cantilever fabrication. This work was supported by IBM and the National Science Foundation-funded Center for Probing the Nanoscale (National Science Foundation Grant Number PHY-0425897). C.L.D. acknowledges support from the Swiss National Science Foundation.
Footnotes
- 1To whom correspondence should be addressed. E-mail: rugar{at}almaden.ibm.com
Author contributions: C.L.D., M.P., H.J.M., and D.R. designed research; C.L.D., M.P., H.J.M., and D.R. performed research; C.T.R. contributed new reagents/analytic tools; C.L.D., M.P., H.J.M., and D.R. analyzed data; and C.L.D., M.P., and D.R. wrote the paper.
The authors declare no conflict of interest.
This article contains supporting information online at www.pnas.org/cgi/content/full/0812068106/DCSupplemental.
- Received August 21, 2008.
- © 2009 by The National Academy of Sciences of the USA
Freely available online through the PNAS open access option.
References
- ↵
- ↵
- Buxton RB
- ↵
- Callaghan PT
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- Rugar D,
- et al.
- ↵
- Zhang Z,
- Hammel PC,
- Wigen PE
- ↵
- ↵
- ↵
- ↵
- ↵
- Garner SR,
- et al.
- ↵
- Mamin HJ,
- Budakian R,
- Chui BW,
- Rugar D
- ↵
- ↵
- ↵
- ↵
- ↵
- Kuehn S,
- Hickman SA,
- Marohn JA
- ↵
- ↵
- Poggio M,
- Degen CL,
- Rettner CT,
- Mamin HJ,
- Rugar D
- ↵
- ↵
- Klug A
- ↵
- Taylor KA,
- Glaeser RM
- ↵
- ↵
- ↵
- ↵
- Müller N,
- Jerschow A
- ↵
- Bertero M,
- Boccacci P
- ↵
- ↵
- Ting M,
- Raich R,
- Hero AO
- ↵
- Mamin HJ,
- Rugar D
- ↵
- ↵
- Haacke EM,
- Brown RW,
- Thompson MR,
- Venkatesan R
Citation Manager Formats
More Articles of This Classification
Physical Sciences
Applied Physical Sciences
Biological Sciences
Microbiology
Related Content
- No related articles found.
Cited by...
- Improving mechanical sensor performance through larger damping
- Nanoscale nuclear magnetic resonance with chemical resolution
- Nonlocally sensing the magnetic states of nanoscale antiferromagnets with an atomic spin sensor
- Nuclear magnetic resonance detection and spectroscopy of single proteins using quantum logic
- Single proteins under a diamond spotlight
- Single-proton spin detection by diamond magnetometry
- Detection of atomic spin labels in a lipid bilayer using a single-spin nanodiamond probe
- Toward Molecular-Scale MRI
- Nuclear Magnetic Resonance Spectroscopy on a (5-Nanometer)3 Sample Volume
- Nanoscale Nuclear Magnetic Resonance with a Nitrogen-Vacancy Spin Sensor
- A nanoscale experiment measuring gravity's role in breaking the unitarity of quantum dynamics
- Scanned-probe detection of electron spin resonance from a nitroxide spin probe


















