Twistoptics in Planar Heterostructures with an Arbitrary Number of Rotated 3D Thin Layers and 2D Conductive Sheets
Pith reviewed 2026-05-10 12:37 UTC · model grok-4.3
The pith
A general analytical model now exists for polaritons in twisted stacks with any number of layers.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We present an analytical model for planar stacks comprising an arbitrary number of finite-thickness anisotropic (biaxial) layers and infinitesimally thin anisotropic conductive sheets. The formalism and its high-momentum and thin-film approximations predict key polaritonic observables, such as wavelength, propagation length, and electromagnetic field distributions.
What carries the argument
Analytical formalism for calculating the dispersion relation in twisted anisotropic multilayer stacks using permittivity tensors rotated by the twist angles and boundary conditions at layer interfaces.
If this is right
- Polariton wavelength and propagation length become calculable for stacks with any number of layers and twists.
- Electromagnetic field distributions can be determined throughout the structure.
- High-momentum approximation yields simple expressions for deeply confined polaritons.
- Thin-film approximation applies when layers are much thinner than the wavelength.
- Open-access scripts enable immediate application of the model.
Where Pith is reading between the lines
- Designers could now systematically vary twist angles across multiple layers to achieve desired canalization or focusing effects without repeated simulations.
- The approach may extend to time-varying or nonlinear responses in twisted structures.
- Connections to twistronics in electronics suggest similar analytical tools could unify descriptions across frequency ranges.
Load-bearing premise
The response of each material layer is accurately represented by a frequency-dependent local anisotropic permittivity tensor, with standard electromagnetic boundary conditions holding at all interfaces.
What would settle it
Measuring the actual polariton wavelength in a fabricated three-layer twisted heterostructure and checking if it matches the model's prediction for the given twist angles and thicknesses.
Figures
read the original abstract
Twistoptics has recently emerged as a branch of nano-optics that explores light propagation in stacks of thin anisotropic layers rotated relative to one another. The concept is particularly relevant for polaritons -- hybrid light-matter quasiparticles -- in van der Waals (vdW) materials, where strong in-plane anisotropy and deep subwavelength confinement make the polaritonic dispersion highly sensitive to interlayer twist angles. This sensitivity enables exotic phenomena such as canalization, i.e., diffraction-free propagation, with potential applications ranging from thermal management to super-resolution imaging. Despite rapid progress, a general analytical framework to describe polariton propagation in twisted planar heterostructures has been missing. Here we present an analytical model for planar stacks comprising an arbitrary number of finite-thickness anisotropic (biaxial) layers and infinitesimally thin anisotropic conductive sheets. The formalism and its high-momentum and thin-film approximations predict key polaritonic observables, such as wavelength, propagation length, and electromagnetic field distributions. We also provide open-access numerical scripts implementing the model to support their practical use. Together, these results provide a general theoretical foundation for twistoptics and should facilitate the discovery and accelerate the implementation of twist-engineered polaritonic phenomena across the electromagnetic spectrum.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents an analytical transfer-matrix formalism for electromagnetic polariton modes in planar heterostructures containing an arbitrary number of finite-thickness rotated biaxial layers together with infinitesimally thin anisotropic conductive sheets. The model is derived from Maxwell's equations with appropriate boundary conditions, supplies closed-form high-momentum and thin-film approximations, and predicts dispersion relations, propagation lengths, and field distributions; open-access numerical scripts implementing the formalism are provided.
Significance. If the central derivations are rigorous, the work supplies a missing general analytical tool for twistoptics that extends beyond existing treatments limited to few layers or isotropic cases. The combination of arbitrary-stack capability, explicit approximations, and reproducible code would enable systematic exploration of canalization and related phenomena across vdW material stacks.
major comments (1)
- [Section on approximations (high-momentum and thin-film limits)] The high-momentum and thin-film approximations are central to the predictive claims, yet their validity domains are stated only qualitatively (accuracy depends on parameters not being too extreme). Explicit bounds in terms of twist angle, layer thickness relative to wavelength, and permittivity anisotropy should be derived and tested against the exact transfer-matrix result.
minor comments (2)
- [Numerical results / validation] The manuscript would benefit from at least one explicit benchmark against a known analytic limit (e.g., zero-twist isotropic stack or single-layer hyperbolic polariton dispersion) with quantitative error metrics.
- [Formalism section] Notation for the rotation matrices and the ordering of layers in the transfer-matrix product should be defined once in a dedicated subsection to avoid ambiguity when the number of layers is arbitrary.
Simulated Author's Rebuttal
We thank the referee for the positive summary and recommendation for minor revision. The single major comment is addressed below; we agree it improves the manuscript and will revise accordingly.
read point-by-point responses
-
Referee: [Section on approximations (high-momentum and thin-film limits)] The high-momentum and thin-film approximations are central to the predictive claims, yet their validity domains are stated only qualitatively (accuracy depends on parameters not being too extreme). Explicit bounds in terms of twist angle, layer thickness relative to wavelength, and permittivity anisotropy should be derived and tested against the exact transfer-matrix result.
Authors: We agree that the current qualitative statements on validity can be strengthened. In the revised manuscript we will add explicit bounds derived from the underlying dispersion relations. For the high-momentum (quasistatic) limit we will state the condition |k| ≫ 2π/λ ⋅ max(√|ε_xx|, √|ε_yy|, √|ε_zz|) together with a quantitative error metric (relative deviation in Re(k) and Im(k) < 5 %). For the thin-film limit we will give d/λ < 0.05 ⋅ min(1/√|ε_⊥|, 1/√|ε_∥|) with the same error threshold. These bounds will be obtained analytically from the small-argument expansion of the transfer-matrix elements and then validated numerically against the exact transfer-matrix result for representative cases: twist angles 0°–90°, d/λ ∈ [0.001, 0.2], and anisotropy ratios |ε_xx/ε_yy| ∈ [1, 20]. The new subsection will include two validation figures (one for each approximation) and a short table summarizing the maximum parameter ranges where the stated error is respected. The open-access scripts will be updated to include an optional “validate_approximation” flag that reproduces these comparisons. revision: yes
Circularity Check
No significant circularity; derivation from standard EM boundary conditions
full rationale
The paper presents a transfer-matrix formalism for arbitrary stacks of rotated biaxial layers and conductive sheets, derived from Maxwell's equations and interface boundary conditions. No load-bearing steps reduce to self-definitions, fitted parameters renamed as predictions, or self-citation chains. The high-momentum and thin-film approximations are standard limits applied to the general solution. Open scripts enable direct verification against isotropic/zero-twist cases, confirming the model is self-contained and falsifiable externally rather than tautological.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Local anisotropic permittivity tensor description of each layer
- standard math Validity of standard electromagnetic boundary conditions at each interface
Reference graph
Works this paper leans on
-
[1]
A. K. Geim and I. V Grigorieva, Van der Waals heterostructures, Nature 499, 419 (2013)
work page 2013
-
[2]
K. S. Novoselov, A. Mishchenko, A. Carvalho, and A. H. Castro Neto, 2D materials and van der Waals heterostructures, Science (1979) 353, aac9439 (2016)
work page 1979
-
[3]
S. Carr, D. Massatt, S. Fang, P. Cazeaux, M. Luskin, and E. Kaxiras, Twistronics: Manipulating the electronic properties of two -dimensional layered structures through their twist angle, Phys Rev B 95, (2017)
work page 2017
-
[4]
Y . Cao, V . Fatemi, S. Fang, K. Watanabe, T. Taniguchi, E. Kaxiras, and P. Jarillo -Herrero, Unconventional superconductivity in magic-angle graphene superlattices, Nature 556, 43 (2018)
work page 2018
-
[5]
Y . Cao et al., Correlated insulator behaviour at half-filling in magic-angle graphene superlattices, Nature 556, 80 (2018)
work page 2018
-
[6]
M. Yankowitz, S. Chen, H. Polshyn , Y . Zhang, K. Watanabe, T. Taniguchi, D. Graf, A. F. Young, and C. R. Dean, Tuning superconductivity in twisted bilayer graphene, Science (1979) 363, 1059 (2019)
work page 1979
-
[7]
Y . Cao, D. Chowdhury, D. Rodan -Legrain, O. Rubies -Bigorda, K. Watanabe, T. Taniguchi, T. Senthil, and P. Jarillo -Herrero, Strange Metal in Magic -Angle Graphene with near Planckian Dissipation, Phys Rev Lett 124, (2020)
work page 2020
-
[8]
G. Hu, A. Krasnok, Y. Mazor, C.-W. Qiu, and A. Alù. Moirè Hyperbolic Metasurfaces. Nano Lett. 2020, 20, 5, 3217–3224
work page 2020
-
[9]
J. Álvarez-Cuervo, M. Obst, S. Dixit et al. Unidirectional ray polaritons in twisted asymmetric stacks. Nat Commun 15, 9042 (2024)
work page 2024
-
[10]
R. Hillenbrand, Y . Abate, M. Liu, X. Chen and D. N. Basov . Visible-to-THz near -field nanoscopy. Nat Rev Mater 10, 285–310 (2025)
work page 2025
-
[11]
S. Dai et al., Tunable Phonon Polaritons in Atomically Thin van der Waals Crystals of Boron Nitride, Science (1979) 343, 1125 (2014)
work page 1979
-
[12]
S. Dai et al., Subdiffractional focusing and guiding of polaritonic rays in a natural hyperbolic material, Nat Commun 6, 6963 (2015)
work page 2015
-
[13]
S. Dai et al., Efficiency of Launching Highly Confined Polaritons by Infrared Light Incident on a Hyperbolic Material, Nano Lett 17, 5285 (2017)
work page 2017
-
[14]
A. J. Giles et al., Ultralow -loss polaritons in isotopically pure boron nitride, Nat Mater 17, 134 (2018). 17
work page 2018
-
[15]
W. Ma et al., In-plane anisotropic and ultra-low-loss polaritons in a natural van der Waals crystal, Nature 562, 557 (2018)
work page 2018
-
[16]
J. Taboada-Gutiérrez et al., Broad spectral tuning of ultra -low-loss polaritons in a van der Waals crystal by intercalation, Nat Mater 19, 964 (2020)
work page 2020
-
[17]
D. N. Basov, M. M. Fogler, and F. J. García de Abajo, Polaritons in van der Waals materials, Science (1979) 354, aag1992 (2016)
work page 1979
-
[18]
T. Low, A. Chaves, J. D. Caldwell, A. Kumar, N. X. Fang, P. Avouris, T. F. Heinz, F. Guinea, L. Martin-Moreno, and F. Koppens, Polaritons in layered two -dimensional materials, Nat Mater 16, 182 (2017)
work page 2017
-
[19]
P. Li et al., Collective near -field coupling and nonlocal phenomena in infrared -phononic metasurfaces for nano-light canalization, Nat Commun 11, 3663 (2020)
work page 2020
-
[20]
A. I. F. Tresguerres-Mata et al., Observation of naturally canalized phonon polaritons in LiV2O5 thin layers, Nat Commun 15, (2024)
work page 2024
-
[21]
G. Hu et al., Topological polaritons and photonic magic angles in twisted α-MoO3 bilayers, Nature 582, 209 (2020)
work page 2020
-
[22]
M. Chen, X. Lin, T. H. Dinh, Z. Zheng, J. Shen, Q. Ma, H. Chen, P. Jarillo -Herrero, and S. Dai, Configurable phonon polaritons in twisted α-MoO3, Nat Mater 19, 1307 (2020)
work page 2020
- [23]
-
[24]
J. Duan, N. Capote-Robayna, J. Taboada-Gutiérrez, G. Álvarez-Pérez, I. Prieto, J. Martín-Sánchez, A. Y . Nikitin, and P. Alonso-González, Twisted Nano-Optics: Manipulating Light at the Nanoscale with Twisted Phonon Polaritonic Slabs, Nano Lett 20, 5323 (2020)
work page 2020
-
[25]
J. Duan et al., Multiple and spectrally robust photonic magic angles in reconfigurable α-MoO3 trilayers, Nat Mater 22, 867 (2023)
work page 2023
-
[26]
N. C. Passler and A. Paarmann, Generalized 4 × 4 matrix formalism for light propagation in anisotropic stratified media: study of surface phonon polaritons in polar dielectric heterostructures, Journal of the Optical Society of America B 34, 2128 (2017)
work page 2017
-
[27]
G. Álvarez-Pérez, A. González-Morán, N. Capote-Robayna, K. V . V oronin, J. Duan, V . S. V olkov, P. Alonso-González, and A. Y . Nikitin, Active Tuning of Highly Anisotropic Phonon Polaritons in Van der Waals Crystal Slabs by Gated Graphene, ACS Photonics 9, 383 (2022)
work page 2022
-
[28]
H. Hu et al., Doping -driven topological polaritons in graphene/ α-MoO3 heterostructures, Nat Nanotechnol 17, 940 (2022)
work page 2022
-
[29]
Z. Zhou et al., Gate -Tuning Hybrid Polaritons in Twisted α-MoO3/Graphene Heterostructures, Nano Lett 23, 11252 (2023)
work page 2023
-
[30]
A. Y . Nikitin, F. Guinea, F. J. Garcia-Vidal, and L. Martin-Moreno, Fields radiated by a nanoemitter in a graphene sheet, Phys Rev B 84, 195446 (2011)
work page 2011
-
[31]
A. Y . Nikitin, F. J. Garcia -Vidal, and L. Martin -Moreno, Analytical expressions for the electromagnetic dyadic green’s function in graphene and thin layers, IEEE Journal on Selected Topics in Quantum Electronics 19, (2013)
work page 2013
-
[32]
Y. Hu, J. Álvarez-Cuervo, E. Terán-García, X. Huang, P. Alonso-González. Modulation of Radiative Heat Transfer at the Nanoscale via Topological Polaritons in Twisted van der Waals Crystals . Nanophotonics: e70000. 18
-
[33]
S. Molesky, Z. Lin, A. Y . Piggott, W. Jin, J. Vucković, and A. W. Rodriguez, Inverse Design in Nanophotonics, Nature Photonics 12, 659–670 (2018)
work page 2018
-
[34]
M. He, J. Nolen, J. Nordlander, A. Cleri, N. Mcllwaine, T. Folland, Y . Tang, B. Landman, J. Maria, and J. Caldwell, Deterministic inverse design of Tamm plasmon thermal emitters with multi - resonant control. Nat. Mater. 20, 1663–1669 (2021)
work page 2021
-
[35]
J. F. Masson, J. S. Biggins, and E. Ringe, Machine learning for nanoplasmonics. Nat. Nanotechnol. 18, 111–123 (2023)
work page 2023
- [36]
-
[37]
N. C. Passler et al., Hyperbolic shear polaritons in low-symmetry crystals, Nature 602, 595 (2022)
work page 2022
-
[38]
G. Hu et al., Real -space nanoimaging of hyperbolic shear polaritons in a monoclinic crystal, Nat Nanotechnol 18, 64 (2023)
work page 2023
-
[39]
J. Matson et al., Controlling the Propagation Asymmetry of Hyperbolic Shear Polaritons in Beta - Gallium Oxide. Nat Commun 14, 5240 (2023)
work page 2023
-
[40]
E. Galiffi et al., Extreme Light Confinement and Control in Low -Symmetry Phonon-Polaritonic Crystals, Nat Rev Mater 9, 9–28 (2024)
work page 2024
-
[41]
G. A. Ermolaev et al., Wandering principal optical axes in van der Waals triclinic materials, Nat Commun 15, (2024)
work page 2024
-
[42]
G. Venturi, A. Mancini, N. Melchioni et al. Visible-frequency hyperbolic plasmon polaritons in a natural van der Waals crystal. Nat Commun 15, 9727 (2024)
work page 2024
-
[43]
Francesco L. Ruta et al. ,Good plasmons in a bad metal.Science 387,786-791(2025)
work page 2025
-
[44]
H. Rubens and E. F. Nichols. Heat Rays Of Great Wave Length. Phys. Rev. (Series I) 4, 314 (1897)
-
[45]
L. D. LANDAU and E. M. LIFSHITZ, CHAPTER I - ELECTROSTATICS OF CONDUCTORS, in Electrodynamics of Continuous Media (Second Edition) , edited by L. D. LANDAU and E. M. LIFSHITZ, V ol. 8 (Pergamon, Amsterdam, 1984), pp. 1–33
work page 1984
-
[46]
G. R. Fowles, Introduction to Modern Optics (Holt McDougal, Evanston, IL, 1975)
work page 1975
-
[47]
G. Álvarez-Pérez, K. V V oronin, V . S. V olkov, P. Alonso-González, and A. Y . Nikitin, Analytical approximations for the dispersion of electromagnetic modes in slabs of biaxial crystals, Phys Rev B 100, 235408 (2019)
work page 2019
-
[48]
See Supplementary Information
-
[49]
A. Y . U. NIKITIN, Graphene Plasmonics , in World Scientific Handbook of Metamaterials and Plasmonics (World Scientific, 2017), pp. 307–338
work page 2017
-
[50]
K. V. V oronin, G. Álvarez-Pérez, C. Lanza, P. Alonso-González, and A. Y . Nikitin, Fundamentals of Polaritons in Strongly Anisotropic Thin Crystal Layers, ACS Photonics 11, 550 (2024)
work page 2024
-
[51]
K. V . V oronin, G. Álvarez -Pérez, A. Tarazaga Martín -Luengo, P. Alonso-González and A. Y . Nikitin. Misalignment between the Directions of Propagation and Decay of Nanoscale -Confined Polaritons Nano Lett. 2025, 25, 37, 13811–13818
work page 2025
-
[52]
Л. А. Фальковский, Оптические свойства графена и полупроводников типа A4B6, Uspekhi Fizicheskih Nauk 178, 923 (2008). 19
work page 2008
-
[53]
J. Chen, M. Badioli, P. Alonso-González, S. Thongrattanasiri, F. Huth, J. Osmond, M. Spasenovic, A. Centeno, A. Pesquera, P. Godignon, A. Zurutuza Elorza, N. Camara, F. J. García de Abajo, R. Hillenbrand & F. H. L. Koppens, Optical nano-imaging of gate-tunable graphene plasmons. Nature 487, 77–81 (2012)
work page 2012
- [54]
-
[55]
C. Lanza, Twistoptics, GitLab repository, available at https://gitlab.com/ChristianLanza/twistoptics (2026)
work page 2026
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.