Abstract
1. INTRODUCTION
At the beginning of the new millennium, a phenomenon of jet-like-structured field enhancement effects originating from light–particle interaction was first known to researchers in laser physics [1–3]. In 2004, this phenomenon was thoroughly studied using numerical analysis on dielectric microcylinders illuminated by a plane wave and formally referred to as the “photonic nanojet (PNJ)” [4]. PNJ is a non-resonant light beam with a subwavelength transverse beam width, high optical intensity, and wavelength-scale propagation distance with weak diffraction that emerges from the shadow-side surface of the irradiated cylindrical or spherical microparticles with dimensions comparable to or slightly larger than the radiation wavelengths [5]. Over the past 16 years, PNJs have been applied in many fields such as enhancement of backscattering, fluorescence, transmission and Raman scattering [6–9], optical forces [10], optical data storage [11], coupled resonator optical waveguides [12], nanopatterning, nanolithography, and nano-marking [13–15], super-resolution imaging [16,17], surface plasmon polaritons [18], scanning probes [19], all-optical switching [20], single nanoparticle detection [21,22], biomolecules/cells trapping and manipulation [23,24], Bloch surface waves [25], optofluidic microcavities [26], biosensors and micromotors [27], and backaction of microparticles [28]. The shape features and field distributions of PNJs, which are usually characterized by the parameters of maximum light intensity (), focal distance () or working distance (), transverse beam width (), and decay length, are dependent on the material attributes and geometrical morphologies of the cylindrical or spherical microparticles, the surrounding medium, and the nature (type, wavelength, polarization, etc.) of incident light [5,29]. There are three basic rules followed between the key parameters (, , and ) of PNJs and the dielectric properties (refractive index, dimension, and shape) of microparticles [4,6,30]. First, the larger refractive index of the microparticle produces PNJs of closer focal distance (or working distance ) and narrower transverse beam width . Second, on the contrary, the larger the curvature radius of the microparticles, the farther (or ) and the wider for the generated PNJs. Third, a tight focus, with small and short (or ), is formed by illuminating the two edge regions of the first light refraction interface (LRI), while a loose focus, with large and long (or ), is formed by illuminating the middle region of the first LRI. Research on PNJ must be based on a full understanding of these three rules. The former two rules, with clear and definite physical mechanisms, have received general attention and been extensively studied by a number of theoretical and experimental works [31,32].
But so far, a systematic discussion on the third rule is lacking. It has been either implicitly or explicitly involved in studies from the last few years. For example, by etching concentric ring structures into the first LRI of a spherical microparticle with an unexplained location, number, and depth, Wu
Ray optics combined with numerical calculation is an effective way to qualitatively and quantitatively analyze the formation mechanism and distribution characteristics of PNJs [29–31,45]. The light scattering of a microparticle with no dissipation can be characterized by its refractive index and size parameters, and approximated as a ray-tracing model upon Snell’s law [31,46]. In terms of Snell’s law, the speeds at which light converges and diverges are closely related to the slopes of the emergent rays, where the slopes have a dependent relationship with the incidence angles of incident light (or the positions of incident light intersecting with the first LRI of the microparticle) [29]. The slope changes of the rays, in the case of no microparticle, are zero for all the propagating photons, while the slope changes are a nonlinear functional curve of incidence angles in the case of a microparticle within the model. In other words, the slopes of the emergent rays for the case of containing a microparticle can be regarded as positive and negative growth of slope values with respect to the case of no slope changes. The speed of light convergence and divergence can be intuitively understood as the length of time that the slopes of the outgoing rays change. This is exactly what natural logarithm means in mathematics: the natural logarithm reveals the time required to reach a certain level of growth [47,48]. In differential calculus and differential geometry, if the curvature of a smooth plane curve changes from being concave to convex, or vice versa, there must be a “boundary point” between the segments of the concave curve and convex curve [49]. We know that mathematically, this point is called the inflection point, which can be solved by finding the second derivative of the obtained logarithmic function curve [50].
Sign up for Photonics Research TOC. Get the latest issue of Photonics Research delivered right to you!Sign up now
In this paper, we present in-depth research on the third rule of PNJs produced by plane-wave-illuminated circular microparticles, i.e., the research on precisely dividing the first LRI of the microparticle into three parts: region of slow change (RSC), region of rapid change (RRC), and region of weak contribution (RWC). The key finding of this research is that the boundary point between RSC and RRC, namely, the inflection point in the text, can be figured out by thoroughly understanding the slope curve of emergent rays with the help of utilizing the mathematical meaning of a natural logarithm function and the second derivative of the function. An infinitely long microcylinder in three dimensions is used as the research model of a circular microparticle. Rigorous two-dimensional (2D) theoretical and numerical studies are carried out by ray optics analyses and finite element method (FEM) simulations. When RSC and RRC are separately irradiated by a plane wave, long PNJs with wide , long (), and short near-field focusing (including one central main lobe and two secondary sidelobes) with narrow , short (), are visually presented. As compared to the case of RWC adding to RRC, the induced relative changes of beam width are, respectively, 0.25% and 0.81% for the central main lobe and secondary sidelobes, while the values of () are the same for RRC and RRC plus RWC. The inflection point position changing with the refractive index contrast (RIC) of the microcylinder relative to its surrounding medium, the corresponding working distance, and relative value of the beam width at different incident angles, and whether the RWC exists in the conditions that the RIC is greater or less than the critical value () are analyzed and discussed with comprehensive ray optics theoretical analyses and full-wave simulation experiments. By illuminating RSC and the composite region (RRC + RWC) with light of different intensities, PNJ and PH are formed in the cases of the two composite regions on both sides of RSC irradiated with plane electromagnetic (EM) waves of equal or unequal intensity. The key parameters of PNJ and PH can be easily modulated by simply adjusting the light intensity ratios between RSC and the composite region for the cases of synchronously or asynchronously changing the light intensity of the two edge regions. The main points of this paper are: (1) finding and determining the inflection point between RRC and RSC; and (2) analyzing the light-focusing characteristics, presenting a new method to modulate PNJs, as well as to generate/modulate PHs, in the inflection point perspective. The results of this study could significantly promote the understanding, utilization, and application of PNJs to a new level.
2. RESULTS AND DISCUSSION
A. Research Approach and Determination of Inflection Point
The research model expressed in a three-dimensional (3D) schematic diagram is shown in Fig. 1(a). A monochromatic plane wave, with a free space wavelength of , unit intensity of , and wave vector propagating along axis is set as the incident light source. The electric and magnetic field components are on the axis and axis of the Cartesian coordinate, respectively. The microcylinder is assumed as a homogeneous, isotropic, and loss-free circular-column structure with infinite length in the direction. FEM and finite-difference time-domain are the two most widely used numerical computation methods for solving the internal and external field distributions of plane-wave-illuminated microparticles [45]; here, the FEM-based COMSOL Multiphysics software package is adopted to perform the 2D full-wave simulations in the later numerical experiments. To fully absorb the outgoing EM waves and eliminate undesired backreflections, a perfectly matched layer is chosen as the absorbing boundary condition on the two lateral sides and one end side of the simulation area. Refractive-index-dependent non-uniform meshes with fine element sizes of less than are used for the whole computational domain.
Figure 1.Schematic illustration of the research model and the obtained inflection point. (a) 3D schematic diagram of a dielectric microcylinder illuminated by an incident plane EM wave. (b) 2D representation of the research model in (a). Green solid lines show the process of a plane wave propagation through the microcylinder. RSC, region of slow change; RRC, region of rapid change; RWC, region of weak contribution;
Figure 1(b) is the schematic diagram of the research model showing a 2D sectional view. and are, respectively, the origin of coordinates and the radius of the microcylinder. The refractive indices of the microcylinder and the environmental medium are and , respectively. RSC and RRC are the two regions, as already mentioned in the Introduction, of light converging and diverging with either a slow or rapid mode. Ray-tracing models based on geometrical optics approximation show positive effects in qualitatively analyzing and calculating the beam variations of light propagating through small particles [46,51]. By way of ray-tracing, the regions of RSC and RRC are marked and roughly divided with the cyan-blue and magenta-pink transmission lines. The green solid lines represent a typical ray-tracing process following the law of refraction as it travels through the dielectric microcylinder. , and , are, respectively, the angles of incidence () and refraction () on both sides of the interfaces when light enters and leaves the microcylinder. Taking the axis of symmetry as the reference line, the range of is (0°, 90°) above the axis and (, 0°) below the axis. is the position of the intersection point between the circular curve of the microcylinder and the ray of emergent light. , , and are the slopes of the ray trajectories before, during, and after entering the microcylinder, respectively. Near the upper and lower vertices of the cylinder along the axis, we define a special region of RWC on each side. The reason that this region is called the weak contribution region is that the emergent light is emitted below the axis and propagates off the axis after the incident light from the upper edge refracted by the microcylinder. Hence, when the emergent light at point is treated as a secondary wave, the light from these points contributes very weakly to the focal spot with a focus on the axis.
The effects of weak contribution, tight focus, and loose focus on the illumination light in RWC, RRC, and RSC can be easily seen from the ray-tracing diagram in Fig. 1(c). The rays represented by the cyan-blue solid lines in RSC focus a little farther away from the microcylinder. In RRC, the region marked by magenta-pink solid lines, the convergent and divergent positions of the light rays are near the rear surface of the microcylinder. For the rays in RWC, the light-gray solid lines in the enlarged view of insertion, light travels as described in the previous definition. The light incident near the upper vertex of the microcylinder has its emitting point below the axis and propagates downward away from the axis. On the contrary, for the incident light near the lower vertex of the microcylinder, the emitting point is above the axis and the propagation direction is upward away from the axis. If only from the perspective of ray optics, the light coming from RWC has no contribution to the focal spot on the axis. However, as the interaction between light and microcylinders in actual conditions is based on wave optics theory, the secondary waves with positions deviating from the axis and propagation directions away from the axis still have a weak contribution to the light focusing of the focal spot on the axis.
According to Snell’s law, the slope of an emergent ray and the coordinate of point in Fig. 1(b) can be written as
Natural logarithms are often used in the mathematics and physics domain to solve the problems related to growth. Even the exact rate of growth is not known, and the cause of the change in growth can be found by taking the log of the function [47,48]. Figure 1(e) is the natural log of the positive part of the slope function in Fig. 1(d). The scope includes part of the RWC, the whole RRC, and the whole RSC with as the boundary. This is a continuous plane algebraic curve, whose outstanding feature is that the curve consists of an upward concave part (magenta-pink solid curve) and a downward convex part (cyan-blue solid curve). According to the concavity and convexity of curves in calculus, the value of the inflection point, marked with green dot in Fig. 1(e), can be obtained by finding the value of incidence angle () where the second derivative of function is equal to zero [49,50]. Figure 1(f) shows the value distributions of . This is a nonlinear curve that goes from positive to negative with increasing from to 0°. It is quite clear that the position of the inflection point is in the regions of RRC and RSC, which is also the position of the boundary point of RRC and RSC. As shown in the inset in the top panel in Fig. 1(f), the value of the incidence angle at the inflection point is .
B. Light-Focusing Analysis Based on Inflection Point Perspective
By using the position of the inflection point obtained, we can precisely delimit the regions of RRC and RSC on the first LRI and execute researches based on this delimitation. Figure 2 shows the comparison of optical field distribution and key parameters of the light-focusing generated from RSC, RRC, and composite region (RRC + RWC) exposed to light alone under the conditions of , , , and . From the previous calculations, the ranges of RSC, RRC, and RWC represented by the interval of incidence angles () are [−44.35°, 44.35°], [−82.82°, −44.35°) ∪ (44.35°, 82.82°], and [−90°, −82.82°) ∪ (82.82°, 90°], respectively. Referring to Refs. [29,30], the key parameters of focal distance (or working distance ) and beam width () are defined as the distance between the center of the microcylinder (vertex of the rear surface of the microcylinder) and the of the focal spot, and the transverse width of the focal spot using full width at half-maximum (FWHM), respectively.
Figure 2.2D optical field distributions of light illuminating the (a-1) RSC, (a-2) RRC, and (a-3) composite region of “RWC + RRC” obtained from FEM-based full-wave simulations. The model parameters are:
When the incident plane wave illuminates the RSC, as shown in Fig. 2(a-1), light focusing with a profile conforming to the conventional PNJ characteristics is formed on the shadow-side surface of the microcylinder. The focused light field contains only one central main lobe. From the normalized spatial distribution of PNJ in the range of to on the axis (cyan-blue solid line), the focal distance is , which is eight times greater than the illumination wavelength . The working distance is . The normalized transverse light field from , to , is drawn in Fig. 2(b-1). A beam width of is obtained by fitting the distribution curve with a Gaussian function. The corresponding distributions of Poynting vector fields and energy flow diagram represented by black conical arrows and purple solid curves in Fig. 2(c-1) indicate that the light beams within RSC propagate through the incoming and exiting interfaces with relatively small deflection angles, which results in the formation of a PNJ with long () and wide .
When the RRC is illuminated, as shown in Fig. 2(a-2), a distinct focused light field with and is formed close to the rear surface of the microcylinder. Such tight focusing consists of a central main lobe and two secondary sidelobes. The spacings between the main lobe and secondary lobes are clear and discernible. As can be seen in Fig. 2(b-2), the peak intensity of the secondary lobe is slightly more than half the . The beam widths of the main lobe and left and right sidelobes (relative to the frame of reference axis) are and , respectively, by using multi-peaks Gaussian fitting, which are far less than the classical diffraction limit (200 nm). The larger deflections of Poynting vector arrows and energy flow streamlines can be clearly seen at the front and rear interfaces of the microcylinder in Fig. 2(c-2), which leads to the production of specially distributed near-field focusing with short () and narrow .
Figures 2(a-3), 2(b-3), and 2(c-3) are, respectively, the light field distribution, transverse FWHM at the focal point along direction, and Poynting vector distributions and energy power flow for the case of the composite region composed by RWC and RRC illuminated by the plane wave. In comparison with the case where only the RRC is irradiated, the optical field profile and the focal distance (working distance ) are almost the same. The differences in the central main lobe and the secondary sidelobes are 0.25% and 0.81% for these two cases. It is rare to see any difference in field and streamline distributions between Fig. 2(c-3) and Fig. 2(c-2). The subtle difference in distribution is found only in the two secondary lobes (Fig. 9 in Appendix A). By comparing the positions of 11 streamlines on both sides of the axis for the cases of “RRC” and “RWC + RRC,” the deviations of the central main lobe within the effective width in the direction are zero, while the deviations of the secondary sidelobes within their effective widths are less than 1.47%. Obviously, the impact of RWC on light focusing is very weak, which is consistent with previous conclusions on ray optics analysis. So, we will use “composite-RRC” (C-RRC) for the region of “RWC RRC” in the remainder of this paper. As the spot intensity of PNJ in the case of RSC irradiated is higher than the spot intensity of near-field focusing in the case of C-RRC irradiated, the focus spot formed by the combined superposition effect, which is equivalent to the whole first LRI of the microcylinder irradiated, is the PNJ.
Figure 3.Inflection point between RSC and RRC of the first LRI for the model of plane-wave-illuminated microcylinder varies with the RIC of microcylinder and environment medium. (b) Position distributions of the emitting points
C. Inflection Point Positions and PNJ Characteristics Under Different RICs
It is known from the previous formulas and functions that the position of the inflection point is determined by the RIC of a dielectric microcylinder and environment medium. The independent variable adopted here is the angle of incidence, which can also be simply substituted by the coordinate of the intersection between the incident light ray and the first LRI of the microcylinder with formula of . Figure 3(a) shows the position of the inflection point along with the RIC increasing from 1.35 to 1.65 with step-size of 0.05. The size of the microcylinder is fixed at . As can be seen from the chart, the position of the inflection point moves towards the symmetry axis of the microcylinder in an approximately linear way as the RIC increases. The linear fitting function calculated by the calculation processing is . The reason for this trend is that the greater the RIC, the faster the deflection rate of the microcylinder to the incident light, and thus the earlier the transformation from RSC to RRC. It should be noted that the RWC does not exist for all RICs. There is a critical value of for RICs, where the grazing incidence at transmits to the right-most vertex of the microcylinder and emits [45]. When , e.g., as shown in the top left panel of Fig. 3(b), the emergent ray propagates outward away from the symmetry axis of the microcylinder with an emitting point above the axis in the case of the incident point () of the plane wave within the region below the blue horizontal line and above the yellow line, while the emitting points locate below the axis and propagate towards the axis for the vast majority of incident light above the blue line and below the axis. In the graph of with respect to , we can clearly see such a characteristic for the position distribution of the emitting points (blue curve). This is similar to the research model described above, in which the contribution to the focused light on the axis is weak. When RIC is equal to the critical value of [top center panel in Fig. 3(b)], with the exception of incident point right at the lower vertex position of the microcylinder, the emitting points of the emergent rays are below the axis (green curve in the figure) for all other incident light. When RIC decreases below , e.g., in the top right panel in Fig. 3(b), none of the incoming rays, after transmission through the microcylinder, will exit with the positions of the emitting points on or above the axis (red curve in the figure). All light emitted from points with negative coordinates and propagating towards the axis is a component of the focal spot on the axis. When the absence of RWC is understood as the value of RWC being null, we can still use C-RRC to represent “RWC RRC” in the following content.
Numerical simulation experiments are carried out to investigate the impact of different RICs on the key parameters of PNJ and near-field focusing. Figures 4(a) and 4(b) show the numerical results with the same parameters and variables as in Fig. 3(a). Because and have the same synchronized trends, we choose only for demonstration. Whether the light shines on the RSC or the C-RRC, the working distance and the transverse beam width of the central main lobe decrease with the increment in RIC. For the case of a plane wave illuminating the RSC, all the exceeds one wavelength, which can be considered as the focused light being the conventionally defined PNJs, and all the beam widths are above the diffraction limit. For the case of a plane wave illuminating the C-RRC, the characteristics and rules of the focused light produced are more complicated. The working distances of the central main lobe are within a wavelength range except (). When , the focal point returns to the surface of the microcylinder. The beam widths of both the central main lobe shown in Fig. 4(b) and the secondary sidelobes not shown here, are far below the diffraction limit. The minimum for , with beyond evanescent decay length () [30], is , which may offer promising prospects for nanoimaging, microprobes, micro/nano-manipulation, etc.
Figure 4.(a) Working distance
By virtue of a ray optics approach, the variation trend of and can be well analyzed and explained. The slope curve of the emergent light () and the position distribution curve of the emitting point () with respect to the coordinate of the incident point () are illustrated in Figs. 4(c) and 4(d). Actually, the two sets of curves in Fig. 4(c) are coming from a continuous set of curves in the range of to 0. For the purpose of making a distinct display of the investigated RSC and RRC, and showing the deviation of the slope curve along the vertical axis with different RIC values, we segment each constituent curve at the slope value corresponding to the inflection point of each RIC value. The position of the inflection point in Fig. 4(c) is the right-most point of each curve in the left set of curves, as well as the left-most point of each curve in the right set of curves. For a better comparison, we show only the part where the slope value is greater than zero. Based on the previous conclusion of critical value (), the set of curves on the left in Fig. 4(c), can cover all values ranging from zero to only for the cases of and . The slope values related to the incident points outside the dividing points for the curves corresponding to are not drawn. It can be seen from the slope curves, associated with either RSC or RRC, that the greater the RIC, the steeper the slope curves of the emergent rays. The steepness of the slope curve determines how fast the light converges and diverges at the rear surface of the microcylinder [29]. Furthermore, the width of the position distribution for both cases of RSC and RRC, as shown in Fig. 4(d), becomes smaller and smaller with the increase in RIC. Taken together, a steeper slope change of the emergent ray and a narrower position distribution of the emitting point produce a PNJ or a near-field focusing of shorter working distance and narrower transverse beam width ; conversely, a smoother slope change of the emergent ray and a wider position distribution of emitting points generate a light focusing of larger and .
D. Modulation of PNJ, Generation and Modulation of PH from Inflection Point
Based on the study of the inflection point discovered, we then propose and demonstrate a unique way for the modulation of light focusing generated from plane-wave-illuminated microparticles. In contrast to the traditional way of realizing modulation by using engineering techniques, such as index engineering (refractive index control [52], RIC variation [29], liquid crystal arrangement [53], material phase transition [54], etc.), geometry engineering (multilayer structure [55], non-spherical structure [56], engineered structure [35], etc.), and non-conventional designs [32], the modulation of light focusing by our approach can be achieved, without any changes in material attributes or geometrical morphologies of the microparticles, by setting and adjusting the different proportions of electric field intensity illuminating on RSC and RRC. In addition, we can also obtain the generation and modulation of PH in this way, with an approach different from the methods reported so far [38,57–59]. Figures 5(a-1) and 5(a-2) show two schematic diagrams of PNJ and PH existing on the shadow-side surface of the microcylinder for two different cases. Taking the inflection points on the first LRI as the cut-off points, the electric field intensities illuminating the C-RRC above, the RSC, and the C-RRC below are set as , , and , respectively. The subscripts u, m, and l denote the electric field applied to the upper, middle, and lower regions of the microcylinder, respectively. The size and refractive indices of the microcylinder and surrounding medium are the same as those in Fig. 2. Cases 1 and 2 represent the illumination conditions of , , and , where is set as the independent variable of the study. The electric field and the illumination wavelength for each region () are the same for cases 1 and 2. As the edge regions above and below are symmetric relative to the axis, the PHs in both the upper and lower directions indicated, respectively, by a red solid spot and light-blue dotted line in Fig. 5(a-2) will be formed when the electric field on one side is modulated.
Figure 5.Schematic diagrams illustrate the existence of (a-1) PNJ and (a-2) PH for the cases of
Figure 5(b) illustrates the modulation results of PNJ key parameters for case 1. The yellow and purple axes are the functions of and with respect to , respectively. During the modulation process, the intensities of and are always the same and change synchronously. The abscissa variable is from zero to ∞, where the leftmost “start” symbol for represents no incident plane wave illuminating the RSC, and the rightmost “stop” symbol for means being infinitely close to the case of the C-RRC above and C-RRC below irradiated by a plane wave. The interval between 0.1 and 500 is divided into four segments, where the step sizes set in the intervals of [0.1, 0.5], [1, 5], [10, 50], and [100, 500] are 0.05, 0.5, 5, and 50, respectively. With the increase in , both and show a nonlinear downward trend. The average rate of change (ARC), calculated by dividing the difference of by the difference of for the two endpoints in each interval, is used to characterize the trend of nonlinear change. The ARCs in the above four intervals, marked in Fig. 4(b), are , , , and for versus function and , , , and for versus function. It is observed that the decline of and probably goes through three stages of fast, slow, and steady. reaches its stable value faster than . The stable value of is almost exactly the same as the for the case of only the C-RRC illuminated by light. This means that there is no need to infinitely increase to modulate the key parameters of the light focusing. In fact, the narrowing of is less than 5% after ( for between 15 and 500). The insets in Fig. 5(b) show the contour plots of light field distributions for equaling 0.3, 1, 3, 30, and 300. When , the resulting PNJ has only one central main lobe. When , the secondary sidelobes appear and become more and more apparent. When is greater than a certain value, e.g., and , the contour distributions of the focusing constituted by the main lobe and secondary lobes change very little.
Case 2 describes the situation where and are not equal and are out of sync. Two examples of , and , are shown in Figs. 6(a) and 6(b), respectively. equals unit intensity. The optical field distributions displayed by heat-camera-color-filled contour plots draw the outline of two typical PHs with bending angles in different directions. The positions of the maximum light intensity (white dot marked in Fig. 6) are at , in Fig. 6(a) and , in Fig. 6(b). With reference to our previous work [58], two trajectory lines connected by several small linear segments, marked as Cambridge-blue and wine-red dotted lines in the two figures, show the optical field distribution along two curved trajectories. Figure 6(c) is a 3D display diagram showing the light intensity of the points on the two trajectories in the plane. The features of bending distribution are clearly visible as well. The normalized transverse optical field distributions at the point along direction for these two samples are shown in Fig. 6(d). The FWHMs obtained by Gaussian function fitting are and , respectively. Obviously, the higher the intensity of light illuminating the edge region of the microcylinder, the higher efficiency of the C-RRC in converging and diverging light, and hence the smaller the and of the produced PH. As can be seen from the results presented here, our method is different from the PH generation by common ways of breaking the symmetry of microparticles in geometrical structures [57] and material composition [58]. PHs generated by these two common ways are mainly from the interference of light with different phase velocities, and the differences in the phase velocity are induced by the different structures or material of a microparticle’s upper and lower parts. Our method proves that PHs can be produced only by setting the intensity of the light irradiating to the C-RRC above and C-RRC below to be different. It also shows that the PH phenomenon does not mean light traveling in curved lines; instead, the lights of different states interfere with each other to form a curved distribution of the light field.
Figure 6.PHs formed at the rear surface of the microcylinder under illumination conditions of (a)
Therefore, the PHs with different profiles can be easily created and effectively modulated by adjusting the values of and . Figure 7 shows the simulation results of working distance and transverse beam width depending on for the cases of and . As can be seen from the top left and bottom left panels in the figure, the value of decreases as increases, while the value of first decreases and then increases. The minimum is at . The reason that goes down in the first place is that as goes up, the edge part representing the light with a faster ability to converge and diverge plays a bigger and bigger role. Consequently, under the conditions that the focal spot intensity arising from the C-RRC illuminated not much more than the focal spot intensity from the RSC is illuminated, the combined effect of these two focusings causes the focal point to get closer and closer to the microcylinder, and the beam waist of the focal spot becomes smaller and smaller. gets bigger when reaches 2.5 and above, because the light focusing related to C-RRC starts to dominate. In the case of boundary illumination (equivalent to the C-RRC illuminated in this paper), the and the distance off the axis for each transverse optical field, in the region from the focal point to the surface of the microcylinder, increase as the distance between the transverse field and microcylinder surface decreases [44]. As goes from 2.75 to three, there is a significant change, that is, a jump in and . When the case of is viewed as the superposition of RSC irradiated and C-RRC above irradiated, the distribution of light intensity on the curved trajectory has two peaks, similar to the wine-red marked curve, 3D curve, and solid curve in Figs. 6(b), 6(c), and 6(d), respectively. Figure 10 in Appendix B clearly shows where these two peaks are located in the plane for , , , and , and the jump in , respectively. The and for the case of , as shown in the right panels in Fig. 7, experience a similar rule of change. Just because the C-RRC below is also illuminated by the plane wave, when the values of and become less and less affected by the increasing light exposure to C-RRC above, the focused field generated from plane-wave-illuminated C-RRC below will cancel out and decelerate this reduction. Thus, the and for each value in the case of are smaller than the case of . The jump also occurs when is larger. There is no further discussion after the jump, because the () is already less than the evanescent decay length ). Similarly, the value of () is less than this decay length at for case as well. And for in both cases, the value is always smaller, as the focus is always moving towards the surface of the microcylinder.
Figure 7.Modulation of PH key parameters (
Bending angle () is a particular characteristic parameter of PH distinct from PNJ [57,58]. At present, there are many definitions of bending angles, which are basically defined around the point, start point, and end point [39,40,42,57–61]. point is also often referred to as the “bend point” and “inflection point” [42,57]. Of course, the inflection point here is not the same as the inflection point discussed in this paper, but refers to the deflection of the PH beam around this point. In fact, the point at which the deflection occurs cannot be named inflection point, because there are several points where the PH beam deflects in some cases. We will give a detailed example in Appendix D. So, we named this point “ point” in this paper due to the uniqueness of point in PH. As for the start point and end point, there are always some problems such as fuzzy definition, unclear description, and no or unclear explanation. Through the reference and synthesis of these methods, we propose a method that can accurately define and calculate the bending angle. The method uses the points with light intensity of , and the farthest point on the contour curve, where each of these three points has a unique and definite position, as the start point, point, and end point. The detailed graphical presentation and text description are available in Appendix C. It is important to note that the bending angle is just one way of describing the curvature characteristic of PH. This is an artificial parameter, not an intrinsic parameter of PH. This is not to say that there can be only one definition of it. To quantify the bending angle, we give a new definition here to uniquely and certainly calculate the value of .
Figure 8 shows the variation of bending angle in both cases as shown in Fig. 7. As in Ref. [58], the angles of deflection in the clockwise and counterclockwise directions are set as negative and positive, respectively. From the previous analysis, the optical field of the main lobe is biased to the C-RRC side with loaded illumination light. The point is above the axis and point, and all the bending angles are negative for the case of . In the process of going from zero to 0.75, the bending angles change very quickly from 0° to the minimum value of . When , the value of is around , and the deviation between the maximum and minimum value is 6.11°. For , the is positive in the range of . In the negative part of , the change rate of is slower than the case of due to the cancellation effect induced by the light illumination of C-RRC below. Because of this cancellation effect, the point overall is closer to the central axis, and, in general, the for is less than the in the case of . The maximum difference between positive and negative bending angles is . The illustration shows two examples of PHs with for , , and for , . The dotted circle marked cases are the generated PNJs in the condition of .
Figure 8.Bending angles (
3. CONCLUSION
In this paper, starting from the third rule of light focusing produced by plane-wave-illuminated microparticles, the inflection point and the PNJ phenomenon on the inflection point perspective have been thoroughly discussed. The inflection point is a special point on the first LRI of the microparticle, which can accurately delimit and divide the regions where light converges and diverges with rapid or slow modes, i.e., the RRC or the RSC. The position of the inflection point can be determined by the sharpness of light focusing reflected by the slope curve of emergent rays, and the practical significance of a natural logarithm function and the second derivative of a function. Based on this partition and the other partition RWC, PNJ with only one central main lobe and near-field focusing with one central main lobe and two secondary sidelobes can be obtained by separately illuminating the RSC, RRC, and “RRC + RWC” of the microcylinder, respectively. The key parameters, i.e., the focal distance (working distance ) and the transverse beam width for RSC are long and wide, while for RRC and “RRC RWC,” they are short and narrow. The deviations of focal distance (working distance) and beam widths of the central lobe and sidelobes are, respectively, 0%, 0.25%, and 0.81% for the latter two cases. With the increase in RIC, the position of the inflection point approaches the central symmetry axis of the microcylinder in an approximately linear way. RWC exists only if RIC is less than the threshold of . The FEM-based simulation results and ray-optics-based theoretical analyses for RICs varying from 1.35 to 1.65 jointly verify that the variations on the working distance and the transverse beam width of the central main lobe follow a downward trend. In the case that and are equal and change simultaneously, the produced PNJ can be modulated, in a nonlinear way, by adjusting the field intensity of the illumination light on RSC and the composite RRC. The modulation ranges are over for and 106.44 nm to 247.27 nm for . In the case of and being out of sync, a new method for the generation of PH can be obtained. When tuning from small to large, the key parameters and bending angle () of the PH can be flexibly modulated. The modulation ranges of the characteristic parameter for ranging from zero to five in the cases of and are 24.09° and 40.1°, respectively. The modulation of PNJ and PH does not require any changes in the morphological structures or material properties of the microparticles. By a thorough study of the inflection point, our work theoretically solves the problem of dividing RSC and RRC on the first LRI of microparticles in the PNJ field, which will enable deeper insights into the nature of light–particle interaction.
APPENDIX A: COMPARISON OF ENERGY FLOW STREAMLINE DISTRIBUTION
To illustrate the subtle difference of the optical field distributions between the cases of plane wave illuminating the RRC and composite region (RWC RRC) of the first LRI, the distributions of energy flow streamlines are compared in the enlarged views. As shown in Figs.?
Figure 9.Distributions of energy flow streamlines of (a-1) RRC and (a-2) “RWC + RRC” on the first LRI illuminated by a plane wave. From the central symmetry axis of the microcylinder to both sides, the 11 streamlines are represented by Arabic numbers 1–11 and marked with small circles. (b) Relative differences in
APPENDIX B: JUMP IN WORKING DISTANCE Wd AND TRANSVERSE BEAM WIDTH ωc
The research results in Fig.?
Figure 10.Contour plots of optical field distributions for the illumination conditions of (a)
APPENDIX C: DEFINITION AND DESCRIPTION OF BENDING ANGLE δb
The key to defining and calculating bending angle is to determine the unique and definite positions of the start point, the point, and the end point. Figure?
Figure 11.(a) Schematic diagram used to define the bending angle of PH. “B,” selected arc line on the boundary between the microcylinder and environmental medium;
APPENDIX D: EXAMPLE OF A PH WITH MULTIPLE DEFLECTION POINTS
To illustrate multiple deflection points at a PH, the following example is established. As shown in Fig.?
Figure 12.Schematic diagram of a plane-wave-illuminated microcylinder made of composite materials. The refractive indices of the upper, middle, and lower parts of the microcylinder are
References
[1] B. S. Luk’yanchuk, Y. W. Zheng, Y. F. Lu. Laser cleaning of solid surface: optical resonance and near-field effects. Proc. SPIE, 4065, 576-587(2000).
[2] M. Mosbacher, H.-J. Münzer, J. Zimmermann, J. Solis, J. Boneberg, P. Leiderer. Optical field enhancement effects in laser-assisted particle removal. Appl. Phys. A, 72, 41-44(2001).
[3] B. Luk’anchuk. Laser Cleaning(2002).
[4] Z. Chen, A. Taflove, V. Backman. Photonic nanojet enhancement of backscattering of light by nanoparticles: a potential novel visible-light ultramicroscopy technique. Opt. Express, 12, 1214-1220(2004).
[5] A. Heifetz, S.-C. Kong, A. V. Sahakian, A. Taflove, V. Backman. Photonic nanojets. J. Comput. Theor. Nanosci., 6, 1979-1992(2009).
[6] X. Li, Z. Chen, A. Taflove, V. Backman. Optical analysis of nanoparticles via enhanced backscattering facilitated by 3-D photonic nanojets. Opt. Express, 13, 526-533(2005).
[7] Y. Li, X. Liu, X. Yang, H. Lei, Y. Zhang, B. Li. Enhancing upconversion fluorescence with a natural bio-microlens. ACS Nano, 11, 10672-10680(2017).
[8] H. S. Patel, P. K. Kushwaha, M. K. Swami, P. K. Gupta. Photonic nanojet assisted enhancement in transmission of light through hollow pyramid shaped near field probes. J. Opt., 17, 055005(2015).
[9] K. J. Yi, H. Wang, Y. F. Lu, Z. Y. Yang. Enhanced Raman scattering by self-assembled silica spherical microparticles. J. Appl. Phys., 101, 063528(2007).
[10] X. Cui, D. Erni, C. Hafner. Optical forces on metallic nanoparticles induced by a photonic nanojet. Opt. Express, 16, 13560-13568(2008).
[11] S.-C. Kong, A. Sahakian, A. Taflove, V. Backman. Photonic nanojet-enabled optical data storage. Opt. Express, 16, 13713-13719(2008).
[12] A. M. Kapitonov, V. N. Astratov. Observation of nanojet-induced modes with small propagation losses in chains of coupled spherical cavities. Opt. Lett., 32, 409-411(2007).
[13] E. McLeod, C. B. Arnold. Subwavelength direct-write nanopatterning using optically trapped microspheres. Nat. Nanotechnol., 3, 413-417(2008).
[14] J. Li, W. Gao, R. Dong, A. Pei, S. Sattayasamitsathit, J. Wang. Nanomotor lithography. Nat. Commun., 5, 5026(2014).
[15] B. Yan, L. Yue, J. N. Monks, X. Yang, D. Xiong, C. Jiang, Z. Wang. Superlensing plano-convex-microsphere (PCM) lens for direct laser nano marking and beyond. Opt. Lett., 45, 1168-1171(2020).
[16] Z. Wang, W. Guo, L. Li, B. Luk’yanchuk, A. Khan, Z. Liu, Z. Chen, M. Hong. Optical virtual imaging at 50 nm lateral resolution with a white-light nanoscope. Nat. Commun., 2, 218(2011).
[17] H. Yang, R. Trouillon, G. Huszka, M. A. M. Gijs. Super-resolution imaging of a dielectric microsphere is governed by the waist of its photonic nanojet. Nano Lett., 16, 4862-4870(2016).
[18] D. Ju, H. Pei, Y. Jiang, X. Sun. Controllable and enhanced nanojet effects excited by surface plasmon polariton. Appl. Phys. Lett., 102, 171109(2013).
[19] M. Duocastella, F. Tantussi, A. Haddadpour, R. P. Zaccaria, A. Jacassi, G. Veronis, A. Diaspro, F. DeAngelis. Combination of scanning probe technology with photonic nanojets. Sci. Rep., 7, 3474(2017).
[20] B. Born, J. Krupa, S. Gagnon, J. Holzman. Integration of photonic nanojets and semiconductor nanoparticles for enhanced all-optical switching. Nat. Commun., 6, 8097(2015).
[21] H. Yang, M. Cornaglia, M. A. M. Gijs. Photonic nanojet array for fast detection of single nanoparticles in a flow. Nano Lett., 15, 1730-1735(2015).
[22] G. Gu, J. Song, M. Chen, X. Peng, H. Liang, J. Qu. Single nanoparticle detection using a photonic nanojet. Nanoscale, 10, 14182-14189(2018).
[23] Y. C. Li, H. B. Xin, H. X. Lei, L. L. Liu, Y. Z. Li, Y. Zhang, B. J. Li. Manipulation and detection of single nanoparticles and biomolecules by a photonic nanojet. Light Sci. Appl., 5, e16176(2016).
[24] Y. Li, H. Xin, X. Liu, Y. Zhang, H. Lei, B. Li. Trapping and detection of nanoparticles and cells using a parallel photonic nanojet array. ACS Nano, 10, 5800-5808(2016).
[25] M. S. Kim, B. V. Lahijani, N. Descharmes, J. Straubel, F. Negredo, C. Rockstuhl, M. Häyrinen, M. Kuittinen, M. Roussey, H. P. Herzig. Subwavelength focusing of Bloch surface waves. ACS Photon., 4, 1477-1483(2017).
[26] X. Wu, Y. Wang, Q. Chen, Y.-C. Chen, X. Li, L. Tong, X. Fan. High-Q, low-mode-volume microsphere-integrated Fabry–Perot cavity for optofluidic lasing applications. Photon. Res., 7, 50-60(2019).
[27] Y. Li, X. Liu, X. Xu, H. Xin, Y. Zhang, B. Li. Red-blood-cell waveguide as a living biosensor and micromotor. Adv. Funct. Mater., 29, 1905568(2019).
[28] Y.-X. Ren, X. Zeng, L.-M. Zhou, C. Kong, H. Mao, C.-W. Qiu, K. K. Tsia, K. K. Y. Wong. Photonic nanojet mediated backaction of dielectric microparticles. ACS Photon., 7, 1483-1490(2020).
[29] G. Gu, R. Zhou, Z. Chen, H. Xu, G. Cai, Z. Cai, M. Hong. Super-long photonic nanojet generated from liquid-filled hollow microcylinder. Opt. Lett., 40, 625-628(2015).
[30] G. Gu, J. Song, H. Liang, M. Zhao, Y. Chen, J. Qu. Overstepping the upper refractive index limit to form ultra-narrow photonic nanojets. Sci. Rep., 7, 5635(2017).
[31] B. S. Luk’yanchunk, R. Paniagua-Dominguez, I. Minin, O. Minin, Z. Wang. Refractive index less than two: photonic nanojets yesterday, today and tomorrow [Invited]. Opt. Mater. Express, 7, 1820-1847(2017).
[32] J. Zhu, L. L. Goddard. All-dielectric concentration of electromagnetic fields at the nanoscale: the role of photonic nanojets. Nanoscale Adv., 1, 4615-4643(2019).
[33] M. X. Wu, B. J. Huang, R. Chen, Y. Yang, J. F. Wu, R. Ji, X. D. Chen, M. H. Hong. Modulation of photonic nanojets generated by microspheres decorated with concentric rings. Opt. Express, 23, 20096-20103(2015).
[34] B. Yan, L. Yue, Z. Wang. Engineering near-field focusing of a microsphere lens with pupil masks. Opt. Commun., 370, 140-144(2016).
[35] M. Wu, R. Chen, J. Soh, Y. Shen, L. Jiao, J. Wu, X. Chen, R. Ji, M. Hong. Super-focusing of center-covered engineered microsphere. Sci. Rep., 6, 31637(2016).
[36] M. Wu, R. Chen, J. Ling, Z. Chen, X. Chen, R. Ji, M. Hong. Creation of a longitudinally polarized photonic nanojet via an engineered microsphere. Opt. Lett., 42, 1444-1447(2017).
[37] L. Yue, B. Yan, J. Monks, R. Dhama, Z. Wang, O. Minin, I. Minin. Intensity-enhanced apodization effect on an axially illuminated circular-column particle-lens. Ann. Phys., 530, 1700384(2018).
[38] C. Liu, H. Chung, O. V. Minin, I. V. Minin. Shaping photonic hook via well-controlled illumination of finite-size graded-index micro-ellipsoid. J. Opt., 22, 085002(2020).
[39] I. V. Minin, O. V. Minin, C.-Y. Liu, H.-D. Wei, Y. E. Geints, A. Karabchevsky. Experimental demonstration of a tunable photonic hook by a partially illuminated dielectric microcylinder. Opt. Lett., 45, 4899-4902(2020).
[40] J. Yang, P. Twardowski, P. Gerard, Y. Duo, J. Fontaine, S. Lecler. Ultra-narrow photonic nanojets through a glass cuboid embedded in a dielectric cylinder. Opt. Express, 26, 3723-3731(2018).
[41] C.-Y. Liu, F.-C. Lin. Geometric effect on photonic nanojet generated by dielectric microcylinders with non-cylindrical cross-sections. Opt. Commun., 380, 287-296(2016).
[42] X. Shen, G. Gu, L. Shao, Z. Peng, J. Hu, S. Bandyopadhyay, Y. Liu, J. Jiang, M. Chen. Twin photonic hooks generated by twin-ellipse microcylinder. IEEE Photon. J., 12, 6500609(2020).
[43] H. Xing, W. Zhou, Y. Wu. Side-lobes-controlled photonic nanojet with a horizontal graded-index microcylinder. Opt. Lett., 43, 4292-4295(2018).
[44] F. Wang, L. Liu, P. Yu, Z. Liu, H. Yu, Y. Wang, W. J. Li. Three-dimensional super-resolution morphology by near-field assisted white-light interferometry. Sci. Rep., 6, 24703(2016).
[45] A. V. Itagi, W. A. Challener. Optics of photonic nanojets. J. Opt. Soc. Am. A, 22, 2847-2858(2005).
[46] A. S. Glassner. An Introduction to Ray Tracing(1989).
[47] M. J. Panik. Growth Curve Modeling: Theory and Applications(2014).
[48] A. B. Shiflet, G. W. Shiflet. Introduction to Computational Science: Modeling and Simulation for the Sciences(2014).
[49] J. Stewart. Calculus: Concepts and Contents(2009).
[50] I. N. Bronshtein, K. A. Semendyayev, G. Musiol, H. Muehlig. Handbook of Mathematics(2007).
[51] H. Guo, Y. Han, X. Weng, Y. Zhao, G. Sui, Y. Wang, S. Zhuang. Near-field focusing of the dielectric microsphere with wavelength scale radius. Opt. Express, 21, 2434-2443(2013).
[52] Y. E. Geints, A. A. Zemlyanov, E. K. Panina. Photonic nanojet calculations in layered radially inhomogennous micrometer-sized spherical particles. J. Opt. Soc. Am. B, 28, 1825-1830(2011).
[53] T. Matsui, K. Tsukuda. Direct imaging of tunable photonic nanojets from a self-assembled liquid crystal microdroplet. Opt. Lett., 42, 4663-4666(2017).
[54] J. Soh, M. Wu, G. Gu, L. Chen, M. Hong. Temperature-controlled photonic nanojet via VO2 coating. Appl. Opt., 55, 3751-3756(2016).
[55] Y. E. Geints, A. A. Zemlyanov, E. K. Panina. Photonic nanojet effect in multilayer micrometre-sized spherical particles. Quantum Electron., 41, 520-525(2011).
[56] L. V. Minin, O. V. Minin, Y. E. Geints. Localized EM and photonic jets from non-spherical and non-symmetrical dielectric mesoscale objects: brief review. Ann. Phys., 527, 491-497(2015).
[57] L. Yue, O. V. Minin, Z. Wang, J. N. Monks, A. S. Shalin, I. V. Minin. Photonic hook: a new curved light beam. Opt. Lett., 43, 771-774(2018).
[58] G. Gu, L. Shao, J. Song, J. Qu, K. Zheng, X. Shen, Z. Peng, J. Hu, X. Chen, M. Chen, Q. Wu. Photonic hooks from Janus microcylinders. Opt. Express, 27, 37771-37780(2019).
[59] C.-Y. Liu, H.-J. Chung, E. Hsuan-Pei. Reflective photonic hook achieved by dielectric-coated concave hemi-cylindrical mirror. J. Opt. Soc. Am. B, 37, 2528-2533(2020).
[60] I. V. Minin, O. V. Minin, G. M. Katyba, N. V. Chernomyrdin, V. N. Kurlov, K. I. Zaytsev, L. Yue, Z. Wang, D. N. Christodoulides. Experimental observation of a photonic hook. Appl. Phys. Lett., 114, 031105(2019).
[61] Y. Geints, I. V. Minin, O. V. Minin. Tailoring ‘photonic hook’ from Janus dielectric microbar. J. Opt., 22, 065606(2020).
[62] A. A. R. Neves. Photonic nanojets in optical tweezers. J. Quant. Spectrosc. Radiat. Transfer, 162, 122-132(2015).
[63] B. Zhang, J. Hao, Z. Shen, H. Wu, K. Zhu, J. Xu, J. Ding. Ultralong photonic nanojet formed by dielectric microtoroid structure. Appl. Opt., 57, 8331-8337(2018).
Set citation alerts for the article
Please enter your email address