pith. sign in

arxiv: 2604.13907 · v1 · submitted 2026-04-15 · ⚛️ physics.optics · cond-mat.mtrl-sci

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

classification ⚛️ physics.optics cond-mat.mtrl-sci
keywords twistopticspolaritonsanisotropic heterostructuresanalytical modelvan der Waals materialselectromagnetic wavestwisted layersconductive sheets
0
0 comments X

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.

The paper introduces an analytical model capable of describing light propagation and polariton modes in planar stacks consisting of any number of rotated thin anisotropic layers and conductive sheets. This generalizes previous limited cases to arbitrary configurations while providing formulas for key quantities like the polariton wavelength and how far it travels before decaying. Sympathetic readers care because it removes the need for custom numerical modeling for each new stack design, potentially speeding up the exploration of twist-controlled phenomena such as canalization for imaging or heat flow applications. The model comes with simplified versions for high-momentum waves and very thin layers, plus computer scripts for direct use.

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

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

  • 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

Figures reproduced from arXiv: 2604.13907 by Aitana Tarazaga Mart\'in-Luengo, Alexey Y. Nikitin, Christian Lanza, Gonzalo \'Alvarez-P\'erez, Javier Mart\'in-S\'anchez, Jos\'e \'Alvarez-Cuervo, Kirill V. Voronin, Pablo Alonso-Gonz\'alez.

Figure 2
Figure 2. Figure 2: b) [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗
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.

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

1 major / 2 minor

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)
  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)
  1. [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.
  2. [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

1 responses · 0 unresolved

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
  1. 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

0 steps flagged

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

0 free parameters · 2 axioms · 0 invented entities

The model rests on standard electromagnetic theory for anisotropic media and interface continuity conditions. No new particles or forces are introduced. Material permittivity tensors are treated as known inputs from experiment or prior literature rather than fitted within the paper.

axioms (2)
  • domain assumption Local anisotropic permittivity tensor description of each layer
    Invoked when defining the constitutive relations for biaxial layers and conductive sheets.
  • standard math Validity of standard electromagnetic boundary conditions at each interface
    Required to match fields across rotated layers.

pith-pipeline@v0.9.0 · 5584 in / 1235 out tokens · 22034 ms · 2026-05-10T12:37:31.725738+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

55 extracted references · 55 canonical work pages

  1. [1]

    A. K. Geim and I. V Grigorieva, Van der Waals heterostructures, Nature 499, 419 (2013)

  2. [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)

  3. [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)

  4. [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)

  5. [5]

    Cao et al., Correlated insulator behaviour at half-filling in magic-angle graphene superlattices, Nature 556, 80 (2018)

    Y . Cao et al., Correlated insulator behaviour at half-filling in magic-angle graphene superlattices, Nature 556, 80 (2018)

  6. [6]

    Yankowitz, S

    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)

  7. [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)

  8. [8]

    G. Hu, A. Krasnok, Y. Mazor, C.-W. Qiu, and A. Alù. Moirè Hyperbolic Metasurfaces. Nano Lett. 2020, 20, 5, 3217–3224

  9. [9]

    Álvarez-Cuervo, M

    J. Álvarez-Cuervo, M. Obst, S. Dixit et al. Unidirectional ray polaritons in twisted asymmetric stacks. Nat Commun 15, 9042 (2024)

  10. [10]

    Hillenbrand, Y

    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)

  11. [11]

    Dai et al., Tunable Phonon Polaritons in Atomically Thin van der Waals Crystals of Boron Nitride, Science (1979) 343, 1125 (2014)

    S. Dai et al., Tunable Phonon Polaritons in Atomically Thin van der Waals Crystals of Boron Nitride, Science (1979) 343, 1125 (2014)

  12. [12]

    Dai et al., Subdiffractional focusing and guiding of polaritonic rays in a natural hyperbolic material, Nat Commun 6, 6963 (2015)

    S. Dai et al., Subdiffractional focusing and guiding of polaritonic rays in a natural hyperbolic material, Nat Commun 6, 6963 (2015)

  13. [13]

    Dai et al., Efficiency of Launching Highly Confined Polaritons by Infrared Light Incident on a Hyperbolic Material, Nano Lett 17, 5285 (2017)

    S. Dai et al., Efficiency of Launching Highly Confined Polaritons by Infrared Light Incident on a Hyperbolic Material, Nano Lett 17, 5285 (2017)

  14. [14]

    A. J. Giles et al., Ultralow -loss polaritons in isotopically pure boron nitride, Nat Mater 17, 134 (2018). 17

  15. [15]

    Ma et al., In-plane anisotropic and ultra-low-loss polaritons in a natural van der Waals crystal, Nature 562, 557 (2018)

    W. Ma et al., In-plane anisotropic and ultra-low-loss polaritons in a natural van der Waals crystal, Nature 562, 557 (2018)

  16. [16]

    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)

    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)

  17. [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)

  18. [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)

  19. [19]

    Li et al., Collective near -field coupling and nonlocal phenomena in infrared -phononic metasurfaces for nano-light canalization, Nat Commun 11, 3663 (2020)

    P. Li et al., Collective near -field coupling and nonlocal phenomena in infrared -phononic metasurfaces for nano-light canalization, Nat Commun 11, 3663 (2020)

  20. [20]

    A. I. F. Tresguerres-Mata et al., Observation of naturally canalized phonon polaritons in LiV2O5 thin layers, Nat Commun 15, (2024)

  21. [21]

    Hu et al., Topological polaritons and photonic magic angles in twisted α-MoO3 bilayers, Nature 582, 209 (2020)

    G. Hu et al., Topological polaritons and photonic magic angles in twisted α-MoO3 bilayers, Nature 582, 209 (2020)

  22. [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)

  23. [23]

    Zheng, F

    Z. Zheng, F. Sun, W. Huang, J. Jiang, R. Zhan, Y . Ke, H. Chen, and S. Deng, Phonon Polaritons in Twisted Double-Layers of Hyperbolic van der Waals Crystals, Nano Lett. 2020 20, 7, 5301–5308

  24. [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)

  25. [25]

    Duan et al., Multiple and spectrally robust photonic magic angles in reconfigurable α-MoO3 trilayers, Nat Mater 22, 867 (2023)

    J. Duan et al., Multiple and spectrally robust photonic magic angles in reconfigurable α-MoO3 trilayers, Nat Mater 22, 867 (2023)

  26. [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)

  27. [27]

    Álvarez-Pérez, A

    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)

  28. [28]

    Hu et al., Doping -driven topological polaritons in graphene/ α-MoO3 heterostructures, Nat Nanotechnol 17, 940 (2022)

    H. Hu et al., Doping -driven topological polaritons in graphene/ α-MoO3 heterostructures, Nat Nanotechnol 17, 940 (2022)

  29. [29]

    Zhou et al., Gate -Tuning Hybrid Polaritons in Twisted α-MoO3/Graphene Heterostructures, Nano Lett 23, 11252 (2023)

    Z. Zhou et al., Gate -Tuning Hybrid Polaritons in Twisted α-MoO3/Graphene Heterostructures, Nano Lett 23, 11252 (2023)

  30. [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)

  31. [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)

  32. [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. [33]

    Molesky, Z

    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)

  34. [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)

  35. [35]

    J. F. Masson, J. S. Biggins, and E. Ringe, Machine learning for nanoplasmonics. Nat. Nanotechnol. 18, 111–123 (2023)

  36. [36]

    Claus, R

    R. Claus, R. CLAUS: Polariton Dispersion and Crystal Optics in Monoclinic Materials Polariton Dispersion and Crystal Optics in Monoclinic Materials, 1978

  37. [37]

    N. C. Passler et al., Hyperbolic shear polaritons in low-symmetry crystals, Nature 602, 595 (2022)

  38. [38]

    Hu et al., Real -space nanoimaging of hyperbolic shear polaritons in a monoclinic crystal, Nat Nanotechnol 18, 64 (2023)

    G. Hu et al., Real -space nanoimaging of hyperbolic shear polaritons in a monoclinic crystal, Nat Nanotechnol 18, 64 (2023)

  39. [39]

    Matson et al., Controlling the Propagation Asymmetry of Hyperbolic Shear Polaritons in Beta - Gallium Oxide

    J. Matson et al., Controlling the Propagation Asymmetry of Hyperbolic Shear Polaritons in Beta - Gallium Oxide. Nat Commun 14, 5240 (2023)

  40. [40]

    Galiffi et al., Extreme Light Confinement and Control in Low -Symmetry Phonon-Polaritonic Crystals, Nat Rev Mater 9, 9–28 (2024)

    E. Galiffi et al., Extreme Light Confinement and Control in Low -Symmetry Phonon-Polaritonic Crystals, Nat Rev Mater 9, 9–28 (2024)

  41. [41]

    G. A. Ermolaev et al., Wandering principal optical axes in van der Waals triclinic materials, Nat Commun 15, (2024)

  42. [42]

    Venturi, A

    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)

  43. [43]

    Ruta et al

    Francesco L. Ruta et al. ,Good plasmons in a bad metal.Science 387,786-791(2025)

  44. [44]

    Rubens and E

    H. Rubens and E. F. Nichols. Heat Rays Of Great Wave Length. Phys. Rev. (Series I) 4, 314 (1897)

  45. [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

  46. [46]

    G. R. Fowles, Introduction to Modern Optics (Holt McDougal, Evanston, IL, 1975)

  47. [47]

    Álvarez-Pérez, K

    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)

  48. [48]

    See Supplementary Information

  49. [49]

    A. Y . U. NIKITIN, Graphene Plasmonics , in World Scientific Handbook of Metamaterials and Plasmonics (World Scientific, 2017), pp. 307–338

  50. [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)

  51. [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

  52. [52]

    Л. А. Фальковский, Оптические свойства графена и полупроводников типа A4B6, Uspekhi Fizicheskih Nauk 178, 923 (2008). 19

  53. [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)

  54. [54]

    Zheng, G

    C. Zheng, G. Hu, X. Liu, X. Kong, L. Wang, and C. W. Qiu, Molding Broadband Dispersion in Twisted Trilayer Hyperbolic Polaritonic Surfaces, ACS Nano 16, 13241 (2022)

  55. [55]

    Lanza, Twistoptics, GitLab repository, available at https://gitlab.com/ChristianLanza/twistoptics (2026)

    C. Lanza, Twistoptics, GitLab repository, available at https://gitlab.com/ChristianLanza/twistoptics (2026)