Loading [MathJax]/jax/element/mml/optable/SuppMathOperators.js

Liu J, Zheng M, Xiong ZJ, Li ZY. 3D dynamic motion of a dielectric micro-sphere within optical tweezers. Opto-Electron Adv 4, 200015 (2021).. doi: 10.29026/oea.2021.200015
Citation: Liu J, Zheng M, Xiong ZJ, Li ZY. 3D dynamic motion of a dielectric micro-sphere within optical tweezers. Opto-Electron Adv 4, 200015 (2021).. doi: 10.29026/oea.2021.200015

Original Article Open Access

3D dynamic motion of a dielectric micro-sphere within optical tweezers

More Information
  • Known as laser trapping, optical tweezers, with nanometer accuracy and pico-newton precision, plays a pivotal role in single bio-molecule measurements and controllable motions of micro-machines. In order to advance the flourishing applications for those achievements, it is necessary to make clear the three-dimensional dynamic process of micro-particles stepping into an optical field. In this paper, we utilize the ray optics method to calculate the optical force and optical torque of a micro-sphere in optical tweezers. With the influence of viscosity force and torque taken into account, we numerically solve and analyze the dynamic process of a dielectric micro-sphere in optical tweezers on the basis of Newton mechanical equations under various conditions of initial positions and velocity vectors of the particle. The particle trajectory over time can demonstrate whether the particle can be successfully trapped into the optical tweezers center and reveal the subtle details of this trapping process. Even in a simple pair of optical tweezers, the dielectric micro-sphere exhibits abundant phases of mechanical motions including acceleration, deceleration, and turning. These studies will be of great help to understand the particle-laser trap interaction in various situations and promote exciting possibilities for exploring novel ways to control the mechanical dynamics of microscale particles.
  • Micro-particles trapped in optical tweezers1, 2 can be utilized as carriers or handles that are biochemically linked to biological molecules. With nanometer accuracy of positioning and pico-newton precision of force measurement, optical tweezers offer an ideal platform to exert and measure the forces and torques on single bio-molecules or micro-particles3-5. In 2018, A. Ashkin was awarded the Nobel Physics Prize for his contributions to invent the optical tweezers technology and promote its applications to biological systems. With the rapid development of optical tweezers technology, simple manipulation of particles is no longer sufficient to support either scientific research or industry demands. For example, the first necessary step for statistically analyzing the trajectory of a trapped particle is to track the particle position accurately. This is often achieved by acquiring a video6, 7, stereomicroscopy8, V-shaped micro-mirrors9 of the particle, or other visualization techniques including position sensitive detection10 and quadrant photodiode detection11, which monitors via detecting the backscattering signal or transmitting signal. In fact, these experimental techniques are limited by sampling rate and in 3D imaging applications. Up to now, it is still very hard to track the 3D trajectory of a micro-particle with high precision. Currently, numerical researches12, 13 are conducted to investigate the optical force exerted on trapped particles. However, the detailed dynamics of a particle stepping into a laser trap14, 15 is still an area rarely studied.

    This situation inspires us to think of an idea to drive a particle into a steady state by a focused laser beam illuminating it. As it is very difficult or even impossible to determine trapping probabilities through real-world experiments due to the large parameter space (including, for example, the initial states of particles position, velocity, etc.), we have developed a computational framework in which trapping probabilities can be determined easily through numerical simulations. Simulation involves calculating the optical forces and torques, non-optical forces and torques, and combining these components to calculate the overall dynamics of the particle.

    A significant concern is how to calculate the optical force and torque of mesoscopic particles. There are two main types of approaches: namely ray optics method and electromagnetic scattering method. In terms of the computational resources and time to solve the dynamic process, the ray optics method has advantages over electromagnetic approaches such as finite-difference time-domain (FDTD) method16 or T-matrix method17, 18. Taking the widely used method of FDTD as an example, the simulation area is needed to be 5 μm × 5 μm × 5 μm for a 3 μm diameter micro-particle. When the micro-particle is far away from the laser trap, the simulation area is even larger. Furthermore, a complete dynamic process involves thousands or even millions of trajectory points, and it is needed to calculate all mechanical quantities (forces and torques) at each point in order to determine the next point of the trajectory by solving electromagnetic fields in a large 3D space. With FDTD or T-matrix method, the time consumption at each point is at the scale of minutes or even hours. Therefore, the solution to a single set of dynamic process costs days or even months. This is a computational burden too heavy which exceeds the capacity of common computers. As for the semi-analytical ray optics method, the time cost is at the scale of millisecond to implement the optical force and torque calculation of each trajectory point, which makes this method a good compromise between clarity and exactness in many situations19, 20. Therefore, we employ the ray optics method to obtain the optical force and torque of the micro-particles in the optical tweezers.

    In this paper, we utilize the ray optics theory to calculate optical forces and optical torques. Then, with the influence of viscosity forces and torques taken into account, we numerically analyze the 3D dynamic process of dielectric micro-spheres in optical tweezers on the basis of Newton mechanical equations. In addition, the supplementary videos show the 3D dynamic process of the micro-sphere particles with different initial states of particle positions and velocities. These studies can be valuable to understand experiments that have been performed, to predict the results of potential experiments, or to explore the mechanical dynamics of trapped particles that are not accessible experimentally. Furthermore, these studies can open up opportunities to explore possible applications based on controlling the dynamics of a particle in the optical tweezers.

    Investigating the dynamic process of a micro-sphere can be divided into three steps. The first step is to calculate the optical force and torque, the second step is to calculate the non-optical force and torque, and the third and final step is to numerically solve the Newton mechanical equations. This section describes the theory and formulas for solving the dynamic process of the micro-sphere in optical tweezers systematically.

    The schematic of the optical tweezers considered in this work and the coordinate system used in our calculation is shown in Fig. 1(a). The center of the optical tweezers, point “O”, is set as the coordinate origin. The collimated Gaussian beam at a wavelength λ=1064nm propagates along the z-axis and is focused by a high numerical aperture (NA=1.4) objective. We implement the procedures and methods developed in ref.21 to calculate the optical force (Foptical=[Fx,Fy,Fz]) and the optical torque (τoptical=[τx,τy,τz]) of a dielectric micro-sphere in optical tweezers when the micro-sphere with radius rs moves around the focus point along the x-, y-, and z-axes, respectively. The initial parameters are as follows. The surrounding medium is water with a refractive index n1=1.33, the dielectric micro-sphere has a refractive index of n2=1.6 and radius rs=1μm, and the laser power is P=10mW. In Figs. 1(b)1(d), we display the calculated optical force and torque as a function of the position along the x− and z− axes of the micro-sphere around the identified equilibrium position. The results analysis about the initial parameters of the micro-sphere in the y-axis can be deduced from that in the x-axis. Due to the symmetry of the x− and y− axis directions, the curve FopticalSy is the same to the curve FopticalSx, and the curve τopticalSy is mirror symmetrical to the curve τopticalSx. In addition, to characterize the optical trapping probability, we need to measure the trap stiffness κ=[κx,κy,κz]. If the particle is not far away from the trap center (so-called Hookean region, |ΔS/rs|), the optical force exerted on the particle is directly proportional to the displacement {\rm{\Delta}} {\boldsymbol{S}} = [{\rm{\Delta}} {S_x}, {\rm{\Delta}} {S_y},{\rm{\Delta}} {S_z}] from the equilibrium position {\boldsymbol{S}} = [{S_x},{S_y},{S_z}] , i.e. {{\boldsymbol{F}}_{{\rm{optical}}}} = {\kappa _x}{\rm{\Delta}} {S_x}\hat x +{\kappa _y}{\rm{\Delta}} {S_y}\hat y + {\kappa _z}{\rm{\Delta}} {S_z}\hat z. Through linear fitting of the data in Hookean region, the trap stiffness coefficients in x− and z− axes are calculated and they are {\kappa _x} = {\text{-}} 12.47\;{\rm{pN/\mu m}} and {\kappa _z} = {\text{-}} 12.93\;{\rm{pN/\mu m}}, respectively. The yellow line in Figs. 1(b) and 1(c) characterizes the trap stiffness. In Fig. 1(d), the micro-sphere only experiences the torque in the y-axis ({\tau _y}) when it moves along the x-axis direction. And it can be inferred that there is no torque acting on the micro-sphere when it moves along the z-axis because of the structure symmetry, as shown in Fig. 1(e).

    Figure 1.  (a) Schematic of basic optical tweezers. A single laser beam is focused by a high numerical aperture objective lens to form a stable optical trap for microscale particle. (b) and (c) illustrate the calculated optical force acting on a micro-sphere when it moves along the x-axis and z-axis, respectively, in the optical tweezers at the laser power P = 10 mW. Here, the two most important features of optical tweezers are shown: the maximum reverse axial force that characterizes the strength of the trap, and spring constants (the slope of the fit line). (d) and (e) display the calculated optical torque acting on a micro-sphere when it moves along the x-axis and z-axis, respectively.

    In general, all the particles experience drag forces coupled with a rapid fluctuation force, which is the result of frequent and numerous collisions with surrounding liquid molecules and the physical origin of famous Brownian motion. The drag force acting on the particle is obtained by the expression

    {{\boldsymbol{F}}_{{\rm{drag}}}} = - \gamma \cdot {\boldsymbol{v}}\;, (1)

    where \gamma is the viscous drag coefficient and {\boldsymbol{v}} is the particle velocity relative to the surrounding liquid. For a spherical bead, \gamma is expressed by the Stokes equation,

    \gamma = 6{{\rm{\pi}}} \eta {r_{\rm{s}}}\;, (2)

    where \eta is the fluid viscosity, 0.89 \times {10^{ - 3}}\;{\rm{Pa}} \cdot {\rm{s}} for water at room temperature and {r_s} represents the radius of the micro-sphere. Meanwhile, the drag torque can be calculated as follows,

    {{\boldsymbol{\tau}} _{{\rm{drag}}}} = 8{\rm{\pi}} r_{\rm{s}}^3\eta {\boldsymbol{\omega}} \;, (3)

    here, {\boldsymbol{\omega}} is the angular velocity of the spherical particle.

    Moreover, the gravity {\boldsymbol{G}} and buoyancy {{\boldsymbol{F}}_{{\rm{buoy}}}} acting on the micro-sphere are written using the density of the micro-sphere and the density of water, {\rho _{\rm{s}}} = 1.03 \times{10^3}\;{\rm{kg/}}{{\rm{m}}^{\rm{3}}} and {\rho _{\rm{w}}} = 1.0 \times {10^3}\;{\rm{kg/}}{{\rm{m}}^{\rm{3}}} respectively, the radius of the micro-sphere {r_{\rm{s}}}, and the gravitational acceleration g,

    {\boldsymbol{G}} = - \frac{4}{3}{\rm{\pi}} r_{\rm{s}}^3{\rho _{\rm{s}}}g\;, (4)
    {{\boldsymbol{F}}_{{\rm{buoy}}}} = \frac{4}{3}{\rm{\pi}} r_{\rm{s}}^3{\rho _{\rm{w}}}g\;. (5)

    It is noted that the direction of the forces {\boldsymbol{G}} and {{\boldsymbol{F}}_{{\rm{buoy}}}} is along the z-axis. Therefore, their resultant forces pass through the center of mass, which provides no torque to the micro-sphere.

    In order to implement numerical analysis over the dynamic process, the Newton mechanical equations are exploited to accurately model the dynamic parameters of the trapped micro-sphere in the optical tweezers. The equation of the micro-sphere translational motion22-24 can be written as

    m\ddot {\boldsymbol{r}} = {{\boldsymbol{F}}_{{\rm{optical}}}}({\boldsymbol{r}},P) + {{\boldsymbol{F}}_{{\rm{drag}}}}({\boldsymbol{v}}) + {\boldsymbol{G}} + {{\boldsymbol{F}}_{{\rm{buoy}}}}\;, (6)

    where {\boldsymbol{r}} is the micro-sphere position with respect to the laser center; m is the mass of the micro-sphere.

    Similarly, the equation of the micro-sphere rotational motion21 can be expressed by

    I\ddot {\bf\textit{ϕ}} = {{\boldsymbol{\tau}} _{{\rm{optical}}}}({\bf\textit{ϕ}} ,P) + {{\boldsymbol{\tau}} _{{\rm{drag}}}}({\boldsymbol{\omega}} )\;, (7)

    here, I = {2 / 5}\;mr_{\rm{s}}^2 represents the moment of inertia of the micro-sphere, and {\bf\textit{ϕ}} is the initial angular of the micro-sphere.

    As Eqs. (6) and (7) cannot be solved analytically, the Runge-Kutta methods, a family of implicit and explicit iterative methods, are used in temporal discretization for the approximate solutions of ordinary differential equations (ODEs). One of the above methods, the fourth-order Runge-Kutta method (RK4 method) is utilized. This method has successfully been employed for solving the ODEs25-27. The second-order differential equation Eqs. (6) and (7) each can be converted to two coupled first-order equations respectively,

    {\dot {\boldsymbol{v}}} = \frac{{{{\boldsymbol{F}}_{{\rm{optical}}}}({\boldsymbol{r}},P) + {{\boldsymbol{F}}_{{\rm{drag}}}}({\boldsymbol{v}}) + {\boldsymbol{G}} + {{\boldsymbol{F}}_{{\rm{buoy}}}}}}{m} = {\boldsymbol{a}}({\boldsymbol{r}},{\boldsymbol{v}})\;, (8)
    \dot {\boldsymbol{r}} = {\boldsymbol{v}}\;. (9)

    And

    {\dot{\boldsymbol{\omega}}} = \frac{{{{\boldsymbol{\tau}} _{{\rm{optical}}}}({\bf\textit{ϕ}} ,P) + {{\boldsymbol{\tau}} _{{\rm{drag}}}}({\boldsymbol{\omega}} )}}{I} = {\boldsymbol{\alpha}} ({\bf\textit{ϕ}} ,{\boldsymbol{\omega}} )\;, (10)
    \dot {\bf\textit{ϕ}} = {\boldsymbol{\omega}} \;. (11)

    Since the time step has a direct effect on the number of iterations required for running a simulation, we have to choose an appropriate time step size. The characteristic time scale of our entire model is given by the relaxation time {m / \gamma } \sim 27\;{\rm{ns}}. In our computations, the numerical integration time step \delta t is set to 10 {\rm{ns}}, such that \delta t < {m / \gamma }. Choosing such a small \delta t provides an opportunity to observe interesting non-equilibrium behavior in the simulation.

    The recursive algorithm for the classical RK4 method can be written as follows. Taking a particular position step {{\boldsymbol{r}}_i} for an example, the mechanical dynamics evolution can be written as

    {{\boldsymbol{f}}_{a1}} = {\boldsymbol{a}}({{\boldsymbol{r}}_i},{{\boldsymbol{v}}_i})\;, (12)
    {{\boldsymbol{f}}_{v1}} = {{\boldsymbol{v}}_i}\;, (13)
    {{\boldsymbol{f}}_{a2}} = {\boldsymbol{a}}({{\boldsymbol{r}}_i} + \delta t/2 \cdot {{\boldsymbol{f}}_{v1}},{{\boldsymbol{v}}_i} + \delta t/2 \cdot {{\boldsymbol{f}}_{a1}})\;, (14)
    {{\boldsymbol{f}}_{v2}} = {{\boldsymbol{v}}_i} + \delta t/2 \cdot {{\boldsymbol{f}}_{a1}}\;, (15)
    {{\boldsymbol{f}}_{a3}} = {\boldsymbol{a}}({{\boldsymbol{r}}_i} + \delta t/2 \cdot {{\boldsymbol{f}}_{v2}},{{\boldsymbol{v}}_i} + \delta t/2 \cdot {{\boldsymbol{f}}_{a2}})\;, (16)
    {{\boldsymbol{f}}_{v3}} = {{\boldsymbol{v}}_i} + \delta t/2 \cdot {{\boldsymbol{f}}_{a2}}\;, (17)
    {{\boldsymbol{f}}_{a4}} = {\boldsymbol{a}}({{\boldsymbol{r}}_i} + \delta t/2 \cdot {{\boldsymbol{f}}_{v3}},{{\boldsymbol{v}}_i} + \delta t/2 \cdot {{\boldsymbol{f}}_{a3}})\;, (18)
    {{\boldsymbol{f}}_{v4}} = {{\boldsymbol{v}}_i} + \delta t/2 \cdot {{\boldsymbol{f}}_{a3}}\;. (19)

    The position and velocity of the micro-sphere in the next integration time interval are solved as

    {{\boldsymbol{r}}_i}_{ + 1} = {{\boldsymbol{r}}_i} + (\delta t/6) \cdot ({{\boldsymbol{f}}_{v1}} + 2{{\boldsymbol{f}}_{v2}} + 2{{\boldsymbol{f}}_{v3}} + {{\boldsymbol{f}}_{v4}})\;, (20)
    {{\boldsymbol{v}}_i}_{ + 1} = {{\boldsymbol{v}}_i} + (\delta t/6) \cdot ({{\boldsymbol{f}}_{a1}} + 2{{\boldsymbol{f}}_{2a}} + 2{{\boldsymbol{f}}_{a3}} + {{\boldsymbol{f}}_{a4}})\;. (21)

    Similarly, the mechanical evolution of the remaining physical values {{\boldsymbol{\omega}} _{i + 1}}, {{\boldsymbol{ϕ}} _{i + 1}} follows the same procedures described above.

    With the above iterative methods, we can obtain the dynamic parameters {\boldsymbol{r}}(t), {\boldsymbol{v}}(t), {\boldsymbol{\omega}} (t), and {\boldsymbol{ϕ}} (t) of the micro-sphere in the optical tweezers with high precision and retrieve every details of the micro-sphere mechanical motion in 3D space. These results can be interpreted in a direct correspondence to the time-space dynamic process. Note in particular that this micro-sphere simulation framework is applicable to either a continuous or a discrete optical force.

    On the basis of the above simulation framework, we set a series of initial parameters for the micro-sphere in the optical tweezers, and go further to theoretically draw a picture of its dynamics of mechanical motion. At first, we choose a special initial position of the micro-sphere close to the laser center. This special initiation indicates a high trapping probability, which can help find out whether the simulation framework agrees with experimental results quantitatively or not.

    The initial parameters are as follows: the initial position {{\boldsymbol{r}}_0} = [0,0,1]\;{\rm{\mu m}}, the initial velocity {{\boldsymbol{v}}_0} = [0,0,0]\;{\rm{\mu m/s}} , the initial angular displacement {{\bf\textit{ϕ}} _0} = [0,0,0]\;{\rm{rad}}, and the initial angular speed {{\boldsymbol{\omega}} _0} = [0,0,0]\;{\rm{rad/s}}. In addition, the gravity and buoyancy of the micro-sphere are {\boldsymbol{G}} = - 4.27 \times {10^{ - 2}}\;{\rm{pN}} and {{\boldsymbol{F}}_{{\rm{buoy}}}} = 4.11 \times {10^{ - 2}}\;{\rm{pN}}, respectively. The simulated numerical results of trajectory, velocity, angular velocity, and optical force of the micro-sphere in the optical tweezers are plotted in Fig. 2 for the laser power of 10 mW. In Fig. 2(a), which illustrates the temporal evolution of the micro-sphere position, it can be seen that the micro-sphere goes straight to the equilibrium position, and finally is tightly trapped in optical tweezers. On the other hand, according to Fig. 2(b), which displays the temporal evolution of micro-sphere velocity, in the beginning, the micro-sphere is in the stationary state. Then it accelerates to a high velocity {v_{\max\_z}} = - 671.79\;{\rm{\mu m/s}} in an extremely short time (below one microsecond), followed by a modest deceleration until a stable capture is reached in a time scale of several milliseconds. Due to the radial symmetrical structure of the micro-sphere and optical tweezers, the micro-sphere only experiences the optical force along the z-axis, {F_z} [shown in Fig. 2(d)]. Moreover, due to the radial symmetrical force exerted on the micro-sphere, no torque exists to supports its rotation [displayed in Fig. 2(c)]. The data in Figs. 2(c) and 2(d) are consistent with the physical analysis and naive expectation. The stable position is [0,0,0.085]. It shows that the trapping potential well center for the micro-sphere particle is in the optical axis, i.e., the z-axis, and slightly shifting forward (0.085 {\rm{\mu m}}) away from the laser beam focus center without the particle, due to the refraction effect of the particle against the laser beam. Besides, the micro-sphere spends about 0.008 second on reaching the stable position, which agrees well with the results reported in ref.28.

    Figure 2.  The dynamic process of the micro-sphere in an optical tweezers. Plots of the (a) position, (b) velocity, (c) angular velocity, and (d) optical force versus time under the initial condition of r0=[0,0,1] μm, v0=[0,0,0] μm/s, ϕ0=[0,0,0] rad, and ω0=[0,0,0] rad/s.

    Having confirmed the accuracy and efficiency of the numerical methods, we proceed to estimate the trapping probabilities by systematically performing simulation at random initial positions and velocities in the parameter space. This greatly helps on judging whether the micro-sphere can be trapped stably and efficiently in the optical tweezers. In the following, we first consider the dynamics of microsphere with different initial positions but with fixed initial zero velocity. We compare the motion trajectories of the micro-sphere when it is set at different initial positions in the optical tweezers. Apart from the initial positions, the rest parameters are as follows, the initial velocity {{\boldsymbol{v}}_0} = [0,0,0]\;{\rm{\mu m/s}}, the angular displacement {{\bf\textit{ϕ}} _0} = [0,0,0]\;{\rm{rad}}, and the angular speed {{\boldsymbol{\omega}} _0} = [0,0,0]\;{\rm{rad/s}}.

    In the first case, the initial position of the micro-sphere is set to be located on the horizontal plane crossing the trapping center and moving along the x-axis, so that the coordinate on the y- and z- axes is given by {r_{0y}} = 0\;{\rm{\mu m}} and {r_{0z}} = 0.085\;{\rm{\mu m}}. Figure 3(a) shows the temporal evolution of the x-directional trajectories of the micro-sphere at different initial positions in the x-axis with {r_{0x}} = 0.5, 1.0, 1.5, 1.8, 1.9, and 2.0 {\rm{\mu m}}, respectively. When the initial position is within the Hookean region (\left| {{r_{0x}}} \right| \leqslant 1.0\;{\rm{\mu m}}), the micro-sphere experiences a strong restoring force, as shown in Fig. 1(b). It rapidly falls into the laser center, then finally locates at the equilibrium position. Interestingly, even though the initial position is out of the Hookean region (1.5\;{\rm{\mu m}} \leqslant \left| {{r_{0x}}} \right| \leqslant 1.9\;{\rm{\mu m}}), the micro-sphere still can be trapped in the optical tweezers. This phenomenon is due to the fact that the micro-sphere experiences the synergetic effect between the x-axis optical force {F_x} and the z-axis optical force {F_z}. At this time, the force {F_z} plays a very important role in dominating the final state of the micro-sphere motion. For clear illustration, we plot the detailed pictures of the trajectory, the optical force and the viscous drag force, and the resultant force {{\boldsymbol{F}}_{{\rm{sum}}}}({{\boldsymbol{F}}_{{\rm{sum}}}} = {{\boldsymbol{F}}_{{\rm{optical}}}} + {{\boldsymbol{F}}_{{\rm{drag}}}} + {\boldsymbol{G}} + {{\boldsymbol{F}}_{{\rm{buoy}}}}, which are the same in all context) of the micro-sphere placed at {r_{0x}} = 1.9 {\rm{\mu m}}, and display them in Figs. 3(b), 3(c) and 3(d), respectively. In Fig. 3(b), the micro-sphere moves towards the laser center in the x-axis direction, while away from that in the z-axis direction in the early stage (T \leqslant 0.015\;{\rm{s}}). In the intermediate stage (0.015\;{\rm{s}} \leqslant T \leqslant 0.045\;{\rm{s}}), the micro-sphere moves slowly in the z-axis direction, and keeps on moving close to the laser center in the x-axis. In the last stage, the micro-sphere overcomes viscous resistance and slowly reaches the Hookean region ({\boldsymbol{r}} = [0.68,0,0.73]\;{\rm{\mu m}}) at the time T = 0.046\;{\rm{s}}, then goes straight to the equilibrium position with a great speed. In Fig. 3(c), the plot of the optical force {{\boldsymbol{F}}_{{\rm{optical}}}} = [{F_x},{F_y},{F_z}] and the drag force {{\boldsymbol{F}}_{{\rm{drag}}}} = [{F_{{\rm{drag}}\_x}},{F_{{\rm{drag}}\_y}},{F_{{\rm{drag}}\_z}}] exerted on the micro-sphere also helps us to understand its dynamic process. In the early stage, the negative force {F_x} pulls the micro-sphere into the trap, while the positive force {F_z} pushes it away. In the intermediate stage, the value of the force {F_z} is low so that the micro-sphere moves slowly in the z-axis. More importantly, the direction of the force {F_z} has changed from the positive z direction to the negative z direction, which indicates the capture of the micro-sphere. In the whole process, the drag force {{\boldsymbol{F}}_{{\rm{drag}}}} opposes its motion, and varies with the velocity of the micro-sphere. On the other hand, in Fig. 3(d), the resultant force {{\boldsymbol{F}}_{{\rm{sum}}}} are helpful in understanding the motion of the micro-sphere. For example, the resultant force is so small in the intermediate stage compared with other forces, and the micro-sphere moves very slowly correspondingly. When {r_{0x}}{\rm{ = }}2.0\;{\rm{\mu m}}, the optical force, drag force, and the resultant force over time are calculated and displayed in Figs. 3(e) and 3(f), respectively. The direction of the force {F_z} always points to the positive z axis, so that it pushes the micro-sphere away from the laser center. Meanwhile, the force {F_x} is too small to pull the micro-sphere back to the trap center. As the micro-sphere moves away from the laser center, the resultant force approaches zero. Therefore, this initial position will not lead to successful trap of the micro-sphere.

    Figure 3.  Mechanical motion dynamics of the micro-sphere at different initial positions in the horizontal trapping plane. (a) The x-directional trajectories of the micro-sphere at different initial positions in the x-axis in the optical tweezers with r0x=0.5, 1.0, 1.5, 1.8, 1.9, and 2.0 μm. (b)−(d) display the temporal evolution of the x- and y- and z- directional trajectory, optical force, drag force, and resultant force, respectively, when the micro-sphere is set at the initial position r0x=1.9 μm. (e) and (f) display the temporal evolution of the optical force and the drag force, and the resultant force, respectively, when r0x=2.0 μm.

    Figure 4 shows the corresponding 3D trajectories of the micro-sphere. The color bar represents the time. This picture can illustrate the temporal evolution of the micro-sphere mechanical motion in 3D space. Moreover, these 3D dynamic processes of the micro-sphere are better visualized in Supporting Information (SI) Video S1. In Fig. 4 and Video S1, it is obvious that the micro-sphere does not follow the shortest straight path to move to the stable position. Besides, the complexity of its trajectory increases as the initial position is located farther away from the laser center. When the micro-sphere is set at {r_{0x}} \le 1.5\;{\rm{\mu m}}, it moves fast along the corresponding trajectory, and takes less than 0.015 second to reach the stable position. While at {r_{0x}} = 1.8\;{\rm{\mu m}} or {\rm{1}}{\rm{.9}}\;{\rm{\mu m}}, its trajectory has been described in Fig. 3(b). As observed in Video S1, the micro-sphere exhibits abundant types of mechanical motions, including acceleration, deceleration, and turning. However, these short dynamic processes are easily overlooked in experiments. With {r_{0x}} = 2.0\;{\rm{\mu m}}, in a period of 0.06 second, the micro-sphere also undergoes many motion phases as in other initial positions, but without a successful capture in the end. If the micro-sphere is utilized as a handle linked to single bio-molecule, these specific motion details can provide key information on the dynamic process of the attached single bio-molecule.

    Figure 4.  The 3D trajectories of the micro-sphere in the optical tweezers at different initial positions in the x-axis with r0x=0.5, 1.0, 1.5, 1.8, 1.9, and 2 μm.

    In the second case, the initial position of the micro-sphere is set to be located at the z-axis with {r_{0x}} = 0\;{\rm{\mu m}}, {r_{0y}} = 0\;{\rm{\mu m}}. Fig. 5(a) shows the temporal evolution of the z-directional trajectories of the micro-sphere at different initial positions in the z-axis with {r_{0z}} = −2.0, −1.5, −1.0, −0.5, 1.0, 1.4, 1.6, 1.7, 1.75, and 1.9 {\rm{\mu m}}, respectively. The corresponding dynamic motions are shown in Supplementary information Video S2. Due to the symmetry of micro-sphere and optical tweezers, the optical force acting on the micro-sphere always has {F_x} = 0 and {F_y} = 0. Therefore, the micro-sphere experiences the optical force along the z-axis only and thus also moves along the z-axis only. In this situation, the forces {\boldsymbol{G}} and {{\boldsymbol{F}}_{{\rm{buoy}}}} also apply a little bit of work. When - 2\;{\rm{\mu m}} \leqslant {r_{0z}} \leqslant 0\;{\rm{\mu m}}, the required equilibrium time, which is the time for the trapped micro-sphere to reach its equilibrium position, is positively correlated with the distance from the trap center. When 0\;{\rm{\mu m}} \leqslant {r_{0z}} \leqslant 1.0\;{\rm{\mu m}}, the micro-sphere experiences a strong propulsion force, rapidly falls into the laser center, and eventually locates at the equilibrium position. When 1.4\;{\rm{\mu m}} \leqslant {r_{0z}} \leqslant 1.75\;{\rm{\mu m}}, the micro-sphere is captured in the end but the situation differs from above. As a reference, the sketches of the optical force and drag force, and the resultant force varying over time are shown in Figs. 5(b) and 5(c) when {r_{0z}} = 1.7\;{\rm{\mu m}}. In the initial stage, the micro-sphere moves slowly, which is primarily due to the small resultant force. While entering into the Hookean region, it rapidly falls into the trap because of the huge restoring optical force. Eventually, it will stop due to the drag force. When {r_{0z}} \geqslant 1.9\;{\rm{\mu m}}, the direction of the optical force {F_z} points to the positive z direction, as shown in Fig. 1(c), and it pushes the micro-sphere away from the laser center. Clearly, this initial position will not lead to successful trap of the micro-sphere.

    Figure 5.  Mechanical motion dynamics of the micro-sphere at different initial positions along the optical axis. (a) The z-directional trajectories of the micro-sphere at different initial positions in the z-axis in the optical tweezers with r0x=−2.0, −1.5, −1.0, −0.5, 1.0, 1.4, 1.6, 1.7, 1.75 and 1.9 μm. (b) and (c) The sketches of the optical force and the drag force, and the resultant force over time, respectively, when the micro-sphere is set at the initial position r0x=1.7 μm.

    In this part, we consider, calculate, and display the motion trajectories of micro-sphere when it is set at different initial velocities. Similar to the past sections, we have listed the initial parameters, the initial position {{\boldsymbol{r}}_0} = [0,0,0]\;{\rm{\mu m}}, the angular displacement {{\bf\textit{ϕ}} _0} = [0,0,0] rad, and the angular speed {{\boldsymbol{\omega}} _0} = [0,0,0]\;{\rm{rad/s}}. In the first case, the micro-sphere is set to be located at the laser trap center and has initial velocity components of {v_{0y}} = 0\;{\rm{\mu m/s}} and {v_{0z}} = 0\;{\rm{\mu m/s}}.

    The corresponding 3D trajectories and dynamic motions of the micro-sphere at different x-directional initial velocities with {v_{0x}} = 100, 1 \times {10^4}, 1 \times {10^6}, 5 \times {10^6}, 7.5 \times {10^6}, 7.75 \times {10^6}, and 8 \times {10^6} μm are shown in Fig. 6 and Supplementary information Video S3. In the initial stage, the micro-sphere with an initial velocity will continue to move and reach a new position because of the inertia. According to the equation {{\boldsymbol{F}}_{{\rm{drag}}}} = \gamma \cdot {\boldsymbol{v}}, the higher speed, the larger viscous force. Thus, the velocities decrease rapidly. Additionally, it moves under a variety of forces, including the optical force, the drag force, the gravity, and the buoyance. In the following processes, its motion is similar to the situations with the different x-axis initial positions in Fig. 4. Moreover, the x-directional initial velocity is essential to the trapping probabilities and the capture time. When {v_{0x}} \leqslant 5 \times {10^6}\;{\rm{\mu m/s}}, it takes less than 0.01 second for the micro-sphere to reach the steady state. While {v_{0x}}=7.75 \times {10^6}\;{\rm{\mu m/s}}, it takes about 0.06 second. When {v_{0x}} \geqslant 8 \times {10^6}\;{\rm{\mu m/s}}, the micro-sphere is hard to be captured. The experience gained via this systematical study is particularly useful in determining the initial velocity of a particle in the liquids required for the particle to be successfully trapped.

    Figure 6.  The 3D trajectories of the micro-sphere at different initial velocities in the x-axis in the optical tweezers with v0x=100, 1×104, 1×106, 5×106, 7.5×106, 7.75×106, and 8×106 μm/s.

    In this part, we go further to investigate the motion trajectories of the micro-sphere in more general initial conditions of positions and velocities. The initial parameters are set as follows: the angular displacement {{\bf\textit{ϕ}} _0} = [0,0,0]\;{\rm{rad}}, the angular speed {{\boldsymbol{\omega}} _0} = [0,0,0]\;{\rm{rad/s}}, the position is located at x-axis ({{\boldsymbol{r}}_0} = [1,0,0]\;{\rm{\mu m}}), while the velocity {{\boldsymbol{v}}_0} is parallel to the z-axis.

    Figure 7(a) shows the 3D trajectories of the micro-sphere at different z-directional initial velocities with {v_{0z}} = 1 \times {10^6}, 5 \times {10^6}, 6 \times {10^6}, and 6.5 \times {10^6} {\rm{\mu m/s}}, respectively. In addition, the initial position is marked with a pentagon. The corresponding dynamic motions can be observed in Supplementary information Video S5. Although the optical force exerted on the micro-sphere at the initial position is {{\boldsymbol{F}}_{{\rm{optical}}}} = [ - 11.9,0,5.3]\;{\rm{pN}}, the micro-sphere still moves upward with the initial z-directional velocities and reaches a new position because of the inertia. When {v_{0z}} = 1 \times {10^6}\;{\rm{\mu m/s}}, the micro-sphere just moves up a little bit. After that, it reaches the stable position rapidly, in less than 0.01 second. When 5 \times {10^6}\;{\rm{\mu m/s}} \leqslant {v_{0z}} \leqslant 6 \times {10^6}\;{\rm{\mu m/s}}, the micro-sphere moves towards the positive x direction while moving downwards concurrently. Later, it makes a large-angle turning. Specifically, the micro-sphere changes its motion direction, shifts in the negative x direction while moving downwards, and approaches the equilibrium position. Taking {v_{0z}} = 6 \times {10^6}\;{\rm{\mu m/s}} as an example, we have calculated and plotted the optical force and drag force, and the resultant force of the micro-sphere in Figs. 7(b) and 7(c), respectively. As is well known, both the magnitude and direction of the optical force {F_x} vary with its location. The key to the turning is that the direction of the force {F_x} is reversed when T = 0.016\;{\rm{s}}. After that, the micro-sphere moves slowly along the trajectory for about 0.02 seconds, and then rushes into the stable position. When {v_{0z}} = 6.5 \times {10^6}\;{\rm{\mu m/s}}, although the micro-sphere moves towards the negative z direction in the early stage, it keeps heading in the positive x direction in the whole process. As a result, it is hard to successfully trap the micro-sphere in this situation.

    Figure 7.  (a) The 3D trajectories of a micro-sphere at different initial velocities in the z-axis with v0x=1×106, 5×106, 6×106, and 6.5×106 μm/s, when the initial position is r0x=[1,0,0] μm. The plots of (b) the optical force and drag force, and (c) the resultant force over time of the micro-sphere when v0=[6.5 × 106,0,0] μm/s, and r0x=[1,0,0] μm.

    In the above discussions, we just list several special cases to describe the micro-sphere trapping probabilities. Nevertheless, there will be various combinations of initial parameters to characterize the motions of the micro-sphere in the optical tweezers. In any of these cases, the theoretical and numerical method proposed in this paper can give a reasonable solution and calculation of the mechanical motion dynamic process, and also explain how and why the micro-sphere moves, whether or not, how and why it can be trapped. It should be noted that the initial angular displacement and angular speed have not been discussed or analyzed in the present work. It is mainly because the optical torque is very small as shown in Fig. 1(d)] while the drag torque is relatively very large ({{\boldsymbol{\tau}} _{{\rm{drag}}}} = 8{\rm{\pi}} r_s^3\eta {\boldsymbol{\omega}} ). Thus, the rotation state of the micro-sphere will not last long and is very hard to be observed in experiments. There is still a critical problem that in the microscopic world, due to the random bombardment of the water molecules, the Brownian motions of particles become increasingly significant with the decreasing size of the particles. In terms of the trapped micro-sphere, the Reynolds number of fluid is rather small, which is much less than 1, therefore the viscous forces rather than Brownian force dominate its dynamic process.

    Through simulating the dynamic process, the change of the relevant physical parameters of the micro-sphere [{\boldsymbol{r}}(t),\;{\boldsymbol{v}}(t),\;{\bf\textit{ϕ}} (t),\;{\boldsymbol{\omega}} (t),\;{{\boldsymbol{F}}_{{\rm{optical}}}}(t),\;{{\boldsymbol{F}}_{{\rm{drag}}}}(t),\;{{\boldsymbol{F}}_{{\rm{sum}}}}(t)] over time are very clear. This can help us understand, and even predict the results of experiments such as trapping probabilities, balance time, the maximum velocity, etc. Meanwhile, the above method can also be utilized in analyzing the dynamic process of a complicated-geometry particle in complicated-structure beams, which creates new possibilities for constructing mixed laser beams to generate various motions for microscale particles.

    In summary, we have proposed a method based on semi-analytical ray optics model to visualize the dynamic process of a dielectric micro-sphere in optical tweezers. This approach allows for a fast calculation of optical force and torque exerted upon micro-sphere by focused laser beam with reasonable accuracy and makes it possible to simulate the whole mechanical motion trajectory and process of the particle under Newton’s law. We discuss the trapping probabilities of the micro-sphere at the different initial positions or velocities along the x- and z- axes respectively. It has been found that, unexpectedly, the particle still has the possibility to be trapped by the optical tweezers even if its position is out of the Hookean region. In contrast, the micro-sphere with a large initial velocity will not necessarily be trapped even if its position is in the Hookean region.

    On the other hand, even in simple optical tweezers, the dielectric micro-sphere exhibits abundant phases of mechanical motions, including acceleration, deceleration, and turning. Solution and analysis of the 3D microscale particle dynamic process within optical trapping is essential for understanding various mechanisms for engineering the light-matter mechanical interactions, the mechanical motion behavior of molecules or micro-particles in liquid, as well as the biophysics and biochemistry of macromolecules and cells. These dynamics will facilitate the investigation of the single-molecular kinematics, the complex motions of diverse particles and so on. In addition, this method opens up a new way to promote optimal external control in various applications, including molecular winding, nanoparticles sorting, and in vivo manipulations. However, this method is still not mature enough. In the future, we shall perfect the theory of the micro-particle dynamics, and take more factors into consideration, such as the Brownian motion of particles, the change of the viscous drag coefficient (due to the action of the cover glass), the physical characteristics of the particles (shape, radius, refractive index), etc.

  • This work is supported by the National Natural Science Foundation of China (Grant No. 11974119 and No. 11804399), and the Guangdong Innovative and Entrepreneurial Research Team Program (Grant No. 2016ZT06C594), the Fundamental Research Funds for the Central Universities, South-Central University for Nationalities (Grant No. CZQ20018), and National Key R&D Program of China (No. 2018YFA 0306200).

  • The authors declare no competing financial interests.

  • [1] Ashkin A, Dziedzic JM, Bjorkholm JE, Chu S. Observation of a single-beam gradient force optical trap for dielectric particles. Opt Lett 11, 288–290 (1986). doi: 10.1364/OL.11.000288

    CrossRef Google Scholar

    [2] Liu J, Li ZY. Controlled mechanical motions of microparticles in optical tweezers. Micromachines (Basel) 9, 232 (2018). doi: 10.3390/mi9050232

    CrossRef Google Scholar

    [3] Guo HL, Li ZY. Optical tweezers technique and its applications. Sci China Phys, Mech Astron 56, 2351–2360 (2013). doi: 10.1007/s11433-013-5355-3

    CrossRef Google Scholar

    [4] Mehta AD, Rief M, Spudich JA, Smith DA, Simmons RM. Single-molecule biomechanics with optical methods. Science 283, 1689–1695 (1999). doi: 10.1126/science.283.5408.1689

    CrossRef Google Scholar

    [5] Liu W, Dong DS, Yang H, Gong QH, Shi KB. Robust and high-speed rotation control in optical tweezers by using polarization synthesis based on heterodyne interference. Opto-Electron Adv 3, 200022 (2020). doi: 10.29026/oea.2020.200022

    CrossRef Google Scholar

    [6] Helgadottir S, Argun A, Volpe G. Digital video microscopy enhanced by deep learning. Optica 6, 506–513 (2019). doi: 10.1364/OPTICA.6.000506

    CrossRef Google Scholar

    [7] Huhle A, Klaue D, Brutzer H, Daldrop P, Joo S et al. Camera-based three-dimensional real-time particle tracking at kHz rates and Ångström accuracy. Nat Commun 6, 5885 (2015). doi: 10.1038/ncomms6885

    CrossRef Google Scholar

    [8] Hay RF, Gibson GM, Lee MP, Padgett MJ, Phillips DB. Four-directional stereo-microscopy for 3D particle tracking with real-time error evaluation. Opt Express 22, 18662–18667 (2014). doi: 10.1364/OE.22.018662

    CrossRef Google Scholar

    [9] Hajjoul H, Mathon J, Viero Y, Bancaud A. Optimized micromirrors for three-dimensional single-particle tracking in living cells. Appl Phys Lett 98, 243701 (2011). doi: 10.1063/1.3599586

    CrossRef Google Scholar

    [10] Huang L, Guo HL, Li KL, Chen YH, Feng BH et al. Three dimensional force detection of gold nanoparticles using backscattered light detection. J Appl Phys 113, 113103 (2013). doi: 10.1063/1.4795272

    CrossRef Google Scholar

    [11] Finer JT, Simmons RM, Spudich JA. Single myosin molecule mechanics: piconewton forces and nanometre steps. Nature 368, 113–119 (1994). doi: 10.1038/368113a0

    CrossRef Google Scholar

    [12] Lehmuskero A, Johansson P, Rubinsztein-Dunlop H, Tong LM, Käll M. Laser trapping of colloidal metal nanoparticles. ACS Nano 9, 3453–3469 (2015). doi: 10.1021/acsnano.5b00286

    CrossRef Google Scholar

    [13] Rahimzadegan A, Fruhnert M, Alaee R, Fernandez-Corbaton I, Rockstuhl C. Optical force and torque on dipolar dual chiral particles. Phys Rev B 94, 125123 (2016). doi: 10.1103/PhysRevB.94.125123

    CrossRef Google Scholar

    [14] Melzer JE, McLeod E. Fundamental limits of optical tweezer nanoparticle manipulation speeds. ACS Nano 12, 2440–2447 (2018). doi: 10.1021/acsnano.7b07914

    CrossRef Google Scholar

    [15] Kim K, Yoon J, Park Y. Simultaneous 3D visualization and position tracking of optically trapped particles using optical diffraction tomography. Optica 2, 343–346 (2015). doi: 10.1364/OPTICA.2.000343

    CrossRef Google Scholar

    [16] Qin JQ, Wang XL, Ding J, Chen J, Fan YX et al. FDTD approach to optical forces of tightly focused vector beams on metal particles. Opt Express 17, 8407–8416 (2009). doi: 10.1364/OE.17.008407

    CrossRef Google Scholar

    [17] Nieminen TA, Loke VLY, Stilgoe AB, Knöner G, Brańczyk AM et al. Optical tweezers computational toolbox. J Opt A: Pure Appl Opt 9, S196–S203 (2007). doi: 10.1088/1464-4258/9/8/S12

    CrossRef Google Scholar

    [18] Stilgoe AB, Mallon MJ, Cao YY, Nieminen TA, Rubinsztein-Dunlop H. Optical tweezers toolbox: better, faster, cheaper; choose all three. Proc SPIE 8458, 84582E (2012). doi: 10.1117/12.929365

    CrossRef Google Scholar

    [19] Zong YW, Liu J, Liu R, Guo HL, Yang MC et al. An optically driven bistable janus rotor with patterned metal coatings. ACS Nano 9, 10844–10851 (2015). doi: 10.1021/acsnano.5b03565

    CrossRef Google Scholar

    [20] Potoček V, Barnett SM. Generalized ray optics and orbital angular momentum carrying beams. New J Phys 17, 103034 (2015). doi: 10.1088/1367-2630/17/10/103034

    CrossRef Google Scholar

    [21] Liu J, Zhang C, Zong YW, Guo HL, Li ZY. Ray-optics model for optical force and torque on a spherical metal-coated Janus microparticle. Photon Res 3, 265–274 (2015). doi: 10.1364/PRJ.3.000265

    CrossRef Google Scholar

    [22] Hajizadeh F, Shao L, Andrén D, Johansson P, Rubinsztein-Dunlop H et al. Brownian fluctuations of an optically rotated nanorod. Optica 4, 746–751 (2017). doi: 10.1364/OPTICA.4.000746

    CrossRef Google Scholar

    [23] Simon A, Libchaber A. Escape and synchronization of a Brownian particle. Phys Rev Lett 68, 3375–3378 (1992). doi: 10.1103/PhysRevLett.68.3375

    CrossRef Google Scholar

    [24] Bui AAM, Stilgoe AB, Lenton ICD, Gibson LJ, Kashchuk AV et al. Theory and practice of simulation of optical tweezers. J Quant Spectrosc Radiat Transf 195, 66–75 (2017). doi: 10.1016/j.jqsrt.2016.12.026

    CrossRef Google Scholar

    [25] Sanderse B, Koren B. Accuracy analysis of explicit Runge–Kutta methods applied to the incompressible Navier–Stokes equations. J Comput Phys 231, 3041–3063 (2012). doi: 10.1016/j.jcp.2011.11.028

    CrossRef Google Scholar

    [26] Zingg DW, Chisholm TT. Runge–Kutta methods for linear ordinary differential equations. Appl Numer Math 31, 227–238 (1999). doi: 10.1016/S0168-9274(98)00129-9

    CrossRef Google Scholar

    [27] Ramos H, Vigo-Aguiar J. A fourth-order Runge–Kutta method based on BDF-type Chebyshev approximations. J Comput Appl Math 204, 124–136 (2007). doi: 10.1016/j.cam.2006.04.033

    CrossRef Google Scholar

    [28] Simmons RM, Finer JT, Chu S, Spudich JA. Quantitative measurements of force and displacement using an optical trap. Biophys J 70, 1813–1822 (1996). doi: 10.1016/S0006-3495(96)79746-1

    CrossRef Google Scholar

  • Related articles

    Ayush Pandey, Zetian Mi
    Opto-Electronic Advances    2024, 7(3)   doi: 10.29026/oea.2024.240011
    Jose Javier Imas, Ignacio R. Matías, Ignacio Del Villar, Aritz Ozcáriz, Carlos Ruiz Zamarreño, Jacques Albert
    Opto-Electronic Advances    2023, 6(10)   doi: 10.29026/oea.2023.230048
    Xinyao Yuan, Quan Xu, Yuanhao Lang, Xiaohan Jiang, Yuehong Xu, Xieyu Chen, Jie Han, Xueqian Zhang, Jiaguang Han, Weili Zhang
    Opto-Electronic Advances    2023, 6(4)   doi: 10.29026/oea.2023.220133
    Ji Tong Wang, Pavel Tonkaev, Kirill Koshelev, Fangxing Lai, Sergey Kruk, Qinghai Song, Yuri Kivshar, Nicolae C. Panoiu
    Opto-Electronic Advances    2024, 7(5)   doi: 10.29026/oea.2024.230186
    Jie Li, Xueguang Lu, Hui Li, Chunyu Song, Qi Tan, Yu He, Jingyu Liu, Li Luo, Tingting Tang, Tingting Liu, Hang Xu, Shuyuan Xiao, Wanxia Huang, Yun Shen, Yan Zhang, Yating Zhang, Jianquan Yao
    Opto-Electronic Advances    2024, 7(10)   doi: 10.29026/oea.2024.240075
    Xiaohu Zhang, Qinmiao Chen, Dongliang Tang, Kaifeng Liu, Haimo Zhang, Lintong Shi, Mengyao He, Yongcai Guo, Shumin Xiao
    Opto-Electronic Advances    2024, 7(5)   doi: 10.29026/oea.2024.230126
    Qichang Ma, Guixin Li
    Opto-Electronic Advances    2024, 7(10)   doi: 10.29026/oea.2024.240166
    Daniele Pirone, Daniele Sirico, Lisa Miccio, Vittorio Bianco, Martina Mugnano, Danila del Giudice, Gianandrea Pasquinelli, Sabrina Valente, Silvia Lemma, Luisa Iommarini, Ivana Kurelac, Pasquale Memmolo, Pietro Ferraro
    Opto-Electronic Advances    2023, 6(1)   doi: 10.29026/oea.2023.220048
  • Video S1: 1Rx.avi
    Video S2: 2Rz.avi
    Video S3: 3Vx.avi
    Video S4: 4Vz.avi
    Video S5: 5RxVz.avi
  • 1.  Liu, J., Feng, Y., Zheng, M. et al. Ultra-thin rapeseed oil-assisted femtosecond laser etching on quartz glass microgrooves. iScience, 2025, 28(2): 111921.
    2.  Yin, G., Wang, Q., Lu, Q. et al. Lamellar beam with similar propagation and imaging characteristics to a Bessel beam. Optics Communications, 2025.
    3.  Zhang, R., Zhao, X., Li, J. et al. Programmable photoacoustic patterning of microparticles in air. Nature Communications, 2024, 15(1): 3250.
    4.  Wang, J., Chen, R., Zang, X. et al. Generation and switching of focused vector and vortex beams by a composite polarization hologram. Optics Express, 2024, 32(24): 43663-43678.
    5.  Bandaru, S., Arora, D., Ganesh, K.M. et al. Recent Advances in Research from Nanoparticle to Nano-Assembly: A Review. Nanomaterials, 2024, 14(17): 1387.
    6.  Shangguan, D., Liu, Y., Chen, L. et al. Modeling and experiment of femtosecond laser processing of micro-holes arrays in quartz. Journal of Applied Physics, 2024, 135(24): 243102.
    7.  Liu, Y., Gholam-Mirzaei, S., Khatri, D. et al. Field-resolved space–time characterization of few-cycle structured light pulses. Optica, 2024, 11(6): 846-851.
    8.  Liu, J., Long, L., Guo, H. et al. Observation of Moon-like Synchronous Revolution and Rotation of Janus Microparticles Trapped in an Annular Optical Trap. ACS Photonics, 2024.
    9.  Miao, W., Long, L., Zhang, F. et al. Characterization of optical forces and torques in polystyrene spheres steered by tightly focused orbital angular momentum beams. Proceedings of SPIE - The International Society for Optical Engineering, 2024.
    10.  Savelyev, D.A.. The Features of the Optical Traps Formation Using Silicon Ring Gratings with Variable Height. Photonics, 2023, 10(11): 1264.
    11.  Liu, H., Xu, J., He, H. et al. Optimization of micropore fabrication on the surface of ultrathick polyimide film based on picosecond UV laser. Journal of Laser Applications, 2023, 35(4): 042032.
    12.  Zhang, F., Camarero, P., Haro-González, P. et al. Optical trapping of optical nanoparticles: Fundamentals and applications. Opto-Electronic Science, 2023, 2(9): 230019.
    13.  Savelyev, D.A., Karpeev, S.V. Development of 3D Microstructures for the Formation of a Set of Optical Traps on the Optical Axis. Photonics, 2023, 10(2): 117.
    14.  Yan, C., Liang, C., Li, C. et al. Preparation of Transmission Quartz Glass Grating via Femtosecond Laser Etching | [飞 秒 激 光 刻 蚀 透 射 式 石 英 玻 璃 光 栅 工 艺 研 究]. Laser and Optoelectronics Progress, 2022, 59(7): 0714006.
    15.  Liu, J., Li, Z. Control of Surface Plasmon Resonance in Silver Nanocubes by CEP-Locked Laser Pulse. Photonics, 2022, 9(2): 53.
    16.  Lu, Y., Xu, Y., Ouyang, X. et al. Cylindrical vector beams reveal radiationless anapole condition in a resonant state. Opto-Electronic Advances, 2022, 5(4): 210014.

    Other cited types(1)

  • Created with Highcharts 5.0.7Amount of accessChart context menuAbstract Views, PDF Downloads StatisticsAbstract ViewsPDF Downloads2024-042024-052024-062024-072024-082024-092024-102024-112024-122025-012025-022025-0301020304050Highcharts.com
    Created with Highcharts 5.0.7Chart context menuAccess Class DistributionFULLTEXT: 87.1 %FULLTEXT: 87.1 %PDF: 12.9 %PDF: 12.9 %FULLTEXTPDFHighcharts.com
    Created with Highcharts 5.0.7Chart context menuAccess Area DistributionOthers: 50.1 %Others: 50.1 %United States: 27.5 %United States: 27.5 %Russia: 8.3 %Russia: 8.3 %France: 3.8 %France: 3.8 %Japan: 2.2 %Japan: 2.2 %India: 1.7 %India: 1.7 %Germany: 0.7 %Germany: 0.7 %Singapore: 0.6 %Singapore: 0.6 %United Kingdom: 0.6 %United Kingdom: 0.6 %Canada: 0.6 %Canada: 0.6 %Saudi Arabia: 0.5 %Saudi Arabia: 0.5 %Brazil: 0.4 %Brazil: 0.4 %Mexico: 0.3 %Mexico: 0.3 %Australia: 0.3 %Australia: 0.3 %Sweden: 0.2 %Sweden: 0.2 %Israel: 0.2 %Israel: 0.2 %Egypt: 0.2 %Egypt: 0.2 %Spain: 0.2 %Spain: 0.2 %Norway: 0.2 %Norway: 0.2 %Romania: 0.1 %Romania: 0.1 %Switzerland: 0.1 %Switzerland: 0.1 %Turkey: 0.1 %Turkey: 0.1 %Denmark: 0.1 %Denmark: 0.1 %Italy: 0.1 %Italy: 0.1 %Ukraine: 0.1 %Ukraine: 0.1 %Cyprus: 0.1 %Cyprus: 0.1 %Thailand: 0.1 %Thailand: 0.1 %Oman: 0.1 %Oman: 0.1 %Netherlands: 0.1 %Netherlands: 0.1 %Poland: 0.1 %Poland: 0.1 %Slovakia: 0.1 %Slovakia: 0.1 %Chile: 0.1 %Chile: 0.1 %Colombia: 0.1 %Colombia: 0.1 %Hungary: 0.1 %Hungary: 0.1 %New Zealand: 0.1 %New Zealand: 0.1 %Bangladesh: 0.1 %Bangladesh: 0.1 %Ethiopia: 0.1 %Ethiopia: 0.1 %Portugal: 0.0 %Portugal: 0.0 %Republic of Korea: 0.0 %Republic of Korea: 0.0 %OthersUnited StatesRussiaFranceJapanIndiaGermanySingaporeUnited KingdomCanadaSaudi ArabiaBrazilMexicoAustraliaSwedenIsraelEgyptSpainNorwayRomaniaSwitzerlandTurkeyDenmarkItalyUkraineCyprusThailandOmanNetherlandsPolandSlovakiaChileColombiaHungaryNew ZealandBangladeshEthiopiaPortugalRepublic of KoreaHighcharts.com
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Figures(7)

Article Metrics

Article views(6017) PDF downloads(887) Cited by(17)

Access History

Other Articles By Authors

Catalog

/

DownLoad:  Full-Size Img  PowerPoint