pith. sign in

arxiv: 2604.05097 · v1 · submitted 2026-04-06 · ❄️ cond-mat.mes-hall

Surface Response, Plasma Modes of coated Multi-Layered anisotropic Semi-Dirac Heterostructures

Teresa Lee , Godfrey Gumbs , Thi Nga Do , Andrii Iurov , Danhong Huang This is my paper

Pith reviewed 2026-05-10 18:52 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall
keywords surface response functionsplasmon dispersionssemi-Dirac heterostructuresanisotropic plasma modesmulti-layer coatingsloss functionsoptical absorption spectra
0
0 comments X

The pith

Closed-form analytical expressions for surface response functions describe plasma modes in coated multi-layered semi-Dirac heterostructures.

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

The paper derives closed-form analytical expressions for the surface response functions of heterostructures made from up to three layers of two-dimensional materials coated with dielectric media. These expressions allow calculation of plasmon dispersions in the long-wavelength limit for single and double layers of anisotropic semi-Dirac materials. A sympathetic reader would care because such analytical results simplify the study of Coulomb-coupled excitations in these structures, which could inform designs for protective coatings with specific optical properties. The work uses Maxwell's equations and linear response theory to model how an external electromagnetic field induces charge density oscillations. It shows anisotropic plasmon behavior and distinct in-phase and out-of-phase modes for multi-layer cases.

Core claim

We derived closed-form analytical expressions for the surface response functions for heterostructures consisting of up to three layered coated two-dimensional materials with dielectric interfaces. For anisotropic semi-Dirac materials, closed-form expressions for plasmon dispersions are obtained in the long wavelength limit, and numerical density plots of loss functions reveal anisotropic behavior. Multi-layer structures exhibit two plasmon branches corresponding to in-phase and out-of-phase oscillations, with optical modes having higher intensity.

What carries the argument

The surface response function derived from Maxwell's equations and linear response theory, which yields the Coulomb coupled plasma excitations from an incident electromagnetic field.

If this is right

  • Analytical expressions simplify computation of plasmon dispersions without numerical approximation in the long wavelength limit.
  • Multi-layer configurations produce distinct optical and acoustic plasmon modes with different intensities.
  • Anisotropic behavior in different momentum directions affects the loss function plots for semi-Dirac materials.
  • Optical absorption spectra can be calculated for plasma modes under external fields with specific polarization and frequency.

Where Pith is reading between the lines

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

  • The derived expressions may apply to additional layer counts or material combinations if charge transfer remains inhibited.
  • Experiments with varying light incidence angles on tilted semi-Dirac materials could test the predicted anisotropy in plasmon responses.
  • Tuning the number of layers might optimize intensity ratios for applications in protective optical coatings.

Load-bearing premise

The dielectric media serves to inhibit charge transfer between layers allowing linear response theory and Maxwell's equations to capture all Coulomb-coupled excitations.

What would settle it

Detection of significant charge transfer across the dielectric film between coated layers would contradict the assumptions underlying the closed-form SRF expressions.

Figures

Figures reproduced from arXiv: 2604.05097 by Andrii Iurov, Danhong Huang, Godfrey Gumbs, Teresa Lee, Thi Nga Do.

Figure 1
Figure 1. Figure 1: shows the energy dispersion derived from the semi-Dirac Hamiltonian when either kx or ky is set equal to zero. The energy is in unit of EF for graphene ℏυF kF . Fig. 1a is the energy dispersion plotted as a function of kx/kF when ky = 0. We define an inverse effective mass parameter a ′ as a/aF , the inverse effective mass in unit of aF , which came from the relation ℏk 2 F aF = ℏυF kF , thus aF = υF /kF o… view at source ↗
Figure 2
Figure 2. Figure 2: shows the energy dispersion when the SOC is present. The magnitude of ∆′ which is defined as ∆/EF is set equal to 0.2. As shown in Fig. 2a, when the energy dispersion is plotted with respect to kx/kF , we see the same magnitude of energy gaps appearing for different magnitudes of inverse effective mass parameter a ′ . Here, since tilting only takes effect in the ky direction, it does not affect the energy … view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Three-dimensional energy dispersion in unit of [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Schematic illustration of a layered structure with two of the electron monolayers serving as outer covers on [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Schematic illustration of a pair of electron monolayers on the surface of a material with dielectric constant [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Schematic representation of a single layer placed on a thick substrate with dielectric constant [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: (Color online) (a) Plasmon dispersion in the long wavelength limit for a monolayer of semi- Dirac material [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: (Color online) With [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: (Color online) Real and imaginary parts of [PITH_FULL_IMAGE:figures/full_fig_p013_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10: (Color online) (a) Plasmon dispersion in the long wavelength limit for two monolayers surrounding a [PITH_FULL_IMAGE:figures/full_fig_p014_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11: (Color online) Density plots for imaginary parts of 1 [PITH_FULL_IMAGE:figures/full_fig_p014_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12: (Color online) Density plots for imaginary parts of 1 [PITH_FULL_IMAGE:figures/full_fig_p015_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13: (Color online) Absorption function in unit of [PITH_FULL_IMAGE:figures/full_fig_p016_13.png] view at source ↗
read the original abstract

We derived closed-form analytical expressions for the surface response functions (SRFs) for heterostructure. We investigate structures consisting of up to three layered, coated heterostructure of two-dimensional (2D) materials with a dielectric medium or vacuum interface. The dielectric media serves to inhibit charge transfer between layers for the case when a pair of 2D layers serve as coatings for a dielectric film. Our results revise the established picture for the dispersion equation for two layers of reduced dimensionality surrounded by dielectric media. An impinging electromagnetic field incident on the surface leads to Coulomb coupled plasma excitations in the structure which are yielded by the SRF. This is achieved by employing Maxwell's equations and linear response theory. We use these results to investigate the plasmonic properties of tilted semi-Dirac materials both analytically and numerically. Closed-form analytical expressions are derived for the plasmon dispersions in the long wavelength limit for single and double layers. We numerically obtain density plots of the loss functions and observe anisotropic behavior in different momentum directions. For the cases when there are two or three layers, we observe two plasmon branches corresponding to in-phase and out-of-phase charge density oscillations, where the in-phase optical modes have higher intensity than the out-of-phase acoustic modes. We calculated the optical absorption spectra for plasma modes in layered semi-Dirac materials produced by an external electromagnetic field carrying an electric polarization and frequency. Possible applications include durable protection coatings providing UV resistance, chemical protection and improving upon traditional ceramic coatings.

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

0 major / 2 minor

Summary. The manuscript claims to derive closed-form analytical expressions for the surface response functions (SRFs) of up to three-layered coated heterostructures of anisotropic semi-Dirac 2D materials using Maxwell's equations and linear response theory. It provides closed-form plasmon dispersion relations in the long-wavelength limit for single and double layers, numerical loss function density plots demonstrating anisotropic behavior, and optical absorption spectra for the plasma modes, with observations of in-phase optical and out-of-phase acoustic modes in multi-layer systems.

Significance. Should the closed-form expressions prove correct upon verification, this work offers significant analytical tools for modeling plasmonic excitations in anisotropic 2D heterostructures, which are often limited to numerical methods. The explicit revision of the two-layer dispersion equation and the distinction between mode intensities provide new insights. The numerical results on loss functions and absorption spectra are concrete contributions that could guide experimental studies in plasmonics.

minor comments (2)
  1. [Numerical Results] The loss function density plots would benefit from clearer labeling of the momentum directions (e.g., along the linear vs. quadratic dispersion axes of the semi-Dirac cone) and the specific parameter values used for the anisotropy and carrier densities.
  2. [Conclusion] The suggested applications in durable protection coatings, UV resistance, and chemical protection are speculative and should either be supported by references to relevant experimental literature or de-emphasized in favor of the core theoretical plasmonic results.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our work on closed-form surface response functions and plasmon dispersions in anisotropic semi-Dirac heterostructures. The recommendation for minor revision is noted. No specific major comments were provided in the report, so we address the overall evaluation below.

Circularity Check

0 steps flagged

Derivation self-contained from Maxwell equations and linear response

full rationale

The paper states that closed-form SRFs and long-wavelength plasmon dispersions are obtained by applying Maxwell's equations together with linear response theory to coated anisotropic semi-Dirac heterostructures (up to three layers). The dielectric-medium assumption that inhibits interlayer charge transfer is a conventional modeling choice for electrostatic treatments and does not reduce any derived expression to a fitted input or prior self-citation. No load-bearing self-citations, self-definitional steps, or renaming of known results are exhibited in the abstract or described derivation chain. The central analytical results therefore remain independent of the paper's own inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Abstract provides insufficient detail to enumerate specific free parameters or invented entities; the central derivations rest on standard assumptions of linear response and dielectric isolation between layers.

axioms (2)
  • domain assumption Linear response theory and Maxwell's equations fully capture the Coulomb-coupled plasma excitations.
    Invoked to obtain the surface response functions and plasmon dispersions from an impinging electromagnetic field.
  • domain assumption Dielectric media completely inhibit charge transfer between 2D layers.
    Stated as the role of the coating to enable independent layer response.

pith-pipeline@v0.9.0 · 5582 in / 1329 out tokens · 61784 ms · 2026-05-10T18:52:10.348311+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

34 extracted references · 34 canonical work pages

  1. [1]

    5: Schematic illustration of a pair of electron monolayers on the surface of a material with dielectric constant ϵ1 and suspended in vacuum

    +s 2c1ϵ1(1−ϵ 2 2)−s 2s1(ϵ2 1 −ϵ 2 2),(11) D(0) 3 =−c 2ϵ2 h 2c1ϵ1 +s 1(1 +ϵ 2 1) i −s 2 h c1ϵ1(ϵ2 2 + 1) +s 1(ϵ2 1 +ϵ 2 2) i =−2c 2ϵ2c1ϵ1 −c 2ϵ2s1(1 +ϵ 2 1)−s 2c1ϵ1(1 +ϵ 2 2)−s 2s1(ϵ2 1 +ϵ 2 2),(12) P3 =c 2ϵ2(c1ϵ1 +s 1) +s 2(c1ϵ1 +s 1ϵ2 2),(13) A3 =c 1ϵ1 +s 1(1−α 0),(14) B3 =c 1ϵ1 −s 1(1 +α 0),(15) γ3 =α 1(c2ϵ2 +s 2) +α 2ϵ2.(16) 7 vacuum FIG. 5: Schematic ...

    work page

  2. [2]

    Plasmons for Single layer on a thick substrate We employ the density matrix ˆρ(r, t) to determine the charge density fluctuationσ(r, ω) . Starting with the non- dissipative equation of motioniℏ∂ˆρ(r, t)/∂t= [ ˆH(t),ˆρ(r, t)] , we separate both the Hamiltonian ˆH(t) = ˆH0 + ˆH1(t) and the density matrix ˆρ(t) = ˆρ0 + ˆρ1(t) into an unperturbed time-indepen...

    work page arXiv

  3. [3]

    Taking the long wavelength limit, then settingD 2s(q, ω) = 0 in Eq

    Plasma excitation for two layers covering a dielectric Similarly, we determine the plasmon modes for two electronic monolayers on the surfaces of a dielectric medium, and all together resting on a thick substrate. Taking the long wavelength limit, then settingD 2s(q, ω) = 0 in Eq. (23), we obtain two plasmon branches in two mutually perpendicular directio...

    work page

  4. [4]

    12: (Color online) Density plots for imaginary parts of 1/D 3(q, ω) for triple monolayers suspended in vacuum

    Plasma excitations for Three monolayers FIG. 12: (Color online) Density plots for imaginary parts of 1/D 3(q, ω) for triple monolayers suspended in vacuum. All three layers have the same tiltingτ= 0.1 and spin-orbit coupling parameter ∆ ′ = 0.1, and to take account of Landau damping, we choseℏδ +/EF = 0.01 in the polarization function. The relative locati...

    work page

  5. [5]

    A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, The electronic properties of graphene, Rev. Mod. Phys81, 109(2009)

    work page 2009

  6. [6]

    Wu, S.-C

    J.-Y. Wu, S.-C. Chen, G. Gumbs, and M.-F. Lin, Feature-rich electronic excitations of silicene in external fields,Phys. Rev. B94, 205427(2016)

    work page 2016

  7. [7]

    A. Kara, H. Enriquez, A. P. Seitsonen, L. C. Lew Yan Voon, S. Vizzini, B. Aufray, and H. Oughaddou, A review on silicene – new candidate for electronics,Surf. Sci. Rep.67, 1(2012)

    work page 2012

  8. [8]

    Zhao, et al.,Rise of silicene: A competitive 2D material,Prog

    J. Zhao, et al.,Rise of silicene: A competitive 2D material,Prog. Mater. Sci.83, 24(2016)

    work page 2016

  9. [9]

    Carvalho, M

    A. Carvalho, M. Wang, X. Zhu, A. S. Rodin, H. Su, and A. H. Castro Neto, Phosphorene from theory to applications, Nat. Rev. Mater.1, 16061(2016)

    work page 2016

  10. [10]

    A. Acun, L. Zhang, P. Bampoulis, M. Farmanbar, A. van Houselt, A. N. Rudenko, M. Lingenfelder, G. Brocks, B. Poelsema, M. I. Katsnelson, and H. J. W. Zandvliet, Germanene: the germanium analogue of graphene,J. Phys.: Condens. Matter 27, 443002(2015)

    work page 2015

  11. [11]

    Sadhukhan and A

    K. Sadhukhan and A. Agarwal, Anisotropic plasmons, Friedel oscillations, and screening in 8−P mmnborophene, Phys. Rev. B96, 035410(2017). 19

    work page 2017

  12. [12]

    Balassis, G

    A. Balassis, G. Gumbs, and O. Roslyak, Polarizability, plasmons, and screening in 1T’-MoS2 with tilted Dirac bands,Phys. Lett. A449, 128353(2022)

    work page 2022

  13. [13]

    Lopez-Bezanilla and P

    A. Lopez-Bezanilla and P. B. Littlewood, Electronic properties of 8−P mmnborophene,Phys. Rev. B93, 241405(R)(2016)

    work page 2016

  14. [14]

    Ross-Harvey, A

    G. Ross-Harvey, A. Iurov, L. Zhemchuzhna, G. Gumbs, D. Huang, and P. Fekete, Dynamical polarization function, anisotropic plasmon modes, and dephasing rates for interacting electrons in semi-Dirac bands,Phys. Rev. B111, 045413(2025)

    work page 2025

  15. [15]

    Yan, C.-Y

    C.-X. Yan, C.-Y. Tan, H. Guo, and H.-R. Chang, Highly anisotropic optical conductivities in two-dimensional tilted semi-Dirac bands,Phys. Rev. B108, 195427(2023)

    work page 2023

  16. [16]

    Kim, S.S

    J. Kim, S.S. Baik, S.H. Ryu, Y. Sohn, S. Park, B.-G. Park, J. Denlinger, Y. Yi, H. J. Choi, and K. S. Kim, Observation of tunable band gap and anisotropic Dirac semimetal state in black phosphorus,Science349, 723–726(2015)

    work page 2015

  17. [17]

    Real, et al., Semi-Dirac transport and anisotropic localization in polariton honeycomb lattices,Phys

    B. Real, et al., Semi-Dirac transport and anisotropic localization in polariton honeycomb lattices,Phys. Rev. Lett.125, 186601(2020)

    work page 2020

  18. [18]

    Dietl, F

    P. Dietl, F. Pi´ echon, and G. Montambaux, New magnetic field dependence of Landau levels in a graphenelike structure, Phys. Lett.100, 236405(2008)

    work page 2008

  19. [19]

    Banerjee , R

    S. Banerjee , R. R. P. Singh, V. Pardo, and W. E. Pickett, Tight-binding modeling and low-energy behavior of the semi-Dirac point,Phys. Rev. Lett.103, 016402(2009)

    work page 2009

  20. [20]

    Pardo, and W

    V. Pardo, and W. E. Pickett, Metal-insulator transition through a semi-Dirac point in oxide nanostructures: VO 2 (001) layers confined within TiO 2.Phys. Rev. B81, 035111(2010)

    work page 2010

  21. [21]

    Delplace, G

    P. Delplace, G. Montambaux, Semi-Dirac point in the Hofstadter spectrum,Phys. Rev. B82, 035438(2010)

    work page 2010

  22. [22]

    Saha, Photoinduced Chern insulating states in semi-Dirac materials,Phys

    K. Saha, Photoinduced Chern insulating states in semi-Dirac materials,Phys. Rev. B94, 081103(R)(2016)

    work page 2016

  23. [23]

    V. N. Kotov, B. Uchoa, and O. P. Sushkov, Coulomb interactions and renormalization of semi-Dirac fermions near a topological Lifshitz transition, Phys. Rev. B103, 045403(2021)

    work page 2021

  24. [24]

    M. M. Elsayed, B. Uchoa, and V. N. Kotov, Coulomb interactions in systems of generalized semi-Dirac fermions,Phys. Rev. B111, 165127(2025)

    work page 2025

  25. [25]

    D. Giri, J. Schliemann, R. A. Molina, and A. Lopez, Tunable plasmon modes and topological transitions in monolayer and bilayer semi-Dirac materials,Phys. Rev. B112165420(2025)

    work page 2025

  26. [26]

    J. K. Jain and P. B. Allen, Dielectric response of a semi-infinite layered electron gas and Raman scattering from its bulk and surface plasmons,Phys. Rev. B32, 997(1985)

    work page 1985

  27. [27]

    Gumbs, A

    G. Gumbs, A. Iurov, J.-Y. Wu, M. F. Lin, and P. Fekete, Plasmon excitations of multi-layer graphene on a conducting substrate,Sci. Rep.6, 21063(2016)

    work page 2016

  28. [28]

    Dahal, G

    D. Dahal, G. Gumbs, and D. Huang, Effect of strain on plasmons, screening, and energy loss in graphene/substrate contacts,Phys. Rev. B98, 045427(2018)

    work page 2018

  29. [29]

    G. F. Giuliani and J. J. Quinn, Charge-density excitations at the surface of a semiconductor superlattice: A new type of surface polariton,Phys. Rev. Lett.51, 919(1983)

    work page 1983

  30. [30]

    Das Sarma, and A

    S. Das Sarma, and A. Madhukar. Collective modes of spatially separated, two-component, two-dimensional plasma in solids,Phys. Rev. B23, 805(1981)

    work page 1981

  31. [31]

    Iurov, L

    A. Iurov, L. Zhemchuzhna, G. Gumbs, and D. Huang: Exploring stable long-lifetime plasmon excitations in a Lieb lattice, Annalen der Physik (2026)

    work page 2026

  32. [32]

    Gumbs, and D

    G. Gumbs, and D. Huang,Properties of Interacting Low Dimensional Systems(Wiley-VCH, 2011)

    work page 2011

  33. [33]

    Solid State Commun.52, 811(1984)

    B.N.J Persson, Inelastic electron scattering from thin metal films. Solid State Commun.52, 811(1984)

    work page 1984

  34. [34]

    N. J. M. Horing, H. C. Tso, and G. Gumbs, Fast-particle energy loss in the vicinity of a two-dimensional plasma,Phys. Rev. B36, 1588(1987)

    work page 1987