Pith. sign in

REVIEW 2 major objections 7 minor 76 references

Large-scale pseudopotential density functional theory calculations using orthogonalized enriched finite element basis

T0 review · 2 major / 7 minor · reviewed 2026-07-10 · glm-5.2

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

Pith's one-line read Atom-centered enrichment cuts DFT cost 5–9×

desk verdict EFE basis delivers real 5-9x speedups for pseudopotential DFT, but the key approximation enabling it is only validated empirically, not analytically. read the letter →

arxiv 2607.06848 v1 pith:7F3GOVC4 submitted 2026-07-07 cond-mat.mtrl-sci

Avirup Sircar , Bikash Kanungo , Sambit Das , Vikram Gavini This is my paper · ORCID
classification cond-mat.mtrl-sci PACS 71.15.Ap71.15.Dx31.15.Ew
keywords basiscalculationselementfinitecomputationalfunctionspseudopotentialreduction
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

The paper proposes that augmenting a standard finite-element basis with compact, atom-centered functions — derived from single-atom Kohn–Sham solutions and then orthogonalized against the finite-element basis — can capture the sharp oscillatory behavior of electronic fields near nuclei without requiring a very fine mesh. This enriched finite element (EFE) basis retains the systematic convergence and boundary-condition flexibility of finite elements while approaching the per-atom efficiency of planewaves. The authors demonstrate that, for pseudopotential DFT calculations on non-periodic systems up to 39,083 electrons, the EFE basis requires 5–7× fewer degrees of freedom, achieves 5–9× lower wall-clock cost, and uses 4–5× less memory than the classical finite-element basis, all while maintaining chemical accuracy (~10⁻⁴ Ha/atom). The framework combines three technical ingredients: orthogonalization of enrichment functions (which block-diagonalizes the overlap matrix), an atom-block-diagonal approximation to the overlap-matrix inverse (which avoids cubic-scaling costs), and a residual-based Chebyshev subspace iteration that tolerates the resulting inexactness in the matrix inverse.

What carries the argument

The argument is carried by three coupled devices: (1) enrichment functions — truncated single-atom Kohn–Sham orbitals orthogonalized against the classical finite-element basis, which block-diagonalize the overlap matrix into classical-classical and enriched-enriched blocks; (2) an atom-block-diagonal approximation to the enriched-enriched overlap inverse, replacing a dense cubic-scaling inversion with O(N_a) independent small inversions; (3) a residual-based reformulation of Chebyshev filtered subspace iteration that remains convergent despite the inexact overlap inverse, because residuals are smaller than eigenvectors and thus more tolerant of matrix-vector product error.

What would settle it

A system where enrichment functions from neighboring atoms overlap significantly — e.g., a close-packed metal with very diffuse semi-core states or a strongly bonded system with short interatomic distances — such that the atom-block-diagonal approximation to the overlap inverse introduces errors the residual Chebyshev filtering cannot absorb, causing either eigenpair residual stagnation or loss of chemical accuracy.

Watch

Extended reading notes

Core claim

The central finding is that compact, orthogonalized atom-centered enrichment functions can be added to a finite-element basis for pseudopotential DFT without destroying conditioning or parallel scalability, and that the resulting basis closes most of the efficiency gap between finite elements and planewaves while preserving systematic convergence. The atom-block-diagonal approximation to the enriched-enriched overlap inverse, combined with residual Chebyshev filtering, is the mechanism that makes this practical: it sidesteps the cubic-scaling exact inverse and still converges to near machine-precision eigenpair residuals.

Load-bearing premise

The atom-block-diagonal approximation to the enriched-enriched overlap inverse neglects overlap between enrichment functions centered on different atoms. The paper shows this works for the tested systems — copper, sodium, platinum, and DNA — but its accuracy for materials with very diffuse valence orbitals or extremely close-packed atoms, where inter-atom enrichment overlap is large, remains the structural dependency on which the entire efficiency argument rests.

Editorial extensions

If this is right

  • If the EFE basis extends to periodic and semi-periodic boundary conditions as the authors intend, it could make finite-element DFT competitive with planewaves across the full range of system types, including metals and insulators, at both small and large scales.
  • The 4–5× memory reduction could push the accessible system-size frontier for real-space DFT from ~10⁵ electrons toward ~5×10⁵ electrons on the same hardware, enabling first-principles treatment of larger dislocation aggregates, quasicrystals, or biomolecular complexes.
  • The residual Chebyshev filtering approach, which tolerates inexact overlap inverses, is a general technique that could be transferred to other enriched or mixed bases where exact overlap inversion is the computational bottleneck.
  • If configurational forces can be computed efficiently in the EFE basis, the method could accelerate ab initio molecular dynamics substantially, since each time step inherits the per-SCF-iteration speedup demonstrated here.
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

2 major / 7 minor

Summary. This manuscript presents an enriched finite element (EFE) basis for pseudopotential Kohn-Sham DFT calculations, augmenting the classical finite element (CFE) basis with compact atom-centered enrichment functions derived from single-atom Kohn-Sham solutions. The enrichment functions are orthogonalized against the CFE basis, yielding a block-diagonal overlap matrix whose inverse is approximated efficiently via an atom-block-diagonal approximation. A residual-based Chebyshev filtering eigensolver is employed to tolerate the resulting inexactness in the overlap matrix inverse. The authors demonstrate systematic convergence to chemical accuracy (~10^-4 Ha/atom) on aspirin and Pt-13, and report 5-7x DoF reduction, 5-9x speedup, and 4-5x memory reduction over CFE (DFT-FE) on benchmark systems up to 39,083 electrons, with good parallel scalability.

Significance. The work addresses a well-known limitation of finite-element DFT—the large number of basis functions per atom compared to plane-waves—by a principled enrichment strategy that preserves systematic convergence. The orthogonalization trick yielding a block-diagonal overlap matrix, combined with the atom-block-diagonal inverse approximation and residual Chebyshev filtering, is a technically coherent and potentially impactful combination. The benchmark systems (Cu and Na nanoparticles, DNA) are non-trivial, and the demonstrated speedups over the highly optimized DFT-FE code are substantial. The framework is falsifiable: convergence tables allow independent verification of accuracy claims, and the performance metrics are clearly reported. The work is a meaningful contribution to large-scale DFT methodology.

major comments (2)
  1. §III.A, Eqs. (19)–(21); §III.B, Fig. 1: The atom-block-diagonal approximation to (M_ee)^{-1} is the key enabler for large-scale EFE calculations, yet its only direct convergence test (Fig. 1) uses a single Cu atom. For a single atom, the atom-block-diagonal approximation is exact (there is only one atom block), so Fig. 1 validates only the GLL quadrature approximation to (M_cc)^{-1}, not the inter-atom enrichment overlap neglect that is the actual concern for multi-atom systems. For close-packed Cu systems with nearest-neighbor distance 4.8 bohr and enrichment cutoff radii up to 10 a.u. (§II.B), enrichment functions from neighboring atoms overlap substantially, making the neglect of off-diagonal atom blocks non-trivial. While the benchmark calculations (Tables III–V) empirically demonstrate convergence to 10^-4 Ha/atom, there is no direct comparison of convergence rates or eigen-residual
  2. §IV.B, Tables III–V: The performance comparisons between EFE and CFE use different discretization parameters (e.g., EFE: h_min=1.25, p=5 vs. CFE: h_min=0.8, p=6 for Cu). Both are stated to achieve ~10^-4 Ha/atom accuracy, but the reported errors are not matched: e.g., Table III Cu 2-shell shows EFE error 1.1e-4 vs. CFE error 8e-5. It would strengthen the claims to either (i) match the error levels more closely by tuning discretization, or (ii) explicitly discuss the sensitivity of the speedup ratios to the target accuracy, since the 5-9x speedup figure depends on both methods being at comparable accuracy.
minor comments (7)
  1. §II.B, Eq. (12): The enrichment function construction involves several free parameters (cutoff smoothness t=0.33, r0 selection criteria, unity integration tolerance 5e-3, radial derivative threshold 5e-3). A brief discussion of how sensitive the results are to these choices, or at minimum a statement that these were determined from numerical experiments on small systems and found to be robust, would help the reader.
  2. Tables I–II: For the Pt-13 system (Table II), the planewave basis achieves similar accuracy with ~3x fewer DoFs than EFE, while for aspirin (Table I) they are comparable. A brief comment on why EFE is less competitive with planewaves for Pt would provide useful context.
  3. Table III, Cu 8-shell: Only per-SCF iteration time is reported (no full ground-state calculation). The text explains this is due to CFE cost, but reporting the number of SCF iterations used for the smaller systems would help the reader extrapolate the Cu 8-shell results to a sense of total wall-time.
  4. §IV.C, Figs. 2–3: The y-axis label 'Relative Speedup' is shown but the ideal scaling line is not drawn. Adding an ideal linear scaling reference line would make the parallel efficiency assessment clearer. Also, the x-axis uses different node ranges for the two figures (2–32 vs. 4–64); using relative speedup with respect to the minimum node count would make the two plots more directly comparable.
  5. §II.B: The enrichment functions are truncated at a maximum r0 of 10 a.u. For systems with very diffuse states (e.g., Rydberg-like states, or weakly bound anions), this cutoff may be insufficient. A brief remark on the applicability limits would be appropriate.
  6. The manuscript is limited to non-periodic systems. While the authors note that extensions to periodic and semi-periodic systems are planned (§V), a brief remark on whether any technical obstacles are anticipated for periodic boundary conditions (e.g., enrichment function overlap across cell boundaries) would be welcome.
  7. §III.B, Eq. (24): The notation and derivation of the residual Chebyshev filtering closely follows Ref. [65]. While the key ideas are summarized, a reader unfamiliar with that reference may find the jump from Eq. (23) to Eq. (24) abrupt. A sentence connecting the recurrence relation to the residual formulation more explicitly would aid readability.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity found; derivation chain is self-contained against external benchmarks

full rationale

The paper's derivation chain is straightforward and non-circular. (1) The EFE basis is constructed by augmenting CFE with atom-centered enrichment functions derived from single-atom Kohn-Sham solutions (Eq. 12) — these are independent inputs, not outputs being predicted. (2) The orthogonalization procedure (Eq. 12) is cited from Ref. [51] (same-group prior work), but it is a mathematical construction whose properties (block-diagonal overlap matrix, Eq. 16) follow directly from the orthogonalization, not from a fitted or assumed result. (3) The atom-block-diagonal approximation to (Mee)^{-1} (Eq. 19-21) is a numerical approximation whose accuracy is empirically validated through benchmark calculations (Tables III-V) against external reference energies from DFT-FE and Quantum Espresso — not a self-referential validation. (4) The residual Chebyshev filtering method is cited from Ref. [65] and applied as an external algorithm; its tolerance to inexact inverses is demonstrated empirically (Fig. 1), not claimed by construction. (5) All performance claims (5-7x DoF reduction, 5-9x speedup, 4-5x memory reduction) are measured against the CFE basis in DFT-FE and planewave basis in QE, both external codes with independently verified accuracy. The self-citations to Refs. [50, 51, 58] provide methodological foundations (orthogonalization, prior all-electron EFE demonstrations) but do not form a load-bearing chain where the present result is defined in terms of its own outputs. The skeptic's concern about Fig. 1 using a single atom (where the atom-block-diagonal approximation is trivially exact) is a correctness/robustness concern, not a circularity issue — the paper does not claim Fig. 1 validates the multi-atom approximation, and the multi-atom accuracy is independently verified through benchmark energy comparisons.

Assumptions & free parameters 4 free parameters · 6 assumptions · 2 invented entities

The paper introduces several hand-tuned parameters for enrichment function construction (t, r0 selection criteria) that are determined empirically from small-system experiments. The core mathematical framework relies on standard DFT and pseudopotential theory. The key ad-hoc assumption is the atom-block-diagonal approximation, which is validated empirically rather than formally proven. No new physical entities or forces are postulated.

free parameters (4)
  • cutoff smoothness parameter t = 0.33
    Fixed value chosen from numerical experiments on small systems (Sec. II.B, Eq. 11). Controls the smoothness of the enrichment function truncation.
  • enrichment function cutoff radius r0 = varies per enrichment function, max 10 a.u.
    Selected per enrichment function based on three criteria: capturing last turning point, integral unity within 5e-3 tolerance, and tail-region location (Sec. II.B).
  • unity integration tolerance = 5e-3
    Tolerance for enrichment function integration to unity, used in r0 selection (Sec. II.B).
  • radial derivative threshold = 5e-3
    Threshold for absolute value of radial derivative at r0 to ensure cutoff is in tail region (Sec. II.B).
assumptions (6)
  • domain assumption Kohn-Sham DFT formulation with local density approximation (Perdew-Wang 92)
    Standard DFT framework adopted throughout (Sec. II.A). Choice of LDA is stated as not affecting main conclusions.
  • domain assumption Norm-conserving ONCV pseudopotentials from PseudoDojo
    All calculations use ONCV PSPs (Sec. IV). Results may not generalize to ultrasoft or PAW pseudopotentials.
  • ad hoc to paper Single-atom Kohn-Sham eigenfunctions are suitable enrichment functions
    The enrichment functions are constructed from isolated atom solutions (Sec. II.B). Authors note this is a convenient choice and do not claim optimality.
  • ad hoc to paper Atom-block-diagonal approximation to Mee inverse is sufficiently accurate when combined with residual Chebyshev filtering
    The approximation in Eq. 19-21 neglects inter-atom enrichment overlap. Justified by Fig. 1 for Cu atom and benchmark results, but not formally proven.
  • standard math Spectral finite elements with GLL quadrature yield diagonal Mcc
    Established in Ref. 36, used to simplify overlap inverse evaluation (Sec. III.A).
  • domain assumption Residual Chebyshev filtering is tolerant to inexact matrix-vector products
    From Ref. 65, invoked in Sec. III.B to justify using approximate overlap inverse in eigensolver.
invented entities (2)
  • Enrichment functions (NE_k,I) independent evidence
    purpose: Compact atom-centered functions augmenting the CFE basis to capture nuclear-region electronic structure
    Derived from single-atom Kohn-Sham solutions with smooth truncation. Their effect is verifiable through systematic convergence benchmarks (Tables I-V) and comparison to CFE and plane-wave results.
  • Atom-block-diagonal overlap inverse approximation (M_tilde_ee)^-1 independent evidence
    purpose: Efficient approximation to the enriched-enriched overlap matrix inverse
    The approximation is validated through benchmark calculations showing chemical accuracy is maintained (Tables III-V) and through the Cu atom residual convergence demonstration (Fig. 1).

how reviews work

0 comments
Cite this review

Pith. "Pith review of Large-scale pseudopotential density functional theory calculations using orthogonalized enriched finite element basis." pith.science (2026). https://pith.science/paper/7F3GOVC4

@misc{pith2026260706848,
  author       = {Pith},
  title        = {Pith review of: Large-scale pseudopotential density functional theory calculations using orthogonalized enriched finite element basis},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/7F3GOVC4}},
  note         = {Machine review of arXiv:2607.06848}
}
abstract

We present an efficient and scalable computational framework for pseudopotential Kohn-Sham density functional theory (KS-DFT) calculations using an enriched finite element (EFE) basis. The EFE basis is formed by augmenting the classical finite element (CFE) basis with compact atom-centered functions, which we term enrichment functions. The key idea is to combine the completeness of a finite element basis with the efficiency of an atom-centered basis. We orthogonalize the enrichment functions with respect to the underlying CFE basis to simultaneously improve the conditioning of the EFE basis and the efficiency of evaluating the inverse of the overlap matrix. To efficiently solve the Kohn-Sham eigenvalue problem, we employ a residual-based Chebyshev subspace iteration approach that is tolerant to approximations in the evaluation of the inverse of the overlap matrix. We demonstrate the accuracy of the framework as compared to the widely available DFT packages. For benchmark non-periodic calculations, ranging up to 39,083 electrons, the EFE basis offers a $5-7\times$ reduction in degrees of freedom over the CFE basis. As a result of this, EFE achieves a $5-9\times$ reduction in computational cost over the CFE basis. The EFE basis also provides a $4-5\times$ reduction in the required memory compared to the CFE basis, thus allowing for optimal utilization of computational resources. Finally, we demonstrate that the EFE basis affords good parallel scalability. Overall, the EFE basis offers a systematically convergent, fast, scalable, resource-efficient basis for pseudopotential DFT calculations.

Figures

Figures reproduced from arXiv: 2607.06848 by the authors.

Figure 1
Figure 1. FIG. 1. Comparison of standard and residual Chebyshev fil [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Breakdown of total wall-time into the various com [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Breakdown of total wall-time into the various compu [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

76 extracted references · 76 canonical work pages

  1. [1]

    results in the following discrete eigenvalue problem HEψ i =ǫiMEψ i, (14) where HE and ME are the discrete Kohn-Sham Hamil- tonian matrix and overlap matrix in the EFE basis; ǫi denotes the i-th discrete Kohn-Sham eigenvalue; and ψ i = { {ψ C i,j }j=1,...,n h; {ψ E i,k,I }k=1,...,n I ;I=1,...,N a } denotes the corresponding eigenvector containing the expa...

    work page

  2. [2]

    (22) The above approximation of ψ i in terms of Vout, includ- ing the projected eigenvalue problem, is often referred to as the Rayleigh-Ritz procedure

    to lie in the space spanned by Vout, that is, ψ i ≈ Voutdi, where di is given by the solution to the following projected eigenvalue problem V† outHEVoutdi =ǫiV† outMEVoutdi. (22) The above approximation of ψ i in terms of Vout, includ- ing the projected eigenvalue problem, is often referred to as the Rayleigh-Ritz procedure. Typically, for a fixed A, we pe...

    work page

  3. [3]

    We still use ME in the Rayleigh-Ritz procedure in Eq

    as an inexpensive approximation to ( ME)− 1 only for the purpose of Chebyshev filter- ing. We still use ME in the Rayleigh-Ritz procedure in Eq. 22. While this offers efficiency, the use of ( ˜ME)− 1 can affect the robustness and accuracy of the Kohn-Sham eigensolve. To illustrate this, we consider the multi-pass Chebyshev filtering for a copper (Cu) atom. As e...

    work page

  4. [4]

    Subsequently, the Chebyshev subspace can be evaluated as Vout = ( ˜ME)− 1 ˜Zm + VinΛ

    as well as the fact that ˜ME(ME)− 1 ≈ I, Zk can be approximated as Zk ≈ ˜Zk = 2σk e HE(˜ME)− 1 ˜Zk− 1 − 2cσk e ˜Zk− 1 − σk− 1σk− 2 ˜Zk− 2 + 2σk e GΛ , (25) where G = HEVin − MEVinΛ is the eigen-residual of the approximate vectors and eigenvalues; ˜Z0 = 0; and ˜Z1 = σ 1 e G . Subsequently, the Chebyshev subspace can be evaluated as Vout = ( ˜ME)− 1 ˜Zm + V...

    work page

  5. [5]

    Hohenberg and W

    P. Hohenberg and W. Kohn, Inhomogeneous electron gas, Physical review 136, B864 (1964)

    work page 1964

  6. [6]

    Kohn and L

    W. Kohn and L. J. Sham, Self-consistent equations in- cluding exchange and correlation effects, Physical review 140, A1133 (1965)

    work page 1965

  7. [7]

    A. J. Cohen, P. Mori-S´ anchez, and W. Yang, Challenges for density functional theory, Chemical reviews 112, 289 (2012)

    work page 2012

  8. [8]

    Burke, Perspective on density functional theory, The Journal of chemical physics 136 (2012)

    K. Burke, Perspective on density functional theory, The Journal of chemical physics 136 (2012)

    work page 2012

Show all 76 references
  1. [9]

    R. O. Jones, Density functional theory: Its origins, ris e to prominence, and future, Reviews of modern physics 87, 897 (2015)

    2015

  2. [10]

    Mardirossian and M

    N. Mardirossian and M. Head-Gordon, Thirty years of density functional theory in computational chemistry: an overview and extensive assessment of 200 density func- tionals, Molecular physics 115, 2315 (2017)

    2017

  3. [11]

    A. M. Teale, T. Helgaker, A. Savin, C. Adamo, B. Aradi, A. V. Arbuznikov, P. W. Ayers, E. J. Baerends, V. Barone, P. Calaminici, et al. , DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science, Physical chemistry chemical physics 24, 287...

    2022

  4. [12]

    W. E. Pickett, Pseudopotential methods in condensed matter applications, Computer Physics Reports 9, 115 (1989)

    1989

  5. [13]

    J. R. Chelikowsky, The pseudopotential-density func- tional method applied to nanostructures, Journal of Physics D: Applied Physics 33, R33 (2000)

    2000

  6. [14]

    Schwerdtfeger, The pseudopotential approximation in electronic structure theory, ChemPhysChem 12, 3143 (2011)

    P. Schwerdtfeger, The pseudopotential approximation in electronic structure theory, ChemPhysChem 12, 3143 (2011)

    2011

  7. [15]

    R. M. Martin, Electronic structure: basic theory and practical methods (Cambridge university press, 2020)

    2020

  8. [16]

    Hamann, M

    D. Hamann, M. Schl¨ uter, and C. Chiang, Norm- conserving pseudopotentials, Phys. Rev. Lett. 43, 1494 (1979)

    1979

  9. [17]

    Vanderbilt, Soft self-consistent pseudopotential s in a generalized eigenvalue formalism, Physical review B 41, 7892 (1990)

    D. Vanderbilt, Soft self-consistent pseudopotential s in a generalized eigenvalue formalism, Physical review B 41, 7892 (1990)

    1990

  10. [18]

    P. E. Bl¨ ochl, Projector augmented-wave method, Physi - cal review B 50, 17953 (1994)

    1994

  11. [19]

    Goedecker, M

    S. Goedecker, M. Teter, and J. Hutter, Separable dual- space gaussian pseudopotentials, Physical Review B 54, 1703 (1996)

    1996

  12. [20]

    D. R. Hamann, Optimized norm-conserving vanderbilt pseudopotentials, Physical Review B—Condensed Mat- ter and Materials Physics 88, 085117 (2013)

    2013

  13. [21]

    S. Carr, S. Fang, and E. Kaxiras, Electronic-structure methods for twisted moir´ e layers, Nature Reviews Mate- rials 5, 748 (2020)

    2020

  14. [22]

    Ismail-Beigi and T

    S. Ismail-Beigi and T. Arias, Ab initio study of screw dislocations in Mo and Ta: a new picture of plasticity in bcc transition metals, Physical Review Letters 84, 1499 (2000)

    2000

  15. [23]

    Shin and E

    I. Shin and E. A. Carter, Possible origin of the discrep- ancy in Peierls stresses of fcc metals: First-principles simulations of dislocation mobility in aluminum, Physi- cal Review B—Condensed Matter and Materials Physics 88, 064106 (2013)

    2013

  16. [24]

    Das and V

    S. Das and V. Gavini, Electronic structure study of scre w dislocation core energetics in aluminum and core energet- ics informed forces in a dislocation aggregate, Journal of the Mechanics and Physics of Solids 104, 115 (2017)

    2017

  17. [25]

    Rodney, L

    D. Rodney, L. Ventelon, E. Clouet, L. Pizzagalli, and F. Willaime, Ab initio modeling of dislocation core prop- erties in metals and semiconductors, Acta Materialia 124, 633 (2017)

    2017

  18. [26]

    Das and V

    S. Das and V. Gavini, Intrinsic ductility enhancement i n Mg alloys elucidated via large-scale ab-initio calculatio ns, Acta Materialia 315, 122377 (2026)

    2026

  19. [27]

    A. Dive, C. Benmore, M. Wilding, S. Martin, S. Beck- man, and S. Banerjee, Molecular dynamics modeling of the structure and Na+-ion transport in Na2S+ SiS2 glassy electrolytes, The Journal of Physical Chemistry B 122, 7597 (2018)

    2018

  20. [28]

    D. J. Cole and N. D. Hine, Applications of large-scale density functional theory in biology, Journal of Physics: Condensed Matter 28, 393001 (2016)

    2016

  21. [29]

    M. C. Payne, M. P. Teter, D. C. Allan, T. A. Arias, and J. D. Joannopoulos, Iterative minimiza- tion techniques for ab initio total-energy calcula- tions: molecular dynamics and conjugate gradients, Rev. Mod. Phys. 64, 1045 (1992)

    1992

  22. [30]

    Kresse and J

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

    1996

  23. [31]

    Giannozzi, S

    P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococ- cioni, I. Dabo, et al. , QUANTUM ESPRESSO: a modu- lar and open-source software project for quantum simula- tions of materials, Journal of physics: Condensed matter 2...

    2009

  24. [32]

    Gonze, J.-M

    X. Gonze, J.-M. Beuken, R. Caracas, F. Detraux, M. Fuchs, G.-M. Rignanese, L. Sindic, M. Verstraete, G. Zerah, F. Jollet, et al. , First-principles computation of material properties: the ABINIT software project, Com- putational Materials Science 25, 478 (2002)

    2002

  25. [33]

    Jensen, Atomic orbital basis sets, Wiley Interdisci - plinary Reviews: Computational Molecular Science 3, 273 (2013)

    F. Jensen, Atomic orbital basis sets, Wiley Interdisci - plinary Reviews: Computational Molecular Science 3, 273 (2013)

    2013

  26. [34]

    T. D. K¨ uhne, M. Iannuzzi, M. Del Ben, V. V. Ry- bkin, P. Seewald, F. Stein, T. Laino, R. Z. Khaliullin, 13 O. Sch¨ utt, F. Schiffmann, et al. , CP2K: An electronic structure and molecular dynamics software package- quickstep: Efficient and accurate electronic structure cal- cul...

    2020

  27. [35]

    J. M. Soler, E. Artacho, J. D. Gale, A. Garc ´ ıa, J. Jun- quera, P. Ordej´ on, and D. S´ anchez-Portal, The SIESTA method for ab initio order-N materials simulation, Jour- nal of physics: Condensed matter 14, 2745 (2002)

    2002

  28. [36]

    Bowler, R

    D. Bowler, R. Choudhury, M. Gillan, and T. Miyazaki, Recent progress with large-scale ab initio calculations: the conquest code, Physica Status Solidi (b) 243, 989 (2006)

    2006

  29. [37]

    S. R. White, J. W. Wilkins, and M. P. Teter, Finite- element method for electronic structure, Physical Review B 39, 5819 (1989)

    1989

  30. [38]

    Tsuchida and M

    E. Tsuchida and M. Tsukada, Electronic-structure cal- culations based on the finite-element method, Physical Review B 52, 5573 (1995)

    1995

  31. [39]

    J. E. Pask, B. M. Klein, P. A. Sterne, and C. Y. Fong, Finite-element methods in electronic-structure theory, Computer Physics Communications 135, 1 (2001)

    2001

  32. [40]

    Motamarri, M

    P. Motamarri, M. R. Nowak, K. Leiter, J. Knap, and V. Gavini, Higher-order adaptive finite-element meth- ods for Kohn–Sham density functional theory, Journal of Computational Physics 253, 308 (2013)

    2013

  33. [41]

    J. R. Chelikowsky, N. Troullier, and Y. Saad, Finite- difference-pseudopotential method: Electronic structure calculations without a basis, Physical review letters 72, 1240 (1994)

    1994

  34. [42]

    Kronik, A

    L. Kronik, A. Makmal, M. L. Tiago, M. Alemany, M. Jain, X. Huang, Y. Saad, and J. R. Chelikowsky, PARSEC–the pseudopotential algorithm for real-space electronic structure calculations: recent advances and novel applications to nano-structures, Physica Status So- lidi (b) 243, ...

    2006

  35. [43]

    Ghosh and P

    S. Ghosh and P. Suryanarayana, SPARC: Accurate and efficient finite-difference formulation and parallel imple- mentation of density functional theory: Isolated clusters , Computer Physics Communications 212, 189 (2017)

    2017

  36. [44]

    S. Das, P. Motamarri, V. Subramanian, D. M. Rogers, and V. Gavini, DFT-FE 1.0: A massively parallel hybrid CPU-GPU density functional theory code using finite- element discretization, Computer Physics Communica- tions 280, 108473 (2022)

    2022

  37. [45]

    S. Das, P. Motamarri, V. Gavini, B. Turcksin, Y. W. Li, and B. Leback, Fast, scalable and accurate finite-element based ab initio calculations using mixed precision com- puting: 46 PFLOPS simulation of a metallic dislocation system, in Proceedings of the international conferenc...

    2019

  38. [46]

    Motamarri, S

    P. Motamarri, S. Das, S. Rudraraju, K. Ghosh, D. Davy- dov, and V. Gavini, DFT-FE: A massively parallel adap- tive finite-element code for large-scale density functiona l theory calculations, Computer Physics Communications 246, 106853 (2020)

    2020

  39. [47]

    Ghosh, H

    K. Ghosh, H. Ma, V. Gavini, and G. Galli, All-electron density functional calculations for electron and nuclear spin interactions in molecules and solids, Physical Review Materials 3, 043801 (2019)

    2019

  40. [48]

    Ghosh, H

    K. Ghosh, H. Ma, M. Onizhuk, V. Gavini, and G. Galli, Spin–spin interactions in defects in solids from mixed all- electron and pseudopotential first-principles calculatio ns, npj Comput Mater 7, 123 (2021)

    2021

  41. [49]

    S. Das, B. Kanungo, V. Subramanian, G. Panigrahi, P. Motamarri, D. Rogers, P. Zimmerman, and V. Gavini, Large-scale materials modeling at quantum accuracy: Ab initio simulations of quasicrystals and interacting ex- tended defects in metallic alloys, in Proceedings of the Inter...

    2023

  42. [50]

    Zhuravel, H

    R. Zhuravel, H. Huang, G. Polycarpou, S. Polydorides, P. Motamarri, L. Katrivas, D. Rotem, J. Sperling, L. A. Zotti, A. B. Kotlyar, et al. , Backbone charge transport in double-stranded DNA, Nature nanotechnology 15, 836 (2020)

    2020

  43. [51]

    W. Baek, S. Das, S. Tan, V. Gavini, and W. Sun, Qua- sicrystal stability and nucleation kinetics from density functional theory, Nature Physics , 1 (2025)

    2025

  44. [52]

    Subramanian, S

    V. Subramanian, S. Das, and V. Gavini, Tucker tensor approach for accelerating fock exchange computations in a real-space finite-element discretization of generalized Kohn–Sham density functional theory, Journal of Chem- ical Theory and Computation 20, 3566 (2024)

    2024

  45. [53]

    Kanungo, P

    B. Kanungo, P. M. Zimmerman, and V. Gavini, Exact exchange-correlation potentials from ground-state elec- tron densities, Nature communications 10, 4497 (2019)

    2019

  46. [54]

    Kanungo and V

    B. Kanungo and V. Gavini, Large-scale all-electron den - sity functional theory calculations using an enriched finite-element basis, Physical Review B 95, 035112 (2017)

    2017

  47. [55]

    N. D. Rufus, B. Kanungo, and V. Gavini, Fast and ro- bust all-electron density functional theory calculations in solids using orthogonalized enriched finite elements, Physical Review B 104, 085112 (2021)

    2021

  48. [56]

    J. P. Perdew and Y. Wang, Accurate and simple ana- lytic representation of the electron-gas correlation ener gy, Physical review B 45, 13244 (1992)

    1992

  49. [57]

    J. Pask, N. Sukumar, and S. Mousavi, Linear scaling solu - tion of the all-electron coulomb problem insolids, Interna - tional Journal for Multiscale Computational Engineering 10 (2012)

    2012

  50. [58]

    J. E. Pask, B. M. Klein, C. Y. Fong, and P. A. Sterne, Real-space local polynomial basis for solid-state electronic-structure calculations: A finite-element ap- proach, Phys. Rev. B 59, 12352 (1999)

    1999

  51. [59]

    J. E. Pask and P. A. Sterne, Finite element meth- ods in ab initio electronic structure calculations, Modelling and Simulation in Materials Science and Engineer ing 13, R71

  52. [60]

    Suryanarayana, V

    P. Suryanarayana, V. Gavini, T. Blesgen, K. Bhat- tacharya, and M. Ortiz, Non-periodic finite-element formulation of Kohn–Sham density functional theory, Journal of the Mechanics and Physics of Solids 58, 256 (2010)

    2010

  53. [61]

    Motamarri, M

    P. Motamarri, M. Iyer, J. Knap, and V. Gavini, Higher-order adaptive finite-element meth- ods for orbital-free density functional theory, Journal of Computational Physics 231, 6596 (2012)

    2012

  54. [62]

    Kanungo, N

    B. Kanungo, N. D. Rufus, and V. Gavini, Efficient all- electron time-dependent density functional theory calcu- lations using an enriched finite element basis, Journal of Chemical Theory and Computation 19, 978 (2023)

    2023

  55. [63]

    J. M. Melenk and I. Babuˇ ska, The partition of unity finite element method: Basic theory and applications, Computer Methods in Applied Mechanics and Engineering 139, 289

  56. [64]

    Babuˇ ska and J

    I. Babuˇ ska and J. M. Melenk, The partition of unity method, International Journal for Numerical Methods in Engineerin g 40, 727 14

  57. [65]

    J. E. Pask and N. Sukumar, Partition of unity finite ele- ment method for quantum mechanical materials calcula- tions, Extreme Mechanics Letters 11, 8 (2017)

    2017

  58. [66]

    Y. Zhou, Y. Saad, M. L. Tiago, and J. R. Chelikowsky, Self-consistent-field calculations using chebyshev-filte red subspace iteration, Journal of Computational Physics 219, 172 (2006)

    2006

  59. [67]

    Y. Zhou, J. R. Chelikowsky, and Y. Saad, Chebyshev- filtered subspace iteration method free of sparse diago- nalization for solving the Kohn–Sham equation, Journal of Computational Physics 274, 770 (2014)

    2014

  60. [68]

    Y. Zhou, Y. Saad, M. L. Tiago, and J. R. Che- likowsky, Parallel self-consistent-field calculations vi a chebyshev-filtered subspace acceleration, Physical Re- view E—Statistical, Nonlinear, and Soft Matter Physics 74, 066704 (2006)

    2006

  61. [69]

    Kodali, K

    N. Kodali, K. Ramakrishnan, and P. Motamarri, Residual-based chebyshev filtered subspace iteration for sparse hermitian eigenvalue problems tolerant to inexact matrix-vector products, Computer Physics Communica- tions 327, 110239 (2026)

    2026

  62. [70]

    M. J. Van Setten, M. Giantomassi, E. Bousquet, M. J. Verstraete, D. R. Hamann, X. Gonze, and G.-M. Rig- nanese, The PseudoDojo: Training and grading a 85 ele- ment optimized norm-conserving pseudopotential table, Computer Physics Communications 226, 39 (2018)

    2018

  63. [71]

    H. F. Walker and P. Ni, Anderson acceleration for fixed- point iterations, SIAM Journal on Numerical Analysis 49, 1715 (2011)

    2011

  64. [72]

    Giannozzi, O

    P. Giannozzi, O. Andreussi, T. Brumme, O. Bunau, M. B. Nardelli, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, M. Cococcioni, et al. , Advanced capabili- ties for materials modelling with Quantum ESPRESSO, Journal of physics: Condensed matter 29, 465901 (2017)

    2017

  65. [73]

    A. Deng, H. Li, M. Bo, Z. Huang, L. Li, C. Yao, and F. Li, Understanding atomic bonding and electronic dis- tributions of a DNA molecule using DFT calculation and BOLS-BC model, Biochemistry and Biophysics Reports 24, 100804 (2020)

    2020

  66. [74]

    Motamarri and V

    P. Motamarri and V. Gavini, Configurational forces in electronic structure calculations using Kohn-Sham den- sity functional theory, Phys. Rev. B 97, 165132 (2018)

    2018

  67. [75]

    N. D. Rufus and V. Gavini, Ionic forces and stress tensor in all-electron density functional the- ory calculations using an enriched finite-element basis, Phys. Rev. B 106, 085108 (2022)

    2022

  68. [76]

    Kanungo and V

    B. Kanungo and V. Gavini, Real time time-dependent density functional theory using higher order finite- element methods, Phys. Rev. B 100, 115148 (2019)

    2019

Pith tools

Reviewed July 10, 2026 · model on record in the stance chip above.