Nothing Special   »   [go: up one dir, main page]

Next Article in Journal
Modeling of High-Power Graded-Index Fiber Amplifiers
Previous Article in Journal
The Rapid Detection of Paclitaxel-Induced Changes in Cervical Cancer Cells Using an Ultrasensitive Biosensor
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tunable Photonic Hook Design Based on Anisotropic Cutting Liquid Crystal Microcylinder

1
School of Physics, Xidian University, Xi’an 710071, China
2
Department of Physics, Government College University, Faisalabad 38000, Pakistan
3
Department of Physics, University of Agriculture, Faisalabad 38000, Pakistan
4
Department of Mathematics, Lahore Campus, COMSATS University Islamabad, Lahore 53713, Pakistan
5
Nondestructuve Testing School, Tomsk Polytechnic University, Lenina 36, Tomsk 634050, Russia
*
Author to whom correspondence should be addressed.
Photonics 2024, 11(8), 736; https://doi.org/10.3390/photonics11080736
Submission received: 26 June 2024 / Revised: 12 July 2024 / Accepted: 18 July 2024 / Published: 7 August 2024
(This article belongs to the Special Issue Vortex Beams: Transmission, Scattering and Application)
Figure 1
<p>The FDTD model that waves illuminate an anisotropic cutting liquid crystal microcylinder in water.</p> ">
Figure 2
<p>Schematic diagram of the photonic hook. The start point, inflection point, and end point are colored green, purple, and red. The inflection point corresponds to the maximum electric field intensity <math display="inline"><semantics> <msub> <mi>I</mi> <mo movablelimits="true" form="prefix">max</mo> </msub> </semantics></math>.</p> ">
Figure 3
<p>The PH is generated by an anisotropic cutting microcylinder with different extraordinary refractive indices irradiated with plane waves. The first and second columns are the distribution figures of electric field intensity. The third and fourth columns represent the field intensities in <span class="html-italic">y</span> and <span class="html-italic">x</span> directions, respectively.</p> ">
Figure 4
<p>Energy distribution of a plane wave incident on the cutting LC microcylinder.</p> ">
Figure 5
<p>Upon incidence of a plane wave, the bending angle, electric field intensity, and FWHM exhibit variations in response to changes in the <math display="inline"><semantics> <msub> <mi>n</mi> <mi mathvariant="normal">e</mi> </msub> </semantics></math>. (<b>a</b>) Bending angle (<b>b</b>) Electric field intensity (<b>c</b>) FWHM.</p> ">
Figure 6
<p>The bimodal phenomenon when a beam illuminates the cutting LC microcylinder.</p> ">
Figure 7
<p>Schematic diagram of PH generation by Gaussian beam illumination on cutting microcylinder particles.</p> ">
Figure 8
<p>The PH generated by anisotropic cutting microcylinder particles with different extraordinary refractive indices irradiated with a Gaussian beam. The first and second columns are the distribution figures of electric field intensity. The third and fourth columns represent the field intensities in <span class="html-italic">y</span> and <span class="html-italic">x</span> directions, respectively.</p> ">
Figure 9
<p>Upon incidence of a Gaussian beam, the bending angle, electric field intensity, and FWHM exhibit variations in response to changes in the <math display="inline"><semantics> <msub> <mi>n</mi> <mi mathvariant="normal">e</mi> </msub> </semantics></math>. (<b>a</b>) Bending angle (<b>b</b>) Electric field intensity (<b>c</b>) FWHM.</p> ">
Figure 10
<p>Schematic diagram of PH generation by Bessel beam illumination on cutting microcylinder particles.</p> ">
Figure 11
<p>The PH is generated by anisotropic cutting microcylinder particles with different extraordinary refractive indices irradiated with a Bessel beam. The first and second columns are the distribution figures of electric field intensity. The third and fourth columns represent the field intensities in <span class="html-italic">y</span> and <span class="html-italic">x</span> directions, respectively.</p> ">
Figure 12
<p>When a Bessel beam is incident, the bending angle, electric field intensity, and FWHM change in response to variations in the <math display="inline"><semantics> <msub> <mi>n</mi> <mi mathvariant="normal">e</mi> </msub> </semantics></math>. (<b>a</b>) Bending angle (<b>b</b>) Electric field intensity (<b>c</b>) FWHM.</p> ">
Review Reports Versions Notes

Abstract

:
The selective control and manipulation of nanoparticles require developing and researching new methods for designing optical tweeters, mainly based on a photonic hooks (PHs) effect. This paper first proposes a tunable PH in which a structured beam illuminates an anisotropic cutting liquid crystal microcylinder based on the Finite-DifferenceTime-Domain (FDTD) method. The PHs generated by plane wave, Gaussian, and Bessel beam are analyzed and compared. The impact of beams and LC particle parameters on the PHs are discussed. Where the influence of the extraordinary refractive index ( n e ) on PHs is emphasized. Our results reveal that introducing birefringence can change the bending direction of PH. Besides, the maximum intensity of the PHs increases as n e increases regardless of the beam type. The PH generated by a plane wave has a higher maximum intensity and smaller FWHM than that generated by the Gaussian and Bessel beams. The smallest FWHM and maximum intensity of the PHs generated by the Gaussian falls between that generated by the plane wave and the Bessel beam. The PH generated by a Bessel beam has the minor maximum intensity and the largest FWHM. Still, it exceeds the diffraction limit and exhibits bending twice due to its self-recovery property. This paper provides a new way to modulate PH. This work offers novel theoretical models and the degree of freedom for the design of PHs, which is beneficial for the selective manipulation of nanoparticles. It has promising applications in Mesotronics and biomedicine.

1. Introduction

The optical tweezer [1] is an experimental technology that uses the interaction of light and matter to achieve particle trapping and manipulation. It was initially designed to facilitate the non-contact and non-destructive manipulation of tiny particles. Much research [2,3,4] has been conducted since Ashkin [5] proposed optical tweezers in the 1980s. Over the years, it played an important role in materials science [6,7], nanophotonics [8,9], and biomedicine [10,11]. Besides, it has been widely used in optical imaging [12], optical storage [13], and particle manipulation [14,15]. In recent years, a new manipulation method based on the photonic nanojet (PNJ) [16] has promoted the development of an optical tweezer. PNJ has the characteristics of much higher intensity than an incident wave, transverse half-height full width of subwavelength, and propagation to non-evanescent field region. It has shown high application value in the fields of super-resolution imaging [17,18], biochemical analysis [19], and nanoscale manipulation [20]. However, as an axial beam, its direction is limited to the axial, making it challenging to achieve selective particle manipulation. Airy beam [21,22] is a non-diffractive beam with self-recovery and self-bending characteristics. It can trap particles from any angle by passing obstacles. However, it is generated by bulky static or dynamic devices with demanding production conditions. It makes particle capture expensive.
In 2016, a curved PNJ was proposed and named Photonic Hook (PH) [23], which opens up new possibilities for achieving selective manipulation of particles. PH has the advantages of the Airy beam and PNJ, with high focusing, a multi-wavelength focal length distance, sub-diffraction full-width at half maximum (FWHM), limited diffraction, and self-bending characteristics [24]. Most importantly, it can detect and influence targets behind obstacles, and the achievement process is more straightforward than that of Airy beams. The research of PH has become a hot issue, and extensive achievements have been produced around it. A study [25] discusses the optical radiation force of PH generated by asymmetric cuboids, finding that PH can manipulate target particles behind obstacles. Following [23] in [26] describe the generation of PH based on a dielectric trapezoidal particle irradiated by a plane wave. Based on this model, a study successfully generated PH [27], confirming the existence of a PH experimentally for the first time. It makes the PH research from theoretical to application and provides strong support for the application of PH. Besides, the proposed concept of photonic hook scalpel [28], which uses an asymmetric fiber to generate PH at the tip, and the related results have broad application prospects in the development of laser surgery. Since then, a study [29] has experimentally generated a plasmonic PH using the dielectric Janus particles, further expanding the field of PH applications. With the further development of the research, a photonic hook generated by a multi-dielectric structure composed of periodically arranged scattering units was proposed [30], which simplified the PH generation process and provided innovative applications in the field of nanolithography and integrated optics. Recently, based on the optical Magnus effect, a novel method for generating PH using the interaction of plane waves and rotating particles has been proposed [31,32]. It illustrates the physical mechanism of PH production from an entirely new perspective. Up to now, the PH has generated considerable attention in diverse domains, such as optical nanomanipulation [33] and high-resolution imaging [34,35,36].
Although extensive research [37,38,39,40] has been carried out on the PNJ generated by the interaction between particles and various structured beams [41,42,43], the research on PH mainly focuses on the scattering of plane waves by particles. Previous studies on PNJs found that the particle and incident beam parameters closely affected the near field. Studying the scattering of particles to the structured beam can provide more modulation degrees of freedom and is expected to achieve selective manipulation of particles. Recently, a study [44] proposed a method to generate a twin photonic hook by irradiating an ellipsoid particle with a Bessel beam. Another research [45] explored the PH produced from a spherical particle illuminated by a Bessel–Gaussian beam. This research concludes that the effective length, intensity, and FWHM of PH can be modulated by varying the beam parameters. Therefore, it is vital for the design of PHs optical tweezers [25,46] to introduce structured beams into the PHs system. It will facilitate the development of a new generation of optical tweezers.
To the best of our comprehension, most of the current research is based on isotropic particles [47,48,49]. Compared with isotropic particles, the shape and refractive index of anisotropic particles [50,51] show asymmetry, which can generate ample torque through the light field to achieve accurate directional control. In addition, the optical anisotropy of anisotropic particles can be designed and controlled for optical detection and directional trap of particles. A study proposed controlled photonic hooks produced by a typical shape composed of uniaxial crystals [52]. It is found that active switching between photon nanojets and PH can be achieved in a variety of structures by changing the polarization state of the incident field. The PHs generation and modulation scheme proposed in this study can be extended to more anisotropic materials. Liquid crystal (LC) is the uniaxial anisotropic material with the lowest symmetry in optics and was used to produce Airy family beams [53,54,55]. Introducing LC dispersion into a polymer or fluid medium results in the spontaneous creation of droplets on a micrometer scale. These droplets exhibit a distinct molecular arrangement [56,57,58] in tangential and radial directions, which express significant advantages in the optical field. It has been verified that the tunable PNJ [59,60,61] is produced by changing the refractive index distribution of LC materials by applying an electric or light field. A conclusion is given that the horizontal and vertical displacements of PNJ depend strongly on the director of the LC. In contrast, LC materials have not been introduced into PH studies.
In this paper, we aim to study the critical properties of the PH for an optical tweeter using the scattering of LC to structured beams and further explore them [62]. A model is proposed using the scattering of a cutting liquid crystal cylinder to beams. The PHs generated by plane wave, Gaussian beam, and Bessel beam are discussed and compared. The characteristics of PHs, including FWHM, bending angle, and effective intensity, are analyzed. Furthermore, we conduct numerical investigations into the scattering characteristics of PH by incorporating the effects of birefringence. The results of this work provide a novel method to design PH, theoretical guidance, and a numerical model for the design of a new generation of optical tweezers.
The overall layout of this article is as follows. Section 2 presents the numerical model of a cutting LC cylinder illuminated by a beam based on the Finite-DifferenceTime-Domain (FDTD) method. The expressions of a plane wave, Gaussian beam, and Bessel beam are provided. The theoretical framework to generate the PHs has also been built. Based on this model, the PHs generated by the scattering of cutting the LC cylinder to plane wave, Gaussian beam, and Bessel beam are calculated in Section 3. The characteristics of PHs generated by various beams are compared. The impact of the beam and particle parameters on PHs is discussed in detail. The anisotropy refractive index of LC is emphasized. In the last section, the findings are summarized, and the potential application areas of this work are indicated. This work paves the way for developing novel optical devices with promising applications in optical manipulation, photoetch, and biomedicine.

2. Materials and Methods

Figure 1 presents the FDTD model that incident wave propagates along the x-axis and illuminates on cutting microcylinder [63] made from an anisotropic liquid crystal. In Figure 1, a perfectly matched layer with scattering boundaries fully absorbs outward waves and configures the polarization as transverse electric (TE). The simulated area has an area of 14 × 14   μ m 2 , the corresponding number of data grids is 700 × 700 . This means that the size of the grid is 0.02 μ m (about λ /32). Thus, the simulation results are sufficiently accurate. Assuming a circular microcylinder with a radius of 2 μ m, infinitely long, is cut at a distance of 1 μ m from the center and filled with LC. During practical experiments, the cutting microcylinder can be securely affixed to a gold base. Gold is chosen for its excellent conductivity and reflectivity [64,65], enabling adequate reflection of incoming plane waves and reducing interference during the propagation process.
Under the theoretical framework of FDTD [66,67,68,69], the Maxwell curl equations for linear, isotropic, non-dispersive, and lossy materials are as follows:
H t = 1 μ × E 1 μ M source   +   σ H
E t = 1 ε × H 1 ε J source   +   σ E
where σ and σ represents the electrical conductivity and effective magnetic loss, respectively. ε and μ are the permittivity and permeability. The current density ( J ) and the effective magnetic current density ( M ) can be considered as independent sources of electric ( J source ) and magnetic ( M source ) field energy.
In the Cartesian coordinate system, the coupled scalar expressions of Equations (1) and (2) are expanded as follows.
H x t = 1 μ E y z E z y M source   x + σ H x
H y t = 1 μ E z x E x z M source   y + σ H y
H z t = 1 μ E x y E y x M source   z + σ H z
E x t = 1 ε H z y H y z J source   x + σ E x
E y t = 1 ε H x z H z x J source   y + σ E y
E z t = 1 ε H y x H x y J source   z + σ E z
Equations (3)–(8) are the partial differential equations that electromagnetic waves interact with general three-dimensional objects in the FDTD numerical algorithm.
As everyone knows, the analytic expression of the incident plane wave in the x direction is expressed as
E ( x ) = E 0 exp ( j k x )
E 0 is the electric field complex amplitude respect with polarization, and E 0 = 1 . k = 2 π n / λ is the wave number, where λ and n are the free space wavelengths and the ambient medium refractive index.
For the fundamental mode Gaussian beam, the analytic expression can be expressed as [66,67]:
E ( ρ , x ) = E 0 w 0 w ( x ) exp ρ 2 w ( x ) 2 exp j k ρ 2 2 R ( x ) exp j k x
where ρ = x 2 + y 2 and z represent the radial and axial distances from the observation point to the beam focus. w ( x ) = w 0 1 + x / X 0 2 is the beam waist radius where the amplitude of the electric field decreases along the beam axis to 1 / e of its value. Where, w 0 is the initial beam waist radius, and R ( x ) = X 0 x X 0 + X 0 x is the curvature radius of the beam wavefront at position x. X 0 = 1 2 k w 0 2 = π w 0 2 λ is the Rayleigh dimension.
For the zero-order Bessel beam, the analytic expression is expressed as follows [70]
E ( x ) = E 0 J 0 k ρ exp j k x
where ρ and φ are the radial and angular coordinates of the field point. k and k are the transverse and axial components of the wave vector k, respectively. x = ρ cos φ and y = ρ sin φ , z is the axial coordinate of the field point. J 0 is the zero-order Bessel function.
To the anisotropic materials, the permittivity ε ˜ ( x , y ) [58,71] is defined as follows
ε ˜ ( x , y ) = ε 0 + Δ ε sin 2 θ ( x , y ) Δ ε sin θ ( x , y ) cos θ ( x , y ) 0 Δ ε sin θ ( x , y ) cos θ ( x , y ) ε 0 + Δ ε cos 2 θ ( x , y ) 0 0 0 ε 0
In Equation (11), Δ ε = ε e ε o is the refractive index difference between the extraordinary and ordinary permittivities. Where n o and n e denote the refractive indices of the material along the ordinary and extraordinary directions, respectively. In addition, angle θ ( x , y ) is the principal axis alignment of the LC molecules concerning the y-axis. To simplify the analysis, the x y plane is set as the isotropic plane, and the specific principal axis direction is the z direction, as shown in Figure 1. Thus, the LC molecules inside the cylindrical structure are tangentially aligned parallel to the surface. This assumption facilitates the introduction of optical anisotropy and allows us to verify its influence on electromagnetic fields.
The scattering models of the plane wave, Gaussian beam, and Bessel beam illuminating an anisotropic particle are established, focusing our attention on the PHs generated by the cutting LC microcylinder. Figure 2 is the schematic diagram of the PH. The start point, inflection point, and end point define it. Locating them is the key to calculating the bending angle. The inflection point is the position exhibiting the highest electromagnetic field intensity. The start point corresponds to the maximum intensity on the particle surface. The position of the extent of optical intensity decay in the PH determines the end point. It is the location where the optical intensity decreases to I max / e , where e is the natural constant. Thus, the bending angle is determined by calculating the angular difference between the vectors extending from the inflection point to the end point and from the inflection point to the start point.

3. Results

In this section, we perform numerical calculations to determine the PH formed by illuminating cutting microcylinder particles with various beam types, including the plane wave, Gaussian beam, and Bessel beam. The research is based on this ideal condition: the liquid crystal is confined to an infinitely thin and rigid film. The numerical results discuss the impact of beam and particle parameters on PHs, focusing on the anisotropic particles. The characteristics of PH, including maximum intensity, bending angle, and FWHM, are examined in detail. We present explicit definitions for the terminology employed in the calculations. FWHM is defined as the specific threshold at which it signifies half of the maximum intensity of the PH. At the threshold level, a horizontal line intersects the intensity curve at two locations, and the horizontal separation between these points corresponds to the FWHM. The electric field intensity of the PH is defined as E 2 E 2 max ( E i n c ) 2 max ( E i n c ) 2 , so the unit is “arb. units” (“a.u.”). Moreover, the refractive index of the environment is set n = 1.33 [72,73], and the wavelength of incident light was 632.8 nm. The ordinary refractive index is n o = 1.5 , but the extraordinary refractive index n e changes from 1.4 to 2.0 . This paper uses TE polarized wave (only the x-direction of the electric field is considered) incident particles as an example to explore the scattering of cutting LC microcylinders to different beams. The research method of the TM case is similar to that of the TE case, but we do not elaborate on it in detail there.

3.1. Plane Wave

Before discussing the PH produced by the interaction between structured beams and cutting LC microcylinders, pay attention to the results generated by the plane wave, as shown in Figure 3. The first column is the electric field intensity a plane wave exerts on the cutting LC microcylinders with different n e . As n e changes, the shape and characteristics of PHs continue to change, including effective length, max intensity, FWHM, and bending angle. The effective length of PHs gradually decreases as n e increases. On the contrary, during this process, the maximum intensity of the PHs gradually increases. This phenomenon is because the particles converge on the plane wave more significantly as n e increases, resulting in greater intensity. However, it is precisely because more beams are concentrated near the particle surface that the effective length of the PH becomes shorter. Therefore, during the experiment, one must choose the suitable refractive index according to the purpose of the experiment to obtain the PH with a longer effective length or greater intensity.
Moreover, compared with the first column, it is found that with the change of n e , the bending direction of the PH also changes, which is meaningful for selective manipulation of particles. To further explore how the characteristic of PH changes with n e , columns 2 to 4 are calculated as auxiliary graphs to the first column. The second column is the three-dimensional figures of PH, where all the leading and side lobes of PH in the x y plane are shown. According to the three-dimensional figures, the energy flow of the electric field near the particle is demonstrated, and the formation process of the PHs is expressed. The third and fourth columns are the electric field intensity of PH in the y and x directions, with the inflection point as the reference point, respectively. With these graphs, we can intuitively observe the change of PH intensity along the axis as n e changes and further explore the influence of n e on the FWHM and bending angle of PHs. In the third column, n e affects the main lobe and the side lobes of PHs. The inflection point of the main lobe is away from the y = 0 axis as n e increases. This shows from the side that the bending degree of PH is gradually increasing. Besides, the FWHM has also become narrower during this process. Therefore, a conclusion is proposed that a PH with a narrower FMHM can be obtained by increasing n e . It has a crucial role in nanoparticle manipulation. Although we can preliminarily obtain the trend of the effective length of the photonic hook with the increase in n e according to the results of the first column, we still calculate the figures of the fourth column to verify this conclusion. Comparing the figures in the fourth column, it is found that the length of PH decreases as n e increases.
To further verify the correctness of conclusions drawn from Figure 3, the PH characteristics are quantitatively calculated in Table 1. Parameter settings are consistent with those shown in Figure 3. The Max Intensity in the second column corresponds to the intensity of the inflection point. Its values are 3.36 , 5.27 , 5.45 , 6.07 , 6.63 , and 7.80 when n e change from 1.4 to 1.9 . On the contrary, the FWHM of PHs are 1.26 λ , 0.66 λ , 0.47 λ , 0.53 λ , 0.43 λ , and 0.37 λ . The FWHM of PHs has narrowed by almost half in this process. This is consistent with our conclusion in Figure 3. The results in Figure 3 do not intuitively yield the bending angle as one of the most important parameters of the PH, so modulation of the bending angle is a big problem for PH optical tweezers. We can not intuitively obtain the bending angle from the results in Figure 3. We give the changing bending angle trend as n e in the fourth column of Table 1. The bending angle is negative from n e = 1.4 to n e = 1.5 . It means that the PHs bend towards the negative direction of the y plane. The bending angle is positive when n e changes from 1.6 to 1.9. Thus, the PHs bend towards the positive direction of the y plane. In this process, the bending angle increases gradually as n e increases. According to Table 1, it can be observed that the FWHM of the PH is less than 0.5 λ , which means it can break the diffraction limit. This indicates that the PH has a high application value.
To further describe the causes and characteristics of PHs generated by a cutting LC cylinder, the energy distribution is calculated in Figure 4. It can be seen that the energy is incident from the left side of the particle, through the convergence of the particles, and mainly distributed in the front of the particle in Figure 4a–d. At the front of the particle, the direction of energy flows from both sides, but the energy distribution on both sides is not uniform. More energy is concentrated under the y < 0 plane. This is the physical mechanism of PH formation determined by the shape of the cutting LC cylinder. Comparing Figure 4a–d, when n e = 1.4 and n e = 1.5, energy is concentrated towards the plane y > 0. When n e = 1.6 and n e = 1.7, energy is focused on the plane y < 0. This corresponds to the results in Figure 3. The direction of energy flow determines the bending direction of the PHs. Besides, the forward energy of the particle is closer to the particle as n e increases. The bending angle of PHs also becomes bigger. It also explains the reason for the narrowing of FWHM: with the rise of n e , the position of the inflection point moves forward correspondingly, and more energy is gathered, which further causes the narrowing of FWHM.
Next, let us reveal the continuous change rule of PH characteristics as n e changes. Figure 5a–c elaborates the variations of bending angle, inflection point intensity, and FWHM of PHs as n e increases, respectively. We are only concerned with the values here. Thus, we ignore the ± of the bending angle in the drawing. When n e < 1.55 , the bending angle is negative. Otherwise, the bending angle is positive. This is because with the increase in n e , the convergence of particles on the beam becomes more assertive, and the focus of PHs moves to the particle surface, resulting in the change of its bending characteristics. We chose the step size of n e to be 0.01 because, under this step size, we can preliminarily obtain the characteristic on PHs and exclude the mutative effect of exceptional n e on PHs as much as possible.
In Figure 5a, it is observed that the overall trend of the bending angle is upward in waves. The maximum value of 26.95° is achieved at n e = 1.94, while the minimum value of 7.35° is achieved at n e = 1.56. The bimodal phenomenon explains the reason for the wavy trend. As shown in Figure 3j,n,r,v, there are two peaks with similar intensity close together. As n e increases, the position of the inflection point (maximum peak) moves back and forth between these two peaks. Thus, the bending angle shows a wave-rising trend as n e increases. In summary, increasing the extraordinary refractive index results in the PH with larger bending angles. In the experiment, we can obtain a smooth trend by suppressing the forward peaks of adjacent peaks. Figure 5b calculates the changing trend of the inflection point electric field intensity as n e . The inflection point intensity gradually increases as n e increases. The minimum value and maximum value are 3.35 and 8.37, corresponding to n e = 1.4 and n e = 1.94. Thus, PHs with strong intensity can be obtained by increasing the extraordinary refractive index in the experiment. It has great potential applications in optical information storage and lithography. Figure 5c gives the impact of n e on the FWHM. It can be seen that the general tendency of FWHM to change with n e is gradually decreasing. However, there is a jump peak from n e = 1.64 to n e = 1.65. The bimodal phenomenon explains it, as shown in Figure 6. When n e is less than 1.64, the inflection point coincides with the forward peak of the adjacent peaks. In contrast, when n e is greater than 1.64, the inflection point coincides with the backward peak of the adjacent peaks. Since the FWHM of the backward peak is always greater than that of the forward peak, the FWHM mutates when n e goes from 1.64 to 1.65. Therefore, within the range of n e from 1.4 to 2, the intensity of the PH gradually increases, the FWHM gradually decreases, and the bending angle fluctuates, but the overall trend is upward. This is beneficial for generating PH with dimensions smaller than the diffraction limit and with practical value.

3.2. Gaussian Beam

Having discussed plane waves, we focus on the most common beam in real life, the Gaussian beams. This section discusses the influence n e exerts on the PH generated by the Gaussian beam. Similar to the analytical process of the plane wave case, the characteristics of PH caused by the Gaussian beam are emphasized, including the bending angle, the electric field intensity, and the FWHM. In the calculation results, a cutting LC microcylinder is illuminated by a Gaussian beam with waist radius w 0 = 2   μ m propagating along the x-axis. Figure 7 is the schematic of the PH that a Gaussian beam interacts with a cutting LC microcylinder. However, the scattering results of particles at different locations in the Gaussian beam are different. In this paper, the particles are placed behind the waist of the Gaussian beam to compare with the plane wave and Bessel beam incident model of cutting liquid crystal cylinder, in which the initial position of the incident light is x = 6 μ m. The scattering of Gaussian beams by particles at different beam positions has been extensively studied, which is not the focus of this paper, and we will build an experimental platform to explore it in the future.
Figure 8 calculates the PHs caused by the Gaussian beam. The first column is the electric field intensity with different n e . Like the plane wave case, the first column shows that a Gaussian beam interacts with the cutting LC microcylinder. A PH appears on the back side of the particle. The second column is the three-dimensional figures of PHs in the x y plane, which are further calculated to elucidate the characteristics of PHs. The third and fourth columns are the electric field intensity of PH in the y and x directions. Comparing with Figure 3 and Figure 8, the intensity of PHs in Figure 8 is less than in Figure 3 under the same n e . The intensity of PHs still increases as n e increases, and the effective length of PHs decreases as n e increases. However, more secondary peaks appear in the x direction when a Gaussian beam illuminates the cutting LC microcylinder compared with the plane wave case comparing with the fourth column of Figure 3 and Figure 8. The intensity of secondary peaks is closer to that of the main peak. This makes it more difficult to choose the inflection point and further causes difficulties in determining the bending angle. Besides, the FWHM of PHs in Figure 8 is wider than that in Figure 3. It is related to the waist radius of the Gaussian beam.
To further compare the characteristics of PH generated by the plane wave and Gaussian beam, Table 2 is given. It can be seen that the intensity of the inflection point is 2.37, 3.35, 3.64, 4.24, 4.36, and 4.92. The trend of FWHM changing with n e is the opposite to the trend of inflection point intensity. The FWHM is 1.36 λ , 0.82 λ , 0.75 λ , 0.65 λ , 0.53 λ , and 0.6 λ . Thus, a conclusion is provided that the incident wave source also determines the characteristics of PHs generated by the cutting LC microcylinder. Compared with the Gaussian beam, the plane wave acting on cutting LC microcylinders is more likely to produce PHs with higher intensity inflection points and narrower FWHM. Of course, compared with the plane wave, the Gaussian beam has more tuning regimes, which is also extremely important for generating tunable PHs. However, the trend of bending angle changing with n e cannot be obtained from Table 2. Thus, Figure 9 is calculated to explore the influence of n e on PHs.
The intensity distribution of the inflection point is shown in Figure 9b, where the minimum value of 2.36 is achieved at n e = 1.4, and the maximum value of 5.26 is achieved at n e = 2.0. It with the conclusions drawn from Figure 9 and Table 2; the intensity of the inflection point increases as n e increases. When it comes to the FWHM, its values decrease initially and then fluctuate around the diffraction limit as n e increases, as shown in Figure 9c. As shown in the third and fourth columns of Figure 8, as well as Table 2, the peak values of transverse light field intensity and axial light field intensity of the PH increase with the increase in n e . As shown in Figure 9c, the FWHM values decrease initially and then fluctuate around the diffraction limit as n e increases. The maximum value 1.35 λ is achieved at n e = 1.4, and the minimum value 0.41 λ is achieved at n e = 1.92. However, the bending angle fluctuates as n e increases. Its minimum value of 1.99° is achieved at n e = 1.41, and the maximum value of 20.72° is achieved at n e = 1.93. The reason for the fluctuation change is the same as the case of the plane wave above, which is caused by the bimodal phenomenon. A more general rule can be obtained in the experiment by suppressing the bimodal phenomenon. Here, we calculate the bending angle according to the definition of the inflection point in this paper.

3.3. Bessel Beam

After discussing the characteristics of PHs, such as how plane waves and Gaussian beams interact with the cutting LC microcylinder, let us focus on Bessel beams. This section discusses the scattering characteristics of the cutting LC microcylinder to the zero-order Bessel beam. The specific principle is shown in Figure 10, where the generated PH is on the right side of the anisotropic cutting microcylinder particle, and the incident beam on the left side is a Bessel beam.
To discuss the scattering characteristics that the Bessel beam exerts on a cutting LC microcylinder, the scattering intensity is calculated in Figure 11. The first to the fourth columns are the electric field intensity, the three-dimensional electric field intensity, and the electric field intensity in the y and x directions, respectively. The intensity of PHs generated by the interaction between the Bessel beam and an LC microcylinder gradually increases as n e increases in Figure 11. The general trend is consistent with plane waves and Gaussian beams producing PHs. However, the intensity of PHs generated by the Bessel beam is less than that produced by the plane wave and Gaussian beam comparing Figure 3, Figure 8 and Figure 11. In the second and fourth columns of Figure 11, the bimodal phenomenon still exists when a Bessel beam illuminates to a cutting LC microcylinder. Also, the adjacent peaks are more intense than Figure 3 and Figure 8. This creates some difficulties in determining the inflection point. Besides, the PH is formed by the forward scattering of particles and the interference of the self-healing Bessel beam. This is why the shape of PHs in the front of the particle is not a whole. Thus, the PH that appears in front of the particle bends twice.
To further elucidate the characteristics of the PHs formed by the Bessel beam acting on particles, Table 3 is given. As we observe in Figure 11, the electric field intensity at the inflection point gradually increases as n e increases. However, the influence of n e on FWHM cannot be judged in Table 3, although FWHM tends to narrow with the rise of n e . Regarding the influence of n e on the bending angle, there is no obvious rule in Table 3. Thus, Figure 12 is calculated to explore how n e affects the characteristics of PHs. Figure 12a–c are the changing trend of bending angle, electric field intensity of inflection point, and FWHM as n e changes. Because of the self-restoring property of the Bessel beam and the influence of bimodal phenomena, the changing trend of the bending angle is chaotic with the increase in n e , as shown in Figure 12a. We only found that the maximum value of 22.02° is achieved at n e = 1.60, and the minimum value of 4.77° is achieved at n e = 1.63. On the contrary, the electric field intensity of the inflection point and FWHM present a more obvious rule with the increase in n e . In Figure 12b, the intensity of the inflection point shows a more robust trend as n e increases, although there is some volatility between n e = 1.84 and n e = 1.90. In Figure 12c, as n e increases, the overall trend of FWHM is gradually narrowing. However, there are some abrupt peaks in the range 1.53 < n e < 1.59 and 1.86 < n e < 1.92. It results from the self-restoring effect of the Bessel beam and the bimodal phenomenon.

4. Conclusions

The FDTD method was employed in this study to compute the PH characteristics (maximum intensity, bending angle, and FWHM) of anisotropic cutting microcylinder particles under different structured beam illuminations (plane wave, Gaussian beam, zero-order Bessel beam). The results show that without changing the beam properties, the PH generated by incident light can be influenced by changing the n e of the particles. The paper analyzes the influence of n e on PH characteristics, which contributes to increasing the degree of freedom in PH modulation. We also observed that introducing birefringence can change the bending direction of PH. The results show that regardless of the beam type, as the n e increases, the maximum intensity of the PH also increases. The overall trend of the FWHM initially decreases and then enters a fluctuating state with no apparent regularity in the bending angle. The PH generated by a plane wave has a higher maximum intensity and smaller FWHM, which is beneficial for applications in super-resolution imaging. The PH generated by a Gaussian beam has a lower maximum intensity than the plane wave but higher than the Bessel beam, and the smallest FWHM falls between the two. However, its length is longer, which is beneficial for applications in optical manipulation. The PH generated by a Bessel beam has the minor maximum intensity and the largest FWHM. Still, it exceeds the diffraction limit and exhibits bending twice due to its self-recovery property. This paper provides a new way to modulate PH. Besides, it also provides insights into the applications of different beams in PH generation and offers novel models for designing PH optical tweezers. The results of this paper have promising application prospects in nanoparticle manipulation, optical detection, and biomedicine.

Author Contributions

Conceptualization, R.L., M.Z., I.V.M. and O.V.M.; methodology, R.L.; software, M.Z. and F.L.; validation, H.T. and M.A. (Muhammad Arfan); formal analysis, H.T. and R.L.; investigation, N.K. and M.A. (Muhammad Asif); data curation, R.Y.; writing—original draft preparation, H.T.; writing—review and editing, R.L., H.T., M.Z., F.L., I.V.M. and O.V.M.; visualization, R.Y.; supervision, R.L., I.V.M. and O.V.M.; funding acquisition, R.L. 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 (62371378, 62071359) and Natural Science Basic Research Program of Shaanxi (Program No. 2023-JC-QN-0074).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Acknowledgments

The research thanks to the partial support of the Tomsk Polytechnic University Development Program.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
FDTDFinite-Difference Time-Domain
PHphotonic hook
LCliquid crystal
FWHMfull width at half maximum
TEtransverse electric
PNJphotonic nanojet

References

  1. Ashkin, A.; Dziedzic, J.M.; Yamane, T. Optical trapping and manipulation of single cells using infrared laser beams. Nature 1987, 330, 769–771. [Google Scholar] [CrossRef] [PubMed]
  2. Moffitt, J.R.; Chemla, Y.R.; Smith, S.B.; Bustamante, C. Recent advances in optical tweezers. Annu. Rev. Biochem. 2008, 77, 205–228. [Google Scholar] [CrossRef] [PubMed]
  3. Bao, Y.; Yu, S.S.; Anderegg, L.; Chae, E.; Ketterle, W.; Ni, K.K.; Doyle, J.M. Dipolar spin-exchange and entanglement between molecules in an optical tweezer array. Science 2023, 382, 1138–1143. [Google Scholar] [CrossRef]
  4. Riccardi, M.; Martin, O.J. Electromagnetic forces and torques: From dielectrophoresis to optical tweezers. Chem. Rev. 2023, 123, 1680–1711. [Google Scholar] [CrossRef] [PubMed]
  5. Essiambre, R.J. Arthur ashkin: Father of the optical tweezers. Proc. Natl. Acad. Sci. USA 2021, 118, e2026827118. [Google Scholar] [CrossRef]
  6. Resnick, A. Use of optical tweezers for colloid science. J. Colloid Interface Sci. 2003, 262, 55–59. [Google Scholar] [CrossRef]
  7. Li, J.; Chen, Z.; Liu, Y.; Kollipara, P.S.; Feng, Y.; Zhang, Z.; Zheng, Y. Opto-refrigerative tweezers. Sci. Adv. 2021, 7, eabh1101. [Google Scholar] [CrossRef]
  8. Ðorđević, T.; Samutpraphoot, P.; Ocola, P.L.; Bernien, H.; Grinkemeyer, B.; Dimitrova, I.; Vuletić, V.; Lukin, M.D. Entanglement transport and a nanophotonic interface for atoms in optical tweezers. Science 2021, 373, 1511–1514. [Google Scholar] [CrossRef]
  9. Kim, M.E.; Chang, T.H.; Fields, B.M.; Chen, C.A.; Hung, C.L. Trapping single atoms on a nanophotonic circuit with configurable tweezer lattices. Nat. Commun. 2019, 10, 1647. [Google Scholar] [CrossRef]
  10. Gordon, R. Biosensing with nanoaperture optical tweezers. Opt. Laser Technol. 2019, 109, 328–335. [Google Scholar] [CrossRef]
  11. LaFratta, C.N. Optical tweezers for medical diagnostics. Anal. Bioanal. Chem. 2013, 405, 5671–5677. [Google Scholar] [CrossRef]
  12. Choi, W.J.; Park, K.S.; Eom, T.J.; Oh, M.K.; Lee, B.H. Tomographic imaging of a suspending single live cell using optical tweezer-combined full-field optical coherence tomography. Opt. Lett. 2012, 37, 2784–2786. [Google Scholar] [CrossRef]
  13. Liang, S.; Cao, Y.; Dai, Y.; Wang, F.; Bai, X.; Song, B.; Zhang, C.; Gan, C.; Arai, F.; Feng, L. A versatile optoelectronic tweezer system for micro-objects manipulation: Transportation, patterning, sorting, rotating and storage. Micromachines 2021, 12, 271. [Google Scholar] [CrossRef]
  14. Ren, Y.; Chen, Q.; He, M.; Zhang, X.; Qi, H.; Yan, Y. Plasmonic optical tweezers for particle manipulation: Principles, methods, and applications. ACS Nano 2021, 15, 6105–6128. [Google Scholar] [CrossRef]
  15. Li, H.; Ren, Y.; Li, Y.; He, M.; Gao, B.; Qi, H. Nanoparticle manipulation using plasmonic optical tweezers based on particle sizes and refractive indices. Opt. Express 2022, 30, 34092–34105. [Google Scholar] [CrossRef]
  16. Chen, Z.; Taflove, A.; Backman, V. Photonic nanojet enhancement of backscattering of light by nanoparticles: A potential novel visible-light ultramicroscopy technique. Opt. Express 2004, 12, 1214–1220. [Google Scholar] [CrossRef] [PubMed]
  17. Yang, H.; Trouillon, R.; Huszka, G.; Gijs, M.A. Super-resolution imaging of a dielectric microsphere is governed by the waist of its photonic nanojet. Nano Lett. 2016, 16, 4862–4870. [Google Scholar] [CrossRef]
  18. Wen, Y.; Yu, H.; Zhao, W.; Wang, F.; Wang, X.; Liu, L.; Li, W.J. Photonic nanojet sub-diffraction nano-fabrication with in situ super-resolution imaging. IEEE Trans. Nanotechnol. 2019, 18, 226–233. [Google Scholar] [CrossRef]
  19. Chen, W.Y.; Liu, Y.Y.; Kong, J.A.N.; Li, L.P.H.; Chen, Y.B.; Cheng, C.H.; Liu, C.Y. Biological cell trapping and manipulation of a photonic nanojet by a specific microcone-shaped optical fiber tip. Opt. Lett. 2023, 48, 1216–1219. [Google Scholar] [CrossRef] [PubMed]
  20. Li, Y.C.; Xin, H.B.; Lei, H.X.; Liu, L.L.; Li, Y.Z.; Zhang, Y.; Li, B.J. Manipulation and detection of single nanoparticles and biomolecules by a photonic nanojet. Light. Sci. Appl. 2016, 5, e16176. [Google Scholar] [CrossRef] [PubMed]
  21. Siviloglou, G.; Broky, J.; Dogariu, A.; Christodoulides, D. Observation of accelerating Airy beams. Phys. Rev. Lett. 2007, 99, 213901. [Google Scholar] [CrossRef] [PubMed]
  22. Efremidis, N.K.; Chen, Z.; Segev, M.; Christodoulides, D.N. Airy beams and accelerating waves: An overview of recent advances. Optica 2019, 6, 686–701. [Google Scholar] [CrossRef]
  23. Minin, I.V.; Minin, O.V. Photonic Jets Formation by Non Spherical Axially and Spatially Asymmetric 3D Dielectric Particles. In Diffractive Optics and Nanophotonics: Resolution Below the Diffraction Limit; Springer: Berlin/Heidelberg, Germany, 2016; pp. 31–54. [Google Scholar]
  24. Minin, O.V.; Minin, I.V. The Photonic Hook: From Optics to Acoustics and Plasmonics; Springer: Cham, Switzerland, 2021. [Google Scholar]
  25. Ang, A.; Karochevsky, A.; Minin, I.; Minin, O.; Sukhov, S.; Shalin, A. Photonic Hook’based optomechanical nanoparticle manipulator. Sci. Rep. 2018, 8, 2029. [Google Scholar] [CrossRef] [PubMed]
  26. Yue, L.; Minin, O.V.; Monks, J.N.; Shalin, A.S.; Minin, I.V. Photonic hook: A new curved light beam. Opt. Lett. 2018, 43, 771–774. [Google Scholar] [CrossRef] [PubMed]
  27. Minin, I.V.; Minin, O.V.; Katyba, G.M.; Chernomyrdin, N.V.; Kurlov, V.N.; Zaytsev, K.I.; Yue, L.; Christodoulides, D. Experimental observation of a photonic hook. Appl. Phys. Lett. 2019, 114, 031105. [Google Scholar] [CrossRef]
  28. 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] [PubMed]
  29. Minin, I.V.; Minin, O.V.; Glinskiy, I.A.; Khabibullin, R.A.; Malureanu, R.; Lavrinenko, A.; Yakubovsky, D.; Volkov, V.; Ponomarev, D. Experimental verification of a plasmonic hook in a dielectric Janus particle. Appl. Phys. Lett. 2021, 118, 131107. [Google Scholar] [CrossRef]
  30. Yang, Y.J.; Zhang, D.L.; Hua, P.R. Array of photonic hooks generated by multi-dielectric structure. Opt. Laser Technol. 2023, 157, 108673. [Google Scholar] [CrossRef]
  31. Tang, H.; Li, R.; Gong, S.; Yang, L.; Yang, L.; Wei, B.; Zhu, Z.; Mitri, F. Curved photonic nanojet generated by a rotating cylinder. Opt. Express 2023, 31, 986–996. [Google Scholar] [CrossRef]
  32. Tang, H.; Shi, Z.; Zhang, Y.; Li, R.; Wei, B.; Gong, S.; Minin, I.V.; Minin, O.V. Unusual optical phenomena inside and near a rotating sphere: The photonic hook and resonance. Opt. Express 2024, 32, 12950–12966. [Google Scholar] [CrossRef]
  33. Baumgartl, J.; Mazilu, M.; Dholakia, K. Optically mediated particle clearing using Airy wavepackets. Nat. Photonics 2008, 2, 675–678. [Google Scholar] [CrossRef]
  34. Minin, I.V.; Minin, O.V.; Li, R.; Tang, H. BaTiO3 microsphere assisted imaging of erythrocytes by light structuring photonic hook. In Proceedings of the International Academic Conference on Optics and Photonics (IACOP 2023), Wuhan, China, 29–31 October 2024; SPIE: Bellingham, WA, USA, 2024; Volume 12972, pp. 73–77. [Google Scholar]
  35. Shang, Q.; Tang, F.; Yu, L.; Oubaha, H.; Caina, D.; Yang, S.; Melinte, S.; Zuo, C.; Wang, Z.; Ye, R. Super-resolution imaging with patchy microspheres. Photonics 2021, 8, 513. [Google Scholar] [CrossRef]
  36. Zou, P.; Xu, C.; Zhi, L.; Melinte, S.; Wang, Z.; Zuo, C.; Ye, R. Photonic hook-assisted contrast-enhanced super-resolution imaging using Janus microspheres. IEEE Photonics Technol. Lett. 2024, 36, 353–356. [Google Scholar] [CrossRef]
  37. Wu, Z.; Han, Y.; Wang, J.; Cui, Z. Generation of Bessel beam sources in FDTD. Opt. Express 2018, 26, 28727–28737. [Google Scholar] [CrossRef] [PubMed]
  38. Mandal, A.; Dantham, V.R. Photonic nanojets generated by single microspheres of various sizes illuminated by resonant and non-resonant focused Gaussian beams of different waists. JOSA B 2020, 37, 977–986. [Google Scholar] [CrossRef]
  39. Tang, H.; Shi, Z.; Zhang, Y.; Li, R.; Wei, B.; Gong, S.; He, X.; Yang, L.; Yan, B.; Sun, H.; et al. Scattering of a radially polarized Bessel beam by a PEMC sphere: Photonic nanojet and bottle beam formation. Appl. Opt. 2023, 62, 9164–9174. [Google Scholar] [CrossRef]
  40. Minin, O.V.; Minin, I.V.; Cao, Y. Time domain self-bending photonic hook beam based on freezing water droplet. Sci. Rep. 2023, 13, 7732. [Google Scholar] [CrossRef]
  41. Arfan, M.; Khaleel, N.; Asif, M. Interaction of Twisted Light Wave with Perfect Electromagnetic Conductor (PEMC) Sphere. Iran. J. Sci. 2023, 47, 1899–1908. [Google Scholar] [CrossRef]
  42. Arfan, M.; Asif, M.; Al-Kenani, A.N. Electromagnetic scattering on a dielectric sphere by vortex Gaussian beam field. Alex. Eng. J. 2024, 90, 112–119. [Google Scholar] [CrossRef]
  43. Arfan, M.; Asif, M.; Alkhoori, H.M. Orbital angular momentum based scattering characteristics for foggy atmosphere. Opt. Quantum Electron. 2024, 56, 630. [Google Scholar] [CrossRef]
  44. Jia, Y.; Li, R.; Zhuang, W.; Liang, J. Photonic nanojet generated by a spheroidal particle illuminated by a vector Bessel beam. Results Opt. 2021, 4, 100106. [Google Scholar] [CrossRef]
  45. Cheng, Y.; Shen, J. Curved photonic jet produced from a spherical particle illuminated by a Bessel–Gaussian beam. J. Quant. Spectrosc. Radiat. Transf. 2021, 272, 107765. [Google Scholar] [CrossRef]
  46. Spector, M.; Ang, A.S.; Minin, O.V.; Minin, I.V.; Karabchevsky, A. Temperature mediated ‘photonic hook’nanoparticle manipulator with pulsed illumination. Nanoscale Adv. 2020, 2, 2595–2601. [Google Scholar] [CrossRef] [PubMed]
  47. Chen, B.; Wei, K.; Cheng, Y.; Su, N.; Xu, Y.; Wu, P. Tunable twin photonic hooks generated by a double-layer fan-shaped microcylinder. Opt. Commun. 2024, 550, 129963. [Google Scholar] [CrossRef]
  48. Minin, I.V.; Minin, O.V.; Liu, C.Y.; Wei, H.D.; Geints, Y.E. Experimental demonstration of a tunable photonic hook by a partially illuminated dielectric microcylinder. Opt. Lett. 2020, 45, 4899–4902. [Google Scholar] [CrossRef]
  49. Zhang, J.; Han, G.; Yang, Z.; Xie, S.; Zhan, K. Photonic hooks generated by a concave micro-cylinder based on structure-constrained functions. Micromachines 2022, 13, 1434. [Google Scholar] [CrossRef]
  50. Lee, K.J.; Yoon, J.; Lahann, J. Recent advances with anisotropic particles. Curr. Opin. Colloid Interface Sci. 2011, 16, 195–202. [Google Scholar] [CrossRef]
  51. Voth, G.A.; Soldati, A. Anisotropic particles in turbulence. Annu. Rev. Fluid Mech. 2017, 49, 249–276. [Google Scholar] [CrossRef]
  52. Li, S.; Xu, J.; Yao, H.; Cheng, H.; Pang, T.; Wang, J.; Chi, T.; Zhang, B.; Lu, Y.; Liu, N. Controllable photonic hooks generated by typical shapes composed of uniaxial crystals. JOSA B 2024, 41, 659–664. [Google Scholar] [CrossRef]
  53. Luo, D.; Dai, H.; Sun, X.; Demir, H. Electrically switchable finite energy Airy beams generated by a liquid crystal cell with patterned electrode. Opt. Commun. 2010, 283, 3846–3849. [Google Scholar] [CrossRef]
  54. Wei, B.Y.; Chen, P.; Hu, W.; Ji, W.; Zheng, L.Y.; Ge, S.J.; Ming, Y.; Chigrinov, V.; Lu, Y.Q. Polarization-controllable Airy beams generated via a photoaligned director-variant liquid crystal mask. Sci. Rep. 2015, 5, 17484. [Google Scholar] [CrossRef]
  55. Wei, B.Y.; Liu, S.; Chen, P.; Qi, S.X.; Zhang, Y.; Hu, W.; Lu, Y.Q.; Zhao, J.L. Vortex Airy beams directly generated via liquid crystal q-Airy-plates. Appl. Phys. Lett. 2018, 112, 121101. [Google Scholar] [CrossRef]
  56. Lavrentovich, O.D. Topological defects in dispersed words and worlds around liquid crystals, or liquid crystal drops. Liq. Cryst. 1998, 24, 117–126. [Google Scholar] [CrossRef]
  57. Humar, M.; Ravnik, M.; Pajk, S.; Muševič, I. Electrically tunable liquid crystal optical microresonators. Nat. Photonics 2009, 3, 595–600. [Google Scholar] [CrossRef]
  58. Humar, M.; Muševič, I. 3D microlasers from self-assembled cholesteric liquid-crystal microdroplets. Opt. Express 2010, 18, 26995–27003. [Google Scholar] [CrossRef] [PubMed]
  59. Eti, N.; Giden, I.; Hayran, Z.; Rezaei, B.; Kurt, H. Manipulation of photonic nanojet using liquid crystals for elliptical and circular core-shell variations. J. Mod. Opt. 2017, 64, 1566–1577. [Google Scholar] [CrossRef]
  60. Liu, C.Y. Tunable photonic nanojet achieved using a core–shell microcylinder with nematic liquid crystal. J. Mod. Opt. 2013, 60, 538–543. [Google Scholar] [CrossRef]
  61. Du, B.; Xia, J.; Wu, J.; Zhao, J.; Zhang, H. Switchable photonic nanojet by electro-switching nematic liquid crystals. Nanomaterials 2019, 9, 72. [Google Scholar] [CrossRef] [PubMed]
  62. Zhang, M.; Li, R.; Wei, B.; Gong, S.; Minin, I.V.; Minin, O.V. Photonic hook generated by a cut anisotropic liquid crystal microcylinder. In Proceedings of the 2023 IEEE XVI International Scientific and Technical Conference Actual Problems of Electronic Instrument Engineering (APEIE), Novosibirsk, Russia, 10–12 November 2023; IEEE: Piscataway, NJ, USA, 2023; pp. 1320–1325. [Google Scholar]
  63. Minin, I.V.; Minin, O.V.; Cao, Y.; Yan, B.; Wang, Z.; Luk’yanchuk, B. Photonic lenses with whispering gallery waves at Janus particles. Opto-Electron. Sci. 2022, 1, 210008. [Google Scholar] [CrossRef]
  64. Guo, Y.; Singh, N.M.; Das, C.M.; Ouyang, Q.; Kang, L.; Li, K.; Coquet, P.; Yong, K.T. Two-dimensional PtSe2 theoretically enhanced Goos-Hänchen shift sensitive plasmonic biosensors. Plasmonics 2020, 15, 1815–1826. [Google Scholar] [CrossRef]
  65. Tao, L.; Deng, S.; Gao, H.; Lv, H.; Wen, X.; Li, M. Experimental investigation of the dielectric constants of thin noble metallic films using a surface plasmon resonance sensor. Sensors 2020, 20, 1505. [Google Scholar] [CrossRef] [PubMed]
  66. Oskooi, A.F.; Roundy, D.; Ibanescu, M.; Bermel, P.; Joannopoulos, J.D.; Johnson, S.G. MEEP: A flexible free-software package for electromagnetic simulations by the FDTD method. Comput. Phys. Commun. 2010, 181, 687–702. [Google Scholar] [CrossRef]
  67. Hammond, A.M.; Oskooi, A.; Chen, M.; Lin, Z.; Johnson, S.G.; Ralph, S.E. High-performance hybrid time/frequency-domain topology optimization for large-scale photonics inverse design. Opt. Express 2022, 30, 4467–4491. [Google Scholar] [CrossRef] [PubMed]
  68. Taflove, A.; Hagness, S.C.; Piket-May, M. Computational electromagnetics: The finite-difference time-domain method. Electr. Eng. Handb. 2005, 3, 15. [Google Scholar]
  69. Taflove, A.; Oskooi, A.; Johnson, S.G. Advances in FDTD Computational Electrodynamics: Photonics and Nanotechnology; Artech House: Norwood, MA, USA, 2013. [Google Scholar]
  70. Yu, Y.Z. Scalar and Vector Analysis of Bessel Beams. J. Quanzhou Norm. Univ. Sci. 2007, 25, 20–24. [Google Scholar]
  71. Matsui, T.; Okajima, A. Finite-difference time-domain analysis of photonic nanojets from liquid-crystal-containing microcylinder. Jpn. J. Appl. Phys. 2013, 53, 01AE04. [Google Scholar] [CrossRef]
  72. Wright, W.; Sonek, G.; Berns, M. Radiation trapping forces on microspheres with optical tweezers. Appl. Phys. Lett. 1993, 63, 715–717. [Google Scholar] [CrossRef]
  73. Gauthier, R.C. Computation of the optical trapping force using an FDTD based technique. Opt. Express 2005, 13, 3707–3718. [Google Scholar] [CrossRef]
Figure 1. The FDTD model that waves illuminate an anisotropic cutting liquid crystal microcylinder in water.
Figure 1. The FDTD model that waves illuminate an anisotropic cutting liquid crystal microcylinder in water.
Photonics 11 00736 g001
Figure 2. Schematic diagram of the photonic hook. The start point, inflection point, and end point are colored green, purple, and red. The inflection point corresponds to the maximum electric field intensity I max .
Figure 2. Schematic diagram of the photonic hook. The start point, inflection point, and end point are colored green, purple, and red. The inflection point corresponds to the maximum electric field intensity I max .
Photonics 11 00736 g002
Figure 3. The PH is generated by an anisotropic cutting microcylinder with different extraordinary refractive indices irradiated with plane waves. The first and second columns are the distribution figures of electric field intensity. The third and fourth columns represent the field intensities in y and x directions, respectively.
Figure 3. The PH is generated by an anisotropic cutting microcylinder with different extraordinary refractive indices irradiated with plane waves. The first and second columns are the distribution figures of electric field intensity. The third and fourth columns represent the field intensities in y and x directions, respectively.
Photonics 11 00736 g003
Figure 4. Energy distribution of a plane wave incident on the cutting LC microcylinder.
Figure 4. Energy distribution of a plane wave incident on the cutting LC microcylinder.
Photonics 11 00736 g004
Figure 5. Upon incidence of a plane wave, the bending angle, electric field intensity, and FWHM exhibit variations in response to changes in the n e . (a) Bending angle (b) Electric field intensity (c) FWHM.
Figure 5. Upon incidence of a plane wave, the bending angle, electric field intensity, and FWHM exhibit variations in response to changes in the n e . (a) Bending angle (b) Electric field intensity (c) FWHM.
Photonics 11 00736 g005
Figure 6. The bimodal phenomenon when a beam illuminates the cutting LC microcylinder.
Figure 6. The bimodal phenomenon when a beam illuminates the cutting LC microcylinder.
Photonics 11 00736 g006
Figure 7. Schematic diagram of PH generation by Gaussian beam illumination on cutting microcylinder particles.
Figure 7. Schematic diagram of PH generation by Gaussian beam illumination on cutting microcylinder particles.
Photonics 11 00736 g007
Figure 8. The PH generated by anisotropic cutting microcylinder particles with different extraordinary refractive indices irradiated with a Gaussian beam. The first and second columns are the distribution figures of electric field intensity. The third and fourth columns represent the field intensities in y and x directions, respectively.
Figure 8. The PH generated by anisotropic cutting microcylinder particles with different extraordinary refractive indices irradiated with a Gaussian beam. The first and second columns are the distribution figures of electric field intensity. The third and fourth columns represent the field intensities in y and x directions, respectively.
Photonics 11 00736 g008
Figure 9. Upon incidence of a Gaussian beam, the bending angle, electric field intensity, and FWHM exhibit variations in response to changes in the n e . (a) Bending angle (b) Electric field intensity (c) FWHM.
Figure 9. Upon incidence of a Gaussian beam, the bending angle, electric field intensity, and FWHM exhibit variations in response to changes in the n e . (a) Bending angle (b) Electric field intensity (c) FWHM.
Photonics 11 00736 g009
Figure 10. Schematic diagram of PH generation by Bessel beam illumination on cutting microcylinder particles.
Figure 10. Schematic diagram of PH generation by Bessel beam illumination on cutting microcylinder particles.
Photonics 11 00736 g010
Figure 11. The PH is generated by anisotropic cutting microcylinder particles with different extraordinary refractive indices irradiated with a Bessel beam. The first and second columns are the distribution figures of electric field intensity. The third and fourth columns represent the field intensities in y and x directions, respectively.
Figure 11. The PH is generated by anisotropic cutting microcylinder particles with different extraordinary refractive indices irradiated with a Bessel beam. The first and second columns are the distribution figures of electric field intensity. The third and fourth columns represent the field intensities in y and x directions, respectively.
Photonics 11 00736 g011
Figure 12. When a Bessel beam is incident, the bending angle, electric field intensity, and FWHM change in response to variations in the n e . (a) Bending angle (b) Electric field intensity (c) FWHM.
Figure 12. When a Bessel beam is incident, the bending angle, electric field intensity, and FWHM change in response to variations in the n e . (a) Bending angle (b) Electric field intensity (c) FWHM.
Photonics 11 00736 g012
Table 1. The comparison of the parameters of the PH generated by cutting microcylinder particles with different extraordinary refractive indices using plane waves.
Table 1. The comparison of the parameters of the PH generated by cutting microcylinder particles with different extraordinary refractive indices using plane waves.
n e Max Intensity (a.u.)FWHMBending Angle
1.43.361.26 λ −11.86°
1.55.270.66 λ −20.91°
1.65.450.47 λ +20.86°
1.76.070.53 λ +18.91°
1.86.630.43 λ +21.55°
1.97.800.37 λ +25.46°
Table 2. The comparison of the parameters of the PH generated by cutting microcylinder particles with different extraordinary refractive indices using a Gaussian beam.
Table 2. The comparison of the parameters of the PH generated by cutting microcylinder particles with different extraordinary refractive indices using a Gaussian beam.
n e Max Intensity (a.u.)FWHMBending Angle
1.42.371.36 λ −7.27°
1.53.350.82 λ −18.11°
1.63.640.75 λ +8.11°
1.74.240.65 λ +15.54°
1.84.360.53 λ +18.16°
1.94.920.60 λ +13.90°
Table 3. The comparison of the parameters of the PH generated by cutting microcylinder particles with different extraordinary refractive indices using a zero-order Bessel beam.
Table 3. The comparison of the parameters of the PH generated by cutting microcylinder particles with different extraordinary refractive indices using a zero-order Bessel beam.
n e Max Intensity (a.u.)FWHMBending Angle
1.41.661.20 λ −18.35°
1.52.070.88 λ −12.51°
1.62.310.90 λ +22.02°
1.72.550.64 λ +11.55°
1.82.930.52 λ +16.19°
1.92.840.63 λ +13.26°
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, R.; Tang, H.; Zhang, M.; Liu, F.; Yang, R.; Khaleel, N.; Arfan, M.; Asif, M.; Minin, I.V.; Minin, O.V. Tunable Photonic Hook Design Based on Anisotropic Cutting Liquid Crystal Microcylinder. Photonics 2024, 11, 736. https://doi.org/10.3390/photonics11080736

AMA Style

Li R, Tang H, Zhang M, Liu F, Yang R, Khaleel N, Arfan M, Asif M, Minin IV, Minin OV. Tunable Photonic Hook Design Based on Anisotropic Cutting Liquid Crystal Microcylinder. Photonics. 2024; 11(8):736. https://doi.org/10.3390/photonics11080736

Chicago/Turabian Style

Li, Renxian, Huan Tang, Mingyu Zhang, Fengbei Liu, Ruiping Yang, Naila Khaleel, Muhammad Arfan, Muhammad Asif, Igor V. Minin, and Oleg V. Minin. 2024. "Tunable Photonic Hook Design Based on Anisotropic Cutting Liquid Crystal Microcylinder" Photonics 11, no. 8: 736. https://doi.org/10.3390/photonics11080736

APA Style

Li, R., Tang, H., Zhang, M., Liu, F., Yang, R., Khaleel, N., Arfan, M., Asif, M., Minin, I. V., & Minin, O. V. (2024). Tunable Photonic Hook Design Based on Anisotropic Cutting Liquid Crystal Microcylinder. Photonics, 11(8), 736. https://doi.org/10.3390/photonics11080736

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