Chinese Optics Letters, 2016, 14 (1): 011405, Published Online: Aug. 6, 2018  

Optically trapping Rayleigh particles by using focused partially coherent multi-rotating elliptical Gaussian beams Download: 971次

Author Affiliations
1 Guangdong Provincial Key Laboratory of Nanophotonic Functional Materials and Devices, South China Normal University, Guangzhou 510631, China
2 CAS Key Laboratory of Geospace Environment, University of Science & Technology of China, Chinese Academy of Sciences, Hefei 230026, China
Abstract
By investigating the cross-spectral density of partially coherent multi-rotating elliptical Gaussian beams (REGBs) that propagate through a focusing optical system, we obtain the radiation force on a Rayleigh particle. The radiation force distribution is studied under different beam indexes, coherence widths, and elliptical ratios of the partially coherent multi REGBs. The transverse and the longitudinal trapping ranges can increase at the focal plane by increasing the beam index or decreasing the coherence width. The range of the trapped particle radii increases as the elliptical ratio increases. Furthermore, we analyze the trapping stability.

In recent years, light beams with rotating characteristics have attracted quite a lot of interest due to their peculiar properties, which present opportunities for research in the fields of science and technology[15" target="_self" style="display: inline;">5]. For example, the ability of rotating beams to rotate particles provides a new degree of control for micrometer-sized particles, and has important applications in optical manipulating and biological specimens rotating. In optics, a more general form of rotating elliptical Gaussian beams (REGBs) are thoroughly studied[1,59" target="_self" style="display: inline;">9]. REGBs, which have an elliptical light spot and whose phase front rotates along the propagating axis of the beam, are obtained from an ordinary Gaussian beam[7].

Ashkin and colleagues demonstrated radiation pressure-trapping particles[1012" target="_self" style="display: inline;">–12] in a series of pioneering Letters, through which the “optical trap” came to be known. The optical force has traditionally been decomposed into two components: one is the scattering force, which is in the direction of the light propagation, and the other is the gradient force, which is in the direction of the spatial light gradient[13]. Optical trapping and manipulation of micrometer-sized particles produced by piconewton-level forces while simultaneously measuring displacement with nanometer-level precision have attracted a very wide range of attention[1418" target="_self" style="display: inline;">18]. Since their invention, optical trapping techniques continue to improve and become better established, and have emerged as a powerful tool with broad-reaching applications in physics and life sciences research[19,20]. Today, many researchers report optical trapping produced by beams such as Gaussian Schell model beams[21], Laguerre–Gaussian beams[22], highly focused Lorentz–Gaussian beams[23], and highly focused, elegant Hermite-cosine-Gaussian beams[24]. In most studies, the incident light beams for optical traps are assumed to be fully coherent[2124]. As we know, any laser field is always partially coherent in practice[2528" target="_self" style="display: inline;">28]. In this Letter, the cross-spectral density (CSD) of partially coherent multi-REGBs propagating through a focusing optical system is derived and used to investigate the radiation forces on a Rayleigh particle. Compared to the work in Ref. [16], the radiation force of the partially coherent multi-REGBs can rotate. To trap oval particles, we can adjust the sizes of the forces in some directions by changing the elliptical ratio of the elliptical beams, which ordinary Gaussian beams cannot do. Furthermore, the trapping stability conditions are analyzed.

Under the free-space rectangular coordinates system, the electric field of the REGB in the z=0 plane can be expressed as[1,5]E0(x,y)=A0exp[(x2wx2+y2wy2+ixyR0)],where A0 is the constant complex amplitude of E0, wx and wy are the initial beam widths on the x- and y-axes, respectively, and R0 specifies the beam rotation. Without loss of generality, the spectral degree of coherence is given by[29]g(r1r2)=exp[(r1r2)2δ2],where r1(x1,y1) and r2(x2,y2) are the original position vectors, and δ is the coherence length. By taking all of the above assumptions into consideration, the CSD of the partially coherent multi-REGBs of the source plane can be written in the Cartesian coordinates as[25]W0(x1,y1,x2,y2;z=0)=1C0m=1M(1)m1m(Mm)E1*E2g(r1r2)=A02C0m=1M(1)m1m(Mm)exp[(x12wx2+y12wy2+ix1y1R0)(x22wx2+y22wy2ix2y2R0)(x12+x222x1x2+y12+y222y1y22mδ2)],where C0=m=1M[(1)m1/m](Mm) is the normalization factor with (Mm) representing the binomial coefficients, and * denotes the complex conjugate.

Considering the propagation of the partially coherent multi-REGBs through an unapertured lens system (see Fig. 1), where z is the distance from the input plane to the output plane, f is the focal length, and z1 is the axial distance from the focal plane to the output plane, the ABCD matrix for the focusing system is (ABCD)=(1z01)(101/f1)=(1z/fz1/f1).

Then, the expression for the CSD of the electric field in the output plane is given by[30]W(x1,y1,x2,y2;z)=(k2πB)2dx1dx2dy1dy2W0×exp{ik2B[(Ax222x2x2+Dx22)+(Ay222y2y2+Dy22)(Ax122x1x1+Dx12)(Ay122y1y1+Dy12)]},where λ is the wavelength of the incident light, k=2π/λ is the wavenumber in free space, and r1(x1,y1) and r2(x2,y2) are the position vectors in the output plane. If we substitute Eqs. (3) and (4) into Eq. (5), performing the related integration, and assume that r1=r2=r(x1=x2=x,y1=y2=y), then the CSD reduces to W(r,r;z)=A02Δ1C0m=1M(1)m1m(Mm)exp(Y1),Δ=(A2+4B2k2X1)(A2+4B2k2X2),X1=(1wx2)2+1mδ2wx2(i2R0)2,X2=(1wy2)2+1mδ2wy2(i2R0)2,Y1=2Δ2[ϵxx2wx22BAkR0(1wx2+1wy2)xy+ϵyy2wy2],ϵx=A2+4B2k2[(1wy2)2+1mδ2wy2wx2wy2(i2R0)2],ϵy=A2+4B2k2[(1wx2)2+1mδ2wx2wy2wx2(i2R0)2].

The intensity distribution and the CSD of the partially coherent multi-REGBs at any point of the output plane are related by the formulas[25]Iout(r;z)=W(r,r;z).

In this Letter, we choose f=0.1mm, λ=1064nm, wx=0.6mm, wy=0.1mm, R0=1m2, A02=1W/mm2, and δ=0.1mm. Figure 2 shows that the beam profile becomes flatter, but the maximum intensity decreases as the beam index M increases (for example, M=10, M=30), which means the potential wells are lower or shallower than their counterparts (for example M=1, M=5).

Fig. 1. Schematic of a focusing optical system. a is the radius of the particle, np is the refractive index of the particle, and nm is the refractive index of the ambient.

下载图片 查看所有图片

Fig. 2. (a) Intensity distribution for different values of M at the focal plane. (b) Two-dimensional intensity distribution for M=30.

下载图片 查看所有图片

We study the radiation force of partially coherent multi-REGBs on a Rayleigh dielectric particle, which is sufficiently small compared with the wavelength (i.e., the particle radius aλ/20) and can be seen as a point dipole in the light field. The scattering force and the gradient force are given by[13,21]F⃗scat(r⃗;z)=e⃗znmα0Iout/c,F⃗grad(r⃗;z)=2πnmβ0Iout/c,where e⃗z is a unity vector along the beam propagation, and α0=(128π5a6/3λ4)[(m021)/(m02+2)]2, β0=a3(m021)/(m02+2), m0=np/nm is the relative refraction index of the particle. We choose a=40nm, np=1.8, and nm=1.33 (water). By substituting Eq. (13) into Eqs. (14) and (15), we can calculate the radiation force induced by partially coherent multi-REGBs on a Rayleigh particle. Both the scattering force and the gradient force will be greatly affected by the beam index M, the coherence width δ, and the elliptical ratio γ (γ=wx/wy); Figs. 35 show this numerically.

Fig. 3. Effect of the beam index M on (a)–(c) Fscat (the scattering force), (d)–(f) Fgrad,t (the transverse gradient force). Fgrad,t=(Fgrad,x2+Fgrad,y2)1/2 while z1<0; Fgrad,t=(Fgrad,x2+Fgrad,y2)1/2 while z1>0. (g) and (h) The effect of the beam index on Fgrad,z (the longitudinal gradient force). The parameters are chosen as y=0 and δ=0.1mm.

下载图片 查看所有图片

Fig. 4. Effect of the coherence width δ on (a)–(c) Fscat, (d)–(f) Fgrad,t, and (g) and (h) Fgrad,z. The parameters are chosen as y=0, M=30.

下载图片 查看所有图片

Fig. 5. Effect of the elliptical ratio γ on (a)–(c) Fscat, (d)–(f) Fgrad,t, and (g) and (h) Fgrad,z. The parameters are chosen as y=0, M=30, δ=0.1mm, and wy=0.1mm.

下载图片 查看所有图片

Figure 3 shows the influence of the beam index M on the radiation force produced by the partially coherent multi-REGBs. From Figs. 3(a), 3(d), and 3(g), we can see that a Rayleigh particle whose refractive index is larger than that of the ambient can be trapped at the focus point by the partially coherent multi-REGBs, because there are stable equilibrium points in Figs. 3(d) and 3(g). From Fig. 3(h), we can see that when M=1, the partially coherent multi-REGBs could not stably trap the particle when x=0.3μm, because the direction of the longitudinal gradient force is same as the scattering force direction. But for larger M (M=5,10,30), the direction of the longitudinal gradient force is against to the direction of the scattering force. So we can conclude that by increasing the beam index M of the partially coherent multi-REGBs, we can increase the transverse trapping range at the focal plane. By comparing the figures in Fig. 3, one can find that both the gradient force and the scattering force decrease as the initial beam index M increases, which means that the trapping stability changes.

Figure 4 shows the influence of the coherence width δ on the radiation force produced by the partially coherent multi-REGBs. From Figs. 4(d) and 4(g), we find that at the focal plane, as the coherence width of the partially coherent multi-REGBs decreases, the positions of the peak values deviate away from the focus point. This means that both the transverse and the longitudinal trapping ranges become larger. However, as is shown in Fig. 4(h), the direction of the longitudinal gradient force will reverse when the coherence width increases to 0.2 mm at x=0.3μm, and the beam will lose the ability to trap the particle. By decreasing the coherence width δ of the partially coherent multi-REGBs, we can increase the trapping ranges. One can see from Fig. 4 that the scattering force and the transverse and the longitudinal gradient forces decrease as the coherence width of the initial partially coherent multi-REGBs decreases. Partially coherent beams have their advantages when trapping.

Figure 5 shows the influence of the elliptical ratio γ on the radiation force produced by the partially coherent multi-REGBs. From Fig. 5, we can see that both the gradient force and the scattering force increase as the elliptical ratio γ increases, which means that the trapping stiffness is higher. The positions of the peak value deviate away from the focus point as γ increases, as shown in the Figs. 5(e) and 5(f). As γ increases, the trapping range broadens, demonstrating the advantage of using elliptical beams. From Figs. 5(a), 5(d), and 5(g), we can see that a Rayleigh particle can be trapped at the focus point by the partially coherent multi-REGBs because there are stable equilibrium points in Figs. 5(d) and 5(g) for different values of γ. For the oval particles, we can adjust the size of forces in all directions by changing the elliptical ratio of elliptical beams, which ordinary Gaussian beams (γ=1) cannot do. So when we choose different values of γ, we can get the different trapping ranges, which is beneficial for the application.

For stable trapping, the gradient force must be large enough to overcome the scattering force R=|Fgrad,z|/|Fscat|1, where the ratio R is called the stability criterion. As is shown in Figs. 35, the scattering force is much smaller than the longitudinal gradient force near the focal plane. Figure 6 shows the two-dimensional force distribution of the partially coherent multi-REGBs near the focus point. In this process, the direction of the scattering force is always along the z direction (white points indicate the +z direction and white forks indicate the z direction), and the direction of the transverse gradient force always points to the center of the beam (arrows). But the direction of the longitudinal force will change when the beams pass through the focus plane. As the direction of the transverse force always points to the center of the beam, the particle can be trapped stably in transverse direction. When z1>0, the particle can be trapped stably in longitudinal direction, as the direction of the scattering force and the center of the longitudinal force are opposite. The partially coherent multi-REGBs rotate while propagating, and the force scope changes. The rotational nature of these beams is interesting for the rotating trapped particle.

Fig. 6. Two-dimensional force distribution of the REGBs near the focus point. (a1)–(a5) are the scattering force distributions. (b1)–(b5) are the transverse gradient force distributions. (c1)–(c5) are the longitudinal gradient forces distributions at positions z1=0.3, 0.1, 0, 0.1, and 0.3 mm.

下载图片 查看所有图片

The Brownian force, which describes the influence of the Brownian motion, can be calculated by the expression[31]|Fb|=12πηaKBT,where η is the viscosity of the ambient, which is 8.0×104Pa·s at T=300K, a is the radius of the Rayleigh particle, and KB is the Boltzmann constant.

We plot the change of the magnitude of all the forces in Fig. 7 with x=Δ/(2ϵx) and f=f+0.1μm, where Fgrad,tm is maximum transverse gradient force, Fgrad,zm is maximum longitudinal gradient force, and Fscatm is maximum scattering force. In Fig. 7(a), the longitudinal gradient force is much bigger than the Brownian motion. In Fig. 7(b), it is found that the particle can be stably trapped when δ>0.005mm, but if δ is too small, the Brownian force will be too strong. In Fig. 7(c), the longitudinal gradient force will be bigger than the transverse gradient force while γ>4, and the longitudinal gradient force is always bigger than the Brownian force. In Fig. 7(d), it is found that the particle can be stably trapped while 14nm<a<50nm, but when a is too small, the disturbance is mainly from the Brownian motion. However, when a is greater than 45 nm, the disturbance of the Brownian force equals the scattering force. The particle can be stably trapped while 46nm<a<50nm for γ=1, 33nm<a<50nm for γ=2, and 15nm<a<50nm for γ=8. So we can conclude that as γ increases, the radius range of the trapped particle increases.

Fig. 7. When z1=0.1μm, the radiation forces Fscatm, Fgrad,tm, Fgrad,zm produced by partially coherent multi-REGBs depend on (a) the initial beam index M with δ=0.1mm and γ=6, (b) the coherence width δ with M=30 and γ=6, (c) the elliptical ratio γ with M=30 and δ=0.1mm, and (d) the particle radius a with M=30, δ=0.1mm, and γ=6. Fb is the Brownian force.

下载图片 查看所有图片

In conclusion, we investigate the CSD of the partially coherent multi-REGBs that propagate through a focusing optical system, and study the radiation force on a Rayleigh dielectric sphere induced by focused partially coherent multi-REGBs. The partially coherent multi-REGBs rotate while propagating, and the force scope changes. The rotational nature of these beams is interesting for rotating trapped particles. We find that we can increase the transverse and the longitudinal trapping ranges at the focal plane by increasing the beam index M, or by decreasing the correlation width δ. As the elliptical ratio γ increases, the radius range of the trapped particle increases. For oval particles, we can adjust the size of forces in some directions by changing the elliptical ratio of the elliptical beams, which ordinary Gaussian beams cannot do. When comparing the Gaussian beams, it is clear that the partially coherent multi-REGBs have better applications in operating flexibility.

References

[1] GoncharenkoA.LogvinY.SamsonA.ShapovalovP., Opt. Commun.81, 225 (1991).OPCOB80030-4018

[2] DengD.PengX.ChenC.ChenB., J. Korean Phys. Soc.66, 774 (2015).KPSJAS0374-4884

[3] DengD., Opt. Commun.285, 3976 (2012).OPCOB80030-4018

[4] AlievaT.AbramochkinE.Asenjo-GarciaA.RazuevaE., Opt. Express18, 3568 (2010).OPEXFF1094-4087

[5] PengX.ChenC.ChenBoDengD., Chin. Phys. B24, 124201 (2015).1674-1056

[6] ArnaudJ.KogelnikH., Appl. Opt.8, 1687 (1969).APOPAI0003-6935

[7] GoncharenkoA., Gaussian Light Beams (Nauka i tehnika, 1977).

[8] LiuD.WangH.WangY.YinH., Opt. Laser Technol73, 12 (2015).

[9] DengD.GuoQ., Phys. Rev. E84, 046604 (2011).

[10] AshkinA.DziedzicJ.BjorkholmJ.ChuS., Opt. Lett.11, 288 (1986).OPLEDP0146-9592

[11] AshkinA.DziedzicJ.YamaneT., Nature330, 769 (1987).

[12] AshkinA., IEEE J. Sel. Top. Quant.6, 841 (2000).IJSQEN1077-260X

[13] HaradaY.AsakuraT., Opt. Commun.124, 529 (1996).OPCOB80030-4018

[14] ShuJ.ChenY.PuJ., J. Opt. Soc. Am. A30, 916 (2013).JOAOD60740-3232

[15] WangS.ChanC., Nat. Commun.5, 1 (2014).NCAOBW2041-1723

[16] LiuX.ZhaoD., Appl. Opt.53, 3976 (2014).APOPAI0003-6935

[17] YuH.SheW., J. Opt. Soc. Am. A32, 90 (2015).JOAOD60740-3232

[18] ZhongM.WangX.ZhouJ.WangZ.LiY., Chin. Opt. Lett.12, 011403 (2014).CJOEE31671-7694

[19] AshkinA.DziedzicJ., Science235, 1517 (1987).SCIEAS0036-8075

[20] LiphardtJ.OnoaB.SmithS.TinocoI.BustamanteC., Science292, 733 (2001).SCIEAS0036-8075

[21] WangL.ZhaoC.WangL.LuX.ZhuS., Opt. Lett.32, 1393 (2007).OPLEDP0146-9592

[22] BhattacharyaM.MeystreP., Phys. Rev. Lett.99, 153603 (2007).PRLTAO0031-9007

[23] JiangY.HuangK.LuX., Opt. Express19, 9708 (2011).OPEXFF1094-4087

[24] LiuZ.ZhaoD., Opt. Express20, 2895 (2012).OPEXFF1094-4087

[25] MandelL.WolfE., Optical Coherence and Quantum Optics (Cambridge University Press, 1995).

[26] ZhaoC.CaiY.LuX.EyyuboğluH., Opt. Express17, 1753 (2009).OPEXFF1094-4087

[27] GoriF.SantarsieroM., Opt. Lett.40, 1587 (2015).OPLEDP0146-9592

[28] ZhangG.LuX.ZhouS.MaW.ZhuJ., Chin. Opt. Lett.12, S11401 (2014).CJOEE31671-7694

[29] WolfE., Introduction to the Theory of Coherence and Polarization of Light (Cambridge University Press, 2007).

[30] SiegmanA., Lasers (University Science Books, 1986).

[31] OkamotoK.KawataS., Phys. Rev. Lett.83, 4534 (1999).PRLTAO0031-9007

Xi Peng, Chidao Chen, Bo Chen, Yulian Peng, Meiling Zhou, Xiangbo Yang, Dongmei Deng. Optically trapping Rayleigh particles by using focused partially coherent multi-rotating elliptical Gaussian beams[J]. Chinese Optics Letters, 2016, 14(1): 011405.

引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!