REVIEW 6 minor 53 references
Separation between long- and short-range part of the Coulomb interaction in low dimensional systems: implications for the macroscopic screening length and the collective charge excitations
T0 review · 0 major / 6 minor · reviewed 2026-08-09 · deepseek-v4-flash
Pith's one-line read In low-dimensional materials, a redefined split of the Coulomb interaction restores a scalar macroscopic dielectric function whose zero crossings identify collective charge excitations and whose static limit yields the screening length.
desk verdict Solid and useful: gives a controlled scalar Dyson equation for the 2D macroscopic response with honest validity conditions, though the 'formally exact LRA' claim is a bit stronger than what the factorization strictly establishes. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central device is the local response approximation (LRA): the ansatz that the microscopic response functions factor as $\chi_{00}(q,r_\perp,r'_\perp;\omega)=f(q,\omega)|\varphi(r_\perp)|^2|\varphi(r'_\perp)|^2$ and the same form for $\bar\chi$, with $\varphi$ the transverse wave function confined to a thickness $\lambda$. This factorization collapses the integral Dyson equation for the head of the response into a scalar equation for $\chi^{\ell\mathrm{D}}_M$ in terms of $\bar\chi^{\ell\mathrm{D}}_M$ and the effective interaction $v^{\ell\mathrm{D}}_{\mathrm{eff}}(q)$. In 2D the ansatz is formally exact in the regime $|q|\lambda\ll1$, where the exponential in the long-range potential can be replaced by one; in 1D the divergence of $K_0$ at zero makes the LRA an approximation and leaves a full microscopic treatment as the only exact route.
What would settle it
Compute the full microscopic response of a 2D material at wave vectors where $|q|\lambda$ is not small (for example, a three-layer system at $q\gtrsim0.1$ a.u.${}^{-1}$) without using the LRA factorization, and compare the loss peaks with the zero crossings of $\Re[1-v^{2\mathrm{D}}_{\mathrm{eff}}\bar\chi^{2\mathrm{D}}_M]$; a mismatch would show the scalar criterion holds only in the strictly long-wavelength regime. In 1D, repeat the LRA-based dielectric calculation for several transverse widths $\lambda$ and check whether the predicted collective-mode energies drift with the cutoff, as the divergence of $K_0$ implies.
Extended reading notes
Core claim
In bulk crystals the macroscopic dielectric function is the head of the inverse dielectric matrix, and collective longitudinal modes appear as zeros of its real part. The paper claims that the same structure survives in one and two dimensions once the long-range Coulomb interaction is redefined as the $\mathbf{G}_\parallel=0$ part, whose partial Fourier transform is $2\pi |q|^{-1}e^{-|q||r_\perp-r'_\perp|}$ in 2D and $2K_0(|q||r_\perp-r'_\perp|)$ in 1D. Under the local response approximation, the low-dimensional macroscopic response $\chi^{\ell\mathrm{D}}_M$ is tied to the modified response $\bar\chi^{\ell\mathrm{D}}_M$ by a scalar Dyson equation, which introduces $\varepsilon^{\ell\mathrm{D}}_M=1-v^{\ell\mathrm{D}}_{\mathrm{eff}}(q)\bar\chi^{\ell\mathrm{D}}_M$. On NbSe$_2$, the prominent 0.72 eV loss feature is identified as an intraband plasmon because $\Re\varepsilon^{2\mathrm{D}}_M$ vanishes there; in BN, the lowest bright exciton is also a zero of $\Re\varepsilon^{2\mathrm{D}}_M$, showing its collective character. The paper further derives $r_0=-\lim_{q\to0}(2\pi/q^2)\bar\chi^{2\mathrm{D}}_M(q,0)$ and obtains 11.2 a.u. for BN and 76.0 a.u. for MoS$_2$.
Load-bearing premise
The argument rests on the local response approximation, which assumes the microscopic charge response can be written as a single function of $q$ and $\omega$ times the transverse density profile $|\varphi(r_\perp)|^2$ for both the full and the modified response; this is exact only in 2D for wavelengths much longer than the material thickness.
Editorial extensions
If this is right
- For 2D semiconductors, the screening length becomes an ab initio quantity through $r_0=-\lim_{q\to0}(2\pi/q^2)\bar\chi^{2\mathrm{D}}_M(q,0)$, with the paper reporting 11.2 a.u. for BN and 76.0 a.u. for MoS$_2$.
- For 2D metals, the metallic screening wave vector follows from $q_0=-\lim_{q\to0}2\pi\bar\chi^{2\mathrm{D}}_M(q,0)$, giving 7.1 a.u.${}^{-1}$ for NbSe$_2$.
- A peak in a 2D loss spectrum can be classified as collective by checking whether $\Re\varepsilon^{2\mathrm{D}}_M(q,\omega)$ vanishes; this identifies the 0.72 eV NbSe$_2$ feature as an intraband plasmon and the lowest bright BN exciton as a collective mode.
- The scalar Dyson structure is independent of the approximation used for $\bar\chi$, so the same criterion applies at RPA, BSE, or higher levels of theory.
- In 1D systems no exact macroscopic dielectric function exists; finite-thickness effects must always be included, and any LRA-based description is only approximate.
Reading between the lines
- If the scalar Dyson equation is used with exchange-correlation kernels beyond RPA/BSE, the zero-crossing criterion could track plasmon softening or exciton-plasmon mixing in 2D materials without phenomenological models; the paper establishes the framework but does not run such calculations.
- The $r_0$ and $q_0$ formulas give a common scale for comparing macroscopic screening across 2D materials; one could test whether these lengths correlate with exciton binding energies or with screening by a surrounding dielectric environment.
- The 1D analysis implies that effective-interaction models of 1D plasmons and excitons must carry an explicit transverse cutoff; a numerical scan of the $|q|\ln(\lambda|q|)$ dispersion as $\lambda$ varies would quantify how strongly the mode energies depend on that regularization.
- For finite-thickness 2D systems the paper's own spectra show small high-frequency residuals, suggesting the LRA can be extended by using a model effective interaction with a finite transverse profile while keeping the scalar Dyson equation; the paper notes this as a regime extension, not a contradiction.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript addresses the long-standing difficulty of defining a macroscopic dielectric function in low-dimensional (ℓD) systems. The authors propose to separate the Coulomb interaction into a long-range part containing only reciprocal vectors parallel to the periodic directions (G‖ = 0) and a short-range part containing all remaining components. Using this split in the microscopic Dyson equation, they derive a scalar Dyson-like equation for the ℓD macroscopic response function χ_M^ℓD (Eq. 19) with an effective long-range interaction v_eff (Eq. 20), and define an ℓD macroscopic dielectric function ε_M^ℓD. Collective charge excitations are then identified with zeros of Re ε_M^ℓD, in analogy with the 3D case. They also express the 2D screening length r0 and the Thomas-Fermi wavevector q0 in terms of the modified response function (Eqs. 23–24), and validate the approach with RPA and BSE calculations for NbSe2, BN, and MoS2, including comparisons between the full microscopic solution and the local response approximation.
Significance. If the derivation stands, the paper provides a general and parameter-free first-principles route to interpret EELS/IXS spectra of low-dimensional materials and to extract macroscopic screening lengths, replacing phenomenological models. The central formal step—the scalar Dyson equation—is controlled in the regime |q|λ << 1, and the numerical comparisons in Figs. 2 and 3 are convincing. The paper is also transparent about the limitations in 1D and for finite-thickness 2D systems. The main strengths are the explicit connection between microscopically accessible response functions and macroscopic spectroscopic quantities, and the consistency of the resulting collective-mode criterion with the bulk analogy.
minor comments (6)
- [Sections III and IV, after Eq. (16) and first paragraph of Section IV] The statement that the LRA is 'formally exact' for |q|λ << 1 is stronger than what is proved: the factorized ansatz of Eqs. (17–18) is not exact in general, although the macroscopic scalar equation (Eq. 19) does follow exactly in that regime because v0^2D is constant in the transverse coordinates. Please revise the wording to distinguish the exact macroscopic closure from the approximate microscopic factorization.
- [Eq. (21)] The Gaussian normalization appears to be incorrect: for ℓ = 2, the normalized transverse wavefunction should be (2/(πλ^2))^{1/4} e^{-z^2/λ^2}, whereas the printed form (2/(πλ))^{1/4} is not normalized, which conflicts with the statement that Eq. (10) implies f = χ_M^ℓD. The same issue affects the ℓ = 1 prefactor. Please correct the prefactors and re-derive the explicit v_eff expressions accordingly.
- [Computational details] The paper does not specify the ground-state code, pseudopotentials, k-grids, cutoffs, or the λ values used for the LRA curves in Figs. 2 and 3; for a journal version, please add a computational methods paragraph or a supplementary material section.
- [Section III, paragraph after Eq. (21)] The statement that arχ_M^ℓD is proportional to |q|^2 at small wavevectors is only valid at finite frequency; for a metal at ω = 0 the static arχ_M^2D tends to a constant, so the equality lim_{q→0} χ_M = lim_{q→0} arχ_M does not hold in that case. Please add the implied frequency qualification.
- [Fig. 3 and surrounding text] The text says the peak at 5.66 eV in -Im[ε^{-1}] and the peak at 5.28 eV in Im[ε] are 'essentially the same'; please clarify that these are the longitudinal and transverse exciton positions related through the long-range interaction, since the 0.38 eV difference is visible in the figure.
- [Throughout] The manuscript contains many typos and grammatical errors, including 'Coulom b' in the title, 'ifferent', 'expres sion', 'well behaves', 'at last at small q', and a malformed double bar in Eq. (18). A careful language edit is needed.
Circularity Check
No significant circularity: the central 2D result follows from the Dyson decomposition plus an acknowledged ansatz, and the screening length is defined, not fitted.
full rationale
The central derivation is self-contained rather than circular. Starting from the exact Dyson decomposition in Section III (Eqs. 13-15), the paper writes a scalar Dyson equation (Eq. 19) only after introducing the Local Response Approximation (Eqs. 17-18) explicitly as an ansatz based on transverse localization of the response; this is an approximation whose validity is discussed (|q|lambda << 1 in 2D, always approximate in 1D), not a hidden restatement of the target result. The collective-excitation criterion (zeros of Re epsilon_2D_M) is then a standard pole condition on the scalar equation, exactly as in bulk, so it is a consequence rather than an input. The screening length r0 in Eq. 23 is obtained by matching Eq. 22 to the asymptotic form 1 + r0 |q|; the resulting expression r0 = lim_{q->0} -2pi/q^2 bar-chi_2D_M(q,0) is a definition/identity evaluated from the computed response, not a parameter fitted to the quantity being predicted. The self-citations (Refs. 41, 44, 47, 51, 52) are used as external benchmarks and standard asymptotic results (plasmon and exciton dispersions, 2D static dielectric behavior) and are not load-bearing for the derivation. The only caveat is that the statement that the LRA is 'formally exact' in 2D is stronger than what the factorization argument proves, but this is a mathematical-rigor concern about an approximation, not circularity.
Assumptions & free parameters
assumptions (4)
- standard math The microscopic charge response function obeys the standard Dyson equation with the Coulomb interaction split into long- and short-range parts.
- domain assumption The low-dimensional macroscopic average is defined by Eq. 6 with Omega_perp -> infinity, which converges because electronic wave functions are confined along the non-periodic directions.
- domain assumption The appropriate long/short-range separation of the Coulomb interaction in low-dimensional systems is Eq. 12: the long-range part includes all G_perp with G_parallel = 0, and the short-range part excludes G_parallel = 0.
- ad hoc to paper The Local Response Approximation ansatz, Eqs. 17-18, factorizes chi and bar-chi as f(q,omega)|phi(r_perp)|^2 |phi(r'_perp)|^2.
Cite this review
Pith. "Pith review of Separation between long- and short-range part of the Coulomb interaction in low dimensional systems: implications for the macroscopic screening length and the collective charge excitations." pith.science (2026). https://pith.science/paper/EBD355LD
@misc{pith2026250203181,
author = {Pith},
title = {Pith review of: Separation between long- and short-range part of the Coulomb interaction in low dimensional systems: implications for the macroscopic screening length and the collective charge excitations},
year = {2026},
howpublished = {\url{https://pith.science/paper/EBD355LD}},
note = {Machine review of arXiv:2502.03181}
}
read the original abstract
Collective charge excitations directly probed in electron energy loss and inelastic X rays scattering spectroscopies play a key role in different fields of condensed matter physics. Being induced by the long-range part of the Coulomb interaction between particles, in standard bulk systems they appear as well-defined features in the spectra associated to the zero crossing of the real part of the macroscopic dielectric function. However, this simple criterion cannot be used to identify collective excitations in low dimensional systems where the macroscopic dielectric function is not a well defined concept. In this work, we discuss how this problem can be traced back to the definition of the long-range Coulomb interaction and we show how the appropriate separation between long- and short-range Coulomb interaction allows one to correctly express the low dimensional macroscopic dielectric function in terms of microscopic quantities accessible in first-principles calculations. This allows disentangling collective charge excitations in low dimensional materials in analogy with what one does in standard bulk systems. In addition, we show how important macroscopic quantities, such as the screening length scale, can be extracted from full first principle calculations. As an illustrative example we perform a study of the screening effects and the collective excitations in prototypical 2D materials including metals (NbSe2) as well as semiconductors (BN and MoS2).
Figures
Reference graph
Works this paper leans on
- [1]
- [2]
- [3]
- [4]
- [5]
-
[6]
C. Lian, Z. A. Ali, and B. M. Wong, Phys. Rev. B 100 , 205423 (2019)
work page 2019
-
[7]
D. Varsano, S. Sorella, D. Sangalli, M. Barborini, S. Corni, E. Molinari, and M. Rontani, Nat. Commun. 8 , 1461 (2017)
work page 2017
-
[8]
M. Barborini, M. Calandra, F. Mauri, L. Wirtz, and P. Cudazzo, Phys. Rev. B 105 , 075122 (2022)
work page 2022
Show all 53 references
-
[9]
Kohn, Phys
W. Kohn, Phys. Rev. Lett. 19 , 439 (1967)
1967
-
[10]
L. V. Keldysh and Y. V. Kopaev, Sov. Phys. Solid State. 6 , 2219 (1965)
1965
-
[11]
J\'erome, T
D. J\'erome, T. M. Rice, and W. Kohn, Phys. Rev. 158 , 462 (1967)
1967
-
[12]
O. V. Dolgov, D. A. Kirzhnits, and E. G. Maksimov, Rev. Mod. Phys. 53 , 81 (1981)
1981
-
[13]
Ichimaru, Rev
S. Ichimaru, Rev. Mod. Phys. 54 , 1017 (1982)
1982
-
[14]
Takada, Journal of Superconductivity 18 , 785 (2005)
Y. Takada, Journal of Superconductivity 18 , 785 (2005)
2005
-
[15]
Takada, Physical Review B 94 , 245106 (2016)
Y. Takada, Physical Review B 94 , 245106 (2016)
2016
- [16]
-
[17]
Davis, H
D. Davis, H. Gutfreund, and W. A. Little, Phys. Rev. B 13 , 4766 (1976)
1976
-
[18]
J. Ihm, M. L. Cohen, and S. F. Tuan, Phys. Rev B 23 , 3258 (1981)
1981
-
[19]
Y. W. Choi, J. Ihm, and M. L. Cohen, Phys. Rev. B 110 , 155127 (2024)
2024
-
[20]
A. O. Govorov, and H. H. Richardson, Nano Today 2 , 30 (2007)
2007
-
[21]
Baffou, and R
G. Baffou, and R. Quidant, Laser Photon. Rev. 7 , 171 (2013)
2013
-
[22]
M. L. Brongersma, N. J. Halas, and P. Nordlander, Nat. Nanotech. 10 , 25 (2015)
2015
-
[23]
Yang, et al., Nano Lett
X. Yang, et al., Nano Lett. 7 , 3798 (2007)
2007
-
[24]
Chen, et al., Nano Lett
Y. Chen, et al., Nano Lett. 11 , 348 (2011)
2011
-
[25]
J. Lee, S. Mubeen, X. Ji, G. D. Stucky, and M. Moskovits, Nano Lett. 12 , 5014 (2012)
2012
-
[26]
X. Zhou, G. Liu, J. Yu, and W. Fan, J. Mater. Chem. 22 , 21337 (2012)
2012
-
[27]
J. G. Smith, J. A. Faucheaux, and P. K. Jain, Nano Today 10 , 67 (2015)
2015
-
[28]
Plasmonics: Fundamental and Applications
S. A. Maier "Plasmonics: Fundamental and Applications" (2007) (New York: Springer)
2007
-
[29]
Fei, et al., Nature 487 , 82 (2012)
Z. Fei, et al., Nature 487 , 82 (2012). F. H. L. Koppens, D. E. Chang, F. J. Garc\'ia de Abajo, Nano Lett. 11 , 3370 (2011)
2012
-
[30]
F. J. Garc\'ia de Abajo, Science 339 , 917 (2013)
2013
-
[31]
Sottile, F
F. Sottile, F. Bruneval, A. G. Marinopoulos, L. K. Dash, S. Botti, V. Olevano, N. Vast, A. Rubio, and L. Reining, Int. J. Quantum Chem. 102 , 684 (2005)
2005
-
[32]
J. M. Pitarke, V. M .Silkin, E. V. Chulkov, and P. M. Echenique, Rep. Prog. Phys. 70 , 1-87 (2007)
2007
-
[33]
Mazzei, and C
S. Mazzei, and C. Giorgetti, Phys. Rev. B 106 , 035431 (2022)
2022
-
[34]
Wiser, Phys
N. Wiser, Phys. Rev. 129 , 62 (1963)
1963
-
[35]
S. L. Adler, Phys. Rev. 126 , 413 (1962)
1962
-
[36]
D. A. Dahland, and L. J. Sham, Phys. Rev. B 17 , 641 (1977)
1977
-
[37]
Hogan, M
C. Hogan, M. Palummo, and R. D. Sole: Comptes Rendus Physique, 10 , 560 (2009)
2009
-
[38]
Sangalli et al., J
D. Sangalli et al., J. Phys.: Condens. Matter 31 325902 (2019)
2019
-
[39]
Mazzei, and C
S. Mazzei, and C. Giorgetti, Phys. Rev. B 107 , 165412 (2023)
2023
-
[40]
Handbook of Mathematical Functions With Formulas, Graphs, and Mathematical Tables
"Handbook of Mathematical Functions With Formulas, Graphs, and Mathematical Tables", edited by M. Abramowitz and I. A. Stegun, NBS Applied Mathematics Series 55 (National Bureau of Standards, Washington, 1964)
1964
-
[41]
Koskelo, G
J. Koskelo, G. Fugallo, M. Hakala, M. Gatti, F. Sottile, and P. Cudazzo, Phys. Rev. B 95 , 035125 (2017)
2017
-
[42]
Cervantes-Villanueva, F
J. Cervantes-Villanueva, F. Paleari, A. Garc\'ia-Crist\'obal, D. Sangalli, and A. Molina-S\'anchez, Phys. Rev. B 109 , 155133 (2024)
2024
-
[43]
Giuliani, and G
G. Giuliani, and G. Vignale, "Quantum Theory of the Electron Liquid (Cambridge: Cambridge University Press) (2005)
2005
-
[44]
Cudazzo, M
P. Cudazzo, M. Gatti, and A. Rubio, New Journal of Physics 15 , 125005 (2013)
2013
-
[45]
Andersen, and K
K. Andersen, and K. S. Thygesen, Phys. Rev. B 88 , 155128 (2013)
2013
-
[46]
F. H. da Jornada, L. Xian, A. Rubio, and S. G. Louie, Nature Comm. 11 , 1013 (2020)
2020
-
[47]
Cudazzo, L
P. Cudazzo, L. Sponza, C. Giorgetti, L. Reining, F. Sottile, and M. Gatti, Phys. Rev. Lett. 116 , 066803 (2016)
2016
-
[48]
D. Y. Qiu, T. Cao, and S. G. Louie, Phys. Rev. Lett. 115 , 176801 (2015)
2015
-
[49]
D. Y. Qiu, G. Cohen, D. Novichkova, and S. Refaely-Abramson, Nano letters 21 , 7644 (2021)
2021
-
[50]
Latini, T
S. Latini, T. Olsen, and K. S. Thygesen, Phys. Rev. B 92 , 245123 (2015)
2015
-
[51]
Cudazzo, I
P. Cudazzo, I. V. Tokatly, and A. Rubio, Phys. Rev. B 84 , 085406 (2011)
2011
-
[52]
Cudazzo, C
P. Cudazzo, C. Attacalite, I. V. Tokatly, and A. Rubio, Phys. Rev. Lett. 104 , 226804 (2010)
2010
-
[53]
T. Ando, A. B. Fowler, and F. Stern, Rev. Mod. Phys. 54 , 437 (1982)
1982
Reviewed August 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.