User:Pranav Rathi/Notebook/OT/2010/12/16/IR Optical Tweezers
Theory of the Optical Trap (physics of the trap
Physics of the optical trap can be explained very well in terms of the size of the particle (micron sphere/bead) relative to the wavelength of the trapping light. Particle size (diameter in comparison to the wavelength helps us deciding which approach (method) is better: Ray optics Or Electric-dipole Approximation.
When the wavelength is larger than the diameter (λ>d: Rayleigh regime, WLDR), electric-dipole method is used, when it is smaller (λ<d: Mie regime, WSDM) ray optics method is used. But there is a discrepancy here: neither method is appropriate, when the particle is in intermediate regime (d≈λ), unfortunately this is our case. To deal it, we will have to extend Mie theory to the case of higher convergent beam, experimental this case is no different. Forces for this case can be fairly approximated thorough ray method.
Ray Optics Method
In the ray optics (geometrical optics) regime, the Gaussian beam can be decomposed into individual rays, each with its own intensity, direction and polarization, which propagates in a straight line of uniform refractive medium, and follow Snell's law at the interfaces. Each ray has the characteristics of a plane wave of zero wavelengths and can change direction (reflection and refraction) at interfaces, and change polarization at dielectric interfaces following Fresnel equations .A graphic view of simple ray model is shown in the figure. The trap consist of an incident parallel beam of Gaussian TEMoo mode and arbitrary polarization enters a high NA numerical aperture objective such the one we are using . The rays are focused one by one to a dimensional less focal point f (original focal point).
- The math for the model is developed on these bases.. All the pictures I am using are from the link.
A ray is incident on the sphere at an angle θ to the surface normal. The linear momentum of light of wavelength λo can be expressed as:
Where p is momentum, E is energy and c is the speed of the light in vacuum. We know the relations c = λof and E = hf so
- where nm is the refractive index of the medium. This is the momentum of a single photon.
From Newton's second law:
Force of N photons:
In Gaussian beam the highest intensity is at the beam axis, where the largest population of photons exist thus highest momentum and force.
In terms of photon flux:
Now force in a scattering medium of refractive index n:
Where P is the incident power and Q is the dimensionless factor. Since force is a function of incident angle and power. Power depends on the Fresnel law of refraction and reflection; this all is taken care by the dimensionless factor. is the incident momentum per second. Where P is the incident power and Q is the dimensionless factor.Since most scattering spheres are in the aqueous medium (water)of some refractive index nm, so the effective refractive index of the particle is . Where np, is sphere refractive index in vacuum. In our case we use Polystyrene beads with refractive index .
The force on the particle can be divided into two: Gradient Force and Scattering Force. Total force is a combination of these two forces. Computation of the total force on the sphere consists of summing the contributions of each beam ray entering the aperture at the radius r with respect to the beam axis and angle β with respect to the y-axis (fig:). The effect of neglecting the finite size of the actual beam focus is negligible for spheres much larger than the wavelength. The point focus of the convergent beam gives the right direction and momentum of the each ray with polarization. The rays then reflect and refract at the surface (interface) of the sphere giving rise to the optical forces. There is a discrepancy here due to the use of highly convergent Gaussian beam. We are ignoring the curvature of the phase-front, that its planner, which is not correct. Gaussian beam has a planner phase-front at the focus which changes to highly curved as beam moves to the far-field (distance larger than the Rayleigh range). Expansion angle (diffraction angle) for highly convergent beam can be as large as 30o; this does not fit with geometrical optics. Another important point is; Gaussian beam propagation formula is strictly correct only for a transverse polarized beam in the limit of small far-field diffraction angle. This formula therefore provides a poor description of a highly convergent beam used in trap.
The proper wave description of a highly convergent beam is much more complex than the Gaussian beam formula. It involves strong axial electric field components at the focus and requires use of the vector wave equations as opposed to the scalar wave equation used for Gaussian beams. So this model does not fit around and at the region of the focal spot, but it is fairly close in the far-field.In our case diffraction angle is , which is larger than the convergence angle φ (65o) inside the water. Inside the bead φ is much larger to be 80o. So for this range this method is applicable.
In Fig: The force due to single ray of power P hitting a dielectric sphere at the in the incident angle θ with incident momentum of . The total(net) force on the sphere is the sum contributions due to reflected ray of power PR and many emergent refracted rays (due to multiple internal reflections) of power PT. The total force is divided into two: Scattered force (FZ,Fs) in the direction of the original incident ray and the gradient force (FY,Fg) orthogonal to the direction of refracted ray. The case shown in the figure is general; particularly the forces are measured in the form of their components in the x, y and z directions acting from the bead center.
Where θ and r are incidence and refraction angles related by Snell's law. R and T are Fresnel reflection and transmission coefficients. Since R and T are polarization dependent the force will be differet for different TE/s and TM/p polarizations. These formulas are sum over all the scattered rays and therefore exact.
Here the polarization orientation is based on the vertical (yz fig: sample plane facing us/the normal to this plane is also a normal into our eyes). In our case the laser is TM/p polarized in the sample plane. So the eqs. 1,2,3 and 4 can be used to computer simulate the trap for our case.
The absolute magnitude of the total force is the vector sub of the components:
The net trapping force is a contribution of the two forces and it is a function of θ and n. For optimum trapping θ is required (see eq. 1 & 2). Optimum θcan be achieved by using a high NA objective with over-field back aperture. In that case incidence angle θ will be equal (close) to convergence angleφ.
- Force along Z-direction
The scattering force along the z-direction can be shown as:
use eq.1 to calculate Fs.
The gradient force along the z-axis can be shown for two cases:
- When the focus is below the sphere center (fig:):
The the zth component of the gradient force is negative and given by:
use eq.2 to calculate Fg. And the total force in the -zdirection:
- When the focus is above the sphere center:
The the zth component of the gradient force is posative and given by:
And the total force in the +zdirection:
Since scattering force is always along the z-direction the equilibrium point is always slightly below the sphere center ( see fig: S≠0). Which mean that the trap center is slightly below the bead center where the backward gradient force just balances the weak forward scattering force. Away from this point the gradient force dominates the scattering force and always has finite net trapping force.
Total force can also be computed for different forcus positions (S) along the z-direction, relative to the sphere center. In such cases incience angle can be computed through:
Where Rb is the bead radius. In ray model forces are independent of the radius Rb=1:
Where ± S is the distance between the focus f (trap canter) and the sphere center O. If focus is below the sphere center S is negative and vice-versa.
- Force along Y-direction:
Fig ; in this case focus of the trapping beam is located along the y-direction. β in the angle between the y-axis and the incidence plane, γ is the convergence angle to the y-axis.α is the angle ray making to the horizontal plane. Knowing α and β we can calculate γ, and using γwe can calculate θ:
WhereS' is the distance between f and O along y-direction.
The net force here particularly depends on the polarization. In eqs. 1& 2 , not only R and T depends on the polarization but also the power P. For the case of an incident beam polarized perpendicular (TE) to the y-axis, for example, one first resolves the polarized electric field E into components E cos β and E sin β perpendicular and parallel to the vertical plane containing the ray. Each of these components can be further resolved into the so-called TM/p and TE/s components parallel and perpendicular to the plane of incidence in terms of these angle μ between the vertical plane and the plane of incidence. By geometry, . The fraction of power lies in the p and s components. This is a general case. In particular case if ray incidence such μ (ray is in the vertical yz-plane) is zero than the polarization can alone be determined on the bases of the vertical plane.
- For TM/p polarization: parallel to the y-axis:
The power in p polarization is:
The power in s polarization is:
- For TE/s polarization: parallel to the y-axis:
The power in s polarization is:
The power in p polarization is:
Knowing θ , PTMand PTE one can compute the gradient and scattering force components for p and s separately using eqs. 1 and 2 and add the results.
- For arbitrary position:
Forces for the arbitrary position of focus is shown below for the first quadrant of the sphere (quadrent of the topright of the fig:). Similar expression are obtained for other quadrants.