- Photonics Research
- Vol. 11, Issue 8, 1437 (2023)
Abstract
1. INTRODUCTION
Nonlinear optical harmonic generation is of great significance in a broad range of technologies and has been attracting much attention in photonics, chemistry, and biosensing. The understanding and manipulation of nonlinear optical properties of nanostructures is a major issue for realizing nonlinear nanophotonic devices. In the past decades, a lot of efforts have been made to enhance the nonlinear frequency conversion efficiency, including the new materials and systems [1–7]. Recently, with the advent of plasmonics, metamaterials, and metasurfaces that go beyond natural materials in many aspects, significant attention has been devoted to the understanding and the observation of nonlinear optical processes in nanophotonics [1–7]. For boosting the second-harmonic generation (SHG) [8–13], as one of the most important nonlinear optical effects, many approaches have been proposed. For example, through exciting resonance modes like surface plasmon polariton resonances [8], Fano resonances [9], anapole modes [10,11], and bound states in the continuum [12], the fundamental fields inside the nonlinear materials can be greatly enhanced, so as to boost the SHG. Notably, epsilon-near-zero (ENZ) media [7,14–17] with a vanishing permittivity are found to exhibit pronounced nonlinear optical properties [7,18–38]. The ENZ media can provide large field enhancement due to the continuity of the normal component of the electric displacement field across the interface [18,19], and at the same time provide unique opportunities for realizing phase-matching conditions [23,38]. Consequently, the ENZ media can give rise to unprecedented strong nonlinear optical responses including the significantly enhanced SHG [18–22,27–29].
In this work, on the contrary, we demonstrate that the ENZ media can be exploited to “turn off” the SHG of a nonlinear particle. We show that when a tiny nonlinear particle is placed very close to a subwavelength ENZ particle, the SHG conversion efficiency can be reduced by more than 2 orders of magnitude as compared to that from the nonlinear particle alone (Fig. 1). The SHG quenching effect attributes to the extraordinary local evanescent fields occurring near the ENZ particle due to evanescent scattering waves. Remarkably, the prohibition of electric fields appears besides the ENZ particle, and therefore, the SHG from the tiny nonlinear particle placed there would be suppressed. We find that this extraordinary local evanescent field exists for both isotropic and anisotropic ENZ particles, irrespective of their shapes. Furthermore, based on the principle of SHG quenching effect, we propose a kind of dynamically controllable optical metasurface integrated with anisotropic ENZ media, which consist of alternative layers of semiconductor material cadmium oxide (CdO) and phase-change material germanium telluride (GeTe). We find that through changing the phase states of GeTe, the SHG from the metasurface can be switched on or off. Our work demonstrates a feasible approach for controlling nonlinear responses with ENZ media, which may find applications in optical switches and modulators.
Figure 1.(a) Schematic illustration of SHG from a tiny nonlinear particle. (b) The SHG is quenched when a linear ENZ particle is placed very close to the tiny nonlinear particle.
2. SHG QUENCHING AND THE UNDERLYING PHYSICS
First, to explore the physical mechanism of the unique SHG quenching effect, we begin with a simple configuration comprising a linear ENZ spherical particle (relative permittivity , relative permeability ) and a tiny nonlinear spherical particle (relative permittivity , relative permeability ) in free space, as illustrated in Fig. 2(a). Both particles are in the deep-subwavelength scale with radii of (ENZ particle) and (nonlinear particle), where is free-space wavelength. The gap (i.e., edge-to-edge distance) between the two particles is . The two particles are illuminated by a plane wave propagating along the direction with electric field polarized in the direction.
Sign up for Photonics Research TOC. Get the latest issue of Photonics Research delivered right to you!Sign up now
Figure 2.(a) Schematic layout of the configuration for exploring the SHG quenching effect. It is composed of a linear ENZ spherical particle accompanied with a tiny nonlinear spherical particle in the deep-subwavelength scale. The green curve denotes the trajectory of the nonlinear particle moving around the ENZ particle. Positions 1 and 4 are close to the poles, positions 2 and 3 are on the equatorial plane, and position 5 is far from the ENZ particle. The two particles are illuminated by a plane wave propagating along the
The SHG of the tiny nonlinear spherical particle can be accounted for by effective surface contribution and characterized by a surface nonlinear current source , where is the nonlinear susceptibility tensor, and is the fundamental frequency electric field inside the nonlinear particle. is the permittivity of free space. For isotropic and centrosymmetric materials, the second-order surface susceptibility tensor can be reduced to only three independent nonzero components , , and [39,40]. Here the symbols and represent the directions perpendicular and parallel to the particle’s surface, respectively.
Figure 2(b) shows the scattering cross section of SHG () from the tiny nonlinear particle when it is successively moved from position 1 to position 5 along the trajectory [green curves in Fig. 2(a)]. Along the trajectory , the edge-to-edge distance is kept unchanged at , and position 5 is at a distance of from the ENZ sphere’s center. Here we assume that the nonlinear susceptibility of the tiny particle is several orders of magnitude larger than that of the ENZ particle, and therefore, the nonlinear response of the ENZ particle is negligible. Such an assumption is reasonable when using appropriate optical materials [30], as we will show in the following practical implementations. The calculation is performed using the finite-element software COMSOL Multiphysics. The solid blue lines and red dots denote the normalized of the cases with and without the ENZ particle, respectively. The latter case is normalized to 1. As expected, in the absence of the ENZ particle, the will not change as varying the position under the illumination of a plane wave. Interestingly, we observe that the presence of the ENZ particle has a great influence on the SHG conversion efficiency of the nonlinear particle. The reaches the maximal value when the nonlinear particle is on the equatorial plane of the ENZ particle (positions 2 and 3), while it is largely reduced when the nonlinear particle is close to the poles of the ENZ particle (positions 1 and 4). Compared with the case without the ENZ particle, we find that the from the nonlinear particle at position 1 or 4 is reduced by more than 1 order of magnitude, indicating that the SHG is quenched. Such an SHG quenching effect would fade away if the nonlinear particle is at a distance comparable with the radius of the ENZ particle (position 5), indicating that the quenching of SHG is a short-range effect. The influence of the geometrical parameters on the SHG quenching effect is elaborated in Appendix A.1.
The underlying physics of such a unique position-dependent characteristic lies on the extraordinary evanescent fields near the ENZ particle. When encountering the ENZ particle, the incident light will be scattered into various directions. If the size of the ENZ particle is comparable to or smaller than the wavelength, strong scattered evanescent waves that are evanescent in the forward and backward directions but capable of transferring energy flux along the perpendicular directions would emerge [41–43]. Such evanescent fields vary dramatically with positions, thus leading to the position-dependent SHG efficiency, as observed in Fig. 2. Particularly, near-zero electric field occurs on two sides of the ENZ particle perpendicular to the direction of incidence and along the direction of polarization, which is also denoted as the “side scattering shadows” [41]. When the tiny nonlinear particle is placed there, the incident light cannot “see” it, and therefore the SHG is suppressed.
The generation of the dramatically varying evanescent fields stems from the continuity boundary condition at the ENZ particle–air interface. Since the ENZ particle is at the deep-subwavelength scale, the fundamental electric field inside it can be considered uniform and can be approximately expressed as [43]. Thus, the continuity of the electric displacement field at the poles of the ENZ particle yields
Equation (3) indicates the dramatically varying fundamental electric field around the ENZ particle and near-zero electric field at the poles. When a tiny nonlinear particle is placed close to the ENZ particle, it will experience significantly different fundamental electric fields at different positions, thus leading to the position-dependent SHG efficiency. Remarkably, the SHG is quenched when the nonlinear particle is located at the poles as the fundamental electric field there is shielded by the ENZ particle. We note that the extraordinary evanescent fields only exist in the regions very close to the ENZ particle, and therefore the quenching of SHG is a short-range effect.
More strict proof can be obtained based on the Mie scattering theory [44] and the electrostatic theory [45]. Under the limit of (here ), the fundamental electric field outside the ENZ sphere () can be derived as (see Appendix A.3 and A.4)
Equation (5) implies that the fundamental electric field is enhanced (i.e., ) on the equator (i.e., ), while is zero (i.e., ) at the poles (i.e., ), coincident with the above analysis.
For numerical verification, in Fig. 3(a) we show the simulated fundamental electric-field amplitude on the ENZ sphere’s surface on the air side in the absence of the nonlinear sphere. The numbers 1–5 represent five different positions corresponding to the positions of the tiny nonlinear sphere in Fig. 2. Apparently, the electric field on the equator (positions 2 and 3) is enhanced compared with that of the incidence, while it is largely reduced to near-zero at the poles (positions 1 and 4). At position 5 which is far away, the influence of the ENZ particle fades away, and the electric field there tends to that of incidence. We note that these extraordinary local electric fields are observed for relatively ideal ENZ particle with . In order to explore the influence of the on the local electric fields, in Fig. 3(b) we plot the normalized electric-field amplitude at position () as a function of . We observe that the electric field (or ) decreases (or increases) quickly as the decreases from unity and becomes stable as long as . Since the evanescent fields only exist very near the ENZ particle, the electric field remains almost unchanged.
Figure 3.(a) Distribution of normalized fundamental electric-field amplitude on the ENZ particle’s surface on the air side in the absence of the nonlinear particle. The numbers 1–5 denote five different positions of the nonlinear particle in Fig.
We note that the electric field is not that close to zero because the positions 1 and 4 have a distance of from the poles of the ENZ sphere. Actually, with reducing , the will tend to be , as implied in Eq. (1). On the other hand, the tends to be 1.5 according to Eq. (2). For verification, in Fig. 3(c) we plot the ratio with the varying ratio ( is kept unchanged at ). Here we set . The results clearly show that the ratio approaches the ideal value (i.e., , red dashed lines) when is considerably small, consistent with the above theory.
However, a practical nonlinear particle cannot be infinitely small. The presence of the nonlinear particle could affect the evanescent fields of the ENZ particle. Therefore, the nonlinear particle should be much smaller than the ENZ particle, such that the original evanescent fields are not largely disturbed, and the nonlinear particle can experience enhanced or weakened local fields instead of the averaged fields. In order to show the influence of the particle’s size, the ratio of SHG scattering cross section as a function of the ratio is plotted in Fig. 3(d). Here is kept unchanged at . The (or ) denotes the SHG scattering cross section when the nonlinear particle is placed at position 1 (or 2), as illustrated by the inset in Fig. 3(d). We see that the ratio increases with the reducing , indicating better SHG quenching performance with smaller .
The above results suggest that in order to achieve excellent performance of SHG quenching effect, the conditions of and shall be satisfied. In addition, it is noteworthy that the SHG quenching effect is robust against the material loss of the ENZ particle. We find that even when the possesses an imaginary part 3 orders of magnitude larger than its real part, the good performance of SHG quenching effect can still be obtained (see Appendix A.2).
3. SHG QUENCHING IN GENERAL SITUATIONS
The extraordinary evanescent waves and the consequent SHG quenching effect are general and could be observed in arbitrary-shaped ENZ particles and anisotropic ENZ particles. As an example, in Fig. 4(a) we present the simulated fundamental electric-field distribution on the surface of a finite-sized isotropic ENZ film (, with dimensions ) on the air side under illumination of a plane wave. Prohibition of electric field above (or below) the cuboid’s upper (or lower) surface is clearly observed. Under this circumstance, a tiny nonlinear particle placed above the ENZ film cannot “see” the incident light, and therefore the SHG will be suppressed. For a demonstration, in Fig. 4(b) we compare the normalized SHG scattering cross section from a nonlinear particle alone (red lines) and that from the nonlinear particle placed above the ENZ film at an edge-to-edge distance of (blue lines) as a function of . Here the nonlinear particle is the same as that in Fig. 2. It is seen that the can be reduced by more than 2 orders of magnitude by the ENZ film with compared with the from the nonlinear particle alone, indicating the suppression of SHG by the ENZ cuboid.
Figure 4.(a) and (c) Distribution of normalized fundamental electric-field amplitude on the surface of an (a) isotropic, (c) anisotropic ENZ film on the air side in the absence of the nonlinear particle. The incident light is polarized in the
Moreover, we consider an anisotropic ENZ film with , , and , which are the , , and components of the relative permittivity tensor, respectively. We see from Fig. 4(c) that the prohibition of electric field above (or below) the cuboid’s upper (or lower) surface still exists irrespective of the existence of anisotropy. This is because the component of the electric field on the air side is proportional to and is therefore near zero, as a result of the continuity of the electric displacement field. The SHG quenching effect is verified in Fig. 4(d). We see that the anisotropic ENZ film with can significantly reduce the SHG scattering cross section from the nonlinear particle to near 2 orders of magnitude (blue lines) compared with the from the nonlinear particle alone (red lines). These results clearly manifest the universality of the SHG quenching effect.
4. DYNAMIC CONTROL OF SHG QUENCHING WITH OPTICAL METASURFACES
The universality of the SHG quenching effect suggests that it could be extended to other systems besides the above particle-scattering system to explore more efficient ways to control the SHG conversion efficiency. In fact, metasurfaces [46–49], as two-dimensional artificial electromagnetic materials structured at the subwavelength scale, provide a versatile platform for manipulating light–matter interactions and realizing significant nonlinear responses [4–6,12,13]. In the following, we would like to show a kind of optical metasurface exhibiting switchable SHG quenching effect.
Figures 5(a) and 5(b) show the schematic graphs of the optical metasurface, consisting of a square array of meta-atoms with a lattice constant of 200 nm on a 2 μm thick silica () substrate. Each meta-atom comprises a large cuboid ( in dimension) accompanied with two small cuboids ( in dimension) on the upper and lower sides, as illustrated by the insets. The top surface of the cuboid is covered by a layer of graphene. In practical realization, such patterned graphene could be fabricated through photolithography [50], electron beam lithography [51], plasma etching [52], laser direct patterning [53], etc. The graphene is known as a unique nonlinear optical material with outstanding optical properties [13,54–56]. Generally, the SHG is forbidden of free-standing graphene due to its centrosymmetry. Interestingly, when the graphene is placed on a substrate with broken inversion symmetry (e.g., the utilized here), SHG would be allowed [13]. Considering the fact that graphene belongs to the symmetry group, its surface second-order nonlinear optical conductivity tensor only has three independent nonzero components, i.e., , , and [13,57,58]. Here, we set , , and , where , , , and is the effective thickness of graphene [13].
Figure 5.(a) and (b) Schematic graphs of an optical metasurface which can be dynamically switched to exhibit (a) high, (b) low SHG conversion efficiency through controlling the phase states of its constituent GeTe. The metasurface is composed of a square array of meta-atoms on a
The central large cuboid is composed of eight alternative layers of semiconductor CdO and phase-change material GeTe. The thickness of each CdO (or GeTe) layer is (or ). The relative permittivity of CdO can be described using the Drude model as [59], where , , and . The GeTe possesses the ability of rapidly switching between the crystalline and amorphous phases in a reversible way, and the phase transition can be achieved optically or through Joule heating and is proven to be reliable, fast, and repeatable [60,61]. Remarkably, the optical properties of GeTe in different phases are significantly different, and thus the optical properties of the CdO–GeTe multilayer can be efficiently manipulated through changing the phase states of GeTe. Due to the deep-subwavelength size, the CdO–GeTe multilayer can be approximately homogenized as an effective anisotropic medium, whose component of effective relative permittivity can be expressed as [62]. In Fig. 5(c), we plot the real part of , i.e., , as a function of wavelength when the GeTe is in the crystalline (blue lines) or the amorphous (red lines) phase. We see that at the wavelength of 6.93 μm in the amorphous phase, in which case the imaginary part of , i.e., , is not very large (see the inset). Under this circumstance, the central CdO–GeTe multilayer in the amorphous phase can be viewed as an effective anisotropic ENZ medium with . However, when the GeTe is switched to the crystalline phase, we have and , indicating the breakdown of the ENZ condition.
This unique property provides us with a route to dynamically control the occurrence and quenching of SHG of the optical metasurface through switching the phases of GeTe, as schematically shown in Figs. 5(a) and 5(b). When the GeTe is in the amorphous phase, the CdO–GeTe multilayer effectively serves as an anisotropic ENZ medium, which can generate near-zero local electric fields near the surfaces normal to the direction (see field distributions in Appendix A.5). The graphene placed there will be hard to “see” the incident light, and therefore the SHG from the graphene is suppressed [Fig. 5(b)]. Interestingly, when the GeTe is switched to the crystalline phase, the ENZ condition breaks down, and the SHG from the graphene emerges [Fig. 5(a)].
For numerical verification, we set the relevant parameters of graphene as , , and at the operating wavelength of 6.93 μm [13]. At this wavelength, the second-order nonlinear susceptibility of graphene () is several orders of magnitude larger than those of CdO () [22], GeTe () [63], and [13], and therefore, here we consider the nonlinear responses of graphene only. In practical implementations, the graphene can be replaced by other optical materials with large nonlinear susceptibility such as gallium arsenide [64]. Figure 5(d) presents the SHG energy from the metasurface with GeTe in the crystalline (green bars) or amorphous (yellow bars) phase under illumination of normally incident light (polarized in the direction). In each phase, the left and right bars correspond to the case with the central CdO–GeTe multilayer considered as an effective medium and the case with the original multilayered structure, respectively, showing good agreement. We note that the SHG energy in the amorphous phase is enlarged by 10 times for better visualization. We see that the SHG energy in the amorphous phase is much lower than that in the crystalline phase, demonstrating the dynamically controllable SHG quenching effect in the optical metasurface.
We note that the pump laser used to enhance the SHG of graphene generally will not inadvertently change the phase of GeTe for two reasons. First, the intensity of the laser pulse that can effectively enhance the SHG of graphene is generally lower than that required to change the phase of GeTe because the transition temperature of GeTe is high [65,66]. Second, the activation energy of GeTe is [67]. This indicates that the wavelength of the laser pulse used to change the phase of GeTe is required to be short, e.g., [65,67], much shorter than the operating wavelength for SHG quenching (i.e., 6.93 μm).
In practical implementation, the constituent materials of the metasurface can be flexibly chosen. For example, graphene could be replaced by other two-dimensional materials with pronounced nonlinear optical properties, such as transition metal dichalcogenides (TMDs) including , , and [68,69]. The CdO–GeTe composite could be replaced by other materials such as indium tin oxide, aluminum-doped zinc oxide, and metal–dielectric composites operating at the ENZ frequency [70,71], but the drawback is that the ENZ wavelength cannot be dynamically controlled.
5. DISCUSSION AND CONCLUSION
The SHG conversion efficiency of a nonlinear material strongly depends on its molecular-scale constituents and their macroscopic ordering [8]. Based on this property, a conventional approach to manipulate the SHG signal is to engineer the local and global symmetries of meta-atoms in metamaterials and metasurfaces [72–74]. In addition, the enhancement of SHG due to surface plasmon polariton resonance in metallic nanostructures is very sensitive to each meta-atom’s geometry, dielectric environment, and polarization of incidence, thus providing another approach to manipulate the SHG signal [75,76]. Here, unlike the previous endeavors, our approach paves a route to control the SHG conversion efficiency through exploiting the near-field effect of a linear material [77], that is, the extraordinary local fundamental electric fields of an ENZ particle, thus further realizing the unique SHG quenching effect. Our approach could be extended to other nonlinear optical processes, including sum- and difference-frequency generation.
It is noteworthy that TMDs, such as , , and , are capable of controlling the SHG conversion efficiency [68,69]. Due to their unique electronic band structures and strong light–matter interactions at the atomic scale, the SHG process of TMDs is sensitive to the number of layers, crystallinity, and interlayer coupling in TMD heterostructures [68,69,78–81]. Interestingly, it was demonstrated that the SHG is quenched in even layered as it belongs to the centrosymmetric space group [78,79]. Fundamentally different from the SHG quenching with TMDs, the SHG quenching by the ENZ media attributes to the macroscopic electromagnetic responses.
Our findings could have an impact on the experimental investigations of optical nonlinearity with ENZ media. A frequently used configuration in experiments is to exploit plasmonic or dielectric nanoantennas placed on a thin ENZ film [25,27–30,37]. Usually, the boost of nonlinear responses induced by the large field enhancement in the ENZ film is reported. Nevertheless, our results in Fig. 4 indicate that the thin ENZ film could significantly suppress the nonlinear responses from the nanoantennas. This breaks the traditional conception that the ENZ media are always responsible for SHG enhancement and enriches the understanding of the experimental observation of optical nonlinearity.
In conclusion, we have demonstrated the SHG quenching effect when a tiny nonlinear particle is placed very close to a subwavelength ENZ particle. The SHG quenching effect is found to originate from the extraordinary near-zero local electric fields occurring near the ENZ particle due to evanescent scattering waves and is proved to be universal with both isotropic and anisotropic ENZ particles, irrespective of their shapes. Based on the principle of SHG quenching effect, we have further demonstrated a kind of dynamically controllable optical metasurface integrated with phase-change material GeTe. Through changing the phase states of GeTe, the SHG from the metasurface can be switched on or off. Our work creates pathways for applications such as nonlinear optical switches and modulators based on the integration of nonlinear materials and ENZ media.
Acknowledgment
Acknowledgment. The authors sincerely thank Professor Yun Lai at Nanjing University for the insights he offered in discussions.
APPENDIX A
In Fig.
Figure 6.(a) Schematic layout of the configuration, which is the same as that in Fig.
Figure 7.(a) Normalized scattering cross section of SHG
Figure 8.Distributions of normalized fundamental electric-field amplitude in one unit cell of the optical metasurface on the plane that contains the graphene. The metasurface is the same as that in Fig.
References
[1] M. Kauranen, A. V. Zayats. Nonlinear plasmonics. Nat. Photonics, 6, 737-748(2012).
[2] N. M. Litchinitser. Nonlinear optics in metamaterials. Adv. Phys. X, 3, 1367628(2018).
[3] D. Smirnova, Y. S. Kivshar. Multipolar nonlinear nanophotonics. Optica, 3, 1241-1255(2016).
[5] G. Li, S. Zhang, T. Zentgraf. Nonlinear photonic metasurfaces. Nat. Rev. Mater., 2, 17010(2017).
[44] H. C. van de Hulst. Light Scattering by Small Particles(1981).
[45] D. J. Griffiths. Introduction to Electrodynamics(1999).
[47] N. Yu, F. Capasso. Flat optics with designer metasurfaces. Nat. Mater., 13, 139-150(2014).
Set citation alerts for the article
Please enter your email address