Photophoretic trapping of multiple particles in tapered-ring optical field

Photophoretic trapping of multiple particles in tapered-ring optical field Fengrui Liu, 1 Zhigang Zhang, 1 Yufeng Wei, 2 Qingchuan Zhang, 1,* Teng Cheng, 1, and Xiaoping Wu 1 1 Key Laboratory of Mechanical Behavior and Design of Materials of Chinese Academy of Sciences, University of Science and Technology of China, Hefei 230027, China 2 Shanghai Aircraft Manufacturing Co., Ltd, Shanghai, 200436, China chteng@ustc.edu.cn * zhangqc@ustc.edu.cn Abstract: We demonstrate the photophoretic trapping of more than several hundreds of absorbing particles by tapered-ring optical traps diffracted from a circular aperture. The experiments with different laser powers show the influence of air flow acting on particles. Three kinds of particles with different densities (about 1~7 g/cm 3 ) and different shapes (spherical, nonspherical) can be trapped. The non-spherical particles (toner particles) disperse in optical field, while the spherical particles (ink droplets and iron particles) arrange as a straight line. More importantly, in the experiments of two counter-propagating tapered-ring beams, the agglomeration of particles is achieved and can help research the dynamics of aerosols. 2014 Optical Society of America OCIS codes: (350.4855) Optical tweezers or optical manipulation; (050.1970) Diffractive optics; (010.1110) Aerosols. References and links 1. A. Ashkin, Acceleration and trapping of particles by radiation pressure, Phys. Rev. Lett. 24(4), 156 159 (1970). 2. M.-C. Zhong, X.-B. Wei, J.-H. Zhou, Z.-Q. Wang, and Y.-M. Li, Trapping red blood cells in living animals using optical tweezers, Nat Commun. 4, 1768 (2013). 3. A. Ashkin and J. M. Dziedzic, Optical trapping and manipulation of viruses and bacteria, Science 235(4795), 1517 1520 (1987). 4. K. Svoboda and S. M. Block, Biological applications of optical forces, Annu. Rev. Biophys. Biomol. Struct. 23(1), 247 285 (1994). 5. H. He, M. E. 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). 6. M. Dienerowitz, M. Mazilu, and K. Dholakia, Optical manipulation of nanoparticles: a review, J. Nanophot. 2(1), 021875 (2008). 7. A. Ashkin, Optical trapping and manipulation of neutral particles using lasers, Proc. Natl. Acad. Sci. U.S.A. 94(10), 4853 4860 (1997). 8. D. G. Grier, A revolution in optical manipulation, Nature 424(6950), 810 816 (2003). 9. E. Higurashi, H. Ukita, H. Tanaka, and O. Ohguchi, Optically induced rotation of anisotropic micro-objects fabricated by surface micromachining, Appl. Phys. Lett. 64(17), 2209 2210 (1994). 10. M. Padgett and R. Bowman, Tweezers with a twist, Nat. Photonics 5(6), 343 348 (2011). 11. K. Dholakia, P. Reece, and M. Gu, Optical micromanipulation, Chem. Soc. Rev. 37(1), 42 55 (2008). 12. A. S. Desyatnikov, V. G. Shvedov, A. V. Rode, W. Z. Krolikowski, and Y. S. Kivshar, Photophoretic manipulation of absorbing aerosol particles with vortex beams: theory versus experiment, Opt. Express 17(10), 8201 8211 (2009). 13. V. G. Shvedov, A. S. Desyatnikov, A. V. Rode, W. Krolikowski, and Yu. S. Kivshar, Optical guiding of absorbing nanoclusters in air, Opt. Express 17(7), 5743 5757 (2009). 14. V. G. Shvedov, C. Hnatovsky, A. V. Rode, and W. Krolikowski, Robust trapping and manipulation of airborne particles with a bottle beam, Opt. Express 19(18), 17350 17356 (2011). 15. F. Liu, Z. Zhang, S. Fu, Y. Wei, T. Cheng, Q. Zhang, and X. Wu, Manipulation of aerosols revolving in taperring optical traps, Opt. Lett. 39(1), 100 103 (2014). 16. V. G. Shvedov, A. V. Rode, Y. V. Izdebskaya, A. S. Desyatnikov, W. Krolikowski, and Yu. S. Kivshar, Giant optical manipulation, Phys. Rev. Lett. 105(11), 118103 (2010). 17. M. Esseling, P. Rose, C. Alpmann, and C. Denz, Photophoretic trampoline Interaction of single airborne absorbing droplets with light, Appl. Phys. Lett. 101(13), 131115 (2012). (C) 2014 OSA 22 September 2014 Vol. 22, No. 19 DOI:10.1364/OE.22.023716 OPTICS EXPRESS 23716

18. Y. L. Pan, S. C. Hill, and M. Coleman, Photophoretic trapping of absorbing particles in air and measurement of their single-particle Raman spectra, Opt. Express 20(5), 5325 5334 (2012). 19. P. Zhang, Z. Zhang, J. Prakash, S. Huang, D. Hernandez, M. Salazar, D. N. Christodoulides, and Z. Chen, Trapping and transporting aerosols with a single optical bottle beam generated by moiré techniques, Opt. Lett. 36(8), 1491 1493 (2011). 20. V. G. Shvedov, A. V. Rode, Y. V. Izdebskaya, A. S. Desyatnikov, W. Krolikowski, and Y. S. Kivshar, Selective trapping of multiple particles by volume speckle field, Opt. Express 18(3), 3137 3142 (2010). 21. J. Lin and Y. Li, Optical trapping and rotation of airborne absorbing particles with a single focused laser beam, Appl. Phys. Lett. 104(10), 101909 (2014). 22. G. Wurm, J. Teiser, and D. Reiss, Greenhouse and thermophoretic effects in dust layers: The missing link for lifting of dust on Mars, Geophys. Res. Lett. 35(10), L10201 (2008). 23. J. B. Wills, K. J. Knox, and J. P. Reid, Optical control and characterisation of aerosol, Chem. Phys. Lett. 481(4-6), 153 165 (2009). 1. Introduction Since the optical manipulation of particles was realized by Ashkin [1] about 40 years ago, it has been applied widely in biology [2 4], physics [5,6] and material science [7 9]. The bestknown example of optical manipulation is trapping and manipulating transparent particles by radiation force [1,10,11]. However, for a light-absorbing particle, the photophoretic force on it is 10 4 times larger than radiation force [12]. It means that the radiation force can be hid by photophoretic force. It is hard for absorbing particles to be manipulated just by radiation force. Fortunately, it has been demonstrated that the photophoretic force can be used to trap and manipulate light-absorbing particles in recent experiments [13 21]. For example, carbon nanoclusters in air can be transported by a dual-vortex optical trap [13]. A single Gaussian beam formed inside the focal volume of a lens with a controlled amount of spherical aberration has been identified that it contains multiple dark traps and can capture airborne absorbing particles [14]. The volume speckle field is used for selective trapping of a few thousand carbon nanoclusters in air [20]. In our previous experiments, a weakly focused beam with tapered-ring optical traps [15] diffracted from a circular aperture was designed to make irregular absorbing particles revolve around optical axis. The laser power is about 50~100 mw. On that condition, the air current with a low speed has weak impact on the particles so that the irregular particles can be driven by tangential photophoretic force to revolve steadily. In this paper, we demonstrate that the tapered-ring optical traps with high laser power (more than 100 mw) can capture more than several hundreds of absorbing particles due to a large number of enclosed dark traps distributed in the optical field. This method provides a novel idea for the optical manipulation of a large number of airborne particles. Rising and descending of many absorbing particles with different laser powers showed the marked influence of air flow acting on particles. The particles with different densities (about 1~7 g/cm 3 ) and different shapes (spherical, non-spherical) could all be optically trapped. More importantly, we institute a two counter-propagating tapered-ring beams system. The agglomeration of particles is achieved and can greatly help research the dynamics of aerosols in air in the future. 2. Simulation and measurement of an optical field In the simulation of a tapered-ring optical field, the incident beam is supposed as a unit amplitude (A = 1) plane wave with λ =532 nm. As shown in Fig. 1(a), the beam passes through both a circular aperture of R = 3 mm diameter and a converging lens of focal distance f = 50mm. A weakly focused optical field forms behind the converging lens. The center of the circular aperture is defined as the origin of the coordinate system. Z axis is parallel to the incident beam. The distance between the converging lens and the circular aperture is very small and can be neglected. The axisymmetrical intensity of the optical field can be given by (derived from the complex amplitude in Ref [15]. (formula (1))) (C) 2014 OSA 22 September 2014 Vol. 22, No. 19 DOI:10.1364/OE.22.023716 OPTICS EXPRESS 23717

2 2 2 2 4π A R jkr 1 1 2π rri 2 2 0 0 (1) I( ri, z) = exp[ ( )] J ( ) rdr λ z 2 z f λz where r is the distance between a point on the diffraction screen and the origin of the coordinate system. The simulated light intensity is shown in Figs. 1(b) and 1(c). Fig. 1. (a) Principle schematic of the simulation of a weakly focused optical field. (b) Simulated 2D light intensity distribution in the y-z plane (x = 0 mm) [15]. (c) Simulated light intensity distribution of one transverse section (z = 48.7 mm). (d) Measured light intensity distribution of one transverse section (z = 48.7 mm). (e) Light intensity distribution of the optical field surrounded by the black square frame in Fig. 1(b). As shown in Fig. 1(b) which shows the 2D light intensity distribution in the y-z plane (x = 0) [15], there are a lot of weakly focused dark stripes distributed in the 2D optical field. Because the optical field is axisymmetric in three-dimensional space, these stripes form multi-layers tapered traps in 3D optical field. Moreover, as shown in Fig. 1(e), a lot of potential barriers (as marked by black dashed ellipses) with high intensity like railings distribute in the passageways of multilayer tapered traps. These barriers separate the tapered traps into a lot of enclosed dark rings. Figures 1(c) and 1(d) show the simulated and measured light intensity distribution of one transverse section (z = 48.7mm). Obviously, the results of experiment and simulation are consistent with each other. The width of each dark ring is about 5~10 micrometers. The particles whose diameters match the dark rings can be captured into these enclosed dark traps when the particles enter the optical field. 3. Reducing laser power To demonstrate that the novel optical traps can capture a large number of light-absorbing particles, we institute the setup as shown in Fig. 2(a). The initial 532nm laser beam waist is about 2 mm and its diameter on the aperture is expanded to about 16 mm. The toner particle, a common light-absorbing particle with the average density about 2 g/cm 3 and the diameter 5~10 μm, is chosen as the sample to be trapped. (C) 2014 OSA 22 September 2014 Vol. 22, No. 19 DOI:10.1364/OE.22.023716 OPTICS EXPRESS 23718

Fig. 2. (a) Principle schematic of the setup trapping particles. L denotes lens (L1, focal length f 1 = 25mm; L2, focal length f 2 = 200 mm; L3, focal length f 3 = 50 mm). C denotes circular aperture with a diameter of 6 mm. (b) (Media 1): The particles fall down from the upper part of the optical field into the lower part when the laser power is reduced. (c) The qualitative analysis about the changing of the particles positions. G denotes the gravity of the particle. S 1 and S 2 denote the upward Stokes drag. P 1 and P 1 denote the photophoretic force. (c1) shows that, with high laser power, the particles are pushed up near the upper part of the optical field because S 1 is larger than G. (c2) shows that, with low power, the particles fall down near the lower part. When laser power is larger than 200 mw, the toner particles are trapped quickly near the upper part of optical field once they are sprayed into the transparent glass cuvette (2cm 2cm 5 cm) as shown in Fig. 2(b1). That is because the air in cuvette absorbs the heat emanated by trapped particles. These particles absorb light and their temperatures increase. Then, the heat is transferred from the particles into air and causes the increasing of the air temperature. The heated air flows upward and produces upward Stokes drag on the toner particles in optical traps. Stokes drag is S = 6πμav, μ is the air viscosity, a is the radius of a particle, v is the relative velocity between air and a particle. In the experiments, some non-trapped particles are lightened when they float across the optical field. We find out the trajectories of these particles by analyzing the recorded pictures and measured the displacement per unit of time. By this method, when the laser power is 200 mw (Fig. 2(b1)), the velocities of 20 toner particles are calculated. To avoid overestimating air velocity, the smallest one (12 mm/s) of the 20 velocities is chosen as the air velocity. Actually, the upward air velocity is larger than the particle s velocity because air needs to overcome particle s gravity if it can drive the particle to rise. The upward Stokes drag on a particle with the radius 5 μm will be 11 5 2 S = 6πμav = 1.96 10 N, where μ = 1.73 10 Nsm. The gravity of the particle is about (C) 2014 OSA 22 September 2014 Vol. 22, No. 19 DOI:10.1364/OE.22.023716 OPTICS EXPRESS 23719

mg 3 6 3 11 = (2.0 10 ) 4 / 3 π (5 10 ) 9.8 = 1. 02 10 N, where the average density of the 3 3 toner particles is ρ = 2.0 10 kg / m. The Stokes drag is larger than the gravity (as analyzed in Fig. 2(c1)) so that the particles are pushed up near the upper part of the optical field. When the laser power is reduced, the temperature of air decreases and then the velocity of air slows down. As shown in Fig. 2(c2), the particles fall down from the upper part to the lower part because the upward Stokes drag decreases that can t offset the gravity G and the photophoretic force P 1. In Fig. 2(b1), with a laser power of 200 mw, almost all of the particles are fixed in the upper part of the optical field (above the white horizontal reference line) and just a few particles are fixed in the lower part (under the reference line). When laser power is reduced, the particles fall down slowly. As shown in Fig. 2(b2), when laser power is reduced to 160 mw, there are about 25 particles falling into the lower part. When laser power is reduced to 80 mw (Fig. 2(b4)), there are only a few particles on the upper part of the optical field. The particles images become blurry from Fig. 2(b1) to Fig. 2(b4) due to the decrease of scatter light from trapped particles when the light intensity decreases. 4. Dynamic behavior of particles with low laser power The number of trapped particles in the optical field is proportional to laser power because the photophoretic force is proportional to the light intensity. When laser power is low, the number of trapped particles is small. Moreover, the particles far from the focus start to bounce slightly. We measured and analyzed the dynamic behavior of trapped particles in the optical field with a low laser power of 30 mw. The light intensity far from the focus is low and the photophoretic force on particles is weak. The impact of air turbulence is relatively strong. Therefore, under air turbulence, only a few particles can be trapped. After about 20 seconds, air turbulence weakens and then more and more particles drifting across the optical field could be trapped easily. Fig. 3. (Media 2) Horizontal views of the moving process of particles. The particles move up gradually into the upper part of optical field. Figure 3 shows the changing of particles positions after the particles are sprayed into the optical field. In Fig. 3(a), some particles are captured in the lower part of optical field after 20 seconds. As time goes on, under irradiation of the laser beam, the temperature of surrounding air increases and the upward Stokes drag becomes stronger. Almost all the particles surrounded by a white ellipse move up into the upper part of optical field under Stokes drag. 5. Trapping of different particles Figure 4 shows the distribution of different kinds of particles observed simultaneously from both horizontal direction (Figs. 4(a1), 4(b1) and 4(c1)) and vertical direction (Figs. 4 (a2), 4(b2) and 4(c2)) with the laser power higher than 400 mw. Figures 4(a1) and 4(a2) show the black toner particles captured by tapered optical traps. Figures 4(b1) and 4(b2) show the trapped spherical black ink droplets with the density about 1 g/cm 3 and the diameter 1~2 μm. Figures 4(c1) and 4(c2) show the spherical iron particles with the density about 7 g/cm 3 and (C) 2014 OSA 22 September 2014 Vol. 22, No. 19 DOI:10.1364/OE.22.023716 OPTICS EXPRESS 23720

the diameter 5 μm. The above experiments show the optical field s ability of trapping different kinds of particles. However, their distribution characteristics are quite different. Almost all the particles are captured in the upper part of optical field (see Figs. 4(a1), 4(b1) and 4(c1)). Observed from vertical direction, the toner particles disperse in the optical field as shown in Fig. 4(a2). In contrast, the trapped black ink droplets and iron particles arrange as a straight line above optical axis (see Figs. 4(b2) and 4(c2)). Fig. 4. (a) - (c): Distribution of three kinds of particles in tapered-ring traps. The left pictures show the results observed from horizontal direction. The right pictures show the results observed simultaneously from vertical direction. (a1) and (a2) show the distribution of toner particles in the optical field. (b1) and (b2) show the distribution of black ink droplets. (c1) and (c2) show the distribution of iron micro-particle. (d) The image of toner particles under the electron microscope. (e) Analysis of the forces acting on trapped toner particles. (f) The image of iron particles under the electron microscope. (g) Analysis of the forces acting on iron particles. We made an analysis (see Figs. 4(e) and 4(g)) for the different distribution characteristics of ink droplets, iron particles and toner particles in the tapered optical traps. As shown in Fig. 4(e), although the upward Stokes drag ( S ) is stronger than the downward gravity (G), the photophoretic force ( P ) composed of PT (tangential component) and PR (radial component) on toner particles is asymmetry and usually has the tangential component ( PT ) because the shape of toner particle is not spherical [15]. Under the tangential component of photophoretic force, the toner particles deviate easily from the medial surface of optical field (see Fig. 4(a2)). The velocity of air in Fig. 4(a) is about 40mm/s, the upward Stokes drag on a toner particle with the average radius a = 5 μm is S = 6πμ av = 6.5 10 11 N. 4 The gravity is mg = πρ a 3 g = 1.0 10 11 N, where ρ = 2.0 103 kg / m3. In Fig. 4(e), the 3 angle θ1 is any angle between 0 and 90. PT and PR can be given by PR = ( S -G ) cos θ1 and PT = ( S -G ) sin θ1. 0, 30, 60 and 90 are chosen as the four representative angles so that we can get a series of PT and PR as shown in Table 1. Table 1. Values of PT and PR when θ1 = 0, 30, 60 and 90 θ1 0 30 60 90 PR PT 5.5 10 11 N 4.8 10 11 N 2.8 10 11 N 0 0 2.8 10 11 N 4.8 10 11 N 5.5 10 11 N #215228 - $15.00 USD Received 11 Jul 2014; revised 6 Sep 2014; accepted 7 Sep 2014; published 19 Sep 2014 (C) 2014 OSA 22 September 2014 Vol. 22, No. 19 DOI:10.1364/OE.22.023716 OPTICS EXPRESS 23721

In contrast, the shapes of iron particle and ink droplet are both spherical. Moreover, the optical field is axisymmetric so that the photophoretic force on them is symmetry (see Fig. 4(e)) and has not the tangential component (only having the radial component ( P R )). The direction of Stokes drag ( S ) is upward so that the spherical particles (iron particles and ink droplets) can be pushed up to the upper part of optical field and arrange in a straight line. Therefore, for ink droplet and iron particle, P= P R =(S-G)cos0 =S-G and P T = 0. In Fig. 4(b), the measured velocity ( v ) of air current is about 35 mm/s, the upward Stokes drag on a 12 ink droplet with the radius a = 0.75 μm will be S = 6πμav = 8.56 10 N. The gravity of the 4 3 14 3 3 ink droplet is mg = πρa g = 1.73 10 N, where ρ = 1.0 10 kg / m. The photophoretic 3 12 force can be express as P = PR = S-G = 8.54 10 N, P T = 0. In Fig. 4(c), the measured velocity ( v ) of air current is about 20 mm/s, the upward Stokes drag on a iron particle with 11 the radius a = 2.5μm will be S = 6πμav = 1.63 10 N. The gravity of the iron particle is 4 3 12 3 3 mg = πρa g = 4.49 10 N, where ρ = 7.0 10 kg / m. The photophoretic force can be 3 11 express as P = P = S-G = 1.18 10 N, P = 0. 6. Dual-beam system R T Based on the above single beam with tapered optical traps, a dual-beam optical traps system is instituted as presented in Fig. 5(a). The counter-propagating beams radiate absorbing particles and generate counter axial photophoretic forces on particles. The particles in optical field are pushed into the junction region of two beams and agglomerate to form micelles. When the toner particles are sprayed into the glass cuvette and pass through the optical field, they gather gradually in the junction region of two beams by photophoretic force (see Figs. 5(b1)-5(b4)), and then agglomerate into micelles. It s easy for about 10-20 particles to agglomerate together. After agglomerating in the junction region, the particles will move as a whole if air turbulence happens. Dozens of experiments are conducted and some micelles often rotate in the dual-beam junction region in almost every experiment. They can keep rotating until outer particles enter the junction region and disturb them heavily. As shown in Fig. 5(c), a micelle composed of more than three particles is rotating clockwise. Its characteristic length is about 25 μm, rotating frequency is 8 Hz and radius of rotation is about 12 μm. Moreover, a particle revolves anticlockwise around a micelle composed of about ten particles as shown in the latter part of Media 3. The revolving frequency is about 0.5Hz, and the revolving track is not circular, but elliptic, semi-major axis of which reaches 60 microns. The particle revolves stably though the tract is not a circular. The rotation of the particles or micelles in the junction region is probably caused by the coupling of the photophoretic force with aerodynamic convection. (C) 2014 OSA 22 September 2014 Vol. 22, No. 19 DOI:10.1364/OE.22.023716 OPTICS EXPRESS 23722

7. Conclusion Fig. 5. (a) Principle schematic of the dual-beam system. (b) Agglomeration process of many particles. (c) (see the first half of Media 3): In the white circle, one micelle is spinning and its rotation plane is parallel to the paper. In addition, the other two micelles are marked by two red ellipses. We demonstrate a novel approach for photophoretic trapping of more than several hundreds of absorbing particles by tapered-ring optical traps with micrometer-size dark rings diffracted from a circular aperture. The experiments with different laser powers show the influence of air flow acting on particles. By changing laser power, we find that, under neither low power nor high power, the heating effect of the laser beam on air cannot be ignored. Therefore, the Stokes drag on particles can produce a demonstrable effect on the distribution of trapped particles. The absorbing particles (toner particle, ink droplet and iron particle) with different densities and different shapes can be trapped. The toner particles disperse in optical field, while the black ink droplets and iron particles are arranged as a straight line above optical axis. Their distribution characteristics are different because the photophoretic forces on particles with different shapes have different directions. The agglomeration of particles is an important direction in the study of aerosols. A counter-propagating beams system is designed to form the junction region, which can make trapped particles agglomerate under the combined effect of dual beams. It allows us to observe the entire process of agglomeration and analysis the dynamic behavior of micelles in real-time. This method provides a novel way for the study of aerosols agglomeration which can be important for greenhouse effect, atmosphere environment [22,23] and even interstellar dusty plasmas. Acknowledgments This work was supported by the National Basic Research Program of China (2011CB302105), the National Natural Science Foundation of China (Grant Nos. 11472266, 11127201, 11102201 and 11332010), and Instrument Developing Project of the Chinese Academy of Sciences (Grant No. YZ201265). (C) 2014 OSA 22 September 2014 Vol. 22, No. 19 DOI:10.1364/OE.22.023716 OPTICS EXPRESS 23723