pith. sign in

arxiv: 1907.03754 · v1 · pith:GIEB5OCSnew · submitted 2019-07-08 · ⚛️ physics.optics

Spectral Numerical Mode Matching Method for 3D Layered Multi-Region Structures

Jie Liu , Na Liu , Qing Huo Liu This is my paper

Pith reviewed 2026-05-25 01:02 UTC · model grok-4.3

classification ⚛️ physics.optics
keywords spectral numerical mode-matching3D layered structuresmetasurfaceslithography modelselectromagnetic simulationreflection and transmission matriceseigenmode expansionsdimensionality reduction
0
0 comments X

The pith

The spectral numerical mode-matching method extends to multiple layers by using reflection and transmission matrices from transverse eigenmode expansions to treat electromagnetic propagation analytically.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper develops an extension of the spectral numerical mode-matching method for simulating electromagnetic fields in three-dimensional structures built from multiple layers. The approach expands fields in terms of eigenmodes in the directions across each layer and then handles the stacking of layers through matrix multiplications that capture reflections and transmissions at each interface. A sympathetic reader would care because this turns a full three-dimensional computation into a collection of simpler two-dimensional eigenvalue problems, which can be much faster when the layers are thick compared with the wavelength. The method is demonstrated on metasurfaces and lithography models where thin conducting or meta surfaces sit at the layer boundaries.

Core claim

The 3D SNMM method is extended from a single interface to multiple layers so that the electromagnetic propagation and scattering in the longitudinal direction is treated analytically through reflection and transmission matrices by using the eigenmode expansions in the transverse directions. Therefore, it effectively reduces the original 3D problem into a series of 2D eigenvalue problems for periodic structures.

What carries the argument

Reflection and transmission matrices obtained from eigenmode expansions in the transverse directions, which analytically handle propagation and scattering between layers.

If this is right

  • The method is especially efficient when the longitudinal layer thicknesses are large compared with wavelength.
  • It accurately characterizes metasurfaces and lithography models.
  • It handles thin surfaces such as good conductor surfaces and metasurfaces deposited at layer interfaces.
  • The semi-analytical nature provides dimensionality reduction to lower computational costs for microwave and optical integrated circuits.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the transverse eigenmodes can be computed efficiently for periodic structures, the overall scaling improves for problems with many thick layers.
  • Similar matrix-based stacking could be applied to other wave problems beyond electromagnetics if eigenmode bases are available.
  • Validation against full 3D solvers would be needed when layer thicknesses approach the wavelength scale.

Load-bearing premise

The structures consist of discrete layers with well-defined interfaces where eigenmode expansions in the transverse directions remain valid and sufficient even when thin surfaces are present.

What would settle it

A direct comparison of SNMM results against a known analytical solution or a converged full-wave 3D simulation for a multi-layer structure with layer thicknesses much larger than the wavelength, checking if the computed fields or scattering parameters match within numerical tolerance.

Figures

Figures reproduced from arXiv: 1907.03754 by Jie Liu, Na Liu, Qing Huo Liu.

Figure 1
Figure 1. Figure 1: The geometry of a multi-region layered doubly periodic structure. [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 1
Figure 1. Figure 1: As shown in Figure 1, we assume that the direction [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Sketch of the air/graphene/dielectric/metal/air (AGDMA) structure for [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Electric field component Ey along the z-axis at (x, y) = (0, 0) for the AGDGMA structure. (a) Real part of Ey. (b) Imaginary part of Ey. structures filled with different media (i.e., different varia￾tions of [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Schematic drawing for the gradient metasurface in the optical [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: Optical frequency case in Figure 4: the field distribution of [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Optical frequency case in Figure 4: the electric field component [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figure 9
Figure 9. Figure 9: Schematic drawing for the EUV lithography. [PITH_FULL_IMAGE:figures/full_fig_p008_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Electric field component Ex along the z-axis at (x, y) = (0, 0) for the lithography model with each layer filled with the homogeneous medium. (a) Real part of Ex. (b) Imaginary part of Ex. 0 50 100 150 200 250 300 350 −2 −1 0 1 2 Re(E x) (V/m) (a) FNMM SNMM 0 50 100 150 200 250 300 350 −2 −1 0 1 2 Im(E x) (V/m) z (nm) (b) FNMM SNMM [PITH_FULL_IMAGE:figures/full_fig_p008_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Electric field component Ex along the z-axis at (x, y) = (0, 0) for the lithography model with patterns DU. (a) Real part of Ex. (b) Imaginary part of Ex. incident to this structure and the direction of polarization is along yˆ, i.e., (θk, φk, φe) = (0, 0, π/2). In COMSOL, the perfect matching layer (PML) absorbing boundary condition is used to truncate the first and the last layers, so that the thickness… view at source ↗
Figure 12
Figure 12. Figure 12: Field distribution of the reflection field [PITH_FULL_IMAGE:figures/full_fig_p009_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Field distribution of the reflection field [PITH_FULL_IMAGE:figures/full_fig_p009_13.png] view at source ↗
read the original abstract

The spectral numerical mode-matching (SNMM) method is developed to simulate the 3D layered multi-region structures. The SNMM method is a semi-analytical solver having the properties of dimensionality reduction to reduce computational costs; it is especially useful for microwave and optical integrated circuits where fabrication is often done in a layered structure. Furthermore, at some layer interfaces, very thin surfaces such as good conductor surfaces and metasurfaces can be deposited to achieve desired properties such as high absorbance and/or anomalous reflection/refraction. In this work, the 3D SNMM method is further extended from a single interface to multiple layers so that the electromagnetic propagation and scattering in the longitudinal direction is treated analytically through reflection and transmission matrices by using the eigenmode expansions in the transverse directions. Therefore, it effectively reduces the original 3D problem into a series of 2D eigenvalue problems for periodic structures. We apply this method to characterize metasurfaces and lithography models, and show that the SNMM method is especially efficient when the longitudinal layer thicknesses are large compared with wavelength. Numerical experiments indicate that the SNMM method is highly efficient and accurate for the metasurfaces and the lithography models.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The paper extends the spectral numerical mode-matching (SNMM) method from single to multiple interfaces in 3D layered structures. Longitudinal propagation and scattering are treated analytically via reflection and transmission matrices constructed from transverse eigenmode expansions, reducing the original 3D problem to a series of 2D eigenvalue problems for periodic structures. The method is applied to metasurfaces and lithography models, with claims of high efficiency and accuracy when longitudinal layer thicknesses are large compared to wavelength.

Significance. If the extension is rigorously validated, the analytic treatment of longitudinal scattering via R/T matrices offers a genuine dimensionality reduction that could be advantageous for layered electromagnetic problems with large layer thicknesses. The applicability to thin metasurfaces would be a notable feature if the interface matching remains complete without additional terms.

major comments (2)
  1. [Abstract and method-extension description] Abstract and method-extension paragraph: the central claim that standard eigenmode expansions and R/T matrices suffice for 'very thin surfaces such as good conductor surfaces and metasurfaces' is load-bearing for the stated applicability. For infinitely thin interfaces, tangential-field discontinuities due to surface currents typically require explicit incorporation into the mode-matching conditions or auxiliary unknowns; no such modification or completeness proof is referenced.
  2. [Numerical experiments and applications] Application sections (metasurface characterization and lithography models): the assertion that 'numerical experiments indicate that the SNMM method is highly efficient and accurate' lacks cited quantitative support (error norms, convergence rates, or comparisons) in the provided description. Without these, the efficiency claim for large layer thicknesses cannot be assessed as load-bearing evidence.
minor comments (1)
  1. [Abstract] The abstract would be strengthened by including at least one concrete performance metric (e.g., CPU time or error level) to support the efficiency statements.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments on the method extension and the need for stronger validation of the claims. We respond to each major comment below.

read point-by-point responses

  1. Referee: [Abstract and method-extension description] Abstract and method-extension paragraph: the central claim that standard eigenmode expansions and R/T matrices suffice for 'very thin surfaces such as good conductor surfaces and metasurfaces' is load-bearing for the stated applicability. For infinitely thin interfaces, tangential-field discontinuities due to surface currents typically require explicit incorporation into the mode-matching conditions or auxiliary unknowns; no such modification or completeness proof is referenced.

    Authors: The transverse eigenmode expansions in SNMM are formulated to satisfy the interface boundary conditions, which for thin metasurfaces and conductors incorporate the tangential-field discontinuities arising from surface currents directly into the 2D eigenvalue problem at each layer. The resulting R/T matrices then carry these modes forward analytically. We will add a clarifying paragraph in the method section describing these interface conditions and noting the completeness of the transverse basis under standard periodic boundary assumptions, together with a reference to the single-interface formulation. revision: yes

  2. Referee: [Numerical experiments and applications] Application sections (metasurface characterization and lithography models): the assertion that 'numerical experiments indicate that the SNMM method is highly efficient and accurate' lacks cited quantitative support (error norms, convergence rates, or comparisons) in the provided description. Without these, the efficiency claim for large layer thicknesses cannot be assessed as load-bearing evidence.

    Authors: The body of the manuscript contains quantitative results in the application sections, including relative error norms, mode-convergence rates, and timing comparisons versus reference solvers that demonstrate the efficiency gain for thick layers. We will revise the abstract to cite these specific metrics and add explicit cross-references to the relevant figures and tables. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation chain is self-contained

full rationale

The paper describes an extension of the SNMM method to multiple layers by constructing reflection/transmission matrices from transverse eigenmode expansions, thereby reducing the 3D problem to 2D eigenvalue problems per layer. No equations or steps in the provided abstract reduce a claimed prediction or result to a fitted parameter, self-citation, or definitional tautology. The central claim rests on the standard validity of eigenmode expansions and interface matching for layered structures (including thin surfaces), which is presented as an independent modeling choice rather than derived from the paper's own outputs. No load-bearing self-citation, ansatz smuggling, or renaming of known results is evident. The method is applied to metasurfaces and lithography models with reported efficiency gains, but these are numerical validations, not circular reductions. This is the expected outcome for a semi-analytical extension paper whose assumptions are stated externally to the derivation.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The method rests on standard electromagnetic theory and numerical techniques for layered media; no free parameters, new entities, or ad-hoc axioms are identifiable from the abstract alone.

axioms (1)
  • domain assumption Eigenmode expansions in transverse directions combined with reflection/transmission matrices accurately capture longitudinal propagation and scattering in layered periodic structures.
    Invoked when the abstract states that propagation is treated analytically through reflection and transmission matrices using eigenmode expansions.

pith-pipeline@v0.9.0 · 5734 in / 1285 out tokens · 23918 ms · 2026-05-25T01:02:34.000662+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

40 extracted references · 40 canonical work pages

  1. [1]

    J. R. Wait, Electromagnetic Wave In Stratified Media. New York: Oxford Univ. Press, 1970

    work page 1970

  2. [2]

    Propagation of electromagnetic waves through geological beds in a geophysical probing environment,

    W. C. Chew and B. Anderson, “Propagation of electromagnetic waves through geological beds in a geophysical probing environment,” Radio Sci., vol. 20, no. 3, pp. 611–621, 1985

    work page 1985

  3. [3]

    A full wave analysis of wave propagation in multiregion, cylindrically stratified media,

    W. C. Chew, Z. P. Nie, Q. H. Liu, and B. Anderson, “A full wave analysis of wave propagation in multiregion, cylindrically stratified media,” in Proc. Int. Symp. Antennas Propag. Soc. , vol. 2, pp. 602–605, 1990

    work page 1990

  4. [4]

    Numerical mode-matching method for the multiregion vertically stratified media,

    Q. H. Liu and W. C. Chew, “Numerical mode-matching method for the multiregion vertically stratified media,” IEEE Trans. Antennas Propag. , vol. 38, no. 4, pp. 498–506, Apr. 1990

    work page 1990

  5. [5]

    Efficient computation of electromagnetic wave in anisotropic orthogonal-plano-cylindrically layered media using the improved numerical mode matching (NMM) method,

    J. W. Dai and Q. H. Liu, “Efficient computation of electromagnetic wave in anisotropic orthogonal-plano-cylindrically layered media using the improved numerical mode matching (NMM) method,” IEEE Trans. Antennas Propag., vol. 63, no. 8, pp. 3569–3578, Aug. 2015

    work page 2015

  6. [6]

    Three-dimensionally periodic dielectric layered structure with omnidirectional photonic band gap,

    S. G. Johnson and J. D. Joannopoulos, “Three-dimensionally periodic dielectric layered structure with omnidirectional photonic band gap,” Appl. Phys. Lett. , vol. 38, no. 4, pp. 498–506, Apr. 1990

    work page 1990

  7. [7]

    A semianalitical spectral element method for the analysis of 3-D layered structures,

    J. F. Chen, B. Zhu, W. X. Zhong, and Q. H. Liu, “A semianalitical spectral element method for the analysis of 3-D layered structures,” IEEE Trans. Microw. Theory Tech., vol. 59, no. 1, pp. 1–8, Jan. 2011

    work page 2011

  8. [8]

    A layered finite element method for electromagnetic analysis of large-scale high-frequency integrated circuits,

    D. Jiao, S. Chakravarty, and C. H. Dai, “A layered finite element method for electromagnetic analysis of large-scale high-frequency integrated circuits,” IEEE Trans. Antennas Propag. , vol. 55, no. 2, pp. 422–432, Feb. 2007

    work page 2007

  9. [9]

    A time-domain layered finite element reduction recovery (LAFE-RR) method for high-frequency VLSI design,

    H. Gan and D. Jiao, “A time-domain layered finite element reduction recovery (LAFE-RR) method for high-frequency VLSI design,” IEEE Trans. Antennas Propag., vol. 55, no. 12, pp. 3620–3629, Dec. 2007

    work page 2007

  10. [10]

    A recovery algorithm of linear complexity in the time-domain layered finite element reduction recovery (LAFE-RR) method for large-scale electromagnetic analysis of high-speed ICs,

    H. Gan and D. Jiao, “A recovery algorithm of linear complexity in the time-domain layered finite element reduction recovery (LAFE-RR) method for large-scale electromagnetic analysis of high-speed ICs,” IEEE Trans. Adv. Packag., vol. 31, no. 3, pp. 612–618, Aug. 2008

    work page 2008

  11. [11]

    Efficient and reliable simulation of multicomponent induction logging response in horizontally stratified inhomogeneous TI formations by numerical mode matching method,

    H. Wang, H. Tao, J. Yao, and Y . Zhang, “Efficient and reliable simulation of multicomponent induction logging response in horizontally stratified inhomogeneous TI formations by numerical mode matching method,” IEEE Trans. Geosci. Remote Sens. , vol. 50, no. 9, pp. 3383–3395, Sep. 2012

    work page 2012

  12. [12]

    Manipulating light absorption of graphene using plasmonic nanoparticles,

    J. F. Zhu, Q. H. Liu, and T. Lin, “Manipulating light absorption of graphene using plasmonic nanoparticles,” Nanoscale, vol. 5, no. 17, pp. 7785-7789, Jul. 2013

    work page 2013

  13. [13]

    Near unity ultraviolet absorption in graphene without patterning,

    J. F. Zhu, S. Yan, N. F. Feng, L. F. Ye, J-Y Ou, and Q. H. Liu, “Near unity ultraviolet absorption in graphene without patterning,” Appl. Phys. Lett., vol. 112, no. 15, 153106, Apr. 2018

    work page 2018

  14. [14]

    J. F. Zhu, C. W. Li, J-Y , Ou, and Q. H. Liu, “Perfect light absorption in graphene by two unpatterned dielectric layers and potential applica- tions,“ Carbon, vol. 142, pp. 430-437, 2019

    work page 2019

  15. [15]

    Simulation of Metasurfaces in finite difference techniques,

    Y . Vahabzadeh, K. Achouri, and C. Caloz, “Simulation of Metasurfaces in finite difference techniques,” IEEE Trans. Antennas Propag., vol. 64, no. 11, pp. 4753-4759, Nov. 2016

    work page 2016

  16. [16]

    High-efficiency broadband cross polarization converter for near-infrared light based on anisotropic plasmonic Meta-surfaces,

    Z. Y . Song, L. Zhang, and Q. H. Liu, “High-efficiency broadband cross polarization converter for near-infrared light based on anisotropic plasmonic Meta-surfaces,” Plasmonics, vol. 11, no. 1, pp. 61-64, Feb. 2016

    work page 2016

  17. [17]

    Spectral element boundary integral method with periodic layered medium dyadic Green’s function for multiscale nano-optical scattering analysis,

    J. Niu, Y . Ren, and Q. H. Liu, “Spectral element boundary integral method with periodic layered medium dyadic Green’s function for multiscale nano-optical scattering analysis,” Opt. Express, vol. 25, no. 20, pp. 24199–24214, Oct. 2017

    work page 2017

  18. [18]

    W. C. Chew, S. Barone, B. Anderson, and C. Hennessy, Diffraction of axisymmetric waves in a borehole by bed boundary discontinuities, Geophysics, vol. 49, no. 10, pp. 1586-1595, 1984

    work page 1984

  19. [19]

    W. C. Chew, Z. P. Nie, Q. H. Liu, and B. Anderson, An efficient solution for the response of electrical well logging tools in a complex environment, IEEE Trans. Geosci. Remote Sens. , vol. 29, no. 2, pp. 308-313, May. 1991

    work page 1991

  20. [20]

    Diffraction of nonaxisymmetric waves in cylindrically layered media by horizontal discontinuities,

    Q. H. Liu and W. C. Chew, “Diffraction of nonaxisymmetric waves in cylindrically layered media by horizontal discontinuities,” Radio Sci. , vol. 27, no. 5, pp. 569–581, Apr. 1992

    work page 1992

  21. [21]

    Electromagnetic field generated by an off-axis source in a cylindrically layered medium with an arbitrary number of horizontal discontinuities,

    Q. H. Liu , “Electromagnetic field generated by an off-axis source in a cylindrically layered medium with an arbitrary number of horizontal discontinuities,” Geophysics, vol. 58, no. 5, pp. 616–625, May. 1993

    work page 1993

  22. [22]

    W. C. Chew, Response of a source on top of a vertically stratified half- space, IEEE Trans. Antennas Propag., vol. AP-33, no. 6, pp. 649C654, Jun. 1985

    work page 1985

  23. [23]

    Numerical mode matching method with perfectly matching layer,

    A. Li, Z. Nie, and Y . Zhao, “Numerical mode matching method with perfectly matching layer,” IEEE Antennas Propag. Soc. Int. Symp. , vol. 4B, pp. 372–375, Jul. 2005

    work page 2005

  24. [24]

    Numerical mode-matching method for tilted-coil antennas in cylindrically layered anisotropic media with multiple horizontal beds,

    Y .-K. Hue and F. L. Teixeira, “Numerical mode-matching method for tilted-coil antennas in cylindrically layered anisotropic media with multiple horizontal beds,” IEEE Trans. Geosci. Remote Sens. , vol. 45, no. 8, pp. 2451-2462, Aug. 2007

    work page 2007

  25. [25]

    Mixed spectral element method for the waveguide problem with Bloch periodic boundary conditions,

    J. Liu, W. Jiang, N. Liu, and Q. H. Liu, “Mixed spectral element method for the waveguide problem with Bloch periodic boundary conditions,” IEEE Trans. Electromagn. Compat., doi: 10.1109/TEMC.2018.2866023, 2018

    work page doi:10.1109/temc.2018.2866023 2018

  26. [26]

    A 3-D spectral-element method using mixed-order curl conforming vector basis functions for electromagnetic fields,

    J. H. Lee, T. Xiao, and Q. H. Liu, “A 3-D spectral-element method using mixed-order curl conforming vector basis functions for electromagnetic fields,”IEEE Trans. Microw. Theory Tech., vol. 54, No. 1, pp. 437-444, Jan. 2006

    work page 2006

  27. [27]

    Spectral element method for band structures of two-dimensional anisotropic photonic crystals,

    M. Luo, Q. H. Liu, and Z. Li, “Spectral element method for band structures of two-dimensional anisotropic photonic crystals,” Phys. Rev. E, vol. 79, 026705, 2009

    work page 2009

  28. [28]

    Spectral element method for band structures of three-dimensional anisotropic photonic crystals,

    M. Luo and Q. H. Liu, “Spectral element method for band structures of three-dimensional anisotropic photonic crystals,” Phys. Rev. E , vol. 80, no. 5, 056702, Nov. 2009

    work page 2009

  29. [29]

    Mixed spectral-element method for 3-D Maxwell’s eigenvalue problem,

    N. Liu, L. E. Tob ´on, Y . M. Zhao, Y . F. Tang, and Q. H. Liu, “Mixed spectral-element method for 3-D Maxwell’s eigenvalue problem,” IEEE Trans. Microw. Theory Tech., vol. 63, no. 2, pp. 3094–3102, Feb. 2015

    work page 2015

  30. [30]

    The mixed spectral-element method for anisotropic, lossy, and open waveguides,

    N. Liu, G. X. Cai, C. H. Zhu, Y . F. Tang, and Q. H. Liu, “The mixed spectral-element method for anisotropic, lossy, and open waveguides,” IEEE Trans. Microw. Theory Tech., vol. 63, no. 10, pp. 3094–3102, Oct. 2015

    work page 2015

  31. [31]

    A spectral-element method for the analysis of 2-D waveguide devices with sharp edges and irregular shapes,

    O. A. Peverini, G. Addamo, G. Virone, R. Tascone, and R. Orta, “A spectral-element method for the analysis of 2-D waveguide devices with sharp edges and irregular shapes,” IEEE Trans. Microw. Theory Tech. , vol. 59, no. 7, pp. 1685–1695, July 2011

    work page 2011

  32. [32]

    3-D numerical mode- matching (NMM) method for resistivity well logging tools,

    G. X. Fan, Q. H. Liu, and S. P. Blanchard, “3-D numerical mode- matching (NMM) method for resistivity well logging tools,”IEEE Trans. Antennas Propag., vol. 48, no. 10, pp. 1544-1552, Oct. 2000

    work page 2000

  33. [33]

    Spectral Numerical Mode Matching Method for Metasurfaces,

    J. Liu, G. X. Cai, J. Yao, N. Liu, and Q. H. Liu, “Spectral Numerical Mode Matching Method for Metasurfaces,” IEEE Trans. Microw. Theory Tech., accepted, 2019

    work page 2019

  34. [34]

    ¨Uber die Quantenmechanik der Elektronen in Kristallgittern,

    F. Bloch, “ ¨Uber die Quantenmechanik der Elektronen in Kristallgittern,” Zeitschrift f ¨ur Physik, vol. 52, no. 7, pp. 555–600, Jul. 1929

    work page 1929

  35. [35]

    A variational analysis of anisotropic inhomogeneous dielectric waveguides,

    W. C. Chew and M. A. Nasir, “A variational analysis of anisotropic inhomogeneous dielectric waveguides,” IEEE Trans. Microw. Theory Tech., vol. 37, no. 4, pp. 661-668, Apr. 1989

    work page 1989

  36. [36]

    Fine structure constant defines visual transparency of graphene,

    R. R. Nair, P. Blake, A. N. Grigorenko, K. S. Novoselov, T. J. Booth, T. Stauber, N. M. R. Peres, and A. K. Geim, “Fine structure constant defines visual transparency of graphene,” Science, vol. 320, no. 5881, pp. 1308, Jun. 2008

    work page 2008

  37. [37]

    Strongly confined gap plasmon modes in graphene sandwiches and graphene-on-silicon,

    Y . Francescato, V . Giannini, and S. A. Maier, “Strongly confined gap plasmon modes in graphene sandwiches and graphene-on-silicon,” New J. Phys., vol. 15, 063020, Jun. 2013

    work page 2013

  38. [38]

    Light propagation with phase discontinuities: generalized laws of reflection and refaction,

    N. F. 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 refaction,” Science, vol. 334, pp. 333-337, Oct. 2011

    work page 2011

  39. [39]

    Gradient- index meta-surface as a bridge linking propagating waves and surface waves,

    S. L. Sun, Q. He, S. Y . Xiao, Q. Xu, X. Li, and L. Zhou, “Gradient- index meta-surface as a bridge linking propagating waves and surface waves,” Nature Mater., vol. 11, no. 5, pp. 426–431, May. 2012

    work page 2012

  40. [40]

    Q. H. Liu, Classnotes for ECE 571 Electromagnetics Theory , Duke University, Durham, North Carolina, Fall 2018

    work page 2018