Pith. sign in

REVIEW 3 major objections 5 minor 82 references

Exchange-Correlation Potentials and Energy Densities through Orbital Averaging and Aufbau Integration

T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

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

Pith's one-line read The paper shows that orbital averaging over Kohn-Sham orbitals derived from FCI densities in a Slater basis recovers exchange-correlation potentials with correct $-1/r$ decay and integer-electron step for He through Ne, and that…

desk verdict Useful benchmark workflow with an imposed asymptote and a KS-orbital assumption that needs more testing. read the letter →

arxiv 2502.10262 v1 pith:YIGY6FEC submitted 2025-02-14 physics.chem-ph

Vaibhav Khanna , Bikash Kanungo , Jeffrey Hatch , Joshua Kammeraad , Paul M. Zimmerman This is my paper · ORCID
classification physics.chem-ph PACS 31.15.E
keywords Kohn-Shaminversionexchange-correlationpotentialenergydensityorbitalaveragingaufbaupathfullconfigurationinteractionSlaterbasisderivativediscontinuity
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 aims to establish a complete route from correlated wavefunction densities to accurate Kohn-Sham exchange-correlation quantities: the potential $v_{\rm xc}$, the integer-electron step in that potential, and a spatially resolved energy density $e_{\rm xc}$. The input densities are full configuration interaction (practically heat-bath CI) results for the atoms He through Ne in a large Slater basis, converted into Kohn-Sham orbitals by balancing density fidelity against kinetic energy. Averaging the Kohn-Sham equations over these orbitals gives potentials with the correct $-1/r$ asymptotic decay and the discontinuous step at integer electron number; integrating those potentials along a piecewise-linear aufbau path between integer densities yields energy densities whose integrals match the CI $E_{\rm xc}$ to better than 0.1%. This matters because it supplies reference data of a kind approximate density functionals currently lack, directly targeting the derivative discontinuity and the spatial shape of the exchange-correlation energy.

What carries the argument

The machinery is the orbital-averaged inversion identity combined with a path integral. The identity replaces the unknown $v_{\rm xc}$ by a density-weighted average over occupied orbitals, so nodes of individual orbitals never divide the expression; the step basis $\{b_i = n_i|\phi_i|^2/\rho\}$ separates the constant shift of the potential from its spatial shape. The q-aufbau path integral $\int_0^1 dq\, v_{\rm xc}([\rho_q];r)\,\Delta\rho_m(r)$ supplies the energy density for each integer interval, and the sum over intervals gives an $e_{\rm xc}$ whose integral recovers $E_{\rm xc}$.

What would settle it

A single decisive check: use the constructed $v_{\rm xc}$ as input to a self-consistent Kohn-Sham calculation in the same Slater basis and compare the output density to the original FCI density. If the potential is the true KS potential, the density should match at least as well as the fitted orbitals do; a noticeably larger mismatch would show the orbital-averaged potential is not the KS potential. A second check: repeat the whole workflow with a substantially larger Slater basis; if the potential changes materially beyond the region of the $\lambda$ variation, the claimed potential is not converged.

Watch

Extended reading notes

Core claim

The paper's central claim is that the orbital-averaged inversion formula $$v_{\rm xc}(r) = \frac{\sum_i n_i\left(\epsilon_i|\phi_i(r)|^2 + \frac{1}{2}\phi_i^*(r)\$nabla^{2}$\phi_i(r)\right)}{\rho(r)} - v_{\rm ext}(r) - v_{\rm H}(r)$$ turns CI-derived Slater-basis orbitals into the physical exchange-correlation potential. The authors obtain the KS eigenvalues by projecting the KS equation onto virtual orbitals, pin the HOMO eigenvalue to the negative ionization energy so that the potential decays as $-1/r$, and isolate the uniform step at integer electron number through a step basis. Along each interval $m-1\to m$, they integrate these potentials over the q-aufbau path $\rho_q(r)=\rho_{m-1}(r)+q(\rho_m(r)-\rho_{m-1}(r))$, producing a spatially resolved $e_{\rm xc}$ that integrates to $E_{\rm xc}$ within $10^{-3}$ relative error for He through Ne.

Load-bearing premise

The load-bearing premise is that the orbitals produced by density fitting with a manually chosen kinetic-energy weight are close enough to true Kohn-Sham orbitals that the orbital-averaged potential is the physical exchange-correlation potential; the density residual alone does not guarantee this.

Editorial extensions

If this is right

  • The recovered potentials give HOMO eigenvalues that match ionization energies to better than $10^{-6}$ Ha for every atom studied, so the reconstructed potential carries the correct long-range physics.
  • The step-basis contribution to $E_{\rm xc}$ follows a quadratic curve in $Z$ ($R^2 = 0.9966$), giving a concrete target for approximate functionals that currently smear the integer-electron discontinuity.
  • Because the q-aufbau energy densities integrate to the CI $E_{\rm xc}$ within 0.1%, they provide spatially resolved reference $e_{\rm xc}$ values that functional approximations can be trained against.
  • The workflow extends to the full sequence of electron counts, so the same construction yields $v_{\rm xc}$ and $e_{\rm xc}$ for cations and fractional occupations, not just neutral atoms.
  • The comparison with a semilocal functional suggests that the energy density from the q-aufbau path is qualitatively similar to the functional's energy density even where the potentials differ strongly, which would justify using such densities as a gauge-matched training target.

Reading between the lines

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

  • The paper does not test molecules; a direct extension would apply the same orbital-averaged workflow to stretched bonds, where the derivative discontinuity should appear as a spatial interatomic step rather than the uniform shift seen in atoms.
  • A stricter reconstruction test would feed the reconstructed $v_{\rm xc}$ back into a self-consistent Kohn-Sham calculation and verify that the resulting density reproduces the FCI density; the current density residual of $\sim 10^{-3}$ per electron does not by itself prove the potential is the KS potential.
  • The quadratic scaling of the step contribution with $Z$ suggests a simple $Z$-dependent correction that could be added to local and semilocal functionals; whether that correction transfers to molecules is an open question.
  • The virial consistency test is reported for only He, Be, and N; extending it to all atoms He through Ne would be a quick check of whether the near-nucleus potential errors seen at larger $\lambda$ are fully absent at the chosen $\lambda$.
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 / 5 minor

Summary. The paper presents a workflow for constructing Kohn-Sham exchange-correlation potentials vxc(r) and energy densities exc(r) from correlated wavefunction densities. The authors use the Rask et al. procedure to obtain KS orbitals from CI densities in a Slater basis, then apply an orbital-averaged inversion (Eq. 13) with eigenvalues determined by projection onto virtual orbitals and a HOMO condition (Eq. 19). A smooth switch to the Slater potential (Eq. 21) handles the low-density tail, and a step-basis decomposition extracts the constant component of vxc. Energy densities are obtained by integrating vxc along a piecewise q-aufbau path. Results are reported for He through Ne at integer and fractional electron counts, including the -1/r asymptotic tail, the integer-electron step, virial tests, and Exc integrals with relative errors below 0.1%.

Significance. If substantiated, this workflow would provide a practical route to benchmark approximate exchange-correlation potentials and energy densities for atoms using correlated densities, complementing methods such as KZG inversion and SlaterRKS. The paper has several concrete strengths: a large Slater basis that avoids the oscillations typical of Gaussian-basis inversion, a direct comparison of the Be potential against an independent KZG reference (Fig. 3), a virial consistency test (Table S2), and a closed-form path-integral construction of exc whose integral agrees with CI Exc to within 0.1% (Table 3). The step-basis decomposition of the constant contribution to vxc is an interesting extension. However, the central claim that the orbital-averaged vxc is the physical Kohn-Sham potential for all atoms He-Ne is only partially supported, because the underlying Rask orbitals are benchmarked directly for only one atom, and two headline features (the -1/r tail and the HOMO eigenvalue) are enforced by construction rather than demonstrated independently.

major comments (3)
  1. [Computing Eigenvalues ϵj (Eq. 19) and Results (Fig. 4)] The claimed asymptotic -1/r decay is not an independent finding: Eq. (19) sets ε_HOMO = -I, which by construction imposes the correct long-range decay of the Kohn-Sham potential, and Eq. (21)-(22) replace the orbital-averaged vxc with the Slater potential wherever ρ(r) < θ = 10^-5, which is exactly the region displayed in Fig. 4. The reported |IP - ε_HOMO| < 10^-6 (Table 3) is therefore a consistency condition of the imposed equation, not a validation of the tail. The authors should either demonstrate the -1/r behavior in a region where the switch (Eq. 21) is inactive and where ε_HOMO has not been fixed by Eq. (19), or explicitly state that the asymptotic tail is enforced by construction rather than fitted.
  2. [Obtaining KS Orbitals from FCI Densities (Eq. 5) and Results (Fig. 3, Table S2)] The orbital-averaged vxc in Eq. (13) inherits the quality of the Rask orbitals, yet the evidence that these orbitals are close to the true KS orbitals is limited to a density residual of order 10^-4-10^-3 per electron (Table 3), a single KZG comparison for Be (Fig. 3), and virial tests for He, Be, and N (Table S2). Density residuals constrain the density-weighted average of the potential but not the pointwise potential; a systematic error in vxc for C, O, F, or Ne that preserves the density and the virial of the tested atoms is not excluded. To support the claim that the OA potential is the physical KS potential across the series, the authors should provide at least one additional independent comparison (e.g., Ne against KZG) or a sensitivity analysis showing that the potential is robust to the choice of λ beyond the density-based arguments.
  3. [Line Integration for Energy Densities (Table 3, Fig. 7)] The <0.1% relative error in Exc obtained by integrating exc is a self-consistency check, not a pointwise validation of the energy density. Since ρ_q, vxc([ρ_q]), and the CI reference Exc are all derived from the same CI densities, the test confirms that the path integral of the reconstructed potential reproduces the input energy, but it does not rule out local errors in exc(r) that cancel under integration. The paper should temper the claim that the energy densities are accurate to 0.1% and, if possible, add a pointwise comparison to a gauge-matched reference or an alternative path to confirm that the spatial structure of exc is meaningful.
minor comments (5)
  1. [Abstract and Computational Details] The abstract repeatedly states that the reference densities come from 'full configuration interaction in a Slater orbital basis,' whereas the Computational Details section describes heat-bath configuration interaction (HBCI) with selection thresholds of 10^-5 Ha (He-B) and 5x10^-5 Ha (C-Ne). Please qualify the abstract to avoid overstating the reference level, or provide evidence that the HBCI wavefunctions are numerically indistinguishable from FCI for the reported quantities.
  2. [Table 3] The table header contains the typo 'Maxium' for 'Maximum'; please correct it. Also, the table lists only the maximum |R| per occupied orbital; providing the range or an error bar would help assess the consistency of the eigenvalue solution.
  3. [Capturing the Step Basis Contribution to vxc] The transformation from {b_i} to {1, b'_2, ..., b'_N} is only stated in words; specifying the linear transformation and the resulting expression for ϵ'_1 explicitly would make the step-basis extraction reproducible without consulting the SI.
  4. [Results and Discussion, Figure 6] The quadratic fit for Estep_xc is reported with R² = 0.9966 but was performed 'while imposing a penalty on positive slopes.' Please report the unpenalized fit and the penalty parameter, or justify why the unconstrained fit is less appropriate.
  5. [Methods and Theory, Eq. (20)] The fractional-occupation interpolation in Eq. (20) is a specific ensemble choice; the text notes this is 'neither unique nor necessary,' but the dependence of the step height and of Estep_xc on this choice should be mentioned, since the step claim relies on this ensemble construction.

Circularity Check

2 steps flagged · score 6.0 of 10

Two headline validations — the −1/r tail and the <1e-6 ionization-energy error — are enforced inputs (Eqs. 19 and 21), although the core density-to-potential inversion and the virial/KZG checks are independent.

  1. fitted input called prediction [Methods and Theory, Eq. 19; Results, Table 3 caption]
    "Here ϵHOMO is set to be equal to −I, (the ionization energy), to ensure that the exchange-correlation potential decays as -1/r at large distances. 42–44 This additional equation is included in the linear system of equations and solved in the same step as equation 18. ... Across the He-Ne series, the error in the predicted ionization energy,|ICI|−| ϵKS HOMO|, was found to be < 10−6."

    Eq. 19 is an explicit constraint fixing ϵ_HOMO = −I in the same least-squares system that determines the other eigenvalues. Therefore the tabulated quantity |I_CI| − |ϵ_HOMO| is not a prediction of the ionization energy: any nonzero value would indicate only that the constraint was not satisfied in the solve. The sentence in the caption calling this a 'predicted ionization energy' presents an input as an output. The same constraint is also the stated justification for the −1/r asymptotic decay, so that part of the physical claim is imported from the input I, not demonstrated by the inversion.

  2. self definitional [Methods and Theory, Eqs. 21–22; Results, Fig. 4 and asymptotic discussion]
    "Whenρ(r) vanishes at large radial distances from nuclei, the orbital-averagedvxc(r), derived through division by ρ(r), may become numerically unstable. However, it is known that at large r, the exchange-correlation potential vxc and the Slater exchange-correlation potential vWF xc,Slater are expected to be identical. 32,33,35 Hence, a smooth transition is introduced between vxc and vWF xc,Slater ... At large r, all computed potentials exhibit the expected asymptotic behavior, decaying as−1/r, in agreement with the known asymptotic of the potential."

    The transition function F(r)=ρ/(ρ+θ) with θ=10^-5 makes the output vxc equal to vWF_xc,Slater wherever ρ ≪ 10^-5. Since this is precisely the region shown in Fig. 4 at large r, the reported −1/r decay is the known decay of the substituted Slater hole potential, not a property independently recovered by the orbital-averaged Eq. 13 inversion. The paper's statement that 'all computed potentials exhibit the expected asymptotic behavior' is thus guaranteed by the construction of Eq. 21, and cannot serve as a validation of the OA potential in that region.

full rationale

The core of the paper is a genuine inversion: the orbital-averaged vxc of Eq. 13 is determined by FCI-derived densities and Rask KS orbitals, not fitted to target vxc or Exc values, and the energy-density line integral integrates the same vxc. The density residuals (Table 3), virial test (Table S2), KZG comparison for Be, and PBE self-consistency check are non-circular evidence. However, two results advertised as successes are imposed by construction. Eq. 19 pins ϵ_HOMO to −I, so the reported IE error is a tautology; Eq. 21 replaces vxc with the Slater potential wherever ρ is below about 10^-5, so the displayed −1/r decay is pre-defined in exactly the asymptotic region used for validation. The fractional-occupation ensemble is acknowledged as non-unique, and the reliance on Rask orbitals from the same group is a reproducibility or correctness concern rather than a circular one. Overall, the central inversion method has independent content, but the paper's headline feature claims are partially circular.

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

The workflow avoids inventing new physical entities but rests on approximations confirmed only indirectly: the density-fitting lambda is a tuned parameter, the fractional occupations are a stated modeling choice, and the -1/r and IP behaviors are enforced by constraint. The most independent evidence is the virial test (He, Be, N: within 4.3 mHa) and the agreement with the KZG potential for Be. No code or raw data are released, which limits independent verification beyond re-implementation.

free parameters (3)
  • lambda (kinetic-energy mixing weight in Rask orbital optimization) = 0.00005 for He-O; 0.0001 for F, Ne
    Chosen by examining the resulting densities and potentials (SI Section S4). The choice is per-element and manual; higher lambda degrades the density fit and produces an unphysical near-nucleus feature. This is a genuine free parameter, not fixed by any theory.
  • fractional p-orbital occupations = (1-delta)/3 on each of px, py, pz
    The paper states this ensemble choice is 'neither unique nor necessary' but is used to maintain spherical symmetry. The step and energy densities depend on the occupation switching.
  • theta (smoothing parameter in the vxc to Slater-potential transition) = 1e-5
    Numerical stabilization constant controlling where the potential transitions to the Slater form at large r (Eq. 22). It does not affect the short-range potential but sets the location of the asymptotic regime.
assumptions (4)
  • domain assumption Existence and accessibility of the KS system: the Rask et al. density-fit procedure (Eq. 5) yields orbitals that are the Kohn-Sham orbitals of the FCI density with a controllable error.
    The entire workflow leans on this. The L1 density errors are ~10^-3 per electron, so the assumption is approximate. If the Rask orbitals were not close to true KS orbitals, the orbital-averaged vxc of Eq. 13 would not be the physical KS potential.
  • standard math The Levy constrained-search definition of the KS kinetic energy and the van Leeuwen-Baerends line integral theorem (Eq. 25) are exact.
    These underpin the potential definition and the energy-density integration. They are standard results in exact DFT and are not questioned here.
  • domain assumption The exact vxc decays as -1/r and the HOMO eigenvalue equals -I (ionization energy).
    This standard asymptotic result is turned into a constraint (Eq. 19): epsilon_HOMO is set to -I by construction. It guarantees the -1/r tail but makes the reported <10^-6 IP error a direct consequence of the constraint.
  • domain assumption Near-FCI HBCI densities with the stated selection thresholds (10^-5 Ha for He-B, 5x10^-5 Ha for C-Ne) are accurate stand-ins for exact ground-state densities.
    Stated in Computational Details. The abstract calls this 'full configuration interaction'; the actual method is heat-bath CI with a tight threshold.
invented entities (1)
  • Step basis and step-basis energy Estep_xc
    purpose: Decomposes vxc into a spatially constant part (the integer-electron step) and a shape part, and extracts the corresponding energy contribution.
    This is an analysis decomposition, not a new physical object. It defines a gauge-like choice that affects the numerical value of Estep_xc, which is then fit quadratically in Z. The fitted curve (Estep_xc = -6.13 Z^2 + 27.78 Z - 30.57, R^2=0.9966) is a descriptive fit to eight data points, not a falsifiable prediction.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Exchange-Correlation Potentials and Energy Densities through Orbital Averaging and Aufbau Integration." pith.science (2026). https://pith.science/paper/YIGY6FEC

@misc{pith2026250210262,
  author       = {Pith},
  title        = {Pith review of: Exchange-Correlation Potentials and Energy Densities through Orbital Averaging and Aufbau Integration},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/YIGY6FEC}},
  note         = {Machine review of arXiv:2502.10262}
}
read the original abstract

Exchange-correlation potentials vxc and energy densities exc are derived for integer and fractional electron counts using an orbital-averaged Kohn-Sham inversion procedure. The reference densities for inversion come from full configuration interaction in a Slater orbital basis. The orbital-averaged potentials accurately capture key features of vxc, including the asymptotic negative one over r decay and the step discontinuity associated with integer electron transitions for the series of atoms He through Ne. Exchange-correlation energy densities exc are produced through an aufbau path integral. The energy densities reach good agreement with total Exc values. By providing full configuration interaction-derived Kohn-Sham quantities, including vxc, exc, and step contributions, this workflow can be instrumental in the development of improved XC functionals that bridge wavefunction-level accuracy with the computational efficiency of density functional theory.

Figures

Figures reproduced from arXiv: 2502.10262 by the authors.

Figure 1
Figure 1. Workflow for computing the exchange correlation potential, step basis contribution, [PITH_FULL_IMAGE:figures/full_fig_p014_1.png] view at source ↗
Figure 2
Figure 2. Difference in the FCI and KS densities for Be, B, N and Ne as a function of the [PITH_FULL_IMAGE:figures/full_fig_p018_2.png] view at source ↗
Figure 3
Figure 3. Comparison of the CI-derived orbital averaged exchange-correlation potential for [PITH_FULL_IMAGE:figures/full_fig_p019_3.png] view at source ↗
Figures from the paper (5 more)
Figure 4
Figure 4. Figure 4: shows the computed exchange-correlation potential for Be, B, N, and Ne as a function of log10 r. At large r, all computed potentials exhibit the expected asymptotic behavior, decaying as −1/r, in agreement with the known asymptotic of the potential. 42,43 This decay of…
Figure 5
Figure 5. Figure 5: Three-dimensional plots of the exchange-correlation potential ( [PITH_FULL_IMAGE:figures/full_fig_p022_5.png]
Figure 6
Figure 6. Figure 6: Variation of the step basis contribution to the exchange-correlation energy with [PITH_FULL_IMAGE:figures/full_fig_p023_6.png]
Figure 7
Figure 7. Figure 7: Plots depicting the variation of a) log10 |exc| and b) log10 ρ for Be, B, N, Ne with the radial distance (r). error |∆Exc/ECI xc | remains consistently below 0.1%, demonstrating the high accuracy of the computed exchange-correlation energy densities for the He to Ne se…
Figure 8
Figure 8. Figure 8: Scatter plots of log10 |∇ρ| versus log10 ρ for Be, B, N, Ne. The colormaps represent log10 |exc| values from a) PBE and b) FCI. Select exc values for Ne have been labeled on both plots. KS inversions using other finite basis sets. By applying this framework to the He-N…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

82 extracted references · 77 canonical work pages

  1. [1]

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

    2014

  2. [2]

    J.; Truhlar, D

    Cramer, C. J.; Truhlar, D. G. Density functional theory for transition metals and transition metal chemistry. Phys. Chem. Chem. Phys. 2009, 11, 10757--10816

    2009

  3. [3]

    T.; Fielicke, A.; Janssens, E.; Lievens, P.; McGrady, J

    Khanna, V.; Singh, R.; Claes, P.; Nguyen, M. T.; Fielicke, A.; Janssens, E.; Lievens, P.; McGrady, J. E. Evolution of Vibrational Spectra in the Manganese–Silicon Clusters Mn2Sin, n = 10, 12, and 13, and Cationic [Mn2Si13]+. J. Phys. Chem. A 2022, 126, 1617--1626

    work page 2022

  4. [4]

    P.; Ruzsinszky, A.; Csonka, G

    Perdew, J. P.; Ruzsinszky, A.; Csonka, G. I.; Vydrov, O. A.; Scuseria, G. E.; Constantin, L. A.; Zhou, X.; Burke, K. Restoring the Density-Gradient Expansion for Exchange in Solids and Surfaces. Phys. Rev. Lett. 2008, 100, 136406

    work page 2008

  5. [5]

    J.; Refson, K.; Probert, M

    Hasnip, P. J.; Refson, K.; Probert, M. I. J.; Yates, J. R.; Clark, S. J.; Pickard, C. J. Density functional theory in the solid state. Philos. Trans. R. Soc. A 2014, 372, 20130270

    work page 2014

  6. [6]

    Khanna, V.; Kanungo, B.; Gavini, V.; Tewari, A.; Zimmerman, P. M. Examining the Impact of Local Constraint Violations on Energy Computations in DFT. J. Comput. Chem. 2025, 46, e70005

    work page 2025

  7. [7]

    Teale, A. M. et al. DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science. Phys. Chem. Chem. Phys. 2022, 24, 28700--28781

    work page 2022

  8. [8]

    Kohn, W.; Sham, L. J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133--A1138

    work page 1965

Show all 82 references
  1. [9]

    Thirty years of density functional theory in computational chemistry: an overview and extensive assessment of 200 density functionals

    Mardirossian, N.; Head-Gordon, M. Thirty years of density functional theory in computational chemistry: an overview and extensive assessment of 200 density functionals. Mol. Phys. 2017, 115, 2315--2372

    2017

  2. [10]

    Direct Method for Optimized Effective Potentials in Density-Functional Theory

    Yang, W.; Wu, Q. Direct Method for Optimized Effective Potentials in Density-Functional Theory. Phys. Rev. Lett. 2002, 89, 143002

    2002

  3. [11]

    I.; Lathiotakis, N

    Gidopoulos, N. I.; Lathiotakis, N. N. Nonanalyticity of the optimized effective potential with finite basis sets. Phys. Rev. A 2012, 85, 052508

    2012

  4. [12]

    Lieb-response

    Gould, T. Toward routine Kohn–Sham inversion using the “Lieb-response” approach. J. Chem. Phys. 2023, 158, 064102

    2023

  5. [13]

    Kohn-Sham potentials in exact density-functional theory at noninteger electron numbers

    Gould, T.; Toulouse, J. Kohn-Sham potentials in exact density-functional theory at noninteger electron numbers. Phys. Rev. A 2014, 90, 050502

    2014

  6. [14]

    Violations of the v -representability condition underlying Kohn-Sham density-functional theory

    Trushin, E.; Erhard, J.; G\"orling, A. Violations of the v -representability condition underlying Kohn-Sham density-functional theory. Phys. Rev. A 2024, 110, L020802

    2024

  7. [15]

    Kohn–Sham inversion for open-shell systems

    Erhard, J.; Trushin, E.; Görling, A. Kohn–Sham inversion for open-shell systems. J. Chem. Phys. 2025, 162, 034116

    2025

  8. [16]

    olkner, A.; Erhard, J.; G\

    Trushin, E.; Fauser, S.; M\"olkner, A.; Erhard, J.; G\"orling, A. Accurate Correlation Potentials from the Self-Consistent Random Phase Approximation. Phys. Rev. Lett. 2025, 134, 016402

    2025

  9. [17]

    A.; Heaton-Burgess, T.; Cohen, A

    Bulat, F. A.; Heaton-Burgess, T.; Cohen, A. J.; Yang, W. Optimized effective potentials from electron densities in finite basis sets. J. Chem. Phys. 2007, 127, 174101

    2007

  10. [18]

    Optimized effective potentials from arbitrary basis sets

    Heaton-Burgess, T.; Yang, W. Optimized effective potentials from arbitrary basis sets. J. Chem. Phys. 2008, 129, 194102

    2008

  11. [19]

    Inverse Kohn–Sham Density Functional Theory: Progress and Challenges

    Shi, Y.; Wasserman, A. Inverse Kohn–Sham Density Functional Theory: Progress and Challenges. J. Phys. Chem. Lett. 2021, 12, 5308--5318

    2021

  12. [20]

    G.; Ospadov, E.; Staroverov, V

    Ryabinkin, I. G.; Ospadov, E.; Staroverov, V. N. Exact exchange-correlation potentials of singlet two-electron systems. J. Chem. Phys. 2017, 147, 164117

    2017

  13. [21]

    Effective potentials in density-functional theory

    Aryasetiawan, F.; Stott, M. Effective potentials in density-functional theory. Phys. Rev. B 1988, 38, 2974--2981

    1988

  14. [22]

    Kohn-Sham potentials and wave functions from electron densities

    Görling, A. Kohn-Sham potentials and wave functions from electron densities. Phys. Rev. A 1992, 46, 3753--3757

    1992

  15. [23]

    Zhao, Q.; Parr, R. G. Quantities \( T_s[n] \) and \( T_c[n] \) in density-functional theory. Phys. Rev. A 1992, 46, 2337--2339

    1992

  16. [24]

    Zhao, Q.; Parr, R. G. Constrained-search method to determine electronic wave functions from electronic densities. J. Chem. Phys. 1993, 98, 543--548

    1993

  17. [25]

    S.; Wasserman, A

    Jensen, D. S.; Wasserman, A. Numerical methods for the inverse problem of density functional theory. Int. J. Quantum Chem. 2017, 118, e25425

    2017

  18. [26]

    Wang, Y.; Parr, R. G. Construction of exact Kohn-Sham orbitals from a given electron density. Phys. Rev. A 1993, 47, R1591--R1593

    1993

  19. [27]

    Algorithm to derive exact exchange-correlation potentials from correlated densities in atoms

    Peirs, K.; Van Neck, D.; Waroquier, M. Algorithm to derive exact exchange-correlation potentials from correlated densities in atoms. Phys. Rev. A 2003, 67, 012505

    2003

  20. [28]

    Exact exchange-correlation potentials from ground-state electron densities

    Kanungo, B.; Zimmerman, P.; Gavini, V. Exact exchange-correlation potentials from ground-state electron densities. Nat. Commun. 2019, 10, 4497

    2019

  21. [29]

    M.; Gavini, V

    Kanungo, B.; Tribedi, S.; Zimmerman, P. M.; Gavini, V. Accelerating inverse Kohn-Sham calculations using reduced density matrices. 2024; https://arxiv.org/abs/2408.02342

    2024 arXiv

  22. [30]

    M.; Gavini, V

    Kanungo, B.; Hatch, J.; Zimmerman, P. M.; Gavini, V. Exact and Model Exchange-Correlation Potentials for Open-Shell Systems. J. Phys. Chem. Lett. 2023, 14, 10039--10045

    2023

  23. [31]

    B.; Bischoff, F

    St \"u ckrath, J. B.; Bischoff, F. A. Reduction of Hartree–Fock Wavefunctions to Kohn–Sham Effective Potentials Using Multiresolution Analysis. J. Chem. Theory Comput. 2021, 17, 1408--1420

    2021

  24. [32]

    G.; Kohut, S

    Ryabinkin, I. G.; Kohut, S. V.; Staroverov, V. N. Reduction of Electronic Wave Functions to Kohn-Sham Effective Potentials. Phys. Rev. Lett. 2015, 115, 083001

    2015

  25. [33]

    W.; Staroverov, V

    Cuevas-Saavedra, R.; Ayers, P. W.; Staroverov, V. N. Kohn–Sham exchange-correlation potentials from second-order reduced density matrices. J. Chem. Phys. 2015, 143, 244116

    2015

  26. [34]

    G.; Staroverov, V

    Ospadov, E.; Ryabinkin, I. G.; Staroverov, V. N. Improved method for generating exchange-correlation potentials from electronic wave functions. J. Chem. Phys. 2017, 146, 084103

    2017

  27. [35]

    Tribedi, S.; Dang, D.-K.; Kanungo, B.; Gavini, V.; Zimmerman, P. M. Exchange correlation potentials from full configuration interaction in a Slater orbital basis. J. Chem. Phys. 2023, 159, 054106

    2023

  28. [36]

    van Leeuwen, R.; Baerends, E. J. Exchange-correlation potential with correct asymptotic behavior. Phys. Rev. A 1994, 49, 2421--2431

    1994

  29. [37]

    V.; van Leeuwen, R.; Baerends, E

    Gritsenko, O. V.; van Leeuwen, R.; Baerends, E. J. Molecular Kohn-Sham exchange-correlation potential from the correlated ab initio electron density. Phys. Rev. A 1995, 52, 1870--1874

    1995

  30. [38]

    A.; Kohut, S

    Kananenka, A. A.; Kohut, S. V.; Gaiduk, A. P.; Ryabinkin, I. G.; Staroverov, V. N. Efficient construction of exchange and correlation potentials by inverting the Kohn–Sham equations. J. Chem. Phys. 2013, 139, 074112

    2013

  31. [39]

    E.; Li, L.; Zimmerman, P

    Rask, A. E.; Li, L.; Zimmerman, P. M. Kohn–Sham Density in a Slater Orbital Basis Set. J. Phys. Chem. A 2024, 128, 3194--3204

    2024

  32. [40]

    Levy, M.; Perdew, J. P. Hellmann-Feynman, virial, and scaling requisites for the exact universal density functionals. Shape of the correlation potential and diamagnetic susceptibility for atoms. Phys. Rev. A 1985, 32, 2010--2021

    1985

  33. [41]

    P.; Chulkov, S

    Gaiduk, A. P.; Chulkov, S. K.; Staroverov, V. N. Reconstruction of Density Functionals from Kohn Sham Potentials by Integration along Density Scaling Paths. J. Chem. Theory Comput. 2009, 5, 699--707

    2009

  34. [42]

    One-electron self-interaction and the asymptotics of the Kohn–Sham potential: an impaired relation

    Schmidt, T.; Kraisler, E.; Kronik, L.; Kümmel, S. One-electron self-interaction and the asymptotics of the Kohn–Sham potential: an impaired relation. Phys. Chem. Chem. Phys. 2014, 16, 14357--14367

    2014

  35. [43]

    W.; Yang, W

    Wu, Q.; Ayers, P. W.; Yang, W. Density-functional theory calculations with correct long-range potentials. J. Chem. Phys. 2003, 119, 2978--2990

    2003

  36. [44]

    Asymptotic Behavior of the Exchange-Correlation Energy Density and the Kohn-Sham Potential in Density Functional Theory: Exact Results and Strategy for Approximations

    Kraisler, E. Asymptotic Behavior of the Exchange-Correlation Energy Density and the Kohn-Sham Potential in Density Functional Theory: Exact Results and Strategy for Approximations. Isr. J. Chem. 2020, 60, 805--822

    2020

  37. [45]

    van Leeuwen, R.; Baerends, E. J. Energy expressions in density-functional theory using line integrals. Phys. Rev. A 1995, 51, 170--178

    1995

  38. [46]

    D.; Staroverov, V

    Elkind, P. D.; Staroverov, V. N. Energy expressions for Kohn–Sham potentials and their relation to the Slater–Janak theorem . J. Chem. Phys. 2012, 136, 124115

    2012

  39. [47]

    P.; Parr, R

    Perdew, J. P.; Parr, R. G.; Levy, M.; Balduz, J. L. Density-Functional Theory for Fractional Particle Number: Derivative Discontinuities of the Energy. Phys. Rev. Lett. 1982, 49, 1691--1694

    1982

  40. [48]

    J.; Mori-Sánchez, P

    Yang, W.; Cohen, A. J.; Mori-Sánchez, P. Derivative discontinuity, bandgap and lowest unoccupied molecular orbital in density functional theory. J. Chem. Phys. 2012, 136, 204111

    2012

  41. [49]

    A.; Tubman, N

    Holmes, A. A.; Tubman, N. M.; Umrigar, C. J. Heat-bath configuration interaction: An efficient selected configuration interaction algorithm inspired by heat-bath sampling. J. Chem. Theory Comput. 2016, 12, 3674--3680

    2016

  42. [50]

    A.; Jeanmairet, G.; Alavi, A.; Umrigar, C

    Sharma, S.; Holmes, A. A.; Jeanmairet, G.; Alavi, A.; Umrigar, C. J. Semistochastic heat-bath configuration interaction method: Selected configuration interaction with semistochastic perturbation theory. J. Chem. Theory Comput. 2017, 13, 1595--1604

    2017

  43. [51]

    A.; Sharma, S.; Umrigar, C

    Li, J.; Otten, M.; Holmes, A. A.; Sharma, S.; Umrigar, C. J. Fast semistochastic heat-bath configuration interaction. J. Chem. Phys. 2018, 149, 214110

    2018

  44. [52]

    A.; Zimmerman, P

    Dang, D.-K.; Kammeraad, J. A.; Zimmerman, P. M. Advances in parallel heat bath configuration interaction. J. Phys. Chem. A 2023, 127, 400--411

    2023

  45. [53]

    D.; Holmes, A

    Chien, A. D.; Holmes, A. A.; Otten, M.; Umrigar, C. J.; Sharma, S.; Zimmerman, P. M. Excited states of methylene, polyenes, and ozone from heat-bath configuration interaction. J. Phys. Chem. A 2018, 122, 2714--2722

    2018

  46. [54]

    V.; Baerends, E

    Lenthe, E. V.; Baerends, E. J. Optimized Slater-type basis sets for the elements 1--118. J. Comput. Chem. 2003, 24, 1142--1156

    2003

  47. [55]

    W.; Zimmerman, P

    Dang, D.-K.; Wilson, L. W.; Zimmerman, P. M. The numerical evaluation of Slater integrals on graphics processing units. J. Comput. Chem. 2022, 43, 1680--1689

    2022

  48. [56]

    On the eigenfunctions of many-particle systems in quantum mechanics

    Kato, T. On the eigenfunctions of many-particle systems in quantum mechanics. Comm. Pure Appl. Math. 1957, 10, 151--177

    1957

  49. [57]

    Handy, N. C. The molecular physics lecture 2004: (i) Density functional theory, (ii) Quantum Monte Carlo. Mol. Phys. 2004, 102, 2399--2409

    2004

  50. [58]

    E.; Knowles, P

    Mura, M. E.; Knowles, P. J. Improved radial grids for quadrature in molecular density‐functional calculations. J. Chem. Phys. 1996, 104, 9848--9858

    1996

  51. [59]

    Quadratures on a sphere

    Lebedev, V. Quadratures on a sphere. USSR Comput. Math. Math. Phys. 1976, 16, 10--24

    1976

  52. [60]

    OpenACC Organization, 2022

    The OpenACC Application Programming Interface, Version 3.3. OpenACC Organization, 2022

    2022

  53. [61]

    J.; Handy, N

    Cohen, A. J.; Handy, N. C. Density functional generalized gradient calculations using Slater basis sets. J. Chem. Phys. 2002, 117, 1470--1478

    2002

  54. [62]

    P.; Burke, K.; Ernzerhof, M

    Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865--3868

    1996

  55. [63]

    J.; Marques, M

    Lehtola, S.; Steigemann, C.; Oliveira, M. J.; Marques, M. A. Recent developments in libxc — A comprehensive library of functionals for density functional theory. SoftwareX 2018, 7, 1--5

    2018

  56. [64]

    E.; Dunning, J., Thom H

    Woon, D. E.; Dunning, J., Thom H. Gaussian basis sets for use in correlated molecular calculations. V. Core‐valence basis sets for boron through neon. J. Chem. Phys. 1995, 103, 4572--4585

    1995

  57. [65]

    Asymptotic near-nucleus structure of the electron-interaction potential in local effective potential theories

    Qian, Z.; Sahni, V. Asymptotic near-nucleus structure of the electron-interaction potential in local effective potential theories. Phys. Rev. A 2007, 75, 032517

    2007

  58. [66]

    Exchange and correlation near the nucleus in density functional theory

    Qian, Z. Exchange and correlation near the nucleus in density functional theory. Phys. Rev. B 2007, 75, 193104

    2007

  59. [67]

    C.; Staroverov, V

    Moore, C. C.; Staroverov, V. N. Jump discontinuities of finite-basis-set exchange–correlation potentials at atomic nuclei. J. Chem. Phys. 2024, 161, 124106

    2024

  60. [68]

    C.; Staroverov, V

    Moore, C. C.; Staroverov, V. N. Are Exact Exchange–Correlation Potentials Continuous at Atomic Nuclei in Molecules? J. Chem. Theory Comput. 2024, 20, 8934--8939

    2024

  61. [69]

    N.; Scuseria, G

    Staroverov, V. N.; Scuseria, G. E.; Tao, J.; Perdew, J. P. Tests of a ladder of density functionals for bulk solids and surfaces. Phys. Rev. B 2004, 69, 075102

    2004

  62. [70]

    P.; Ryabinkin, I

    Gaiduk, A. P.; Ryabinkin, I. G.; Staroverov, V. N. Removal of Basis-Set Artifacts in Kohn–Sham Potentials Recovered from Electron Densities. J. Chem. Theory Comput. 2013, 9, 3959--3964

    2013

  63. [71]

    J.; Mori-Sánchez, P.; Yang, W

    Cohen, A. J.; Mori-Sánchez, P.; Yang, W. Fractional spins and static correlation error in density functional theory. J. Chem. Phys. 2008, 129, 121104

    2008

  64. [72]

    The Kohn--Sham gap, the fundamental gap, and the optical gap: The physical meaning of occupied and virtual Kohn--Sham orbital energies

    Baerends, E.; Gritsenko, O.; van Meer, R. The Kohn--Sham gap, the fundamental gap, and the optical gap: The physical meaning of occupied and virtual Kohn--Sham orbital energies. Phys. Chem. Chem. Phys. 2013, 15, 16408--16425

    2013

  65. [73]

    Hellgren, M.; Gross, E. K. U. Discontinuities of the exchange-correlation kernel and charge-transfer excitations in time-dependent density-functional theory. Phys. Rev. A 2012, 85, 022514

    2012

  66. [74]

    Delocalization Errors in Density Functional Theory Are Essentially Quadratic in Fractional Occupation Number

    Hait, D.; Head-Gordon, M. Delocalization Errors in Density Functional Theory Are Essentially Quadratic in Fractional Occupation Number. J. Phys. Chem. Lett. 2018, 9, 6280--6288

    2018

  67. [75]

    P.; Staroverov, V

    Gaiduk, A. P.; Staroverov, V. N. Construction of integrable model Kohn-Sham potentials by analysis of the structure of functional derivatives. Phys. Rev. A 2011, 83, 012509

    2011

  68. [76]

    M.; Gavini, V

    Kanungo, B.; Zimmerman, P. M.; Gavini, V. A Comparison of Exact and Model Exchange–Correlation Potentials for Molecules. J. Phys. Chem. Lett. 2021, 12, 12012--12019

    2021

  69. [77]

    J.; Yang, W

    Mori-S\'anchez, P.; Cohen, A. J.; Yang, W. Localization and Delocalization Errors in Density Functional Theory and Implications for Band-Gap Prediction. Phys. Rev. Lett. 2008, 100, 146401

    2008

  70. [78]

    Chan, G. K.-L. A fresh look at ensembles: Derivative discontinuities in density functional theory. J. Chem. Phys. 1999, 110, 4710--4723

    1999

  71. [79]

    Baerends, E. J. Chemical potential , derivative discontinuity , fractional electrons , jump of the Kohn–Sham potential , atoms as thermodynamic open systems , and other (mis)conceptions of the density functional theory of electrons in molecules. Phys. Chem. Chem. Phys. 2022, 2...

    2022

  72. [80]

    J.; Mori-Sánchez, P.; Yang, W

    Cohen, A. J.; Mori-Sánchez, P.; Yang, W. Challenges for Density Functional Theory. Chem. Rev. 2012, 112, 289--320

    2012

  73. [81]

    step basis

    Hodgson, M. J. P.; Kraisler, E.; Schild, A.; Gross, E. K. U. How Interatomic Steps in the Exact Kohn–Sham Potential Relate to Derivative Discontinuities of the Energy. J. Phys. Chem. Lett. 2017, 8, 5974--5980 mcitethebibliography main.tex000066400000000000000000001612131475365...

    2010

  74. [82]

    4A v +a,Ďh -&Ăw,QcA]AX0vE P -&.ԡ<sĝs,;J A @DDDDDNK J o

    Kramida, A.; Ralchenko, Y.; Reader, J.; Team, N. A. NIST Atomic Spectra Database (version 5.12) . https://physics.nist.gov/asd, 2024; National Institute of Standards and Technology, Gaithersburg, MD mcitethebibliography si.bib0000664000000000000000000000745014753654340010670 0...

    2024

Pith tools

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