Localized surface plasmons (LSPs)^[1,2] are collective charge oscillations confined within the interface between the noble metal nanoparticles and the surrounding dielectric medium, and the resonance is sensitively related with the structural geometry, material, and environment^[3,4]. It has been proved that LSP nanostructures with spatial asymmetry own different coupling behaviors for polarized light of opposite handedness and exhibit artificially designed chiral responses^[4–8].

Nanogap structures, such as dimers, bowties, and pairs of semicircles and rods, can confine the electromagnetic field into the narrow gaps or at the sides of particles with different incident polarizations^[9,10]. Due to this polarization-sensitive resonance, the chiral responsibility of nanogap structures is studied. Until now, planar nanogap structures, including dolmen structures^[11], split rings^[12], and vortex gaps^[13], were developed for giant chiral response.

Recently, metal–insulator–metal (MIM) three-dimensional (3D) nanostructures have attracted extensive concerns in the researches of plasmonics. It has been demonstrated that, by stacking layer by layer, the vertically aligned plasmonic structures can induce complicated charge distribution within 3D structures^[14–16]. Furthermore, the twistedly stacked MIM structures can lead to the 3D-twisted charge distribution and realize higher-order or hybridized modes for chiral response^[17–19].

In this work, we studied the chiral response with 3D-stacked nanogap structures. This chiral MIM nanostructure consists of two vertically and twistedly aligned nanodisks with a nanogap, resulting in hybridized charge distribution between two layers. Excited with right-handed circular-polarized (RCP) and left-handed circular-polarized (LCP) plane waves, the stacked structure twisted at 60° exhibits plasmonic coupling behaviors with/without gap modes, respectively, showing a giant chiral response (ΔT) of 60% at the wavelength of 1550 nm. This chiral response can be further tuned by controlling the particle size and the insulator layer of the structure.

To understand the hybridization mechanism, plasmonic modes in nanogap structures are investigated firstly. It has been previously proved that nanogap structures own different LSP coupling performances corresponding to the alignment between the polarization and gap^[10,20]. When the incident polarization is perpendicular to the gap, the near-field intensity is confined within the gap area, called as gap mode [shown in Fig. 1(a), panel (1)]. When the incident polarization is parallel with the gap, the charge distribution and the near-field intensity are located at the sides of the components, which is similar to LSP coupling behavior of single particles (SPs), named the SP mode here [shown in Fig. 1(a), panel (2)]. By twistedly stacking the two layers, the charge distribution in each gap layer is impacted and distorted by each other, leading to the hybridization of plasmonic resonance [shown in Fig. 1(a), panel (3)].

Figures 1(b) and 1(c) illustrate the proposed periodic MIM nanostructure. The 3D finite-difference time-domain (FDTD) method is carried out using commercial software (FDTD Solution, Lumerical Inc.) to examine the charge distribution, near-field profile, and chiral response. For normal incidence, normalized LCP and RCP plane waves with electric field E0=1V/m are selected to produce the chiral response, ranging from 900 nm to 2000 nm. A mesh size of 2 nm (along the x, y, and z directions) is set to ensure the accuracy of electric/magnetic-field calculations. To optimize the circular dichroism (CD) signal at the wavelength of 1550 nm, which is one of the most important wavelengths for optical communication, the thickness and radius of the Au nanodisks with the gap size of 50 nm are set at 50 nm and 220 nm, respectively. The optical constants of Au are taken from Johnson and Christy (900−2000nm)^[21]. Silica with the refractive index of 1.42 is selected as the substrate. The thickness of the insulator layer is set as 90 nm. The refractive indices for the surrounding medium, insulator layer, and gaps are also set as 1.42 to achieve a uniform environment for LSP coupling. The impact of the refractive index n of the insulator layer is further investigated by varying the value from 1.38 to 1.46. Periodic boundary conditions in the x and y directions are applied with the period Px,Py=800nm, and absorbing boundary conditions with a perfect matching layer (PML) are applied along the ±z directions. The twisted angle θ between the two layers is varied from 0°to90° with a step of 10° to investigate chiral response, as illustrated in Fig. 1(b). The distributions of charge and electrical and magnetic near-field are further investigated with varied twisted angles, and the corresponding chiral response is valued by CD^[22]: ΔT=TRCP−TLCP.(1)

The transmission spectra for LCP and RCP plane waves with different twisted angles are investigated firstly. As shown in Figs. 2(a) and 2(b), the LCP and RCP transmission spectra are the same for the twisted angle at 0° and 90°, since there is no structural asymmetry for these structures. For the twisted angle at 60°, the absorption peak around 1550 nm [marked in Fig. 2(a)] is discovered in the transmission spectrum for the LCP plane wave, and the transmittance is only 21.3%. In contrast, the transmittance for the RCP plane wave at the same wavelength is 83.5%. The transmission difference proves that there is a giant chiral response for this twisted MIM nanostructure.

Figure 2(c) illustrates the trends (marked as yellow dash lines) of transmission by twisting the structure in details. For the LCP plane wave, the absorption peak is red-shifted originally from 1152 nm (marked as red star) to 1785 nm with decreased θ from 90° to 0°. Meanwhile, the same absorption for the RCP plane wave is inverted to enhanced transmission with decreased θ. The giant transmission difference with the twisted angle at 60° at the wavelength of 1550 nm is noted as green stars in Fig. 2(c). Calculated CD signals (ΔT) impacted by varied structural parameters are concluded in Fig. S1 of Supplementary Material. Absorption spectra also perform chiral responsivity, shown in Fig. S2 of Supplementary Material.

The transmission difference indicates that the chiral response is impacted by the twisted angle. Correspondingly, charge distributions of the MIM structures twisted from 0° to 90° are analyzed at the wavelength of 1550 nm. The results with the twisted angle of 0°, 30°, 60°, and 90° are shown in Figs. 3(a)–3(d), respectively. For easy understanding, simplified charge oscillations with the same conditions are illustrated in Figs. 3(e)–3(h). With the twisted angle at 0°, the charge distributions of both LCP and RCP excitations are located at the sides of the particles, showing the typical SP mode [shown in Figs. 3(a) and 3(e)]. With the twisted angle at 90°, the SP mode exists in the bottom layer, and the gap mode exists in the top layer for both LCP and RCP excitations [shown in Figs. 3(b) and 3(f)]. By increasing θ from 0° to 30°, the SP mode still exists in the top layer, but the intensity of charge distribution for the RCP condition is weakened, causing the difference of the LSP coupling effect and transmittance [shown in Figs. 3(c) and 3(g)]. Noticeably, for θ at 60°, the top layer for the LCP plane wave shows the SP mode, while for the RCP condition it diversely shows the gap mode [shown in Fig. 3(d) and marked in the inset of Fig. 3(h)]. For the LCP condition, the charge distribution of the two layers is close to an antibonding-type coupling between two layers [marked in Fig. 3(d)], which is previously convinced in MIM structures^[23], while the SP mode of the bottom layer for the RCP condition performs ambiguously with lower intensity. Thus, the antibonding-type coupling for the RCP condition is much weaker than that for the LCP condition, leading to lower absorption and higher transmittance. Corresponding to this difference in coupling modes, the maximum transmission difference achieves 60% for θ at 60°.

The optical chirality C for monochromatic electromagnetic radiation in free space is expressed as^[24]C≡ε02E·∇×E+12μ0H·∇×H,(2)where E, H, ε0, and μ0 are the electric field, magnetic field, permittivity, and permeability, respectively. It has been proved that chiral metamaterials with spatial asymmetry own different near-field coupling behaviors with different polarizations of incident waves and can enhance the macroscopic chiroptical response^[25,26]. For chiral structures, the interference between electric dipole and magnetic dipole impacted by structures plays a critical role in the chiral response. The mixed electric dipole and magnetic dipole can be expressed asp=αE−iGH,(3)m=χH+iGE,(4)where α is the electric polarizability, χ is the magnetic susceptibility, and G is the mixed electric–magnetic polarizability, presenting a pseudo-scalar for chiral structures. E and H are the electric/magnetic local fields at the structures^[27,28]. Thus, the coupling between the spin angular momentum of the incident light and the structural chirality induces different electric/magnetic near-field distributions and impacts the interference between electric dipole and magnetic dipole, leading to the chirality-dependent optical response. Figures 4(a)–4(d) show |H| near-field distributions of the structures twisted with 0°, 30°, 60°, and 90° for LCP and RCP plane waves at 1550 nm, respectively. For the structures twisted with 0° [Fig. 4(a)] and 90° [Fig. 4(b)], |H| near-field distributions are symmetric between LCP and RCP conditions, and the 90°-twisted one owns higher |H| field intensity, which could be due to the greater coupling strength of magnetic resonance. For the structure twisted with 30° and 60°, the |H| field distributions perform differently for LCP and RCP plane waves due to the increased structural chirality, as shown in Figs. 4(c) and 4(d). The greatest difference for |H| near-field distributions is realized for the twisted angle of 60°, and the corresponding ΔT is over 60%. Such a coupling difference induced by structural chirality is also discovered in |E| near-field distributions shown in Fig. S3 of Supplementary Material. Correspondingly, Fig. 5(a) proves that the CD signal (ΔT) is maximum when the nanogaps are twisted at 60°, leading to the giant chiral response at 1550 nm, as marked in Fig. 5(b).

Finally, the tunability of chiral behavior of this structure is investigated. The size of the nanodisk impacts the chiral response, as shown in Fig. 6(a). For the structure twisted at 60°, the maximum CD signal (ΔT) shows red shift from 1538 to 1626 nm by increasing the radius of the nanodisk from 130 nm to 260 nm. The impact of the refractive index n of the insulator layer is also investigated in Fig. 6(b). The chiral response can be tuned from 1538 nm to 1626 nm when the n increases from 1.38 to 1.46, showing critical sensitivity to the change of the medium index.

We theoretically study the chiral response of twistedly stacked nanogap structures. With the twisted angle of 60°, LSP coupling performance with/without gap modes for RCP and LCP excitation at 1550 nm is discovered. The electrical and magnetic near-field distributions of this twisted MIM structure corresponding to different polarizations are further validated. As a result, giant chiral response (ΔT=60%) at the wavelength of 1550 nm is realized, which shows critical sensitivity to the change of the nanodisk size and the medium index. Our result reveals the hybridization of LSP coupling behavior in twistedly stacked nanogap structures, which perform as a kind of ideal chiral metamaterial.