1. Introduction
The operation of the OT, part of the idea that light exerts pressure on material objects, phenomenon called radiation pressure. The first speculations on the existence of the phenomenon go back to the XVII century with the astronomical observations of Johannes Kepler, approximately in the year 1619, and later with theoretical contributions of James Clerk Maxwell in 1873; the first observations and preliminary documentation of the radiation pressure observed in 1901 by E.F. Nichols y G.F. Hull. However, with since of quantum mechanics, with Max Planck and Albert Einstein and with the invention of the laser, that is in 1970 it was possible to obtain the experimental proof where the laser light has the capacity of moving and physically holding microscopic objects by Arthur Ashkin [1]. With the invention of OT, a window been opened for the application it several areas like Atomic Physics, Nanotechnology, Genetics, Biology and Microbiology, because with OT it is possible to individually analyze the microorganisms, allowing to measure mechanical forces and elastic properties [2-3]. In addition to its ability to hold and move microscopic objects without direct contact and exerting small forces around pN. However, conventional OT with Gaussian capture beam and many other limitations: The particles are trapped in the zone of high intensity of the beam and therefore, susceptible to optical damage of the sample (opticution) [4], due to absorption heating, plus of attracting multiple particles in the same trap, requiring the use of diluted samples and lasers with specific wavelengths.
Currently, the use of optical vortex trapping beams has been increasing and many applications have been found [5-9]. Compared to regular Gaussian beams, doughnut-shaped modes, such as TEM*01, are suitable to avoid laser-induced heating and optical damage since the intensity profile drops to zero at the optical axis [6,9].
The prediction of the optical forces that are present in optical trapping is a current problem in continuous evolution. Now, with the use of complex beam like TEM*01 makes the problem even more complex. In this paper, we obtained math expressions for the radiation forces in the Rayleigh’s regime, in order to predict the behavior of a dielectric sphere, under illumination of TEM*01 mode linearly polarized.
2. Thoretical analysis
To explain the behavior of the radiation forces, present in an OT when capturing a micrometer object, there are some methods; in particular, the so-called Approximate Methods, these methods have proven to be useful in the field of OT due to the theoretical and phenomenological analysis they provide. One of these approximate methods is Rayleigh's Dipolar Theory [10], which may be lower in the sense of some approximations, which is smaller than the wavelength of the capture beam. This method consists in consider light as electromagnetic waves and analyzer the behavior of the radiation sphere when it interacts with the capture radiation. As the size of the sphere is considered small, the electromagnetic field incident on the sphere is considered homogeneos, allowing that the dielectric sphere it is behavior as a puntual dielectric dipole.
By electric field action an electric dipole moment is induced in the sphere that has the form [10-12]:
Where n m is the refractive index of the medium, Є0 is the electric permittivity of vacuum, α is the radius of sphere, m=n p /n m is the refractive index of the sphere (n m ) relative to environment (n m ) and E(r,t) is the electric field of the incident radiation. The first force considering is a consequence of the interaction of the electric and magnetic fields of the light beam on the sphere and is exactly the Lorentz's force.
In order to find the mathematical form of the force, it is necessary to calculate the electric and magnetic force on the sphere, taking as an observation point the beam waist and thus obtaining an expression for the force because of the interaction of both fields on the regions positively and negatively charged in the dielectric sphere, adding up both contributions and expressing in terms of electric dipole momento p, it is possible to say that:
Substituting the equation (1) in (2), applying vector identities and because of the electric and magnetic field changes rapidly in time, the temporary average of them is calculated (see Appendix), obtaining:
Where C is the light velocity of vacuum and ∇I(r) is the Gradient of the intensity light beam; thus, the equation (3) is called Gradient Force.
If is have in account the description of the Laguerre beams [14-15], and considering the parameters beam, is possible writing the linearly polarized electric field in direction as;
where r, Ф, z, are cylindrical coordinates, E 0 is the amplitude of electric field and the parameters of the beam: W(Z) is radius of the beam, R(Z) is the radius of curvature and therefore, it is the description of how the wavefront envolves the propagation along the axis, k is the wave number in the medium and ξ(z) is the Gouy's phase change which is the delay of the phase of beam relative to a wave plane [16].
The capture beam that considered is the doughnut-shaped mode type TEM*01. Substituting in the equation (4) the radial index as as p = 0, this is obtained the electric field for this beam, which is related with its radial distribution and azimuth index l = 1 is his helical phase. Later we calculate the intensity of the light beam [17], and replacing in equation (3) we find the shape of the Gradient force in the three normal coordinates, defined in this way
Where, P is the power of the laser beam and W 0 is the radius of the beam waist.
The second force that is takes in account arises from the scattering of the light when the electromagnetic wave incise on the sphere. The dipole oscillates and acts as a source of secondary emission, that emit electromagnetic waves in all directions, due to the oscillating nature of the electric and magnetic fields, and due at size significantly small of the sphere with respect at the wavelenght of the light the dipole oscillates synchronically with the field. As the medium that around of the sphere is considered homogeneous, the resulting propagations is in the direction of the incident wave.
If takes in account the interchanged of the lineal momentum of the beam at the dipole, the transverse section of scattering and the intensity of capture beam, the shaped of the Scattering Force is obtained [10-12]:
Where σ S is denomined cross section of scattering, being this a perpendicular plane to the propagation vector of the scattered wave and for a dielectric sphere that scatters light of isotropic form, has the form [18]:
Taking the form of intensity of the capture beam TEM*01 and cross section of the scattering, equation (9); is possible writing of the Scattering force in function of normalized coordinates as:
3. Results
Through the equations (5), (6), (7) and (10) it is possible to analyze the behavior of the Gradient force and the Scattering force in the capture of a sphere of radius α = 5nm , with refractive index n p = 1,593, which is immersed in a medium of refractive index n m = 1,332, using a Gaussian Laguerre capture beam polarized linearly doughnut-shaped mode type TEM*01 on the axe X with a power P = 100Mw and a wavelength λ = 514,5nm and a beam waist w 0 = 5µm.
The results obtained for the component in x of the Gradient force as is observed in the Fig. 1. In addition to the transverse profile in x of the intensity of the beam TEM*01 (dotted line). Intensity profile has been given an adequate maximum value, in order to give clarity to the behavior of said force in the beam intensity zones.
In the Fig. 1, is possible to observe that for values -1< the Gradient force take positive values, which indicating that it is exerts a force that sends the sphere towards the zone of greater intensity and for values when the sphere it outside of the greater intensity zone but close dark central area of the beam, the Gradient force take negative values, indicating that it exerts a force that sends the sphere towards the zone of greater intensity . The same behavior is evidenced for zones where . Finally, it is possible to conclude that there is a restorative behavior of the transverse nature x of the Gradient force, that is present in the location of the dielectric sphere in the zones of greater intensity of the beam TEM*01. Furthermore, it should be noted that as shown in Fig. 1, the Gradient Force is not present in the dark zone of the beam, because the force depends on the intensity of the beam.
The behavior of the transverse component and of the Gradient Force will be the same that in the component x, as is possible to observe in the equations (5) and (6), because of symmetry. In Fig. 2 it is possible to obtain a view in the transverse plane xy, where the origin of each arrow indicates the position of the sphere, the length of the arrow represents the magnitude of the Gradient force and the direction of the arrow represents the direction of the force respectively. In order to facilitate the analysis, the intensity of the beam has been recorded in the same figure and where it is possible to corroborate the previous restorative behavior in Fig. 1 of the transverse components of the Gradient force. In other words, a potential well is generated in the region of greater intensity of the beam.
In the Fig. 3, is possible to observe the behavior of the transverse component of the Gradient Force on the front plane xz and a cut of the intensity beam in its same plane. In the Fig. 3 is evident its restorative behavior at the zones in x where is to get the greater intensity beam. Nevertheless, is possible to observe that is force is more intensity in the Rayleigh’s length , which is the distance of propagation of beam in where the beam not diverge significantly and is taken since of beam waist.
In Fig. 4, is possible to observe the behavior of the longitudinal component z of the Gradient Force, in a front plane xz. This force positioned at sphere in the greater intensity beam that it is location in the region of beam waist .This force is more intensity in the Rayleigh’s length and is minor in a distance more far of this region. Thus, is possible to observe that on the axes of beam there not is contribution of the longitudinal component of the Gradient Force, because there is absence of radiation in the dark zone of the beam.
In Fig. 5, is observe the behavior of Scattering Force on a sphere that is positioned in (Fig. 7) and in some positions on the transverse axes We observe that the values in the Scattering Force are always positive, permitting that sphere is accelerated in the direction of propagation beam. Is possible too to observe in Fig. 5 that the longitudinal component of the Gradient Force for , take negative values; which permited to transport sphere in opossite direction of propagation, for the zone of greater intensity, which localization is .
This same behavior occur for . Neverthless, is possible to observe that for and for , the Gradient force push to the sphere in the direction of propagation beam. This behavior is because of that the Gradient force has restorative nature at to the zones more near of greater intensity of the beam, and as the sphere is encountered very close of dark zone (without radiation), the intensity zone more close is encountered in the direction propagation of beam. This phenomenon can be observed with more clarity Fig. 4 in and for . 1,5.
In the same Fig. 5, it can be evidenced the longitudinal component Z of the Total force, that is the addition of longitudinal of the Gradient force and of the Scattering force. For , the longitudinal component of the Total force in , have negative values that indicate for this regions that the dielectric sphere is trapped and directed at the beam waist. This behavior is equal for Nevertheless, for and for the sphere is accelerated in the direction propagation of beam, phenomenon that is because of the Gradient force, explained a priori.
A qualitative analysis of the behavior of the Gradient force in the 3D space, can obtained in Fig. 6. In this figure, the hollow hyperboloid represent the zones of greater intensity of the mode TEM*01. Is possible evident the restorative nature of the Gradient force, at to regions of greater intensity. Nevertheless, is possible to observe that arrows that represent at the Gradient force only have components in the plane xy. This is because of that transverse components in x and y(Fig.1), are of order 10-7 while the longitudinal component z, is of order 10-10 (Fig. 5), and therefore, thousand times smaller than the transverse component.
This means that the vectors have an inclination with respect to the plane xy of θ ≈ 0,0405 that can not be appreciated in Fig. 6. As it is possible to predict, the sphere will experience a force that will position it in the regions more intense in the transverse plane xy and that will be much greater than the longitudinal effect.
In more detail, in Fig. 7 which corresponds to a frontal plane xz, it is possible to observe the accelerating behavior in the direction of propagation, which has the Scattering force on the sphere, and its direct proportionality with the intensity of beam. In addition, it should be specified that there is no Scattering force in the dark zone of the beam, which is due to the absence of radiation.
The Fig. 8 correspond to the 3D Scattering force and permitted to observe that for each position of the space where is encountered the sphere this force is responsible of to accelerate to the sphere in the direction of propagation of beam. The maximum magnitude of this force is encountered in the greater intensity zones of beam.
In Fig. 9 it is possible to observe the longitudinal component of the Total force, corresponding to addiction of the Scattering force and the longitudinal component of the Gradient force. For when the Total force change of direction, and have a magnitude that must face at to the Gradient force and scattering force, that is minor that in 0 which both components have the same direction.
The Total force in space will be the Gradient force in transverse plane xy and addition of its longitudinal component with the Scattering force which is possible to see in Fig. 10. In this figure is evident that under the parameters: refraction index, radius of the sphere, wavelength, beam waist and power of beam, is possible to observe a capture of the sphere in the ring of beam that is encountered in beam waist, establishing a potential well in this zone where the greater intensity is found, it should be noted that no Total force is present in the dark zone of the beam, as previously demonstrated. Finally, we could also observe that the arrows that represent the Total force are only on the xy plane, a phenomenon that is due to the difference in the order of the longitudinal component with respect to the transverse component.
4. Interface
To facilitate the analysis in the radiation forces, establishing different parameters, a graphical user interface was developed in MATLAB. This interface implements the equations developed previously with the approximations of Rayleigh Scattering Regime. This method uses equations of electrodynamics to have a theoretical, phenomenological idea and a notion closer what gives rise to the operation of the OT.
The developed interface has two options: 2D Analysis (Fig. 11) and 3D Analysis (Fig. 12), which allow observing the behavior of Gradient, Scattering and Total forces on some transverse and frontal planes, and in the 3D space respectively.
The interface has a panel for inputting parameters by the user which are: refractive index of the medium surrounding the sphere, refractive index of the sphere, radius of the sphere (nm), wavelength of the beam capture (nm), capture beam power in (mW), beam waist(µm).
5. Conclusions
In this work, we presented math expressions for Gradient and Scattering forces exerted on a dielectric sphere under illumination for TEM*01 mode linearly polarized under Rayleigh approach.
Our results under this approximation, can be firstly used by the understanding of behavior in trapping microbead of optical force profile generated by complex beams, e.g. the TEM*01 mode beam.
To show and study the behavior of radiation forces we realized graphical comparisons in the longitudinal, transverse components and vector field. In addition, we developed a GUIDE in MATLAB to study the radiation forces easier when a user changes the trapping parameters.