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
Beyond Rayleigh's criterion: A resolution measure with application to singlemolecule microscopy

Edited by Robert H. Austin, Princeton University, Princeton, NJ, and approved January 4, 2006 (received for review September 14, 2005)
Abstract
Rayleigh's criterion is extensively used in optical microscopy for determining the resolution of microscopes. This criterion imposes a resolution limit that has long been held as an impediment for studying nanoscale biological phenomenon through an optical microscope. However, it is well known that Rayleigh's criterion is based on intuitive notions. For example, Rayleigh's criterion is formulated in a deterministic setting that neglects the photon statistics of the acquired data. Hence it does not take into account the number of detected photons, which, in turn, raises concern over the use of Rayleigh's criterion in photoncounting techniques such as singlemolecule microscopy. Here, we reexamine the resolution problem by adopting a stochastic framework and present a resolution measure that overcomes the limitations of Rayleigh's criterion. This resolution measure predicts that the resolution of optical microscopes is not limited and that it can be improved by increasing the number of detected photons. Experimental verification of the resolution measure is carried out by imaging singlemolecule pairs with different distances of separation. The resolution measure provides a quantitative tool for designing and evaluating singlemolecule experiments that probe biomolecular interactions.
According to Rayleigh's criterion, the resolution of an optical microscope is defined as the minimum distance between two point sources such that their presence can be distinguished in the image (1). It is widely believed that the resolution limit imposed by this criterion precludes the singlemolecule study of molecular interactions at distances of <200 nm. Despite the widespread use of Rayleigh's criterion, it is well known that this criterion is based on heuristic notions (2). Formulated within a deterministic framework, Rayleigh's criterion neglects the stochastic nature of the photondetection process and hence does not consider the total photon count in the acquired data. This formulation is not surprising, because Rayleigh's criterion was developed at a time when the unaided human eye was used as the detector. Therefore, Rayleigh's criterion is not well adapted to current microscope setups, in which highly sensitive photoncounting cameras are used. Recent singlemolecule experiments have shown that Rayleigh's criterion can, in fact, be surpassed in an optical microscope setup (3–5). Thus Rayleigh's resolution limit is an inadequate performance criterion for current quantitative imaging techniques.
By adopting a stochastic framework, we propose a resolution measure that overcomes the shortcomings of Rayleigh's criterion and provides a quantitative measure of a microscope's ability to determine the distance between two point sources. Unlike Rayleigh's criterion our resolution measure predicts that the resolution of a microscope can be improved by increasing the number of photons collected from the point sources. The resolution measure is given in terms of quantities such as the expected number of detected photons, numerical aperture of the objective lens, and wavelength of the detected photons. We also investigate how the resolution measure is influenced by various experimental factors that affect the acquired data such as pixelation of the detector. Experimental results are presented by estimating distances from images of closely spaced single molecules. These results show that distances well below Rayleigh's resolution limit can be determined with an accuracy as specified by our resolution measure.
Results
Fundamental Resolution Measure (FREM).
Our approach to the derivation of the resolution measure is to obtain a bound/limit to the accuracy with which the distance between two point sources can be estimated based on the acquired data (see Materials and Methods). Analogous to Rayleigh's criterion, we consider an optical microscope setup that images two identical, selfluminous, infocus point sources emitting unpolarized, incoherent light. The analytical expression of the FREM for this imaging condition is given by
where λ denotes the emission wavelength of the detected photons, n
_{a} denotes the numerical aperture of the objective lens, Λ_{0} denotes the photon detection rate (intensity) per point source, [t
_{0}, t] denotes the acquisition time interval, and Γ_{0}(d) is given by
with J
_{n} denoting the nth order Bessel function of the first kind, α = 2πn
_{a}/λ; r
_{01} =
The stochastic framework used to obtain the FREM models the photon emission (detection) process as a random process (shot noise process). The spatial locations at which the photons hit the detector are assumed to be randomly distributed according to the image profiles of the point sources/single molecules. This framework considers an optical microscope setup in which the detector provides the time points and the spatial coordinates of every detected photon without adding any extraneous noise. For any imaging condition, this setup can be thought of as an idealization of current imaging detectors in which the presence of finitesized pixels and measurement noise deteriorates the acquired data. Thus the resolution measure derived within this framework provides a result that is fundamental for the given imaging condition. In the present context, the FREM is obtained for imaging conditions analogous to those of Rayleigh's criterion. Hence the spatial distribution of the detected photons from each point source is described by the Airy profile (1).
The FREM predicts how accurately the distance d between two point sources can be resolved. A small numerical value for the FREM predicts a high accuracy in determining d, whereas a large numerical value of the FREM predicts a low accuracy in determining d. Fig. 1 A shows the behavior of the FREM as a function of the distance of separation between a pair of GFP molecules (λ=520 nm) and for a pair of Cy5 molecules (λ=690 nm) that are imaged with an objective lens of numerical aperture 1.45. In Fig. 1 A, it is assumed that the expected photon count is the same for both fluorophores. For the GFP molecules, Rayleigh's criterion predicts the smallest resolvable distance to be ≈220 nm (≈0.61 λ/n _{a}). In contrast, Fig. 1 A shows that the FREM has a small numerical value for distances ranging from 50 to 220 nm, which are well below Rayleigh's criterion. For distances <50 nm, however, the FREM deteriorates (i.e., increases) significantly with decreasing distance of separation (see Fig. 1 A Inset). In particular, as the distance of separation decreases to zero, the FREM becomes infinitely large, because the term Γ_{0}(d) (Eq. 2 ), which appears in the denominator of the FREM, tends to zero. An analogous behavior of the FREM is also seen for the Cy5 molecules. Note that the numerical value of the FREM for the Cy5 molecules is consistently larger than that of the GFP molecules for the same expected photon count per fluorophore. For example, the FREM predicts that a distance of 10 nm between two GFP molecules can be determined with an accuracy not better than ±5.7 nm when the expected photon count per GFP molecule is 3,000. On the other hand, for the same distance of separation and photon count per molecule, the FREM predicts an accuracy not better than ±9.1 nm for the Cy5 molecules. In the case of the Cy5 molecules, however, the numerical value of the FREM is comparable to the distance of separation itself. Because the FREM exhibits an inverse square root dependence on the expected number of detected photons, this deterioration can be compensated for by increasing the expected number of detected photons, as shown in Fig. 1 B. Thus in the above example, if we increase the expected photon count per Cy5 dye molecule to 10 ^{4} , then the FREM predicts that a distance of 10 nm can be determined with an accuracy not better than ±5 nm.
Practical Resolution Measure (PREM).
The FREM provides the bestcase scenario for a microscope setup, where experimental factors that potentially deteriorate the acquired data were not taken into account. We next investigate how the resolution measure is affected by such experimental factors. Here we obtain an analytical expression for the resolution measure that takes into account these experimental factors. We refer to this result as the PREM. The PREM can be thought of as an extension to the FREM. For instance, the PREM takes into account the presence of additive noise sources, namely Poisson and Gaussian noise. Poisson noise is used to model spurious photons in the acquired image, which, for example, arise because of the autofluorescence of the sample and dark current of the detector (6). Gaussian noise is used to model measurement noise in the acquired data, which, for example, arise during the readout process in the detector (6). The additive Poisson noise considered here is distinct from the shot noise, which describes the statistics of the photondetection process from the single molecules and is already accounted for by the FREM. Aside from these extraneous noise sources, the PREM also takes into account the effect of pixelation of the detector.
Fig. 2 A shows the behavior of the PREM as a function of the distance between two Cy5 molecules in the presence and absence of noise sources for a pixelated detector. Fig. 2 A also shows the FREM for reference. Note that even in the absence of extraneous noise sources the numerical value of the PREM is consistently greater than that of the FREM because of the pixelation of the detector. Moreover, in the presence of noise sources this behavior of the PREM becomes more pronounced. In particular, for very small distances (≤50 nm), the numerical value of the PREM is at least three to five times greater than that of the FREM (see Fig. 2 A Inset). Analogous to Fig. 1 B, the deterioration of the PREM at very small distances can be compensated for by collecting more photons from the point sources (Fig. 2 B). In contrast, for distances ranging from 100 to 250 nm, which are below Rayleigh's criterion (≈0.61 λ/n _{a} ≈ 290 nm), the numerical value of the PREM approaches that of the FREM even in the presence of noise, as shown in Fig. 2 A. As an application of these results, consider a practical scenario in which we require distances ranging from 50 to 200 nm to be resolved between two Cy5 molecules with an accuracy of at least 5 nm. From Fig. 2 B, we know that to estimate a distance of 50 nm with an accuracy not better than 5 nm the PREM predicts the expected number of detected photons per single molecule to be at least 15,000. On the other hand, from Fig. 2 A we see that to estimate a distance of 200 nm with similar accuracy the PREM predicts the expected number of detected photons per single molecule to be at least 2,500. Hence to resolve distances ranging from 50 to 200 nm between two Cy5 molecules with an accuracy not better than 5 nm on average at least 15,000 photons must be collected per single molecule.
Improving the Resolution Measure by Using Additional Spatial Information.
It is shown in Figs. 1 B and 2 B that the resolution measure can be improved by increasing the number of detected photons from each of the point sources. In singlemolecule experiments increasing the photon count is not always possible, because the fluorophores may photobleach. However, for a singlemolecule pair that exhibits a doublestep photobleaching behavior, additional information can be obtained from the photons collected from the fluorophore that remains after the first photobleaching event (3). Fig. 2 C shows the behavior of the resolution measure for a pair of Cy5 molecules spaced 10 nm apart by taking into account the number of photons collected before and after the first photobleaching event. Here the resolution measure is determined for an imaging condition with numerical values analogous to those used in Fig. 2 A. From Fig. 2 A we see that the resolution measure predicts an accuracy not smaller than ±5.8 nm to determine a distance of 10 nm, when on average 5,000 photons are collected from each fluorophore before and after the first photobleaching event. This result is in contrast to the case when the additional information obtained after the first photobleaching event is not used. In this case, the resolution measure predicts an accuracy not smaller than ±43.0 nm to resolve a distance of 10 nm for the same photon count per fluorophore.
Experimental Verification.
The resolution measure provides a bound/limit to the smallest possible standard deviation of any unbiased estimator of the distance between the point sources. To verify whether this bound can be attained in experiments, images of closely spaced Cy5 molecules were collected, and their distances of separation were estimated by using the maximumlikelihood estimator. According to Rayleigh's criterion, the minimum resolvable distance is given by 0.61 λ/n _{a}, which, in the present case is ≈290 nm. Table 1lists the results of distance estimation and the predicted resolution measure for two pairs of single molecules. One of the singlemolecule pairs has a mean distance of separation of 293 nm (data analysis 1), which is close to Rayleigh's criterion, and the other singlemolecule pair has a mean distance of separation of 207 nm (data analysis 2), which is below Rayleigh's criterion. From Table 1 we see that for each data analysis the standard deviation of the maximumlikelihood estimates of the distance comes close to the resolution measure. Note that the numerical values of the standard deviations are themselves only estimates based on the acquired data. With larger data sets, the agreement with the resolution measure is expected to increase further.
The above data sets were also analyzed by estimating the distances of separation through the global maximumlikelihood estimator (Fig. 3, which is published as supporting information on the PNAS web site), which used the additional spatial information available in the acquired data, i.e., the images collected after the first photobleaching event. Table 1 lists the results of the distance estimates (data analyses 3 and 4) for the Cy5 singlemolecule pairs analyzed above. From Table 1 we see that the accuracy of the distance estimates obtained by using the additional spatial information is consistently smaller than the accuracy obtained when the additional information is not used (data analyses 1 and 2). Table 1 also shows an analogous behavior of the resolution measure for the corresponding data sets. For example, in the case when the additional spatial information is not used, the standard deviation of the distance estimates for the Cy5 singlemolecule pair with a mean distance of separation of 207 nm is equal to ±10.1 nm (data analysis 2). On the other hand, for the same singlemolecule pair, when additional spatial information is used the standard deviation of the distance estimates is equal to ±3.8 nm (data analysis 4). Table 1 also lists the mean endtoend distance estimates of the DNA molecular ruler, which were determined by using the global estimation approach (see Fig. 4, which is published as supporting information on the PNAS web site). From Table 1 it can be seen that for each DNA data set the standard deviation of the global maximumlikelihood estimator comes close to the resolution measure. For example, in data analysis 6 the standard deviation of the distance estimates is equal to ±8.7 nm and the resolution measure predicts an accuracy not smaller than ±7.4 nm to resolve the distance of 12 nm. We note that in the same data set, if the additional spatial information is not used, then the resolution measure predicts an accuracy not smaller than ±52.1 nm to resolve a distance of 12 nm.
Generalization of the FREM.
The FREM given in Eq. 1 was derived for imaging conditions analogous to those of Rayleigh's criterion, which considered two equal intensity, infocus point sources that emit unpolarized, incoherent light. However, in several applications these conditions are not met, for example, when using polarized illumination and detection (7). We now consider a situation where the point sources can potentially have unequal intensities that vary as a function of time, and where the image profiles of the point sources can be distinct. The expression for the generalized FREM (gFREM) is given by where q _{1} and q _{2} denote the image functions of the point sources and Λ_{1} and Λ_{2} denote the intensities of the point sources. In many situations the image of the point source significantly differs from the Airy profile, for example, because of the defocus in the objective lens (8), the different orientations of the pointsource emission dipole (7, 9), or the aberrations present in the imaging setup (10). Moreover, depending on the nature of illumination, the intensity of the point sources can be unequal when their emission dipole orientations are different (7, 9). We note that Eq. 3 provides a general expression for the FREM that is applicable to a wide variety of imaging conditions, including the abovementioned scenarios. If in Eq. 3 we set the intensities to be constant and identical, i.e., Λ_{1}(τ)=Λ_{2}(τ)=Λ_{0}, τ ≥ t _{0}, and assume the image functions to be given by the Airy profile, then we immediately obtain the expression for the FREM given in Eq. 1 . Analogous to the gFREM, an expression has also been derived for the generalized PREM (see Materials and Methods).
The stochastic framework used to derive the gFREM models the photon detection process for each point source as a Poisson process (shot noise process). Recently, there have been reports of the generation of nonclassical states of light from a fluorescent light source in which the photon statistics deviate from the classical shot noise process (11–13). To take into account this deviation, a further generalization of the FREM can be obtained by modeling the photon detection process as a general counting process (ref. 14 and Supporting Text, which is published as supporting information on the PNAS web site).
Discussion
The advent of singlemolecule microscopy has generated significant interest in studying nanoscale biomolecular interactions. Classically, fluorescence resonance energy transferbased methods have been used to probe interactions in the distance range of 1 to 10 nm (15). It is widely believed that Rayleigh's criterion precludes the resolution of two single molecules at distances of <200 nm, which leaves a gap in the distance range of 10 to 200 nm that is vital for the study of many biological processes with an optical microscope. It has been suggested that Rayleigh's resolution limit can be superseded if the distance between two point sources is determined by curve fitting the image with the sum of two pointsource image profiles (16). In fact, by adopting this approach several groups have shown that Rayleigh's limit can be surpassed in experiments (3–5).
The FREM (Eq. 1 ) is a resolution measure that overcomes several deficiencies of Rayleigh's criterion. It gives a bound for the accuracy with which the distance between two point sources can be estimated when the acquired data are not affected by deteriorating experimental factors. An important property of the FREM is that it provides a quantitative assessment of how the optical characteristics of the experimental setup and the photon budget influence the resolution performance in determining a particular distance of separation. Fig. 1 A shows that the numerical value of the FREM for a pair of GFP molecules is consistently smaller than that for a pair of Cy5 molecules when the expected photon count per fluorophore is 3,000 in both cases. For example, to resolve distances of 8, 50, and 200 nm between a pair of GFP molecules, the FREM predicts an accuracy not smaller than ±6.4, ±2.7, and ±1.9 nm, respectively. On the other hand, for a pair of Cy5 molecules with the same expected photon count per fluorophore and distances, the FREM predicts an accuracy not smaller than ±10.1, ±4.0, and ±2.7 nm, respectively. For the 50 and 200nm distances, the numerical values of the FREM for the GFP and Cy5 molecules are significantly smaller than the corresponding actual distances. This observation implies that the FREM predicts a relatively high accuracy in resolving distances ranging from 50 to 200 nm between single molecules. For the 8nm distance, however, the numerical value of the FREM for the GFP and Cy5 molecules are either comparable to or greater than the actual distance, which suggests that even in the bestcase scenario, i.e., in the absence of deteriorating experimental factors, distances of <8 nm are difficult to resolve between the GFP/Cy5 molecules, unless a higher than average number of photons is detected. In singlemolecule experiments typically 3,000 photons can be collected before a GFP molecule irreversibly photobleaches (see e.g., refs. 17 and 18), whereas for Cy5 molecules typically >10 ^{4} photons can be collected (see ref. 19). Thus, to resolve distances of 8, 50, and 200 nm with an expected photon count of 10 ^{4} per Cy5 molecule, the FREM predicts an accuracy not smaller than ±5.5, ±2.2, and ±1.5 nm, respectively, i.e., a ≈2fold improvement from the case when on average 3,000 photons are collected per Cy5 molecule. Note that the localization accuracy of a single molecule, i.e., the accuracy with which the position of a single molecule can be determined also depends on the number of collected photons (20, 21).
The PREM derived here extends the results of the FREM by illustrating how the resolution measure is deteriorated by experimental factors such as pixelation of the detector and extraneous noise sources. The PREM, i.e., the bound on the accuracy with which the distance can be estimated, for typical imaging conditions (see Fig. 2 A) is given by ±31.6, ±5.3, and ±2.2 nm for the case when the Cy5 single molecules are 8, 50, and 200 nm apart, respectively, and the expected photon count per fluorophore is 10 ^{4} . For the same distances, if the expected photon count is 3,000 per Cy5 molecule, the PREM is significantly higher at ±76.7, ±12.5, and ±4.5 nm, respectively. Similarly, for a pair of GFP molecules with an expected photon count of 3,000 per molecule, the PREM predicts an accuracy not smaller than ±42.0, ±7.4, and ±3.0 nm to resolve distances of 8, 50, and 200 nm, respectively. This calculation shows that especially for small distances the predicted resolution measure is probably not acceptable for many applications. This deterioration in the limit of the accuracy with which the distance can be measured is because the acquired data are a discretized version of the actual image, and the presence of extraneous noise sources corrupts the acquired data (e.g., scattered photons, noise in the acquisition electronics). Moreover, a comparison with the FREM illustrates that control of the noise sources is also of great importance to improve the accuracy of the estimated distance parameter in a practical scenario (Fig. 2 A).
The above results suggest that the distance between two single molecules can be estimated with a reasonable level of accuracy, depending on the photon count, certainly for distances >50 nm, but possibly also for smaller distances. For distances ≈10 nm, however, the predicted resolution measures are typically worse. To be able to resolve such distances the number of detected photons would have to be increased substantially, which is typically not possible because of photobleaching. For instance, if a PREM of 7 nm is to be achieved for a Cy5 singlemolecule pair spaced 10 nm apart, then under the noisy imaging conditions of Fig. 2 A at least 1.5 × 10 ^{5} photons need to be detected per Cy5 molecule, which underscores the importance of the development of brighter and more photostable fluorescent markers to carry out such studies (7, 22).
In ref. 23 it was shown that GFP singlemolecule pairs typically photobleach together. However, in the case where twostep photobleaching occurs, additional information can be used by imaging the remaining single molecule (3). The development of the FREM and the PREM discussed so far was based on the case when both point sources/single molecules do not bleach during the acquisition of the image. The approach can, however, also be applied to investigate the accuracy that can be achieved when additional information is taken into account from the remaining single molecule that did not bleach in the first photobleaching step. For example, consider a pair of GFP molecules that are 8, 15, 25, or 50 nm apart. Assume that 2,000 photons are collected from the single molecules before and after the first photobleaching event. If the photons collected before and after the first photobleaching event are taken into account, and assuming the experimental conditions for pixelation and noise sources of Fig. 2 A, then an accuracy of no better than ±8.1, ±7.9, ±7.6, and ±6.7 nm, respectively, can be expected for these distances. Under the present assumptions on the expected number of detected photons, this calculation shows that for distances up to 15 nm probably unreasonably large errors would be incurred in the estimation. However, for distances >25 nm an error level of <30% could be expected.
The resolution measure provides a bound to the accuracy/standard deviation with which the distance between two point sources can be estimated. This was experimentally verified by imaging closely spaced Cy5 molecules and estimating their distance of separation from the acquired data. Here, the maximumlikelihood estimator was used, because it possesses favorable properties for estimating parameters (24). In general, other estimation algorithms can also be used for determining the distance of separation between two point sources. However, a question that arises is which of the different estimation algorithms is the most suitable for estimating the distance of separation. In such a scenario, the knowledge of the resolution measure becomes crucial, because it can be used as a standard to compare the performance of the different estimation algorithms. The experimental results presented in Table 1 show that the standard deviations of the maximumlikelihood distance estimates come close to the bound predicted by the resolution measure, thereby validating the choice of this estimator.
The FREM was derived for imaging conditions that were analogous to those assumed in Rayleigh's criterion. In some singlemolecule experiments, however, the conditions are different to those assumed in the derivation of Rayleigh's criterion and which formed the basis for the derivation of the FREM. Whereas the FREM assumes that the image of a point source is given by an Airy profile, the gFREM was derived so that more complex image profiles can be analyzed. Such profiles could arise, for example, because of outoffocus conditions (8), the presence of aberrations (10), or the use of polarized illumination (7). We note that the gFREM can also be used to calculate the resolution measure for determining the distance of separation between any two (distinct) objects such as cellular organelles, provided the intensities and the image functions of the objects are known. The resolution limit of Rayleigh's criterion is thought to arise because of the finite width of the central peak of the pointsource image, i.e., the point spread function (2). This approach led to the development of microscopy techniques such as 4Pi confocal microscopy (25), stimulated emission depletion microscopy (26), and image interference microscopy (27) in which the width of the central peak of the point spread function is smaller than that of the conventional optical microscope. These techniques have reported improvement in resolving features that are typically unresolvable in conventional optical microscopes. However, for some of the techniques it was reported that this improvement was achieved at a severe cost to the signal (i.e., number of collected photons) (28). This observation illustrates the importance of considering the photon/light budget when discussing resolution performance, especially in fluorescence imaging applications that typically use photobleachable fluorophores. In the present context, the expression for the FREM/gFREM explicitly shows the tradeoff between the intensities of the point sources, which determine the photon budget, and the image functions of the point sources, which determine the point spread function shape.
Materials and Methods
SingleMolecule Microscopy.
A molecular ruler consisting of a 30bp DNA duplex was used. The oligonucleotide 5′ATC TCG GTG CGT AAT ACT CAC GGG CAG GAC3′ (ref. 29 and D. Holowka, personal communication) and its complementary sequence (both labeled with Cy5 at the 5′ end) were purchased from Synthegen (Houston). The oligonucleotides were annealed in 200 mM Tris·HCl, 10 mM MgCl_{2}, pH 8.0 and stored at 4°C. Molecular modeling by the vendor of the DNA duplex labeled with Cy5 dye at both ends predicts the distance between the two Cy5 dyes to be 12 nm. Cy5 dye purchased from GE Healthcare (Piscataway, NJ) was used for calibration purposes. To image single DNA molecules, a cleaned glassbottomed dish (MatTek, Ashland, MA) was coated with a layer of polylysine (0.01% solution) followed by the fluorescent sample at a concentration of 1 pM. Imaging experiments were carried out in a custom setup that was built on a Zeiss Axiovert S100 fluorescence microscope. The setup consisted of a 643nm laser (Research ElectroOptics, Boulder, CO), a cooled chargecoupled device camera (ORCAER, Hamamatsu, Bridgewater, NJ), and an αPlanFLUAR (numerical aperture 1.45, ×100) Zeiss objective lens. The sample was illuminated with circularly polarized light in widefield mode, the camera was operated in 2 × 2 binning mode, and the exposure time was either 1 or 3 s.
Theory.
Let q
_{1}, q
_{2} be the image functions of the two point sources, i.e., their images at the origin of the coordinate system assuming unit magnification and normalization ∫_{ℝ2}
q
_{i}(x, y)dxdy=1, i=1, 2. Let Λ_{1}(τ), Λ_{2}(τ), τ ≥ t
_{0}, be the intensities of the Poisson processes that model the time points of the detected photons. If the two single molecules are located a distance d apart, the impact locations of the photons on the detector plane are randomly distributed with density
where M is the magnification and the parameter to be estimated is θ=d. For this data model the Fisher information matrix is given by (14, 30)
where κ_{θ}(r, t)=(Λ_{1}(τ) + Λ_{2}(τ))f
_{θ,τ}(r), r∈ℝ^{2}, τ ≥ t
_{0}. According to the CramerRao lower bound (14, 24), for the standard deviation of any unbiased estimator θ of θ we have
The data model that describes the image acquisition with a pixelated detector with N _{p} pixels during the integration time interval [t _{0}, t] is given by independent random variables ℐ_{θ,k}=S _{θ,k} + B_{k} + W_{k} , where S _{θ,k} (B_{k} ) is a Poisson random variable with mean μθ(k, t) (β(k, t)) that models the detected photon count at the kth pixel from the point sources (spurious sources), and W_{k} is an independent Gaussian random variable with mean η_{k} and variance σ_{w,k} ^{2} that models the readout noise at the kth pixel, k=1, …, N_{p} . The general expression for the Fisher information matrix I(θ) is given by (21, 30) where and v _{θ}(k, t)=μ_{θ}(k, t) + β(k, t), k=1, …, N _{p}, θ ∈ Θ.
To obtain the Fisher information matrix for the resolution problem, the general expression for I(θ) is evaluated with where C_{k} denotes the kth pixel, (x _{0}, y _{0}) denotes the location of one of the point sources, φ denotes the angle subtended by the line joining the two point sources, and θ=(x _{0}, y _{0}, d, φ) (Fig. 5, which is published as supporting information on the PNAS web site).
According to the CramerRao lower bound, the bound on the standard deviation of any unbiased estimator of the distance d is
Evaluating the gFREM (generalized PREM) for the special case of identical photon detection rates Λ_{1}(τ)=Λ_{2}(τ)=Λ_{0}, τ ≥ t _{0}, and assuming that q _{1}, q _{2} are given by the Airy profile, we obtain the FREM (PREM). See Supporting Text for more details and discussion.
Data Analysis.
The singlemolecule images were analyzed and processed by using matlab (Mathworks, Natick, MA). The singlemolecule pairs were identified by their doublestep photobleaching behavior. Before estimating the unknown parameters, the expected number of detected photons of the single molecules, the background photon count, and the parameter α were first independently determined from the singlemolecule images. The parameters were then estimated by the method of maximum likelihood (see Supporting Text).
Acknowledgments
This research was supported in part by National Institutes of Health Grants R01 GM071048 and R01 AI50747.
Footnotes
 ^{§}To whom correspondence should be addressed. Email: ober{at}utdallas.edu

Author contributions: S.R., E.S.W., and R.J.O. designed research; S.R. performed research; S.R., E.S.W., and R.J.O. contributed new reagents/analytic tools; S.R. analyzed data; and S.R., E.S.W., and R.J.O. wrote the paper.

Conflict of interest statement: No conflicts declared.

This paper was submitted directly (Track II) to the PNAS office.
 Abbreviations:
 FREM,
 fundamental resolution measure;
 gFREM,
 generalized FREM;
 PREM,
 practical resolution measure
 © 2006 by The National Academy of Sciences of the USA
References

↵
 Born M. ,
 Wolf E.

↵
 Inoue S. ,
 Spring K. R.

↵
 Gordon M. P. ,
 Ha T. ,
 Selvin P. R.

↵
 Qu X. ,
 Wu D. ,
 Mets L. ,
 Scherer N. F.
 ↵

↵
 Snyder D. L. ,
 Helstrom C. W. ,
 Lanterman A. D. ,
 Faisal M. ,
 White R. L.
 ↵
 ↵
 ↵
 ↵
 ↵
 ↵

↵
 Tinnefeld P. ,
 Müller C. ,
 Sauer M.

↵
 Snyder D. L. ,
 Miller M. I.
 ↵

↵
 Hobbs P. C. D.
 ↵

↵
 Mashanov G. I. ,
 Tacon D. ,
 Peckham M. ,
 Molloy J. E.
 ↵
 ↵
 ↵
 ↵
 ↵

↵
 Zacks S.

↵
 Hell S. ,
 Stelzer E. H. K.

↵
 Klar T. A. ,
 Jakobs S. ,
 Dyba M. ,
 Egner A. ,
 Hell S. W.
 ↵
 ↵

↵
 Paar J. M. ,
 Harris N. T. ,
 Holowka D. ,
 Baird B.
 ↵
Citation Manager Formats
More Articles of This Classification
Physical Sciences
Engineering
Biological Sciences
Related Content
Cited by...
 Subdiffraction imaging on standard microscopes through photobleaching microscopy with nonlinear processing
 Nanometerscale mapping and singlemolecule detection with colorcoded nanoparticle probes
 New directions in singlemolecule imaging and analysis
 Macromolecularscale resolution in biological fluorescence microscopy
 Using photon statistics to boost microscopy resolution