Dynamic radiation force of a pulsed Gaussian beam acting on a Rayleigh dielectric sphere

We investigate the dynamic evolution of the radiation forces produced by the pulsed Gaussian beams acting on a Rayleigh dielectric sphere. We derive the analytical expressions for the scattering force and all components of the ponderomotive force induced by the pulsed Gaussian beams. Our analysis shows that the radiation force, for both the transverse and longitudinal components, can be greatly enhanced as the pulse duration decreases. It is further found that for the pulse with long pulse duration, it can be used for the stable trapping and manipulating the particle, while for the pulse with short pulse duration it may be used for guiding and moving the small dielectric particle. Finally we discuss the stability conditions of the effective manipulating the particle by the pulsed beam.

Usually optical trapping or tweezers in many experiments are constructed by using the CW laser. It is well known that the CW laser, with the power of a few milliwatts, can only produce the radiation force with an order of a few pN to manipulate the micro-sized particles. Recently, Ambardekar et al. [10] and Deng et al. [11] used a pulsed laser to generate the large gradient force (up to 100pN) within a short duration (~45µs) for overcoming the adhesive interaction between the particles and the surface. Little et al. [12] and Agate et al. [13] have made a comparison between the femtosecond and the CW optical tweezers, and have pointed out that femtosecond optical tweezers are as effective as CW optical tweezers.
However, there is no any investigation about "how does the pulse duration affect on the radiation force acting on a Rayleigh particle?" The radiation forces as electromagnetic forces have been discussed in detail [14,15]. In this paper, we have derived the analytical expressions for the longitudinal and transverse components of the radiation force under the Rayleigh approximation. Based on the derived formula, we analyze the effect of the pulse duration on the physical manipulating impact on the micro-sized particles. At last, we show the stability conditions for effective manipulating the particles.

Radiation force induced by a pulsed Gaussian beam
As shown in Fig. 1, a paraxial pulsed Gaussian beam is radiating on a dielectric sphere. The polarization direction of the pulsed electric field is assumed to be along the x direction. The expression for the electric field of the paraxial pulsed Gaussian beam is given by where 0 w is the beam waist at the plane 0 = z , ρ is the radial coordinate, x is the unit vector of the polarization along the x direction, λ π / 2 = k is the wave number, τ is the pulse duration, and 0 ω is the carrier frequency. Eq. (1) is applicable under the narrow-spectral condition that the pulse spectrum width (proportional to τ / 1 ) is much smaller than 0 ω . For the fixed input energy U of a single pulsed beam, the constant 0 E is determined is the light speed in vacuum, 0 ε and 0 µ are the dielectric constant and permeability in the vacuum, respectively, and 2 n is the refractive index of the surrounding medium. The corresponding magnetic field under the paraxial approximation can be given by are normalized temporal spatial coordinates, and ẑ is the unit vector in the beam propagating direction.
For the sake of simplicity, we assume that the radius a of the particle is much smaller than the wavelength of the pulse, i. e., λ << a . In this case, the Rayleigh approximation is applicable and the dielectric particle can be seen as a point dipole. Under this approximation, for the CW laser, the radiation force includes the scattering force scat F v and gradient force arises from the inhomogeneous field distribution [14,17]; for the pulse, grad F v is a part of the ponderomotive force. In a dilute medium (without dipole-dipole interaction), the ponderomotive force p F v is simply the Lorentz force [14,[16][17]: where is the polarizability of a spherical particle in the Rayleigh regime, and 2 1 / n n m = (here 1 n is the refractive index of the particle). Obviously, for the CW laser, the second term in Eq. (4) is always zero. In our case, for the short pulse, except for grad F v , the second term t F v of Eq. (4) may also play an important role for manipulating the small particles. Substituting Eqs. (1) and (2) into Eq. (4), we can obtain all components of the ponderomotive force as follows: where ρ is the unit vector of the radial direction. For scat F v , it is proportional to light intensity and is along the z + direction [17], is the cross section of the radiation pressure of the spherical particles in the Rayleigh regime. From Eqs. (5) and (6), we find that the magnitude of the radiation forces, especially for the longitudinal components is greatly affected by the pulse duration τ . From Eqs. (5a) and (5b), both the transverse and longitudinal components of the gradient force act as the restoring forces pointed towards the pulse center for the particle with 1 > m although the magnitude of these forces will be changed with different τ . In the following, we find that whether there is a stable trapping strongly depends on the duration τ because t F v is strongly affected by τ . For smaller τ , the second term of Eq. (5c) can be greatly increased in manner of 2 / 1 τ .

Numerical simulation and analysis
Without loss of generality, in all our simulations we choose the following parameters [18]: (for example, the small glass bead and water), nm, and the input pulse power is fixed to be  First, let us see the temporal evolution of the radiation force acting on the small particle. Figures 2(a), 2(b), and 2(c) plot the dynamic evolution of the transverse radiation force. It shows that there is a trapping region in the transverse direction, and the maximal transverse radiation force is enhanced as τ decreases [see Fig. 2(d)]. Figures 2(e), 2(f) and 2(g) plot the dynamic evolution of the longitudinal radiation force, , i.e., all forces along the z axis. As τ becomes small, the maximal | | z F v also increases greatly [see Fig.   2(h)], and the dynamic behavior of z F v changes dramatically from stably to unstably along the z axis. Therefore, for the pulse with large τ , the force is determined by the gradient force, which forms a stable trapping region at the pulse center; in this situation, the gradient force is much larger than scat F v and t F v , thus it is very similar to the case of the CW laser [2,[17][18][19][20]. For the pulse with small τ , although the transverse radiation force provides the restoring force to pull the particle back to the stable point on the z axis, the longitudinal radiation force will make the particle accelerate and decelerate along the z + direction.   Now we turn to discuss how the pulse duration affects on the distribution of the radiation force. Figure 3 shows the dynamic distributions of the transverse and the longitudinal radiation forces acting on the Rayleigh dielectric sphere for the pulse with 1 = τ ps at different times. It is clear that for the pulse with larger τ , the radiation force has a stable trapping region, which is located near the focused center of the pulse and the trapping is enhanced as the pulse arrives at the focused center 0 = z . Therefore in this case the pulse may trap and manipulate the micro-sized particles as effectively as the CW laser [12,13]. With the decreasing of the pulse duration, see Fig. 4, the stable trapping region becomes small and the trapping time becomes short, although the transverse force still provides the restoring force in the transverse direction. For the pulse with very short duration, from Fig. 5, although the pulse cannot be used as an optical trap to manipulate the small particle, it may be used to accelerate the particle because (i) during the first half pulse, the particle will be confined around the z axis and be accelerated along the z + axis, and (ii) during the other half pulse, although the particle will be decelerate it still moves along the z + axis. It also should be noticed that the magnitude of the radiation force are greatly enhanced for the short pulse (also see Fig. 2). In this sense, the pulsed beam may be used to accelerate and levitate the particle for overcoming the van der Walls force from the surface of the vessel [10,11].

Conditions for the stable manipulation
In the above discussion, our analysis shows that the radiation force of the pulsed beam may be used to manipulate the Rayleigh dielectric sphere. In order to know whether the particle could be effectively controlled by the pulse, first we estimate the potential well induced by the radiation force. The potential well must be deep enough to overcome the kinetic energy of the particle due to thermal fluctuation. The condition can be given by using the Boltzmann factor [2,8]: Second, we consider the effect of the random diffusion of the small particle in the surrounding medium. As we know that during the interval of two pulses, the radiation force vanishes, and the particle moves under the gravity and the Brownian force. If the particle doesn't move out the trapping region during the pulse interval, then the radiation force of the second pulse can drag it back to the stable trapping region or manipulate the particle again. In our case, the density of the small glass bead is assumed to be N , which is much smaller than the radiation force, therefore the effect of the gravity can be ignored for the small particle with 5 = a nm. When a particle moves freely inside the ambient, it suffers the Brownian motion. From the Stokes-Einstein relation, the diffusion coefficient is nm. If the pulse repetition rate is larger than 100Hz, the particle diffuses into an area about 12 10 55 . 0 − × m 2 during the pulse interval. This diffusion area is smaller than the acting range of the radiation force. Therefore the particle could be effectively confined within the manipulating region of the pulsed beam.

Conclusion
In summary, we have theoretically investigated the dynamic radiation force produced by the pulsed Gaussian beam acting on a Rayleigh dielectric particle. We find that the radiation force of the pulse can increase greatly as the pulse duration decreases. For the pulsed beam with large pulse duration, it can be used for stable trapping and manipulating the particle; for the pulsed beam with short pulse duration it may be used for guiding and moving the small dielectric particle. Finally the conditions for effective trapping and manipulating the particle have been analyzed.