Pith. sign in

REVIEW 2 major objections 5 minor 98 references

Unified open-boundary electrostatics in real-space density functional theory

T0 review · 2 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

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

Pith's one-line read One local variational functional with Green's-function Dirichlet values yields the exact electrostatic potential, energy, forces, and stress for isolated, 1D, and 2D periodic systems in DFT, with or without an applied uniform field.

desk verdict Unified open-boundary electrostatics in real-space DFT is a genuine, well-tested advance, but the field-dependent stress tensor (Eq. 24) is asserted without derivation and looks wrong on tensorial grounds. read the letter →

arxiv 2608.08474 v1 pith:ODG4LLFC submitted 2026-08-09 cond-mat.mtrl-sci

Rajat Kumar , David Codony , Phanish Suryanarayana , Abhiraj Sharma This is my paper · ORCID
classification cond-mat.mtrl-sci PACS 71.15.Mb
keywords real-spacedensityfunctionaltheoryopen-boundaryelectrostaticsDirichletboundaryconditionsPoissonequationuniformelectricfieldmultipoleexpansionstaticpolarizabilitypiezoelectriccoefficients
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 claims that the open-boundary electrostatics of isolated molecules, 1D-periodic wires, and 2D-periodic slabs can be handled in real-space density functional theory by a single local variational functional, Eq. (3), whose maximizer is the exact electrostatic potential of the total charge density plus any applied uniform electric field along the open directions. The required Dirichlet boundary values are written analytically from Green's functions: spherical multipoles for isolated systems, cylindrical multipoles with modified Bessel functions for wires, and a dipole step with in-plane Fourier terms for slabs. If the claim is right, the vacuum region only needs to be large enough for the electron density to decay, not large enough to suppress multipole-image errors, so real-space DFT would need substantially less vacuum than plane-wave supercell calculations at the same accuracy. The paper verifies exponential convergence with vacuum size, agreement with a plane-wave reference, stresses that match numerical derivatives of the energy, and polarizabilities and piezoelectric coefficients that agree with published values.

What carries the argument

The load-bearing object is the local variational functional of Eq. (3), a maximum over $\phi \in \mathcal{U} \subset H^1(\Omega)$, the affine space of functions periodic along the periodic directions and fixed to $\phi_0(\mathbf{x}) + \mathbf{x}\cdot\mathbf{E}$ on the open boundaries, with boundary surface integrals that absorb the slow or growing decay of the potential along the open directions. This functional replaces the nonlocal Coulomb kernel with a differential operator whose inversion reproduces the long-range interaction exactly, so stationarity of the functional guarantees consistency between the energy and its derivatives. The Dirichlet values $\phi_0$ come from the Green's functions of each geometry: spherical harmonics for the 0D case, cylindrical multipole moments together with zeroth-order modified Bessel functions $K_0$ for the 1D case, and a dipole step term with exponentially screened in-plane Fourier components for the 2D case, each valid where the charge density has vanished and convergent for a charge-neutral system.

What would settle it

Run the method on a neutral molecule at progressively smaller vacuum sizes until the electron density at the boundary is no longer negligible, and compare the energy and forces against a well-converged plane-wave calculation: the exponential-convergence curve should flatten or diverge at the point where the density support reaches the boundary. Equivalently, apply the formulation to a charged molecule, for which the zeroth multipole term is dropped: the computed energy will drift from the exact value as the cell grows, exposing the neutrality condition as the binding assumption.

Watch

Extended reading notes

Core claim

The central discovery is a maximization formulation of electrostatics that makes no assumption that the potential or its gradient vanishes on the open boundaries. The paper defines an affine space of functions $\phi$ that are periodic along the periodic directions and take prescribed values $\phi_0(\mathbf{x}) + \mathbf{x}\cdot\mathbf{E}$ on the open faces, where $\phi_0$ is the potential of the total charge density alone. Eq. (3) expresses the electrostatic energy as the maximum over this space of a local functional containing surface integrals on the open boundaries; the maximizer solves the Poisson equation $-\frac{1}{4\pi}\nabla^2\phi = \rho + b$ for the total charge density, and substituting it back yields the closed-form energy, Eq. (5). The Dirichlet data $\phi_0$ are derived by a Green's-function ansatz for each dimensionality: the spherical multipole series of Eq. (9) for isolated systems, the cylindrical multipole and Bessel-function series of Eq. (13) for 1D periodic systems, and the dipole step plus exponentially screened in-plane Fourier series of Eq. (17) for 2D periodic systems. The force and stress follow by differentiation, with the electric field entering only through $\phi$; the stress, Eq. (24), is the zero-field stress plus an isotropic diagonal term, and because no existing code provides stresses for such systems they are validated against numerical derivatives of the energy.

Load-bearing premise

The load-bearing premise is that the total charge density (electrons plus ionic pseudocharge) is exactly zero on and beyond the open boundary faces and that the system is charge neutral, so the multipole and Fourier series for the Dirichlet values converge; if the vacuum is so small that the density has not decayed, an atom crosses the boundary, or the system carries a net charge, the boundary values are no longer exact and the claimed exactness collapses.

Editorial extensions

If this is right

  • One framework covers isolated, 1D-periodic, and 2D-periodic systems, with an applied uniform electric field entering entirely through the electrostatic potential rather than through an explicit change to the Kohn–Sham Hamiltonian.
  • Energy, forces, and stresses converge exponentially with vacuum size, reaching roughly $10^{-6}$ Ha/atom and $10^{-5}$ Ha/Bohr at about 10 Bohr of vacuum, several times less vacuum than the plane-wave reference calculations needed in the test cases.
  • The stress tensor is available for systems with a non-vanishing potential on the open boundaries, so cell relaxation and equation-of-state studies need no finite-difference energy derivatives for such systems.
  • Static polarizabilities and piezoelectric coefficients of polar low-dimensional systems follow from the same machinery without separate correction schemes, matching published values.
  • The additional cost of evaluating the Dirichlet boundary values is negligible relative to the Poisson solve, so the parallel scalability of the real-space solver is retained.
  • The boundary conditions also accommodate fixed-potential electrodes and could be matched to bulk boundary conditions for semi-infinite surface calculations.

Reading between the lines

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

  • Charged systems are the natural next stress test: the authors explicitly defer them to future work, and the omission of the $\ell = 0$ term in Eq. (9) suggests a neutralizing-background term would need to be added to $\phi_0$ before the formulation extends to net-charged species.
  • The built-in consistency check between the energy derivative and the dipole moment could be turned into an automated diagnostic for the truncation parameters, flagging when $\ell_{\max}$, $m_{\max}$, or $Q^{\max}_{mn}$ is too small for the system's multipole content.
  • Because the electric field enters only through $\phi$, the first-order variation of the boundary values should slot naturally into real-space density functional perturbation theory, a path the authors name as future work.
  • The Green's-function ansatz is the only part that changes with geometry, so the machinery could likely be re-derived for cyclic and helical symmetry, giving open-boundary electrostatics for bent and twisted nanostructures.
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

2 major / 5 minor

Summary. The manuscript develops a unified real-space electrostatics framework for Kohn–Sham density functional theory covering isolated, 1D-periodic, and 2D-periodic systems with open boundaries and a uniform applied electric field along the open directions. The central construction is the local variational functional of Eq. (3), whose stationarity is claimed to yield the Poisson equation for the total (electron plus pseudocharge) density, with Dirichlet values supplied by the analytical Green's-function/multipole expressions in Eqs. (9), (13), and (17). The paper further derives ground-state energy, atomic forces, and stress, implements them in SPARC, validates exponential vacuum convergence and agreement with Quantum ESPRESSO, and applies the framework to static polarizabilities and piezoelectric coefficients.

Significance. If the central claim holds, the framework provides a unified, parameter-free treatment of open-boundary electrostatics with applied fields, eliminating multipole-image errors and substantially reducing the vacuum required in real-space DFT calculations, while also supplying a stress tensor not currently available in other implementations. The strengths of the paper are its explicit variational derivation, the absence of fitted parameters in the boundary values, the exponential-convergence validation, the direct comparisons with an established plane-wave code, and the open data repository. These features make the work a potentially important contribution to real-space electronic structure methods. The main caveat is the stress expression in Eq. (24), which is asserted without derivation and appears inconsistent with the energy functional from which it is supposed to follow; because stress is part of the paper's central claim, this issue is load-bearing.

major comments (2)
  1. [Sec. 3, Eq. (24)] The field-dependent stress contribution is stated without derivation and its tensor structure is inconsistent with the energy functional from which it is supposed to follow. In Eqs. (5) and (22), the applied-field part of the electrostatic energy is the dipole coupling ∫_Ω (ρ+b) x·E. Under an affine strain ε, the coordinates transform as x → (1+ε)x, so the first-order change of this term is E_i μ_j ε_ij, with μ_j = ∫ (ρ+b) x_j. The corresponding stress contribution is the symmetric part of (1/|Ω|) E_i μ_j, which is anisotropic and carries no factor of 1/2. Equation (24) instead asserts an isotropic half-dipole term, δ_αβ (1/(2|Ω|)) ∫ (ρ+b) x·E. The numerical validation in Sec. 5.2 is too weak to distinguish the two: the errors are quoted as percentages of unspecified stress magnitudes, and the field-induced change is small relative to the zero-field stress. The authors should either derive Eq. (24) from Eq. (22), including all surface, boundary-value, and pseudocharge terms, or replace it; as written, the claimed exactness of the stress is unsupported.
  2. [Sec. 3, Eq. (23)] The force expression is obtained by differentiating the ground-state energy and invoking the Hellmann–Feynman theorem, but the admissible space of potentials in Eq. (3) has Dirichlet values φ0(R) + x·E that depend on atomic positions through φ0. For a variational problem with R-dependent constraints, the stationarity condition alone does not eliminate the derivative of the constraint. The manuscript does not show that the boundary terms involving ∂φ0/∂R_I vanish, and the asserted expression in Eq. (23) has no such terms. Please provide this step or state the conditions under which the contribution of the R-dependent Dirichlet values cancels. The agreement with Quantum ESPRESSO in Table 1 is reassuring, but it does not replace the missing derivation.
minor comments (5)
  1. [Sec. 2, after Eq. (6)] The phrase 'the the' should be corrected to 'the'.
  2. [Sec. 2.2, Eq. (13)] The symbols x2D and θ are used in Eq. (13) and defined only afterward; please move the definitions before the equation or restructure the sentence for readability.
  3. [Sec. 3, Eq. (22)] The quantities b̃ and V_c appear in Eq. (22) but are defined only in the following text; please define them before first use.
  4. [Sec. 5.2] The stress validation reports percentage errors without giving the reference stress magnitudes. Please report the actual zero-field stresses and the field-induced changes so that the reader can assess the significance of the 0.1% level agreement.
  5. [Fig. 4] The convergence plots show multiple truncation-parameter curves, but the caption and legends do not identify which curve corresponds to which truncation parameter. Please add this information for reproducibility.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity found; the boundary-value electrostatics is derived from the charge density and anchored to external benchmarks, with Eq. (24) as a derivation gap rather than a circular reduction.

full rationale

The central derivation is self-contained. Equation (2) is recast into the local variational problem Eq. (3), whose Euler–Lagrange equation is the Poisson equation Eq. (4); the affine Dirichlet space uses phi0, which is not fitted but computed from the instantaneous total charge density rho+b via Green's-function multipole expansions in Eqs. (8)–(17). The closed-form energy Eq. (5), the Hellmann–Feynman force expression Eq. (23), and the zero-field stress taken from Ref. [75] form a derivation chain, not a circular reduction: no target energy, force, polarizability, or piezoelectric coefficient is used to determine the Dirichlet values or the functional. External anchors are independent of the fitted values: Quantum ESPRESSO comparisons for energy, forces, and polarization density (Table 1), literature values for polarizability and piezoelectric coefficients, and finite-difference energy derivatives for stress. The only notable weakness is Eq. (24): the field-dependent stress is asserted as 'extending the zero-field result of Sharma et al. [75]' with no derivation, and its isotropic half-dipole form is in tension with the full dipole coupling in Eq. (5)/(22). This is an omitted proof and a correctness risk, not circularity: the stress is not made true by construction, nor is it a fitted parameter renamed as a prediction. Ref. [75] is a self-citation, but it supplies a parameter-free prior zero-field result, and the novel field-dependent term is not derived from it, so the self-citation does not make the central claim circular.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

No new physical entities are invented. The electrostatic boundary values are derived from standard Green's functions and multipole expansions; the only adjustable numbers are series truncation orders and numerical settings, all of which are convergence-tested rather than fitted to target data.

free parameters (3)
  • lmax = 6 in production, 9 in convergence reference
    Spherical multipole truncation for the 0D boundary values, Eq. (9); convergence is checked in Fig. 4a, so this is a numerical setting, not a fitted physical parameter.
  • mmax, nmax = (3,0) in production, (4,2) in convergence reference
    Cylindrical multipole and axial Fourier truncation for the 1D boundary values, Eq. (13); convergence is checked in Fig. 4b.
  • Qmax_mn = 0 in production, 0.11 Bohr^-1 in convergence reference
    In-plane Fourier truncation for the 2D boundary values, Eq. (17); convergence is checked in Fig. 4c.
assumptions (5)
  • domain assumption System is charge neutral: total charge integral of (rho + b) is zero.
    Invoked after Eq. (7) and in each Green's function derivation to drop the l=0 multipole and to make the boundary potential convergent; the paper states that extension to charged systems is future work.
  • domain assumption Total charge density has compact support strictly inside Omega and vanishes on and beyond the open boundaries dOmega_alpha.
    Stated in Sec. 2 opening; required for the multipole and Fourier expansions of phi0 on the boundary and for the claimed domain-reduction benefit.
  • domain assumption The applied uniform field is produced by a distant source electronically disconnected from the system, so the potential in Omega is x dot E.
    Sec. 2, first paragraph; this is the standard way to treat a uniform field in a finite cell.
  • domain assumption Pseudopotential approximation with local pseudocharges b_I; point nuclei are replaced by smooth pseudocharge densities.
    Used throughout Eqs. (1) through (2) and in Sec. 4; the central claim is within pseudopotential DFT.
  • standard math The Green's functions used for the screened Poisson equations are the standard fundamental solutions: log, K0, and exponential decay.
    Eqs. (8), (12), and (16); these are classical results, cited to Lennard-Jones and Dent [67].

how reviews work

0 comments
Cite this review

Pith. "Pith review of Unified open-boundary electrostatics in real-space density functional theory." pith.science (2026). https://pith.science/paper/ODG4LLFC

@misc{pith2026260808474,
  author       = {Pith},
  title        = {Pith review of: Unified open-boundary electrostatics in real-space density functional theory},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ODG4LLFC}},
  note         = {Machine review of arXiv:2608.08474}
}
read the original abstract

We present an electrostatic formulation in real-space density functional theory that provides a systematic and unified treatment of the open-boundary electrostatics of isolated and partially periodic systems, including in the presence of an applied uniform electric field along the open (finite) directions. Specifically, we formulate a local electrostatic energy functional whose stationarity yields the Poisson equation for the electrostatic potential, subject to periodic and Dirichlet boundary conditions along the periodic and open directions, respectively. Using a Green's function approach, we derive analytical expressions for the Dirichlet values arising from the total charge density of the system. We also derive the expressions for the energy, atomic forces, and stress tensor. We implement the resulting expressions within the large-scale parallel real-space SPARC electronic structure code. Using representative examples, we verify the accuracy and efficiency of the framework, demonstrating exponential convergence of the computed quantities with vacuum size and excellent agreement with established plane-wave codes while requiring significantly less vacuum at comparable accuracy. Since no existing implementation provides the stresses for such systems, we instead verify them against numerical derivatives of the energy, finding close agreement. Finally, we apply the framework to compute static polarizabilities and piezoelectric coefficients, obtaining very good agreement with values reported in the literature.

Figures

Figures reproduced from arXiv: 2608.08474 by the authors.

Figure 1
Figure 1. Illustration of a 2-atom unit cell Ω centered at the origin and defined by the vectors 𝐋1 , 𝐋2 , and 𝐋3 . The pair of boundary faces of Ω spanned by {𝐋𝛽 ∶ 𝛽 ≠ 𝛼} is denoted by 𝜕Ω𝛼 , 𝛼, 𝛽 ∈ {1, 2, 3}. Consider a 𝑑-dimensional (where 𝑑 ∈ {0, 1, 2, 3}) charge-neutral system spanned by the vectors {𝐋𝛼 ∶ 𝛼 ∈ 𝑝 ⊆ {1, 2, 3}, |𝑝 | = 𝑑} in the periodic directions and defined by the vectors {𝐋𝛼 ∶ 𝛼 ∈ 𝑜 = {1, 2, 3} ⧵ 𝑝 , |… view at source ↗
Figure 2
Figure 2. Schematic representation of the electrostatic potential 𝜙̄(𝑧), averaged over the plane perpendicular to the open 𝑧-direction, for isolated (red), 1D periodic (green), and 2D periodic (blue) systems. The black line represents the linear potential due to an applied uniform electric field, and the two vertical gray lines mark the boundary planes (𝜕Ω3 ) outside which the average electron density ̄𝜌(𝑧) (dashed curve) van… view at source ↗
Figure 3
Figure 3. Illustration showing the atomic arrangement and the unit cell with vacuum for (a) the 4,5-diaminophthalonitrile molecule, (b) the polycarbonitrile wire, and (c) the molybdenum sulfoselenide Janus monolayer. The top and bottom rows show projections onto the 𝑥–𝑦 and 𝑦–𝑧 planes, respectively. Here, 𝑥𝑣 , 𝑦𝑣 , and 𝑧𝑣 denote the vacuum size in the 𝑥, 𝑦, and 𝑧 directions, respectively. Atomic species are indicated in the l… view at source ↗
Figures from the paper (5 more)
Figure 4
Figure 4. Figure 4: Convergence of the energy (solid blue), forces (dashed red), and stresses (dotted green) with respect to vacuum size for (a) the 4,5-diaminophthalonitrile molecule, (b) the polycarbonitrile wire, and (c) the molybdenum sulfoselenide monolayer, at different truncation p…
Figure 5
Figure 5. Figure 5: a shows the change in computed energy (markers) with electric field together with its quadratic fit (curves), and Fig. 5b shows the change in computed dipole moment (markers) with electric field together with the derivative of the fitted energy with respect to the elec…
Figure 6
Figure 6. Figure 6: Induced polarization density as a function of the applied electric field for (a) the 4,5-diaminophthalonitrile molecule, (b) the polycarbonitrile wire, and (c) the molybdenum sulfoselenide monolayer. Markers are computed values and solid lines are linear fits; the slop…
Figure 7
Figure 7. Figure 7: Polarization density as a function of applied strain for (a) the polycarbonitrile wire and (b) the molybdenum sulfoselenide monolayer. Markers are computed values and solid lines are linear fits; the slopes give the piezoelectric coefficients 𝑒3𝑗𝑗, with 𝑅2 indicating t…
Figure 8
Figure 8. Figure 8: Convergence of the energy (solid blue), forces (dashed red), and stresses (dotted green) with respect to vacuum size for (a) a long polycarbonitrile wire, and (b) a large molybdenum sulfoselenide monolayer, both with randomly perturbed atomic positions. References [1] …

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

98 extracted references · 78 canonical work pages

  1. [1]

    Hohenberg, W

    P. Hohenberg, W. Kohn, Inhomogeneous electron gas, Phys. Rev. 136 (3B) (1964) B864

    1964

  2. [2]

    W. Kohn, L. J. Sham, Self-consistent equations including exchange and correlation effects, Phys. Rev. 140 (4A) (1965) A1133

    1965

  3. [3]

    Burke, Perspective on density functional theory, J

    K. Burke, Perspective on density functional theory, J. Chem. Phys. 136 (2012) 150901

    2012

  4. [4]

    A. D. Becke, Perspective: Fifty years of density-functional theory in chemical physics, J. Chem. Phys. 140 (2014) 18A301

    2014

  5. [5]

    W. E. Pickett, Pseudopotential methods in condensed matter applications, Comput. Phys. Rep. 9 (3) (1989) 115–197

    1989

  6. [6]

    Martin, Electronic Structure: Basic theory and practical methods, Cambridge University Press, 2004

    R. Martin, Electronic Structure: Basic theory and practical methods, Cambridge University Press, 2004

    2004

  7. [7]

    Kresse, J

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

    1996

  8. [8]

    (2005) 567–570

    S.J.Clark,M.D.Segall,C.J.Pickard,P.J.Hasnip,M.J.Probert,K.Refson,M.C.Payne,Firstprinciplesmethodsusingcastep,Z.Kristallogr. (2005) 567–570

    2005

Show all 98 references
  1. [9]

    X.Gonze,J.-M.Beuken,R.Caracas,F.Detraux,M.Fuchs,G.-M.Rignanese,L.Sindic,M.Verstraete,G.Zerah,F.Jollet,etal.,First-principles computation of material properties: the abinit software project, Comput. Mater. Sci. 25 (3) (2002) 478–492

    2002

  2. [10]

    Giannozzi, S

    P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococcioni, I. Dabo, et al., Quantum espresso:amodularandopen-sourcesoftwareprojectforquantumsimulationsofmaterials,J.Phys.Condens.Matter21(39)(2009)395502

    2009

  3. [11]

    Ismail-Beigi, T

    S. Ismail-Beigi, T. Arias, New algebraic formulation of density functional calculation, Comput. Phys. Commun. 128 (1-2) (2000) 1–45

    2000

  4. [12]

    Gygi, Architecture of qbox: A scalable first-principles molecular dynamics code, IBM J

    F. Gygi, Architecture of qbox: A scalable first-principles molecular dynamics code, IBM J. Res. Dev. 52 (1.2) (2008) 137–144. Rajat et al.: Preprint submitted to arXiv Page 18 of 21 Unified open-boundary electrostatics in real-space DFT

    2008

  5. [13]

    Valiev, E

    M. Valiev, E. J. Bylaska, N. Govind, K. Kowalski, T. P. Straatsma, H. J. J. Van Dam, D. Wang, J. Nieplocha, E. Aprà, T. L. Windus, et al., Nwchem:Acomprehensiveandscalableopen-sourcesolutionforlargescalemolecularsimulations,Comput.Phys.Commun.181(9)(2010) 1477–1489

    2010

  6. [14]

    P. P. Ewald, Die berechnung optischer und elektrostatischer gitterpotentiale, Ann. Phys. (Leipzig) 369 (3) (1921) 253–287

    1921

  7. [15]

    J. Ihm, A. Zunger, M. L. Cohen, Momentum-space formalism for the total energy of solids, J. Phys. C: Solid State Phys. 12 (21) (1979) 4409–4422

    1979

  8. [16]

    Makov, M

    G. Makov, M. C. Payne, Periodic boundary conditions in ab initio calculations, Phys. Rev. B. 51 (7) (1995) 4014

    1995

  9. [17]

    K. Kunc, R. Resta, External fields in the self-consistent theory of electronic states: a new method for direct evaluation of macroscopic and microscopic dielectric response, Phys. Rev. Lett. 51 (8) (1983) 686

    1983

  10. [18]

    Neugebauer, M

    J. Neugebauer, M. Scheffler, Adsorbate-substrate and adsorbate-adsorbate interactions of na and k adlayers on al (111), Phys. Rev. B 46 (24) (1992) 16067

    1992

  11. [19]

    Bengtsson, Dipole correction for surface supercell calculations, Phys

    L. Bengtsson, Dipole correction for surface supercell calculations, Phys. Rev. B. 59 (19) (1999) 12301

    1999

  12. [20]

    Meyer, D

    B. Meyer, D. Vanderbilt, Ab initio study of BaTiO3 and PbTiO3 surfaces in external electric fields, Phys. Rev. B 63 (20) (2001) 205426

    2001

  13. [21]

    C. A. Rozzi, D. Varsano, A. Marini, E. K. Gross, A. Rubio, Exact coulomb cutoff technique for supercell calculations, Phys. Rev. B. 73 (20) (2006) 205119

    2006

  14. [22]

    Ismail-Beigi, Truncation of periodic image interactions for confined systems, Phys

    S. Ismail-Beigi, Truncation of periodic image interactions for confined systems, Phys. Rev. B 73 (23) (2006) 233103

    2006

  15. [23]

    I.Dabo,B.Kozinsky,N.E.Singh-Miller,N.Marzari,Electrostaticsinperiodicboundaryconditionsandreal-spacecorrections,Phys.Rev.B 77 (11) (2008) 115139

    2008

  16. [24]

    G. J. Martyna, M. E. Tuckerman, A reciprocal space based method for treating long range interactions in ab initio and force-field-based calculations in clusters, J. Chem. Phys. 110 (6) (1999) 2810

    1999

  17. [25]

    R. N. Barnett, U. Landman, Born-oppenheimer molecular-dynamics simulations of finite systems: Structure and dynamics of (h 2 o) 2, Phys. Rev. B 48 (4) (1993) 2081

    1993

  18. [26]

    T.Sohier,M.Calandra,F.Mauri,Densityfunctionalperturbationtheoryforgatedtwo-dimensionalheterostructures:Theoreticaldevelopments and application to flexural phonons in graphene, Phys. Rev. B 96 (7) (2017) 075448

    2017

  19. [27]

    Rivano, N

    N. Rivano, N. Marzari, T. Sohier, Density functional perturbation theory for one-dimensional systems: Implementation and relevance for phonons and electron-phonon interactions, Phys. Rev. B 109 (24) (2024) 245426

    2024

  20. [28]

    C.Freysoldt,P.Eggert,P.Rinke,A.Schindlmayr,M.Scheffler,Screeningintwodimensions:Gwcalculationsforsurfacesandthinfilmsusing the repeated-slab approach, Phys. Rev. B: Condens. Matter Mater. Phys. 77 (23) (2008) 235428

    2008

  21. [29]

    A. D. Becke, Basis-set-free density-functional quantum chemistry, Int. J. Quantum Chem. 36 (S23) (1989) 599–609

    1989

  22. [30]

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

    1989

  23. [31]

    J. R. Chelikowsky, N. Troullier, Y. Saad, Finite-difference-pseudopotential method: Electronic structure calculations without a basis, Phys. Rev. Lett. 72 (8) (1994) 1240–1243

    1994

  24. [32]

    A.P.Seitsonen,M.J.Puska,R.M.Nieminen,Real-spaceelectronic-structurecalculations:Combinationofthefinite-differenceandconjugate- gradient methods, Phys. Rev. B 51 (20) (1995) 14057

    1995

  25. [33]

    Tsuchida, M

    E. Tsuchida, M. Tsukada, Electronic-structure calculations based on the finite-element method, Phys. Rev. B 52 (8) (1995) 5573

    1995

  26. [34]

    E.Briggs,D.Sullivan,J.Bernholc,Real-spacemultigrid-basedapproachtolarge-scaleelectronicstructurecalculations,Phys.Rev.B54(20) (1996) 14362

    1996

  27. [35]

    Fattebert, Finite difference schemes and block rayleigh quotient iteration for electronic structure calculations on composite grids, J

    J.-L. Fattebert, Finite difference schemes and block rayleigh quotient iteration for electronic structure calculations on composite grids, J. Comput. Phys. 149 (1) (1999) 75–94

    1999

  28. [36]

    T. A. Arias, Multiresolution analysis of electronic structure: semicardinal and wavelet bases, Rev. Mod. Phys. 71 (1) (1999) 267–311

    1999

  29. [37]

    Shimojo, R

    F. Shimojo, R. K. Kalia, A. Nakano, P. Vashishta, Linear-scaling density-functional-theory calculations of electronic structure based on real- space grids: design, analysis, and scalability test of parallel algorithms, Comput. Phys. Commun. 140 (3) (2001) 303–314

    2001

  30. [38]

    Skylaris, P

    C.-K. Skylaris, P. D. Haynes, A. A. Mostofi, M. C. Payne, Introducing onetep: Linear-scaling density functional simulations on parallel computers, J. Chem. Phys. 122 (8) (2005) 084119

    2005

  31. [39]

    J. E. Pask, P. A. Sterne, Finite element methods in ab initio electronic structure calculations, Model. Simul. Mater. Sci. Eng. 13 (2005) R71–R96

    2005

  32. [40]

    Bowler, R

    D. Bowler, R. Choudhury, M. Gillan, T. Miyazaki, Recent progress with large-scale ab initio calculations: the conquest code, Phys. Status Solidi B. 243 (5) (2006) 989–1000

    2006

  33. [41]

    Status Solidi B

    A.Castro,H.Appel,M.Oliveira,C.A.Rozzi,X.Andrade,F.Lorenzen,M.A.Marques,E.Gross,A.Rubio,Octopus:atoolfortheapplication of time-dependent density functional theory, Phys. Status Solidi B. 243 (11) (2006) 2465–2488

    2006

  34. [42]

    Genovese, A

    L. Genovese, A. Neelov, S. Goedecker, T. Deutsch, S. A. Ghasemi, A. Willand, D. Caliste, O. Zilberberg, M. Rayson, A. Bergman, et al., Daubechies wavelets as a basis set for density functional pseudopotential calculations, The J. Chem. Phys. 129 (1) (2008) 014109

    2008

  35. [43]

    Iwata, D

    J.-I. Iwata, D. Takahashi, A. Oshiyama, T. Boku, K. Shiraishi, S. Okada, K. Yabana, A massively-parallel electronic-structure calculations based on real-space density functional theory, J. Comput. Phys. 229 (6) (2010) 2339–2363

    2010

  36. [44]

    P.Suryanarayana,V.Gavini,T.Blesgen,K.Bhattacharya,M.Ortiz,Non-periodicfinite-elementformulationofkohn–shamdensityfunctional theory, J. Mech. Phys. Solids 58 (2) (2010) 256–280

    2010

  37. [45]

    P.Suryanarayana,K.Bhattacharya,M.Ortiz,Amesh-freeconvexapproximationschemeforkohn–shamdensityfunctionaltheory,J.Comput. Phys. 230 (13) (2011) 5226–5238

    2011

  38. [46]

    L. Lin, J. Lu, L. Ying, et al., Adaptive local basis set for kohn–sham density functional theory in a discontinuous galerkin framework I: Total energy calculation, J. Comput. Phys. 231 (4) (2012) 2140–2154

    2012

  39. [47]

    Ghosh, P

    S. Ghosh, P. Suryanarayana, Sparc: Accurate and efficient finite-difference formulation and parallel implementation of density functional theory: Isolated clusters, Comput. Phys. Commun. 212 (2017) 189–204. Rajat et al.: Preprint submitted to arXiv Page 19 of 21 Unified open-b...

    2017

  40. [48]

    Ghosh, P

    S. Ghosh, P. Suryanarayana, Sparc: Accurate and efficient finite-difference formulation and parallel implementation of density functional theory: Extended systems, Comput. Phys. Commun. 216 (2017) 109–125

    2017

  41. [49]

    Q.Xu,P.Suryanarayana,J.E.Pask,Discretediscontinuousbasisprojectionmethodforlarge-scaleelectronicstructurecalculations,J.Chem. Phys. 149 (9) (2018) 094104

    2018

  42. [50]

    P.Motamarri,S.Das,S.Rudraraju,K.Ghosh,D.Davydov,V.Gavini,Dft-fe–amassivelyparalleladaptivefinite-elementcodeforlarge-scale density functional theory calculations, Comput. Phys. Commun. 246 (2020) 106853

    2020

  43. [51]

    T. L. Beck, Real-space mesh techniques in density-functional theory, Rev. Mod. Phys. 72 (4) (2000) 1041–1080

    2000

  44. [52]

    Y. Saad, J. R. Chelikowsky, S. M. Shontz, Numerical methods for electronic structure calculations of materials, SIAM Rev. 52 (1) (2010) 3–54

    2010

  45. [53]

    P. P. Pratapa, P. Suryanarayana, J. E. Pask, Spectral quadrature method for accurate o (n) electronic structure calculations of metals and insulators, Comput. Phys. Commun. (2015)

    2015

  46. [54]

    P.Suryanarayana,P.P.Pratapa,A.Sharma,J.E.Pask,Sqdft:Spectralquadraturemethodforlarge-scaleparallelo(n)kohn–shamcalculations at high temperature, Comput. Phys. Commun. 224 (2018) 288–298

    2018

  47. [55]

    Gavini, S

    V. Gavini, S. Baroni, V. Blum, D. R. Bowler, A. Buccheri, J. R. Chelikowsky, S. Das, W. Dawson, P. Delugas, M. Dogan, et al., Roadmap on electronic structure codes in the exascale era, Model. Simul. Mater. Sci. Eng. 31 (6) (2023) 063301

    2023

  48. [56]

    A. S. Banerjee, P. Suryanarayana, Cyclic density functional theory: A route to the first principles simulation of bending in nanostructures, J. Mech. Phys. Solids 96 (2016) 605–631

    2016

  49. [57]

    Ghosh, A

    S. Ghosh, A. S. Banerjee, P. Suryanarayana, Symmetry-adapted real-space density functional theory for cylindrical geometries: Application to large group-iv nanotubes, Phys. Rev. B 100 (12) (2019) 125143

    2019

  50. [58]

    Sharma, P

    A. Sharma, P. Suryanarayana, Real-space density functional theory adapted to cyclic and helical symmetry: Application to torsional deformation of carbon nanotubes, Phys. Rev. B. 103 (3) (2021) 035101

    2021

  51. [59]

    Gavini, J

    V. Gavini, J. Knap, K. Bhattacharya, M. Ortiz, Non-periodic finite-element formulation of orbital-free density functional theory, J. Mech. Phys. Solids 55 (4) (2007) 669 – 696

    2007

  52. [60]

    Alemany, M

    M. Alemany, M. Jain, L. Kronik, J. R. Chelikowsky, Real-space pseudopotential method for computing the electronic properties of periodic systems, Phys. Rev. B 69 (7) (2004) 075101

    2004

  53. [61]

    Hirose, T

    K. Hirose, T. Ono, Y. Fujimoto, S. Tsukamoto, First-principles calculations in real-space formalism (2005)

    2005

  54. [62]

    J. Han, M. L. Tiago, T.-L. Chan, J. R. Chelikowsky, Real space method for the electronic structure of one-dimensional periodic systems, J. Chem. Phys. 129 (14) (2008) 144109

    2008

  55. [63]

    Natan, A

    A. Natan, A. Benjamini, D. Naveh, L. Kronik, M. L. Tiago, S. P. Beckman, J. R. Chelikowsky, Real-space pseudopotential method for first principles calculations of general periodic and partially periodic systems, Phys. Rev. B—Condensed Matter and Materials Physics 78 (7) (2008) 075109

    2008

  56. [64]

    K.Ramakrishnan,G.SaiGautam,P.Motamarri,Real-spacemethodsforabinitiomodelingofsurfacesandinterfacesunderexternalpotential bias, J. Chem. Theory Comput. 21 (14) (2025) 7087–7101

    2025

  57. [65]

    Q. Xu, A. Sharma, B. Comer, H. Huang, E. Chow, A. J. Medford, J. E. Pask, P. Suryanarayana, Sparc: Simulation package for ab-initio real-space calculations, SoftwareX 15 (2021) 100709

    2021

  58. [66]

    Impacts 20 (2024) 100649

    B.Zhang,X.Jing,Q.Xu,S.Kumar,A.Sharma,L.Erlandson,S.J.Sahoo,E.Chow,A.J.Medford,J.E.Pask,etal.,Sparcv2.0.0:Spin-orbit coupling, dispersion interactions, and advanced exchange–correlation functionals, Softw. Impacts 20 (2024) 100649

    2024

  59. [67]

    Lennard-Jones, B

    J. Lennard-Jones, B. M. Dent, Cohesion at a crystal surface, Trans. Faraday Soc. 24 (1928) 92–108

    1928

  60. [68]

    Abramowitz, I

    M. Abramowitz, I. A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, National Bureau of Standards, Washington, D.C., 1964

    1964

  61. [69]

    Bhowmik, A

    S. Bhowmik, A. J. Medford, P. Suryanarayana, Bulk boundary condition for surface calculations in density functional theory, arXiv preprint arXiv:2607.07894 (2026)

    2026 arXiv

  62. [70]

    N. D. Mermin, Thermal properties of the inhomogeneous electron gas, Phys. Rev. 137 (5A) (1965) A1441

    1965

  63. [71]

    Kleinman, D

    L. Kleinman, D. Bylander, Efficacious form for model pseudopotentials, Phys. Rev. Lett. 48 (20) (1982) 1425

    1982

  64. [72]

    Harris, Simplified method for calculating the energy of weakly interacting fragments, Phys

    J. Harris, Simplified method for calculating the energy of weakly interacting fragments, Phys. Rev. B 31 (4) (1985) 1770

    1985

  65. [73]

    W. M. C. Foulkes, R. Haydock, Tight-binding models and density-functional theory, Phys. Rev. B 39 (17) (1989) 12520

    1989

  66. [74]

    R. P. Feynman, Forces in molecules, Phys. Rev. 56 (4) (1939) 340

    1939

  67. [75]

    A.Sharma,P.Suryanarayana,Onthecalculationofthestresstensorinreal-spacekohn-shamdensityfunctionaltheory,J.Chem.Phys.149(19) (2018) 194104

    2018

  68. [76]

    Sharma, A

    A. Sharma, A. Metere, P. Suryanarayana, L. Erlandson, E. Chow, J. E. Pask, Gpu acceleration of local and semilocal density functional calculations in the sparc electronic structure code, J. Chem. Phys. 158 (20) (2023)

    2023

  69. [77]

    X. Jing, A. Sharma, J. E. Pask, P. Suryanarayana, Gpu acceleration of hybrid functional calculations in the sparc electronic structure code, J. Chem. Phys. 162 (18) (2025)

    2025

  70. [78]

    Suryanarayana, D

    P. Suryanarayana, D. Phanish, Augmented lagrangian formulation of orbital-free density functional theory, J. Comput. Phys. 275 (2014) 524–538

    2014

  71. [79]

    Y. Zhou, Y. Saad, M. L. Tiago, J. R. Chelikowsky, Self-consistent-field calculations using chebyshev-filtered subspace iteration, J. Comput. Phys. 219 (1) (2006) 172–184

    2006

  72. [80]

    Y.Zhou,Y.Saad,M.L.Tiago,J.R.Chelikowsky,Parallelself-consistent-fieldcalculationsviachebyshev-filteredsubspaceacceleration,Phys. Rev. E 74 (6) (2006) 066704

    2006

  73. [81]

    P. P. Pratapa, P. Suryanarayana, Restarted pulay mixing for efficient and robust acceleration of fixed-point iterations, Chem. Phys. Lett. 635 (2015) 69–74

    2015

  74. [82]

    A. S. Banerjee, P. Suryanarayana, J. E. Pask, Periodic pulay method for robust and efficient convergence acceleration of self-consistent field iterations, Chem. Phys. Lett. 647 (2016) 31 – 35. Rajat et al.: Preprint submitted to arXiv Page 20 of 21 Unified open-boundary electr...

    2016

  75. [83]

    S.Kumar,Q.Xu,P.Suryanarayana,Onpreconditioningtheself-consistentfielditerationinreal-spacedensityfunctionaltheory,Chem.Phys. Lett. 739 (2020) 136983

    2020

  76. [84]

    P.Suryanarayana,P.P.Pratapa,J.E.Pask,Alternatinganderson–richardsonmethod:Anefficientalternativetopreconditionedkrylovmethods for large, sparse linear systems, Comput. Phys. Commun. 234 (2019) 278–285

    2019

  77. [85]

    P.P.Pratapa,P.Suryanarayana,J.E.Pask,Andersonaccelerationofthejacobiiterativemethod:Anefficientalternativetokrylovmethodsfor large, sparse linear systems, J. Comput. Phys. 306 (2016) 43–54

    2016

  78. [86]

    W. H. Press, Numerical recipes 3rd edition: The art of scientific computing, Cambridge university press, 2007

    2007

  79. [87]

    H. J. Monkhorst, J. D. Pack, Special points for brillouin-zone integrations, Phys. Rev. B 13 (12) (1976) 5188

    1976

  80. [88]

    Hamann, Optimized norm-conserving vanderbilt pseudopotentials, Phys

    D. Hamann, Optimized norm-conserving vanderbilt pseudopotentials, Phys. Rev. B 88 (8) (2013) 085117

    2013

  81. [89]

    M.F.Shojaei,J.E.Pask,A.J.Medford,P.Suryanarayana,Softandtransferablepseudopotentialsfrommulti-objectiveoptimization,Comput. Phys. Commun. 283 (2023) 108594

    2023

  82. [90]

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

    1996

  83. [91]

    A. C. Riis-Jensen, T. Deilmann, T. Olsen, K. S. Thygesen, Classifying the electronic and optical properties of janus monolayers, ACS nano 13 (11) (2019) 13354–13364

    2019

  84. [92]

    L.Dong,J.Lou,V.B.Shenoy,Largein-planeandverticalpiezoelectricityinjanustransitionmetaldichalchogenides,ACSnano11(8)(2017) 8242–8248

    2017

  85. [93]

    D.Codony,I.Arias,P.Suryanarayana,Transversalflexoelectriccoefficientfornanostructuresatfinitedeformationsfromfirstprinciples,Phys. Rev. Mater. 5 (3) (2021) L030801

    2021

  86. [94]

    Kumar, D

    S. Kumar, D. Codony, I. Arias, P. Suryanarayana, Flexoelectricity in atomic monolayers from first principles, Nanoscale 13 (3) (2021) 1600– 1607

    2021

  87. [95]

    Kumar, P

    S. Kumar, P. Suryanarayana, Bending moduli for forty-four select atomic monolayers from first principles, Nanotechnology 31 (43) (2020) 43LT01

    2020

  88. [96]

    Bhardwaj, A

    A. Bhardwaj, A. Sharma, P. Suryanarayana, Torsional strain engineering of transition metal dichalcogenide nanotubes: an ab initio study, Nanotechnology 32 (47) (2021) 47LT01

    2021

  89. [97]

    Sharma, P

    A. Sharma, P. Suryanarayana, Calculation of phonons in real-space density functional theory, Phys. Rev. E 108 (4) (2023) 045302

    2023

  90. [98]

    Sharma, P

    A. Sharma, P. Suryanarayana, Cyclic-and helical-symmetry-adapted phonon formalism within density functional perturbation theory, Phys. Rev. B 113 (20) (2026) 205116. Rajat et al.: Preprint submitted to arXiv Page 21 of 21

    2026

Pith tools

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