Switching and measuring a force of 25 femtoNewtons with an optical trap

An experiment is described, where an ultra-weak force of 25 femtoNewtons (fN = 10 N) displaces a 0.53μm latex bead captured in a soft optical trap. The displacement is measured by exploiting the phase anomaly of the trapping beam, resulting in a small phase shift of the scattered light. The 25fN force stems from the radiation pressure of a nearly collimated second laser beam, which is switched by an accousto-optical modulator. To the best of my knowledge, this is the smallest switchable force which has ever been directly measured. © 2005 Optical Society of America OCIS Codes: (140.7010) Trapping; (300.2530) Fluorescence, laser-induced; (290.2200) Extinction; (290.4020) Mie theory; (260.3160) Interference


Introduction
Although photons are massless, they can exert a pressure on matter, by transferring their momentum h/λ, which is the inverse incident wavelength quantized by Planck's constant. Although one photon carries a tiny momentum of 6.626 × 10 −28 Ns at a wavelength of λ = 1 µm, this momentum sums up to a gigantic radiation pressure in the center of hot stars. This pressure works against the gravitational forces, which would induce a collapse of the star. In principle the photon momentum can be exploited for an interplanetary transport of spaceships with solar sails. Large, photon reflecting surfaces (sails) with low mass, are slowly, but almost frictionless accelerated by the photon flux emitted by the sun. This situation changes completely in the earth's atmosphere, where collisions with air molecules dissipate the kinetic energy of the objects accelerated by radiation pressure. The question arises by how far these forces can be scaled down to laboratory situations.
Already 1970 Ashkin [1] dragged micron-sized latex spheres with an Argon laser and measured radiation forces in the picoNewton (pN) range by finding the equilibrium of optical and friction forces. However, the smaller and lighter the optically accelerated particles, the more pronounced becomes the viscous drag force γ⋅v and thermal random forces F th (leading to Brownian motion). For example, a 1µm latex bead can take up a kinetic or thermal energy of 8 k B T (k B T ≈ 4⋅10 -21 Nm at T = 300 K, k B = Boltzmann constant) after a few seconds [2]. By moving over a distance of e.g. Δx = 30nm against an external potential W(x), the effective thermal force can easily become F th = ∂ x 2 W(x) ⋅ Δx = 1 pN. The gravity G = m⋅g of this bead is only 10 fN. A molecular motor transporting this 1µm bead along fibers in a living cell exerts forces of 3-5 pN in this Brownian storm of thermal forces. Forces in the pN regime have been successfully measured with optical tweezers (see review of Neuman and Block [3]), whereas direct force measurements in the low fN regime are very difficult to perform. With a direct measurement of a force F mes I refer to the equilibration with a single calibrated force F cal = F mes .
Where do we find fN forces? For instance, the entropic force that arises when the entropy is reduced by stretching e.g. a DNA polymer with a persistence length of L p = 50 nm is F ent = k B T/L p = 80 fN. A 1µm small bacterium has to overcome a friction force of F fr = γ⋅v = 95 fN at a propagation speed of v = 10 µm/s in water. The gravitational attraction between two 3 mm glass spheres close to each other is about 10 fN. The repulsive electrostatic force of a 0.5 µm latex sphere in de-ionized water is about 20 fN in a distance of 1.5 µm to the surface due to the auto-dissociation of water. In his famous experiments to prove the quantization of electric charge [4], Millikan switched a 100-200 fN small electric force F E = E⋅n⋅q to counterbalance the gravitational force m⋅g on 3-4 µm sized oil drops, loaded with n elementary charges q. And coming back to the radiation pressure, the daily sun light, defined by the solar intensity I sol = 1370 W/m 2 exerts a force F rad = I 0 ⋅n/c ⋅Q pr ⋅(a 2 π) of only 20 fN on a small absorbing dust particle of 2a = 75µm in diameter in the atmosphere. I 0 = I sol is the intensity incident on the particle's cross-section (a 2 π), c/n is the speed of light in a medium with refractive index n and Q pr is the pressure efficiency, which is 1 for an absorbing particle. I have setup an arbitrary collection of ultra-small forces to give the reader a feeling for magnitudes in the low femtoNewton range. Force measurements in the low femtoNewton range are mentioned regularly in the literature [5][6][7][8][9], but attention has to be paid, whether these studies present a direct measurement or an estimation of the apparatus' resolution for force measurements.
So far no report has been found in the literature where such small forces could be controlled and measured directly. In this paper I demonstrate how a force of 25 fN on a small probe is switched and measured directly with help of an ultra-weak optical trap and an interferometric position detection system. The force can be switched on and off, such that no extrapolation of a system inherent force distribution is necessary. Only two parameters need to be determined for the measurement: the stiffness κ z of the optical trap [10], which is obtained by signal autocorrelations 〈 S z (t)⋅S z (t+τ)〉 of the probe's position z(t), and secondly, the shift Δz of the probe's mean position, measured before and after the femtoNewton force was switched on. The measured force is then F z = κ z ⋅Δz. The displacement Δz is obtained via the tiny phase change of the trapping beam scattered at the probe.

Instrumentation and experiment
Which femtoNewton force can be quickly switched on and off and leaves the system in thermal equilibrium in both the on-and off-state? A conservative force, which can be easily controlled, results from the radiation pressure of a collimated laser beam on a non-absorbing sphere. A 488nm Argon laser beam was switched with a accousto-optical modulator (AOM) leading to an optical power of P = 11 mW at the dielectric probe, which was a 0.535 µm latex bead with refractive index n s = 1.591+0.0003i [11].
The ultra-low force measurement was performed with a photonic force microscope (PFM), which is an optical tweezers related system with a 3D interferometric probe tracking system. Technical details of the PFM are described elsewhere [12]. The experimental setup is illustrated by Fig. 1. An unexpanded 488nm laser beam is weakly focused into the back focal plane (BFP) of the trapping lens. This results in a nearly collimated beam with waist s in the focal plane of the trapping lens. A trapped latex bead with radius a = 268 nm (Bangs Labs

Fluorescence collection
Inc.) is moved in axial direction by the radiation pressure of the 488nm laser beam, until the radiation force of the 488nm beam equilibrates with the trapping force of the 1064nm beam. The trapping force F trap is obtained by calibrating both the optical trap and the detection system. The axial force F trap = κ z ⋅ Δz is calibrated by the axial stiffness κ z = γ/τ z and the detection system is calibrated by the axial signal sensitivity g z = σ sz / σ z = σ sz ⋅(κ z /k B T) 1/2 of the particle's position signal S z (z) = g z ⋅z, which is linear to the position z over a reasonable distance [12]. The mean square displacement of the trapped particle σ z 2 :=〈σ z (t) 2 〉 follows the equipartition theorem κ z ⋅σ z 2 = k B T. The resulting half width σ sz of the signal histogram ( Fig.  2(b)) is obtained by a Gaussian fit exp[-(S z /σ sz ) 2 ]. The autocorrelation time τ z is the time the particle needs to relax inside the optical trap in axial direction and is provided by the signal autocorrelations is the viscous drag far from an interface (η viscosity of water). Therefore the relation between the particle's position z(t) and it's position signal S z (z(t)) is Therefore the radiation force from the 488nm laser is provided by the force equilibrium The axial displacement of the particle Δz is proportional to the difference ΔS z in the mean position signals on off on off The position signal S z (t) = 〈 S〉 ± σ sz (t) fluctuates by σ sz (t) = g z ⋅σ z (t) around the mean value 〈 S〉 due to thermal motion (see Fig. 2(a)). The integration times were ΔT on ≈ ΔT off ≈ 4s at a sampling rate of 50 kHz. The temperature was T = 21°C. The signals S z are acquired by integrating the interference intensity over the full area of the QPD [12], the signal sensitivity g z (NA det ) can be increased by decreasing the numerical aperture of the detection lens NA det [13,14], which we reduced to NA det = 0.3. S z is the integrated interference intensity between the incident field E i and the scattered field E s (z) which varies with the position of the probe z.
Close to the focus, the phase difference Δφ(Δz) = φ i -φ s (Δz) ~ Δz is proportional to the particle's displacement Δz due to the Gouy phase shift [15,16]. The Gouy phase shift leads to the reduced momentum carried by a focused beam in axial direction. Via a phase shift of π it leads to a stretching of the mean wavelength in the focal area. Once the detection system is calibrated, i.e. g z = ΔS z / Δz has been determined, an alternative way to gain the radiation force is to subtract the off-state potential V off (z) from the on-state potential V on (z) = V off (z) + V rad (z). The potentials V(z) are related to the probability density p(z) of the histograms via Boltzmann statistics p(z) ~ exp[-V(z)/k B T]. The radiation force then is the derivative of the radiation potential, i.e. F rad = ∂ / ∂z V rad (z). This method is illustrated in Fig. 2(c). To estimate the number of photons transferring their momentum on to the probe, the intensity I 0 of the 488nm laser beam at the probe was determined. As illustrated in Fig. 1, a diluted fluorescein solution (mass ratio 2⋅10 -5 ) was used to measure the beam waist s of the beam. After fluorescence excitation, the beam profile was imaged onto a CCD camera (see Fig. 3). The axial extent of the focus, i.e. the extent of the detection plane (focal plane) is about 1µm, whereas the Rayleigh length of the beam is much longer, ensuring that no artifact is introduced by the integration of the detection lens. By fitting a Gaussian beam profile to the fluorescence intensity I(r) = I 0 ⋅exp(-r 2 /s 2 ), a half width of s = 7µm at I 0 /e was obtained (Fig. 3  left). The trapped probe was well centered relative to the 488nm beam, as indicated in Fig. 3 right. Here the fluorescence increase due to coherent scattering behind the sphere is 70%, integrated over the depth of field of the trapping lens.

Results and discussion
The total power P 488 = 11 mW in the focal plane is 62% of the power measured behind the coverslip without the trapping lens. P is distributed over the whole focal plane according to From Eq. (6) we find I 0 = P 488 / (πs 2 ) = 71.5 MW/m 2 , which corresponds to a power of P geo = I 0 ⋅a 2 π = 16.1 µW distributed over the cross section of the sphere, which is C geo = a 2 π = 0.226 µm 2 . In other words P geo /(hν) ≈ 5.3⋅10 13 photons per second hit the sphere's cross section.
The resulting radiation force can be measured by using eq.(3). Totally 11 measurements were performed at two strengths of the optical trap (λ trap = 1064nm). At a very low trapping laser power P trap = 1.5mW in the focal plane, a radiation force of F rad = 25.6 ± 2.2 fN was obtained (N=5 measurements). This force is the product of a mean trap stiffness of κ z = 0.37 pN/µm and a mean displacement Δz = 69nm. At a second, still low laser power P trap = 4.4mW nearly the same radiation force of F rad = 22.3 ± 5.7 fN was measured (N=6), which is the product of a mean trap stiffness of κ z = 1.12 pN/µm and a mean displacement Δz = 20nm. The consistency and the relative small errors underline the reliability of the measurements. The uncertainty in the force measurement due to thermal fluctuations is F ± ΔF = κ z ⋅(Δz ± σ z ) with Fig. 3. The Gaussian profile of the weakly focused laser beam is visible due to fluorescence emission in the fluorescein solution (left image). The intensity profile reveals a beam waist of 2s = 14µ m at a 1/e decay. Right: A bead trapped in the center of visible laser beam produces an increase in electric field density behind the bead due to scattering. In this region the fluorescence intensity is increased by 70% as indicated by the profile and the photo. Integration time was 20ms. σ z 2 = k B T/κ z ). The uncertainty ΔF increases with (k B T κ z ) 1/2 at a constant external force F rad and can be reduced by using a weak optical trap, i.e. a small trapping stiffness κ z , and, of course, by a low temperature T. Since resolution is a question of definition, I chose that ΔF ≤ F or σ z ≤ Δz, which was approximately fulfilled in this experiment (Δz ≈ σ z ≈ 70nm). By defining a different resolution limit or by achieving an even smaller trap stiffness, forces smaller than 25fN can be measured.
As a further verification, a similar radiation force F rad (< 15% deviation) was determined by the gradient of the potential difference as illustrated in Fig. 2 right. Furthermore, the correctness of the detector calibration g z via the relaxation time τ z (eq.(2)) was tested by moving a fixed bead over a distance Δz through the focus by a piezo-stage. The ratio g z = ΔS z / Δz provided nearly the same value.
According to Mie-theory, a sphere with index n s = n R +i⋅n I = 1.591+0.0003i and radius a = 268nm has the scattering function S(θ=0) = 7.715-10.548i leading to a scattering efficiency Q sca = 4/(k⋅a) 2 ⋅Re{S(0)} = 1.4653, k = 2π n/λ, λ = 488nm. This would mean that 147% of the incident photons are scattered at the sphere (extinction paradox) [17]. However, due to diffraction, i.e. the redistribution of optical energy and momentum, more photons are scattered in forward direction (θ=0) due to interference with the unscattered field. The resulting pressure efficiency without absorption is Q pr -Q abs = 0.180, where ( ) Here n is the refractive index of the medium (water) and Q abs is the absorption efficiency. T 1 (θ) and T 2 (θ) are the angle-dependent Mie-scattering functions parallel and perpendicular to the plane of incidence [18]. Using Eq. (1), we obtain a theoretical radiation pressure of F rad = 12.8fN, if Q abs = 0. This theoretical value is only half as large as the measured value for the radiation force. It is therefore likely that absorption takes place at the sphere. The momentum uptake ~ Q abs due to absorption and incoherent scattering must be about the same as due to coherent scattering ~(Q pr -Q abs ). Assuming in the simplest case a dense monolayer of fluorophores, we expect a number N of absorbing fluorophores on the front and back side of the bead with geometrical cross-section C geo = a 2 π. With a geometrical cross-section r 2 π of the fluorophores, the number of molecules is N = 2a 2 π / r 2 π = 2.26⋅10 6 . The diameter of a fluorescein fluorophore is about 2r = 0.5nm and the molar extinction coefficient of fluorescein is about C ext,fl = 72000 (mol/l) -1 cm -1 (see Sigma-Genosys, www.fisheroligos.com/tec_flc.htm), which is C ext,fl = 72000 m 2 /N A = 0.012 (nm) 2 (N A = Lohschmidt number, C sca,fl = 0 due to incoherent, isotropic scattering). This results in a total absorption cross-section of C abs = N⋅C ext,fl = 0.034 µm 2 or an absorption efficiency Q abs = C abs / C geo = 0.15, which is about the same as Q pr -Q abs = 0.18 due to coherent scattering at the sphere. Again, according to this estimation, a latex sphere with diameter 535nm takes up the same amount of light momentum as a 0.5nm thin monolayer of fluorescein molecules around the sphere. The emitted fluorescence of the dye monolayer is small in comparison to the increased excitation and fluorescence behind the sphere due to coherent scattering of the excitation light.
Some years ago Matsuo et al. [6] tried to exploit the difference between the measured radiation pressure efficiency on a particle and the computed scattering efficiency to estimate the absorption of the particle. However, in this study unique measurements were complicated due to a more complex scattering at the 10-fold larger spheres and the close proximity of two interfaces.
Finally I want to remark that also the absorbance of the 488nm light by fluorophores inside the fluorescein droplet should have an effect on the final momentum transfer, because it mainly reduces I 0 at the sphere. An estimate of this reduction is complicated due to the diffusion of fluorophores through the beam.
Although different ratios of momentum uptake due to resonant and non-resonant scattering at the (coated) sphere may result in a radiation force of 25fN, the force switching and measurement remains unaffected by this ratio.

Conclusions
The smallest switchable force ever was measured by using an ultra-weak optical trap. Different measurements manifest a force of about 25 fN resulting from the radiation pressure of a 11mW laser beam (λ = 488nm) on a 0.53µm latex sphere. About half of this force results from momentum transfer due to coherent scattering, the other half is likely due to absorption of a dense monolayer of fluorescein molecules around the sphere. With this study I try to set a new milestone concerning force measurements in the sub-picoNewton regime and I encourage other researchers to underbid this value with careful measurements.