REVIEW 2 major objections 6 minor 54 references
Minimizing propagated density errors of atomic core-electron for simultaneously accurate bandgaps and lattice constants in closed-shell Copper semiconductors
T0 review · 2 major / 6 minor · reviewed 2026-07-09 · glm-5.2
Pith's one-line read Copper semiconductor errors traced to core-electron density in pseudopotentials
desk verdict Practical fix for Cu semiconductor bandgaps and lattice constants via HF core pseudopotentials; causal mechanism plausible but indirectly validated read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central mechanism is the quantum screening effect of core electrons on nuclear attraction. In pseudopotential construction, the core electrons' density determines how much of the nuclear Coulomb attraction is effectively cancelled by wavefunction orthogonality (quantum screening), leaving a residual pseudopotential for the valence electrons. When (semi-)local functionals over-delocalize the core density, this screening is overestimated, weakening the effective nuclear attraction and pushing semi-core 3d electrons to higher energies. The mHF pseudopotential (HF exchange + PBE correlation) removes self-interaction from the core, restoring correct screening, lowering the 3d level, and deep-
What would settle it
Construct an independent, high-level all-electron reference density for Cu-containing solids (e.g., via quantum Monte Carlo or diffusion Monte Carlo) and show that the mHF-core density does not agree better with it than the LDA-core density. Alternatively, find a class of Cu semiconductors where mHF@LDA systematically fails on bandgaps or lattice constants, which would contradict the claim of class-wide correction.
Extended reading notes
Core claim
The primary cause of bandgap and lattice-constant errors in Cu-containing closed-shell semiconductors is not a valence-electron functional deficiency but an error in the atomic core-electron density that is encoded into the pseudopotential and propagated into the solid. Self-interaction in (semi-)local functionals over-delocalizes Cu core electrons, enhancing their screening of the nuclear charge and raising the 3d energy level. Correcting the core density via a modified Hartree-Fock pseudopotential, while retaining a local functional for valence electrons, fixes both bandgaps and lattice constants simultaneously across an entire material class.
Load-bearing premise
The entire improvement rests on the claim that the modified Hartree-Fock pseudopotential faithfully encodes a more accurate Cu core-electron density. This is validated indirectly—through better atomic d-s splitting and better solid-state bandgaps and lattices—but no independent high-level benchmark of the solid-state density itself is provided, leaving open the possibility that mHF introduces a compensating error that happens to improve the target properties rather than truly
Editorial extensions
If this is right
- If the core-density propagation mechanism is general, the same mHF-core-plus-local-valence strategy should improve bandgaps and lattice constants for other transition-metal semiconductors where d-band misplacement is the dominant error source.
- The finding reframes pseudopotentials from mere computational accelerators into physical correction tools: the choice of functional used to construct the pseudo-core directly controls solid-state accuracy, opening a design space for pseudopotential engineering.
- The d-s energy splitting in the isolated atom (4.35 eV by mHF vs. 0.74 eV by LDA vs. 5.04 eV experiment) serves as a cheap diagnostic for whether a given pseudopotential will produce correct bandgaps in the resulting solid.
- The approach suggests that an ideal all-electron exchange-correlation functional would use position-dependent exact-exchange mixing: full HF in the core region, transitioning to local-density behavior in the valence region.
- Corrected Cu-3d energy positions would propagate into improved predictions for d-band center, catalytic activity, defect levels, and exciton binding energies in Cu-based materials.
Reading between the lines
- If the core-density error mechanism is as dominant as claimed for Cu, one would expect a monotonic relationship between the atomic d-s splitting error and the solid-state bandgap error across different transition metals, which could be tested systematically.
- The partial-core variant (mHF-17, including only 1s2s2p in the pseudo-core) yielding intermediate corrections suggests a tunable dial: the fraction of core electrons treated with exact exchange could be optimized per element, though this risks fitting rather than physics.
- The argument that LDA is more appropriate for valence electrons because they resemble a free-electron gas may break down for strongly correlated or Mott-insulating Cu compounds where valence electrons are themselves localized, potentially limiting the method's transferability.
- An independent all-electron coupled-cluster density benchmark for the Cu atom, compared directly against LDA, PBE, and mHF densities, would provide a stronger test of whether mHF is genuinely closer to the true density or merely produces a compensating error.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript investigates how errors in Cu core-electron density, encoded in pseudopotentials and propagated into solid-state DFT calculations, affect bandgaps and lattice constants in Cu-containing closed-shell semiconductors. The authors show that (semi-)local functionals (LDA/PBE) over-delocalize the Cu core density due to self-interaction, leading to over-screening of the nuclear attraction and too-high Cu 3d energy levels. They propose using modified Hartree-Fock (mHF: HF exchange + PBE correlation) pseudopotentials for the Cu core while retaining (semi-)local functionals for valence electrons. Across 50+ Cu-containing semiconductors, this mHF@LDA approach yields MRE of ~29% for bandgaps (vs ~80% for LDA@LDA) and ~0.9% for lattice constants, with no erroneous metals. The physical mechanism is traced through atomic d-s splitting (0.74 eV for LDA vs 4.35 eV for mHF vs 5.04 eV experimental), density difference analysis, and pseudo-core boundary tests.
Significance. The paper addresses a practically important and under-studied problem: how atomic density errors propagate through pseudopotentials into solid-state DFT. The statistical evidence across 50+ compounds is a clear strength, as is the simultaneous benchmarking of two independent physical quantities (bandgap and lattice constant). The real-space partitioning philosophy—different functionals for core vs. valence regions, mediated through pseudopotentials—is conceptually clean and immediately implementable. The pseudo-core boundary test (Table I, mHF-17 vs mHF-11) provides a falsifiable check on the mechanism rather than a fitted parameter. The all-electron LDA comparison (Section III.D) usefully rules out pseudopotential construction artifacts as the source of improvement.
major comments (2)
- [Section I and Section III.B] The central causal claim—that LDA/PBE self-interaction produces an over-delocalized core density, which over-screens nuclear attraction and raises the 3d level—is validated only indirectly. The atomic d-s splitting (4.35 eV for mHF vs 5.04 eV experimental, vs 0.74 eV for LDA) is a single scalar observable; it cannot distinguish 'mHF core density is more accurate' from 'mHF systematically lowers the 3d level for reasons that include but are not limited to corrected core density.' The paper itself notes (Section I) that high-level wavefunction methods provide 'nearly exact electron density distribution' for atoms, yet no comparison of the full radial density ρ(r) against a CCSD(T) or QMC benchmark is performed. Such a comparison would directly confirm whether the mHF core density is closer to the true density, transforming the current plausible narrative into a demonstrated causal chain. I
- [Section III.B, Eq. (1)] The mHF functional (HF exchange + PBE correlation) differs from LDA/PBE in multiple respects simultaneously—exchange treatment, correlation treatment, and the resulting screening. Attributing the improvement solely to 'corrected core-electron density' is therefore underdetermined. Can the authors provide any disentangling test? For example, comparing against a core density from PBE0 or other hybrid functionals, or decomposing the density change into exchange-driven vs. correlation-driven components, would help isolate the mechanism. The pseudo-core boundary test (Table I) partially addresses this by showing that including 3s3p in the mHF core matters, which is consistent with the self-interaction mechanism, but it does not rule out alternative explanations for why exact exchange improves things.
minor comments (6)
- [Figure 1(a)] The caption states that for erroneously metallic systems, 'the gap is defined between Cu-4s and -3d states at Γ point,' but it is unclear how these negative gaps are treated in the MRE calculation. The caption mentions 'marginal statistical impact' but the precise treatment (excluded? included as negative?) should be stated explicitly.
- [Figure 1(c)] The histogram binning convention ('histograms are at (0, ±10%, ±30%, etc.); each covers a range of ±10%') is somewhat confusing. A clearer statement of the bin edges would help readers interpret the probability distributions.
- [Section III.C] The term 'quantum screening' is introduced to describe the pseudopotential's effective repulsion from orthogonality requirements. This terminology may cause confusion with classical electrostatic screening. A brief clarifying remark distinguishing the two would improve readability.
- [Table I] The notation 'mHF-17' and 'mHF-11' refers to the number of valence electrons, but this is not immediately obvious. A footnote or parenthetical clarifying the convention would help.
- [Section III.D] The all-electron LDA calculations using the Elk code are mentioned for Cu2S and CuCl, but the computational parameters (k-mesh, basis set convergence) for these all-electron calculations are not provided. Including these details would aid reproducibility.
- [Reference 14] Reference 14 (Ye et al., J. Chem. Theory Comput. 2025) appears to be by some of the present authors and may contain related methodology. The relationship to the present work should be clarified to ensure novelty is transparent.
Circularity Check
No significant circularity: the mHF pseudopotential is constructed from an independent functional form (HF exchange + PBE correlation), not fitted to the target solid-state properties; validation is against external experimental benchmarks.
full rationale
The paper's central claim is that using mHF pseudopotentials for the Cu core eliminates propagated density errors, improving bandgaps and lattice constants. The derivation chain is: (1) LDA/PBE self-interaction over-delocalizes core density → (2) over-screening of nuclear attraction → (3) too-high Cu 3d energy → (4) bandgap underestimation. Step (1) is supported by atomic d-s splitting (0.74 eV LDA vs. 4.35 eV mHF vs. 5.04 eV experiment), an external benchmark. Step (4) is validated against experimental bandgaps and lattice constants for 50+ compounds, also external. The mHF functional (Eq. 1: E_xc^mHF = E_x^HF + E_c^PBE) is a fixed functional form, not fitted to solid-state targets. The choice of pseudo-core electrons (1s2s2p3s3p, 11-electron) is tested against alternatives (Table I, mHF-17 vs. mHF-11), showing sensitivity but not circularity — the 11-electron choice is justified by the physical argument that all inner-shell electrons with self-interaction error should be included, and the results are compared against experiment, not against a fitted target. The paper cites prior work by some of the same authors (Refs. 27, 28, 36) for the mHF construction and pseudopotential methodology, but these citations provide methodological background, not the load-bearing validation. The central claim's support comes from external experimental data, not from self-citation. The reader's concern about indirect validation (d-s splitting as a single scalar, no CCSD(T) density comparison) is a correctness risk, not a circularity issue — the paper does not define its inputs in terms of its outputs. No step in the derivation chain reduces to its own inputs by construction.
Assumptions & free parameters
free parameters (2)
- Pseudo core boundary (11-electron: 1s2s2p3s3p core) =
11 electrons in core
- mHF functional form (HF exchange + PBE correlation) =
E_xc = E_x^HF + E_c^PBE
assumptions (3)
- domain assumption Core electrons in solids retain their atomic characteristics and can be treated as fixed background charges encoded in pseudopotentials.
- domain assumption Simultaneous accuracy in bandgap and lattice constants across a class of materials implies correct electron density.
- domain assumption The mHF functional (HF exchange + PBE correlation) provides a more accurate description of localized core-electron density than LDA or PBE.
Cite this review
Pith. "Pith review of Minimizing propagated density errors of atomic core-electron for simultaneously accurate bandgaps and lattice constants in closed-shell Copper semiconductors." pith.science (2026). https://pith.science/paper/GBBLI54D
@misc{pith2026260707037,
author = {Pith},
title = {Pith review of: Minimizing propagated density errors of atomic core-electron for simultaneously accurate bandgaps and lattice constants in closed-shell Copper semiconductors},
year = {2026},
howpublished = {\url{https://pith.science/paper/GBBLI54D}},
note = {Machine review of arXiv:2607.07037}
}
abstract
Density functional theory struggles to accurately determine electron density of atoms, whose error is inevitably encoded into the pseudopotential and propagated into solid-state calculations. However, little is known about how this affects accuracy nor how to remedy it. In this work, through a systematic study of the effect of Cu atomic density on bandgap and lattice constants of over 50 Cu-containing simple closed-shell semiconductors, we find that core-electron density can drastically affect nuclear attraction to valence electrons and subsequent charge distribution and energy position of Cu 3$d$ electrons. The error can be eliminated at its source by employing modified Hartree-Fock pseudopotentials for Cu core while retaining (semi-)local functionals for valence electrons. This real-space partitioning approach leads to simultaneous high-accuracy in bandgap and lattice constants across the entire material class.
Figures
Reference graph
Works this paper leans on
-
[1]
Huang, Bing and von Rudorff, Guido Falk and von Lilienfeld, O. Anatole , month = jul, year =. The central role of density functional theory in the. Science , publisher =. doi:10.1126/science.abn3445 , abstract =
-
[2]
Density functional theory:. Rev. Mod. Phys. , author =. 2015 , pages =. doi:10.1103/RevModPhys.87.897 , language =
-
[3]
A conundrum for density functional theory , volume =
Hammes-Schiffer, Sharon , month = jan, year =. A conundrum for density functional theory , volume =. Science , publisher =. doi:10.1126/science.aal3442 , number =
-
[4]
Inhomogeneous electron gas.Phys
Hohenberg, P. and Kohn, W. , month = nov, year =. Inhomogeneous. Phys. Rev. , publisher =. doi:10.1103/PhysRev.136.B864 , abstract =
-
[5]
Kohn, W. and Sham, L. J. , month = nov, year =. Self-. Phys. Rev. , publisher =. doi:10.1103/PhysRev.140.A1133 , abstract =
-
[6]
The top 100 papers , volume =. Nature News , author =. 2014 , note =. doi:10.1038/514550a , abstract =
-
[7]
Density functional theory is straying from the path toward the exact functional , volume =. Science , author =. 2017 , pages =. doi:10.1126/science.aah5975 , abstract =
-
[8]
Quantum chemistry:. Nat. Rev. Chem. , author =. 2017 , pages =. doi:10.1038/s41570-017-0019 , language =
Show all 54 references
-
[9]
WIREs Comput Mol Sci , author =
The uniform electron gas , volume =. WIREs Comput Mol Sci , author =. 2016 , pages =. doi:10.1002/wcms.1257 , abstract =
2016 doi
-
[10]
Comment on "
Zhang, Yingkai and Yang, Weitao , month = jan, year =. Comment on ". Phys. Rev. Lett. , publisher =. doi:10.1103/PhysRevLett.80.890 , abstract =
-
[11]
A. J. Chem. Phys. , author =. 1935 , pages =. doi:10.1063/1.1749559 , number =
1935 doi
-
[12]
Fifty years of pseudopotentials , volume =. Int. J. Quantum Chem. , author =. 1983 , pages =. doi:10.1002/qua.560240862 , abstract =
1983 doi
-
[13]
Pseudopotentials, an. J. Chem. Theory Comput. , author =. 2025 , pages =. doi:10.1021/acs.jctc.5c01373 , abstract =
2025 doi
-
[14]
J Chem Phys , author =
The densities produced by the density functional theory: comparison to full configuration interaction , volume =. J Chem Phys , author =. 2008 , pages =. doi:10.1063/1.2821123 , abstract =
2008 doi
-
[15]
and Levy, Mel , month = nov, year =
Perdew, John P. and Levy, Mel , month = nov, year =. Physical. Phys. Rev. Lett. , publisher =. doi:10.1103/PhysRevLett.51.1884 , abstract =
-
[16]
Sham, L. J. and Schlüter, M. , month = nov, year =. Density-. Phys. Rev. Lett. , publisher =. doi:10.1103/PhysRevLett.51.1888 , abstract =
-
[17]
Electronic band structure of cuprous delafossite
Zhang, Min-Ye and Jiang, Hong , month = mar, year =. Electronic band structure of cuprous delafossite. Phys. Rev. B , publisher =. doi:10.1103/PhysRevB.111.125204 , abstract =
-
[18]
J PHYS-CONDENS MAT , author =
Density functional theory study of the electronic structure of fluorite. J PHYS-CONDENS MAT , author =. 2013 , pages =. doi:10.1088/0953-8984/25/12/125503 , abstract =
2013 doi
-
[19]
Structural, electronic, and optical properties of. Appl. Phys. Lett. , author =. 2013 , pages =. doi:10.1063/1.4824770 , abstract =
2013 doi
-
[20]
Lukashev, Pavel and Lambrecht, Walter R. L. and Kotani, Takao and van Schilfgaarde, Mark , month = nov, year =. Electronic and crystal structure of. Phys. Rev. B , publisher =. doi:10.1103/PhysRevB.76.195202 , abstract =
-
[21]
Quasiparticle band structures of
Gao, Weiwei and Xia, Weiyi and Wu, Yabei and Ren, Wei and Gao, Xiang and Zhang, Peihong , month = jul, year =. Quasiparticle band structures of. Phys. Rev. B , publisher =. doi:10.1103/PhysRevB.98.045108 , abstract =
-
[22]
Uncertainty of. J. Chem. Theory Comput. , author =. 2024 , pages =. doi:10.1021/acs.jctc.4c01036 , abstract =
2024 doi
-
[23]
Putting error bars on density functional theory , volume =. Sci. Rep. , author =. 2024 , pages =. doi:10.1038/s41598-024-69194-w , abstract =
2024 doi
-
[24]
Calculation of the lattice constant of solids with semilocal functionals , volume =
Haas, Philipp and Tran, Fabien and Blaha, Peter , month = feb, year =. Calculation of the lattice constant of solids with semilocal functionals , volume =. Phys. Rev. B , publisher =. doi:10.1103/PhysRevB.79.085104 , abstract =
-
[25]
and Meek, Terry L
Mann, Joseph B. and Meek, Terry L. and Knight, Eugene T. and Capitani, Joseph F. and Allen, Leland C. , month = may, year =. Configuration. J. Am. Chem. Soc. , publisher =. doi:10.1021/ja9928677 , abstract =
-
[26]
A way to identify whether a. J. Chem. Phys. , author =. 2024 , pages =. doi:10.1063/5.0196706 , abstract =
2024 doi
-
[27]
Accuracy trade-off between one-electron and excitonic spectra of cuprous halides in first-principles calculations , volume =. J. Chem. Phys. , author =. 2021 , pages =. doi:10.1063/5.0043999 , abstract =
2021 doi
-
[28]
Effect of. Phys. Chem. Chem. Phys. , author =. 2018 , pages =. doi:10.1039/C8CP00990B , abstract =
2018 doi
-
[29]
Advanced capabilities for materials modelling with
Giannozzi, P and Andreussi, O and Brumme, T and Bunau, O and Buongiorno Nardelli, M and Calandra, M and Car, R and Cavazzoni, C and Ceresoli, D and Cococcioni, M and Colonna, N and Carnimeo, I and Dal Corso, A and de Gironcoli, S and Delugas, P and DiStasio, R A and Ferretti, ...
-
[30]
and Rappe, Andrew M
Yang, Jing and Tan, Liang Z. and Rappe, Andrew M. , month = feb, year =. Hybrid functional pseudopotentials , volume =. Phys. Rev. B , publisher =. doi:10.1103/PhysRevB.97.085130 , abstract =
-
[31]
Understanding the origin of bandgap problem in transition and post-transition metal oxides , volume =. J. Chem. Phys. , author =. 2019 , pages =. doi:10.1063/1.5111188 , abstract =
2019 doi
-
[32]
Ceperley, D. M. and Alder, B. J. , month = aug, year =. Ground. Phys. Rev. Lett. , publisher =. doi:10.1103/PhysRevLett.45.566 , abstract =
-
[33]
Perdew, J. P. and Zunger, Alex , month = may, year =. Self-interaction correction to density-functional approximations for many-electron systems , volume =. Phys. Rev. B , publisher =. doi:10.1103/PhysRevB.23.5048 , abstract =
-
[34]
and Burke, Kieron and Ernzerhof, Matthias , month = oct, year =
Perdew, John P. and Burke, Kieron and Ernzerhof, Matthias , month = oct, year =. Generalized. Phys. Rev. Lett. , publisher =. doi:10.1103/PhysRevLett.77.3865 , abstract =
-
[35]
, month = oct, year =
Kohn, Walter and Mattsson, Ann E. , month = oct, year =. Edge. Phys. Rev. Lett. , publisher =. doi:10.1103/PhysRevLett.81.3487 , abstract =
-
[36]
An energy functional for surfaces , volume =. J. Chem. Phys. , author =. 2001 , pages =. doi:10.1063/1.1396649 , abstract =
2001 doi
-
[37]
Perdew, J. P. and Burke, K. and Ernzerhof, M. , month = jan, year =. Perdew,. Phys. Rev. Lett. , publisher =. doi:10.1103/PhysRevLett.80.891 , abstract =
-
[38]
OPIUM Documentation , howpublished =
-
[39]
Advancing. J. Chem. Phys. , author =. 2024 , pages =. doi:10.1063/5.0232711 , abstract =
2024 doi
-
[40]
Lambert, D. S. and O’Regan, D. D. , month = mar, year =. Use of. Phys. Rev. Res. , publisher =. doi:10.1103/PhysRevResearch.5.013160 , abstract =
-
[41]
Band gaps of crystalline solids from. Proc. Natl. Acad. Sci. U.S.A. , author =. 2021 , pages =. doi:10.1073/pnas.2104556118 , abstract =
2021 doi
-
[42]
Interpretations of ground-state symmetry breaking and strong correlation in wavefunction and density functional theories , volume =. Proc. Natl. Acad. Sci. U.S.A. , author =. 2021 , pages =. doi:10.1073/pnas.2017850118 , abstract =
2021 doi
-
[43]
Inapplicability of exact constraints and a minimal two-parameter generalization to the. Phys. Rev. B , author =. 2016 , note =. doi:10.1103/PhysRevB.94.220104 , abstract =
2016 doi
-
[44]
The Journal of Chemical Physics , author =
Electronic structure of mononuclear. The Journal of Chemical Physics , author =. 2021 , pages =. doi:10.1063/5.0054439 , abstract =
2021 doi
-
[45]
Point defects in. Sci. Technol. Adv. Mater. , author =. 2011 , pages =. doi:10.1088/1468-6996/12/3/034302 , abstract =
2011 doi
-
[46]
and Lany, Stephan and Chang, Young Jun and Rotenberg, Eli and Zunger, Alex and Toney, Michael F
Lim, Linda Y. and Lany, Stephan and Chang, Young Jun and Rotenberg, Eli and Zunger, Alex and Toney, Michael F. , month = dec, year =. Angle-resolved photoemission and quasiparticle calculation of. Phys. Rev. B , publisher =. doi:10.1103/PhysRevB.86.235113 , abstract =
-
[47]
Hamann, D. R. and Schl. Norm-Conserving Pseudopotentials , journal =. 1979 , doi =
1979
-
[48]
Computer Physics Communications , volume =
Optimization algorithm for the generation of ONCV pseudopotentials , author =. Computer Physics Communications , volume =. 2015 , doi =
2015
-
[49]
, title =
Hellmann, H. , title =. The Journal of Chemical Physics , volume =. 1935 , doi =
1935
-
[50]
, title =
Cohen, Marvin L. , title =. International Journal of Quantum Chemistry , volume =. 1983 , doi =
1983
-
[51]
Melius, C. F. and Goddard, III, W. A. , title =. Physical Review A , volume =. 1974 , doi =
1974
-
[52]
Dewhurst, J. K. and Sharma, S. and Nordstr. The Elk Code , howpublished =
-
[53]
Booth, George H. and Gr. Towards an Exact Description of Electronic Wavefunctions in Real Solids , journal =. 2013 , doi =
2013
-
[54]
Natural Orbitals for Wave Function Based Correlated Calculations Using a Plane Wave Basis Set , journal =
Gr. Natural Orbitals for Wave Function Based Correlated Calculations Using a Plane Wave Basis Set , journal =. 2011 , doi =
2011
Reviewed July 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.