Local stress distribution on the surface of a spherical cell in an optical stretcher

We calculate stress distribution on the surface of a spherical cell trapped by two counter-propagating beams in an optical stretcher in the ray optics regime. We explain the apparition of peaks in the stress distribution, which were not revealed in the earlier published results. We consider the divergence of the incident beams from the fibers, and express the stress distribution as a function of fiber-to-cell distance. In an appendix, we show that the local scattering stress is perpendicular to the spherical refractive surface regardless of incident angle, polarization, the reflectance and transmittance at the surface. Our results may serve as a guideline for the optimization of experimental parameters in optical stretchers. ©2006 Optical Society of America OCIS codes: (140.7010) Trapping; (000.0000) General. References and Links 1. A. Constable, Jinha Kim, J. Mervis, F. Zarinetchi, and M. Prentiss, "Demonstration of a fiber-optical lightforce trap," Opt. Lett. 18, 1867-1869 (1993). 2. J. Guck, R. Ananthakrishnan, H. Mahmood, T. J. Moon, C.C. Cunningham, and J. Käs, "The Optical Stretcher: A Novel Laser Tool to Micromanipulate Cells," Biophys. J. 81, 767-784 (2001). 3. M. Wei, , K. Yang, A. Karmenyan, and A. Chiou, "Three-dimensional optical force field on a Chinese hamster ovary cell in a fiber-optical dual-beam trap," Opt. Express 14, 3056-3064 (2006). 4. Sleep, J., D. Wilson, R. Simmons, and W. Gratzer, "Elasticity of the red cell membrane and its relation to hemolytic disorders: an optical tweezers study," Biophys. J. 77, 3085–3095 (1999). 5. S. Hénon, G. Lenormand, A. Richert, and F. Gallet, "A New Determination of the Shear Modulus of the Human Erythrocyte Membrane Using Optical Tweezers," Biophys. J. 76, 1145-1151, (1999). 6. Y.P. Liu, Chuan Li, A.C.K. Lai, "Experimental study on the deformation of erythrocytes under optically trapping and stretching," Mater. Sci. Eng. A 423, 128–133 (2006). 7. A. L Weisenhornt, M. Khorsandit, S. Kasast, V. Gotzost and H.-J. Butt, "Deformation and height anomaly of soft surfaces studied with an AFM," Nanotechnology 4, 106-113 (1993). 8. R. M. Hochmuth, "Micropipette aspiration of living cells," J. Biomech. 33, 15–22 (2000). 9. A. Ashkin, "Forces of a single-beam gradient laser trap on a dielectric sphere in the ray optics regime," Biophys. J. 61, 569-582 (1992). 10. G. Roosen, "A theoretical and experimental study of the stable equilibrium positions of spheres levitated by two horizontal laser beams," Opt. Commun. 21, 189-194 (1997). 11. P. J. Rodrigo, I. R. Perch-Nielsen, and J. Glückstad, "Three-dimensional forces in GPC-based counterpropagating-beam traps," Opt. Express 14, 5812-5822 (2006). 12. K.C. Neuman, "Characterization of Photodamage to Escherichia coli in Optical Traps," Biophys. J. 77, 2856–2863 (1999).


Introduction
The ability of cells to sense an external tension and to react to the physical stress is intimately related to their functions and properties.In this context, morphological deformations of yeast cells [1], red blood cells (RBCs) [2], and Chinese hamster ovary (CHO) cells [3] under external optical forces have been studied with dual-beam fiber-optical stretchers consisting of two non-focused counter-propagating laser beams from two well-aligned single-mode optical fibers.In an optical stretcher, the cell is confined on the common optical axis of the two beams by the transverse gradient forces, and is stabilized at a point on the optical axis where the two laser beam scattering forces balance each other.As the index of refraction inside the cell is usually greater than that of the surrounding medium, the changes in the photon momentum due to Fresnel reflection and refraction at the interface tend to stretch the cell.The optical stretcher can potentially allow us to measure the elasticity of cells, to sort cells by their membrane viscoelastic properties and to detect differences (or malignancy) in cell's structures for early cancer detection.In these applications, the optical stretcher applies forces that are less localized, smoother and with less disturbance to the cells, in comparison with other approaches such as optical tweezers with beads attached to the cells [4][5][6], atomic force microscope [7] or micropipette [8] so that the elasticity of the whole cell is measured in a more natural environment.
To study the cell's response to external forces in an optical stretcher, the precise stress distribution on the cell surface needs to be determined.The studies on the photonics radiation pressure of a laser trap on a dielectric sphere in the ray optics regime by Ashkin and many authors [9][10][11] are well known.Most works model the trapped particle as a tiny rigid nondeformable sphere.Consequently, only the total force applied to the trapped object is of interest.However, in the dual-beam optical stretcher applications, we are interested in the local stress distribution on the cell's surface and in the concomitant cell's deformation [2].Guck gave the first comprehensive but approximate theory on the photonics stress distribution in the dual-beam fiber-optic stretcher on a spherical deformable red blood cell.In this paper, we present a more precise theory of the scattering stress distribution on the surface of a spherical cell.Our results show local peaks in the stress profile that were not previously outlined.The peaks result from the focusing power of the spherical cell acting as a thick convergent lens.An analytical proof is given in the appendix to show that the local scattering stress is perpendicular to the spherical surface, regardless of the incident angle, polarization of the incident beam, as well as the reflectance and the transmittance at the cell surface.In addition, we consider the divergence of the beams from the fibers, and express the stress distribution as a function of fiber-to-cell distance, which is directly measurable in experiments, for a given fiber numerical aperture.

Theory
Consider a spherical cell trapped in a dual-beam fiber-optical stretcher.Cells do not usually have a spherical shape.However, with red blood cells (RBCs) one can dilute the buffer to have the cell osmotically swollen to a spherical shape.Then, the RBC is modeled as a spherical elastic membrane shell filled with homogeneous and isotropic fluid [2].Experimentally, near-infrared light is often used to minimize the damage to the cell caused by light absorption [12].We use λ=1.064μm in our calculation.For a RBC cell of radius around ρ=3μm, the ray optics regime criterion (2πρ/λ>>1) is satisfied.Any near-infrared wavelength can be used in our theory.Only the calculated results are expected to change accordingly.
When a laser beam encounters a dielectric interface, scattering forces are generated by the change in photon's momentum P=n x E/c [9], where n x is the refractive index of the buffer medium, E is the beam energy and c the speed of light.We denote the momentum of the incident, transmitted and reflected rays by i P , t P and r P , and their directional unit vectors by k a , t a and r a respectively.Then, according to the law of momentum conservation, ( ) , the stress σ applied to the cell's refractive surface is expressed as: where P is the laser beam power, A the area covered by the beam, n = n 2 /n 1 , with n 1 and n 2 being the index of the medium surrounding the cell and inside the cell, respectively, T and R are the Fresnel transmittance and reflectance, respectively, and Q is the dimensionless momentum transfer vector defined by Eq. ( 1).The problem is thus reduced to finding Q .For a stable trap [10], the ratio of the full width at half maximum (FWHM) of the laser beam w over the radius of the cell ρ, (w/ρ) > 1.In Ref. [2] the incident rays are assumed parallel to the optical axis (denoted as the x-axis in Fig. 1).We consider a diverging Gaussian beam from the fiber hitting the front surface (left half) of the sphere.For a given distance from the fiber end to the cell center D, there is a unique relationship between the incident point defined by polar angle φ 1 and the incident angle ε.In fact, as shown in Fig. 1   , and that of the transmitted ray is ε + φ 1 -2β, as shown in Fig. 1.As the reflectance is on the order of 10 -3 at normal incidence for the refractive index n 1 =1.33 for the buffer, and n 2 ≈1.38 for the cells, the third and subsequent reflections inside the cell would have relatively weak power and result in weak stress, which can be neglected.Once all angles are found, one can deduce Q .For the front surface, we have: For the back surface we have: where exp[-2ρ 2 sin 2 (φ 1 )/w 2 ] is a Gaussian beam correction factor.One should note that we made no assumption thus far about the state of polarization, which can affect the reflectance and the transmittance.
It is interesting to prove that the scattering stress is always perpendicular to the spherical refraction surface regardless of the incident angle and no matter whether the rays hit the surface from outside or inside the sphere; both front Q and back Q are perpendicular to the spherical surface.In the appendix we analytically prove, using Eq. ( 1) and Eqs.(2-5), that: Note that this proof is independent of the Fresnel reflectance R and transmittance T and therefore of the incident beam polarization.With the proof of the orthogonality we can write: Equations (7) were given in Ref. [2] without Gaussian beam correction and have been implied without proof that the stress is perpendicular to the spherical surface.For the sake of simplicity, a random polarization is usually considered [2,9] and the average reflectance and transmittance of the parallel and perpendicular polarizations are used.In that case, no parameter in Eqs.(4-7) varies with the meridional angle, so that the three dimensional trapping system is rotationally symmetric around the x-axis and can be analyzed in the x-y plane only.
It is interesting to examine the output polar angle, , where n 1 < n 2 , as a function of the incidence angle ε.For small ε , φ 2 increases with ε monotonically.Then, the increase of φ 2 is slowed down and finally φ 2 decreases with increasing ε for ε > ε ~, as shown in Fig. 2, where ε ~, corresponding to a maximum output polar angle φ 2 ( ε ~), can be computed by the derivative of φ 2 with respect to ε.Consequently, there is an upper limit of the output polar angle φ 2 that depends on the indices n 1 and n 2 , the ratio w/ρ, the fiber NA and the cell radius ρ.As an example, for NA=0.11,n 1 =1.335, n 2 =1.37,D=39.9μm (w/ρ=1.1)and ρ=3μm, φ 2 =71°.There is no incident ray at the front surface, whose refracted ray can hit the back surface at a position of polar angle higher than the upper limit.Below the upper limit there is a range of φ 2 , where the same output position φ 2 can be reached by two different incident angles, as shown in Fig. 2(a).This range is 65°≤φ 2 ≤71° in our example.In the dual-beam optical stretcher the two counter-propagating laser beams generate two stress distributions, which are added up.In the example shown in Fig. 2 added up.One can therefore separate both the front and back surfaces into 4 regions.In the region, one incident ray from outside and one ray from inside the cell hit the same point on the surface.This region is limited by the polar angle between 0°<φ<65° in our example.In the second region, three rays, one incident ray from outside and two rays from inside of the cell hit the same point at the surface.This region is limited by the polar angle between 65°<φ<71° in our example.In the third region only one external ray hits the cell and exerts a stress on the cell surface, and no other ray can hit this point from inside the cell.The third region is limited by the polar angle 71°<φ<87° in our example.The fourth region is limited by the highest position that a ray can hit on the surface for a given fiber NA and fiber-to-cell distance D. Thus, no stress is applied for φ>87° in the example.In the first and second regions, to calculate the stress applied to the front surface at a position φ 1 , we need to find the incidence angle of the ray (ε 2 ) and φ 2 (ε 2 ) coming in the -x direction, which will make a polar angle of φ 1 =(2β-φ 2 ) to the -x direction when hitting the position φ 1 .In other words, we need to solve equation to find φ 2 for a given φ 1.Then we compute Eqs.(2-3) with φ 1 , Eqs. (4-5) with φ 2 and we add the results using Eq. ( 8): A sketch of Q back and Q front is depicted in Fig. 3, where the front surface is denoted as the first surface and the back surface as the second for the incident beam in +x direction, and vice versa for the incident beam propagating in -x direction.Indeed, in the region 0°<φ <65° the stress is the addition of contributions of the two rays.In the region 65°<φ <71°, the stress is the addition of contributions of the three rays and is thus much more intense.In the region φ >71° only the rays hitting the surface from outside of the cell contribute.Finally, in the region φ >87° the stress profile at the first surface is cut to zero.The total stress is the sum of those applied at the first and second surfaces, shown by the two curves in Fig. 3.The stress distribution, which is symmetric in each quadrant of Fig. 1, is shown in Fig. 4 for different distances D from the fibers end to the cell center.We can see the peaks in the stress distribution located at about 60°, 120°, 240° and 300° positions, which were not revealed in the earlier results [2].We anticipate that this difference may influence the predicted cell deformation.At small incident angles, close to 0° and 180°, for a given cell radius ρ the stress strength is proportional to intensity of the input beam, which decreases as 1/D 2 .However the Gaussian beam correction factor exp[-2ρ 2 sin 2 (φ 1 )/w 2 ] increases as the beam size w or the distance D increases.This is why at greater incident angles, close to 90° and 270°, the stress profiles are less sensitive to the distance.We see in Fig. 4 that the surface region near 90°, where no stress is applied, gets smaller as the distance D increases, and that the width of the peak stays almost the same at all distances.

Conclusion
Neglecting the subsequent reflections inside the spherical cell, we have expressed the optical stress distribution as a function of the fiber-to-cell distance, which is directly measurable in experiments for a given fiber NA.We have shown that the focusing power of the spherical cell concentrates the refracted rays to a smaller area on the second interface, resulting in peaks on the stress distribution around certain angular positions.In addition, we have demonstrated that the optical stress is perpendicular to the spherical surface independently on the incident angle, polarization of the incident beam, reflectance and transmittance at the cell surface.

.
The polar angle φ 1 is the incident angle when the beam is parallel to the x-axis.

Fig. 1 .
Fig. 1.Incident, reflected and transmitted rays on a spherical object.The refraction angle β is determined by the Snell's law β ε sin sin 2 1 n n = .After the first refraction, the transmitted ray hits the back surface (right-half) of the sphere from the inside of the cell at a point defined by polar angle φ 2 =2β -φ 1 .Τhe angle of the reflected ray to the xaxis is ) 3 ( 1 φ β π − −

Fig. 2 .
Fig. 2. (a) Position of the output ray on the back surface as a function of incidence angle.Intersection of the horizontal lines and curves are the solutions for two incident rays at a same output position.In the example shown in (b), at incident angles ε =78.2° and 84.3° the rays pass through the cell and hit the back surface at the same point with a polar angle of φ 2 =70°.
(b) on the back surface and at the polar angle position φ 2 =70° the cell is also hit by a third incident ray from the counter-propagating laser beam.At that position, the contributions of the three rays should be #74957 -$15.00USD Received 11 September 2006; revised 20 November 2006; accepted 24 November 2006 (C) 2006 OSA