Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Analytical description of STED microscopy performance

Open Access Open Access

Abstract

Stimulated emission depletion (STED) resolves fluorescent features that are closer than the far-field optical diffraction limit by applying a spatially modulated light field keeping all but one of these features dark consecutively. For estimating the efficiency of transient fluorophore darkening, we developed analytical equations considering the spatio-temporal intensity profile of the STED beam. These equations provide a quick analysis and optimization of the resolution and contrast to be gained under various conditions, such as continuous wave or pulsed STED beams having different pulse durations. Particular emphasis is placed on fluorescence fluctuation methods such as correlation spectroscopy (FCS) using STED.

©2010 Optical Society of America

1. Introduction

Stimulated emission depletion (STED) overcomes the diffraction limit in fluorescence microscopy by ensuring that a (large) part of fluorophores in the diffraction-limited focal region is non-fluorescent (dark) for a brief period of time in which the remaining fluorophores are detected. As a result, a subset of fluorescent (bright) fluorophores is separated from a subset of dark neighbors within the diffraction region. Transient fluorescence deactivation is attained by a beam of light, called the STED beam, inducing stimulated emission from the fluorescent state of the dye [13]. Fluorophores that are exposed to an intensity I STED of the STED beam that is well beyond a threshold I s (at which the rate of stimulated emission equals the spontaneous decay) spend virtually all time in the ground state even when illuminated by excitation light; their capability to emit fluorescence is ‘switched off’.

In the focal plane of a typical STED microscope an excitation light spot is overlaid with a STED light spot featuring a central intensity minimum (“zero”). In regions with high STED intensity I STED the spontaneous decay of fluorophores is largely suppressed, but this is clearly not the case for those fluorophores residing at or close to the central minimum. The more intense the STED beam, the narrower the region becomes in which the fluorophores are allowed to emit. For most fluorophores the lifetime of the fluorescent state is in the order of nanoseconds and the cross section for stimulated emission σ STED is about 10−17–10−16cm2. Hence a STED beam providing a time-averaged focal intensity ISTED>>Is~5MW/cm2 turns the fluorophores off even when they are flooded by excitation photons.

To keep the fluorophores dark using a moderate average laser power, most STED microscopes (‘nanoscopes’) employ a train of light pulses for both excitation and STED [13]. Numerical calculations of the STED efficiency have been performed elsewhere [1,2,4,5]; in the following, we will show how to optimize the performance of STED microscopy by tuning the pulse width of the STED beam. For this purpose, we calculate the fluorescence deactivation (‘off-switching’) efficiency by using approximations suitable for obtaining an analytical expression for the STED nanoscopy performance.

2. STED nanoscopy

STED efficiency

Figure 1 shows the Jablonski diagram of a typical organic fluorophore with the electronic ground singlet state S0, the electronic excited singlet state S1, the lowest triplet or other dark state T1, and vibrational sub-levels as well as the most relevant transitions between these states. After excitation from S0 and vibrational relaxation, the fluorophore is found in the base level of the excited state S1, from where it spontaneously decays to S0 or crosses to T1. If spontaneous decay takes place, a photon is emitted with a fluorescence quantum yield qfl1. On the other hand, the excited fluorophore can be sent down to a higher vibrational sub-level S0* with a red-shifted light field inducing stimulated emission. If the stimulated emission rate is larger than the spontaneous decay, the fluorophore remains dark even under excitation irradiation; the residence in the excited state S1 is essentially disallowed by the presence of the STED beam. Interruption of the STED beam turns the dye on instantly, because S1 is allowed again. Keeping alternating features dark for a brief period of time enables the separation of features that are mutually closer than the diffraction limit. It is essential to realize that in the concept of STED microscopy, the sole purpose of the STED beam is to ensure that not all but just one of the features located within the diffraction region emits signal.

 figure: Fig. 1

Fig. 1 Simplified Jablonski diagram with electronic and vibrational energy levels of a typical fluorophore.

Download Full Size | PDF

For estimating the fluorescence ability of a fluorophore under STED, we assume that the fluorophore is initially in S1 due to a short excitation pulse. Furthermore, we assume a rectangular STED pulse lasting for τ STED, and a stimulated emission rate kSTED=ζkS1, with k S1 denoting the spontaneous decay rate and ζ=ISTED/Is being the saturation factor. ζ scales with the ratio of the pulse intensity I STED of the STED light and the saturation (threshold) intensity I s [6]. The saturation intensity here defines the STED intensity at which the rate kSTED=σSTEDΦSTED of stimulated emission equals the spontaneous decay rate k S1 (i.e., ζ = 1). We define ΦSTED=ISTEDλSTED/(hc) as the photon flux, λ STED is the STED wavelength in free space and hc=1.991025Jm. Consequently, Is=kS1hc/(λSTEDσSTED) is a molecule-specific constant given by the cross section σ STED of stimulated emission and k S1. During the STED pulse, i.e. for 0 < t < τ STED, the relevant rate equations for the populations P S1 of S1 and PS0* of S0*, respectively, are

tPS1=kS1PS1kSTED(PS1PS0*)tPS0*=kvibPS0*+kSTED(PS1PS0*)
where k vib is the vibrational relaxation rate of S0* to S0. We normalize the populations to one molecule, i.e. PS1(0)=1; all other states are initially unpopulated. Because kvib>>kS1 and as long as τSTED>>1/kvib, the equilibration of the population PS0* is fast, such that the excited state population can be well approximated by PS1(0<t<τSTED)exp(kS1t(1+γ)), where γ=ζkvib/(ζkS1+kvib) is the effective saturation factor that also accounts for the generally undesired excitation from S0* to S1 by the STED beam itself. Thereafter, the excited state population spontaneously decays at a rate k S1 and is simply given by PS1(τSTED<t<T)exp(kS1(γτSTED+t)), where T is the period between two STED pulse peaks in a pulse train, i.e. the inverse of the pulse repetition rate. The fluorescence probability F per molecule is then given by

F(γ)=qflkS10TPS1(t)dt=qfl(1+γexp(kS1τSTED(1+γ))1+γexp(kS1(γτSTED+T))).

The probability of spontaneous decay η ps is now obtained as F(γ)/F(0), which equals the fraction to which the fluorescence is suppressed by the STED pulse of intensity I STED (which is proportional to ζ and γ), and where F(0)=qfl(1exp(kS1T)) is the fluorescence probability of the undisturbed molecule.

ηps=(1+γexp(kS1τSTED(1+γ))1+γexp(kS1(γτSTED+T)))(1exp(kS1T))1
This is the main equation describing STED analytically.

Here, k S1 τ STED and γkS1τSTEDkSTEDτSTED=ζkS1τSTED=ISTEDσSTEDτSTED characterize the average number of spontaneous and the average number of possible stimulated de-excitation events during the STED pulse τ STED, respectively [7]. k STED τ STED is a measure of the strength of the applied STED irradiation. For a constant pulse repetition rate 1/T, it scales linearly with the average STED power PSTED.

If the period T is long compared to the excited state lifetime, k S1 T is large and Eq. (3) simplifies to the first term, i.e.

ηps1+γexp(kS1τSTED(1+γ))1+γ.

For large kS1τSTED(1+γ), Eq. (4) becomes then

ηps1/(1+γ).

This is the case for long STED pulses, i.e. for the case of applying pulsed excitation in conjunction with continuous wave (CW) STED (τ STED = T). The equation also holds for large k S1, i.e. for fluorophores with short fluorescence lifetimes, in which case I s is large and the I STED necessary to achieve a certain fluorescence inhibition becomes very large as well. Likewise Eq. (5) also applies if large I STED, i.e. large ζ and γ are used. This condition also defines the maximum achievable fluorescence inhibition ηps(kSTED)=1/(1+kvib/kS1), which is determined by the finite vibrational relaxation rate k vib. For typical fluorophores and measurement conditions, the ratio k vib/k S1 ~103–104 (k S1 ~1/ns compared to k vib > 1/ps), hence η ps < 1‰ approaches almost complete fluorescence suppression. An infinite relaxation rate results in ηps1/(1+ζ)=1/(1+ISTED/Is)0 for ISTED. In most experiments however, ηps(ISTED)>0 because of a small though non-zero probability of S0 → S1 excitation by the STED light.

Equation (3) is applicable for pulsed excitation; STED can be either pulsed or CW. In case of CW excitation and CW STED, the ongoing excitation continuously works against the inhibition of the fluorescent state occupation, i.e. against the off-switching of the fluorophore. Neglecting the triplet state, the probability of spontaneous decay η CW writes as

ηCWkex+kS1kex+(1+γ)kS1.

This is the main equation for CW STED, which derives from the ratio of the average time periods a fluorophore requires for a spontaneous decay from S1 with and without STED. Without STED, this is the time for a full excitation–de-excitation cycle kex1+kS1-1. With STED, more excitations by a factor of γ are required, that is a time of (1+γ)kex1+kS1-1 per spontaneous decay. For low excitation avoiding saturation of the S0 → S1 transition (k ex << k S1), Eq. (6) simplifies to the expression of Eq. (5). In this case of CW STED the same fluorescence inhibition is achieved regardless of applying pulsed or CW excitation. More detailed calculations taking into account the pulse delay as well as the excitation power and triplet state dependent response are summarized in the annexes.

Previous work has given different descriptions of the fluorescence inhibition efficiency η ps by STED [6]. Assuming kvib and regarding only the residual fluorescence emitted after the STED pulse, η ps was calculated from the single rate equation PS1/t=kS1PS1kSTEDPS1 as the remaining population of S1 after the pulse duration τ STED.

ηps*=exp(kS1τSTEDγ)=exp(ln(2)ISTED/Is*)
The saturation intensity Is*=ln(2)hc/(λSTEDσSTEDτSTED) had been introduced as the STED intensity to inhibit half of the S1 population, i.e. half of the spontaneously emitted fluorescence (ηps*(Is*)=0.5) after the STED pulse. Equation (3) differs mainly in the aspects that it accounts for k vib and the pulse period T and that it includes the spontaneous fluorescence emitted during the STED pulse.

Pulse length dependence

Figure 2 outlines the estimated STED performance, i.e. the probability of spontaneous decay η ps versus ζ (i.e. pulse STED intensity) calculated for k S1 = 1/3.9ns, k vib = 5/ps and 80MHz pulse rate. These values are taken from the organic dye Atto647N in a supported lipid bilayer [8] and a mode-locked Ti:Sapphire laser. In Fig. 2a, the thin line shows the lower limit 1/(1+ζ), which is achieved for pulsed excitation but continuous wave (CW) STED and infinite k vib. The bold lines indicate the performance for pulsed STED (Eq. (3)) showing an ever more pronounced knee, the shorter the STED pulses become. This is because the STED intensity has to counterbalance the short excited state depletion time before the suppression becomes efficient. Regardless of the pulse width, the suppression efficiency saturates at a lower bound (here η ps ≈5·10−5) due to the vibrational relaxation rate k vib, which sets the upper limit for the effective stimulated depopulation rate of S1 (Eq. (5)). In addition, the thin dotted lines of Fig. 2a and the lines of Fig. 2b trace η ps versus ζ and τ STED, respectively, for constant average STED powers that potentially stimulate the S1 de-excitation kSTEDτSTED=ζkS1τSTED equals 1–16 times per pulse.

 figure: Fig. 2

Fig. 2 Probability of spontaneous decay η ps as a function of the saturation factor ζ and the pulse width τ STED of the STED light (as labeled) for k S1 = 1/3.9ns, k vib = 5/ps and T = 1/80MHz. (a) The thin line shows the lower limit achieved for CW STED and infinitely fast vibrational relaxation. Thin dashed lines show the dependence of η ps on variable pulse lengths τ STED for fixed numbers of potential de-excitations k STED τ STED, i.e. for fixed average STED powers. (b) This dependence is shown as a function of τ STED. The inset outlines the pulse scheme used in the calculations: a Dirac excitation pulse followed by a rectangular STED pulse.

Download Full Size | PDF

Increasing ζ or τ STED along the k STED τ STED isolines entails an according decrease of τ STED and ζ, respectively. For very long pulses, i.e. small ζ, the suppression of the fluorescent state is quite inefficient because spontaneous emission is as fast or faster. For very short pulses on the other hand, i.e. large ζ, the suppression becomes inefficient again because of the limited S1 de-excitation rate stemming from the finite vibrational relaxation rate k vib. Figure 2 reveals that a pulse width of τ STED ≈20–30ps would exploit the available STED power maximally and would yield the steepest decrease of η ps with increasing ζ. Both characteristics are preferable in STED imaging as outlined further on.

In many practical cases however, it turned out that a pulse width of τ STED ≈100–150 ps leads to better results. The following reasons can be cited. Firstly, longer pulses leave more time for the excited molecules to rotate, which helps minimizing polarization effects [9]. Secondly, if pulsed diode lasers are used for excitation [10], the excitation pulse width of about 50–70ps and the synchronization jitter can be better countered with longer STED pulses. Thirdly, it has been shown that longer pulses are more advantageous with respect to photobleaching of the fluorophores [11], since these phenomena usually scale with higher orders of the intensity [12]. Strongly related with the previous argument is the fact that longer pulses disfavor multiphoton excitation processes, including plain multiphoton excitation of relaxed molecules from the ground state [13].

Spatial resolution and STED imaging

In the focal plane of a typical STED microscope a diffraction-limited excitation light spot is overlaid with a STED light spot featuring a local intensity minimum at the center of the excitation spot. An example is the doughnut-shaped intensity pattern recalled in Fig. 3 [14]. Consequently, fluorescence is inhibited everywhere but at the focal center and the more intense the STED beam, the more fluorophores and fluorescent features are switched off by the STED light and the narrower the region becomes in which the fluorophores are allowed to fluoresce. Figure 4 exemplifies calculated cross-sections through effective STED point-spread-functions (PSF) of such a confocal epi-fluorescence microscope with an oil immersion objective of numerical aperture (NA) of 1.40 and a fluorescent dye with a lifetime of 3.9ns [15,16]. The effective STED PSF describes the region in which the dyes are allowed to emit and subsequently be registered by the detector. More precisely, it gives the normalized probability with which the fluorophore is able to contribute to the measured signal at a given location in the sample. The excitation wavelength λ ex is 635nm and the STED wavelength λ STED is 770nm. Both beams overfill the objective aperture and are circularly polarized. The STED beam passes through a 2π phase ramp (vortex plate) for creating a doughnut with a central “zero”, i.e. a minimum with a relative intensity of ε times the intensity at the doughnut crest [6]. The detection pinhole has a projected diameter of 500nm in the sample space, i.e. 0.9 × the Airy disc diameter. At low excitation intensity, the detected fluorescence brightness profile is the product of the excitation probability times the emission probability η times the detection efficiency.

 figure: Fig. 3

Fig. 3 Cross-sections through the foci of the excitation and the STED light spots (a) and of the fluorescence collection efficiency (b) used in the calculation examples. Insets show the spot profiles in the xy-planes.

Download Full Size | PDF

 figure: Fig. 4

Fig. 4 Effective PSF of a STED nanoscope with (a) pulsed (τ STED = 50ps) and (b) CW STED (τ STED = T = 12.5ns): lateral cross-sections of the detected fluorescence brightness without STED (bold solid lines) and with average STED powers required for achieving FWHM diameters of 72nm (thin solid lines) and 33nm (pointed red lines). Dash-dotted lines show the profiles with an imperfect “zero” and dotted lines show them normalized again to their peak brightness. For values below 10% of the peak brightness, the logarithmic scale allows a better distinction of the differences.

Download Full Size | PDF

Without STED, the detected fluorescence (λ fl = 670nm) has a brightness profile with a full width at half-maximum (FWHM) diameter of 232nm. Applying STED pulses or CW STED irradiation, the FWHM diameter is narrowed with increasing STED power in all cases and follows a similar dependence regardless of the STED modality as summarized in Fig. 5 .

 figure: Fig. 5

Fig. 5 FWHM diameters versus STED power. For k STED τ STED > 500 in case of 50ps STED pulses with ε = 1%, the increase of the FWHM indicates loss of signal due to a too strong reduction of the central peak brightness versus the contributions from the focal periphery.

Download Full Size | PDF

For obtaining the same FWHM diameter, the CW STED power has to be increased about five fold (~1.5k S1 T) as compared to the average STED power in the pulsed case [17]. As predicted [18] and experimentally verified [9], the inverse of the FWHM diameters of the STED foci follow a square root law. In the pulsed mode, the FWHM is about FWHM0/1+kSTEDτSTED with the diffraction limited confocal FWHM0 ≈0.52λ ex/NA and where k STED τ STED is the number of potential stimulated emission events per molecule at the doughnut crest (if there were as many excitations). With 50ps long STED pulses and k STED τ STED = 50 at the crest, a FWHM of about 33nm is reached in the outlined example. In the vicinity of the “zero”, changing the pulse width τ STED does not significantly alter the profile as shown by the CW PSF, where the STED “pulse width” equals the excitation pulse period T = 12.5ns.

In Fig. 4, the dash-dotted lines show the effect of an imperfect “zero” with a residual intensity amounting to ε = 1% of the crest intensity. In the outlined examples, a vortex plate is used for creating the doughnut, which is rather tolerant to many experimental imperfections such as astigmatic or spherical wave front aberrations or lateral misalignments. Using a vortex plate, an imperfect zero typically stems from wrongly polarized light (c.f [19]. for examples) and/or from asymmetric wave front aberrations. For instance, elliptically polarized light with an electric field varying by ~20% peak-to-peak results in ε = 1%. In this case, the peak brightness drops significantly, in particular at high STED power, but the diameter is hardly affected for short STED pulses. However, with CW STED the resolution deteriorates (with respect to the case of ε = 0) to a value of FWHM = 43nm as shown by the normalized profiles. Note, that along with the reduction of the fluorescence at the central peak, suppression of fluorescence is not complete in the focal periphery. While the contribution of this peripheral fluorescence is in principle low, it is less for short STED pulses and for ε = 0, i.e. the effective STED PSF is sharper and offers a better contrast. The improved performance in the case of short STED pulses results from the steep dependence of η ps with increasing ζ (Fig. 2).

3. STED nanoscopy for fluorescence fluctuation spectroscopy

So far, we measured the resolution by the FWHM of the central peak. Unfortunately, this measure may become inappropriate for very high resolution, i.e. small FWHM values. Because the emission of peripheral fluorophores cannot be suppressed completely, their contribution to the total signal becomes increasingly important in this case (Fig. 4). As long as STED nanoscopy is applied for imaging comparatively sparse features, the effective PSF in real space, respectively its optical transfer function (OTF) in Fourier space, can be usually exploited for reducing undesired image artifacts like blur or ghost images. Within the limits of the signal-to-noise ratio (SNR) of the image(s), image deconvolution allows improving the image quality by re-attributing faint peripheral signals to their original features. However, such image restoration methods are hardly (if ever) applicable for spectroscopic measurements, such as fluorescence correlation spectroscopy (FCS) [20,21] for instance, essentially because the fluorescence signal from randomly distributed and/or mobile fluorophores is collected and the spatial information (neighborhood) has to be ignored in favor of sufficient time resolution, i.e. the assumption of sparse and spatio-temporally resolved features fails. For example, previous STED-FCS measurements of fluorophores diffusing in open measurement volumes of nanoscale dimensions showed a decrease in signal-to-background ratio (SBR) with increasing resolution due to a growing contribution of peripheral fluorescence [8,22]. Besides increasing the noise in the FCS data, a poor SBR damps the cor-relation amplitude, impeding the determination of concentration and brightness values [23].

In order to quantify the performance of STED in FCS, we introduce a general measure based on the SBR within a homogeneous sample. Note that a spatially homogenous sample represents a well-defined worst case, i.e. non-structured and hence non-resolvable case in imaging, whereas it is virtually the default case in FCS (when the random molecular trajectories are averaged over time). The SBR shall denote the contrast between the signals from the brightest focal region(s) compared to the background from the dimmer region(s), where we apply an arbitrary threshold to separate the bright from the dim region(s). For 50ps pulsed and CW STED irradiation, Fig. 6 exemplifies the SBR versus the radii of the bright signal area for a threshold set at 10% of the maximum fluorescence intensity in the focal center. Without STED, the signal area has a radius of about 203nm and nearly 90% of the detected fluorescence is emitted within this area (SBR ≈8.6). With increasing STED power, the SBR decreases along with the radius of the signal area.

 figure: Fig. 6

Fig. 6 (a) SBR at 10% of maximum threshold. (b) RHS diameters (thin lines: FWHM).

Download Full Size | PDF

The decrease of the SBR is much more pronounced for the CW than for the pulsed STED mode. With a perfect “zero” (ε = 0), the SBR remains >3 for signal area radii >5nm in pulsed mode and 118nm with CW STED, respectively. This significant advantage of pulsed STED was experimentally observed when comparing STED-FCS measurements in pulsed and CW STED mode [8]. With ε = 1% residual intensity in the “zero”, this radius increases sharply from 5nm to 30nm in the pulsed mode, because the bright central region is no longer able to surpass three-fold the background from the entire periphery. In this case, increasing the STED power further can even lead to a complete loss of signal, which is indicated by the return point on the curve in Fig. 6a, where the diameter increases because the first order side lobe (Airy ring) becomes part of the signal area. For the same imperfect “zero” in CW STED, the diameter hardly increases to about 119nm because the peripheral regions always contribute notable background (Fig. 4b).

With the SBR set up here, an effective size of the STED focal region for FCS may be obtained as the area or volume of the brightest focal region(s). Analogous to the FWHM, requiring a SBR = 1 yields a region of half-signal (RHS) from where 50% of the detected signal originates. Figure 6b compares the FWHM diameter with the RHS diameter 2RHS/π, if the outlined system is used to investigate a two-dimensional sample, such as fluorophores diffusing in a membrane sheet in the focal plane. For low to medium STED power, the RHS diameter follows the FWHM with < 20% difference. If k STED τ STED increases beyond the maximum contrast between the “zero” and the crest region, the RHS diameter deviates and eventually grows again, which accounts for a too dim signal from the tiny effective STED PSF. For the outlined example, the minimal RHS diameters are about 4nm (50ps) and 29nm (50ps with ε = 1%) at k STED τ STED ≈104 and 140, respectively.

As pointed out and demonstrated, STED does not require pulsed operation [17]. In the outlined case of a fluorophore with 3.9ns lifetime and a STED laser with 80MHz pulse rate, pulsed excitation but CW STED would just require about five times the average pulsed STED power in order to obtain equivalent FWHM diameters [24,25]. For FCS however, it became clear that the CW STED performs significantly worse because the fluorescence in the peripheral region is less efficiently suppressed. This is indicated by minimal CW RHS diameters of 40nm with perfect and 93nm with ε = 1% “zero”, respectively. With strong CW excitation, the situation may even become worse because fluorophores in the periphery cycle more often than those in the center (c.f. annexes).

4. Conclusions

In summary, we outlined a straightforward analytical model that allows a rather accurate estimation of the inhibition of spontaneous fluorescence by STED and the resulting brightness profiles in the effective STED PSF, i.e. in the region in which the fluorophore remains bright. This region scales down in diameter with a factor proportional to the square root of the number of stimulated emission events per molecule that are in principle possible (maximum statistical average) at the focal intensity maximum of the STED beam. We introduced a combined resolution and contrast measure, the region of half-signal (RHS), stating that the collected signal from within the tiny region in which the dye is allowed to fluoresce should preferably be at least as large as all background contributions from elsewhere, in particular from the focal periphery. For nearly perfect “zeros”, the FWHM diameter continuously shrinks with increasing STED power, whereas the RHS diameter decreases only as long as the fluorophores in the active region generate enough signal to overcome the background from their neighbors. The achievable RHS diameters are smaller in the pulsed than in the CW STED mode, because the pulsed mode better suppresses the excited state in the focal periphery for a given number of stimulating photons.

However, it is worth noting that remedies can be found. If pulsed excitation is used, time-gated detection improves the performance if the early spontaneous emission of the fluorophore is discarded [26]. Likewise, it is important to understand that the outlined RHS ruler fully applies only to the most unfavorable case, i.e. a dense sample with many fluorophores in the focal periphery, each contributing. Whereas this assumption holds for most spectroscopic measurements such as FCS, which typically disregard spatial information in favor of time resolution, it represents a rare worst case in imaging, where the effective PSF (or the corresponding optical transfer function) can be fully exploited for improving the image based on the spatial information discriminating the neighboring features. In both cases, our simple model allows optimizing the experimental conditions for STED nanoscopy as well as taking into account the effective STED PSF for an improved analysis of the measurement data.

Annexes

The probability of spontaneous decay η is given by the ratio of the fluorescence F(γ) under STED versus the fluorescence F(0) without STED. The emission rate per fluorophore k fl is the product of the fluorescence quantum yield qfl and the rate of spontaneous decays k from the excited state S1: Fkfl = qflk.

Pulsed excitation

The derivation of Eq. (3) and its simplifications are outlined in section 2 and summarized as case 2–4 in Table 1 (cumulative conditions). If the STED pulse follows the excitation pulse with a delay τ, 1exp(kS1τ) spontaneous decays take place before the start of the STED pulse and only the remaining fraction exp(kS1τ) of molecules are affected (case 1).

Tables Icon

Table 1. STED efficiency ηps = nS1(γ)/nS1(0) with instantaneous excitation and pulsed or continuous STED irradiation.

As long as the fluorescence response follows the excitation intensity linearly, the fluorescence emission is proportional to the product of the excitation intensity and the spontaneous emission probability. However, for larger excitation intensities, the ground state S0 of the fluorophore becomes depopulated during the excitation pulse (saturated excitation) and/or because of triplet state accumulation. For the sake of simplicity, we keep assuming τex<<1/kS1 (instantaneous excitation) and assume in addition kS1>>kT1 (slow triplet decay).

Let PX be the population of the molecular state X just before the excitation pulse and PX+ just thereafter. Let n S0 be the number of excitations from the ground state S0 during the excitation pulse. The singlet state populations are then described by the following relations:

PS0=1PS1PT1PS1+=PS1+nS0PS0PS1=exp(kS1(γτSTED+T))PS1+
The triplet state T1 is populated via intersystem crossing, who’s number of transitions is the product of the intersystem crossing probability q isc and the number n S1 of spontaneous decays from S1. The triplet population relaxes then with the rate k T1 to the ground state S0.
PT1+PT1+qiscnS1PS1+PT1exp(kT1T)PT1+
At the dynamic equilibrium, Eq. (9) yields the triplet population PT1qiscnS1PS1+/(exp(kT1T)1) and insertion into Eq. (8) yields then a relation for the population of the excited singlet state S1.
PS1+=nS0+(1nS0)exp(kS1(γτSTED+T))PS1+qiscnS1nS0PS1+exp(kT1T)1
Solving Eq. (10) for PS1+ and noting that the rate of spontaneous decays k ps is given by nS1PS1+/T, we finally obtain

kps=nS1nS0T(1(1nS0)exp(kS1(γτSTED+T))+qiscnS1nS0exp(kT1T)1)1.

Equation (11) includes the effect of the triplet state accumulation, which is most pronounced at the local “zero” of the STED beam, where the spontaneous decays are most frequent. Therefore, if a fluorophore shows a significant triplet state population, the accumulation in the dark triplet state will lead to a significant reduction of the fluorophores brightness in the vicinity of the “zero”, much as if there would be some residual STED light.

In addition to the triplet state accumulation, saturated excitation occurs at higher excitation rates k ex because of the depopulation of the ground state S0 during the excitation pulse. The number of excitations n S0 is approximately given by nS0=1exp(kexτex), where we neglect spontaneous decays during the pulse according to our assumptions. However, for high excitation rates, the effective excitation rate kex*=kexkvib/(kex+kvib) should be used. Furthermore, if the STED pulse overlaps the excitation pulse, the (net) number of excitations is reduced to the dynamic equilibrium of the fluorescence excitation and suppression rates and reaches approximately

nS0*=kex*kex*+kS1(1+γ)(1exp((kex*+kS1(1+γ))τex)).

Validity of the effective STED rate

For the sake of simplicity and physical insight, we solved the rate Eqs. (1) at the dynamic equilibrium and found the effective STED rate kSTED*=γkS1=kSTEDkvib/(kSTED+kvib). On the other hand, the exact solution of the time-varying population PS1(t) reads as

PS1(0<t<τSTED)=kS1kvib+k12k1exp(k2+k12t)kS1kvibk12k1exp(k2k12t),
where k1=(kvibkS1)2+4kSTED2 and k2=kvib+kS1+2kSTED. Thereafter, P S1(t) decays with the rate k S1 for the remainder of the pulse period T as outlined before. The integration of q fl k S1 P S1(t) for 0 < t < T is straightforward and yields the probability of fluorescence F(γ), from which we again calculate η ps = F(γ)/F(0). Equation (13) can be simplified by noting that kvib±kS1kvib. However, it is not obvious to get Eq. (3) as the STED rate k STED in the doughnut crest may approach or even surpass k vib, which prohibits any further simplification. Therefore, Fig. 7 compares the exact results with those obtained by Eq. (3) (Fig. 2) showing the excellent agreement at all relevant conditions in the outlined examples. Notable differences occur only at very short STED pulses, at which the condition τSTED>>1/kvib is violated and our assumption of a dynamic (pseudo-)equilibrium during the STED pulse fails.

 figure: Fig. 7

Fig. 7 Probability of spontaneous decay η ps obtained with Eq. (3) compared to the exact results based on Eq. (13) for the very same parameters as described in Fig. 2. The exact results are overlaid as thin lines and are in excellent agreement except at STED pulse widths τ STED < 10ps.

Download Full Size | PDF

Continuous wave excitation

Continuous wave (CW) excitation requires CW STED, which simplifies the calculations because pulse synchronization and timing is not required. The calculation of the spontaneous decay rate k CW and the probability of spontaneous decay η CW are straightforward. Let kSTED* be the effective STED rate kSTED*=γkS1=kSTEDkvib/(kSTED+kvib) as defined before, and let qisc*=qisckS1/(kS1+kSTED*)=qisc/(1+γ) be the corresponding intersystem crossing yield. At the dynamic equilibrium, the average cycle time τ CW is then given by the sum of the occupation times of the singlet and triplet states.

τCW=1kex*+1kS1+kSTED*+qisc*kT1=1kex*+1kS1(1+γ)+qisckT1(1+γ)
The second term is the average time τS1* the molecule spends in the excited state S1. Thus, the average population of S1 becomes PS1=τS1*/τCW and the spontaneous decay rate is kCW=kS1PS1. We finally obtain
kCW=(1+γkex*+1kS1+qisckT1)1andηCW=(1kex*+1kS1+qisckT1)kCW.
If the triplet state population is small, i.e. qisc/kT1<<1/kS1, ηCW(kex*+kS1)/(kex*+(1+γ)kS1); and if in addition the excitation rate is small, i.e. kex*<<kS1, we obtain again ηCW1/(1+γ) as for the case 4 in Table 1. This convergence at low CW excitation rate can be expected because individual excitations become well separated in time and, therefore, the fluorophores undergo STED as if a pulsed excitation scheme were used.

Acknowledgments

Marcel Leutenegger is a fellow of the Alexander von Humboldt foundation.

References and links

1. S. W. Hell and J. Wichmann, “Breaking the diffraction resolution limit by stimulated emission: stimulated-emission-depletion fluorescence microscopy,” Opt. Lett. 19(11), 780–782 (1994). [CrossRef]   [PubMed]  

2. T. A. Klar, S. Jakobs, M. Dyba, A. Egner, and S. W. Hell, “Fluorescence microscopy with diffraction resolution barrier broken by stimulated emission,” Proc. Natl. Acad. Sci. U.S.A. 97(15), 8206–8210 (2000). [CrossRef]   [PubMed]  

3. S. W. Hell, “Far-field optical nanoscopy,” Science 316(5828), 1153–1158 (2007). [CrossRef]   [PubMed]  

4. A. Thomas, Progress in Stimulated Emission Depletion Microscopy (Shaker, Achen 2001); Ph.D. thesis, Heidelberg (2001).

5. S. W. Hell, “Increasing the Resolution of Far-Field Fluorescence Microscopy by Point-Spread-Function Engineering,” in Topics In Fluorescence Spectroscopy 5: Nonlinear and Two-Photon-Induced Fluorescence, J. Lakowicz, ed. (Plenum Press, New York, 1997).

6. B. Harke, J. Keller, C. K. Ullal, V. Westphal, A. Schönle, and S. W. Hell, “Resolution scaling in STED microscopy,” Opt. Express 16(6), 4154–4162 (2008). [CrossRef]   [PubMed]  

7. The actual number of de-excitation events is limited by the number of excitations per pulse.

8. C. Ringemann, B. Harke, C. von Middendorff, R. Medda, A. Honigmann, R. Wagner, M. Leutenegger, A. Schönle, S. W. Hell, and C. Eggeling, “Exploring single-molecule dynamics with fluorescence nanoscopy,” N. J. Phys. 11(10), 103054 (2009). [CrossRef]  

9. V. Westphal and S. W. Hell, “Nanoscale resolution in the focal plane of an optical microscope,” Phys. Rev. Lett. 94(14), 143903 (2005). [CrossRef]   [PubMed]  

10. V. Westphal, C. M. Blanca, M. Dyba, L. Kastrup, and S. W. Hell, “Laser-diode-stimulated emission depletion microscopy,” Appl. Phys. Lett. 82(18), 3125–3127 (2003). [CrossRef]  

11. M. Dyba and S. W. Hell, “Photostability of a fluorescent marker under pulsed excited-state depletion through stimulated emission,” Appl. Opt. 42(25), 5123–5129 (2003). [CrossRef]   [PubMed]  

12. C. Eggeling, J. Widengren, R. Rigler, and C. A. M. Seidel, “Photobleaching of fluorescent dyes under conditions used for single-molecule detection: Evidence of two-step photolysis,” Anal. Chem. 70(13), 2651–2659 (1998). [CrossRef]   [PubMed]  

13. W. Denk, J. H. Strickler, and W. W. Webb, “Two-photon laser scanning fluorescence microscopy,” Science 248(4951), 73–76 (1990). [CrossRef]   [PubMed]  

14. K. I. Willig, J. Keller, M. Bossi, and S. W. Hell, “STED microscopy resolves nanoparticle assemblies,” N. J. Phys. 8(6), 106 (2006). [CrossRef]  

15. M. Leutenegger, R. Rao, R. A. Leitgeb, and T. Lasser, “Fast focus field calculations,” Opt. Express 14(23), 11277–11291 (2006). [CrossRef]   [PubMed]  

16. M. Leutenegger and T. Lasser, “Detection efficiency in total internal reflection fluorescence microscopy,” Opt. Express 16(12), 8519–8531 (2008). [CrossRef]   [PubMed]  

17. K. I. Willig, B. Harke, R. Medda, and S. W. Hell, “STED microscopy with continuous wave beams,” Nat. Methods 4(11), 915–918 (2007). [CrossRef]   [PubMed]  

18. S. W. Hell, “Toward fluorescence nanoscopy,” Nat. Biotechnol. 21(11), 1347–1355 (2003). [CrossRef]   [PubMed]  

19. X. Hao, C. Kuang, T. Wang, and X. Liu, “Effects of polarization on the de-excitation dark focal spot in STED microscopy,” J. Opt. 12(11), 115707 (2010). [CrossRef]  

20. D. Magde, E. Elson, and W. Webb, “Thermodynamic fluctuations in a reacting system - measurement by fluorescence correlation spectroscopy,” Phys. Rev. Lett. 29(11), 705–708 (1972). [CrossRef]  

21. R. Rigler, and E. S. Elson, Fluorescence Correlation Spectroscopy: Theory and Applications, Springer Ser. Chem. Phys. 65, Berlin Heidelberg, Germany (2001).

22. L. Kastrup, H. Blom, C. Eggeling, and S. W. Hell, “Fluorescence fluctuation spectroscopy in subdiffraction focal volumes,” Phys. Rev. Lett. 94(17), 178104 (2005). [CrossRef]   [PubMed]  

23. D. E. Koppel, “Statistical Accuracy in Fluorescence Correlation Spectroscopy,” Phys. Rev. A 10(6), 1938–1945 (1974). [CrossRef]  

24. In case of a long-lived excited state, i.e. ≈12ns in a nano-diamond nitrogen vacancy center [25], the required CW power is only about twice the average power for pulsed operation.

25. K. Y. Han, K. I. Willig, E. Rittweger, F. Jelezko, C. Eggeling, and S. W. Hell, “Three-dimensional stimulated emission depletion microscopy of nitrogen-vacancy centers in diamond using continuous-wave light,” Nano Lett. 9(9), 3323–3329 (2009). [CrossRef]   [PubMed]  

26. E. Auksorius, B. R. Boruah, C. Dunsby, P. M. P. Lanigan, G. Kennedy, M. A. A. Neil, and P. M. W. French, “Stimulated emission depletion microscopy with a supercontinuum source and fluorescence lifetime imaging,” Opt. Lett. 33(2), 113–115 (2008). [CrossRef]   [PubMed]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (7)

Fig. 1
Fig. 1 Simplified Jablonski diagram with electronic and vibrational energy levels of a typical fluorophore.
Fig. 2
Fig. 2 Probability of spontaneous decay η ps as a function of the saturation factor ζ and the pulse width τ STED of the STED light (as labeled) for k S1 = 1/3.9ns, k vib = 5/ps and T = 1/80MHz. (a) The thin line shows the lower limit achieved for CW STED and infinitely fast vibrational relaxation. Thin dashed lines show the dependence of η ps on variable pulse lengths τ STED for fixed numbers of potential de-excitations k STED τ STED, i.e. for fixed average STED powers. (b) This dependence is shown as a function of τ STED. The inset outlines the pulse scheme used in the calculations: a Dirac excitation pulse followed by a rectangular STED pulse.
Fig. 3
Fig. 3 Cross-sections through the foci of the excitation and the STED light spots (a) and of the fluorescence collection efficiency (b) used in the calculation examples. Insets show the spot profiles in the xy-planes.
Fig. 4
Fig. 4 Effective PSF of a STED nanoscope with (a) pulsed (τ STED = 50ps) and (b) CW STED (τ STED = T = 12.5ns): lateral cross-sections of the detected fluorescence brightness without STED (bold solid lines) and with average STED powers required for achieving FWHM diameters of 72nm (thin solid lines) and 33nm (pointed red lines). Dash-dotted lines show the profiles with an imperfect “zero” and dotted lines show them normalized again to their peak brightness. For values below 10% of the peak brightness, the logarithmic scale allows a better distinction of the differences.
Fig. 5
Fig. 5 FWHM diameters versus STED power. For k STED τ STED > 500 in case of 50ps STED pulses with ε = 1%, the increase of the FWHM indicates loss of signal due to a too strong reduction of the central peak brightness versus the contributions from the focal periphery.
Fig. 6
Fig. 6 (a) SBR at 10% of maximum threshold. (b) RHS diameters (thin lines: FWHM).
Fig. 7
Fig. 7 Probability of spontaneous decay η ps obtained with Eq. (3) compared to the exact results based on Eq. (13) for the very same parameters as described in Fig. 2. The exact results are overlaid as thin lines and are in excellent agreement except at STED pulse widths τ STED < 10ps.

Tables (1)

Tables Icon

Table 1 STED efficiency η ps = n S1(γ)/n S1(0) with instantaneous excitation and pulsed or continuous STED irradiation.

Equations (15)

Equations on this page are rendered with MathJax. Learn more.

t P S1 = k S1 P S1 k STED ( P S1 P S 0 * ) t P S0 * = k vib P S0 * + k STED ( P S1 P S0 * )
F ( γ ) = q fl k S1 0 T P S1 ( t ) d t = q fl ( 1 + γ exp ( k S1 τ STED ( 1 + γ ) ) 1 + γ exp ( k S1 ( γ τ STED + T ) ) ) .
η ps = ( 1 + γ exp ( k S1 τ STED ( 1 + γ ) ) 1 + γ exp ( k S1 ( γ τ STED + T ) ) ) ( 1 exp ( k S1 T ) ) 1
η ps 1 + γ exp ( k S1 τ STED ( 1 + γ ) ) 1 + γ .
η ps 1 / ( 1 + γ ) .
η CW k ex + k S1 k ex + ( 1 + γ ) k S1 .
η ps * = exp ( k S1 τ STED γ ) = exp ( ln ( 2 ) I STED / I s * )
P S0 = 1 P S1 P T1 P S1 + = P S1 + n S0 P S0 P S1 = exp ( k S1 ( γ τ STED + T ) ) P S1 +
P T1 + P T1 + q isc n S1 P S1 + P T1 exp ( k T1 T ) P T1 +
P S1 + = n S0 + ( 1 n S0 ) exp ( k S1 ( γ τ STED + T ) ) P S1 + q isc n S1 n S0 P S1 + exp ( k T1 T ) 1
k ps = n S1 n S0 T ( 1 ( 1 n S0 ) exp ( k S1 ( γ τ STED + T ) ) + q isc n S1 n S0 exp ( k T1 T ) 1 ) 1
n S0 * = k ex * k ex * + k S1 ( 1 + γ ) ( 1 exp ( ( k ex * + k S1 ( 1 + γ ) ) τ ex ) ) .
P S1 ( 0 < t < τ STED ) = k S1 k vib + k 1 2 k 1 exp ( k 2 + k 1 2 t ) k S1 k vib k 1 2 k 1 exp ( k 2 k 1 2 t ) ,
τ CW = 1 k ex * + 1 k S1 + k STED * + q isc * k T1 = 1 k ex * + 1 k S1 ( 1 + γ ) + q isc k T1 ( 1 + γ )
k CW = ( 1 + γ k ex * + 1 k S1 + q isc k T1 ) 1 and η CW = ( 1 k ex * + 1 k S1 + q isc k T1 ) k CW .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.