Abstract
Metasurfaces, with their superior capability in manipulating the optical wavefront at the subwavelength scale and low manufacturing complexity, have shown great potential for planar photonics and novel optical devices. However, vector field simulation of metasurfaces is so far limited to periodic-structured metasurfaces containing a small number of meta-atoms in the unit cell by using full-wave numerical methods. Here, focusing on achiral meta-atoms only with electric polarizability and thickness far less than the wavelength of light, and ignoring the coupling between meta-atoms, we propose a general phenomenological method to analytically model the metasurfaces based on the assumption that the meta-atoms possess localized resonances with Lorentz-Drude forms, whose exact form can be retrieved from the full wave simulation of a single element. Applied to phase modulated geometric metasurfaces constituted by identical meta-atoms with different orientations, our analytical results show good agreement with full-wave numerical simulations. The proposed theory provides an efficient method to model and design optical devices based on metasurfaces.
© 2016 Optical Society of America
1. Introduction
Metallic metasurface is a thin layer of structured surface consisting of an array of planar metallic antennas with thickness far less than the wavelength of light. By tailoring the orientation and anisotropic resonances of meta-atoms, metasurfaces can be designed to control the phase, amplitude and polarization of light at the nanoscale, which has opened door for various planar photonic and optical devices with low manufacturing complexity [1–3]. They have been employed for ultra-thin wave plates [4–6], phase gradient plates [7–9], planar lenses [10,11], optical vortex plates [7,12] and holograms [13–15]. Among various types of metasurfaces, a special type of metasurfaces that manipulate the phase of light by the orientation of the antennas have shown great potential for practical applications due to the simplicity and robustness of the phase control. The underlying mechanism of the phase control is a Pancharatnam-Berry phase in the polarization state, which leads to strong spin–orbital interaction of light [16–19]. This concept has been used to realize spin controlled directional coupler for surface waves [20,21] and broadband holograms [22,23].
Up to date, theoretical modeling of metallic metasurfaces is mainly limited to full-wave numerical simulations based on the Finite-Difference Time-Domain (FDTD) method and the Finite Element Method (FEM). So far, most reported analytical models for metasurfaces, such as homogenization with a Lorentz-formed polarizability [23], electric surface impedance and magnetic surface impedance method have only been applied to uniform metasurfaces [6,24–26]. Very recently, an analytical model based on lumped impedance was applied to the modeling of Huygens’ metasurfaces [27,28]. For metasurfaces used as metalens, holograms and vortex plates that consist of a large number of different meta-atoms, there have been no efficient theoretical methods that could be used to accurately predict their optical performance. In this paper, we propose a general analytical theory to model geometric metasurfaces with arbitrarily distributed phase profiles. The method is shown to agree very well with, but much more efficient than the full-wave numerical simulations.
2. Analytical theory of metasurface
Without loss of generality, we consider a metasurface sandwiched between the homogenous and isotropic cover and substrate, as depicted in Fig. 1. Where, () denote tangential component of electric fields propagating along direction at the interface between the cover (substrate) and metasurface.
Because the metasurface has a negligible thickness h (), the overall tangential component of the electric fields is assumed to be continuous across the metasurface. Therefore the electric fields on the interfaces of the cover and substrate satisfy
For the sake of simplicity, we only consider metasurfaces made up of achiral meta-atoms and neglect the surface current induced by the magnetic fields of light. The tangential components of magnetic fields at the interfaces of the cover and substrate thus satisfy [29]
Where, is the electric susceptibility tensor of the metasurface. Due to their deep subwavelength sizes, assuming that the coupling between meta-atoms can be ignored, each antenna can be considered as an electric dipole located at. The susceptibility tensor can be written as,In Eq. (3), is the local polarizability tensor of the meta-atom located at. Now we consider a geometric phase metasurface consisting of an array of anisotropic antennas of identical geometry but different orientations as schematically shown in Fig. 2(a), where the Pancharatnam-Berry phase for circularly incident beam is determined by the orientation angle of each antenna.
The local polarizability tensor can be expressed as
where is the antenna orientation angle, and are two orthogonal components of the local polarizability tensor of the anisotropic antenna. Note that when is 0, the polarizability is reduced to a diagonalized tensor [, 0; 0 ]. The frequency responses of the two orthogonal components can be assumed to have Lorentz-Drude forms [30]. That iswhere is a parameter indicating how strong the meta-atom couples with the incident electric fields, c is the light speed in vacuum, and are the resonance frequency and the damping rate of the meta-atom, respectively.We can express susceptibility tensor and the electric field as the sum of various spatial frequency components through a Fourier expansion
where S = LxLy is the area of the metasurface. Lx and Ly are the smallest periods of the metasurface along the direction x and y. And, the discrete wavenumbersandare equal toand , respectively. Combining Eqs. (1)-(3) and Eq. (6), the Fourier expansions of electric and magnetic fields at the interfaces of the cover and substrate satisfy In Eq. (8), is the admittance tensor of the substrate, that relates the electric field and magnetic field in the media, and can be expressed asWhere,,andare the permittivity, magnetic permeability and refractive index of the substrate. The expression for the admittance tensor of the cover is given by replacing ‘S’ by ‘C’ in Eq. (9). A general expression of spatial frequency components of scattered fields produced by the metsurface can be obtained from Eqs. (7)-(8). That isWhere, I is identity matrix. Matrix and are defined asIt has been shown that the diffraction efficiency of metasurface can be dramatically enhanced through a three layer design, where the plasmonic antenna layer is on top of a homogeneous dielectric substrate layer of finite thickness d and a thick ground metal layer such that the metausrface operate in reflection mode shown Fig. 2(b). It has been shown that, with suitably designed thickness of the substrate layer, the interplay between the antenna resonance and the Fabry-Perot resonance can lead to very high efficiency over a broad bandwidth [23]. We next extend our analytical analysis to such a multilayer design. Based on Eqs. (7)-(8), the electric and magnetic fields at the interfaces of the cover and metasurface satisfy the following set of equations,
where is an effective admittance that take into account the reflection by the metal ground plane and is expressed as,Where,In the above expressions, is the complex reflection coefficient of TE (TM) polarized light at the surface of the ground metal plane. For a normally incident light, Eq. (14) can be written asEquation (15) shows that the effective admittance of the combined dielectric layer and the ground plane is dispersive as from the circuit model [31], which helps to cancel the intrinsic dispersion of the metal dipoles. Thus, with a proper thickness of the dielectric layer, a broadband and high-efficient optical device could be realized.3. Retrieval of Lorentzian parameters
For an antenna of a specific geometry, the Lorentzian parameters can be retrieved by fitting the full wave simulation of a simple periodic array of the antennas with our analytical theory. Here CST microwave studio software is used for the simulation. We consider a metasurface consisting of a periodic array of identical gold nanorods with the same orientation in a square lattice with period , as shown in Fig. 3(a). The semi-infinite cover and substrate are air and SiO2 with refractive indexes of 1 and 1.45, respectively. The relative permittivity of gold is modeled by Drude model. That is
For wavelength much greater than the period, we only need to consider the contribution of the zero-order Fourier component of the polarizability tensor. For a normally incident light from the cover side, combining Eqs. (4), (6)-(8), the zero-order transmitted and reflected electric fields, denoted by subscript 0, can be written aswhere , and are defined in Eq. (4). For the periodic-structured metasurface, all the meta-atom are identical. Consequently, the subscript mn shown in Eq. (4) is omitted. Rewriting Eq. (14) in the circularly polarized bases, we haveEquation (17) shows that only cross-circular polarizations in transmitted light can gain different additional phases for the normally incident beam with circular polarizations. The additional phases are non-dispersive geometric phases determined by the orientation angle of the meta-atoms.
Figures 3(b)-3(e) shows the transmission amplitude and phase of the X (Y) -polarized light normally incident from cover onto the metasurface that is schematically shown in Fig. 3(a). We fit the curves based on Eq. (5) and Eq. (16). The retrieved parameters of the nanorod are
The fitting results agree very well with the CST simulation. These fitting parameters will be used to model more complex metasurfaces consisting of antennas of spatially varying orientations for controlling the wavefront of light in the subsequent sections.4. Verification of the analytical theory with full-wave numerical simulations
By feeding the retrieved parameters obtained from Section III into the analytical model, we can analytically model metasurfaces consisting of antennas of arbitrary orientations for manipulating the wavefront of incident beam. The analytical modeling results are compared to full wave simulations for several representative metasurfaces.
4.1 Phase gradient metasurface
We first consider a metasurface with one-dimensional geometric-phase gradient. In a periodic unit cell, the phase gradient metasurface comprises N identical meta-atoms, whose orientations angle φ exhibits a constant gradientalong one direction. Figure 4(a) shows the simulated phase gradient metasurface with N equal to 8, made up of the same Au nanorods as that in Fig. 3(a). The periods of the array along X axis and Y axis are and, respectively. From Eq. (4) and Eq. (6), we can obtain an approximate expression of the susceptibility tensor of the phase gradient metasurface by only keeping the first order spatial harmonic terms. That is
where is the orientations angle of the first nanorod chosen in the periodic unit cell.In Eq. (18), we only consider the −1, 0 and + 1 order diffractive light for the phase gradient metasurface, whose wave vectors are given by ,andfor an incident light with wave vectors . Combining Eq. (10), Eq. (11) and Eq. (18), we can obtain the expressions of the transmitted and reflected electric fields for the incident light from the cover,
where subscripts −1, 0 and + 1 denote diffraction orders of the transmitted and reflected electric fields from the metasurface, respectively, andFor a right circularly polarized light normally incident on the metasurface shown on Fig. 4(a) from the cover, using the retrieved parameters of the nanorod, we can obtain the reflectivity and transmittivity of the zero-order and first-order diffractive light from Eqs. (19.1)-(19.3). Figures 4(b)-4(c) show that the analytical results agree well with those obtained from the CST simulation.
Furthermore, to check the applicability of Eqs. (19.1)-(19.3) for the obliquely incident light, we simulate a RCP with an incident angle of 30° from air onto the metasurface as shown in Fig. 4(a). Figure 5 shows that the analytical results still agree well with those calculated from CST simulation for the obliquely incident beam.
4.2 Reflective metasurface
The configuration of the reflective metasurface is shown in Fig. 6(a), which consists of an array of antennas atop a homogeneous dielectric layer and a ground plane. For simplicity, we confine ourselves to the discussion of incident light at normal direction. The admittance tensor of the combined dielectric layer and the ground plane is given by Eq. (15). For Au ground plane, that is
Where, is the reflection coefficient of the Au plate. For a reflective metasurface consisting of antennas with identical orientation, we can obtain the zero-order reflected electric fields by replacing the refractive index of the substrate with in Eq. (16). Figure 6(b) shows the analytical and simulated reflectivity spectra of co-polarization and cross polarization for normally incident RCP onto the three layer metasurface. It is obvious that the metasurface can work as a reflective half-wave plate with the wavelength from 1100nm to 1400nm.
Next, we consider a reflective metasurface shown in Fig. 7(a) with a linear phase gradient. By replacing with in Eq. (19), the analytical expressions of zero-order and first-order reflective electric fields from the metasurface can be obtained. Figure 7(b) shows the zero-order and first-order reflectivity obtained from both the analytical theory and CST simulation for a right circularly polarized light normally incident onto the phase-gradient metasurface. Again, the analytical results agree well with the simulated results.
4.3 One-dimensional phase modulated metasurface
In order to further validate the applicability of the analytical theory for more complex phase profiles, we apply it to a metasurface with a carefully designed nonlinear phase gradient for selectively exciting a number of diffraction orders. The metasurface unit cell consists of 13 antennas with orientation angles with respect to the x axis designed as [0.9π, 0.06π, 0.21π, 0.34π, 0.39π, 0.37π, 0.47π, 0.61π, 0.76π, 0.93π, 0.14π, 0.46π, 0.72π]. As shown in Fig. 8, three diffractive orders from −3 to −1 are selectively excited over the whole displayed wavelength range from 970nm to 1700nm. Meanwhile, the zeroth order is very low across the wavelength between 1200 nm and 1400 nm. Both numerical simulation and analytical modeling are in good agreement. All the calculation and simulation prove the accuracy of the analytical theory, and it can be applied to the modeling of more complicated metasurface holograms.
5. Conclusion
In summary, we have derived a general analytical theory of metasurfaces based on the assumptions that the tangential electric fields on both sides of the metasurface are continuous and the meta-atoms possess localized resonances with Lorentz-Drude forms. By retrieving the six parameters used to describe the anisotropic polarizability tensor of a meta-atom, optical performance of geometric metasurfaces constituted by the same meta-atoms but different orientations can be accurately calculated by the analytical formulae. As the analytical results show excellent agreement with full-wave numerical simulations, the analytical theory provides an efficient method to design and model optical devices based on metasurfaces of complex phase profiles. By introducing the magnetic polarizabilities and chirality into meta-atoms, the phenomenological modeling is possible to be extended to Huygens' metasurfaces and bianisotropic metasurfaces.
Acknowledgment
This work was supported by EPSRC (Grant No. EP/J018473/1), National Natural Science Foundation of China (NSFC) (Grant Nos. 11374367, 61328503), and Leverhulme Trust (Grant No. RPG-2012-674).
References and links
1. A. V. Kildishev, A. Boltasseva, and V. M. Shalaev, “Planar Photonics with Metasurfaces,” Science 339(6125), 1232009 (2013). [CrossRef] [PubMed]
2. N. Yu and F. Capasso, “Flat optics with designer metasurfaces,” Nat. Mater. 13(2), 139–150 (2014). [CrossRef] [PubMed]
3. N. Meinzer, W. L. Barnes, and I. R. Hooper, “Plasmonic meta-atoms and metasurfaces,” Nat. Photonics 8(12), 889–898 (2014). [CrossRef]
4. M. A. Kats, P. Genevet, G. Aoust, N. Yu, R. Blanchard, F. Aieta, Z. Gaburro, and F. Capasso, “Giant birefringence in optical antenna arrays with widely tailorable optical anisotropy,” Proc. Natl. Acad. Sci. U.S.A. 109(31), 12364–12368 (2012). [CrossRef]
5. S. Jiang, X. Xiong, Y. Hu, Y. Hu, G. Ma, R. Peng, C. Sun, and M. Wang, “Controlling the Polarization State of Light with a Dispersion-Free Metastructure,” Phys. Rev. X 4(2), 021026 (2014). [CrossRef]
6. B. Yang, W. M. Ye, X. D. Yuan, Z. H. Zhu, and C. Zeng, “Design of ultrathin plasmonic quarter-wave plate based on period coupling,” Opt. Lett. 38(5), 679–681 (2013). [CrossRef] [PubMed]
7. N. Yu, P. Genevet, M. A. Kats, F. Aieta, J. P. Tetienne, F. Capasso, and Z. Gaburro, “Light propagation with phase discontinuities: generalized laws of reflection and refraction,” Science 334(6054), 333–337 (2011). [CrossRef] [PubMed]
8. X. Ni, N. K. Emani, A. V. Kildishev, A. Boltasseva, and V. M. Shalaev, “Broadband Light Bending with Plasmonic Nanoantennas,” Science 335(6067), 427 (2012). [CrossRef] [PubMed]
9. S. Sun, Q. He, S. Xiao, Q. Xu, X. Li, and L. Zhou, “Gradient-index meta-surfaces as a bridge linking propagating waves and surface waves,” Nat. Mater. 11(5), 426–431 (2012). [CrossRef] [PubMed]
10. X. Ni, S. Ishii, A. V. Kildishev, and V. M. Shalaev, “Ultra-thin, planar, Babinet inverted plasmonic metalenses,” Light Sci. Appl. 2(4), e72 (2013). [CrossRef]
11. A. Pors, M. G. Nielsen, R. L. Eriksen, and S. I. Bozhevolnyi, “Broadband Focusing Flat Mirrors Based on Plasmonic Gradient Metasurfaces,” Nano Lett. 13(2), 829–834 (2013). [CrossRef] [PubMed]
12. P. Genevet, N. Yu, F. Aieta, J. Lin, M. A. Kats, R. Blanchard, M. O. Scully, Z. Gaburro, and F. Capasso, “Ultra-thin plasmonic optical vortex plate based on phase discontinuities,” Appl. Phys. Lett. 100(1), 013101 (2012). [CrossRef]
13. J. Lin, P. Genevet, M. A. Kats, N. Antoniou, and F. Capasso, “Nanostructured holograms for broadband manipulation of vector beams,” Nano Lett. 13(9), 4269–4274 (2013). [CrossRef] [PubMed]
14. X. Ni, A. V. Kildishev, and V. M. Shalaev, “Metasurface holograms for visible light,” Nat. Commun. 4, 2807 (2013). [CrossRef]
15. W. T. Chen, K. Y. Yang, C. M. Wang, Y. W. Huang, G. Sun, I. D. Chiang, C. Y. Liao, W. L. Hsu, H. T. Lin, S. Sun, L. Zhou, A. Q. Liu, and D. P. Tsai, “High-efficiency broadband meta-hologram with polarization-controlled dual images,” Nano Lett. 14(1), 225–230 (2014). [CrossRef] [PubMed]
16. G. Li, M. Kang, S. Chen, S. Zhang, E. Y. Pun, K. W. Cheah, and J. Li, “Spin-Enabled Plasmonic Metasurfaces for Manipulating Orbital Angular Momentum of Light,” Nano Lett. 13(9), 4148–4151 (2013). [CrossRef] [PubMed]
17. X. Yin, Z. Ye, J. Rho, Y. Wang, and X. Zhang, “Photonic spin hall effect at metasurfaces,” Science 339(6126), 1405–1407 (2013). [CrossRef] [PubMed]
18. N. Shitrit, I. Yulevich, E. Maguid, D. Ozeri, D. Veksler, V. Kleiner, and E. Hasman, “Spin-optical metamaterial route to spin-controlled photonics,” Science 340(6133), 724–726 (2013). [CrossRef] [PubMed]
19. E. Karimi, S. A. Schulz, I. De Leon, H. Qassim, J. Upham, and R. W. Boyd, “Generating optical orbital angular momentum at visible wavelengths using a plasmonic metasurface,” Light Sci. Appl. 3(5), e167 (2014). [CrossRef]
20. L. Huang, X. Chen, B. Bai, Q. Tan, G. Jin, T. Zentgraf, and S. Zhang, “Helicity dependent directional surface plasmon polariton excitation using a metasurface with interfacial phase discontinuity,” Light Sci. Appl. 2(3), e70 (2013). [CrossRef]
21. J. Lin, J. P. Mueller, Q. Wang, G. Yuan, N. Antoniou, X. C. Yuan, and F. Capasso, “Polarization-controlled tunable directional coupling of surface plasmon polaritons,” Science 340(6130), 331–334 (2013). [CrossRef] [PubMed]
22. L. Huang, X. Chen, H. Mühlenbernd, H. Zhang, S. Chen, B. Bai, Q. Tan, G. Jin, K. Cheah, C. Qiu, J. Li, T. Zentgraf, and S. Zhang, “Three-dimensional optical holography using a plasmonic metasurface,” Nat. Commun. 4, 2808 (2013). [CrossRef]
23. G. Zheng, H. Mühlenbernd, M. Kenney, G. Li, T. Zentgraf, and S. Zhang, “Metasurface holograms reaching 80% efficiency,” Nat. Nanotechnol. 10(4), 308–312 (2015). [CrossRef] [PubMed]
24. Y. Zhao, N. Engheta, and A. Alù, “Homogenization of plasmonic metasurfaces modeled as transmission-line loads,” Metamaterials (Amst.) 5(2-3), 90–96 (2011). [CrossRef]
25. Y. Zhao and A. Alù, “Manipulating light polarization with ultrathin plasmonic metasurfaces,” Phys. Rev. B 84(20), 205428 (2011). [CrossRef]
26. C. Pfeiffer, C. Zhang, V. Ray, L. J. Guo, and A. Grbic, “High performance bianisotropic metasurfaces: Asymmetric transmission of light,” Phys. Rev. Lett. 113(2), 023902 (2014). [CrossRef] [PubMed]
27. A. Epstein and G. V. Eleftheriades, “Floquet-Bloch analysis of refracting Huygens metasurfaces,” Phys. Rev. B 90(23), 235127 (2014). [CrossRef]
28. A. Epstein and G. V. Eleftheriades, “Huygens’ metasurfaces via the equivalence principle: design and applications,” J. Opt. Soc. Am. B 33(2), A31 (2016). [CrossRef]
29. E. F. Kuester, M. A. Mohamed, M. Piket-May, and C. L. Holloway, “Averaged Transition Conditions for Electromagnetic Fields at a Metafilm,” IEEE Trans. Antenn. Propag. 51(10), 2641–2651 (2003). [CrossRef]
30. D. J. Griffiths, Introduction to Electrodynamics, 3rd ed. (PrenticeHall International, New Jersey, 1999)
31. A. M. Patel and A. Grbic, “A Printed Leaky-Wave Antenna Based on a Sinusoidally-Modulated Reactance Surface,” IEEE Trans. Antenn. Propag. 59(6), 2087–2096 (2011). [CrossRef]