Robust tunable plasmon induced transparency in coupled-resonance finite array of metasurface nanostructure

Robust and dynamically polarization-controlled tunable plasmon induced transparency (PIT) resonance in designed finite-array nanostructures metasurface is demonstrated, where sharp resonance is guaranteed by design and protected against large geometrical imperfections even for micro-zone sub-array. By employing the explicit analysis of near-field characteristic in the reciprocal-space based on the momentum matching, and the far-field radiation features with point-scattering approach in real-space sparked from Huygens’s principles, the physics of interference resonance for plane-wave optical transmission and reflection of the metasurface is theoretically and thoroughly investigated. The distinctive polarization-selective and Q-tunable PIT shows robust features to performance degradations in traditional PIT system caused by inadvertent fabrication flaws or geometry asymmetry-variations, which paves way for the development of reconfigurable and flexible metasurface and, additionally, opens new avenues in robust and multifunctional controllable nanophotonics device design and applications.

Plasmon-induced-transparency (PIT) is a widespread wave scattering phenomenon associated with a peculiar coupled interference effect. Since its first demonstration 1 in the plasmonic nanostructures, PIT has been widely and thoroughly studied both theoretically and experimentally 2-4 due to a large variety of pivotal and potential applications in sensing, filter, switches, slow-light, and nonlinearity. However, the study of PIT with reconfigurable dynamic tunability poses several challenges. Many works reported tunable PIT resonance by employing the geometrical asymmetry/variations, which is inherently static and the control of channels for the PIT windows are rarely studied 3 . Favorably, phase-change-materials, graphene, magnetic-materials, and nonlinear dielectric nanostructures have auspiciously emerged as a promising alternative and platform for dynamic control of PIT, and are expected to complement or even replace plasmonic nanostructures for a wide range of potential applications 5 .
Recently, carefully designed plasmonic nano-antennas taking advantages of the polarization-selective excitation and evolution of the split-ring-resonator like unit have been proposed to shed light on dynamic PIT construction, not only providing a notable sequential polarization-controlled switchable PIT channels, but also flexibility and configurability 3 . Nevertheless, these devices perform structure-parameters sensitivity, which conversely challenges the state-of-art of fabrication and owns inferior tolerance to fabrication errors and flaws, leading to deficiency of robustness of the PIT resonance for practical applications. Tunable PIT with dynamic manipulation properties is unresolved for PIT with immunity against fabrication imperfections/flaws and geometrical variations. Topological optics is inspiring a new perspective on design of nanophotonics metasurface and devices 6,7 . The manipulation of photonics by employing nanoscale photonics structures with exotic properties has been well demonstrated with fascinating designs [6][7][8] . Due to their simple construction and efficient tunability, metasurface using nanoslits or slot antennas topology as building blocks has been identified as a versatile platform to study topological nature and manipulation of photonics for helicity dependent directional surface plasmon polaritons (SPP) excitation 9 , flexible control of in-plane SPP shaping and focusing or far-field radiation 10,11 , chirality of asymmetric transmission 12 , parity-time symmetry breaking 13 , and spin optics in metasurface 14 .
In this study, we introduce a class of compact slits nano-antennas array in metallic film that lift all these limitations and show their feasibility to be dynamically controlled by incidence polarization, where robust PIT with tunable quality factors is demonstrated. The physics of the generation, evolution, and annihilation of the PIT, and the immunity of the proposed system to the asymmetry geometry parameters variation and micro-zone-array are investigated and illustrated by exploring the dispersion characteristics in the real-space and the reciprocalspace. The simulated results are verified by the point-scattering approach model theoretical predictions, where the interference of the field in the slits and the far-field radiation features are revealed. The results pave the way for reconfigurable, flexible multifunctional metasurface and robust nanophotonics device design and applications.
Proposed metasurface and theoretical calculation method. The proposed topology consists of periodic rectangular nanoslits antenna array in gold film on substrate. Figure 1a illustrates the schematic of plasmonic nanostructures. A 50-nm-thick gold film was deposited on a glass substrate. The incidence is denoted as linear polarized light with an angle of polarization (AOP) α. Figure 1b shows the scheme of the unit cell slit's orientation and the corresponding coordinate system. The unit cell with perforated specifically arranged resemble slot antennas (L = 200 nm, W = 40 nm) is shown in Fig. 1c. The periods P x , P y are set as 600 nm equally along both x and y directions.
The geometrical parameters shift S and gap G are denoted in Fig. 1c, respectively. The rotation coordinate system denotes as u-v and the translational coordinate system as x-y are shown in the Fig. 1b. The coordinate vector system for the proposed metasurface is shown in Fig. 1c. For the case of a rectangular-shaped nano-slot in metal film with a high aspect ratio, the resonance and amount of SPP generation are significantly dependent on the orientation of the incident polarization [15][16][17][18] . Different from the metallic nanobars where the plasmonic resonance mode is excited with polarization along the long-arm, the SPP can only be excited when the polarization direction of the incident light is perpendicular to the metal slits.
The physics and model of optical transmission of periodic subwavelength holes array (SHA) nanostructures have been vastly and thoroughly studied, where several model such as the coupled-mode equation 19 , hybridwave model 20 , and the ab initio theory of Fano-formula 21 were proposed. Meanwhile, the resonance role of SPP modes was studied by reciprocal-space momentum-match conditions. For single-interface SPP, the dispersion relation can be described as: where ε d and ε m are the relative permittivity of the dielectric and the metal. The dispersion relation of the SPP Bloch mode generated under the Bragg coupling condition can be obtained by the momentum match: where ⇀′′ k in is the component of wave vector k parallel to the array surface, ⇀ u x and ⇀ u y are the unit vector in plane. Instead of straightforward description of the transmission, the reciprocal-space method relies on a priori definition of resonance as unperturbed smooth surface mode without the nanostructures. In 2005, Genet et al. presented a real-space Huygens description instead of the more conventional reciprocal-space description of interfering resonance features in the SHA 22 . However, the anisotropic designed holes arrays increasingly utilized in the metasurface for flexible control and modulation of electromagnetic wave are seldom treated. Recently, the SPP far-field radiation in an anisotropic SHA metasurface is demonstrated 11,23 . In the point-scattering model, the unit cell of the hole openings are treated as resonance point-scatter. In Fig. 1b, the angle ϕ between the v axis (also referred to as the normal direction of a nanoslit) and y axis is defined as the rotating angle of the nanoslit. By setting the first, lower-left nanoslit's center as the origin O of the whole system, the contribution of all the point-scatter to the transmission spectra can be processed based on the Huygens' Principles.
The disposition vector of the center of the first nanoslit in the jth (the mth row and the nth column) supercell The unit vector for the position vector is u j = ⇀ R j / ⇀ R j . In one unit cell, the relative coordinate of the center of the jth nanoslit's center relative to the first nanoslit's center in the same supercell is set as ⇀′ r s (Fig. 1c). The absolute coordinate of the jth nanoslit's center can be obtained as x sin φ. The full polarization behavior of the SPP is contained in the tensorial u j ⊗ u j nature of the elementary scattering matrix. Assuming that the scattering matrix of the far-field is spherical so that the intensity of the transmission for the metasurface can be calculated 11,22,23 , where the spectrum for m = 1, n = 1 are defined as a normalized coefficient T 0 . Under point-scattering limit, the shape factor S (|K|) can be replaced with a constant. This relates the present formulation to the two-dimensional Green function of the surface-wave with no need to investigate the full three-dimensional Green function half-space problem 22 . For a given hole, the incoming field is converted into two-dimensional surface waves that propagate away from the hole as a spherical (Huygens) wave. Based on the point-scattering-approach, the planewave optical transmission of the designed metasurface for arbitrary angles and polarizations can be evaluated and predicated theoretically.
The intensity of transmission spectra for both TM/TE plane-wave at normal incident for the symmetric/asymmetric structure are theoretically calculated using Eq.  23 . To verify the theory prediction, the full-field electromagnetic calculations were performed using the finite difference time domain (FDTD) method. The permittivity of Au was taken from the experiment data 24 , while the refractive index of the glass substrate was fixed at 1.45. The calculated reflectance spectra are shown in the Fig. 2b. Clear transparency windows with narrow peaks are observed at 904 nm wavelength, for both the system with and without broken-symmetry. Excellent agreement between the theory and the simulation is observed.

Discussions
Optical response and tunability of the nanostructure. The PIT is maintained and obtained in the system without geometry symmetry-breaking, which shows extraordinary immunity to the symmetry-breaking employed and usually required in the traditional plasmonic nanostructure using the resemble resonance nanobar/nanoslit unit. The in-plane interference of the surface wave by the scattering of light from the arranged nanoslits is attributed to the variation of the reflection spectra for different gaps of the nanoslits.
The results in the Fig. 3a show clear emerging and vanishing of the PIT peak for different gap values, where the PIT peak is shrinking and fading for increased gaps below 20 nm, while the notable regeneration of the PIT with broaden and enlarged peak is observed for further increased gaps. The result shows quite unique trends featured with unusual fading and revival of the PIT resonance when increasing the gaps G. Hence, once the distance for varied coupling is increased or decreased, the PIT resonances position moves. The oscillatory behavior as a function of varied gaps is also observed as theoretically predicted from the surface-wave mediated interference model, where the first-order interference between SPP and these out-of-plane component at the subwavelength www.nature.com/scientificreports/ scattering sites formed by the slits openings gives rise to the observed spectra 23 . The PIT-like resonance is maintained for vast range of gaps, which shows relatively robust performance to fabrication flaws or imperfection.
To further verify this, the rotation angles (ϕ 1 , ϕ 2 , ϕ 3 , ϕ 4 ) of the slits are slightly and randomly changed, where by theory prediction, besides of the shrink ratio values, the unchanged fixed peak position centered at 904 nm is still observed agreeing with the simulation results. At the same time, changing the distance (Shift S) between the nanoslits vertical to the incidence polarization will not significantly affect the position of the PIT resonance due to the fact that the slot-antenna's polarizationexcitation-selective property. Figure 3b shows the reflectance spectrum at normal incidence as a function of varied shifts of the nanoslit 1 (see Fig. 1c), where the PIT resonance effect shows slight changes for shifts below 40 nm, and the transparency window shows negligible and almost unchanged features for increased shifts. By carefully designing utilizing the polarization-selective bound-charge-oscillator like slot antenna unit 15,17 , the singular features of the PIT resonance in the system is reported for the first time, to the best of our knowledge, that the topologically robust PIT resonance is demonstrated in metasurface for both symmetric and asymmetric unit. Most works indicated that asymmetry was crucial to realize PIT when the interaction of dipole resonance and the plasmon resonance modes in the nanostructure is considered. Here we propose structure that is not limited by the asymmetry condition and PIT resonance with robustness features shows immunity to the asymmetry geometry variation.
For varied incidence polarization orientation, the tunability of the system with fascinating sequentially switchable PIT resonance peaks with tuned quality factor (Q-factor) and contrast is obtained, where dynamic PIT with tunable quality can be efficiently tailored by adjusting the incidence polarization. The result in Fig. 4a exhibits clear evolution of the PIT resonance, where broaden resonance associated with almost fixed resonance position is demonstrated for increased angles of polarization. The resonance quality factors of the system with minimum 9 nm linewidth can be gradually tuned by controlling the AOP. For asymmetry structures, the PIT resonance is also observed from the results in Fig. 4b. Q-tunable robust PIT in designed metasurface immune to the geometry symmetry/asymmetry is obtained.  www.nature.com/scientificreports/ Furthermore, there has been plenty of designed controllable light beams with tailored wave-front and vector vortex are realized using coded metasurface or spatial light modulator (SLM), the fruitful dynamic PIT resonance with striking features inspired by manipulating the incidence polarization shows remarkable potential for dynamic metasurface devices design and applications.
Dispersion manipulation and theoretical calculation of Bloch surface mode supported. To gain clear systematic properties of the system, the topological properties are investigated by exploring the dispersion of the nanostructure. The results in Figs. 5 and 6 show the dispersion calculated band for TM-polarization and TE-polarization for both k x , and k y momentum using the reciprocal-space theory (Eq. 2).
The PIT-like resonance is revealed origins from the directional excitation coupled Bloch-surface-wave. It is crucial to dynamically control PIT resonances while maintain robust profile features. However, typically emerged