Pith. sign in

REVIEW 3 major objections 4 minor 33 references

Alkali Intercalation of Moire Heterostructures for Low-Loss Plasmonics

T0 review · 3 major / 4 minor · reviewed 2026-08-15 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Sodium intercalation makes a Moiré hBN heterostructure host plasmons that evade first-order electron-phonon decay, with lifetimes set by electron-electron scattering near 10^14 s^-1.

desk verdict Predicts a concrete candidate for first-order lossless plasmons in sodium-intercalated Moiré hBN, but the e-e decay that sets the real lifetime rests on an unvalidated approximation. read the letter →

arxiv 2505.10225 v1 pith:HORNKI2N submitted 2025-05-15 cond-mat.mtrl-sci cond-mat.mes-hall

Ali Ghorashi , Nicholas Rivera , Ravishankar Sundararaman , Efthimios Kaxiras , John Joannopoulos , Marin Soljačić This is my paper · ORCID
classification cond-mat.mtrl-scicond-mat.mes-hall PACS 73.20.Mf71.15.Mb63.22.-m
keywords Moiréheterostructureshexagonalboronnitridealkaliintercalationplasmonslow-lossplasmonicselectron-phononinteractionflatbandselectron-electronscattering
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper proposes a concrete material platform—sodium atoms intercalated into a √7×√7 Moiré heterostructure of hexagonal boron nitride—that should support plasmons with deep subwavelength confinement that are immune to the usual decay channel of one-phonon emission in a frequency window near 1 eV. If the calculations hold, this would be a practical route to low-loss plasmonics in a structurally stable, nonmagnetic, fabricable system, moving past earlier proposals that were ferromagnetic, unstable, or difficult to make. The key design principle is an isolated, moderately flat band at the Fermi level: with bandwidth 0.9 eV and a maximum phonon energy of 0.18 eV, plasmons above 1.08 eV cannot lose energy by emitting one phonon within that band. The paper shows that higher-order electron-phonon decay is negligible in that window, but that electron-electron interactions cause plasmons to decay to lower-frequency plasmons at a rate around $10^{14}$ $s^{-1}$, which sets the practical lifetime.

What carries the argument

The central object is the isolated, moderately flat electronic band at the Fermi level created by sodium intercalation into the √7×√7 Moiré pattern of hBN. The argument rests on two threshold conditions: a plasmon with energy ℏω is protected from intraband one-phonon decay when ℏω > W + ℏω_ph (with bandwidth W = 0.9 eV and maximum phonon energy ℏω_ph ≈ 0.18 eV), and from interband one-phonon decay when ℏω < |E_{v,c} − E_Fermi| − ℏω_ph. The quantitative machinery includes a three-layer dielectric model for the plasmon dispersion (Eq. 1), the general self-energy-based decay-rate formula of Allen (Eq. 2) that accounts for finite bandwidth and higher-order electron-phonon processes, and the plasmon-pole approximation for the electron-electron self-energy (Eq. 5) that yields the dominant plasmon-plasmon scattering rate. The same framework is validated against the one-phonon formula (Eq. 4) and against graphene, where all three methods agree.

What would settle it

Compute the electron-electron self-energy without the plasmon-pole approximation (for example, with a full RPA dielectric function or a quantum Monte Carlo approach) and compare the resulting plasmon decay rate in the 1.08–1.16 eV window to the paper's ~$10^{14}$ $s^{-1}$ prediction; a rate differing by an order of magnitude would invalidate the central lifetime claim. Alternatively, if a fabricated sodium-intercalated √7×√7 hBN sample shows first-order phonon-assisted plasmon decay below 1.16 eV in optical or near-field spectroscopic measurements, the band-structure picture would be wrong.

Watch

Extended reading notes

Core claim

The central claim is that in the √7×√7 sodium-intercalated hBN Moiré heterostructure, the Moiré potential traps the sodium atoms and produces an isolated band of width 0.9 eV at the Fermi level, with the nearest other bands far enough away that plasmons in the 1.08–1.16 eV window cannot decay through any one-phonon-assisted intraband or interband transition. For hole dopings of n = 0.1 and n = 0.2, the paper finds a first-order lossless window from 1.08 eV up to 1.16 eV, above which interband one-phonon processes turn on. Using self-energy-based decay formulas, the authors compute higher-order electron-phonon decay rates near $10^{7}$ $s^{-1}$ in that window, and then, using a plasmon-pole approximation for the electron-electron self-energy, compute a dominant plasmon-to-plasmon scattering rate of about $10^{14}$ $s^{-1}$. The conclusion is that, in a clean sample, the plasmon lifetime in this material is controlled not by phonons but by intrinsic electron-electron scattering, and that the material is dynamically stable with no imaginary phonon frequencies.

Load-bearing premise

The prediction that electron-electron scattering sets the plasma decay rate at about $10^{14}$ $s^{-1}$ rests on the plasmon-pole approximation for the electron-electron self-energy; if that approximation is quantitatively wrong, the lifetime and the practical value of the material would change.

Editorial extensions

If this is right

  • Plasmons in the 1.08–1.16 eV window of the sodium-intercalated heterostructure would propagate with decay rates set by electron-electron scattering (~10^14 s^-1), yielding quality factors orders of magnitude higher than conventional two-dimensional metals at comparable confinement.
  • The design rule 'isolated flat band at the Fermi level plus small maximum phonon energy' becomes a systematic search principle: among the five alkalis tested, only sodium produces a sufficiently isolated band, and the √7×√7 angle is the one that avoids ferromagnetism.
  • The decay-rate evaluation protocol (combining finite-bandwidth self-energy formulas with one-phonon checks) provides a template for predicting plasmon lifetimes in other bandwidth-limited two-dimensional metals.
  • The lossless window is bounded: plasmons above 1.16 eV acquire interband one-phonon decay channels, so applications must deliberately operate inside the 1.08–1.16 eV range.
  • If realized, the material would offer near-infrared subwavelength plasmonics in a stable, nonmagnetic van der Waals heterostructure, a regime previously inaccessible in hBN-based systems.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The paper's own caveat is that the dominant decay rate rests on the plasmon-pole approximation; if a full random-phase-approximation or diagrammatic Monte Carlo treatment shifts that rate significantly, the practical lifetime and the claim of 'low-loss' would change, though the first-order phonon immunity would remain.
  • The Moiré-potential trapping mechanism suggests a broader design space: varying the intercalant size, Moiré angle, and substrate layers could tune the bandwidth and the gap to neighbouring bands, potentially opening lossless windows at other frequencies.
  • The 1.08–1.16 eV window sits near the telecom and near-infrared region, so a successful experimental realization could enable low-loss nanophotonic devices, but the fabrication of large-angle Moiré hBN with controlled sodium doping remains an untested experimental challenge.
  • The authors mention impurity and two-plasmon scattering as additional channels; the predicted 10^14 s^-1 rate therefore represents an upper bound on lifetime in ideal clean samples, and real devices would likely show additional losses.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. This manuscript proposes a class of Moiré heterostructures of hexagonal boron nitride intercalated with alkali atoms as platforms for low-loss plasmons. The authors screen five alkali species and two Moiré angles, and identify sodium intercalated in a √7×√7 hBN Moiré structure as the only candidate with an isolated, moderately flat band at the Fermi level. Using DFT for band structures and phonons and a three-layer model for the plasmon dispersion, they find that for hole dopings n=0.1 and n=0.2, plasmons in the energy window 1.08–1.16 eV are immune to first-order electron-phonon decay. They then compute higher-order electron-phonon decay rates via Eq. (2) and find rates around 10^7 s^-1 in that window. Finally, using the plasmon-pole approximation for the electron-electron self-energy, Eq. (5), they report that plasmon-to-plasmon scattering dominates, with a rate around 10^14 s^-1. The paper concludes that the electron-phonon channel is effectively suppressed but the electron-electron channel sets the practical lifetime.

Significance. If the quantitative predictions hold, the paper provides a concrete, structurally stable, nonmagnetic material in which a specific plasmon frequency window is immune to first-order electron-phonon decay, addressing a limitation of previous defect-lattice proposals in hBN. The kinematic argument for the lossless window is transparent and rests on computed band and phonon properties rather than on parameters fitted to the target decay rates. The comparison of Eqs. (2), (3), and (4) offers a useful internal consistency check for the electron-phonon channel, and the prediction of a 1.08–1.16 eV lossless window is falsifiable. The main weakness is that the dominant loss channel, the electron-electron plasmon-plasmon scattering rate of about 10^14 s^-1, is obtained solely from the plasmon-pole approximation and is not validated by an independent calculation or a quantitative error estimate. Because this channel governs the practical lifetime, the central low-loss claim is not yet fully established.

major comments (3)
  1. [Plasmon to Plasmon Scattering, Eq. (5), Fig. 4] The dominant decay rate of about 10^14 s^-1 is computed by inserting the plasmon-pole electron-electron self-energy, Eq. (5), into Eq. (2). This is the only calculation of the electron-electron channel in the paper, and the Outlook explicitly states that the prediction rests on the validity of the plasmon-pole approximation. Since this channel dominates the lifetime in the 1.08–1.16 eV window, the central practical claim of controlled losses is not yet quantitatively established. I recommend adding an independent validation, such as a full RPA calculation of the loss function, a comparison with the f-sum rule, or a quantitative assessment of the plasmon-pole approximation error for the relevant doping and frequency range. A qualitative caveat in the Outlook is not sufficient for a reported rate of 10^14 s^-1.
  2. [Plasmonic Losses, Eq. (2), Footnote [23]] The manuscript uses only the imaginary part of the electron self-energy and states in footnote [23] that the real part contributes to a change in the effective bandwidth. The first-order lossless window, however, is defined by the computed bandwidth (0.9 eV), the maximum phonon energy, and the band-edge position relative to the Fermi level. A real-part correction of a few tens of meV could shift the 1.08 eV lower edge or the 1.16 eV upper edge for n=0.1 and materially change the predicted window. The authors should quantify the real-part contribution or argue with specific numbers that it is negligible.
  3. [Introduction, Fig. 2, lossless window] For n=0.1, the first-order lossless window spans only 1.08 to 1.16 eV, a margin of about 80 meV. The paper reports no error bars or convergence tests for the DFT bandwidth, the maximum phonon frequency, or the band-edge energies that set both edges of this window. Because the window is the manuscript's core prediction, its robustness to computational parameters (k-point sampling, Wannier interpolation, exchange-correlation functional, and van der Waals correction) should be demonstrated. A short convergence or sensitivity study would substantially strengthen the claim.
minor comments (4)
  1. [Introduction] There are typographical errors, including 'subsitutional' instead of 'substitutional' and inconsistent spelling of 'Moiré' as 'Moire'.
  2. [Eq. (1)] The typeset form of Eq. (1) in the manuscript is difficult to parse; please ensure the three-layer plasmon dispersion equation is rendered clearly and that all symbols are defined at first use.
  3. [Fig. 2(c)] The caption for Fig. 2(c) should state explicitly what the pink shaded region represents and what the solid arrows denote, since the current description is ambiguous.
  4. [References] Reference [33] is incomplete; it should include the full author list and a proper citation for the JJDFTX.jl repository.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the low-loss window and decay rates follow from DFT-computed band and phonon properties and standard many-body formulas; no fitted parameter is renamed as a prediction.

full rationale

The paper's central claims are derived from independently computed inputs. The first-order lossless window (1.08–1.16 eV) is obtained from the DFT/Wannier-computed isolated-band bandwidth W = 0.9 eV and the computed maximum phonon energy ≈0.18 eV via the energy-conservation condition ℏω > W + ℏω_ph for intraband one-phonon decay, with an upper bound from the interband threshold; this is a kinematic consequence, not a fit to the reported decay rates. The electron-phonon decay rates are evaluated with the standard lifetime formula of Eq. (2), attributed to Allen [21], and the one-phonon golden-rule expression of Eq. (4); the paper independently determines the validity region ω_0 by comparing Eq. (2) and Eq. (4) in Fig. 3, including a graphene benchmark. The electron-electron rate is computed by inserting the explicitly labeled plasmon-pole approximation, Eq. (5), into Eq. (2); the authors candidly state in the Outlook that this prediction rests on the validity of that approximation. That is an approximation caveat and a correctness risk, not a circular reduction: the self-energy parameters are the computed plasmon dispersion and occupation factors, not the target lifetime. Self-citations [9] and [24] provide methodological context for one-phonon equivalence and an alternative correlation-bubble treatment, but the central claim does not depend solely on those citations because the band structure, phonons, and dielectric response are computed here and the kinematic criteria are stated and applied directly. No parameter is fitted to the predicted loss rates, and no derived quantity is defined in terms of the quantity it is used to predict. The derivation chain is therefore self-contained and non-circular.

Assumptions & free parameters 1 free parameters · 4 assumptions · 0 invented entities

No new particles, mediators, forces, or conserved quantities are introduced. The proposed material structure is a prediction, not an invented fundamental entity, and its experimental evidence is absent by design.

free parameters (1)
  • Hole doping n = n = 0.1 and n = 0.2 (lossless window); n = 0.3 to 0.5 also examined
    Five filling factors are scanned; the lossless plasmon window exists only for n=0.1 and n=0.2. Doping is an experimental control, not fitted to the target result, but the headline prediction is conditional on it.
assumptions (4)
  • domain assumption The isolated band at the Fermi level is the only relevant electronic band for the plasmonic decay channels considered; other bands contribute only through static dielectric screening of hBN.
    Used in the three-layer conductivity model Eq. (1) and in restricting the phonon-assisted decay sums to the isolated band, as stated below Eq. (4).
  • domain assumption A plasmon with energy between W+ω_ph and E_c−E_F−ω_ph cannot decay via a one-phonon intraband or interband process; this defines the first-order lossless window.
    This criterion is taken from prior work [9] by the same group and is used to identify the 1.08 to 1.16 eV window.
  • domain assumption DFT-PBE with D3 corrections and the chosen pseudopotentials gives quantitatively reliable band structures, Fermi-level band isolation, and phonon stability for the screened Moiré structures.
    All electronic and phononic inputs come from this method; no convergence tests or experimental benchmarks are shown.
  • domain assumption The plasmon-pole approximation for the electron-electron self-energy, Eq. (5), yields quantitatively reliable electron-electron decay rates.
    This approximation underlies the reported 10^14 s^-1 decay rate; the Outlook explicitly says the prediction rests on its validity.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Alkali Intercalation of Moire Heterostructures for Low-Loss Plasmonics." pith.science (2026). https://pith.science/paper/HORNKI2N

@misc{pith2026250510225,
  author       = {Pith},
  title        = {Pith review of: Alkali Intercalation of Moire Heterostructures for Low-Loss Plasmonics},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/HORNKI2N}},
  note         = {Machine review of arXiv:2505.10225}
}
abstract

Two-dimensional metals generically support gapless plasmons with wavelengths well below the wavelength of free-space radiation at the same frequency. Typically, however, this substantial confinement of electromagnetic energy is associated with commensurately high losses, and mitigating such losses may only be achieved through judicious band structure engineering near the Fermi level. In a clean system, an isolated, moderately flat, band at the Fermi level with sufficiently high carrier density can support a plasmon that is immune to propagation losses up to some order in the electron-phonon interaction. However, proposed materials that satisfy these criteria have been ferromagnetic, structurally unstable, or otherwise difficult to fabricate. Here, we propose a class of band structure engineered materials that evade these typical pitfalls -- Moire heterostructures of hexagonal boron nitride intercalated with alkali atoms. We find that only sodium atoms engender a sufficiently isolated band with plasmons lossless at first order in the electron-phonon interaction. We calculate higher order electron-phonon losses and find that at frequencies of about $1$eV the electron-phonon decay mechanism is negligible -- leading to a contribution to the decay rate of about 10^7 Hz in a small frequency range. We next calculate losses from the electron-electron interaction and find that this is the dominant process -- leading plasmons to decay to lower frequency plasmons at a rate of around 10^14 Hz.

Figures

Figures reproduced from arXiv: 2505.10225 by the authors.

Figure 1
Figure 1. FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 4
Figure 4. FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figures from the paper (1 more)
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p004_3.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

33 extracted references · 28 canonical work pages

  1. [9]

    Ghorashi, N

    A. Ghorashi, N. Rivera, B. Shi, R. Sundararaman, E. Kaxi- ras, J. Joannopoulos, and M. Solja ˇci´c, Highly confined, low- loss plasmonics based on two-dimensional solid-state defect lattices, Physical Review Materials8, L011001 (2024)

    work page 2024

  2. [23]

    In principle, the real part contributes to a change in the effective bandwidth

    In this work, we consider only the imaginary part of the self energy. In principle, the real part contributes to a change in the effective bandwidth

    work page

  3. [1]

    Pines,Theory of quantum liquids: normal Fermi liquids (CRC Press, 2018)

    D. Pines,Theory of quantum liquids: normal Fermi liquids (CRC Press, 2018)

    work page 2018

  4. [2]

    Bohm and D

    D. Bohm and D. Pines, A collective description of electron in- teractions: Iii. coulomb interactions in a degenerate electron gas, Physical Review92, 609 (1953)

    work page 1953

  5. [3]

    Jablan, M

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

    work page 2013

  6. [4]

    H. A. Atwater and A. Polman, Plasmonics for improved photo- voltaic devices, Nature materials9, 205 (2010)

    work page 2010

  7. [5]

    Wokaun, J

    A. Wokaun, J. Bergman, J. Heritage, A. Glass, P. Liao, and D. Olson, Surface second-harmonic generation from metal is- land films and microlithographic structures, Physical Review B 24, 849 (1981)

    work page 1981

  8. [6]

    S. S. Yee, J. Homola, and G. Gauglitz, Surface plas- mon resonance (spr) optical sensors, current technology and applications-preface (1999)

    work page 1999

Show all 33 references
  1. [7]

    Jablan, H

    M. Jablan, H. Buljan, and M. Solja ˇci´c, Plasmonics in graphene at infrared frequencies, Physical Review B—Condensed Matter and Materials Physics80, 245435 (2009)

    2009

  2. [8]

    Stauber and N

    T. Stauber and N. Peres, Effect of holstein phonons on the elec- tronic properties of graphene, Journal of Physics: Condensed Matter20, 055002 (2008)

    2008

  3. [10]

    J. B. Khurgin and G. Sun, In search of the elusive lossless metal, Applied Physics Letters96(2010)

    2010

  4. [11]

    J. B. Khurgin and A. Boltasseva, Reflecting upon the losses in plasmonics and metamaterials, MRS bulletin37, 768 (2012)

    2012

  5. [12]

    M. N. Gjerding, M. Pandey, and K. S. Thygesen, Band struc- ture engineered layered metals for low-loss plasmonics, Nature communications8, 15133 (2017)

    2017

  6. [13]

    X. Hu, Z. Wu, Z. Li, Q. Xu, K. Chen, W. Hu, Z. Ni, K. Jin, H. Weng, and L. Lu, High-throughput search for lossless met- als, Physical Review Materials6, 065203 (2022)

    2022

  7. [14]

    Beida, E

    W. Beida, E. S ¸as ¸ıo˘glu, M. Tas, C. Friedrich, S. Bl ¨ugel, I. Mer- tig, and I. Galanakis, Correlation effects in two-dimensional mx 2 and ma 2 z 4 (m= nb, ta; x= s, se, te; a= si, ge; z= n, p) cold metals: Implications for device applications, Physical Review Materials9, 0...

    2025

  8. [15]

    Papaj and C

    M. Papaj and C. Lewandowski, Plasmonic nonreciprocity driven by band hybridization in moir ´e materials, Physical Re- view Letters125, 066801 (2020)

    2020

  9. [16]

    Lewandowski and L

    C. Lewandowski and L. Levitov, Intrinsically undamped plas- mon modes in narrow electron bands, Proceedings of the Na- tional Academy of Sciences116, 20869 (2019)

    2019

  10. [17]

    Lon ˇcari´c, Z

    I. Lon ˇcari´c, Z. Rukelj, V . M. Silkin, and V . Despoja, Strong two- dimensional plasmon in li-intercalated hexagonal boron-nitride film with low damping, npj 2D Materials and Applications2, 33 (2018)

    2018

  11. [18]

    Rukelj, Dynamical conductivity of lithium-intercalated hexagonal boron nitride films: A memory function approach, Physical Review B102, 205108 (2020)

    Z. Rukelj, Dynamical conductivity of lithium-intercalated hexagonal boron nitride films: A memory function approach, Physical Review B102, 205108 (2020)

    2020

  12. [19]

    Cassabois, P

    G. Cassabois, P. Valvin, and B. Gil, Hexagonal boron nitride is an indirect bandgap semiconductor, Nature photonics10, 262 (2016)

    2016

  13. [20]

    K. S. Thygesen, Calculating excitons, plasmons, and quasipar- ticles in 2d materials and van der waals heterostructures, 2D Materials4, 022004 (2017)

    2017

  14. [21]

    P. B. Allen, Electron self-energy and generalized drude for- mula for infrared conductivity of metals, Physical Review B 92, 054305 (2015)

    2015

  15. [22]

    Giustino, Electron-phonon interactions from first principles, Reviews of Modern Physics89, 015003 (2017)

    F. Giustino, Electron-phonon interactions from first principles, Reviews of Modern Physics89, 015003 (2017)

    2017

  16. [24]

    Sundararaman, T

    R. Sundararaman, T. Christensen, Y . Ping, N. Rivera, J. D. Joannopoulos, M. Solja ˇci´c, and P. Narang, Plasmonics in ar- gentene, Physical Review Materials4, 074011 (2020)

    2020

  17. [25]

    Jablan and D

    M. Jablan and D. E. Chang, Multiplasmon absorption in graphene, Physical review letters114, 236801 (2015)

    2015

  18. [26]

    Sundararaman, K

    R. Sundararaman, K. Letchworth-Weaver, K. A. Schwarz, D. Gunceler, Y . Ozhabes, and T. A. Arias, Jdftx: Software for 6 joint density-functional theory, SoftwareX6, 278 (2017)

    2017

  19. [27]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Physical review letters77, 3865 (1996)

    1996

  20. [28]

    Schlipf and F

    M. Schlipf and F. Gygi, Optimization algorithm for the genera- tion of oncv pseudopotentials, Computer Physics Communica- tions196, 36 (2015)

    2015

  21. [29]

    Sundararaman and T

    R. Sundararaman and T. Arias, Regularization of the coulomb singularity in exact exchange by wigner-seitz truncated interac- tions: Towards chemical accuracy in nontrivial systems, Phys- ical Review B—Condensed Matter and Materials Physics87, 165122 (2013)

    2013

  22. [30]

    Grimme, J

    S. Grimme, J. Antony, S. Ehrlich, and H. Krieg, A consistent and accurate ab initio parametrization of density functional dis- persion correction (dft-d) for the 94 elements h-pu, The Journal of chemical physics132(2010)

    2010

  23. [31]

    Marzari and D

    N. Marzari and D. Vanderbilt, Maximally localized generalized wannier functions for composite energy bands, Physical review B56, 12847 (1997)

    1997

  24. [32]

    Giustino, M

    F. Giustino, M. L. Cohen, and S. G. Louie, Electron- phonon interaction using wannier functions, Physical Re- view B—Condensed Matter and Materials Physics76, 165108 (2007)

    2007

  25. [33]

    JJDFTX.jl,https://github.com/AliGhorashiCMT/ JJDFTX.jl, accessed 26 March, 2025

    2025

Pith tools

Reviewed August 15, 2026 · model on record in the stance chip above.