## Power spectral analysis for optical trap stiffness calibration from high-speed camera position detection with limited bandwidth

Optics Express, Vol. 18, Issue 8, pp. 7670-7677 (2010)

http://dx.doi.org/10.1364/OE.18.007670

Acrobat PDF (1826 KB)

### Abstract

The use of camera imaging enables trap calibration for multiple particles simultaneously. For stiff traps, however, blur from image integration time affects the detected particle positions significantly. In this paper we use power spectral analysis to calibrate stiff optical traps, taking the effects of blur, aliasing and position detection error into account, as put forward by Wong and Halvorsen [Opt. Express **14**, 12517 (2006)]. We find agreement with simultaneously obtained photodiode data and the expected relation of corner frequency *f _{c}
* with laser power, up to

*f*= 3.6 kHz for a Nyquist frequency of 1.25 kHz. Spectral analysis enables easy identification of the contribution of noise. We demonstrate the utility of our approach with simultaneous calibration of multiple holographic optical traps.

_{c}© 2010 OSA

## 1. Introduction

*f*>

_{c}*f*. To do so, we fit the power spectra taking blur, aliasing and position detection error into account, to obtain

_{Nyq}*f*for a trapped particle at different laser powers. We compare the results to those from power spectra obtained simultaneously using a position sensitive photodiode (PSD) and find good agreement. To demonstrate the use of this technique to calibrate traps that display spatial modulations and drift, we apply it to calibrate multiple holographic optical traps (HOTs) simultaneously, finding the expected linear relation with laser power. Our work demonstrates, for the first time, the viability of using camera power spectral analysis for calibration of very stiff optical traps, with

_{c}*f*up to at least 3 ×

_{c}*f*.

_{Nyq}## 2. The effect of blur, aliasing, and detection error on measured power spectra

*S*of the

_{xx}*x*position of a particle trapped in a harmonic potential well as a function of the frequency

*f*is described by a Lorentzian [5

5. K. Berg-Sørensen and H. Flyvbjerg, “Power spectrum analysis for optical tweezers,” Rev. Sci. Instrum. **75**(3), 594–612 (2004). [CrossRef]

*k*is the Boltzmann constant and

_{B}*T*is the absolute temperature. In Fig. 1 this function is plotted for three different corner frequencies

*f*(solid black lines).

_{c}*W*,

*e.g.*camera imaging, the blur due to this integration time affects the acquired spectrum and has to be taken into account [2

2. W. P. Wong and K. Halvorsen, “The effect of integration time on fluctuation measurements: calibrating an optical trap in the presence of motion blur,” Opt. Express **14**(25), 12517–12531 (2006), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-14-25-12517. [CrossRef] [PubMed]

*S*. With the Nyquist frequency equal to half the sampling frequency,

*f*= 0.5 ×

_{Nyq}*f*, we find:

_{s}*f*with respect to

_{c}*f*. For the blurred spectrum, however, the effect of aliasing is limited for

_{Nyq}*f*= 1/

_{s}*W*(red dash-dotted lines), but increases for 1/W >

*f*(solid blue lines: 1/

_{s}*W*= 8 ×

*f*). The relation among the parameters

_{Nyq}*f*,

_{Nyq}*f*and

_{c}*W*affects the extent to which an acquired power spectrum differs from a pure Lorentzian. As can be seen from Fig. 1, the greater

*f*and 1/

_{c}*W*are compared to

*f*, the larger the deviation of the spectrum from Eq. (1), suggesting that it becomes more important to use the adjusted Lorentzian (Eq. (2)) to obtain the correct

_{Nyq}*f*.

_{c}*f*, an error with variance

*ε*

^{2}will introduce a noise level of

*ε*

^{2}/2

*f*to the power spectrum. Taking the effects of blur, aliasing and tracking error into account, the measured double-sided power spectral density is described by:

_{Nyq}## 3. Experimental setup and fitting of measured power spectra

3. A. van der Horst and N. R. Forde, “Calibration of dynamic holographic optical tweezers for force measurements on biomaterials,” Opt. Express **16**(25), 20987–21003 (2008), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-16-25-20987. [CrossRef] [PubMed]

^{2}) to image polystyrene spheres (Spherotech) trapped in water. Particles were trapped in the center of our sample cell (~40 μm thick), which consisted of two cover slips closed off with candle wax. Due to our current illumination, we limited the acquisition rate to 2500 frames/s. Image correlation was used for particle position detection [3

3. A. van der Horst and N. R. Forde, “Calibration of dynamic holographic optical tweezers for force measurements on biomaterials,” Opt. Express **16**(25), 20987–21003 (2008), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-16-25-20987. [CrossRef] [PubMed]

*f*, distinctive noise peaks are eliminated. For the camera data, we average over 5 to 15 power spectra and fit to the measured frequency range (

_{c}*f*= 0,

*f*= 1250 Hz). For the position sensitive photodiode (

_{Nyq}*f*= 16384 Hz) we average over >100 spectra. We fit only up to

_{Nyq}*f*= 5 kHz, due to the decreased sensitivity of our silicon-based PSD (OSI Optoelectronics, DL-10) above ~5 kHz for 1064 nm light [7

7. E. J. G. Peterman, M. A. van Dijk, L. C. Kapitein, and C. F. Schmidt, “Extending the bandwidth of optical-tweezers interferometry,” Rev. Sci. Instrum. **74**(7), 3246–3249 (2003). [CrossRef]

*f*>10 kHz. As a consequence, the noise level due to the detection error,

_{Nyq}*ε*

^{2}/2

*f*, is usually not significant, even for high trap stiffness. In addition, the detection integration time is on the order of (sub)μs, making the effect of blur negligible. Because of the effective low-pass filtering by our photodiode, we do not include aliasing and therefore fit our single-sided PSD power spectra, measured in V

_{Nyq}^{2}/Hz, with:in which

*β*is an additional prefactor to account for the (unknown) conversion from Volts to nm. From Fig. 1(a) it can also be seen that for

*f*<<

_{c}*f*, the corner frequency can accurately be obtained from the shape of the spectrum, even when

_{Nyq}*γ*or

*β*are not known. When fitting PSD power spectra using Eq. (5),

*f*and

_{c}*γ*/

*β*are free fitting parameters for each spectrum.

*W*is on the order of ms and blur has to be taken into account. In addition, with a bandwidth of typically no more than several kHz, both the noise level and aliasing cannot

*a priori*be neglected. Thus, we use

*f*,

_{c}*ε*

^{2}and/or

*γ*as free fitting parameters (see below) to fit the camera power spectra with the single-sided version of Eq. (4). Even though

*γ*can be calculated, using it as a fitting parameter allows us to account for possible errors in the estimate of particle radius and medium viscosity. The conversion from pixels to μm is constant (and known) for all of our measurements, so no free conversion parameter like

*β*is necessary. We stop the summation in Eq. (4) at

*n*= ± 3, which we found to be sufficient to account for the effect of blur.

*f*>>

_{c}*f*(Fig. 1(c)), the power spectrum shows so little distinct curvature that

_{Nyq}*f*cannot easily be obtained from the shape of the spectrum. In addition, for higher trap stiffness, the motion of the particle in the trap decreases and the apparent motion due to the detection error

_{c}*ε*

^{2}becomes increasingly more significant. Therefore, at these higher trap stiffnesses,

*f*can only be determined when both

_{c}*γ*and

*ε*

^{2}are known. The viscous drag coefficient

*γ*can be obtained from a power spectrum for which

*f*<<

_{c}*f*(Fig. 1(a)). In this regime, the relative contribution of

_{Nyq}*ε*

^{2}is very limited. The plateau region of the power spectrum at low frequencies and the Brownian motion at high frequencies are clearly visible and fitting will yield both

*f*and

_{c}*γ*. Alternatively,

*γ*can be obtained by fitting to only the high-frequency region, since for

*f*>4 ×

*f*, the tail of the power spectrum is practically independent of

_{c}*f*. It is clear that the range of

_{c}*f*over which this calibration method gives valid results will largely be determined by the accuracy to which γ and

_{c}*ε*

^{2}are known.

*f*can be determined by least-squares analysis of power spectra in which data points are statistically independent and normally distributed, criteria that are relatively easy to meet for PSD-derived data [5

_{c}5. K. Berg-Sørensen and H. Flyvbjerg, “Power spectrum analysis for optical tweezers,” Rev. Sci. Instrum. **75**(3), 594–612 (2004). [CrossRef]

*f*means that measurements need to be made over orders of magnitude longer time to achieve a comparable number of data points. Instead, we average our camera data (25,000 to 75,000 points for each laser power) over 5 to 15 spectra to reduce noise, but do not apply additional blocking.

_{Nyq}## 4. Experimental results and discussion

*ε*

_{x}^{2}/

*f*= 7.2 × 10

_{Nyq}^{−4}nm

^{2}/Hz (

*ε*

_{x}^{2}= 0.90 nm

^{2}; Fig. 2(b)) and

*ε*

_{y}^{2}/

*f*= 14.5 × 10

_{Nyq}^{−4}nm

^{2}/Hz (

*ε*

_{y}^{2}= 1.8 nm

^{2}). Large peaks in the power spectra, small or absent in power spectra of trapped particles, are from mechanical noise coupled in through the sample surface. As

*ε*

^{2}depends on the illumination and on

*W*, the images were taken under the same conditions as used in the following experiment (

*W*= 0.4 ms).

*f*and

_{c}*γ*/

*β*as independent fitting parameters for each spectrum. The normalized PSD power spectra in

*x*and their fits (red lines) for three laser powers are shown in Fig. 3(a) .

*f*>4 ×

*f*) region of a power spectrum depends very little on

_{c}*f*. Therefore, to obtain

_{c}*γ*, we first fit this region of spectra at low laser power (6.1 mW) with the error fixed at

*ε*

_{x}^{2}= 0.90 nm

^{2}and

*ε*

_{y}^{2}= 1.8 nm

^{2}. (The parameter

*f*was fixed at an initially estimated value obtained from fitting to the whole frequency range.) We found

_{c}*γ*= 1.93 × 10

_{x}^{−8}kg/s and

*γ*= 1.96 × 10

_{y}^{−8}kg/s, which are close to the expected value at 20 °C of 1.98 × 10

^{−8}kg/s. Next, we fixed

*γ*at these measured values and fit spectra at all laser powers with

*f*as the only free fitting parameter. Figure 3(b) shows three camera power spectra and their fits. In Fig. 3(c), a camera power spectrum for 105 mW and its fit (red) are shown, together with a pure Lorentzian fit (Eq. (5); blue), with

_{c}*f*and

_{c}*γ*as free fitting parameters. Both curves seem to fit the data well, yet the corner frequencies found differ greatly (

*f*= 774 and 583 Hz, respectively), demonstrating the significant differences in obtained corner frequencies that result from using Eq. (4) vs. Equation (5) to fit power spectra of stiffer traps.

_{c}*f*from camera images and from PSD data for all measured laser powers. We find agreement between PSD and camera results (within 10%), up to

_{c}*f*2.9 ×

_{c}=*f*(at 416 mW). At this high trap stiffness (

_{Nyq}*κ*= 432 pN/µm), the height of the plateau is 23.7 × 10

_{x}^{−4}nm

^{2}/Hz. Therefore, 30% of the particle’s

*apparent*motion is due to the detection error

*ε*

^{2}, making it quite remarkable that an accurate

*f*can be obtained. Also plotted are the results from pure Lorentzian fits to the camera spectra. As can be seen, for higher laser power,

_{c}*f*is increasingly underestimated when using Eq. (5) rather than Eq. (4) to fit power spectra from the camera.

_{c}*γ*to be constant for all laser powers. However, due to the absorption of laser light by the solvent and resultant heating, the viscosity is expected to decrease and therefore

*f*to increase with laser power [8

_{c}8. E. J. G. Peterman, F. Gittes, and C. F. Schmidt, “Laser-induced heating in optical traps,” Biophys. J. **84**(2), 1308–1316 (2003). [CrossRef] [PubMed]

9. H. Mao, J. R. Arias-Gonzalez, S. B. Smith, I. Tinoco Jr, and C. Bustamante, “Temperature control methods in a laser tweezers system,” Biophys. J. **89**(2), 1308–1316 (2005). [CrossRef] [PubMed]

*κ*is independent of temperature,

*f*is expected not to be [9

_{c}9. H. Mao, J. R. Arias-Gonzalez, S. B. Smith, I. Tinoco Jr, and C. Bustamante, “Temperature control methods in a laser tweezers system,” Biophys. J. **89**(2), 1308–1316 (2005). [CrossRef] [PubMed]

8. E. J. G. Peterman, F. Gittes, and C. F. Schmidt, “Laser-induced heating in optical traps,” Biophys. J. **84**(2), 1308–1316 (2003). [CrossRef] [PubMed]

*ΔT*<3K, resulting in a decrease in

*γ*of <7%. Because

*γ*/

*β*is a free fitting parameter, able to capture temperature changes, our fits should correctly return

*f*, which we expect to increase from heating effects by <7%. For the camera data, we expect to observe a smaller increase in fit values of

_{c}*f*with increasing power, because the use of a fixed

_{c}*γ*obtained at low laser power results in underestimating the true

*f*. The lines in Fig. 4 are linear fits to the

_{c}*x*and

*y*camera results at low laser power (up to 53 mW) for which we expect the possible effect of heating to be minimal. We see that for laser powers >100mW,

*f*lies above this linear fit for both

_{c}*x*and

*y*. Further investigation is needed to confirm whether this is indeed a result of a decrease in

*γ*due to heating.

*ε*

^{2}can be found by fitting a spectrum whose shape is predominantly determined by the corner frequency

*f*, but to which the detection error

_{c}*ε*

^{2}nonetheless has a significant contribution. As mentioned before,

*γ*can be found from the high-frequency region of a power spectrum at relatively low laser power. Now, by fitting two power spectra iteratively, one at low and one at high laser power, we can obtain both

*ε*

^{2}and

*γ*. Because the power spectra at the highest laser power show very little curvature (416 mW, see Fig. 3(b)), making it hard to independently find

*f*, we used the spectra for 6.1 and 216 mW and found

_{c}*ε*

_{x}^{2}= 0.83 nm

^{2}(

*γ*= 1.93 × 10

_{x}^{−8}kg/s) and

*ε*

_{y}^{2}= 1.0 nm

^{2}(

*γ*= 1.95 × 10

_{x}^{−8}kg/s), which are within a factor of two of the results for the stuck particle.

*i.e.*, iteratively fitting two power spectra to determine

*ε*

^{2}and

*γ*and subsequently fitting all power spectra with

*f*as the only free fitting parameter. For the 2.10-μm-diameter particle we used the power spectra for 63 and 806 mW and found

_{c}*ε*

_{x}^{2}= 0.65 nm

^{2}(

*γ*= 2.07 × 10

_{x}^{−8}kg/s) and

*ε*

_{y}^{2}= 1.1 nm

^{2}(

*γ*= 1.98 × 10

_{y}^{−8}kg/s). For the 3.17-μm-diameter particle we used the power spectra for 122 and 921 mW, which gave

*ε*

_{x}^{2}= 0.73 nm

^{2}(

*γ*= 3.18 × 10

_{x}^{−8}kg/s) and

*ε*

_{y}^{2}= 2.3 nm

^{2}(

*γ*= 2.96 × 10

_{y}^{−8}kg/s). These values for

*γ*compare well to the expected values of

*γ*(

*γ*= 2.99 × 10

^{−8}kg/s for a 3.17-μm particle at 20 °C). In Fig. 5(b),

*f*is plotted as a function of the laser power for

_{c}*x*and

*y*, showing a linear relation for both particles. We see no apparent effect of heating.

## 5. Conclusions

*κ*for a trapped particle at different laser powers, we obtained the corner frequency

*f*by fitting our camera power spectra with an adjusted Lorentzian (Eq. (4)) that accounts for motion blur, aliasing, and position detection error

_{c}*ε*

^{2}, as given in [2

2. W. P. Wong and K. Halvorsen, “The effect of integration time on fluctuation measurements: calibrating an optical trap in the presence of motion blur,” Opt. Express **14**(25), 12517–12531 (2006), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-14-25-12517. [CrossRef] [PubMed]

*γ*and

*ε*

^{2}(which can be obtained independently or from a pair of power spectra measured at high and low laser powers),

*f*can be determined for any camera power spectrum up to at least

_{c}*f*= 3 ×

_{c}*f*. In addition, we showed that fitting a pure rather than an adjusted Lorentzian would increasingly underestimate

_{Nyq}*f*for higher laser power.

_{c}*f*with laser power. The ability to calibrate multiple very stiff HOT traps simultaneously will find its application in the use of such traps for quantitative force measurements [10

_{c}10. A. Farré, A. van der Horst, G. A. Blab, B. P. B. Downing, and N. R. Forde, “Stretching single DNA molecules to demonstrate high-force capabilities of holographic optical tweezers,” J. Biophoton. **3**(4), 224-233 (2010). [CrossRef]

## Acknowledgments

## References and links

1. | K. C. Neuman and S. M. Block, “Optical trapping,” Rev. Sci. Instrum. |

2. | W. P. Wong and K. Halvorsen, “The effect of integration time on fluctuation measurements: calibrating an optical trap in the presence of motion blur,” Opt. Express |

3. | A. van der Horst and N. R. Forde, “Calibration of dynamic holographic optical tweezers for force measurements on biomaterials,” Opt. Express |

4. | A. van der Horst, B. P. B. Downing, and N. R. Forde, “Position and intensity modulations in holographic optical traps created by a Liquid Crystal Spatial Light Modulator,” in |

5. | K. Berg-Sørensen and H. Flyvbjerg, “Power spectrum analysis for optical tweezers,” Rev. Sci. Instrum. |

6. | O. Otto, C. Gutsche, F. Kremer, and U. F. Keyser, “Optical tweezers with 2.5 kHz bandwidth video detection for single-colloid electrophoresis,” Rev. Sci. Instrum. |

7. | E. J. G. Peterman, M. A. van Dijk, L. C. Kapitein, and C. F. Schmidt, “Extending the bandwidth of optical-tweezers interferometry,” Rev. Sci. Instrum. |

8. | E. J. G. Peterman, F. Gittes, and C. F. Schmidt, “Laser-induced heating in optical traps,” Biophys. J. |

9. | H. Mao, J. R. Arias-Gonzalez, S. B. Smith, I. Tinoco Jr, and C. Bustamante, “Temperature control methods in a laser tweezers system,” Biophys. J. |

10. | A. Farré, A. van der Horst, G. A. Blab, B. P. B. Downing, and N. R. Forde, “Stretching single DNA molecules to demonstrate high-force capabilities of holographic optical tweezers,” J. Biophoton. |

**OCIS Codes**

(140.7010) Lasers and laser optics : Laser trapping

(170.4520) Medical optics and biotechnology : Optical confinement and manipulation

(350.4855) Other areas of optics : Optical tweezers or optical manipulation

**ToC Category:**

Optical Trapping and Manipulation

**History**

Original Manuscript: February 18, 2010

Revised Manuscript: March 19, 2010

Manuscript Accepted: March 19, 2010

Published: March 29, 2010

**Virtual Issues**

Vol. 5, Iss. 8 *Virtual Journal for Biomedical Optics*

**Citation**

Astrid van der Horst and Nancy R. Forde, "Power spectral analysis for optical trap stiffness calibration from high-speed camera position detection with limited bandwidth," Opt. Express **18**, 7670-7677 (2010)

http://www.opticsinfobase.org/vjbo/abstract.cfm?URI=oe-18-8-7670

Sort: Year | Journal | Reset

### References

- K. C. Neuman and S. M. Block, “Optical trapping,” Rev. Sci. Instrum. 75(9), 2787–2809 (2004). [CrossRef]
- W. P. Wong and K. Halvorsen, “The effect of integration time on fluctuation measurements: calibrating an optical trap in the presence of motion blur,” Opt. Express 14(25), 12517–12531 (2006), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-14-25-12517 . [CrossRef] [PubMed]
- A. van der Horst and N. R. Forde, “Calibration of dynamic holographic optical tweezers for force measurements on biomaterials,” Opt. Express 16(25), 20987–21003 (2008), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-16-25-20987 . [CrossRef] [PubMed]
- A. van der Horst, B. P. B. Downing, and N. R. Forde, “Position and intensity modulations in holographic optical traps created by a Liquid Crystal Spatial Light Modulator,” in Optical Trapping Applications, Vol. 1 of 2009 OSA Technical Digest (CD) (Optical Society of America, 2009), paper OMB3. http://www.opticsinfobase.org/oe/abstract.cfm?URI=OTA-2009-OMB3
- K. Berg-Sørensen and H. Flyvbjerg, “Power spectrum analysis for optical tweezers,” Rev. Sci. Instrum. 75(3), 594–612 (2004). [CrossRef]
- O. Otto, C. Gutsche, F. Kremer, and U. F. Keyser, “Optical tweezers with 2.5 kHz bandwidth video detection for single-colloid electrophoresis,” Rev. Sci. Instrum. 79(2), 023710 (2008). [CrossRef] [PubMed]
- E. J. G. Peterman, M. A. van Dijk, L. C. Kapitein, and C. F. Schmidt, “Extending the bandwidth of optical-tweezers interferometry,” Rev. Sci. Instrum. 74(7), 3246–3249 (2003). [CrossRef]
- E. J. G. Peterman, F. Gittes, and C. F. Schmidt, “Laser-induced heating in optical traps,” Biophys. J. 84(2), 1308–1316 (2003). [CrossRef] [PubMed]
- H. Mao, J. R. Arias-Gonzalez, S. B. Smith, I. Tinoco, and C. Bustamante, “Temperature control methods in a laser tweezers system,” Biophys. J. 89(2), 1308–1316 (2005). [CrossRef] [PubMed]
- A. Farré, A. van der Horst, G. A. Blab, B. P. B. Downing, and N. R. Forde, “Stretching single DNA molecules to demonstrate high-force capabilities of holographic optical tweezers,” J. Biophoton. 3(4), 224-233 (2010). [CrossRef]

## Cited By |
Alert me when this paper is cited |

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

« Previous Article | Next Article »

OSA is a member of CrossRef.