## Abstract

We optically manipulated a metal particle in borosilicate glass. The glass in the neighborhood of the laser-heated metal particle softened; hence, the metal particle was able to migrate in the glass. In this letter, the driving force of the metal particle toward the light source in the glass provided by laser illumination was investigated. The variation in the surface tension of the glass at the interface between the glass and the metal particle induced by the temperature gradient was calculated via a numerical temperature calculation. It was found that the temperature at the laser-illuminated surface of a stainless-steel particle with a radius of 40 μm was ~320 K higher than that on the nonilluminated side. The force applied to the metal particle from the surrounding glass was calculated to be ~100 μN, which was approximately equal to the viscous resistance force. In addition, the experimental and numerically calculated speeds of the moving particle, which was measured while varying the laser power, are discussed.

© 2013 OSA

## Corrections

Hirofumi Hidai, Makoto Matsushita, Souta Matsusaka, Akira Chiba, and Noboru Morita, "Moving force of metal particle migration induced by laser irradiation in borosilicate glass: erratum," Opt. Express**22**, 25194-25195 (2014)

https://www.osapublishing.org/oe/abstract.cfm?uri=oe-22-21-25194

## 1. Introduction

Surface tension generally decreases with increasing temperature [1]. The variation of surface tension at a droplet surface under a temperature gradient induces a driving force [2]. Using the driving force, the manipulation of a droplet has been demonstrated by local heating with a laser [1]. Light can move matter by other types of driving force, for example, optical tweezers [3], a photophoretic force [4], fluid flow generated by local melting, and a volume change due to a phase change [5].

We optically manipulated a metal particle in borosilicate glass [6]. The glass in the neighborhood of the laser-heated metal particle softened; hence, the metal particle was able to migrate in the softened glass. The metal particle was obtained by melting a metal foil deposited on the back of the glass by illuminating the foil with a laser through the glass. However, the driving force of particle migration has not been identified.

In this report, we reveal the driving force of metal particle migration toward the light source. The variation in the surface tension of the glass at the interface between the glass and the metal particle induced by the temperature gradient was calculated via a numerical temperature calculation and compared with the moving speed obtained experimentally.

## 2. Numerical calculation

First, the temperature distribution of the melted metal particle and the surrounding glass was calculated numerically. Second, the surface tension of the metal particle was estimated from the temperature distribution. Third, the total pressure applied by the surface tension was clarified to be the moving force.

Figure 1 shows a schematic image of the simulation model. Two-dimensional axisymmetric spherical coordinates were used because both the laser beam and the metal particle were symmetrical about the optical axis (z-axis in Fig. 1). The flow around and within a spherical fluid have been derived by Hadamard and Rybczynski under the assumption that the inertia terms are set identically equal to zero everywhere in the flow field [7].

Here, ${v}_{r}$ and ${v}_{\theta}$ are the flow velocities in the radial and circumferential directions, respectively, ${v}_{0}$ is the velocity of the metal particle, ${r}_{0}$ is the radius of the metal particle, and $\mu $ and $\mu \text{'}$ are the viscosities of the surrounding glass and metal, respectively.

Thermal radiation from the surface of the metal particle and heat conduction were considered. The temperature $T(t,r,\theta )$ at time *t* at location $(r,\theta )$ was calculated utilizing the following two-dimensional heat conduction equation in the spherical coordinate system:

*c*is the specific heat, $\rho $ is the density, $k$ is the thermal conductivity, ${Q}_{las}(t,r,\theta )$ is the heat supplied by the laser, and ${Q}_{rad}(t,r,\theta )$ is the heat loss by radiation.

Regarding the laser parameters, a cylindrical Gaussian energy distribution was set. We assumed that the laser energy was completely absorbed on the surface because the absorption depths of metals are small. The heat supplied by the laser is expressed as

where *R* is the reflectivity of the metal particle. $I\left(r,\theta \right)$ is the laser intensity, expressed as

where *P* is the laser power and *w* is the spot radius. The heat loss of radiative heat transfer is expressed as

*ε*is the emissivity,

*β*is the Stefan-Boltzmann constant, and ${T}_{0}$ is the environmental temperature.

The initial and boundary conditions were set as

where ${T}_{0}$ is the environmental temperature and ${r}_{ib}$ is the value of $r$ at the boundary of the calculated region 160 μm from the center of the metal particle, where constant-temperature conditions were applied. Heat Eqs. (5)-(8) were solved by an explicit finite difference method, and flow Eqs. (1)-(4) were included in the calculation using the following values: ${v}_{0}=50$ μm/s, ${r}_{0}=40$ μm, $\mu =8.6\times {10}^{-3}$ Pa∙s, $\mu \text{'}=6.0\times 10$ Pa∙s, $R=0.6$, $w=\text{150}$μm,${T}_{0}=293$ K, $\epsilon =0.4$, $\beta =5.67\times {10}^{-8}$ W/(m^{2}∙K

^{4}), $\Delta r=2$ μm, and$\Delta \theta =7.5$°. The properties of the glass and melted metal are $c=1.15\times {10}^{3}$ J/(kg∙K) and $7.24\times {10}^{2}$ J/(kg∙K);$\rho =2.2\times {10}^{3}$ kg/m

^{3}and $7.93\times {10}^{3}$ kg/m

^{3}; and $k=2$ W/(m∙K) and 49 W/(m∙K), respectively. ${v}_{0}$ is set to 50 μm/s in this calculation under all conditions. The temperature increased upon laser illumination and then became constant after a certain time.

The steady temperature $T\left({t}_{c},{r}_{0},\theta \right)$ at the surface of the metal particle was obtained by the above calculation. The surface tension of the glass $\sigma \left(T\right)$ varies in accordance with the temperature. The pressure applied on the metal particle by the surface tension of the glass was calculated using the Young–Laplace equation

Therefore, the total applied force ${F}_{sur}$ on the metal particle from the surrounding glass with a temperature distribution was expressed asOn the other hand, the viscous resistance force ${F}_{res}$ generated by the movement of the metal particle was expressed under the assumption that the viscosity of the glass was constant.

Here,*V*is the velocity of particle migration. The temperature dependences of the surface tension $\sigma \left(T\right)$ N/m and viscosity $\mu \left(T\right)$ Pa∙s were approximated asandfrom the data for Glass No. 5 in Fig. 7 in [8], and Fig. 1 in [9], respectively.

The velocity of the particle was calculated using

where *m* is the mass of the particle.

## 3. Experimental procedure

The laser illumination setup and side-view microscope system used to observe particle migration were similar to those described in [6]. Briefly, borosilicate glass with a thickness of 10 mm (Pyrex®, Corning 7440, Corning Inc.) was used. SUS304 foil (753173 Nilaco Corp.) with a thickness of 10 μm was sandwiched between the glass sample and another glass plate. A Nd:YAG laser beam (Awave1064-50W50K, Advanced Optowave Corp.) with a wavelength of 1064 nm and continuous wave oscillation was focused on the foil through the sample glass by a convex lens with a focal length of 170 mm. The beam mode of the laser was TEM_{00} (M^{2} < 1.3). The theoretical beam radius was calculated to be ~150 μm from the specifications of the laser oscillator by assuming that M^{2} was 1.3. The laser power was varied in the range of 17.3 – 26.7 W. The laser focus was fixed during laser irradiation. Side-view shadowgraph images under white-light illumination were obtained to monitor the process in situ. Scattered laser light and thermal emission were filtered by a band-pass filter (10BPF10-440, Newport) placed in front of the CCD camera used to obtain images of the particle.

## 4. Results and discussion

After the start of laser irradiation, a black particle moved backward (toward the light source) from the SUS304 foil. The particle migration stopped when laser illumination ceased. The particle started moving again when the laser illumination was restarted, as shown in Media 1. Figure 2 shows snapshots from Media 1 taken 8 s (a) and 13 s (b) after the start of laser illumination. The particle moved ~360 μm.

The black particle was investigated with an energy-dispersive X-ray spectroscopy system (EDX, Link ISIS 300, Oxford Instruments Plc.) installed with a field-emission scanning electron microscopy system (FE-SEM, JSM-6340, JEOL Ltd.) after cutting to expose the particle. As a result, the elements comprising stainless steel, *i.e.,* iron, nickel, and chrome, were detected from the particle. Therefore, the black particle was an implanted stainless-steel particle.

The trajectory of metal particle migration shown in Fig. 2 could be observed because heating and quenching by laser illumination caused a change in the refractive index [6]. It is noteworthy that the melting of the stainless-steel particle was confirmed, because two particles spontaneously fused upon contact with one another.

The temporal temperature change at the hottest point (*r* = 40 μm and *θ* = 0°) was calculated and is shown in Fig. 3. The temperature increases rapidly until ~0.02 s, then gradually increases, then becomes constant ~0.17 s after the start of laser illumination. Figure 4(a) shows the distribution of the temperature at a power of 21.3 W after the temperature became constant. Arrows show the flow within and around the particle. The color scale indicates the temperature. Figure 4(b) shows the temperature on the surface of the particle. The temperature at$\theta =0$° (with the optical axis corresponding to the normal line) is the highest (~2270 K) and the temperature on the other side ($\theta =180$°) is the lowest (~1950 K) on the surface of the particle. The temperature difference on the surface is ~320 K.

The force applied to the metal particle from the surrounding glass ${F}_{sur}$ was calculated using Eqs. (10) and (12) to be ~100 μN. The temperature and viscosity of the surrounding glass (at *r* = 54 μm and *θ* = 90°) were calculated to be 1440 K and 3.8 × 10^{3} Pa∙s, respectively, from Eq. (13). The migration speed of the particle was calculated to be 34 μm/s by assuming the viscous resistance force ${F}_{res}$ to be equal to the force applied from the surrounding glass ${F}_{sur}$. The experimental particle speed was measured to be ~72 μm/s and corresponded reasonably well with the calculated speed (34 μm/s).

Figure 5 shows the speed of the moving particle after the speed became constant, which was measured for various laser powers, and the calculated speed obtained by the method described above. The migration speed increased with the power, and the maximum speed was ~100 μm/s. When the power was higher than 26.7 W, the glass itself absorbed the light, a phenomenon that was also reported in [10, 11], and the particle did not migrate. The calculated particle speed also increases with increasing laser power, as shown in Fig. 5. However, the variation was larger than that in the experimental results, and the maximum speed was calculated to be 700 μm/s at a power of 26.7 W. The reason for this was considered to be as follows. The viscous resistance force was calculated under the assumption that the particle migrated in a medium with constant viscosity. However, the viscosity increased with increasing distance from the particle, hence only the glass around the particle was fluid. The viscous resistance force was calculated to be smaller than the actual force, particularly when particle speed was high. Therefore, the calculated particle speed was higher than the actual speed.

Figure 6 shows the velocity of the particle, which was measured by plotting the position of the particle in Media 1 every second. The particle starts moving ~1 s after the start of illumination, the velocity increases gradually, and then the velocity becomes constant ~8 s after the start of illumination.

The speed of the particle was calculated using Eq. (14) by assuming that the stationary particle was accelerated immediately upon laser illumination with a constant force of 100 μN. The speed of the particle exceeded 90% of its terminal speed in less than 10 ns. Therefore, the speed of the particle became constant within ~0.17 s after the start of laser illumination, considering that the temperature distribution became constant after ~0.17 s as mentioned before.

The calculated duration of the laser illumination after which a particle reached a constant speed was shorter than the experimental result. This is considered to be because the calculated duration after which a constant temperature was reached was much shorter than the experimental duration, since the dependence of the material properties was not considered. In particular, the reflectivity of the solid metal is greater than that of the melted metal; therefore, the calculated temperature increased faster than the experimental temperature.

Finally, the migration force was considered. As discussed in [6], the illumination of light generates other types of force. The force generated by the surface tension gradient was calculated to be ~100 μN. In contrast, for example, the dipole force is on the order of nano-Newtons [3]. Therefore, the force generated by surface tension gradient was much larger than the dipole force; thus, the force generated by the surface tension is believed to be the dominant force causing the particle migration in glass.

## 5. Conclusion

The driving force of a metal particle toward a light source in glass upon laser illumination was investigated. The variation in the surface tension of the glass at the interface, induced by the temperature gradient, was calculated via a numerical temperature calculation. As a result, the temperature difference on the surface of a stainless-steel particle with a radius of 40 μm was calculated to be ~320 K. The force applied to the metal particle from the surrounding glass was calculated to be ~100 μN, which was approximately equal to the viscous resistance force. In addition, the experimental and numerically calculated speeds of the moving particle were discussed. Surface tension was the dominant force driving the particle toward the light source.

## Acknowledgment

Support from Japan Society for the Promotion of Science through a Grant-in-Aid for Scientific Research (24656096) is gratefully acknowledged.

## References and links

**1. **K. T. Kotz, K. A. Noble, and G. W. Faris, “Optical microfluidics,” Appl. Phys. Lett. **85**(13), 2658–2660 (2004). [CrossRef]

**2. **N. O. Young, J. S. Goldstein, and M. J. Block, “The motion of bubbles in a vertical temperature gradient,” J. Fluid Mech. **6**(03), 350–356 (1959). [CrossRef]

**3. **A. Ashkin, “Acceleration and trapping of particles by radiation pressure,” Phys. Rev. Lett. **24**(4), 156–159 (1970). [CrossRef]

**4. **M. Sitarski and M. Kerker, “Monte Carlo simulation of photophoresis of submicron aerosol particles,” J. Atmos. Sci. **41**(14), 2250–2262 (1984). [CrossRef]

**5. **F. M. Weinert and D. Braun, “Optically driven fluid flow along arbitrary microscale patterns using thermoviscous expansion,” J. Appl. Phys. **104**(10), 104701 (2008). [CrossRef]

**6. **H. Hidai, T. Yamazaki, S. Itoh, K. Hiromatsu, and H. Tokura, “Metal particle manipulation by laser irradiation in borosilicate glass,” Opt. Express **18**(19), 20313–20320 (2010). [CrossRef] [PubMed]

**7. **S. Maeda, “Kihou no seisei to undou” Chem. Eng. J. **31**, 438–443 (1967) (in Japanese).

**8. **C. W. Parmelee, K. C. Lyon, and C. G. Harman, “The surface tensions of molten glass,” Univ. Illinois Eng. Exp. Station Bull. **36**, 5–50 (1939).

**9. **. Liu, J. Shang, J. Tang, and Q.-A. Huang, “Micromachining of Pyrex 7740 glass by silicon molding and vacuum anodic bonding,” J. Microelectromech. Syst. **20**(4), 909–915 (2011). [CrossRef]

**10. **H. Hidai, M. Yoshioka, K. Hiromatsu, and H. Tokura, “Glass modification by continuous-wave laser backside irradiation (CW-LBI),” Appl. Phys., A Mater. Sci. Process. **96**(4), 869–872 (2009). [CrossRef]

**11. **H. Hidai, M. Yoshioka, K. Hiromatsu, and H. Tokura, “Structural changes in silica glass by continuous-wave laser backside irradiation,” J. Am. Ceram. Soc. **93**, 1597–1601 (2010).