Next Article in Journal
The Cost of Improving the Precision of the Variational Quantum Eigensolver for Quantum Chemistry
Previous Article in Journal
An Innovative Formulation Based on Nanostructured Lipid Carriers for Imatinib Delivery: Pre-Formulation, Cellular Uptake and Cytotoxicity Studies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optical Force on a Metal Nanorod Exerted by a Photonic Jet

1
School of Physics and Optoelectronic Engineering, Xidian University, Xi’an 710071, China
2
Collaborative Innovation Center of Information Sensing and Understanding, Xidian University, Xi’an 710071, China
3
School of Nondestructive Testing, Tomsk Polytechnic University, 634050 Tomsk, Russia
4
China Research Institute of Radiowave Propagation, Qingdao 266000, China
*
Author to whom correspondence should be addressed.
Nanomaterials 2022, 12(2), 251; https://doi.org/10.3390/nano12020251
Submission received: 30 November 2021 / Revised: 29 December 2021 / Accepted: 4 January 2022 / Published: 13 January 2022
(This article belongs to the Special Issue Nanomaterial-Based Sharp Focusing and Application)

Abstract

:
In this article, we study the optical force exerted on nanorods. In recent years, the capture of micro-nanoparticles has been a frontier topic in optics. A Photonic Jet (PJ) is an emerging subwavelength beam with excellent application prospects. This paper studies the optical force exerted by photonic jets generated by a plane wave illuminating a Generalized Luneburg Lens (GLLs) on nanorods. In the framework of the dipole approximation, the optical force on the nanorods is studied. The electric field of the photonic jet is calculated by the open-source software package DDSCAT developed based on the Discrete Dipole Approximation (DDA). In this paper, the effects of the nanorods’ orientation and dielectric constant on the transverse force Fx and longitudinal force Fy are analyzed. Numerical results show that the maximum value of the positive force and the negative force are equal and appear alternately at the position of the photonic jet. Therefore, to capture anisotropic nanoscale-geometries (nanorods), it is necessary to adjust the position of GLLs continuously. It is worth emphasizing that manipulations with nanorods will make it possible to create new materials at the nanoscale.

1. Introduction

Optical capture by manipulating neutral atoms by lasers was first reported by Ashkin et al. in 1970 [1]. At present, optical manipulation can achieve the capture of living cells and organelles [2], DNA rotation [3], and chromosome surgery [4] with low damage, which is very important in modern medicine and biology. Later, optical capture was gradually applied to the manipulation of various nanomaterials. Nanomaterials and nanotechnology are already indispensable and important members of modern technology. For example, quantum dots play an important role in fluorescence detection, bioluminescence detection and other fields due to their unique photoelectric properties [5]. Semiconductor nanowires [6] have excellent performance in electronics and photonics due to their tunable direct bandgap and high carrier mobility. Anisotropic nanoscale-geometries, including nanorods [7], can be used in molecular imaging and photothermal cancer therapy. The capture and manipulation of nanomaterials are significant in the optoelectronic industry. Different arrangements of nanomaterials have different optical and physical properties [8], and optical force is an important way to realize the complex arrangement structure of nanomaterials. In 2006, Pelton et al. [9] reported three-dimensional optical capture of gold nanorods and enhanced optical force by longitudinal surface plasmon resonance. In 2008, Selhuber-Unkel [10] performed a quantitative analysis of the optical force on gold nanorods. The results showed that the interaction between an optical trap and a nanorod is related to the particle’s polarizability. In 2010, Tong et al. [11] used a linearly polarized near-infrared laser to arrange and rotate nanostructures, including nanorods, optically. In 2011, Ruijgrok et al. [12] quantitatively measured the torque exerted on the nanorod. In 2014, Liaw et al. [13] studied the polarizability of gold nanorods and analyzed the optical torque under the Maxwell stress tensor. In 2017, Fick et al. [14] captured nanorods using optical fiber nanotweezers. In 2020, Huang et al. [15] studied the capture of gold nanorods by a plasmonic tweezer. By adjusting the mode of local surface plasmon resonance, the force direction can be changed. Capturing nanorods and arranging them according to different patterns can enable nanoarrays to have new functions. In addition, the dipole approximation is a common and accurate method to study optical force. In this article, we use the dipole approximation method to numerically simulate the optical force exerted by the photonic jet on the nanorods.
When a beam of light illuminates a particle (its size is close to the wavelength), the backlit side of the particle will generate a strong and narrow beam. Scholars vividly named it a photonic jet. Photonic jet has been used in many fields. For example, optical data storage [16], superresolution imaging [17], and Raman signal enhancement [18], etc. At the same time, a photonic jet, as a highly focused beam, can also be used as a light source to capture and manipulate particles [19]. Its application in cell surgery [20,21,22] proves this point. In this paper, we will generate PJ by illuminating a mesoscale Generalized Luneburg Lens with a plane wave, which, however, does not reduce the generality of the problem, and the results can be generalized to other particle-lens combinations. GLLs is the particle model mentioned in Mao et al.’s [23] paper. The refractive index of this model is rotationally symmetrical and can generate an ultra-narrow or ultra-long PJ. This PJ is of great significance in the optical trapping of nanoparticles. We can change PJ significantly by changing the wavelength, focal length, and radius of GLLs, and then analyze the influence of PJ on the optical force.
PJ and optical force are studied in the framework of Discrete Dipole Approximation [24]. The basic idea of DDA comes from DeVoe’s writings. DDA has studied scattering by first dividing particles of any shape into dipole arrays (polarizable) whose size is negligible compared to the wavelength [25]. Then the interaction between the dipoles and the incident field is analyzed to solve the whole field after the scattering. Until today, DDA has been developed into a numerical simulation method of the electromagnetic field with high precision and fast operation speed. In particular, DDSCAT [26] is an open-source software package of DDA with simple operation methods and accurate calculation results. DDSCAT is especially useful for a particle with a special shape. The coordinate system for the dipole array inside the particle in DDSCAT is straightforward to understand and modify. When calculating the scattering of a particle with an irregular shape, we only need to find the surface function and combine it with the dipole coordinate system [27]. Therefore, DDSCAT is software suitable for calculating the scattered field of GLLs. Dipole approximation [28] is a common method for numerical simulating the optical force exerted on nanostructures. In this calculation model, the nanostructure is regarded as a dipole (the size is almost negligible). In the calculation of optical force, the shape of the nanostructure can be expressed by its polarizability. The polarizability of a nanoparticle with a more complex shape is often a tensor [29]. The expression of the polarizability of the nanorod in this article will be shown in the second section.
The rest of this paper is distributed as follows. The second section introduces the theory of DDA and DDSCAT to calculate the near-field scattering, and then analyzes the theory of force on nanorod in the dipole approximation framework. In the third section, the numerical results of transverse and longitudinal forces on nanorods and the effects of wavelength, focal length, and radius of GLLs on the optical forces are discussed. The fourth section summarizes the full article.

2. Materials and Methods

2.1. Discrete Dipole Approximation

In the above, we briefly introduced the method of DDA to calculate near-field scattering. In this section, we will analyze the principle of DDA in detail. The electric dipole moment of a polarizable (the polarizability is α i ) dipole is P i = α i E i [30]. The total electric field at r i in the calculation area can be expressed as [31]:
E i = E i n c , i + E s c a , i
where E i n c , i and E s c a , i are the incident electric field and scattered electric field at r i . A i j P j represents the scattering field generated by the dipole (at r j ) at r i . Therefore, Equation (1) can be expressed as:
E i = E i n c , i j i A i j P j
where, A i j [32] represents the influence of the dipole at position r j on the dipole at position r i . A i j = exp ( i k r i j ) r i j × [ k 2 ( r ^ i j r ^ i j 1 3 ) + i k r i j 1 r i j 2 ( 3 r ^ i j r ^ i j 1 3 ) ] , i j , where, k is the wave vector, r i j is the distance between point i and point j , and r ^ i j is the unit vector from i to j . 1 3 is the identity matrix. r ^ i j r ^ i j and 1 3 are shown in Equations (3) and (4) [33,34]. r x , r y , and r z are the coordinate components of r ^ i j along the x , y , and z directions, respectively. j i A i j P j represents the scattered field at that point:
r ^ i j r ^ i j = [ r x 2 r x r y r x r z r y r x r y 2 r y r z r z r x r z r y r z 2 ]
1 3 = [ 1 0 0 0 1 0 0 0 1 ]
Now, we can get the electric dipole moment at any position [35,36]:
P i = α i ( E i n c , i j i A i j P j )
E i n c , i = j i A i j P j + α i 1 P i
where [37]:
α i α C M 1 + ( α C M / d 3 ) [ ( b 1 + m 2 b 2 + m 2 b 3 S ) ( k d ) 2 ( 2 / 3 ) i ( k d ) 3 ] , b 1 = 1.891531 , b 2 = 0.1648469 , b 3 = 1.7700004 , S i = 1 3 ( a ^ i e ^ i ) 2 ,
α i C M = 3 d 3 4 π n i 1 n i + 2
where, α i C M is the Clausius-Mossotti polarizabilities (as shown in Equation (8)), m and n i are the refractive index and the dielectric constant at any position. d is the distance between dipoles, b 1 , b 2 , and b 3 are constants, k is the wave vector, a ^ and e ^ are the unit vectors of the incident direction and the polarization direction [38].
In this paper, n i = n 0 [ 1 + f G L L s 2 ( r 0 / R ) 2 ] 1 / 2 / f G L L s , where, n 0 = 1 . f G L L s are the focal length normalized radii of GLLs. It is unitless. r o is the radial coordinate ( o = 1, 2, 3......, 30) and R is the maximum radius of the GLLs. Due to the limitation of memory and computing speed, we only divide the GLLs into 30 layers, and we will verify the correctness of this model in Section 3.1. The design method and correctness verification of GLLs in DDSCAT have been completed in our previous work [39], and will not be described in this paper. Defining A j j α j 1 , Equation (6) has the following form [40]:
E inc ,   i = j   target   A i j P j
In DDSCAT, the inner (original target sites j) and outer (vacuum sites i) fields of a particle will be represented separately as [41]:
E = E i n c + E s c a t = { α j 1 P j   original   target   sites   j E i n c , i j t arg e t A i j P j   vacuum   sites   i
Finally, the electric dipole moment outside the particle is 0.

2.2. Optical Force on a Nanorod

In the framework of the dipole approximation, the polarizability of arbitrarily oriented nanorod is put into the optical force equation as follows [42]:
α z = V ε d ε m ε d ε d
α t = 2 V ε d ε m ε d ε m + ε d
where, V is the volume of nanorods, ε m and ε d are the dielectric functions of metal and medium, respectively. The polarizability tensor of the particle is [29]:
α ^ = | α z 0 0 0 α t 0 0 0 α t |
The coordinate system inside the particle is ( ξ , η , ζ ). The schematic diagram of the polarization distribution and orientation of the nanorods is shown in Figure 1. α z and α t represent the longitudinal (parallel to the axis) and transverse (parallel to the bottom surface) polarizability of the nanorod, respectively. Both the ξ and η axes are parallel to the bottom surface of the nanorod, so the polarizability along these two directions is the same. The rotation matrix connecting the two systems is [29,43]:
R ^ ( θ 0 , φ 0 ) = ( cos φ 0 cos θ 0 sin φ 0 cos φ 0 sin θ 0 sin φ 0 cos θ 0 cos φ 0 sin φ 0 sin θ 0 sin θ 0 0 cos θ 0 )
The polarizability tensor of the rotated nanorod in the system of coordinates xyz is then done by the matrix product [44]:
A ^ = R ^ 1 α ^ R ^
and the nanorod dipole moment p excited by the total electric field E at the nanorod position r = ( x , y , z ) is:
p = A ^ ( θ 0 , φ 0 ) E ( x , y , z )
Finally, the optical force on the nanorods is [45]:
F ξ = 1 2 { p ξ E * } , ξ = x , y , z
By rewriting p as:
p = [ p 1 p 2 p 3 ]
Equation (17) can be expanded to:
F x = 1 2 { p 1 E x * x + p 2 E y * x + p 3 E z * x } F y = 1 2 { p 1 E x * y + p 2 E y * y + p 3 E z * y } F z = 1 2 { p 1 E x * z + p 2 E y * z + p 3 E z * z }
In this paper, we will analyze the F x , longitudinal force F y and F z exerted by the photonic jet on a nanorod.

3. Results

This section will analyze the simulation results. We use DDSCAT to generate a GLLs containing 30 layers with different refractive indices. The internal structure of GLLs is shown in Figure 2. The thickness of each layer is almost the same, and the refractive index of each layer can be calculated by n i = n 0 [ 1 + f G L L s 2 ( r 0 / R ) 2 ] 1 / 2 / f G L L s . In all numerical simulations, the maximum radius R and focal length f G L L s of GLLs are 2 μm and 1.2 a.u., and the wavelength will be introduced separately in each section. The center of mass of GLLs is at ( 0 , 0 , 0 ) , and the plane wave always propagates along + x and is polarized along the y -direction. The reason for choosing these parameters is to ensure correctness (compared with Mie theoretical calculation results) and a better jet effect. The nanoparticles are gold nanorods, and their polarizability is given in Equations (11)–(13). As the calculation framework in this paper is dipole approximation, the volume of nanoparticles must meet specific requirements, and the effective radius of nanorods must meet certain requirements 2 π a e f f 1 . In this paper, we only discuss the effect of orientation and dielectric constant of nanorod on F x , F y , and F z . Both electric field (a.u.) and optical force (in N) are displayed in the x o y plane and the positions of GLLs are marked with white circles. When the electric field is shown separately, the value is | E | / | E 0 | . For convenience, we will use | E | to represent the value. Because the photonic jet is a superposition of the incident and scattering fields on the shadow side of the particle, we must know the ratio of the field in the area of the photonic jet to the incident field. There are two other important parameters for PJ, which are the focal length f and the full width at half maxima (FWHM) at the focal point. These two parameters describe the distance from the focal point of the PJ to the outer surface of the GLLs and the width where the intensity is half the intensity of the focal point, respectively. They respectively represent the horizontal and vertical scales of the PJ. When the FWHM is less than half of the wavelength, the PJ may break the diffraction limit, so these two scales tend to be measured in wavelength λ .
In the following, we first compare the photonic jet of GLLs under the DDA framework with that under the Mie theory framework. Here, we will only show the comparison results of photonic jet with wavelength 0.5 μm and 0.6328 μm. Then we analyzed the influence of the orientation of the nanorods on the optical force, and finally analyzed the influence of the wavelength on the optical force.
This section may be divided by subheadings. It should provide a concise and precise description of the experimental results, their interpretation, as well as the experimental conclusions that can be drawn.

3.1. Numerical Validation

In this section, we calculate the photonic jet generated by a plane wave irradiating a GLLs respectively through DDA and Mie theory and compare the results. The schematic diagram of the PJ irradiating nanorod is shown in Figure 3. DDSCAT and Jan Schäfer [46,47,48,49,50]’s Mie theory software were used, and the two results were normalized by I norm   = ( I i I min ) / ( I max I min ) [51], where, I i , I min and I max are the intensity at any position, the maximum value, and the minimum value, respectively. The error of all comparisons is within one-thousandth. As shown in Figure 4, the wavelengths are 0.5 (Figure 4a) and 0.6328 μm (Figure 4b) respectively, the maximum radius of GLLs is 2 μm (i.e., Mie size parameter q 20), and the focal length is 1.2 a.u..

3.2. Orientation of Nanorods

In this section, we will analyze the influence of nanorod’s orientation on the optical force. As mentioned earlier, the focal length f G L L s and the maximum radius of GLLs are 1.2 and 2 μm. The incident wavelength is 0.6328 μm, the dielectric constant of the gold nanorod is −11.740 + 1.2611i, and the volume is 4.2 × 10−3 μm3. We will first show the electric field of the photonic jet under this set of parameters, as shown in Figure 5. The intensity of the focal point of the photonic jet is 51.7 a.u., the focal length f is 0.62 λ, and the F W H M is 0.5 λ. Next, we will use two subsections to discuss the influence of the orientation φ 0 and θ 0 of the nanorods on the optical force.
The optical force in this paper has several discrete regions. For convenience, we define the independent region of optical force as discrete optical force, as shown in Figure 6. Figure 6a,b show two kinds of discrete optical forces, respectively.

3.2.1. φ 0 Orientation of Nanorods

In this section, we will study the influence of the change of φ 0 on the optical force ( θ 0 = 90 ° ). We respectively show F x and F y when φ 0 = 0 ° ,   90 ° ,   99 ° ,   and   105 ° . The reason why we choose these angles is because of the unique properties of the discrete optical force at these angles. There are significant differences in the size and shape of these discrete optical forces. Of course, we have calculated the optical forces at all angles. However, the distribution patterns of these optical forces at these four angles is the most representative. We will explain their unique properties one by one. When φ 0 is one of these four cases, the optical force (especially F y ) has a unique and typical distribution pattern of positive and negative forces. When φ 0 = 0 ° , F y is almost perfectly symmetric with respect to y = 0, including the magnitude and scope of positive and negative forces. When φ 0 = 90 ° , the F y in the y > 0 area is significantly smaller than the F y in the y < 0 area (reversed when φ 0 = 105 ° ). When φ 0 = 99 ° , the scope of positive force increases significantly, and the negative force converges on the particle surface. Still, the maximum value of the negative force is higher than that of the positive force. Because of the uniqueness of the optical forces in these four angles, we decided to show them, as shown in Figure 7a–d. When φ 0 is equal to other values, the distribution pattern of positive and negative forces is always similar to the above four cases, so it will not be shown. In increasing φ 0 from 0 to 180 degrees, the change of F x is mainly reflected in the intensity of the force. At the position of the PJ, the optical force is symmetrical about y = 0, and the positive force and the negative force alternately appear. The maximum value of the positive force and the negative force are almost equal, and both have a clear action area. Therefore, when capturing nanorods with different orientations in the light propagation direction ( x ), the negative force can always be exerted on the nanorods by accurately moving the GLLs to achieve the purpose of capture. Fy has different directions of action on both sides of y = 0, and the positive and negative forces also appear alternately along + x . Therefore, the nanorods oscillate repeatedly on both sides of y = 0. However, when φ 0 45 ° , the positive and negative forces on both sides of y = 0 appear asymmetry. The force on one side is more potent than the other side. Such particles may always move to one side in the y-direction. At this time, we can also move GLLs along the x-direction so that the nanorods always receive the same magnitude and opposite optical forces on both sides of y. This allows the nanorods to be in a state of dynamic equilibrium and to be stably captured in the y-direction. Figure 8 shows the variation of the maximum and minimum values of F x , F y , and F z with φ 0 . As shown in Figure 8a,b we calculate the change rules of the maximum and minimum values of F x , F y , and F z in the process of φ 0 increasing from 0 ° to 180 ° . Since the light propagates in the x direction and θ 0 = 90 ° , the radius of the bottom surface of the nanorod is small, so the intensity of F z is relatively tiny compared to F x and F y . Through research, it can be found that the maximum and minimum values of F x and F y have precisely the same trend with the increase of φ 0 . The optical force decreases first, reaching a minimum at 39 ° . Then increase, get a local maximum at 69 ° , and drop immediately. When φ 0 = 99 ° , it comes the local minimum value and increases. Reach the maximum value when φ 0 = 159 ° . In summary, different φ 0 will cause the force on nanorods to be very different, so we can propose targeted solutions for particles with different φ 0 to achieve stable capture.

3.2.2. θ 0 Orientation of Nanorods

In this section, we will discuss the influence of θ 0 on optical forces ( φ 0 = 45 ° ). As shown in Figure 9a–d, we offer F x and F y at θ 0 = 0 ° ,   30 ° ,   60 ° ,   and   180 ° , respectively. Similarly, these angles were chosen to demonstrate the unique properties of discrete optical forces and the distribution patterns of positive and negative forces. The distribution pattern of positive and negative forces determines the direction of force and the way to capture the particles. When θ 0 = 0 ° , the discrete optical force of F y is almost perfectly symmetric with respect to y = 0. When θ 0 = 30 ° , the discrete optical force of F y in the y < 0 area is smaller than the F y in the y > 0 area. When θ 0 = 60 ° , the scope of negative force increases significantly, and the positive force converges on the particle surface. Still, the maximum value of the positive force is higher than that of the negative force. When θ 0 = 180 ° , the range of negative force is greater than that of positive force, but the maximum value of negative force is smaller. Because of the uniqueness of the optical forces in these four angles, we decided to show them, as shown in Figure 9a–d. When θ 0 is other values, the distribution pattern of optical force is similar to the above four situations, and the repeated content will not be shown in this article. θ 0 in 0° − 42°, F x and F y are almost all symmetric with respect to y = 0, with only variations in the intensity of the force. Still, the range of F y changes significantly around 30 ° , 60 ° and 180 ° , as shown in Figure 9b–d. F x is not symmetric about y = 0. The range of the F x ’s negative force is significantly increased. The wide distribution of negative forces allows particles to be captured more easily. Figure 10 shows the variation of the maximum and minimum values of F x , F y , and F z with θ 0 . As shown in Figure 10a,b, we studied the maximum values of the positive and negative forces of F x , F y , and F z in the range of 0 ° θ 0 180 ° . The numerical results show that the maximum values of F x and F y decrease first, then increase and finally decrease again with the increase of θ 0 . When θ 0 = 120 ° , F z is a local minimum value. Otherwise, the trend of F z increasing with θ 0 is basically the same as that of F x . The maximum value of positive force and negative force are almost equal. Under the same parameters, the area where the maximum positive force and the negative force are located in the photon ejection is fixed. By moving the GLLs, the nanorods can be accurately captured.

3.3. Dielectric Constant

In this section, we discuss the effects of the real ( ε 1 ) and imaginary ( ε 2 ) parts of the nanorod’s dielectric constant on the optical force, as shown in Figure 11. The f G L L s and R of GLLs are 1.2 and 2 μm. The wavelength is 0.6328 μm, and the volume of the nanorod is 4.2 × 10−3 μm3. In Figure 11a,b we show the maximum and minimum optical force changes with the real part ( ε 1 ) of the dielectric constant when the imaginary part ε 2 of the nanorod’s dielectric constant is 1, respectively. From the figures, we can find that the optical force first increases and then decreases with the increase of ε 1 . When ε 1 = 2 , the optical force has the local maximum value. From Figure 11c,d, we can find that as ε 2 increases, the optical force gradually increases.

4. Conclusions

In this paper, we studied the optical forces F x and F y exerted by a photonic jet on nanorods under the framework of the dipole approximation. The PJ is generated by a plane wave irradiating a GLLs, and the numerical result is calculated by DDSCAT. We divided GLLs into 30 layers of concentric rings with different refractive indices. We compare the calculation results of DDSCAT with the results of the Mie theory to validate the PJ’s electric field. Then we use the PJ to illuminate the nanorods and study the optical forces Fx and Fy on the nanorods. Of course, the size of nanorods must satisfy the Rayleigh approximation. The effects of the orientation and dielectric constant of nanorods on the optical force are investigated. The angle between the axis of the nanorod and the z-axis is represented by θ 0 , and the angle between the projection of the axis on the xoy plane and the x-axis is defined by φ 0 . We first discuss φ 0 ’s effect on F x and F y by making θ 0 = 90 ° (the axis of the nanorods is in the xoy plane). We find that F x is always symmetric about the y-axis, and F y is not symmetric about the y-axis, only around φ 0 = 90 ° and φ 0 = 105 ° . In contrast, the intensity of F z can be ignored since the maximum value of the positive and negative optical forces are almost equal. It is necessary to constantly adjust the position of the GLLs when capturing these oriented nanorods so that the nanorods are always in the area of negative or positive forces. Then we investigate the effect of θ 0 on optical force when φ 0 = 45 ° . The numerical results show that F x is asymmetric about the y-axis near θ 0 = 45 ° and 60 ° , and the area of the negative force of F y increases significantly (including θ 0 = 180 ° ), and the positive force converges on the surface of GLLs. Similarly, the maximum values of the positive and negative optical forces are equal. Due to the wider distribution of negative forces at certain angles, it is easier to adjust GLLs to capture nanorods. Finally, we discussed the effect of the dielectric constant of the nanorods on the optical force. Under the same conditions, a larger dielectric constant generally leads to a powerful optical force. In this paper, the optical force exerted by the photonic jet on elongated particles in the form of nanorods instead of spherical nanoparticles with different orientations and materials (dielectric constant) is studied. The numerical results describe the direction and magnitude of the optical force, which makes it possible to directionally manipulate the nanorods. Note that the presented results will be valid not only for GLLs, but also for other particles [52] that form a photonic jet with similar characteristics. These results are expected to provide theoretically support for the manipulation of nanorods and the arrangement of nanoarrays.

Author Contributions

Conceptualization, R.L.; Data curation, L.L.; Funding acquisition, S.G. and R.L.; Investigation, L.L.; Supervision, S.G., I.V.M. and O.V.M.; Visualization, B.W.; Writing—original draft, B.W.; Writing—review and editing, R.L., I.V.M. and O.V.M.; Initiated the work, I.V.M. and O.V.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China under Grant (61971385) and the National Natural Science Foundation of China (Grant No: 61771375).

Data Availability Statement

Some or all data, models, or code generated or used during the study are available in a repository or online in accordance with funder data retention policies (Provide full citations that include URLs or DOIs.)

Acknowledgments

I.V.M. and O.V.M. acknowledge support from the Tomsk Polytechnic University Development Program, Russia.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
PJPhotonic jet
GLLsGeneralized Luneburg Lens
DDADiscrete Dipole Approximation
FWHMFull Width at Half Maxima

References

  1. Ashkin, A. Atomic-beam deflection by resonance-radiation pressure. Phys. Rev. Lett. 1970, 25, 1321–1324. [Google Scholar] [CrossRef]
  2. Ambardekar, A.A.; Li, Y.Q. Optical levitation and manipulation of stuck particles with pulsed optical tweezers. Opt. Lett. 2005, 301, 797–1799. [Google Scholar] [CrossRef] [PubMed]
  3. Oroszi, L.; Galajda, P.; Kirei, H.; Bottka, S.; Ormos, P. Direct measurement of torque in an optical trap and its application to double-strand dna. Phys. Rev. Lett. 2006, 97, 058301. [Google Scholar] [CrossRef]
  4. Ponelies, N.; Scheef, J.; Harim, A.; Leitz, G.; Greulich, K. Laser micromanipulators for biotechnology and genome research. J. Biotechnol. 1994, 35, 109–120. [Google Scholar] [CrossRef]
  5. Yao, J.; Li, L.; Li, P.; Yang, M. Quantum dots: From fluorescence to chemiluminescence, bioluminescence, electrochemiluminescence, and electrochemistry. Nanoscale 2017, 9, 13364–13383. [Google Scholar] [CrossRef] [PubMed]
  6. Fang, M.; Han, N.; Wang, F.; Yang, Z.X.; Yip, S.; Dong, G.; Hou, J.J.; Chueh, Y.; Ho, J.C. Cheminform abstract: Iii—v nanowires: Synthesis, property manipulations, and device applications. ChemInform 2015, 46. [Google Scholar] [CrossRef]
  7. Huang, X.; El-Sayed, I.H.; Qian, W.; El-Sayed, M.A. Cancer cell imaging and photothermal therapy in the near-infrared region by using gold nanorods. J. Am. Chem. Soc. 2006, 1282, 115–2120. [Google Scholar] [CrossRef]
  8. Maragò, O.M.; Jones, P.H.; Gucciardi, P.G.; Volpe, G.; Ferrari, A.C. Optical trapping and manipulation of nanostructures. Nat. Nanotechnol. 2013, 88, 07–819. [Google Scholar] [CrossRef] [Green Version]
  9. Pelton, M.; Liu, M.; Kim, H.Y.; Smith, G.; Guyot-Sionnest, P.; Scherer, N.F. Optical trapping and alignment of single gold nanorods by using plasmon resonances. Opt. Lett. 2006, 312, 2075–2077. [Google Scholar] [CrossRef]
  10. Selhuber-Unkel, C.; Zins, I.; Schubert, O.; Sönnichsen, C.; Oddershede, L.B. Quantitative optical trapping of single gold nanorods. Nano Lett. 2008, 82, 998–3003. [Google Scholar] [CrossRef]
  11. Tong, L.; Miljkovi´c, V.D.; Käll, M. Alignment, rotation, and spinning of single plasmonic nanoparticles and nanowires using polarization dependent optical forces. Nano Lett. 2010, 102, 68–273. [Google Scholar] [CrossRef]
  12. Ruijgrok, P.V.; Verhart, N.R.; Zijlstra, P.; Tchebotareva, A.L.; Orrit, M. Brownian fluctuations and heating of an optically aligned gold nanorod. Phys. Rev. Lett. 2011, 107, 037401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Liaw, J.W.; Lo, W.J.; Kuo, M.K. Wavelength-dependent longitudinal polarizability of gold nanorod on optical torques. Opt. Express 2014, 221, 0858–10867. [Google Scholar] [CrossRef] [PubMed]
  14. Fick, J.; Leménager, G.; Thiriet, M.; Lallid, K.; Gacoin, T.; Valdivia-Valero, F.; Colas des Francs, G. Trapping of rare earth-doped nanorods with high aspect ratios using optical fiber-tip nano-tweezers. In Proceedings of the 2017 Conference on Lasers and Electro-Optics Europe European Quantum Electronics Conference, Munich, Germany, 25–29 June 2017. [Google Scholar]
  15. Huang, W.H.; Li, S.F.; Xu, H.T.; Xiang, Z.X.; Long, Y.B.; Deng, H.D. Tunable optical forces enhanced by plasmonic modes hybridization in optical trapping of gold nanorods with plasmonic nanocavity. Opt. Express 2018, 266, 202–6213. [Google Scholar] [CrossRef]
  16. Kong, S.C.; Sahakian, A.; Taflove, A.; Backman, V. Photonic nanojet-enabled optical data storage. Opt. Express 2008, 161, 3713–13719. [Google Scholar] [CrossRef] [PubMed]
  17. Chen, L.; Zhou, Y.; Li, Y.; Hong, M. Microsphere enhanced optical imaging and patterning: From physics to applications. Appl. Phys. Rev. 2019, 60, 21304. [Google Scholar] [CrossRef]
  18. Upputuri, P.K.; Wu, Z.; Gong, L.; Ong, C.K.; Wang, H. Super-resolution coherent anti-stokes raman scattering microscopy with photonic nanojets. Opt. Express 2014, 221, 2890–12899. [Google Scholar] [CrossRef]
  19. Minin, I.; Minin, O. Recent trends in optical manipulation inspired by mesoscale photonics and diffraction optics. J. Biomed. Photonics Eng. 2020, 6, 020301. [Google Scholar] [CrossRef]
  20. Minin, I.V.; Minin, O.V.; Liu, Y.Y.; Tuchin, V.V.; Liu, C.Y. Concept of photonic hook scalpel generated by shaped fiber tip with asymmetric radiation. J. Biophotonics 2021, 14, e202000342. [Google Scholar] [CrossRef]
  21. Ferrari, H.; Renner, G.; Luebrecht, S.; Love, G. High-frequency jet ventilation: Applications for endoscopy and surgery of the airway. South. Med. J. 1986, 799, 41–943. [Google Scholar] [CrossRef]
  22. Darafsheh, A.; Fardad, A.; Fried, N.M.; Antoszyk, A.N.; Ying, H.S.; Astratov, V.N. Contact focusing multimodal microprobes for ultraprecise laser tissue surgery. Opt. Express 2011, 193, 440–3448. [Google Scholar] [CrossRef]
  23. Mao, X.; Yang, Y.; Dai, H.; Luo, D.; Yao, B.; Yan, S. Tunable photonic nanojet formed by Generalized Luneburg Lens. Opt. Express 2015, 232, 6426–26433. [Google Scholar] [CrossRef]
  24. DeVoe, H. Optical properties of molecular aggregates. I. classical model of electronic absorption and refraction. J. Chem. Phys. 1964, 413, 93–400. [Google Scholar] [CrossRef]
  25. DeVoe, H. Optical properties of molecular aggregates. II. Classical theory of the refraction, absorption, and optical activity of solutions and crystals. J. Chem. Phys. 1965, 433, 199–3208. [Google Scholar] [CrossRef]
  26. Yon, J.; Liu, F.; Morán, J.; Fuentes, A. Impact of the primary particle polydispersity on the radiative properties of soot aggregates. Proc. Combust. Inst. 2019, 371, 151–1159. [Google Scholar] [CrossRef]
  27. Yurkin, M.A.; de Kanter, D.; Hoekstra, A.G. Accuracy of the discrete dipole approximation for simulation of optical properties of gold nanoparticles. J. Nanophotonics 2010, 4, 041585. [Google Scholar]
  28. Gordon, J.P. Radiation forces and momenta in dielectric media. Phys. Rev. A 1973, 8, 14–21. [Google Scholar] [CrossRef]
  29. Li, M.; Yan, S.; Yao, B.; Lei, M.; Yang, Y.; Min, J.; Dan, D. Trapping of rayleigh spheroidal particles by highly focused radially polarized beams. J. Opt. Soc. Am. B 2015, 324, 468–472. [Google Scholar] [CrossRef]
  30. Draine, B.T.; Flatau, P.J. Discrete-dipole approximation for scattering calculations. J. Opt. Soc. Am. A 1994, 111, 491–1499. [Google Scholar] [CrossRef]
  31. Yang, W.; Schatz, G.C.; Van Duyne, R.P. Discrete dipole approximation for calculating extinction and raman intensities for small particles with arbitrary shapes. J. Chem. Phys. 1995, 1038, 869–875. [Google Scholar] [CrossRef] [Green Version]
  32. Hoekstra, A.G.; Frijlink, M.; Waters, L.; Sloot, P. Radiation forces in the discrete-dipole approximation. JOSA A 2001, 181, 944–1953. [Google Scholar] [CrossRef] [Green Version]
  33. Amendola, V. Surface plasmon resonance of silver and gold nanoparticles in the proximity of graphene studied using the discrete dipole approximation method. Phys. Chem. Chem. Phys. 2016, 18, 2230–2241. [Google Scholar] [CrossRef] [PubMed]
  34. Jensen, T.; Kelly, L.; Lazarides, A.; Schatz, G. Electrodynamics of noble metal nanoparticles and nanoparticle clusters. J. Clust. Sci. 1999, 102, 95–317. [Google Scholar]
  35. Collinge, M.J.; Draine, B.T. Discrete-dipole approximation with polarizabilities that account for both finite wavelength and target geometry. JOSA A 2004, 212, 2023–2028. [Google Scholar] [CrossRef] [Green Version]
  36. Yurkin, M.A.; Maltsev, V.P.; Hoekstra, A.G. The discrete dipole approximation for simulation of light scattering by particles much larger than the wavelength. J. Quant. Spectrosc. Radiat. Transf. 2007, 106, 546–557. [Google Scholar] [CrossRef] [Green Version]
  37. Loke, V.L.; Mengüç, M.P.; Nieminen, T.A. Discrete-dipole approximation with surface interaction: Computational toolbox for matlab. J. Quant. Spectrosc. Radiat. Transf. 2011, 1121, 1711–1725. [Google Scholar] [CrossRef]
  38. Draine, B.T.; Goodman, J. Beyond Clausius-Mossotti: Wave Propagation on a Polarizable Point Lattice and the Discrete Dipole Approximation. Astrophys. J. 1993, 405, 685. [Google Scholar] [CrossRef]
  39. Wei, B.; Xu, Q.; Li, R.; Zhang, S.; Gong, S.; Sun, H.; Song, N. Optical torque on a rayleigh particle by photonic jet. J. Quant. Spectrosc. Radiat. Transf. 2021, 272, 107775. [Google Scholar] [CrossRef]
  40. Chaumet, P.C.; Rahmani, A. Efficient iterative solution of the discrete dipole approximation for magnetodielectric scatterers. Opt. Lett. 2009, 349, 917–919. [Google Scholar] [CrossRef] [PubMed]
  41. Yurkin, M.; Hoekstra, A. The discrete dipole approximation: An overview and recent developments. J. Quant. Spectrosc. Radiat. Transf. 2007, 1065, 558–589. [Google Scholar] [CrossRef] [Green Version]
  42. Spector, M.; Ang, A.S.; Minin, O.V.; Minin, I.V.; Karabchevsky, A. Photonic hook formation in near-infrared with Mxene Ti3C2 nanoparticles. Nanoscale Adv. 2020, 2, 5312–5318. [Google Scholar] [CrossRef]
  43. Li, M.; Yan, S.; Yao, B.; Liang, Y.; Han, G.; Zhang, P. Optical trapping force and torque on spheroidal rayleigh particles with arbitrary spatial orientations. J. Opt. Soc. Am. A 2016, 331, 341–1347. [Google Scholar] [CrossRef]
  44. Trojek, J.; Chvátal, L.; Zemánek, P. Optical alignment and confinement of an ellipsoidal nanorod in optical tweezers: A theoretical study. J. Opt. Soc. Am. A 2012, 291, 224–1236. [Google Scholar] [CrossRef]
  45. Mitri, F.G. Adjustable vector Airy light-sheet single optical tweezers: Negative radiation forces on a subwavelength spheroid and spin torque reversal. Eur. Phys. J. D 2018, 722, 1–14. [Google Scholar] [CrossRef]
  46. Schäfer, J.P. Implementierung und Anwendung Analytischer und Numerischer Verfahren zur Lösung der Maxwellgleichungen für die Untersuchung der Lichtausbreitung in Biologischem Gewebe. Ph.D. Thesis, Univerität Ulm, Ulm, Germany, 2011. [Google Scholar]
  47. Schäfer, J.; Lee, S.C.; Kienle, A. Calculation of the near fields for the scattering of electromagnetic waves by multiple infinite cylinders at perpendicular incidence. J. Quant. Spectrosc. Radiat. Transf. 2012, 1132, 2113–2123. [Google Scholar] [CrossRef]
  48. Bohren, C.F.; Huffman, D.R. Absorption and Scattering of Light by Small Particlesl; John Wiley & Sons: Hoboken, NJ, USA, 2008. [Google Scholar]
  49. Lee, S. Dependent scattering of an obliquely incident plane wave by a collection of parallel cylinders. J. Appl. Phys. 1990, 684, 952–4957. [Google Scholar] [CrossRef]
  50. Kerker, M. Electromagnetic scattering (book reviews: The scattering of light and other electromagnetic radiation). Science 1970, 167, 861. [Google Scholar]
  51. Dobler, C.P.; Rosoff, J. SPSS for Introductory Statistics: Use and Interpretation/SPSS for Intermediate Statistics: Use and Interpretation. Am. Stat. 2005, 59, 352. [Google Scholar]
  52. Minin, O.V.; Minin, I.V. Optical phenomena in mesoscale dielectric particles. Photonics 2021, 8, 591. [Google Scholar] [CrossRef]
Figure 1. (a,b) respectively show the nanorod’s coordinate system ( ξ , η , ζ ) and the angle between the axis ( ξ ) of the nanorod and the ( x , y , z ) coordinate system ( θ 0 , φ 0 ). θ 0 represents the angle between the axis of the nanorod and the z-axis, and φ 0 represents the angle between the projection of the axis of the nanorod on the xoy plane and the x-axis.
Figure 1. (a,b) respectively show the nanorod’s coordinate system ( ξ , η , ζ ) and the angle between the axis ( ξ ) of the nanorod and the ( x , y , z ) coordinate system ( θ 0 , φ 0 ). θ 0 represents the angle between the axis of the nanorod and the z-axis, and φ 0 represents the angle between the projection of the axis of the nanorod on the xoy plane and the x-axis.
Nanomaterials 12 00251 g001
Figure 2. Internal structure diagram of GLLs. We have marked 28, 29, and 30 layers of concentric rings with different refractive indices in the figure. We use r1r28 to represent concentric rings from 1 to 28 layers, R represents the maximum radius of GLLs, and fGLLs represents the focal length of GLLs.
Figure 2. Internal structure diagram of GLLs. We have marked 28, 29, and 30 layers of concentric rings with different refractive indices in the figure. We use r1r28 to represent concentric rings from 1 to 28 layers, R represents the maximum radius of GLLs, and fGLLs represents the focal length of GLLs.
Nanomaterials 12 00251 g002
Figure 3. PJ generated by a plane wave illuminating a GLLs. In the figure, plane waves, GLLs, photonic jets, and nanorods and some of their parameters (such as the maximum radius R and the focal length fGLLs of GLLs, the focal point of the photonic jet, the focal length f and FWHM) are marked respectively. GLLs are composed of 30 layers of concentric rings with different refractive indexes centered on the center of the sphere.
Figure 3. PJ generated by a plane wave illuminating a GLLs. In the figure, plane waves, GLLs, photonic jets, and nanorods and some of their parameters (such as the maximum radius R and the focal length fGLLs of GLLs, the focal point of the photonic jet, the focal length f and FWHM) are marked respectively. GLLs are composed of 30 layers of concentric rings with different refractive indexes centered on the center of the sphere.
Nanomaterials 12 00251 g003
Figure 4. Comparison of theoretical calculation results between DDSCAT and Mie. (a,b) describe the comparison results of the two methods when the wavelength is 0.5 μm and 0.6328 μm, respectively. The first picture in each row is the intensity distribution along the x-axis through the focal point of the photonic jet. The second picture is also the intensity distribution but along the y-axis through the focal point.
Figure 4. Comparison of theoretical calculation results between DDSCAT and Mie. (a,b) describe the comparison results of the two methods when the wavelength is 0.5 μm and 0.6328 μm, respectively. The first picture in each row is the intensity distribution along the x-axis through the focal point of the photonic jet. The second picture is also the intensity distribution but along the y-axis through the focal point.
Nanomaterials 12 00251 g004aNanomaterials 12 00251 g004b
Figure 5. PJ generated by a plane wave with a wavelength of 0.6328 μm irradiating a GLLs.
Figure 5. PJ generated by a plane wave with a wavelength of 0.6328 μm irradiating a GLLs.
Nanomaterials 12 00251 g005
Figure 6. (a,b) The different kinds of discrete optical forces.
Figure 6. (a,b) The different kinds of discrete optical forces.
Nanomaterials 12 00251 g006
Figure 7. Optical force exerted on a nanorod with different φ0. Panels (ad) show the optical forces Fx and Fy when φ 0 = 0 ° ,   90 ° ,   99 ° ,   and   105 ° , respectively.
Figure 7. Optical force exerted on a nanorod with different φ0. Panels (ad) show the optical forces Fx and Fy when φ 0 = 0 ° ,   90 ° ,   99 ° ,   and   105 ° , respectively.
Nanomaterials 12 00251 g007aNanomaterials 12 00251 g007b
Figure 8. The maximum and minimum values of optical force varies with φ 0 . (a,b) show the change rules of the maximum and minimum values of F x , F y , and F z in the process of φ 0 increasing from 0 ° to 180 ° .
Figure 8. The maximum and minimum values of optical force varies with φ 0 . (a,b) show the change rules of the maximum and minimum values of F x , F y , and F z in the process of φ 0 increasing from 0 ° to 180 ° .
Nanomaterials 12 00251 g008
Figure 9. Optical force exerted on a nanorod with different θ0. Panels (ad) show the optical forces Fx and Fy when θ 0 = 0 ° ,   30 ° ,   60 ° ,   and   180 ° , respectively.
Figure 9. Optical force exerted on a nanorod with different θ0. Panels (ad) show the optical forces Fx and Fy when θ 0 = 0 ° ,   30 ° ,   60 ° ,   and   180 ° , respectively.
Nanomaterials 12 00251 g009aNanomaterials 12 00251 g009b
Figure 10. The maximum and minimum values of optical force varies with θ 0 . (a,b) show the change rules of the maximum and minimum values of F x , F y , and F z in the process of θ0 increasing from 0° to 180°.
Figure 10. The maximum and minimum values of optical force varies with θ 0 . (a,b) show the change rules of the maximum and minimum values of F x , F y , and F z in the process of θ0 increasing from 0° to 180°.
Nanomaterials 12 00251 g010
Figure 11. The maximum and minimum values of optical force varies with εm. Panels (ad) describe the effects of the real and imaginary parts of the nanorod’s dielectric constant on the maximum and minimum values of optical force, respectively.
Figure 11. The maximum and minimum values of optical force varies with εm. Panels (ad) describe the effects of the real and imaginary parts of the nanorod’s dielectric constant on the maximum and minimum values of optical force, respectively.
Nanomaterials 12 00251 g011
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wei, B.; Gong, S.; Li, R.; Minin, I.V.; Minin, O.V.; Lin, L. Optical Force on a Metal Nanorod Exerted by a Photonic Jet. Nanomaterials 2022, 12, 251. https://doi.org/10.3390/nano12020251

AMA Style

Wei B, Gong S, Li R, Minin IV, Minin OV, Lin L. Optical Force on a Metal Nanorod Exerted by a Photonic Jet. Nanomaterials. 2022; 12(2):251. https://doi.org/10.3390/nano12020251

Chicago/Turabian Style

Wei, Bojian, Shuhong Gong, Renxian Li, Igor V. Minin, Oleg V. Minin, and Leke Lin. 2022. "Optical Force on a Metal Nanorod Exerted by a Photonic Jet" Nanomaterials 12, no. 2: 251. https://doi.org/10.3390/nano12020251

APA Style

Wei, B., Gong, S., Li, R., Minin, I. V., Minin, O. V., & Lin, L. (2022). Optical Force on a Metal Nanorod Exerted by a Photonic Jet. Nanomaterials, 12(2), 251. https://doi.org/10.3390/nano12020251

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop