1. Introduction

Single-molecule-fluorescence imaging has been a powerful tool in particle localization and tracking studies [1, 2, 3, 4]. In single-molecule imaging, the fluorescence intensity profile of a point light source is called a point spread function (PSF). While the PSF is described by an Airy function, it is, in practice, approximated by a Gaussian. A Gaussian fit to the PSF of a stationary single fluorophore has two fitting parameters: centroid and standard deviation (SD). The centroid is the center of the PSF, and the SD is its width. While the centroid determines the lateral position of the particle in the imaging plane, the standard deviation determines its axial position [5, 6, 7] and orientation [8, 9, 10], as well as the spatial resolution of the fluorescence imaging system [11].

Error analysis provides the precision for a physical measurement, and is essential for validation of the method used. While error analysis of single-molecule PSF centroid measurements has provided the precision for lateral localization measurements [1], which enabled differentiation of various biological mechanisms (such as the walking mechanisms of myosin V on actin [2]), PSF SD measurement error analysis will provide the precision in the following applications: (1) single-molecule axial position measurements, where the SD of a single molecule’s PSF increases with the defocusing distance [5, 6, 7]; (2) single-molecule orientation measurements, where at different molecular orientations with respect to the imaging plane, the molecule exhibits an elliptical-shaped Gaussian PSF with a SD that changes in both lateral directions depending on its orientation [8, 9, 10]; and (3) characterization of a single-molecule-fluorescence imaging system, where the measured SD of an imaged single fluorophore determines whether the imaging system is diffraction-limited [11].

In contrast to the precision of centroid measurements which has been extensively investigated and applied to many systems [1, 2, 12], the precision of SD measurements of single-molecule PSFs has not been evaluated. Here we report SD measurement error studies of immobile-single-molecule PSFs using analytical calculation, numerical integration, simulation, and experimental measurements. As with centroid analysis, the precision of SD measurements is affected by the experimental settings of a finite number of photons per PSF N, the standard deviation of the background noise σ_b, and the camera’s finite pixel size a. We have obtained an analytical expression for the PSF SD measurement error as a function of these parameters. Our SD measurements have achieved nanometer resolution for a wide range of experimental conditions. This expression for the SD measurement error will provide confidence in determining a particle’s axial position and molecular orientation from measurements using a single-molecule imaging system of known resolution.

2. Theory

2.1. Formulating SD measurement error, ∆s, by χ² minimization

In this article, we derive the SD error for a PSF, which is a collection of N photons, from a common distribution emitted by a point light source. We include the additional experimental effects of photon count fluctuation per PSF, background noise, and camera pixelation in our study.

We utilized the method developed by Bobroff [13] and subsequently used for centroid error analysis by Thompson, Larson, and Webb [1] to derive the error associated with SD measurements of single fluorophores. The approach uses Chi-square statistics to estimate the error associated with fitting of experimental data to expected theoretical values. In order to maintain consistency in notation for single-molecule tracking studies, we will retain many of the same notations used in Ref. [1]. Below, we derive the analytical solution to the PSF SD error as a function of N, a, and σ_b (b in prior studies) beginning with one dimension and extending to two dimensions.

In 1D least squares fitting of the intensity profile of an immobile single fluorophore, χ ²(s) is proportional to the sum of squared errors between the observed photon count at pixel i, y_i, and the expected photon count N_i(x,s), of a PSF. Here x and s are the measured position and SD of the PSF, respectively, while x ₀ and s ₀ are the true location and the theoretical SD of the molecule:

where σ_i,photon is the expected photon count uncertainty at pixel i without accounting for photon-to-camera count conversion (described in the following section). In this article, we emphasize the SD error and assume that the location measurement errors are negligible, i.e. x = x ₀. For simplicity, N_i(x ₀,s) is denoted as N_i in this article unless otherwise specified.

There are two sources for σ_i,photon at pixel i: one is the Poisson-distributed photon shot noise of the PSF where the variance is the mean expected photon count of the pixel, N_i, and the other is the SD of the background noise, σ_b, expressed in photons. The variances of the two sources add to yield

The deviation of s from s ₀, ∆s = s-s ₀, is obtained by setting dχ ²(s)/ds to 0, expanding N_i about s ₀, and keeping the first order term in ∆s. Appendix A shows the detailed derivation of ∆s from dχ ²(s)/ds = 0. The mean squared value of ∆s is

The root mean square of ∆s, ∆s_rms, is the PSF SD error that we calculate in this article.

2.2. Modifying σ_i,photon to include camera count conversion effects

When an EMCCD (Electron Multiplying Charge Coupled Device) camera is used in imaging single fluorescent molecules, the detected pixel reading is in camera counts. In converting from camera counts to photon counts, an additional variance in σ_i,photon appears. Below we derive the uncertainty in photon counts,σ_i, to use in place of σ_i,photon in Eq. (3) for experiments where EMCCD camera count conversions are involved.

An EMCCD camera amplifies the detected photons by converting each photon to a distribution of photoelectrons through many multiplication stages. At the final stage, one photon yields a distribution of camera counts (equivalent to the last stage photoelectron counts) with a distribution function f(n ^*) [14],

where n ^* is the camera counts in the distribution and M is the photon multiplication factor of the camera. Here we use^* to denote camera counts in order to differentiate from photon counts. The n ^* distribution has a mean of M and a variance of M ².

At pixel i, the PSF photon count distribution is described by a Poisson distribution with the variance being equal to the mean. Each photon at the pixel contributes two terms to the pixel’s camera count variance: the mean photon shot noise variance M ² (variance of a single photon count, which is one, multiplied by the square of the multiplication factor), and the photon-to-camera count conversion variance M ². The total camera count variance contributed by one photon is 2M ²; therefore, a mean of N_i photons yields a camera count variance of 2N_iM ². This variance agrees with the expression in Ref. [15] where the variance in camera counts σ_out,camera ², is related to the variance in photon counts σ_in,photon ², by an excess noise factor F ²,

for EMCCD cameras with a large number of multiplication stages.

Fluorescence from buffer, diffusing molecules in the solution, and camera counts from electronic readout and thermal noise constitutes the total background photon count at pixel i, with a variance of σ_b ² and a mean of 〈b〉. The total background variance in camera counts is the sum of the background count variance σ_b ² M ², and the variance introduced by the average number of background photons, 〈b〉, each with a variance of M ²: (σ_b ² + 〈b〉)M ².

Summing the PSF and the background contributions, the total camera count variance at pixel i is

When expressed in photon counts,

Revising Eq. (3) with the modified σ_i we have

2.3. Expressing ∆s in photon counts

To evaluate Eq. (8) in 1D, we use a normalized Gaussian distribution

where we set the location of the PSF to be at x ₀ = 0 for simplicity and without loss in generality. We approximate the pixel summation in Eq. (8) by an integral going from negative to positive infinity, and we estimate 〈(∆s)²〉 at the two extrema of σ_i ²: the high photon count regime where σ_b ² + 〈b〉 can be neglected, and the high background noise regime where 2N_i, can be neglected. In the high photon count regime,

and in the high background noise regime,

An alternative derivation of Eq. (10) is presented in Ref. [16], although the photon-to-camera count conversion variance was not included and thus 〈(∆s)²〉 = s ₀ ²/2N. The total 1D 〈(∆s)²〉 is the sum of Eqs. (10) and (11) (without the pixelation effect discussed below)

The method of approximating 〈(∆s)²〉 by summing these results for both extema of σ_i ² is validated by numerical calculation results shown in Fig. 2, and is in accordance with Ref. [1]. We now calculate the effect of camera pixelation on 〈(∆s)²〉. Each photon in a PSF is associated with two variances with respect to the centroid. One is the mean variance of the PSF, s ₀ ², and the other is due to the fact that each photon is further binned into a pixel that has an intensity profile described by a uniform distribution with a width corresponding to the pixel size a. The variance of this distribution is a ²/12. Thus, the total variance of a photon due to pixelation is the sum of the two,

Under experimental conditions, the measured s should be (s ₀ ² + a ²/12)^1/2 and for theoretical formulations, the expected SD of a PSF should include the pixelation effect. We have verified that s ₀ ^†2 increases with a according to Eq. (13) by simulation. Plugging Eq. (13) into Eq. (12) we have for 1D

Extending the 1D 〈(∆s)²〉 calculation to 2D where s_x,y, which for the remainder of this article, represents the SD in either the x or y direction of the imaging plane, and s _0x and s _0y are the theoretical SD values in the x and y directions, respectively,

The derivation of Eq. (15) is provided in Appendix B.

A more accurate estimation of 〈(∆s_x,y)²〉 can be obtained by numerically integrating Eq. (8), incorporating the transition region between the high photon count and the high background noise regimes. The numerical integration results are shown in Fig. 2 to be consistently higher than the analytical calculation results by ≈ 15%.

3. Methods

3.1. Experimental setup

Single-molecule imaging was performed using a Nikon Eclipse TE2000-S inverted microscope (Nikon, Melville, NY) attached to an iXon back-illuminated EMCCD camera (DV897ECS-BV, Andor Technology, Belfast, Northern Ireland). Prism-type Total Internal Reflection Fluorescence (TIRF) microscopy was used to excite the fluorophores with a linearly polarized 532 nm laser line (I70C-SPECTRUM Argon/Krypton laser, Coherent Inc., Santa Clara, CA) focused to a 40 μm × 20 μm region on fused-silica surfaces (Hoya Corporation USA, San Jose, CA). The incident angle at the fused-silica water interface was 68° – 71° with respect to the normal. The laser was pulsed with illumination intervals between 1 ms and 500 ms and excitation intensity between 0.3 kW/cm² and 2.6 kW/cm². By combining laser power and pulsing interval variations we obtained 50 to 3000 photons per PSF. A Nikon 100X TIRF objective (Nikon, 1.49 NA, oil immersion) was used in combination with a 2X expansion lens, giving a pixel size of 79 nm.

At focus, the PSF image generated by a point light source with a mean emission wavelength of 580 nm and symmetric polarization has a full width at half-maximum (FWHM) of ≈ λ/2NA = 580 nm/2.9 ≈ 200 nm and theoretical s ₀ = FWHM/2.35 ≈ 85 nm. Including the pixelation effect [Eq. (13)], the measured PSF SD s _0x,0y ^†, for our imaging system should be 88 nm. Due to random fluctuations in the emission polarization direction of streptavidin-Cy3 molecules attached to surfaces [8] and variations in focus between each measurement, we observed a range of s _0x,0y ^† values from 90 nm to 140 nm.

Single streptavidin-Cy3 molecules (SA1010, Invitrogen, Carlsbad, CA; 530/10 excitation, 580/60 emission) were immobilized on fused-silica surfaces by depositing 6 μl of 0.04 nM streptavidin-Cy3 powder dissolved in 0.5X TBE buffer (45 mM Tris, 45 mM Boric Acid, 1 mM EDTA, pH 7.0). A coverslip flattened the droplet and its edges were sealed with nail polish. The fused-silica chips were cleaned using oxygen plasma before use. We inspected for possible surface fluorescence contaminations by imaging the TBE buffer alone; no impurities were found on either the fused-silica surface or in the buffer. The immobilization of the adsorbed molecules was verified by centroid vs time measurements.

3.2. Data acquisition and selection

Typical movies were obtained by synchronizing the onset of camera exposure with laser illumination for different intervals. The gain levels of the camera were adjusted such that none of the pixels of a PSF reached the saturation level of the camera. For the initial step, streptavidin-Cy3 monomers were first selected in ImageJ (NIH, Bethesda, MD) by examining the fluorescence time traces of the molecules for a single bleaching step [17]. For a selected monomer, the intensity values for 25 × 25 pixels centered at the molecule were recorded. The center 15 × 15 pixels of the PSF were used for 2D Gaussian fitting with peripheral pixels used for background analysis.

The intensity values of the selected molecules were first converted to photon counts (see the following section) and then fitted to the following 2D Gaussian function using a least squares curve fitting algorithm (lsqcurvefit) provided by MATLAB (The Mathworks, Natick, MA):

where f ₀ was the amplitude and 〈b〉 was the mean background value. A background pixel’s total count is the sum of the floor, electronic readout noise, and background fluorescence counts. For the 〈b〉 in this article, the floor value, determined by the lowest background pixel value, has already been subtracted. With this fitting, the PSF’s SD values in both the x and y directions, its measured location (x ₀,y ₀), and the image’s mean background value were obtained.

The selected streptavidin-Cy3 monomers were further characterized to satisfy the following conditions used for SD error analysis. (1) No stage drift detected by using centroid vs time measurements. Stage drift introduces additional blur to each single-molecule PSF and thus affects the measured SD values. (2) A minimum of 75 valid PSF images, each with a photon count N that fluctuated less than 20% from the experimental mean 〈N〉, of the monomer. The PSF N count restriction is necessary for precise SD error analysis at N by using a statistically sufficient number of PSFs with consistent N. (3) PSFs with signal-to-noise ratios (I ₀/√I ₀ + σ_b ²) larger than 2.5, where I ₀ is the peak PSF photon count (total photon count minus 〈b〉) and σ_b ² is the background variance in photons. (4) Mean 〈s_x〉 and 〈s_y〉 obtained by Gaussian fitting of the s_x and s_y distributions of all valid images did not differ by more than 10 nm, or ±5% of the mean SD value to minimize polarization effects of Cy3. (5) The mean SD values 〈s_x,y〉 were between 95 nm and 135 nm to minimize defocusing effects. These constraints on s_x and s_y are necessary for obtaining the expression for ∆s_rms, as a function of N, with minimal variations in the other parameters.

To convert from a pixel’s camera count to photons, the camera count value was divided by M. In order to obtain M for each experimental setting, the center nine pixel values of the PSF were evaluated if the molecule’s average signal-to-noise ratio was greater than 3. According to Eq. (7),

where 〈N_i ^*〉 and σ_i ^* are the Gaussian fitted mean and standard deviation of the measured camera count distribution of pixel i, respectively. Here 〈N_i ^*〉 is the mean camera count that includes background fluorescence and electronic noise counts.

3.3. PSF and background simulations

Single-fluorescent-molecule PSFs were generated using the Gaussian random number generator in MATLAB. For Fig. 2, the s _0x,0y of each simulated PSF was determined by the experimental means 〈s_x,y〉. The observed fluctuation in the number of photons N, was incorporated. The generated photons of each PSF were binned into 15 × 15 pixels with a pixel size of 79 nm. Then each photon count in a pixel was converted into camera count using Eq. 4 with a M value of one. Random background photons at each pixel were generated using the corresponding experimental background distribution function. Although the exact experimental background distributions were used for the simulations, the numerical integrations and analytical calculations were computed using the theoretical variance and the mean of all background counts, σ_b ² and 〈b〉, respectively, rather than their fitted values. The background counts are primarily drawn from two types of distributions: a full Gaussian with a high mean or a truncated Gaussian with a low mean, depending on the background fluorescence level of each specific experiment. The final simulated PSFs with background noise were fitted to a 2D Gaussian [Eq. (16)] to obtain the centroid and SD values of the PSF.

For each simulated ∆s_x,y,rms data point, 1000 iterations (2000 iterations for Fig. 3) were performed and the Gaussian fitted SDs of the s_x,y distributions were the simulated ∆s_x,y,rms results.

4. Results

We report our study of 2D ∆s_x,rms using four different methods: (1) experimental measurements, (2) simulations, (3) numerical integrations of Eq. (8), and (4) analytical calculations using Eq. (15).

Figure 1A shows a set of single streptavidin-Cy3 molecule images with an increasing number of detected photons N. These molecules have similar mean SD 〈s_x〉 values of 110 nm, 111 nm, and 107 nm, respectively. In order to demonstrate the decreasing SD error with increasing N, each representative image was chosen such that the 2D SD value was the sum of the mean SD 〈s_x〉, and one standard deviation of the molecule’s s_x distribution ∆s_x,rms (SD_image = 〈s_x〉 + ∆s_x,rms). To clearly illustrate the change in the SD error, which is measured as the PSF SD minus 〈s_x〉, the 1D intensity profiles of the PSFs are plotted in Fig. 1B as opposed to their 2D intensity profiles for clarity. The 1D intensity values were obtained by averaging transverse pixel intensity values of the PSF at each longitudinal pixel i. The measured 2D SD_image values deviate from their respective means, 〈s_x〉 values, by 10.3 nm, 7.2 nm, and 2.7 nm. As expected, the 2D SD error decreases with increasing N.

 

Fig. 1. (A) Representative images with increasing N of 151, 393, and 1891 photons of single streptavidin-Cy3 molecules. (B) 1D intensity profiles (circles) of the molecules in (A) and their Gaussian fits (lines). The respective 1D SD values are 195.4 nm, 140.5 nm, and 110.9 nm. The scale bar is 500 nm.

Download Full Size | PDF

Figure 2 shows ∆s_x,rms obtained by using experimental measurements, simulations, numerical integrations, and analytical calculations. Each experimental ∆s_x,rms data point is the standard deviation of the s_x distribution for a single streptavidin-Cy3 monomer. A simulation was performed for each experimental data point. The parameters were based upon experimental results including fluctuations in a PSF’s total detected photons, background distribution, and the s _0x,0y values determined by the mean experimental 〈s_x,y〉 after subtracting for the pixelation effect [Eq. (13)]. The finite bandwidth of the emission filter was also taken into consideration by simulating each photon as being drawn from a PSF whose width is varied according to a Gaussian distribution centered about s _0x,0y (with SD of 2 nm). Numerical integrations and analytical calculations used the same 〈N〉, s _0x,0y, σ_b, and 〈b〉 as those in the corresponding experimental data point. For all N, the numerically integrated ∆s_x,rms results are ≈ 15% higher than the theoretical results while the experimental results are ≈ 57% higher. The causes of these discrepancies are discussed in the following section. The simulations agree well with the experimental results.

The above results are for our pixel size of 79 nm. For different experimental settings the pixel size will vary and affect ∆s_x,rms. Figure 3 shows ∆s_x,rms vs a/s ₀ studied by simulations using s _0x = s _0y = s ₀ = 120 nm, N = 500 photons, σ_b = 1 photon, and 〈b〉 = 4 photons. The generated photons of each PSF were binned into 19 × 19 pixels and subsequently converted into camera counts following the same procedure described above for Fig. 2. As a/s ₀ increases, there is an initial decline in ∆s_x,rms until rising at a/s ₀ ≈ 0.73. Beyond a/s ₀ ≈ 0.73, ∆s_x,rms increases slightly and then continues the decline again at a/s ₀ ≈ 1.18. This decline after a/s ₀ ≈ 1.18 disagrees with theory, which suggests an increase in ∆s_x,rms beyond the theoretical minimum of ${\left(a/s_0\right)}^4=\frac{144}{\frac{9N}{4\pi \left({\sigma }_b^2+〈b〉\right)}+1}$ at a/s₀ = 1.18 (vertical dashed line). The overall decreasing ∆s_x,rms trend after the theoretical minimum occurs because when the pixel size increases, the measured PSF SD is increasingly affected by the width of the pixel and approaches the SD of the pixel; thus, variations among measured SD values decrease. Eventually, at sufficiently large pixel sizes where the whole PSF is contained within one pixel, the measured SD will be the SD of the pixel, inferred by the top-hat distribution function, and the measured SD error will be zero. The analytical calculation does not take this large pixelation effect into consideration; consequently, these results and those of the simulations begin to rapidly diverge.

 

Fig. 2. Comparing ∆s_x,rms vs N obtained by using four different methods: experimental measurements (solid squares), simulations (circles), numerical integrations (crosses), and analytical calculations (dashed line). Each experimental ∆s_x,rms data point is the SD from the Gaussian fit to the s_x distribution of a single streptavidin-Cy3 monomer. For each data point, its experimental N and background distributions were used for simulation, and its experimental 〈N〉, 〈s_x,y〉, σ_b, and 〈b〉 values were used for the numerical integrations and analytical calculations. The experimental data are on average 57% higher than the analytical calculation data.

Download Full Size | PDF

 

Fig. 3. ∆s_x,rms vs a/s ₀ studied by simulations. In these simulations, there were no fluctuations in N or s ₀. The vertical dashed line at a/s ₀ = 1.18 is where the theoretical ∆s_x,rms minimum occurs, determined by differentiating Eq. (15) with respect to a.

Download Full Size | PDF

The simulated local ∆s_x,rms minimum occurs at a/s ₀ = 0.73, rather than at the theoretical minimum of a/s ₀ = 1.18 due to the pixel size effect described above. Our experimental settings of a = 79 nm and s ₀ = 120 nm yield a/s ₀ = 0.66 and is close to the simulated ∆s_x,rms minimum. We have also performed additional simulations using different parameter sets where the theoretical minimum always preceded the continued decline in ∆s_x,rms. According to Fig. 3 and our other simulations, a good a/s ₀ range for future studies should be between ≈ 0.5 and 1, as is usually the case. Future ∆s_x,rms studies using different pixel sizes should take this discrepancy into account.

Note that the simulated ∆s_x,rms minimum at a/s ₀ = 0.73 is different from the theoretical ∆x_rms minimum at a/s ₀ = 0.88 described in Ref. [1], and our theoretical ∆s_x,rms minimum at a/s ₀ = 1.10 calculated from Eq. 18, using our set of parameter values. Future studies should take this difference into consideration by selecting an optimal pixel size.

5. Discussion and Extensions

Here we discuss four issues: (1) causes for discrepancies between results obtained using different methods; (2) modifications to the centroid measurement error developed by Thompson, Larson, and Webb [1] to include the EMCCD camera photon conversion effects; (3) relation between SD error and the error of the measured quantities associated with each of the aforementioned applications; and (4) methods to determine the SD error ∆s_x,y,rms, for dimeric fluorophores and mobile molecules in future studies.

5.1. Causes for discrepancies

Numerical integration results are consistently higher than the analytical results by 15%, while simulation (and experimental) results are higher than analytical results by 57% for all N. There are four reasons for these discrepancies:

In summary, the analytical calculation of the SD measurement error expressed in Eq. (15) is a reasonable approximation for a large range of experimental parameters. When the 57% difference is corrected for, the expression is in excellent agreement with our experimental results. Future studies using this formula should be aware of the limitations and be sure to include a 57% difference from underestimation of the true error for similar a/s ₀ values.

5.2. Modifications to centroid error analysis

The PSF centroid error expression developed by Thompson, Larson, and Webb [1] did not take the photon-to-camera count conversion variance into consideration. Additionally, the theoretical standard deviation s ₀, should be modified to include the pixelation effect (s ₀ ² + a ²/12)^1/2, with respect to both directions. We have modified the PSF centroid measurement error to be

This theoretical expression for the centroid measurement error underestimates the experimental results by 42%.

5.3. Interpreting ∆s_x,y,rms in SD measurement applications

With regards to the three applications of SD measurements discussed in this article, the SD measurement error of a single image can be translated into the precisions associated with each of the application’s measured quantities. For the direct translation, SD error is the uncertainty of an imaging system’s measured resolution; for the indirect translations, the precisions for axial localization and molecular orientation measurements can be expressed by SD error.

For characterization of an imaging system’s spatial resolution, if the system is diffraction-limited, the SD of the imaging system should fall between the measured SD of the PSF from a visible point light source ± the SD error [11, 18]. Thus, SD measurement error directly provides the precision for quantifying an imaging system’s resolution.

For axial localization studies, it has been previously shown that the SD of the PSF of a molecule located at a distance z away from the focal plane can be expressed as [5, 6, 7]

where D ≈ 400 nm is the imaging depth of a typical single-molecule imaging system. Consequently, by error propagation, the precision in the SD measurement of a single image, ∆s_rms, can be used to determine the localization error associated with the molecule’s axial position, ∆z:

For molecular orientation studies, the polarized PSF, for a range of orientations, has an elliptical intensity profile that can be fit by a 2D Gaussian with different standard deviations in the x and y directions [8, 9, 10]. When an expression relating s_x and s_y to the orientation is developed, the error in measuring s_x and s_y, once again by error propagation, can be used to calculate an error associated with the reported orientation of the molecule.

5.4. ∆s_x,y,rms calculation for future SD measurement applications

In addition to stationary single molecules, SD measurements can be used in future applications to investigate molecules such as stationary dimers or moving fluorophores. We are currently exploring these two areas of interest: (1) two sub-diffraction limit, separated molecules labeled with identical fluorophores that exhibits a combined PSF with a SD that increases with their separation [19]; (2) a moving molecule (i.e. a freely-diffusing fluorophore) which produces a blurred image given a finite exposure time, whereby the measured SD of the resulting intensity profile can be used to study the dynamic properties of the molecule, particularly its diffusion coefficient [20]. With modification, the method for estimating the SD error of stationary single molecules in this article can be extended to these two cases. For these studies, the N_i distribution function at each pixel may be different from the Gaussian assumption for stationary molecules in Eq. (9). A new N_i distribution function for each specific case can be obtained and a new σ_i ² formula [Eq. (7)] can be derived. Using the new N_i distribution function and σ_i ², the SD error for these cases can be obtained following the same procedure outlined in the theory (Sec. 2).

6. Conclusion

In this article we report the precision analysis for SD measurements of single-fluorescent-molecule intensity profiles. Our analytical expression of the PSF SD error allows for proper quantification of the precision associated with determination of the imaging system’s resolution and both axial localization and molecular orientation measurements of single molecules. Furthermore, we propose additional studies to characterize multiple fluorophores and examine the diffusive properties of mobile molecules by evaluating the measured SDs of their corresponding intensity profile to known precision. When our theoretical framework is extended to these studies, SD analysis will be advanced into a powerful tool for single-molecule-fluorescence imaging studies.