Main

For many decades, diffraction precluded far-field optical imaging with resolution better than half the wavelength of light. This limit has been overcome with the advent of STED microscopy, which is a lens-based fluorescence microscope with a resolution conceptually no longer limited by diffraction13,16,17. In brief, a STED microscope separates two or more adjacent objects by ensuring that only one of them fluoresces at a time. Sequential object signalling is accomplished by transiently switching the objects off by stimulated emission.

STED microscopy has so far provided a resolution of 16–80 nm, which has been used for mapping proteins in cells18, and for the real-time imaging of colloidal particles or synaptic vesicles in living neurons19. The fluorescent markers used were organic molecules or fluorescent proteins having a singlet ground state (S0) and a first excited state (S1), each with vibrational sublevels. Once excited, these molecules equilibrate to the lowest vibrational sublevel of the S1 state, from where they can return to a vibrational sublevel of the S0 state by fluorescence emission after an average period τfl. Alternatively, the so-called STED beam switches the molecules off by forcing them down to the S0 state. This occurs if the stimulated emission rate σ ISTEDλSTED/hc is much larger than that of the spontaneous decay, 1/τfl, where σ is the cross-section for stimulated emission, h is Planck's constant, c is the speed of light in vacuum, and ISTED is the STED beam intensity. The wavelength λSTED is chosen such that the fluorophore is pumped into a rapidly decaying higher vibrational sublevel of the S0 state, which makes de-excitation predominant. If we neglect polarization, the ability of the molecule to fluoresce decreases in proportion to exp(− ISTED/IS) where IS is a fluorophore characteristic17. Hence, for ISTED ≥ 3IS the fluorophore is essentially switched off.

Organic fluorophores follow such switching, but have the disadvantage of suffering from bleaching and blinking, making solid-state-based fluorophores very appealing. It is also interesting, however, to explore applications of far-field optical nanoscopy to transparent solids. Luminescent transitions arising from defects in diamond, particularly nitrogen-vacancy centres composed of a substitution nitrogen next to a charged vacancy (Fig. 1a), have attracted much attention2,20. Absorption at 532 nm elevates them from their ground (3A) state to an excited (3E) state of τfl = 11.6 ns lifetime, which emits at 600–850 nm with virtually no bleaching or blinking. Phononic coupling yields a pronounced Stokes shift, but the involved states are now spin triplets. Centres with magnetic quantum number ms = 0 emit more strongly, because their ms = ±1 counterparts have an increased tendency to convert into a dark singlet state (1E), ultimately yielding ground-state centres with ms = 0, a feature that allows for the orientation and detection of spin by illumination2,20. When used as scanning tips, nitrogen-vacancy centres in nanosized diamonds provide optical21 and magnetic imaging9,10, but the prospect of reading and controlling individual defects at the nanoscale with conventional lenses is very attractive.

Figure 1: Stimulated emission depletion of the excited state ‘switches off’ a charged nitrogen-vacancy colour centre by light.
figure 1

a, Energy diagram sketching the ground (3A), fluorescent (3E) and the dark singlet state (1E) along with excitation (Exc), emission (Em) and stimulated emission (STED) transitions in nitrogen-vacancy centres. b, The steep (‘rectangular’) fluorescence depletion curve showing the abrupt decrease of the population of 3E with intensity ISTED facilitates a nearly perfect implementation of STED. The inset representation confirms the exponential suppression of the excited state. For ISTED  > 20 MW cm−2, the centre is deprived of its ability to fluoresce (switched off).

We therefore investigated STED on nitrogen-vacancy defects within diamonds of type IIa grown by chemical vapour deposition. Excitation at 532 nm was achieved using a 60-ps pulsed laser diode, which was triggered by the pulses of the STED beam provided by a laser system operating at λSTED = 775 ± 1 nm and emitting 3.2 ns pulses at 8 MHz. The synchronized pulses were combined using dichroic mirrors and collinearly directed through a 1.46 numerical aperture (NA) oil immersion lens, which also imaged the fluorescence onto the confocally arranged aperture of a counting avalanche photodiode.

The phononic coupling in the 3A state enables quenching of excited nitrogen-vacancy centres by stimulated emission. This is demonstrated in Fig. 1b, which shows the fluorescence of a single centre as a function of ISTED, here defined as the peak intensity in the focal plane. The time-averaged power of the excitation light was 40 µW. A steep exponential drop of the fluorescence with ISTED is observed, as one would anticipate for perfect STED (ref. 13). Extracting IS = 6.6 MW cm−2 yields σ = 1.2 × 10−17 cm2, which agrees with the value determined from the emission spectrum. The main outcome, however, is that the nitrogen-vacancy centres can be sharply switched off by light (Fig. 1b).

To make use of this nearly ‘rectangular’ depletion curve13 for nanoscopy, we produced a doughnut-shaped focal intensity distribution22  (x,y) with a central minimum (0,0) ≈ 0. The minimum is ideally zero, but in practice it amounts to a fraction ɛ 1 of the maximum intensity at the doughnut crest:  = ɛ max[((x,y)]. Overlapping the Airy disk of the excitation beam with the doughnut and enforcing IS switched off all the centres covered by the Airy disk except for those at the doughnut minimum where (x,y) ≤ 3IS. Therefore, the Airy disk can be ignored when calculating the spot in which the centres can still be ‘on’, that is, the effective point-spread function (PSF) of the system. It is described by the function exp(−(x,y)/IS) having a full-width at half-maximum (FWHM) of Δx,y ≈ λSTED /√(C/IS), with the constant C depending on the shape of the minimum22. Although the resolution Δ x,y scales with the wavelength, adjusting IS squeezes the PSF continuously, indicating that the wavelength has surrendered its limiting role.

We imaged by scanning the sample with respect to the beams using a three-dimensional piezo stage (Physik Instrumente, P-733.3DD) with 0.1 nm positional accuracy and 0.03% linearity. Figure 2 compares a confocal image (Fig. 2a) with the corresponding STED image (Fig. 2b), taken at = 1.4 GW cm−2. Although the confocal reference lacks any detail, its STED counterpart resolves each nitrogen-vacancy defect. The similar brightness and sparse occurrence of the dots shown in the crystal allow us to conclude that they represent single nitrogen-vacancy defects. The images consist of 75 × 65 pixels of 4 × 4 nm in size, recorded with a dwell time of 5 ms. To reduce statistical error (number of counts given in the colour table), the displayed and analysed STED image comprises the sum of 97 consecutively recorded images. Representing the effective PSF, the images of individual centres reveal Δx = 18.3 nm and Δy = 16.0 nm as the resolution under the conditions used. Note that the resolution gain derives just from the optical preparation of the states.

Figure 2: Stimulated emission depletion microscopy reveals densely packed nitrogen-vacancy centres in diamond.
figure 2

a,b, Confocal (a) and STED (b) images from the same crystal region. c, The individual centres resolved in b automatically yield the effective PSF of the STED recording whose y-profile exhibits a FWHM of Δy = 16.0 nm. d, The coordinate of each centre can be calculated with 0.14 nm precision. Comparing a with d highlights the dramatic gain in information resulting from the unique increase in resolution. e,f, Applying  = 3.7 GW cm−2 shrinks a confocal spot of 223 nm diameter (FWHM) down to 8 nm. Note that the increase in resolution is a purely physical phenomenon.

The STED recording inherently resolves individual nitrogen-vacancy centres because they are all further apart from each other than Δx or Δy. An important insight is that, when scanned over the sample together with the excitation beam, the doughnut-shaped STED beam allows only a single nitrogen-vacancy centre to fall within the doughnut minimum; all others are switched off or simply not excited. Thus STED nitrogen-vacancy centres are resolved individually and hence registered sequentially in time.

Once all the centres are resolved, increasing the resolution further is no longer required. Their position can be calculated from the centroid of their subdiffraction-sized spot23. The precision of this calculation depends on the effective PSF, the number of collected photons and the background24,25. We measured a spatially homogeneous background of 1.6 counts ms−1 throughout the imaged region, as indicated in the images. Under these conditions and the assumption (here not made or asserted) that scanning errors are negligible, the xy-coordinates of each centre can be derived with 1–2 Å precision; that is, by three orders of magnitude below the diffraction limit (Fig. 2d). The information provided in Fig. 2d should be contrasted with that contained in the confocal recording (Fig. 2a), which fails to render any information about the number and location of the colour centres. The distance between the centres in the upper left corner projected in the xy-plane is 67.9 ± 0.1 nm. By extending STED into the third dimension, one can also incorporate depth information. Thus, provided that care is exerted to avoid scanning distortions and drift, STED microscopy is able to map nitrogen-vacancy centres with nanometric resolution and subnanometre precision.

Figure 2e shows a confocal image of an isolated nitrogen-vacancy centre, but near its middle, the STED beam is momentarily engaged. Applying  = 3.7 GW cm−2 produced a spot of 8 nm FWHM. In this case, the peak is also reduced by 3.5-fold, which is due to the nonzero doughnut minimum (0,0) = ɛ suppressing the peak by exp(–ɛ/IS). Clearly, such peak reductions become important at large and therefore at high resolution. Evaluation of the measurement using Fig. 1b yields ɛ = 7 × 10−4, which is due to doughnut imperfections. Note that the effective focal area decreases by a factor of 777.

Importantly, the outstanding photostability of nitrogen-vacancy centres1 is maintained under the intense STED beams, as made evident in consecutive STED recordings (Fig. 3) in which both the brightness and the tenfold resolution increase remain unchanged. Demonstrating optical nanoscopy without photobleaching, Fig. 3 highlights that with diamond nanoparticles as markers11,12, STED will enable the use of (biological) fluorescence nanoscopy with countless recordings.

Figure 3: STED microscopy without photobleaching.
figure 3

a, Consecutive confocal and STED recordings (as indicated) of the same crystal region. b, The total fluorescence signal in each image frame, showing virtually no change in brightness and in resolution (here <25 nm) with increasing number of recordings. The recording time for all images was the same, causing the STED images to be darker approximately by the (here 100-fold) reduction in focal area. Importantly, the peak brightness of single centres is virtually equal in the STED and confocal recordings. Note the absence of blinking. The experiment demonstrates far-field optical nanoscopy recording without photobleaching.

Recording images similar to those in Figs 2 and  3 is impractical with near-field optics, because, being located more than 2–3 µm beneath the surface of the crystal, the colour centres cannot be accessed by a tip. Besides being noninvasive and able to access the interior of the crystal, compared to near-field optics, STED microscopy is easier to put into practice, control and operate, and provides substantially higher resolution. Likewise, alternative nanoscopy strategies that sparsely switch individual emitters (from a dark state to a state emitting a bunch of photons allowing their localization26,27,28,29) are hard to implement for nitrogen-vacancy colour centres, because at room temperature, blinking or metastable switching between a bright and a dark state is nearly absent. Hence, nitrogen-vacancy centres are ideal for STED and STED is ideal for their imaging. This is likely to apply to many other colour centres in crystals.

Next we used a focal STED beam spot with a y-oriented central zero-line with steep peaks along the x-axis16. Increasing to 8.6 GW cm−2 squeezes the FWHM of the PSF from Δx = 223 nm down to 6 nm (Fig. 4a). The increase in resolution follows the anticipated inverse square-root law15,16. An exemplary super-resolved PSF of a single nitrogen-vacancy centre is shown in Fig. 4b and has a measured Δx = 5.8 ± 0.8 nm. This value is 133-fold below the responsible wavelength and represents a new regime of obtainable all-physics-based resolution using conventional optics.

Figure 4: Maximizing the resolving power by optimizing STED in a single direction.
figure 4

a, The measured decrease of the lateral FWHM Δx from 223 nm down to 6 nm scales inversely with the square-root of the intensity of the STED beam (red line). Note that all data were taken from the same single nitrogen-vacancy defect. b, A line profile of the effective PSF along the x-axis, demonstrating an unprecedented all-physics-based far-field optical resolving power of 5.8 ± 0.8 nm.

Despite the leaps and bounds STED microscopy has made, its limits have not yet been reached. Applying a shorter λSTED entailing a larger σ should decrease Δxy even further. Implemented with a deep zero (ɛ 1 × 10−4), the initially desired nearly ‘rectangular’ depletion curve13 should enable a further ‘spot sharpening’ by more than a factor of 10, which may surpass the dimensions of the electron cloud and, leaving the dipole approximation behind, open up the fascinating prospect of sensing the (size of the) cloud and its immediate environment. It would also be desirable to investigate STED on nitrogen-vacancy centres subject to magnetic and microwave fields. In particular, reading out nitrogen-vacancy centres and hence local magnetic fields on the nanoscale should be easier with STED and conventional lenses than with nanosized optical tips. Thus, STED microscopy may greatly simplify magnetic imaging at the nanoscale9,10. Likewise it will be interesting to apply STED nanoscopy to other colour centres and crystals. Finally, the achieved and future gains in all-physics-based resolution should stimulate novel explorations in applied and fundamental (quantum) optics.

Methods

The cylindrical doughnut focal spot used for generating the data of Figs 2 and  3 was produced by inserting a glass waveplate (RPC Photonics), inducing an azimuthal linear helical phase ramp from 0 to 2π on the initially flat wavefront of the STED beam. In these experiments both the excitation and the STED beam were circularly polarized. The bandwidth of the STED beam was <10 GHz.

The data of Fig. 4 were produced with linearly polarized excitation and STED beams. A phase shift of π, resulting in the y-oriented line-shaped valley, was produced by inserting a home-built phase-shifter consisting of tilted optical flats. In this case, the polarization of both beams was oriented along the direction of the line-shaped valley. Images in Fig. 2a,b were taken from a diamond sample with a higher density of nitrogen-vacancy centres than the crystal used for the measurements shown in Figs 2e and  4a,b.

The intensity values refer to spatial and temporal peak intensities assuming aberration-free focusing and rectangular laser pulses, with the transmission of the objective lens 70% at 775 nm. A 95% transmission was assumed at the interface between the immersion oil and the diamond crystal. In all recordings, background noise was reduced to 120 counts s−1 by gated detecting within a 40-ns window after the STED pulse.

The stimulated emission cross-section was calculated using σ ≡ hc/(λSTEDISτp) under the approximation that the change in population of the excited state of the nitrogen-vacancy centre was negligible during the duration τp of the STED pulse. The cross-section was also derived from the emission spectrum E(λ) of a single defect normalized to the quantum yield3,30 (Q = 0.7 = ∫E(λ) dλ) using σ = λ4E(λ)/8π cn2τfl, with the speed of light in vacuum c and the refractive index of diamond n = 2.417.

The STED image of Fig. 2b was generated by summing 97 STED images of the same region. Alignment was ensured by introducing subpixel translational shifts to maximize the spatial correlation between the images. This corrects exclusively translational drift between consecutive images. The central positions of the defects in the registered images were stochastically distributed with a variance close to the expected localization accuracy, supporting the notion that no further substantial corrections were required in this case. The centroids were determined by Gaussian fits.

Data points within Fig. 4a were fitted with the general expression22 Δx ≈ xc /√(1 + xc2a2 /(IS ln 2)), where xc is the FWHM for  = 0 and a is the steepness of the STED intensity pattern ISTED(x) ≡ 4a2x2, giving a = 1.76 λ/NA. The measured STED PSF independently yielded a steepness a = 1.75 λ/NA. The line profile of Fig. 4b represents eight consecutive x-scans over the same nitrogen-vacancy defect, revealing an average measured FWHM of the scans of 5.8 nm, with a standard deviation of 0.8 nm.