Pith. sign in

REVIEW 3 major objections 4 minor 44 references

A Recursive Hybrid Tetrahedron Method for Brillouin-zone Integration

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

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

Pith's one-line read A recursive tetrahedron refinement scheme turns Brillouin-zone integrals into weighted sums on the original k-grid.

desk verdict A genuinely useful recursive weight-collection extension of the hybrid tetrahedron method, with the main caveat being the honestly flagged band-sorting problem. read the letter →

arxiv 2411.17162 v1 pith:O6VRCIMX submitted 2024-11-26 cond-mat.mtrl-sci

Kun Dong , Yihao Lin , Xiaoqiang Liu , Jiechao Feng , Ji Feng This is my paper · ORCID
classification cond-mat.mtrl-sci PACS 71.15.-m
keywords Brillouin-zoneintegrationtetrahedronmethodquadraticinterpolationweighted-sumlinearresponsefunctionsLindhardfunctionKohn-Shamsusceptibilitysingularintegrands
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 recursive extension of the hybrid tetrahedron method that computes Brillouin-zone integrals as weighted sums over the original k-grid. The weights are generated by repeatedly dividing each quadratic tetrahedron into eight smaller tetrahedra using quadratic interpolation, then collecting linear-tetrahedron integration weights from the finest level back to the initial grid. The payoff is that refinement depth can be increased without storing the refined grid, and integrands carrying step, Dirac-delta, and pole singularities—the kind that appear in response and spectral functions—can be integrated in one framework. The authors demonstrate the convergence on the Lindhard function, the transverse spin susceptibility of the honeycomb-lattice Hubbard model, and the Kohn-Sham susceptibility of fcc cobalt.

What carries the argument

The load-bearing object is the quadratic tetrahedron: a tetrahedron with values at its four vertices and six edge midpoints, which determines a unique quadratic interpolant. Dividing it into eight equal-volume subordinate quadratic tetrahedra and adding the new edge-midpoint values as weighted linear combinations of the parent values makes the refinement repeatable. Repeating this division and then applying analytic linear-tetrahedron integration on the finest grid, with all weights collected back through the linear maps $Q^{(p,p-1)}$, yields the main weighted-sum identity Eq. (6).

What would settle it

Take a two-band model with a band crossing, compute a susceptibility with this method on a coarse grid using the code's sorted eigenvalues, and compare against a dense-grid reference computed with continuously labeled bands; if increasing the number of refinements does not shrink the error toward the dense-grid reference, the assumption of interpolable, continuous bands has failed.

Watch

Extended reading notes

Core claim

On its own terms, the paper claims that the main identity, Eq. (6), holds: after any number of tetrahedron refinements, a Brillouin-zone integral can be expressed as $\sum_m w_m F(k_m^{(0)})$, with the weights $w_m$ assembled by quadratic-interpolation maps from the finest grid back to the initial grid. The same recursive weight collection works when the singular weight function $W(k)$ contains a step function, a Dirac delta, a pole $1/D(k)$, or several of these at once, because the singular factors are evaluated on quadratically interpolated band eigenvalues and denominators and then treated by analytic linear-tetrahedron rules on the finest level. As a consequence, response and spectral functions, whose integrands contain denominators that can vanish, can be computed in the zero-broadening, zero-temperature limit rather than with a smearing parameter.

Load-bearing premise

The method assumes that band eigenvalues and denominator functions are smooth enough on the Brillouin zone that quadratic interpolation is faithful; the authors note that in practice eigenvalues from electronic-structure codes are sorted by size and become discontinuous at band crossings, and an effective band sorting algorithm is still needed.

Editorial extensions

If this is right

  • With enough refinements, the method's accuracy approaches that of direct quadratic tetrahedron integration while using only the initial coarse k-grid values.
  • The same weighted-sum machinery handles step, delta, and pole singularities, so response functions can be computed without artificial broadening or smearing.
  • The integral weights depend only on the band structure and can be computed once and reused, for example inside a self-consistent density-functional perturbation theory loop.
  • In the demonstrated cases, refinements remove artificial peaks and jaggedness produced by broadening methods and reduce the overestimated magnon energies and Goldstone gap on coarse grids.

Reading between the lines

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

  • Because Eq. (6) only requires a linear relation between refined and initial grid values, the recursive weight-collection idea could in principle be applied with interpolation schemes other than the quadratic tetrahedron, such as adaptive or higher-order local interpolants.
  • The authors' flagged band-sorting problem suggests a natural paired extension: if eigenvalues are relabeled to be continuous across band crossings, the method should extend cleanly to degenerate multi-band response functions.
  • The appendix's two-dimensional triangle formulas make the same recursive treatment available for surface and 2D-material calculations, although the paper only demonstrates the 3D case.
  • Reusing the stored weights could make self-consistent linear-response loops much cheaper at equal accuracy, since the expensive interpolation and analytic-integration steps would not be repeated on every iteration.
Share X Bluesky LinkedIn Reddit HN

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 paper proposes a recursive extension of the hybrid tetrahedron method for Brillouin-zone integration. The key result, Eq. (6), expresses a generic integral ∫ W(k)F(k) dk as a weighted sum of F on the initial k-grid, where the weights are obtained by propagating linear tetrahedron weights through the composite quadratic-interpolation map Q^(N,0). The weight function W(k) is assumed to depend on k through smooth functions (band eigenvalues ε(k), denominator D(k)) that can be quadratically interpolated; the method then handles step, delta, and 1/D singularities, including combinations such as Θ(ε_F−ε(k))/D(k). The Appendix provides closed-form linear tetrahedron weights for W=1/D with all limiting cases, as well as the standard step-function weights. Numerical demonstrations include the Lindhard function (compared to the exact result), the RPA transverse susceptibility of the honeycomb-lattice Hubbard model, and the Kohn-Sham susceptibility of fcc Co within TDDFT. A Julia implementation is released.

Significance. The central mathematical rearrangement in Eq. (6) is clean and correct, and the recursive weight-collection idea is a genuine practical improvement: it avoids storing exponentially many refined-grid values and makes the integral weights reusable for multiple integrands. The Lindhard-function test is a strong, falsifiable benchmark: the method converges systematically toward the exact analytic result with refinement, and the reported mean absolute error decreases rapidly with the number of refinements. The Appendix's closed-form weights for 1/D with all limiting cases is a useful reference for implementers, and the release of BZIntegral.jl is a concrete reproducibility asset. However, the paper's own Section V identifies an unresolved band-crossing problem that limits the method's practical reach to systems with smooth, continuously labeled bands; this must be addressed or clearly scoped in the manuscript before the broad claims in the abstract can be considered established.

major comments (3)
  1. [Section V and Eqs. (15), (18)] Section V explicitly concedes that "improper band labeling can cause uncertain discontinuities in band eigenvalue εnk" and that "an effective band sorting algorithm is needed for the method to work at its full potential." This is load-bearing for the advertised response-function applications, because Eqs. (15) and (18) mix band indices n,n′ and involve eigenvector-dependent matrix elements that are generically discontinuous at band crossings, while Table I and Eq. (6) rely on quadratic interpolation of smooth functions. The fcc Co benchmark in Fig. 5(c) is converged against the method's own 30×30×30, nr=1 result, so it does not independently validate accuracy in the presence of crossings. I request that the authors either supply a band-sorting algorithm, benchmark a crossing system against an independent exact or very-dense-grid reference, or explicitly restrict the claim of practical applicability to cases with continuous band labels.
  2. [Section III B, Eq. (9)] The treatment of weight functions containing products of step functions, Eq. (8) and Eq. (10), applies identity (9) and then invokes the linear tetrahedron rule on Θ(−x1x2), with the justification that this is "acceptable" when the final tetrahedra are sufficiently tiny. The manuscript provides no error estimate and no numerical test for this case, despite the abstract's claim of "simultaneously handling multiple singularities." I ask for a test of at least one product-of-step weight (e.g., a joint density of states or Eq. (10b)) or a clear statement that this part of the advertised functionality is heuristic and not benchmarked.
  3. [Appendix VII C, final paragraph] The final paragraph states that the 1/D weights "still hold after changing the arguments of logarithms from absolute values of D(k) to complex numbers D(k) itself," but the derivation preceding it treats real D with possible sign changes and absolute values. Since the response-function applications, Eqs. (13), (15), and (18), all use complex denominators with +iη, the analytic continuation of the closed-form weights to complex D should be justified, including the handling of branch cuts when D is near zero. Without this, a central ingredient of the flagship applications rests on an unproven assertion.
minor comments (4)
  1. [Section IV] Typographical errors should be corrected: "respnse" (Section IV, first paragraph), "coloser" (Section IV B), "mehtod" (Section IV B), and "studys" (Section IV B).
  2. [Fig. 3(c,d) caption] The caption's "log-log linear" phrase should be clarified: the mean absolute error appears to scale as a power of Δk^3 for fixed nr, and the exponential decrease with nr is a separate claim; please specify the axes and the functional form being reported.
  3. [Eq. (12)] Eq. (12) introduces Q^(N,0)_ij with subscript order that is consistent with Eq. (6), but the index convention is not explicitly stated when the composite map is first defined in Eq. (4); please define the direction of the map (from refined to initial) at its introduction.
  4. [Appendix VII C] The limiting-case notation "a&b → c" and "a&b&c → d" is not self-explanatory; please define these as simultaneous limits (e.g., a→c and b→c) so that the formulas are unambiguous.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: Eq. (6) is a linear-map bookkeeping identity, the method fits no parameters, and the self-citations are ancillary rather than load-bearing.

full rationale

Walking the derivation chain, the central result Eq. (6) is an algebraic identity rather than a fitted prediction: the composite quadratic-interpolation map Q^(N,0) expresses refined-grid values as linear combinations of initial-grid values, and the linear tetrahedron weights are functionals of W; substituting the former into the latter yields a weighted sum on the initial grid. No parameter is fitted to any benchmark, and the singular weight W is not interpolated directly but evaluated through interpolants of the smooth functions epsilon and D, an approximation whose validity is explicitly conditional on smoothness. The Lindhard test is compared against the exact analytical expression, providing an external check. The Hubbard-model and fcc-Co demonstrations use the method's own converged high-resolution result as a reference, which limits their independent confirmatory force but is a convergence study, not circularity. The only self-references are the companion Julia implementation [35] and the authors' earlier DFPT implementation [30], used for standard first-order perturbation equations and code infrastructure; these are not load-bearing for the derivation of Eq. (6). Section V explicitly concedes that 'an effective band sorting algorithm that restores the continuity in εnk is needed,' which is an honest robustness limitation rather than a circular step. Overall, the derivation is self-contained and no claimed prediction reduces by construction to its inputs.

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

The method introduces no new physical entities or fitted parameters. It relies on standard interpolation mathematics, a domain assumption about smooth band structure, and one paper-specific approximation for composite step-function weights. The main open assumption is the band-sorting continuity requirement.

assumptions (4)
  • standard math A unique quadratic form is determined by values at 10 points in R3 (tetrahedron vertices and edge midpoints).
    Invoked in Section II with footnote [11] to justify quadratic interpolation within each tetrahedron.
  • domain assumption The BZ is a parallelepiped tiled by a regular grid with open boundary conditions; crystal symmetry can recover the full BZ.
    Assumed in Section III A and discussed in Section V; the method does not assume periodic boundary conditions on the working grid.
  • domain assumption Band eigenvalues ε_n(k) and denominators D(k) are smooth enough for quadratic interpolation after proper band sorting.
    Stated in Section V as an 'unsettled issue': sorted band indices are discontinuous at band crossings, which can cause overfitting. The method relies on local smoothness of the interpolated functions.
  • ad hoc to paper For composite weight functions with products of step functions, applying the identity Θ(x1)Θ(x2)=1/2[Θ(x1)+Θ(x2)-Θ(-x1x2)] and linear tetrahedron rules is accurate for sufficiently fine tetrahedra.
    Introduced in Section III B after Eq. (9); the paper notes this is acceptable only when tetrahedra are sufficiently small, making it an approximation specific to this work.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A Recursive Hybrid Tetrahedron Method for Brillouin-zone Integration." pith.science (2026). https://pith.science/paper/O6VRCIMX

@misc{pith2026241117162,
  author       = {Pith},
  title        = {Pith review of: A Recursive Hybrid Tetrahedron Method for Brillouin-zone Integration},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/O6VRCIMX}},
  note         = {Machine review of arXiv:2411.17162}
}
read the original abstract

A recursive extension of the hybrid tetrahedron method for Brillouin-zone integration is proposed, allowing iterative tetrahedron refinement and significantly reducing the error from the linear tetrahedron method. The Brillouin-zone integral is expressed as a weighted sum on the initial grid, with integral weights collected recursively from the finest grid. Our method is capable of simultaneously handling multiple singularities in the integrand and thus may provide practical solutions to various Brillouin-zone integral tasks encountered in realistic calculations, including the computation of response and spectral function with superior sampling convergence. We demonstrate its effectiveness through numerical calculations of the density response functions of two model Hamiltonians and one real material system, the face-centered cubic cobalt.

Figures

Figures reproduced from arXiv: 2411.17162 by the authors.

Figure 1
Figure 1. FIG. 1. Overview of the recursive hybrid tetrahedron method. [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. (a) Quadratic tetrahedron: a tetrahedron with func [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. (a),(b)Real and imaginary part of [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: FIG. 4. (a), (b)Im [PITH_FULL_IMAGE:figures/full_fig_p008_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5. Im [PITH_FULL_IMAGE:figures/full_fig_p010_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. How to split tetrahedron. (a) A linear tetrahedron. Its vertices are labeled according to order of band eigenvalues, and [PITH_FULL_IMAGE:figures/full_fig_p011_6.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

44 extracted references · 33 canonical work pages

  1. [1]

    a → b, b ̸= c or d ϕ(a, b, c, d) → b log |b| 2(b − c)(b − d)

    work page

  2. [2]

    b → c, c ̸= a or d ϕ(a, b, c, d) + ϕ(a, c, b, d) = 1 36(a − c)3(c − d)2 h (a − c) 4a3 + 10ac(c − d) + c2(6c − d) − a2(8c + d) − 6a2 2a2 − 6ac + 6c2 + ad − 3cd log |a| + 6c2(2ac − 3ad + cd) log|c| i

    work page

  3. [3]

    a&b → c, a ̸= d ϕ(a, b, c, d) + ϕ(a, c, b, d) = − −c + d + (2c + d) log|c| 6(c − d)2

    work page

  4. [4]

    b&c → d, a ̸= d ϕ(a, b, c, d) + ϕ(a, c, d, b) + ϕ(a, d, b, c) = 2a3 + 3a2d − 6ad2 + d3 − 6a2d log |a| + 6a2d log |d| 12(a − d)4 14

    work page

  5. [5]

    The calculation of integral weights is susceptible to precision loss close to any of these limiting cases

    a&b&c → d ϕ(a, b, c, d) + ϕ(a, c, d, b) + ϕ(a, d, b, c) = 1 24d Case 1 can come simultaneously with case 2, for example, when a → d and b → c which means both ϕ(a, d, b, c) and ϕ(a, b, c, d) +ϕ(a, c, b, d) have limits to take care of. The calculation of integral weights is susceptible to precision loss close to any of these limiting cases. If any pair of ...

    work page

  6. [6]

    P. E. Bl¨ ochl, O. Jepsen, and O. K. Andersen, Im- proved tetrahedron method for Brillouin-zone integra- tions, Phys. Rev. B 49, 16223 (1994)

    1994

  7. [7]

    a → b, b ̸= c ϕ(a, b, c) = 1 4b − 4c

    work page

  8. [8]

    [9] for readers’ convenience

    This table is adapted from Ref. [9] for readers’ convenience. Index Value Index Value 11 6(7)+3(1)-(4) 24 4(5)+4(10)+2(7)-(1)-(4) 12 6(7)+3(4)-(1) 25 4(7)+4(10)+2(5)-(1)-(2) 13 6(9)+3(4)-(3) 26 4(10)+4(9)+2(8)-(3)-(2) 14 6(9)+3(3)-(4) 27 4(10)+4(8)+2(9)-(4)-(3) 15 6(8)+3(3)-(2) 28 4(9)+4(8)+2(10)-(2)-(4) 16 6(8)+3(2)-(3) 29 4(5)+4(8)+2(6)-(1)-(3) 17 6(5)+...

    work page

Show all 44 references
  1. [9]

    b → c, a ̸= c ϕ(a, b, c) + ϕ(a, c, b) = a2 − c2 − 2ac log |a| + 2ac log |c| 2(a − c)2

  2. [10]

    If D(k) is complex, the above formulas still hold after changing the arguments of logarithms from absolute values of D(k) to complex numbers D(k) itself

    a&b → c ϕ(a, b, c) + ϕ(a, c, b) = 1 6c As in the 3D case, we should take care of the situation when any pair of denominators are close to each other. If D(k) is complex, the above formulas still hold after changing the arguments of logarithms from absolute values of D(k) to co...

  3. [11]

    Methfessel and A

    M. Methfessel and A. T. Paxton, High-precision sampling for Brillouin-zone integration in metals, Phys. Rev. B 40, 3616 (1989). 15

    1989

  4. [12]

    Jepson and O.K

    O. Jepson and O.K. Anderson, The electronic structure of h.c.p. Ytterbium, Solid State Communications 9, 1763 (1971)

    1971

  5. [13]

    Lehmann and M

    G. Lehmann and M. Taut, On the Numerical Calculation of the Density of States and Related Properties, physica status solidi (b) 54, 469 (1972)

    1972

  6. [14]

    Wiesenekker, G

    G. Wiesenekker, G. te Velde, and E. J. Baerends, Ana- lytic quadratic integration over the two-dimensional Bril- louin zone, Journal of Physics C: Solid State Physics 21, 4263 (1988)

    1988

  7. [15]

    Wiesenekker and E

    G. Wiesenekker and E. J. Baerends, Quadratic integra- tion over the three-dimensional Brillouin zone, Journal of Physics: Condensed Matter 3, 6721 (1991)

    1991

  8. [16]

    Lambin and J

    Ph. Lambin and J. P. Vigneron, Computation of crystal Green’s functions in the complex-energy plane with the use of the analytical tetrahedron method, Phys. Rev. B 29, 3430 (1984)

    1984

  9. [17]

    Kawamura, Y

    M. Kawamura, Y. Gohda, and S. Tsuneyuki, Improved tetrahedron method for the Brillouin-zone integration ap- plicable to response functions, Phys. Rev. B 89, 094515 (2014)

    2014

  10. [18]

    A. H. MacDonald, S. H. Vosko, and P. T. Coleridge, Ex- tensions of the tetrahedron method for evaluating spec- tral properties of solids, Journal of Physics C: Solid State Physics 12, 2991 (1979)

    1979

  11. [19]

    Kohn and L

    W. Kohn and L. J. Sham, Self-Consistent Equations In- cluding Exchange and Correlation Effects, Phys. Rev. 140, A1133 (1965)

    1965

  12. [20]

    In general, a unique quadratic form can be determined given values at 10 points in R3, if no more than 3 points are colinear and no more than 6 points are coplanar

  13. [21]

    Kuneˇ s and R

    J. Kuneˇ s and R. Laskowski, Magnetic ground state and Fermi surface of bcc Eu, Phys. Rev. B 70, 174415 (2004)

    2004

  14. [22]

    M. D. Johannes and I. I. Mazin, Fermi surface nesting and the origin of charge density waves in metals, Phys. Rev. B 77, 165135 (2008)

    2008

  15. [23]

    R. S. Markiewicz, Bridging k and q space in the cuprates: Comparing angle-resolved photoemission and STM re- sults, Phys. Rev. B 69, 214517 (2004)

    2004

  16. [24]

    McElroy, G.-H

    K. McElroy, G.-H. Gweon, S. Y. Zhou, J. Graf, S. Uchida, H. Eisaki, H. Takagi, T. Sasagawa, D.-H. Lee, and A. Lanzara, Elastic Scattering Susceptibility of the High Temperature Superconductor Bi 2Sr2CaCu2O8+δ: A Comparison between Real and Momentum Space Pho- toemission Spectr...

    2006

  17. [25]

    Simon, F

    L. Simon, F. Vonau, and D. Aubel, A phenomenological approach of joint density of states for the determination of band structure in the case of a semi-metal studied by FT-STS, Journal of Physics: Condensed Matter 19, 355009 (2007)

    2007

  18. [26]

    Kubo, Statistical-Mechanical Theory of Irreversible Processes

    R. Kubo, Statistical-Mechanical Theory of Irreversible Processes. I. General Theory and Simple Applications to Magnetic and Conduction Problems, Journal of the Physical Society of Japan 12, 570 (1957)

    1957

  19. [27]

    R. B. Muniz and D. L. Mills, Theory of spin excitations in Fe(110) monolayers, Phys. Rev. B 66, 174417 (2002)

    2002

  20. [28]

    Karlsson and F

    K. Karlsson and F. Aryasetiawan, Spin-wave excitation spectra of nickel and iron, Phys. Rev. B 62, 3006 (2000)

    2000

  21. [29]

    Lindhard, On the properties of a gas of charged par- ticles, Kgl

    J. Lindhard, On the properties of a gas of charged par- ticles, Kgl. Danske Videnskab. Selskab Mat.-Fys. Medd. 28 (1954)

    1954

  22. [30]

    Kresse and J

    G. Kresse and J. Furthm¨ uller, Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set, Phys. Rev. B 54, 11169 (1996)

    1996

  23. [31]

    Giuliani and G

    G. Giuliani and G. Vignale, Quantum theory of the elec- tron liquid (Cambridge university press, 2005)

    2005

  24. [32]

    Hubbard, Electron correlations in narrow energy bands, Proceedings of the Royal Society of London

    J. Hubbard, Electron correlations in narrow energy bands, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences 276, 238 (1963)

    1963

  25. [33]

    Kanamori, Electron Correlation and Ferromagnetism of Transition Metals, Progress of Theoretical Physics 30, 275 (1963)

    J. Kanamori, Electron Correlation and Ferromagnetism of Transition Metals, Progress of Theoretical Physics 30, 275 (1963)

    1963

  26. [34]

    M. C. Gutzwiller, Effect of Correlation on the Ferromag- netism of Transition Metals, Phys. Rev. Lett. 10, 159 (1963)

    1963

  27. [35]

    Tasaki, The Hubbard model - an introduction and selected rigorous results, Journal of Physics: Condensed Matter 10, 4353 (1998)

    H. Tasaki, The Hubbard model - an introduction and selected rigorous results, Journal of Physics: Condensed Matter 10, 4353 (1998)

    1998

  28. [36]

    H. Tang, M. Plihal, and D. Mills, Theory of the spin dynamics of bulk Fe and ultrathin Fe(100) films, Journal of Magnetism and Magnetic Materials 187, 23 (1998)

    1998

  29. [37]

    Raczkowski, R

    M. Raczkowski, R. Peters, T. T. Ph` ung, N. Takemori, F. F. Assaad, A. Honecker, and J. Vahedi, Hubbard model on the honeycomb lattice: From static and dynam- ical mean-field theories to lattice quantum Monte Carlo simulations, Phys. Rev. B 101, 125103 (2020)

    2020

  30. [38]

    Sorella and E

    S. Sorella and E. Tosatti, Semi-Metal-Insulator Transi- tion of the Hubbard Model in the Honeycomb Lattice, Europhysics Letters 19, 699 (1992)

    1992

  31. [39]

    X. Liu, Y. Lin, and J. Feng, Implementation of the den- sity functional perturbation theory for generalized sus- ceptibility in the projector augmented wave framework, Phys. Rev. B 108, 094405 (2023)

    2023

  32. [40]

    Runge and E

    E. Runge and E. K. U. Gross, Density-Functional Theory for Time-Dependent Systems, Phys. Rev. Lett. 52, 997 (1984)

    1984

  33. [41]

    K. Cao, H. Lambert, P. G. Radaelli, and F. Giustino, Ab initio calculation of spin fluctuation spectra using time- dependent density functional perturbation theory, plane waves, and pseudopotentials, Phys. Rev. B 97, 024420 (2018)

    2018

  34. [42]

    R. W. Wyckoff, Cubic closest packed, ccp, structure, Crystal structures 1, 7 (1963)

    1963

  35. [43]

    Landau, On the vibrations of the electronic plasma, Journal of Experimental and Theoretical Physics 16, 574 (1946)

    L. Landau, On the vibrations of the electronic plasma, Journal of Experimental and Theoretical Physics 16, 574 (1946)

    1946

  36. [44]

    Yihao Lin, BZIntegral.jl, https://github.com/ SelimLin/BZIntegral.jl (2024), [Online; accessed November-2024]

    2024

Pith tools

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