The radiation force of circular Airy beams (CAB) on a dielectric Rayleigh particle is investigated in this paper. Our results show that the CAB can be used to trap the particle whose refractive index is larger than the ambient at different positions along the beam axis. Comparing with the Gaussian beam under the same conditions, the longitudinal and the transverse gradient force of CAB on the Rayleigh particle are increased, and the particle can be trapped more stable. Our analyses also demonstrate that the trapping properties of CAB can be modulated by controlling corresponding parameters of CAB.
© 2013 Optical Society of America
The circular Airy beam (CAB) has attract intensive interest due to its unique abruptly autofocusing property [1–3]. This beam abruptly focuses its energy right before the focal point while maintaining a low intensity profile until the very point. This feature is very useful in biomedical treatment, for only the intended area would be affected in the propagation direction. Up to now, the CAB with optical vortices [4, 5], and the radially polarized CAB  are proposed and observed in experiment. In optical trapping, the CAB may produce a greater gradient force on the particle, because the abruptly autofocusing property means a great intensity gradient in the focal region. The CAB has been used to trap and guide the micro particles in experiment . But the advantages or the particularities of the optical trapping through CAB have not been detailedly analyzed.
Nowadays, the optical tweezers have been widely applied in many fields, since they are first invented in 1986 by Ashkin . It now becomes an important tool to investigate the biological cells, DNA molecules, neutral atoms, and other particles [9–13]. The conventional optical tweezers are generally constructed with fundamental Gaussian beams (GB). Because of the gradient force of the Gaussian beam, the particle would be stably trapped at the beam waist . Many researches show that other beams also can be used in optical manipulation . For example, the Laguerre Gaussian beam can be used to rotate the particle ; the Bessel Gaussian beam can trap particles in multiple planes ; the Airy beam can be used to transversely clear the particles . Besides, the radiation force of the cylindrical vector beam , Gaussian Schell model beam , Lorentz-Gaussian beam  and other beams [22–24] on the micro particles have also been studied in theory or in experiment. To our best knowledge, the radiation force of CAB on the Rayleigh particle has not been studied before.
In this paper, we investigate the propagation of the CAB under paraxial condition and non-paraxial condition. Because in most situations the light beams must be small enough to trap particles. So it is necessary to know whether the paraxial propagation expression is satisfied when the beam size is very small. The distributions of radiation forces on the Rayleigh particle demonstrate the unique trapping properties of CAB. Similar trapping can also be achieved with generalized abruptly autofocusing beams . By analyzing the optical trap of CAB with different parameters and comparing it with that of the GB, some interesting and useful results are found in our investigations.
2. Non-paraxial propagation and paraxial propagation of CAB
The electric field of the circular airy beam (CAB) at the input plane is defined as :
At first, we’d like to investigate the non-paraxial effect on the abruptly autofocusing property of the CAB. For simplicity, we assume the incident light is linear polarized in the x-direction. The Rayleigh-Sommerfeld diffraction formula can be applied to analyze the propagation of the CAB in the free space. The electric field can be expressed as :27]:26]:28]. The numerical computation of Eq. (3)-(5) is complicated, and some numerical methods could be used [29, 30]. But under paraxial approximation, we can obtain an analytic expression for the propagation of the CAB at the beam center (r = 0) which is useful in calculating the distribution of the radiation force along the beam axis.Eq. (6) can be extended towithout affecting the result. So the electric field in the beam center can be expressed as follows:32],
Now we can compare the non-paraxial results calculated by Eq. (2) with the paraxial results calculated by Eq. (10). In our calculation, we assume a = 0.08, = 1064 nm, r0 = 10s, , the incident power is 1W. The parameter s is usually considered as a proportional coefficient. It does not affect the beam shape or the propagation characteristics of the CAB under the paraxial approximation,as we can see in Fig. 1. As Fig. 1(a) shows, when s = 5, the paraxial result calculated by Eq. (10) and the non-paraxial result calculated by Eq. (2) agree quite closely. The CAB shows abruptly autofocusing property at about, and the intensity reaches its maximum value at. When s is too small, the paraxial result does not agree with the non-paraxial result, as we can see in Fig. 1(b). When s = 1, the non-paraxial result shows that the abruptly autofocusing property of the CAB along the beam axis becomes weak, which indicates that the beam size of the abruptly autofocusing Airy beam is generally limited in practice.
2. Radiation force of CAB exerted on a Rayleigh particle
As is well known, the Rayleigh dielectric particle can be considered as a point dipole in the light fields.And the polarisabilityis :14, 34]:
In our calculation, we choose s = 5. So the paraxial expression of Eq. (10) can be used to calculate the longitudinal radiation force of the CAB at the beam center. However, Eq. (2) should also be used when we want to calculate the transverse radiation force of CAB. We assume R = 20 nm, the refractive index of the surrounding medium nm = 1.33 (i.e., water). Two kinds of particles are discussed below, whose refractive index is np = 1.59 (i.e., glass micro particle) or np = 1.00 (i.e., air bubble), so,.
Figure 2 shows the distributions of the longitudinal and transverse radiation force when np = 1.59. In order to understand the advantages of the CAB over the GB in optical trapping, we assume the initial radius (i.e., r0) and the focal length (i.e., zf, the distance between the initial plane and the waist plane) of the GB are the same as the CAB. The scattering force Fs is directed to the propagation direction, and the gradient force Fg is directed to the equilibrium point in Fig. 2. From Figs. 2(a)-2(c), we can see that the particle could be trapped at the focus point (zf = 975) by the longitudinal gradient force. And there are several stable equilibrium points in the distribution of the sum of the longitudinal gradient force and scattering force (Fig. 2(c)), which indicates that the CAB can trap particles in multiple planes. Typically speaking, for the first two equilibrium points za = 981, zc = 1076, and the zero point between them, zb = 1038, the particle can be longitudinally trapped at za and zc, but cannot be trapped at zb. Note that the position of za is slightly shifted forward from the focus point, because of the influence of the scattering force. Comparing with GB, the longitudinal gradient force of CAB is increased, the scattering force of CAB is decreased. In our case, the GB cannot trap the particle at the beam waist, because the backward gradient force is not large enough to overcome the forward scattering force (Fig. 2(c)). Figures 2(d)-2(f) show that the particle can be transversely trapped at za and zc, but cannot be transverse trapped at zb by CAB. The transverse gradient force decreases as z increases. The gradient force generated from CAB is also larger than that generated from GB at the same position. However, the transverse gradient force cannot attract particles at the beam center when z = zb, which is different from GB. Comparing with GB, the gradient force of CAB is larger and the trapping region of CAB is smaller, as is shown in Figs. 2(a) and 2(d), which indicates that the stiffness of the optical trap formed by CAB is greater. So the particle will be trapped more stably by CAB.
The distribution of the radiation force for the particle with np = 1.00 are shown in Fig. 3. In contrary to the case of np = 1.59, this particle could be longitudinally trapped at zb = 1038, but could not be trapped at za and zc. Although there is equilibrium point at the beam center in Fig. 3(b), yet the gradient force is too small to overcome the Brownian motion, as we will show below. So the particle could not be trapped at zb.
The radiation force can be modulated by controlling corresponding parameters of the CAB. Figure 4(a) shows that when a becomes smaller, the longitudinal radiation force (i.e., the sum of the longitudinal gradient force and the scattering force) for each optical trap will become larger. So the number of effective optical traps formed along the beam axis will be increased. However, the trapping positions are not changed with a. The transverse gradient force also increases with a, as Fig. 4(b) shows. Figure 4(c) shows that the parameter r0 is related with positions of the optical traps, the trapping positions would be moved backward as r0 decreases. The transverse gradient force also increases as r0 decreases (Fig. 4(d)). However, r0 should not be too small, because the CAB possesses the abruptly autofocusing property only when r0 is large enough .
The distribution of radiation force will be stretched or compressed by varying s, because s is a proportional coefficient. When s increases, the positions of the optical traps will be moved forward, the longitudinal radiation force for each optical trap is decreased and the corresponding trapping region is broadened obviously, as Fig. 4(e) shows. The transverse gradient force also increases as s decreases (Fig. 4(f)). However, s should not be too small because the abruptly autofocusing property of CAB will disappear as Fig. 1(b) shows.
3. Analysis of trapping stability
There are several necessary conditions for stably trapping in optical traps by CAB. First the backward longitudinal gradient force must be large enough to overcome the forward scattering force, which is satisfied as we can see in Figs. 2(c), 3(a) and 4. Second, the longitudinal gradient force must be large enough to overcome the influences of the buoyancy and the gravity. In our investigations, the difference between the buoyancy and the gravity for the glass micro particle or the air bubble is about 10−7 pN, which is smaller than maximum value of longitudinal radiation force.Figs. 2 and 3. For the glass particle trapped at za, Rb for the longitudinal gradient force and the transverse gradient force is about and , respectively. For the glass particle trapped at zc, Rb for the longitudinal gradient force and the transverse gradient force is aboutand, respectively. These values are far less than 1. So the glass particle could be stably trapped at za and zc. But for the air bubble trapped at zb, Rb for the longitudinal gradient force is and for the transverse gradient force is 0.59, which approaches 1. So this particle cannot be stably trapped by the transverse gradient force of CAB, because of the influence of the Brownian motion.
In conclusions, we have studied the radiation force of CAB exerted on a dielectric spherical particle in the Rayleigh scattering regime. The numerical results show that the CAB can be used to trap the particle whose refractive index is larger than the surrounding medium. The particle can be trapped at different positions along the beam axis, at the focus or the position some distance behind the focus. Comparing with the conventional optical trap with GB, the longitudinal gradient force and the transverse gradient force are both increased by using the CAB, when the initial radius, the focus length, and the incident power of the two beams are the same. Because of the unique abruptly autofocusing property of CAB, the trap stiffness of the optical trap formed by CAB is larger than that of the GB. For the particle whose refractive index is smaller than the surrounding medium, the longitudinal trapping by CAB is stable, but the transverse trapping is not stable because of the influence of the Brownian motion. The trapping properties of CAB can be controlled by corresponding parameters of the CAB. The value of the longitudinal radiation force can be controlled by varying a or s, the positions of the optical traps along the beam axis can be controlled by varying r0 or s, and the trapping range can be effectively controlled by s. In most cases, the decrease of the beam size is helpful to increase the gradient force. But for the optical trap formed by CAB, the abruptly autofocusing property will disappear when the beam size is too small. We believe our investigation results are useful in optical micromanipulation or other fields.
This work was supported by the National Basic Research Program of China (Grant No.2012CB921602), National Nature Science Foundation of China (Grant No.10974177, 10874012) and the program of International S&T Cooperation of China (Grant No.2010DFA04690).
References and links
3. I. Chremmos, P. Zhang, J. Prakash, N. K. Efremidis, D. N. Christodoulides, and Z. G. Chen, “Fourier-space generation of abruptly autofocusing beams and optical bottle beams,” Opt. Lett. 36(18), 3675–3677 (2011). [CrossRef] [PubMed]
7. P. Zhang, J. Prakash, Z. Zhang, M. S. Mills, N. K. Efremidis, D. N. Christodoulides, and Z. G. Chen, “Trapping and guiding microparticles with morphing autofocusing Airy beams,” Opt. Lett. 36(15), 2883–2885 (2011). [CrossRef] [PubMed]
10. L. Tskhovrebova, J. Trinick, J. A. Sleep, and R. M. Simmons, “Elasticity and unfolding of single molecules of the giant muscle protein titin,” Nature 387(6630), 308–312 (1997). [CrossRef] [PubMed]
13. P. Zemanek and C. J. Foot, “Atomic dipole trap formed by blue detuned strong Gaussian standing wave,” Opt. Commun. 146(1-6), 119–123 (1998). [CrossRef]
14. Y. Harada and T. Asakura, “Radiation forces on a dielectric sphere in the Rayleigh scattering regime,” Opt. Commun. 124(5-6), 529–541 (1996). [CrossRef]
16. H. He, M. E. J. Friese, N. R. Heckenberg, and H. Rubinsztein-Dunlop, “Direct Observation of Transfer of Angular Momentum to Absorptive Particles from a Laser Beam with a Phase Singularity,” Phys. Rev. Lett. 75(5), 826–829 (1995). [CrossRef] [PubMed]
17. V. Garcés-Chávez, D. McGloin, H. Melville, W. Sibbett, and K. Dholakia, “Simultaneous micromanipulation in multiple planes using a self-reconstructing light beam,” Nature 419(6903), 145–147 (2002). [CrossRef] [PubMed]
18. J. Baumgartl, M. Mazilu, and K. Dholakia, “Optically mediated particle clearing using Airy wavepackets,” Nat. Photonics 2(11), 675–678 (2008). [CrossRef]
19. Q. W. Zhan, “Radiation forces on a dielectric sphere produced by highly focused cylindrical vector beams,” J. Opt. A, Pure Appl. Opt. 5(3), 229–232 (2003). [CrossRef]
20. L. G. Wang, C. L. Zhao, L. Q. Wang, X. H. Lu, and S. Y. Zhu, “Effect of spatial coherence on radiation forces acting on a Rayleigh dielectric sphere,” Opt. Lett. 32(11), 1393–1395 (2007). [CrossRef] [PubMed]
22. C. L. Zhao, Y. J. Cai, X. H. Lu, and H. T. Eyyuboğlu, “Radiation force of coherent and partially coherent flat-topped beams on a Rayleigh particle,” Opt. Express 17(3), 1753–1765 (2009). [CrossRef] [PubMed]
24. Z. R. Liu and D. M. Zhao, “Radiation forces acting on a Rayleigh dielectric sphere produced by highly focused elegant Hermite-cosine-Gaussian beams,” Opt. Express 20(3), 2895–2904 (2012). [CrossRef] [PubMed]
26. W. H. Carter, “Electromagnetic Field of a Gaussian Beam with an Elliptical Cross Section,” J. Opt. Soc. Am. 62(10), 1195–1201 (1972). [CrossRef]
27. J. W. Goodman, Introduction to Fourier Optics, (Roberts & Company Publishers, 2005).
28. F. W. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark, The NIST Handbook of Mathematical Functions, (Cambridge University Press, 2010).
30. M. Guizar-Sicairos and J. C. Gutiérrez-Vega, “Computation of quasi-discrete Hankel transforms of integer order for propagating optical wave fields,” J. Opt. Soc. Am. A 21(1), 53–58 (2004). [CrossRef] [PubMed]
31. S. A. Collins Jr., “Lens-System Diffraction Integral Written in Terms of Matrix Optics,” J. Opt. Soc. Am. 60(9), 1168–1177 (1970). [CrossRef]
32. O. Vallée and M. Soares, Airy Functions and Applications to Physics, (Imperial College, 2004).
33. B. T. Draine, “The Discrete-Dipole Approximation and Its Application to Interstellar Graphite Grains,” Astrophys. J. 333, 848–872 (1988). [CrossRef]
34. M. Dienerowitz, M. Mazilu, and K. Dholakia, “Optical manipulation of nanoparticles: a review ,” J. Nanophotonics 2, (2008).