Robust enhancement of high-harmonic generation from all-dielectric metasurfaces enabled by polarization-insensitive bound states in the continuum

The emerging all-dielectric platform exhibits high-quality ($Q$) resonances governed by the physics of bound states in the continuum (BIC) that drives highly efficient nonlinear optical processes. Here we demonstrate the robust enhancement of third-(THG) and fifth-harmonic generation (FHG) from all-dielectric metasurfaces composed of four silicon nanodisks. Through the symmetry breaking, the genuine BIC transforms into the high-$Q$ quasi-BIC resonance with tight field confinement for record high THG efficiency of $3.9\times10^{-4}$ W$^{-2}$ and FHG efficiency of $4.8\times10^{-10}$ W$^{-4}$ using a moderate pump intensity of 1 GW/cm$^{2}$. Moreover, the quasi-BIC and the resonantly enhanced harmonics exhibit polarization-insensitive characteristics due to the special $C_{4}$ arrangement of meta-atoms. Our results suggest the way for smart design of efficient and robust nonlinear nanophotonic devices.

Nonlinear frequency conversion, such as harmonic generation that greatly broadens the accessible spectrum, are ubiquitous in applications from our daily life to cutting-edge science and technology 1 . Over the past decade, the emergence of artificially structured nanoresonators, i.e., metamaterials and metasurfaces, has revolutionized our perception in nonlinear optics 2 . Recently, all-dielectric metasurfaces provide a paradigm for efficient and flexible harmonic generation in nonlinear nanoscale optics due to the high index, low loss, and resonant field enhancement 3,4 . For example, the thirdharmonic generation (THG) was studied in dielectric metasurfaces made from centrosymmetric semiconductors such as silicon and germanium, and the secondharmonic generation (SHG) was reported in the III-V semiconductor metasurfaces. In particular, the ability to support both electric and magnetic Mie-type resonances is beneficial to enhance nonlinear responses. Various dielectric metasurfaces, supporting the magnetic dipole resonance [5][6][7] , toroidal dipole resonance 8 , and nanoradiating anapole resonance [9][10][11] , have been explored to increase the field confinement in tight volumes, resulting in a considerable increase of nonlinear conversion efficiency.
Most recently, a special optical resonant mode originating from the physics of bound states in the continuum (BIC) gains increasing interest in boosting the harmonic generation from dielectric metasurfaces. Optical BIC underpins the existence of perfectly confined states embedded into the radiation spectrum and thus suffers from inaccessibility of external excitations 12 . It was reported that dielectric metasurfaces with broken in-plane inversion symmetry can support ultrahigh quality (Q) factor resonances arising from distortion of symmetry-protected BICs [13][14][15] . Such quasi-BIC resonances with large field confinement have been employed to significantly enhance nonlinear responses such as THG and SHG by orders of magnitude on dielectric metasurface platform [16][17][18][19][20] . Es-pecially, efficient high-harmonic generation (HHG) up to 11th order was also experimentally observed from silicon metasurfaces empowered by the symmetry-protected BIC 21 . However, a direct consequence of these quasi-BICs excited by introducing structural symmetry perturbation is the sensitivity to polarization angle of the pump beam, which severely hinders the robustness in their practical applicability.
In this letter, we demonstrate the robust enhancement of the harmonic signals from all-dielectric metasurfaces that is not dependent on the incidence polarization. The nonlinear metasurfaces with in-plane symmetry perturbation are designed using the symmetry selection rules such that the quasi-BIC and optical response are polarization independent under normal incidence. The high Q factor resonance with strongly enhanced local fields leads to significant enhancement of the THG and fifthharmonic generation (FHG) signals, and the nonlinear harmonics show robustness to the polarization angles. Our results have important implications for novel nonlinear optical devices and restore a flexible design strategy for polarization-independent harmonic signals boosted by symmetry-protected BICs.
We consider dielectric metasurfaces comprising complex meta-atoms in the form of four silicon nanodisks on a silica substrate under the normal incident plane wave, as shown in Fig. 1(a). The silicon is adopted for the odd order nonlinear processes including THG, direct FHG, and cascaded FHG. The asymmetric nanodisks are designed to show a twofold rotational symmetry, i.e. C 2 symmetry. We perform a numerical analysis of the linear response using the finite element method solver in COMSOL Multiphysics in the frequency domain. The metasurfaces are considered as an infinite structure with perfect periodicity. The refractive index of silicon is taken from experiments and that of silica semi-infinite substrate is set as 1.45 22 . The square lattice is designed with period P = 1400 nm, nanodisk height h = 200 nm, center-tocenter separation between nanodisks d = 700 nm, and nanodisk radius r 1 varies while r 2 is fixed as 272 nm. Figure 1(b) shows the transmission spectra of the metasurfaces with nanodisk radius r 1 ranging from 228 nm to 272 nm under x-polarized normally incident plane waves. When their radii are perfectly identical, i.e. r 1 = r 2 , the metasurfaces support a genuine BIC marked by a circle, since the symmetry is incompatible with modes of free space. As expected, there is no dip in the transmission spectrum for the genuine BIC case, implying a completely decoupling of the mode from external radiation. The energy coupling of such symmetryprotected BIC system necessitates a certain radiation leakage channel, which is achieved through varying r 1 in this case. In the broken-symmetry metasurfaces, the resonance manifests itself as a narrow dip in the transmission spectrum. In Fig. 1(c), we plot the evolution of the Q factor of quasi-BIC resonance as a function of the asymmetry parameter α. Here α is defined by the difference in nanodisk cross sectional area ∆S/S. The radiative Q factor depends on α in the form of the in-verse quadratic law, Q = Q 0 α −2 , where Q 0 is a constant determined by the metasurface design independent of α [13][14][15] . As it can be seen, the Q factor experiences a dramatic decrease as α increases. Consequently, the nonradiating symmetry-protected BIC with infinite Q factor is effectively transformed into the radiating quasi-BICs with ultrahigh Q factors as the symmetry perturbation is introduced in the metasurfaces.
We use the quasi-BIC resonance to achieve the local field enhancement and highly efficient harmonics.
Here the broken-symmetry metasurfaces with r 1 =250 nm (α = 0.15) are adopted and the quasi-BIC resonant wavelength is designed at 1556 nm. In Fig. 2(a), the transmission spectrum of the metasurfaces demonstrates a quasi-BIC resonance with a Fano asymmetric line shape, which can be fitted by T = |a 1 + ia 2 + b ω−ω0+iγ |, where a 1 , a 2 , and b are constant real numbers, ω 0 is the resonant frequency, γ is the overall damping rate of the resonance. Then the Q factor can be evaluated by ω 0 /2γ as 1504. We further perform the multipolar decomposition of the optical linear response to uncover the nature of the quasi-BIC resonance. It is observed in Fig. 2(b) that toroidal dipole (TD) and magnetic quadrupole (MQ) are excited in the metasurfaces, with predominant contributions to quasi-BIC resonance. Such resonant feature can also be clearly observed from the near-field distribution in the inset. Note that the simulated enhancement of the nearfield amplitude at resonance can be up to 40 times inside the structures.
Then we concentrate on the nonlinear responses in the metasurfaces, starting with the THG process from the above-designed configuration. To calculate the harmonic generation from the structure, two coupled electromagnetic models within the undepleted pump approximation are exploited. First, the linear response at pump wavelength is solved for the local field distributions. Second, the induced nonlinear polarization inside the meta-atoms is employed as a source for the electromagnetic simulation at the harmonic wavelength to obtain the harmonic field. For the THG process, the nonlinear polarization in the silicon nnaodisks can be computed by P 3ω = 3ε 0 χ (3) (E·E)E ω , where ε 0 = 8.8542×10 −12 F/m is the vacuum permittivity, χ (3) = 2.45 × 10 −19 m 2 /V 2 is the third-order susceptibility of silicon, and E ω is the electric field at the pump wavelength. The THG conversion efficiency is normalized as η TH = P TH /(P FF ) 3 , where P TH is the total radiated power at the third harmonic and P FF is the pump power incident on the resonator at the fundamental frequency.
In Fig. 2(c), the THG efficiency reaches its peak of 3.9 × 10 −4 W −2 at the pump wavelength of 1556 nm corresponding to the position of the quasi-BIC resonance. Compared to wavelengths detuned from the resonance, the peak efficiency at quasi-BIC represents an enhancement of more than six orders of magnitude. Such a record high value of THG efficiency significantly outperforms the previously reported performance of the existing silicon resonators under normal excitation, as listed in Ta- ble 1. Essentially, since THG is predominantly a volume phenomenon that responds mainly to the electric component of the incident wave, the electric energy inside the structure at the pump wavelength can be used to interpret the THG efficiency characteristics 1 . As reported on the right axis of the Fig. 2(c), the electric energy within the pump wavelength regime exhibits a peak at 1556 nm. This reasonable agreement intuitively illustrates the large enhancement of THG efficiency at the quasi-BIC resonance resulting from the tight energy confinement. To analyze the dependence of the THG efficiency with the symmetry perturbation of the metasurfaces, the normalized efficiencies with respect to the asymmetry parameter α are exhibited in Fig. 2(d). The results show the expected trend of THG efficiency decreasing with the increase of α, since large deviation from the symmetric unit cell facilitates the radiation leakage of the symmetry-protected BIC. This behavior reveals the controllable THG efficiency by the magnitude of the symmetry perturbation, and thus the THG performance can be improved by fabricating small defects with specialized process.
The proposed dielectric metasurfaces that host a high Q factor quasi-BIC resonance provides a promising solution to generate high-harmonic signal on the subwavelength scale at moderate driving intensities. Here we consider the potential in boosting FHG process under normally incident plane wave. In a conventional way, a direct process involved with the fifth-order susceptibility χ (5) is used to obtain the FH signal. Accordingly, this direct FHG process is simulated through the nonlinear polarization P 5ω = 3ε 0 χ (5) (E · E) 2 E ω , where χ (5) of silicon is estimated as 5.6 × 10 −39 m 4 /V 4 . On the other hand, the FHG can be achieved by a cascaded process of the THG followed by degenerate four wave mixing (dFWM), as illustrated in Fig. 3(a). The cascaded process shows advantages of higher efficiency and extremely short time constant. In the simulation, the polarization of the cascade THG is defined as , where E ω and E 3ω are the electric field at fundamental frequency ω and TH frequency 3ω, respectively. The calculated FHG efficiencies of the two pathways are depicted in Fig. 3(b). At the pump intensity of 1 GW/cm 2 , the cascaded effects show the peak efficiency of 4.8 × 10 −10 W −4 at the quasi-BIC resonance wavelength, which is higher than that of the direct process of 1.5 × 10 −11 W −4 due to the relation between χ (5) and χ (3) by χ (5) 3. (a) The schematic of FHG via direct process, or through cascaded THG and dFWM. (b) The normalized conversion efficiency of FHG via direct and cascaded pathway for x-polarized normal incidence.
As the most intriguing part, we demonstrate that the proposed metasurfaces support polarization-insensitive quasi-BIC resonance and its optical response in either linear or nonlinear regime keeps unchanged under polarization variation. In most previous works that concentrated on BIC-enhanced harmonic generation, the quasi-BIC resonance can only be excited by a specific polariza- tion light and the significant harmonic signals can only be expected in this situation, following the selection rule of the so-called symmetry-protected BIC. Figs. 4(a)-(d) present the linear transmission spectra and the normalized nonlinear harmonic signals of the metasurfaces with various pump polarization angles. It is observed that the transmission spectrum does not change as the polarization angle varies from 0 to 2π, revealing the polarizationinsensitive feature of the quasi-BIC here, which leads to the THG and FHG enhancement robust to the polarization angles because of the direct dependence of nonlinear response with linear response. Such robustness can be unveiled by the symmetry in arranging the metaatoms, as illustrated in Fig. 4(e). As highlighted by the solid square, the unit cell shows C 2 symmetry. In the meanwhile, the unit cell constructed by neighboring nanodisks with the same radius r 1 possesses the C 4 symmetry, and this kind of C 4 symmetry also applies to the unit cell with nanodisks of radius r 2 , shown by the dashed squares. Since the C 4 symmetry is required for the polarization insensitivity, the special arrangement of meta-atoms in the proposed metasurfaces enables the polarization-insensitive optical responses in both the linear and nonlinear regimes.
In summary, we propose the robust BIC enhanced of THG and FHG in the dielectric metasurfaces composed of four silicon nanodisks. Through the symmetry perturbation in the unit cell, the quasi-BIC with high Q factor and tight field confinement are observed. Boosted by the quasi-BIC, the metasurfaces produces a high THG efficiency of 3.9 × 10 −4 W −2 at the pump wavelength of 1556 nm and outperforms the reported silicon nanodisk resonators. We also show that the FHG can be enhanced with efficiency of up to 4.8 × 10 −10 W −4 at a moderate pump intensity. Most importantly, the quasi-BIC resonance and the resultant nonlinear harmonics are independent on pump polarization, showing great prospects in designing metasurface and metadevices for nonlinear applications.