• Photonics Research
  • Vol. 9, Issue 4, 484 (2021)
Shuxia Zhao1, Lei Shao2, Jianfang Wang3, Hai-Qing Lin2、4、*, and Wei Zhang1、2、5、*
Author Affiliations
  • 1Institute of Applied Physics and Computational Mathematics, Beijing 100088, China
  • 2Beijing Computational Science Research Center, Beijing 100193, China
  • 3Department of Physics, The Chinese University of Hong Kong, Shatin, Hong Kong SAR, China
  • 4e-mail: haiqing0@csrc.ac.cn
  • 5e-mail: zhang_wei@iapcm.ac.cn
  • show less
    DOI: 10.1364/PRJ.416015 Cite this Article Set citation alerts
    Shuxia Zhao, Lei Shao, Jianfang Wang, Hai-Qing Lin, Wei Zhang. Chirality-selective transparency induced by lattice resonance in bilayer metasurfaces[J]. Photonics Research, 2021, 9(4): 484 Copy Citation Text show less

    Abstract

    Chiral optical responses of bilayer metasurfaces made of twisted metallic nanorods are investigated in detail with focus on the collective effect due to lattice resonance (LR). Using an analytical approach based on the coupled dipole method (supported by full wave simulation), we find optical chirality is dramatically increased by the coupling between localized surface plasmon resonances and LR. The collective effect results in significant chiral signal even for metasurfaces made of achiral unit cells. The interlayer coupling generally destroys the Wood’s anomaly and the associated transparency. While making use of Pancharatnam–Berry (PB) phase and propagation phase, one can modulate the optical activity effectively and achieve chirality-selective transparency induced by LR in a designed structure with a g-factor of absorption as high as 1.99 (close to the upper limit of 2). Our studies not only reveal a new mechanism of modulating chiral optical response by combination effects from PB phase, propagation phase, and LR, but also give a quantitative relationship between the geometry configuration and chiral optical properties, thus providing helpful guidance for device design.

    1. INTRODUCTION

    Chirality is of great importance in fundamental science, material design, biomedicine, and so on. Chiral optics associated with natural chiral molecules is limited by the fixed geometrical structure and weak optical signal mainly in the ultraviolet (UV) range. Plasmonic clusters supporting localized surface plasmon resonances (LSPRs) interact with visible light strongly, which can be controlled by the shape, size, position, and permittivity of the objects [1]. Plasmonic nanostructures can be precisely designed with complex geometric patterns and thus provide many opportunities for exploring light–matter/structure interaction, in particular the chiral nature of photonics. In the past decade, many types of chiral clusters have been explored to study circular dichroism (CD), such as helical particles [2,3], tetramers [4], hybrid L-shaped resonators [5], U-shaped antennas [6,7], hybrid rod-sphere structures [8], various dimers [914], and others [15,16]. Using intrinsic chiral clusters as building blocks, chiral metamaterials/metasurfaces [1723] can discriminate right- and left-circularly polarized (RCP and LCP) light, and they can even realize selective transmission of different circularly polarized light [2426]. Furthermore, 2D planar metasurfaces with extrinsic chirality can realize asymmetric transmission/reflection [2730] as well as modulation of light’s wavefront, direction, and polarization [3134], which is promising in developing metalenses. The responses of metasurfaces of stacked layers are discussed in the view of photonic systems [19,20,35,36]. Besides designing nanostructures, a structured light field (superchiral light) [37,38] is proposed to increase light–material interaction, which is promising in amplification of chiral molecules’ CD signal [39].

    Collective effects {such as lattice resonance (LR) [4044]} and interference effects (such as Fano resonance [45] and electromagnetic induced transparency) modulated by electromagnetic (EM) field phase play important roles in controlling the optical properties. However, the study of the combination of those effects in the chiral metasurface area is in the preliminary stage [22,4651]. Many important issues need systematic studies. For example, it is known that the geometric configuration, such as the relative orientation/position of the substructure, affects the collective resonance [43,50,52], phases of the EM field {Pancharatnam–Berry (PB) [31,53,54] and propagation phase}, and mirror symmetry breaking (related to chiral response) [4,14,52]. Much effort is needed to investigate the physical consequence of the interplay of these effects, the resulting (chiral) optical properties, the applications in optics modulation, and device design.

    Here we take an archetypical chiroptical nanorod dimer [9,55] as the building block of the metasurfaces to explore the combination effects of collective resonance and interference on the chiral optics. Using an analytical approach based on a coupled dipole model (CDM) and finite-difference time-domain (FDTD) simulation, we find that the optical chirality of the metasurface is dramatically increased by the coupling between LSPRs and LR. The collective effect results in significant chiral signal, even for metasurfaces made of achiral unit cell-orthogonal nanorod dimers. By using phase (PB phase and propagation phase) modulated LR, different chiral responses can be obtained; in particular, chirality-selective transparency [100% transmission and a g-factor of absorption up to 1.99 (close to the upper limit of 2)] associated with Wood’s anomaly is realized in a specially designed structure, though in general the interlayer coupling destroys the Wood’s anomaly. Our studies reveal a new mechanism of tuning the chiral optical properties based on collective effect assisted by phase modulation (PB phase and propagation phase), and the analytical results (verified by numerical simulation) clearly reveal the quantitative relationship between the geometry configuration and chiral optical responses, providing modern photonic applications such as circular polarizers, optical communication, and quantum information processing.

    The paper is organized as follows. First we introduce the analytical formula to calculate the transmittance of metasurfaces. Second, we discuss the chirality of twisted nanorods and LR in monolayer metasurfaces. Third, bilayer metasurfaces’ extraordinary chiral transmission induced by LR is discussed and the general relationship between geometry and chiral optical response is given. Finally, we discuss the extraordinary chiral transmission of cases with higher-order resonance, square lattices, and the optical activity beyond the coupled dipole limit.

    2. THEORETICAL FORMULA

    A. Coupled Dipole Model

    Scheme of the metasurface made of arrays of twisted nanorod dimers.

    Figure 1.Scheme of the metasurface made of arrays of twisted nanorod dimers.

    For simplicity, the silver nanorods are modeled as ellipsoids, whose long- and short-axis radii are 80 and 30 nm. The longitudinal and transverse LSPRs’ resonance wavelengths are at 590 and 355 nm. In the wavelength range 450–700 nm, the dominant mode is due to longitudinal LSPR. Thus, the component of the dipole polarizability tensor along the long axis of the nanorod is considered. Higher-order multipole modes and transverse dipole modes are neglected. Furthermore, the distance between two rods in the nearby/same unit cell is 3 times larger than the long/short axis radius. As a result, the CDM [1,56] is used to describe the interaction between the nanorods: 1αPnA=F^A·(EnA+mnG^nmAAPmA+mG^nmABPmB),1αPnB=F^B·(EnB+mnG^nmBBPmB+mG^nmBAPmA),where PnA/B and EnA/B are the induced electric dipole moments and incident electric field at the position of rnA/B=(nxΛx,nyΛy,0/z0), nx and ny are integers, and A and B represent the upper and lower layer. F^A/BPnA/BPnA/B/|PnA/B|2 is the projection tensor. G^nmAA=G^(|rnArmA|), and the Green’s tensor is G^(r)=(k2U^+)eikr4πr,where k is the light momentum in free space and U^ is the 3×3 unit tensor.

    The optical response of each nanorod can be described by the polarizability using either an analytical [57] or semi-analytical [58] method. Here we directly use the multipole decomposition method [5961] to give a precise description. This method is promising for describing nanoparticles with arbitrary shape and size. The dipole polarizability can be obtained as D/E, where E is the electric field along the long axis of the nanorod and D is the total electric dipole moment calculated by the FDTD method (see subsection 2.B on FDTD simulation). The multipole decomposition analysis indicates that the total response of each nanorod is dominated by the total electric dipole moment, and higher multipoles are negligible.

    Because of the periodicity in the xy plane of the metasurface, the Bloch theorem guarantees the solution of Eq. (1) must have the form PnA/B=PA/Bexp(ik||r||n) and the incident electric field EnA/B=EA/Bexp(ik||r||n), including a phase factor involving the incident light’s momentum component k|| parallel to a 2D metasurface. Specifically, k|| is zero for normal incident light, i.e., EnA/B=EA/B=E0exp(ikzA/B), where E0 is the electric field amplitude vector. We consider two cases of (PA,PB)L/R for left-/right-circularly polarized (LCP/RCP) light with electric amplitude vector E0;L/R=[12,±i2,0]E0. After some calculation, we have PL/RA=22ei(kz0±ϕ)H12+e±iθH(ϕ)H(θ)H(ϕ)H122E0e^A,PL/RB=22e±iθH12+ei(kz0±ϕ)H(θ)H(θ)H(ϕ)H122E0e^B,where e^A/B is the unit vector along the long axis of the nanorod in the upper (A) layer and lower (B) layer. H12=G^xxABcosθcosϕ+G^yyABsinθsinϕ, H(θ)=1α(G^xxAAcos2θ+G^yyAAsin2θ), and G^AA=G^BB=m0G^mnAA is the interaction tensor for nanorods in the same layer. G^AB=G^BA=mG^mnAB describes the interaction of nanorods in different layers. With the help of the Weyl identity eikrr=i2π1kzqei(q||r||+kzq|z|)dq||,with kzq=k2q||2, and q=(q||,kzq) an arbitrary 3D vector in k-space, we obtain the interaction tensor as G^=j=0G^(rrj)=i8π2j=0k2U^qqkzqei[q||(rrj)+kzq|zzj|]dq||,where r is the observation point. Using the relationship j=0eiq||r||j=(2π)2ALδ(q||L),with A the area of lattice unit cell and L the 2D reciprocal-lattice vectors, the sum of Green’s functions in Eq. (3) can be written as G^AA=limz0(i2ALk2U^LLkzleikzl|z|i2πk2U^qqkzqeikzq|z|dq||),G^AB=i2ALk2U^LLkzleikzl|z0|,where kzl=k2L2 and L=(L,kzl) is a 3D vector. Then we get the conditions to make G^ divergent, namely, k=|L|. Furthermore, the Green’s function of the far field (rfar) due to those dipoles is G^far=i2ALk2U^LLkzleiLr+ikzl|zzj|.

    If the wavelength is larger than the lattice constant (k<min{Lx,Ly}), then all terms of the sums are evanescent, except that for the L(0,0)=(Lx=02πΛx,Ly=02πΛy), representing the electromagnetic wave propagating in the z direction. If k>min{Lx,Ly}, the transmission has higher-order (m,n) grating modes, whose directions are the same as L(m,n)=(m2πΛx,n2πΛy,kzl). Finally, the far field of the (0,0) order grating mode [rfar=(0,0,z) and z>z0] is E(0,0)A/B=j=0G(rfarrjA/B)PjA/B=G^(0,0);A/BfarPA/B=ik2Aeik(zzA/B)PA/B.

    After solving the electric field, we can calculate the (0,0) order transmittance as T(0,0)=P(0,0)Pinc=|E0eikz+E(0,0)A+E(0,0)B|2|E0eikz|2=1σA+k24A2|E0|2|PA+eikz0PB|2.

    Note that the extinction cross section (σ) in each unit cell is of the form σ=k|E0|2Im(E0A*PA+E0B*PB).

    Similarly, the reflectance can be calculated as R(0,0)=P(0,0)(R)Pinc=|E(0,0)A+E(0,0)B|2|E0eikz|2=k24A2|E0|2|PA+eikz0PB|2.

    The absorption of the metasurface is defined as Ab=1TR. Then the absorption g-factor g=2(AbLAbR)/(AbL+AbR), and AbL/R is the absorption in the presence of LCP/RCP. The analytical approach based on the CDM gives a clear physical picture, and the quantitative results are supported by the FDTD simulation as seen below.

    B. FDTD Simulation

    In FDTD simulation, the circularly polarized light is generated by the x- and y-direction polarized plane waves with ±π/2 phase difference. For a nanorod dimer, perfectly match layer (PML) boundary conditions are used in all boundaries of three dimensions. For metasurfaces, the PML boundary conditions are used in the z direction, while periodic conditions are used in the x and y directions. The mesh cell is cube with size of 4 nm, and silver’s refractive index in the visible light range is obtained from Ref. [62].

    To calculate the polarizability (α) of one single rod used in CDM, we use the multipole decomposition method [5961] based on FDTD simulation of one rod: p=P(r)dr,m=iω2[r×P(r)]dr,T=iω10{2r2P(r)[r·P(r)]r}dr,for electric, magnetic, and toroidal dipole moments, where r belongs to the total ellipsoid. P(r)=[ϵ(r)ϵ0]E(r), where E(r) can be directly obtained from the FDTD method. The total electric dipole moment used in this paper is D=p+ikcT,with c the velocity of light. The polarizability is calculated as α=D/E, where D and E are the complex amplitudes of the dipole moment and incident electric field in the long-axis direction.

    3. RESULTS AND DISCUSSIONS

    A. Optical Responses of a Single Nanorod Dimer and a Monolayer Metasurface: Localized Surface Plasmon Resonance and Lattice Resonance

    Extinction cross sections of a dimer with twist angle π/4 and separation z0=200 nm for LCP and RCP light. Solid/dashed lines are the results calculated by CDM/FDTD.

    Figure 2.Extinction cross sections of a dimer with twist angle π/4 and separation z0=200  nm for LCP and RCP light. Solid/dashed lines are the results calculated by CDM/FDTD.

    For comparison with new features of chiral optics due to LR in bilayer metasurfaces, we discuss the optical response of a monolayer metasurface made of silver nanorods with a rotational angle θ with respect to the x axis. The CDM gives 1αPnA=F^A·(EnA+mnG^nmAAPmA).

    In the presence of normally incident LCP/RCP light, the effective dipole moment is PL/R=2e±iθ2H(θ)E0(cosθe^x+sinθe^y).

    Transmittances of monolayer metasurfaces with θ=0,π/6,π/3,π/2 for LCP and RCP light.

    Figure 3.Transmittances of monolayer metasurfaces with θ=0,π/6,π/3,π/2 for LCP and RCP light.

    B. Chiral Optical Response of Bilayer Metasurface: The Combination Effects from Lattice Resonance, PB Phase, and Propagation Phase

    1. Chirality-Selective Transparency

    As for bilayer metasurfaces, nanorods on the two layers may have different rotation angles θ and ϕ. So metasurface’s intrinsic chirality leads to different transmission for LCP and RCP light. The transmittances of a typical bilayer metasurface are displayed in Figs. 4(a) and 4(c) for those from FDTD simulation and CDM calculation. Different responses to LCP and RCP light are clearly seen, and there is very good agreement between the results from FDTD simulation and those based on CDM. Interestingly, a quite different response can be found as shown in Figs. 4(b) and 4(d) (for FDTD and CDM). Here one sees chirality-selective transparency. To understand the physics involved, it is helpful to rewrite Eq. (3) in the basis suitable for circularly polarized light. Then one has PL/RA=E02Δ(θ,ϕ){[ei(kz0±ϕθ)H12+ei(±θθ)H(ϕ)]e^L+[ei(kz0±ϕ+θ)H12+ei(±θ+θ)H(ϕ)]e^R},PL/RB=E02Δ(θ,ϕ){[ei(±θϕ)H12+ei(kz0±ϕϕ)H(θ)]e^L+[ei(±θ+ϕ)H12+ei(kz0±ϕ+ϕ)H(θ)]e^R},where Δ(θ,ϕ)H(θ)H(ϕ)H122(θ,ϕ) and e^L/R=22(e^x±ie^y). In our modulation of optics by metasurfaces, two important factors are phase and LR. We first take a close look at various phases in the above equations. We take PRB=E02Δ(θ,ϕ){[ei(θ+ϕ)H12+ei(kz02ϕ)H(θ)]e^L+[ei(ϕθ)H12+eikz0H(θ)]e^R]}as an example to explain the origin of the phases. Obviously, kz0 is the propagation phase for nanorod B. For normal incidence, there is additional PB phase 2ϕ for the left-circular component (in the second term). The term proportional to H12 is due to the interaction from nanorod A. Transforming from coordinates associated with nanorod A to those of B leads to a phase of ϕθ (in the third term). Also additional PB phase 2ϕ for the left-circular component leads to ϕθ2ϕ=(ϕ+θ) (in the first term). The modulation of the PB phase and propagation phase provides a useful method to tune the chiral optics as shown below. The case of θ+ϕ=π.

    In this case, H(θ)=H(ϕ). (PLA+eikz0PLB)·e^R=(PRA+eikz0PRB)·e^L. Under the additional condition of z0=λ/2, (PLA+eikz0PLB)·e^L=(PRA+eikz0PRB)·e^R. Then identical transmittances for LCP and RCP light for all θ=πϕ are obtained, which has been verified by FDTD simulation (not shown here).

    Transmittance of a bilayer metasurface with parameters z0=200 nm, Λx=600 nm, Λy=300 nm. (a), (c) θ=π/3 and ϕ=π/2; (b), (d) θ=π/6 and ϕ=5π/6. (a), (b) Results based on FDTD simulation; (c), (d) those based on the CDM.

    Figure 4.Transmittance of a bilayer metasurface with parameters z0=200  nm, Λx=600nm, Λy=300  nm. (a), (c) θ=π/3 and ϕ=π/2; (b), (d) θ=π/6 and ϕ=5π/6. (a), (b) Results based on FDTD simulation; (c), (d) those based on the CDM.

    We then discuss LR and the related Wood’s anomaly. From Eqs. (3) and (19), one can see that the quadratic terms of G^xxAA or G^yyAA in the denominator of (PL,PR) are zero. So the denominator and numerator only have linear terms of G^xxAA and G^yyAA. If k=|L(±1,0)|, only G^yyAA is divergent. In general, the divergent terms in the denominator and numerator cancel each other out (“/”), leading to the disappearance of Wood’s anomaly as seen in Figs. 4(a) and 4(c). Compared with the results of monolayer as shown in Fig. 3, the interlayer coupling generally destroys the Wood’s anomaly and related transparency, while in some specific conditions, the divergent terms in the numerator may be cancelled by the phase modulation. The coefficient G^yyAA in numerator for LCP(+) and RCP(-) is ei(kz0±ϕ)sinθsinϕe±iθsin2ϕ.

    Under the condition θ=πϕ=(π±kz0)/2, the above coefficient for LCP or RCP (not both) vanishes, which leads to the recovery of the transparency induced by Wood’s anomaly for one of the circularly polarized lights. It is the chirality-selective transparency as seen in Fig. 4 (for θ=π/6, ϕ=5π/6, and z0=200  nm satisfying the condition). This transparency is obtained by modulating the PB phase and propagation phase [see Eq. (19)]. The nice agreement between the results from FDTD and CDM further verifies the analysis based on CDM.

    From Eq. (21), one sees that the condition for chirality-selective transparency depends on the geometric parameters such as the twist angle and the separation between the nanorods. The phase-modulated lattice resonance and related chirality-selective transparency are the general effect. Changing the dimer length or the refractive index of the surrounding medium leads to a shift of the frequency of the LSPR of a single ellipsoid. The main physics of chirality-selective transparency remains unchanged. The case of ϕθ=π/2.

    (a), (b) Transmittance of a metasurface with achiral building blocks [i.e., unit cell with nanorod dimer of twist angle π/2 (θ=π/4 and ϕ=3π/4)] based on FDTD simulation. (a) Λx=600 nm,Λy=350 nm,z0=150 nm (the condition of chirality-selective transparency is satisfied); (b) Λx=400 nm,Λy=350 nm,z0=150 nm (the condition of chirality-selective transparency is violated). (c)–(f) The electric field distribution (|E|) for nanorods A and B at wavelength of 600 nm. (c), (d) The distribution at the x–y cross-section plane with z=0 (rod A)/150 nm (rod B) for LCP light; (e), (f) the distribution at the x–y cross-section plane with z=0/150 nm for RCP light. The other parameters are the same as those in (a). (g) The g-factor corresponding to (a). (h) The transmittance versus the long-axis radius of the nanorods at a wavelength of 600 nm.

    Figure 5.(a), (b) Transmittance of a metasurface with achiral building blocks [i.e., unit cell with nanorod dimer of twist angle π/2 (θ=π/4 and ϕ=3π/4)] based on FDTD simulation. (a) Λx=600  nm,Λy=350  nm,z0=150  nm (the condition of chirality-selective transparency is satisfied); (b) Λx=400  nm,Λy=350  nm,z0=150  nm (the condition of chirality-selective transparency is violated). (c)–(f) The electric field distribution (|E|) for nanorods A and B at wavelength of 600 nm. (c), (d) The distribution at the xy cross-section plane with z=0 (rod A)/150 nm (rod B) for LCP light; (e), (f) the distribution at the xy cross-section plane with z=0/150  nm for RCP light. The other parameters are the same as those in (a). (g) The g-factor corresponding to (a). (h) The transmittance versus the long-axis radius of the nanorods at a wavelength of 600 nm.

    In general, the phase states of transmission/reflection field depend on the detailed geometric parameters and wavelength, and the transmission and reflection waves are elliptically polarized. Using Eqs. (3), (10), and (12), one can see that at the lattice resonance (λ=600  nm), transmission and reflection waves are circularly polarized with equal electric field amplitude in the x/y direction for transmitted and reflected far fields. For the RCP incident field, transparency appears, and the transmitted wave possesses phase difference (phase of the electric field in the x direction minus that in the y direction) π/2. For the LCP incident field, the transmission and reflection waves possess phase difference π/2 and π/2.

    2. Higher-Order Lattice Resonance

    Transmittance of LCP/RCP light associated with higher-order resonance based on FDTD simulation. (a) θ=π/4, ϕ=3π/4, Λx=1200 nm and Λy=300 nm, z0=150 nm. k=|L(±2,0)|. (b) θ=π/4, ϕ=3π/4, Λx=931 nm and Λy=781 nm, z0=150 nm. k=|L(±1,±1)|.

    Figure 6.Transmittance of LCP/RCP light associated with higher-order resonance based on FDTD simulation. (a) θ=π/4, ϕ=3π/4, Λx=1200  nm and Λy=300  nm, z0=150  nm. k=|L(±2,0)|. (b) θ=π/4, ϕ=3π/4, Λx=931  nm and Λy=781  nm, z0=150  nm. k=|L(±1,±1)|.

    The conditions for the chirality-selective transparency related to higher-order lattice resonance (k=|L(0,±n)|) are ei(kz0±ϕ)cosθcosϕe±iθcos2ϕ=0,which leads to θ=πϕ=±kz0/2.

    Here we would like to point out that the chirality-selective transparency associated with k=|L(±n,±m)| (n,m0) is absent. Unlike the cases of LR under the conditions (k=|L(±n,0)|)/(k=|L(0,±n)|), both G^xxAA and G^yyAA are divergent for k=|L(±n,±m)| (n,m0). The coefficient of the term (G^xxAAG^yyAA, G^xxABG^yyAB) in the denominator is 2cosθcosϕsinθsinϕ(cos2θsin2ϕ+cos2ϕsin2θ),which is zero only for θ=ϕ. If θϕ, the transparency associated with Wood’s anomaly appears for LCP/RCP light as shown in Fig. 6(b). If θ=ϕ, the metasurface shows an achiral response. The small dip/peak at 465.5 nm in Fig. 6(b) is LR with k=|L(±2,0)|. This wavelength is far from that for local plasmonic resonance, so the energy transmitted is very large. Because Eq. (21) is not satisfied, 100% transparency cannot be achieved.

    3. Bilayer Metasurface with Square Lattice

    The metasurface with a square lattice (Λx=Λy) possesses higher lattice symmetry than that with a rectangular lattice, which brings about new/different features in optical response. For the metasurface with a square lattice, G^xxAA=G^yyAA and G^xxAB=G^yyAB, and the effective dipole moments are PL/RA=22ei(kz0±ϕ)G^xxABcos(θϕ)+e±iθ(1αG^xxAA)(1αG^xxAA)2[G^xxABcos(θϕ)]2E0e^A,PL/RB=22e±iθG^xxABcos(θϕ)+ei(kz0±ϕ)(1αG^xxAA)(1αG^xxAA)2[G^xxABcos(θϕ)]2E0e^B.

    There is no chirality-selective transparency because the necessary condition of transparency θ=ϕ eliminates the chirality of the structure.

    In the following, we consider the case of ϕθ=π/2. The term H12=G^xxABcos(θϕ) vanishes, indicating the effective decoupling of the upper and lower metasurfaces. Therefore, the bilayer metasurface shows achiral optical response like that of the monolayer metasurface. Also, the transparency due to Wood’s anomaly appears for both LCP and RCP light. From Eq. (24), we find that the scattering far field EL/RfarPL/RA+eikz0PL/RB1He^L/R, with H independent of θ or ϕ. In general, for the metasurface with a rectangular lattice in the presence of normally incident LCP/RCP field, the transmitted field contains both LCP and RCP components [see Eq. (19)]. For the metasurface with a square lattice in the presence of normally incident LCP/RCP field, the transmitted field contains only an LCP/RCP component. Moreover, the transmittance is independent of θ or ϕ.

    4. Beyond Coupled Dipole Approximation

    Transmittance of LCP/RCP light with Λx=600 nm and Λy=350 nm calculated by the FDTD method. (a) z0=90 nm, θ=7π/20, and ϕ=13π/20. (b) z0=75 nm, θ=3π/8, and ϕ=5π/8.

    Figure 7.Transmittance of LCP/RCP light with Λx=600  nm and Λy=350  nm calculated by the FDTD method. (a) z0=90  nm, θ=7π/20, and ϕ=13π/20. (b) z0=75  nm, θ=3π/8, and ϕ=5π/8.

    4. CONCLUSION

    We theoretically explore the chiral optics of bilayer metasurfaces made of twisted nanorods, focusing on the collective effect due to LR. Through detailed analytical calculation based on CDM and FDTD simulations, we find that combination effects from LR and phase modulation (PB phase, propagation phase) can bring about novel chiral optical responses, including chirality-selective transparency (recovery of Wood’s anomaly) and chiral response for metasurfaces with achiral unit cells. The theoretical results deepen our understanding of light–matter interaction at the nanometer scale. In particular, the analytical results (supported by numerical simulation) give a quantitative relationship between local geometric structure, lattice structure, and their (chiral) optical properties; for example, the condition for chirality-selective transparency θ=πϕ=π±kz02. Those special properties of metasurfaces could play important roles in optical communication, circular dichroism spectroscopy, and quantum information processing.

    References

    [10] X. Yin, M. Schäferling, B. Metzger, H. Giessen. Interpreting chiral nanophotonic spectra: the plasmonic Born–Kuhn model. Nano Lett., 13, 6238-6243(2013).

    Shuxia Zhao, Lei Shao, Jianfang Wang, Hai-Qing Lin, Wei Zhang. Chirality-selective transparency induced by lattice resonance in bilayer metasurfaces[J]. Photonics Research, 2021, 9(4): 484
    Download Citation