| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |

* Division of Physics and Astronomy, Vrije Universiteit, Amsterdam, The Netherlands, and
Department of Physics, Washington State University, Pullman, Washington 99164-2814, USA
Correspondence: Address reprint requests to Erwin J. G. Peterman, Division of Physics and Astronomy, Vrije Universiteit, De Boelelaan 1081 HV, Amsterdam, The Netherlands. Tel.: +31-20-4447987; Fax: +31-20-4447899; E-mail: erwinp{at}nat.vu.nl.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
5 K/W heating in the focus (Wurlitzer et al., 2001Here we report a more direct method to observe temperature changes in a focused laser beam which is based on the analysis of the thermal motion of a trapped bead. Trapping a polystyrene or silica bead is commonly used in optical tweezers experiments. In addition, we will present a model that takes into account more accurately the entire spatial profile of the focused beam in a low numerical aperture (NA) approximation. By comparing data taken in water and glycerol, we show that light absorption by and dissipation in the solvent is the primary determinant of the temperature change, rather than heat absorbed by the trapped particle. The (cooling) effect of the sample cell wall is part of the model and is also demonstrated in the experiments. Presenting data for two cases of solvents with rather different heat conductivities demonstrates the applicability of the model. The model can be used for all other cases (with not too high absorption) as long as heat absorption and conductivity are known, and the results are largely independent of trapped particle properties, provided again that the absorption is not too high. Goals of the paper are to present a correct physical picture of heating in optical tweezers and to provide a practical model for experimenters to use to calculate potential heating effects in their particular situations.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
To visualize the beads, the sample was illuminated with the 546-nm line of an Hg arc lamp (100 W, Zeiss). The transmitted light was collected with an objective (Zeiss, NeoFluar 100x, oil immersion, NA 1.3) and imaged onto a camera (VT1000, Dage-MTI).
The laser-beam path consisted of two 3x beam expanders, a combination of a
/2 plate and a Glan-laser polarizer (to adjust laser intensity), a 1:1 telescope (to position the trap focus in x-, y-, and z-directions), and beam steering mirrors. The laser beam was focused in the sample by the objective. The near-infrared and visible beam paths were separated using dichroic mirrors. Transmitted laser light was collected by the condenser (Zeiss, oil immersion, NA 1.4). Displacement of a bead trapped in the laser focus was measured in two dimensions, normal to the optical axis of the microscope (x, y) by imaging the back-focal plane of the condenser onto a quadrant photodiode (SPOT-9DMI, UDT Sensors) (Gittes and Schmidt, 1998a
). The signals of this photodiode were amplified and anti-alias filtered using custom-built electronics, digitized with a PC board (AT-MIO16X, National Instruments or AD16/ChicoPlus, Innovative Integration), and further processed with custom software (LabView, National Instruments).
To exactly determine the laser intensity in the objective focus, the transmission of the objective was measured by replacing the condenser with an identical objective and measuring the transmitted light. We measured the transmission of this particular objective at 1064 nm as 62 ± 2%, in agreement with values published by others for the same brand and type of objective (Liu et al., 1995
; Neuman et al., 1999
; Svoboda and Block, 1994
).
Optical absorption of water and glycerol
The optical absorption of glycerol and water at 1064 nm was determined by measuring the transmitted intensity of a laser beam passing through a cuvette filled with a variable path length of the respective liquid. The cuvette was made from an acrylic glass cylinder, oriented vertically, open on the top and closed with a cover glass at the bottom. Liquid was pumped in with a syringe connected to the side of the cylinder, close to the bottom. Another coverslip was placed inside the cylinder floating on top of the liquid layer, in order to obtain a flat meniscus. The water and glycerol where filtered through a 0.2-µm syringe filter and care was taken to avoid air bubbles. The height of the liquid column was measured with a precision ruler. Laser light (from a diode-pumped Nd:YVO4; Topaz 106C, Spectra Physics, Mountain View, CA) was sent through the cuvette and the transmitted light intensity was measured with a laser power meter (Newport 1815-C with detector head 818T-10) as a function of the variable liquid path length in the cuvette as shown in Fig. 2. The extinction coefficients could be extracted from these data by fitting a single exponential function. The extinction coefficients were found to be 14.2 m-1 for water (similar to the
12 m-1 obtained from the graph in Hale and Querry (1973)
and 21.4 m-1 for glycerol.
|
, Clarkson University, NY) were used. Unless stated differently, single beads were trapped
10 µm above the coverslip-liquid interface to avoid thermal and hydrodynamic surface effects. The effect of heating in the focus of the trapping beam was measured in two ways. First, power spectra of the Brownian motion of a bead in the trap were measured as function of the laser power, and second, displacements of trapped beads were measured upon exertion of viscous drag (by moving the sample chamber using the piezo-actuated stage), also as a function of laser power.
| RESULTS |
|---|
|
|
|---|
![]() | (1) |
= 3
x
x d is the Stokes' drag of the bead with diameter d in a solvent with viscosity
, kB is the Boltzmann constant, T is the absolute temperature,
is the trap stiffness, and f0 is the corner frequency of the spectrum. Some typical power spectra of a 502-nm-diameter polystyrene bead in glycerol are shown in Fig. 3. If the laser power only affected the trap stiffness
(see Eq. 1), the corner frequency f0 would increase with increasing power. For high frequencies, the f02 term in the denominator of Eq. 1 becomes negligible, and the spectral density would become independent of trap stiffness and thus independent of laser power. If, on the other hand, the temperature, T, in the vicinity of the bead were to change with power, the spectral density would change also at high frequencies. Such an effect is evident in the data (Fig. 3) from the fact that the slanting parts of the spectra (f > f0) do not overlap.
|
From Eq. 1, the laserpower dependence of the two fitted parameters (f0 and S0 x f02) can be expressed as
![]() | (2) |
![]() | (3) |
![]() |
The power dependence of P/f0 and 1/(S0 x f02) for a polystyrene bead with diameter 502 nm in glycerol and a silica bead with diameter 444 nm in water are shown in Fig. 4. Also shown are fits of the data by Eqs. 2 and 3. The fit results for 502-nm (the experiments shown in Figs. 24
), 840-nm, and 2200-nm polystyrene beads in glycerol and 444-nm (the experiments shown in Fig. 2) and 500-nm silica beads in water are listed in Table 1. These experiments indicate that in the case of glycerol the temperature increase is 41.1 ± 0.7 K/W laser light for the 502-nm beads, 43.8 ± 2.2 K/W for the 840-nm beads, and 34.2 ± 0.1 K/W for the 2200-nm beads. There is only a small effect coming from the size of the trapped bead. In water the temperature increase is substantially smaller than in glycerol, namely 7.7 ± 1.2 K/W for 500-nm beads and 8.1 ± 2.1 K/W for 444-nm beads. This is due to the lower absorption of 1064-nm light by water and the higher thermal conductivity of water (see below).
|
|
|
A value for the temperature coefficient B in glycerol was obtained in an independent way (for a distance of the bead from the coverslip of 10 µm). The displacement of a trapped bead with respect to the center of the trap as a function of the laser power was measured while moving the sample stage (thus the surrounding glycerol) back and forth with constant speed. This displacement (
x) is proportional to the force exerted on the bead and is equal to
![]() | (4) |
![]() | (5) |
Displacements (multiplied with the laser power) as a function of laser power are shown in Fig. 5 C. Also shown is a fit of the data by Eq. 5. From this fit a temperature factor of 36.6 ± 2.0 K/W was obtained (see Table 1). For the 444-nm beads a temperature factor of 39.2 ± 2.2 K/W was measured. The values obtained for the laser-induced temperature increase obtained from the fluid-drag experiments are in agreement with those obtained from the power spectra (Table 1).
| THEORETICAL MODEL |
|---|
|
|
|---|
|
![]() | (6) |
is the extinction coefficient of the solvent, defined by I(x) = I0 x e-
x.
The heat generated by the absorption will dissipate. The dissipating heat flow
obeys the local differential equation
![]() | (7) |
1.4 W/(m x K); Weast, 1973
is the temperature deviation from ambient temperature due to heating at position
. In steady state the heat dissipated is equal to the heat generated:
J =
Q/
t. Using Eqs. (6) and (7), this becomes:
![]() | (8) |
, a model for the intensity profile of the laser is needed. Here a more accurate approximation of the intensity profile is considered than in a model published before (Liu et al., 1995
, and
, with
= 0 along the incoming optical axis). To describe the beam intensity I, suppose that
![]() | (9) |
<
/2, outgoing for
>
/2). At high NA, f(
) could be quite a broad function of angle. At large distances, this incoming or outgoing intensity falls off as an inverse-square law,
![]() | (10) |
) is normalized over incoming solid angle,
![]() | (11) |
=sin
d
d
. (The attenuation in Eq. 6 is assumed to be small, so that the outgoing power has been set equal to the incoming power.) To clarify Eq. 9, we derive the function f(
) and the constant a for the case of a Gaussian beam, which is a low NA approximation for a focused laser. Manipulating well-known expressions (Siegman, 1986
![]() | (12) |
0 is the equivalent angular radius of a uniformly illuminated aperture with the same total power. In terms of
0 one finds that the intensity of the Gaussian beam through the focus (along the
= 0 axis) has the form
![]() | (13) |
/2
02. The NA in this approximation is just
0
sin
0, so a plausible extrapolation to NA
1 gives a =
/2
in Eq. 9.
We are interested in solving Eq. 8 for the temperature in the focus, T(0). The Green's function for Eq. 8 satisfies
![]() | (14) |
Here,
is the Dirac delta function. The Green's function with its argument
' set to zero and obeying the boundary condition G = 0 at some large radial distance R, is
![]() | (15) |
The G = 0 boundary condition might represent, for example, distant glass surfaces of higher heat conductivity (and lower attenuation). In terms of G(r), the solution
T(r) of Eq. 8 at the origin is
![]() | (16) |
V = r2 dr d
,
![]() | (17) |
/2
to get
![]() | (18) |
= 1064 nm), close to the experimentally determined values. It is important to note that the experimentally found difference between glycerol and water is reproduced by the model, indicating that indeed the determining parameters for the heat generation are absorption of light, as governed by the extinction coefficients,
, and dissipation of the heat, as governed by the thermal conductivity, C. It should be noted that the dependence of the temperature difference on R is not very strong: if R was 100 µm, the temperature increase would be 21 K/W for water and 66 K/W for glycerol. Our experiments show that, indeed, the temperature increase is larger if the distance to the cover glass is larger. In general, in diffusion problems the shortest length scale dominates the effect, which is reflected in the logarithmic distance dependence. In reality the cutoff (i.e., the glass coverslip) is close only on one side. To obtain readily solvable equations, a symmetric cutoff R was assumed (as if the bead were in the center of a chamber, with varying thickness 2R). In Fig. 4 E the data is compared to the distance dependence of Eq. 8. Within the error margins of the experiment this simplified model describes the data well without any fitting parameters.
In the preceding discussion we ignored the influence of the trapped particle. To obtain a more accurate description of the heating effect, a correction term can be added to Eq. 18 that incorporates the difference in absorption and thermal conductivity of the bead compared to the solvent. We represent these effects as an effective change in the coefficient b over the volume of the bead:
![]() | (19) |
Our experiments on fixed, trapped beads gave the same results as the drag experiments (in which the bead was held in place and the solvent was moved) suggesting that the heat equilibration is faster than the movement in these experiments. To estimate time scales for the equilibration time for any given temperature distribution, we need not only C (with units J/sxmxK) but also the heat capacity per volume cV (for water 4.2x106 J/m3xK, for glycerol 3.0x106 J/m3xK). On dimensional grounds, the equilibration time to a given distance r must be on the order of r2xcV/C. Thus, for both water and glycerol, temperature equilibration out to 10 µm takes roughly 1 ms. The time-dependent heating model by Schönle and Hell (1998)
does not apply to our situation, because they assume an infinite sample size, which causes the infinite equilibration times they report. Temperature kinetics should thus not be relevant even if the solvent is moved fairly quickly. In our drag experiments the sample is moved with a speed of 910 nm/s, slow enough to allow temperature equilibration.
| DISCUSSION |
|---|
|
|
|---|
This value for the laser-induced heating is comparable to the values measured before (5 K/W, Wurlitzer et al., 2001
;
1015 K/W; Liu et al., 1995
, 1996
). We think that our measurements are likely to be more accurate since well-defined solid particles (polystyrene, silica) were trapped, in contrast to some of the previous experiments where vesicles and cells were used. We have shown both experimentally and in our model that the effect of the trapped (polystyrene, silica) particles is only minor; the main contribution to the heating is light absorption by the surrounding solvent. The distance from the glass-solvent interface, however, has a substantial effect on the heating. We have shown that the heating increases substantially when the bead is trapped further away from the interface.
It remains to be discussed for which kind of experiments these results are relevant. Most optical trapping experiments use water (or watery solutions) as solvent and laser powers on the order of 100 mW (at 1064 nm). This leads to a temperature increase of only
0.8 K in the focus, which, depending on the solvent, may change the viscosity more dramatically. In many cases the Lorentzian fit to a power spectrum of a trapped bead is used for the calibration of the trap and detector response (Gittes and Schmidt, 1998b
). If the heating effect is not taken into consideration using this calibration method, the trap stiffnesswhich is proportional to the estimated viscosity (
) times the measured corner frequency (f0)is overestimated 2% (10%) when a laser power of 100 mW (500 mW) is used. Here we assume a temperature increase of 8 K/W and a base temperature of 294.55 K. The detector response (in m/V) is proportional to the temperature (T) divided by the viscosity (
) and the zero-frequency intercept of the power spectrum (S0 f02) and is in the same circumstances underestimated 2% and (11%). Consequently, heating effects due to laser-light absorption by the solvent in optical trapping experiments even in watery solution have a small but measurable effect, and should be taken into consideration, especially when laser powers higher than
100 mW are used.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
, Clarkson University, New York, for the supply of silica particles. This research was supported by the Foundation for Fundamental Research on Matter (FOM). E.P. is supported by a Postdocs Universitaire Loopbaan Stimulerings Programma (PULS)-fellowship from the Research Council for Earth and Life Sciences (ALW) with financial aid from the Dutch Organization for Scientific Research (NWO).
Submitted on November 7, 2001; accepted for publication October 16, 2002.
| REFERENCES |
|---|
|
|
|---|
Ashkin, A. 1997. Optical trapping and manipulation of neutral particles using lasers. Proc. Natl. Acad. Sci. USA. 94:48534860.
Ashkin, A., J. M. Dziedzic, and T. Yamane. 1987. Optical trapping and manipulation of single cells using infrared laser beams. Nature. 330:769771.[Medline]
Block, S. M. 1990. Optical tweezers: a new tool for biophysics. In Noninvasive Techniques in Cell Biology. S. Grinstein, K. Foskett, editors. Wiley-Liss, New York, New York. 375401.
Daubert, T. E., and R. P. Danner. 1989. Physical and Thermodynamic Properties of Pure Chemicals: Data Compilation. Hemisphere Pub. Corp., New York, New York.
Finer, J. T., R. M. Simmons, and J. A. Spudich. 1994. Single Myosin molecule mechanicspiconewton forces and nanometer steps. Nature. 368:113119.[Medline]
Gittes, F., and C. F. Schmidt. 1996. Microscopic approaches to dynamics and structure of biological motors. Curr. Opin. Solid State Mat. Sci. 1:412424.
Gittes, F., and C. F. Schmidt. 1998a. Signals and noise in micromechanical measurements. Methods Cell Biol. 55:129156.[Medline]
Gittes, F., and C. F. Schmidt. 1998b. Thermal noise limitations on micromechanical experiments. Eur. Biophys. J. Biophys. Lett. 27:7581.
Hale, G. M., and M. R. Querry. 1973. Optical constants of water in the 200-nm to 200-µm wavelength region. Appl. Optics. 12:555563.
Kellermayer, M. S. Z., S. B. Smith, H. L. Granzier, and C. Bustamante. 1997. Folding-unfolding transitions in single titin molecules characterized with laser tweezers. Science. 276:11121116.
Liu, Y., D. K. Cheng, G. J. Sonek, M. W. Berns, C. F. Chapman, and B. J. Tromberg. 1995. Evidence for localized cell heating induced by infrared optical tweezers. Biophys. J. 68:21372144.
Liu, Y., G. J. Sonek, M. W. Berns, and B. J. Tromberg. 1996. Physiological monitoring of optically trapped cells: assessing the effects of confinement by 1064-nm laser tweezers using microfluorometry. Biophys. J. 71:21582167.
Lytle, J. D., G. W. Wilkerson, and J. G. Jaramillo. 1979. Wideband optical transmission properties of seven thermoplastics. Appl. Optics. 18:18421846.
MacKintosh, F. C., and C. F. Schmidt. 1999. Microrheology. Curr. Opin. Colloid Interface Sci. 4:300307.
Mehta, A. D., M. Rief, J. A. Spudich, D. A. Smith, and R. M. Simmons. 1999. Single-molecule biomechanics with optical methods. Science. 283:16891695.
Molloy, J. E., and M. J. Padgett. 2002. Lights, action: optical tweezers. Contemp. Phys. 43:241258.
Neuman, K. C., E. H. Chadd, G. F. Liou, K. Bergman, and S. M. Block. 1999. Characterization of photodamage to Escherichia coli in optical traps. Biophys. J. 77:28562863.
Schönle, A., and S. W. Hell. 1998. Heating by absorption in the focus of an objective lens. Optics Lett. 23:325327.[Medline]
Siegman, A. E. 1986. Lasers. University Science Books, Mill Valley, California.
Smith, S. B., Y. J. Cui, and C. Bustamante. 1996. Overstretching B-DNA: the elastic response of individual double-stranded and single-stranded DNA molecules. Science. 271:795799.[Abstract]
Svoboda, K., and S. M. Block. 1994. Biological applications of optical forces. Annu. Rev. Biophys. Biomolec. Struct. 23:247285.[Medline]
Svoboda, K., C. F. Schmidt, B. J. Schnapp, and S. M. Block. 1993. Direct observation of kinesin stepping by optical trapping interferometry. Nature. 365:721727.[Medline]
Tskhovrebova, L., J. Trinick, J. A. Sleep, and R. M. Simmons. 1997. Elasticity and unfolding of single molecules of the giant muscle protein titin. Nature. 387:308312.[Medline]
Weast, R. C., editor. 1973. CRC Handbook of Chemistry and Physics, 53rd ed. CRC Press, Boca Raton, Florida.
Wurlitzer, S., C. Lautz, M. Liley, C. Duschl, and T. M. Fischer. 2001. Micromanipulation of Langmuir-monolayers with optical tweezers. J. Phys. Chem. B. 105:182187.
This article has been cited by other articles:
![]() |
Y. Zhang, G. Sun, S. Lu, N. Li, and M. Long Low Spring Constant Regulates P-Selectin-PSGL-1 Bond Rupture Biophys. J., December 1, 2008; 95(11): 5439 - 5448. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. B. Rasmussen, L. B. Oddershede, and H. Siegumfeldt Optical Tweezers Cause Physiological Damage to Escherichia coli and Listeria Bacteria Appl. Envir. Microbiol., April 15, 2008; 74(8): 2441 - 2446. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. V. Desai, T. G. Bishop, L. Vicci, E. T. O'Brien Sr., R. M. Taylor 2nd, and R. Superfine Agnostic Particle Tracking for Three-Dimensional Motion of Cellular Granules and Membrane-Tethered Bead Dynamics Biophys. J., March 15, 2008; 94(6): 2374 - 2384. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. E. Wallin, A. Salmi, and R. Tuma Step Length Measurement--Theory and Simulation for Tethered Bead Constant-Force Single Molecule Assay Biophys. J., August 1, 2007; 93(3): 795 - 805. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Nassoy How Accurate Are Ultrasensitive Biophysical Force Probes? Biophys. J., July 15, 2007; 93(2): 361 - 362. [Full Text] [PDF] |
||||
![]() |
M. Andersson, E. Fallman, B. E. Uhlin, and O. Axner Dynamic Force Spectroscopy of E. coli P Pili Biophys. J., October 1, 2006; 91(7): 2717 - 2725. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Zondervan, F. Kulzer, H. van der Meer, J. A. J. M. Disselhorst, and M. Orrit Laser-Driven Microsecond Temperature Cycles Analyzed by Fluorescence Polarization Microscopy Biophys. J., April 15, 2006; 90(8): 2958 - 2969. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Titushkin and M. Cho Distinct Membrane Mechanical Properties of Human Mesenchymal Stem Cells Determined Using Laser Optical Tweezers Biophys. J., April 1, 2006; 90(7): 2582 - 2591. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. A. Abbondanzieri, J. W. Shaevitz, and S. M. Block Picocalorimetry of Transcription by RNA Polymerase Biophys. J., December 1, 2005; 89(6): L61 - L63. [Abstract] [Full Text] [PDF] |
||||
![]() |
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., August 1, 2005; 89(2): 1308 - 1316. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. H. Werner, H. Cai, R. A. Keller, and P. M. Goodwin Exonuclease I Hydrolyzes DNA with a Distribution of Rates Biophys. J., February 1, 2005; 88(2): 1403 - 1412. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. F. Tolic-Norrelykke, A. M. Engh, R. Landick, and J. Gelles Diversity in the Rates of Transcript Elongation by Single RNA Polymerase Molecules J. Biol. Chem., January 30, 2004; 279(5): 3292 - 3299. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |