pith. sign in

arxiv: 2606.12274 · v1 · pith:NPALI7PQnew · submitted 2026-06-10 · ❄️ cond-mat.mtrl-sci

Plasmonic properties and correlation energies from a compact multipole representation of the dielectric response in 2D metals

Dario A. Leon , Claudia Cardoso , Kristian Berland This is my paper

Pith reviewed 2026-06-27 09:04 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords 2D metalsdielectric responseplasmonic modescorrelation energiesmultipole-Padé approximantsdynamical screeningBrillouin zoneab initio calculations
0
0 comments X

The pith

A compact multipole representation of the inverse dielectric function captures the full momentum-dependent response of 2D metals using few dispersive plasmon modes.

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

The work generalizes multipole-Padé approximants to include momentum dependence over the entire Brillouin zone of two-dimensional metals through a symmetry-conserving anisotropic form of the inverse dielectric function. This compact analytic expression is built from ab initio data for seven metals that span different electronic regimes. The resulting representation reproduces the dielectric response accurately with only a small number of collective modes. It also delivers reliable values for correlation energies that involve dynamical screening. The construction therefore connects detailed numerical calculations directly to simpler analytical treatments of screening effects.

Core claim

Multipole-Padé approximants are generalized to two-dimensional metals by constructing a symmetry-conserving, anisotropic representation of the inverse dielectric function that includes momentum dependence over the full Brillouin zone. This analytic form enables efficient and accurate evaluation of dynamical screening quantities. For seven two-dimensional metals, a small number of dispersive plasmonic modes suffices to describe the dielectric response accurately across the Brillouin zone and to yield accurate correlation energies.

What carries the argument

symmetry-conserving anisotropic multipole-Padé representation of the inverse dielectric function, which encodes the full momentum-dependent response using a small number of dispersive collective modes

If this is right

  • Dynamical screening quantities can be evaluated efficiently and accurately over the full Brillouin zone.
  • Spectral features arising from plasmon modes are reproduced to high accuracy.
  • Correlation energies computed from the representation match those obtained from the complete ab initio response.
  • A direct connection is established between ab initio data and analytical models of screening in 2D systems.

Where Pith is reading between the lines

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

  • The same compact form could be applied to response functions beyond the dielectric function in other low-dimensional materials.
  • The reduced number of modes may allow analytic progress on many-body problems that currently require full numerical integration over momentum.
  • Device-scale modeling of screening in 2D heterostructures could become feasible if the parameters are precomputed once per material.

Load-bearing premise

The inverse dielectric function of 2D metals admits an accurate symmetry-conserving anisotropic multipole-Padé representation whose parameters can be fixed from ab initio calculations without losing essential momentum dependence.

What would settle it

Direct numerical comparison, for any one of the seven metals, between the dielectric function or the correlation energy obtained from the multipole-Padé approximant and the same quantities obtained from the full ab initio dielectric function, revealing large systematic deviations.

Figures

Figures reproduced from arXiv: 2606.12274 by Claudia Cardoso, Dario A. Leon, Kristian Berland.

Figure 1
Figure 1. Figure 1: FIG. 1. Spectral properties of the dielectric response of monolayer Na obtained with MPA( [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Comparison of the [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Convergence of the RPA correlation energy [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
read the original abstract

Multipole-Pad\'e approximants provide a compact representation of dynamical response functions in terms of a small number of collective modes. Here, we generalize this framework to incorporate momentum dependence across the full Brillouin zone of 2D metals by constructing a symmetry-conserving, anisotropic representation of the inverse dielectric function. This analytic form enables efficient and accurate evaluation of quantities involving dynamical screening, including spectral features and correlation energies. We construct such compact representations for a set of seven two dimensional metals spanning distinct electronic regimes, and show that a small number of dispersive plasmonic modes suffices to accurately describe the dielectric response across the full Brillouin zone, while also yielding accurate correlation energies. The proposed representation therefore establishes a direct bridge between {\it ab initio} calculations and analytical models of screening, opening new avenues for applications in condensed matter systems.

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 generalizes multipole-Padé approximants to a symmetry-conserving anisotropic representation of the inverse dielectric function in 2D metals that incorporates momentum dependence over the full Brillouin zone. It constructs such compact forms for seven 2D metals spanning different electronic regimes and reports that a small number of dispersive plasmonic modes suffices to reproduce the dielectric response while also producing accurate correlation energies, thereby bridging ab initio calculations and analytical screening models.

Significance. If the numerical validations hold, the work supplies a practical analytic bridge between first-principles dielectric data and model treatments of dynamical screening. The explicit multi-material construction and the emphasis on compactness and symmetry conservation are strengths that could enable efficient calculations of plasmonic spectra and correlation energies in 2D systems.

major comments (1)
  1. [Results section (validation for the seven materials)] The central claim of accuracy rests on fitting the multipole parameters to external ab initio dielectric data; the manuscript should therefore report quantitative error metrics (e.g., mean absolute deviation or maximum relative error in the inverse dielectric function) across the full Brillouin zone for each of the seven materials, together with direct comparisons of the resulting correlation energies to independent benchmarks.
minor comments (2)
  1. [Abstract] The abstract states that the representation is constructed from ab initio data but does not indicate the typical number of poles/residues retained or the momentum sampling used; adding this information would clarify the compactness claim.
  2. [Theory section] Notation for the anisotropic multipole-Padé form (poles, residues, and symmetry constraints) should be defined explicitly in the main text before the applications to specific materials.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive assessment and the recommendation for minor revision. We address the single major comment below by agreeing to strengthen the quantitative validation as suggested.

read point-by-point responses

  1. Referee: [Results section (validation for the seven materials)] The central claim of accuracy rests on fitting the multipole parameters to external ab initio dielectric data; the manuscript should therefore report quantitative error metrics (e.g., mean absolute deviation or maximum relative error in the inverse dielectric function) across the full Brillouin zone for each of the seven materials, together with direct comparisons of the resulting correlation energies to independent benchmarks.

    Authors: We agree that explicit quantitative error metrics would provide a more rigorous and transparent validation of the multipole-Padé representation. Although the manuscript demonstrates agreement through figures comparing the approximant to ab initio data across the Brillouin zone and reports accurate correlation energies, we did not include tabulated error measures. In the revised manuscript we will add, for each of the seven materials, tables (or supplementary tables) reporting the mean absolute deviation and maximum relative error of the inverse dielectric function over a dense grid of momenta spanning the full Brillouin zone. We will also include direct numerical comparisons of the resulting correlation energies against independent ab initio benchmarks. revision: yes

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper constructs symmetry-conserving anisotropic multipole-Padé representations of the inverse dielectric function by fitting parameters to external ab initio data for seven 2D metals. The central claim—that a small number of dispersive modes suffices for accurate description across the Brillouin zone and accurate correlation energies—is an empirical validation result obtained by comparing the fitted representation's outputs to the same ab initio benchmarks. No step reduces a claimed prediction to its own fitted inputs by construction, no uniqueness theorem is imported from self-citations, and no ansatz is smuggled via prior work. The derivation chain is self-contained as a data-driven compact representation whose accuracy is externally falsifiable.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the existence of a compact multipole-Padé form that faithfully reproduces ab initio dielectric data; the parameters of that form are determined by fitting and therefore constitute free parameters. No new physical entities are introduced.

free parameters (1)
  • multipole poles and residues
    A small number of frequency- and momentum-dependent parameters are adjusted to match the computed inverse dielectric function for each material.
axioms (1)
  • domain assumption The inverse dielectric function admits an accurate representation by a finite multipole-Padé approximant that preserves crystal symmetry and anisotropy across the Brillouin zone.
    This is the modeling premise that enables the compact analytic form.

pith-pipeline@v0.9.1-grok · 5680 in / 1273 out tokens · 20504 ms · 2026-06-27T09:04:50.700952+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

35 extracted references

  1. [1]

    E. H. Hwang and S. Das Sarma, Dielectric function, screening, and plasmons in two-dimensional graphene, Phys. Rev. B75, 205418 (2007)

    2007

  2. [2]

    Pfn¨ ur, C

    H. Pfn¨ ur, C. Tegenkamp, and L. Vattuone, Plasmons in one and two dimensions, inSpringer Handbook of Surface Science, edited by M. Rocca, T. S. Rahman, and L. Vat- tuone (Springer International Publishing, Cham, 2020) pp. 557–584

    2020

  3. [3]

    F. H. da Jornada, L. Xian, A. Rubio, and S. G. Louie, Universal slow plasmons and giant field enhancement in atomically thin quasi-two-dimensional metals, Nature Communications11, 1013 (2020)

    2020

  4. [4]

    Muniain and V

    U. Muniain and V. M. Silkin, Impact of the energy dis- persion anisotropy on the plasmonic structure in a two- dimensional electron system, Phys. Chem. Chem. Phys. 24, 17885 (2022)

    2022

  5. [5]

    Cardoso, Z

    C. Cardoso, Z. Kandemir, P. D’Amico, G. Sesti, K. m. c. S ¸endur, M. V. Miloˇ sevi´ c, and C. Sevik, Many-body ef- fects and excitonic corrections in the optical response of two-dimensional metallic mxenes, Phys. Rev. B113, 125131 (2026)

    2026

  6. [6]

    Jablan, M

    M. Jablan, M. Soljaˇ ci´ c, and H. Buljan, Plasmons in graphene: Fundamental properties and potential appli- cations, Proceedings of the IEEE101, 1689 (2013)

    2013

  7. [7]

    Y. Chen, Z. Fan, Z. Zhang, W. Niu, C. Li, N. Yang, B. Chen, and H. Zhang, Two-dimensional metal nanoma- terials: Synthesis, properties, and applications, Chemical Reviews118, 6409 (2018), pMID: 29927583

    2018

  8. [8]

    Fan, N.-H

    Y. Fan, N.-H. Shen, F. Zhang, Q. Zhao, H. Wu, Q. Fu, Z. Wei, H. Li, and C. M. Soukoulis, Graphene plasmonics: A platform for 2d optics, Advanced Optical Materials7, 1800537 (2019)

    2019

  9. [9]

    T. Wang, M. Park, Q. Yu, J. Zhang, and Y. Yang, Sta- bility and synthesis of 2d metals and alloys: a review, Materials Today Advances8, 100092 (2020)

    2020

  10. [10]

    Yousaf, M

    A. Yousaf, M. S. Gilliam, S. L. Y. Chang, M. Augustin, Y. Guo, F. Tahir, M. Wang, A. Schwindt, X. S. Chu, D. O. Li, S. Kale, A. Debnath, Y. Liu, M. D. Green, E. J. G. Santos, A. A. Green, and Q. H. Wang, Exfo- liation of quasi-two-dimensional nanosheets of metal di- borides, The Journal of Physical Chemistry C125, 6787 (2021)

    2021

  11. [11]

    H. T. B. Do, M. Zhao, P. Li, Y. W. Soh, J. Rangaraj, B. Liu, T. Jiang, X. Zhang, J. Lu, P. Song, J. Teng, and M. Bosman, Slow and highly confined plasmons observed in atomically thin tas2, Nature Communications16, 5801 (2025)

    2025

  12. [12]

    R. M. Martin, L. Reining, and D. M. Ceperley,Interact- ing Electrons(Cambridge University Press, Cambridge, 2016)

    2016

  13. [13]

    D. A. Leon and K. Berland, Momentum- and frequency- resolved collective electronic excitations in solids: in- sights from spectroscopy and first-principles calculations, Journal of Physics: Condensed Matter38, 183001 (2026)

    2026

  14. [14]

    D. A. Leon, C. Cardoso, T. Chiarotti, D. Varsano, E. Molinari, and A. Ferretti, Frequency dependence in gwmade simple using a multipole approximation, Phys. Rev. B104, 115157 (2021)

    2021

  15. [15]

    D. A. Leon, A. Ferretti, D. Varsano, E. Molinari, and C. Cardoso, Efficient full frequency gw for metals using a multipole approach for the dielectric screening, Phys. Rev. B107, 155130 (2023)

    2023

  16. [16]

    D. A. Leon, K. Berland, and C. Cardoso, Spectral prop- erties from an efficient analytical representation of the 6 gwself-energy within a multipole approximation, Phys. Rev. B111, 195147 (2025)

    2025

  17. [17]

    D. A. Leon, C. Cardoso, and K. Berland, Bulk plasmons in elemental metals, Phys. Rev. B113, 125138 (2026)

    2026

  18. [18]

    H. Wang, X. Huang, J. Lin, J. Cui, Y. Chen, C. Zhu, F. Liu, Q. Zeng, J. Zhou, P. Yu, X. Wang, H. He, S. H. Tsang, W. Gao, K. Suenaga, F. Ma, C. Yang, L. Lu, T. Yu, E. H. T. Teo, G. Liu, and Z. Liu, High- quality monolayer superconductor nbse2 grown by chem- ical vapour deposition, Nature Communications8, 394 (2017)

    2017

  19. [19]

    S. Guan, S. A. Yang, L. Zhu, J. Hu, and Y. Yao, Electronic, dielectric and plasmonic properties of two- dimensional electride materials x2n (x=ca, sr): A first- principles study, Scientific Reports5, 12285 (2015)

    2015

  20. [20]

    Cudazzo and M

    P. Cudazzo and M. Gatti, Collective charge excitations of the two-dimensional electride ca 2N, Phys. Rev. B96, 125131 (2017)

    2017

  21. [21]

    Eshuis, J

    H. Eshuis, J. E. Bates, and F. Furche, Electron correla- tion methods based on the random phase approximation, Theoretical Chemistry Accounts131, 1084 (2012)

    2012

  22. [22]

    X. Ren, P. Rinke, C. Joas, and M. Scheffler, Random- phase approximation and its applications in computa- tional chemistry and materials science, Journal of Mate- rials Science47, 7447 (2012)

    2012

  23. [23]

    G. D. Mahan,Many-Particle Physics(Springer New York, NY, 2013)

    2013

  24. [24]

    Giannozzi, S

    P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococ- cioni, I. Dabo, A. D. Corso, S. de Gironcoli, S. Fabris, G. Fratesi, R. Gebauer, U. Gerstmann, C. Gougoussis, A. Kokalj, M. Lazzeri, L. Martin-Samos, N. Marzari, F. Mauri, R. Mazzarello, S. Paolini, A. Pasquarello, L. Paulatto, C. Sbraccia, S. Sc...

    2009

  25. [25]

    Giannozzi, O

    P. Giannozzi, O. Andreussi, T. Brumme, O. Bunau, M. B. Nardelli, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, M. Cococcioni, N. Colonna, I. Carnimeo, A. D. Corso, S. de Gironcoli, P. Delugas, R. A. DiStasio, A. Ferretti, A. Floris, G. Fratesi, G. Fugallo, R. Gebauer, U. Gerstmann, F. Giustino, T. Gorni, J. Jia, M. Kawa- mura, H.-Y. Ko, A. Kokalj, E. K¨...

    2017

  26. [26]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett. 77, 3865 (1996)

    1996

  27. [27]

    D. R. Hamann, Optimized norm-conserving vanderbilt pseudopotentials, Phys. Rev. B88, 085117 (2013)

    2013

  28. [28]

    Marini, C

    A. Marini, C. Hogan, M. Gr¨ uning, and D. Varsano, yambo: An ab initio tool for excited state calculations, Comput. Phys. Commun.180, 1392 (2009)

    2009

  29. [29]

    Sangalli, A

    D. Sangalli, A. Ferretti, H. Miranda, C. Attaccalite, I. Marri, E. Cannuccia, P. Melo, M. Marsili, F. Paleari, A. Marrazzo, G. Prandini, P. Bonf` a, M. O. Atambo, F. Affinito, M. Palummo, A. Molina-S´ anchez, C. Hogan, M. Gr¨ uning, D. Varsano, and A. Marini, Many-body per- turbation theory calculations using the yambo code, J. Phys.: Condens. Matter31, 3...

    2019

  30. [30]

    Mounet, M

    N. Mounet, M. Gibertini, P. Schwaller, D. Campi, A. Merkys, A. Marrazzo, T. Sohier, I. E. Castelli, A. Ce- pellotti, G. Pizzi, and N. Marzari, Two-dimensional ma- terials from high-throughput computational exfoliation of experimentally known compounds, Nature Nanotech- nology13, 246 (2018)

    2018

  31. [31]

    Campi, N

    D. Campi, N. Mounet, M. Gibertini, G. Pizzi, and N. Marzari, Expansion of the materials cloud 2d database, ACS Nano17, 11268 (2023), pMID: 37310789

    2023

  32. [32]

    See Supplemental Materials for a detailed description

  33. [33]

    Guandalini, D

    A. Guandalini, D. A. Leon, P. D’Amico, C. Cardoso, A. Ferretti, M. Rontani, and D. Varsano, Efficient GW calculations via interpolation of the screened interac- tion in momentum and frequency space: The case of graphene, Phys. Rev. B109, 075120 (2024)

    2024

  34. [34]

    Sesti, A

    G. Sesti, A. Guandalini, A. Ferretti, P. D’Amico, C. Car- doso, M. Rontani, and D. Varsano, Efficient gw calcula- tions for metals from an accurate ab initio polarizability, https://arxiv.org/abs/2508.06930(2025)

    arXiv 2025

  35. [35]

    D. A. Leon, (2026),https://github.com/ DarioALeonValido/MPAq_2D-metals

    2026