Recognition: no theorem link
Systematically Improvable Numerical Atomic Orbital Basis Using Contracted Truncated Spherical Waves
Pith reviewed 2026-05-15 11:33 UTC · model grok-4.3
The pith
Contracting truncated spherical waves produces numerical atomic orbital bases that are systematically improvable and transferable across molecules and solids.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Contracting truncated spherical waves to minimize the kinetic trace in the residual space generates NAO bases that inherit systematic improvability from prior methods while achieving superior transferability, because the spherical-wave expansion connects reference states to the basis without spurious long-range couplings and supports inclusion of extensive unoccupied states.
What carries the argument
The contraction step that minimizes the trace of the kinetic operator in the residual space when expanding in truncated spherical waves, which directly produces the numerical atomic orbital functions.
If this is right
- Total energies, bond lengths, and atomization energies of molecules reach satisfactory accuracy with modest basis sizes.
- Lattice constants, cohesive energies, and band gaps of bulk solids converge reliably using the same NAO construction.
- Energy-level alignments between occupied and unoccupied states improve when unoccupied reference states are included in the contraction.
- The absence of periodic-image interactions allows the same basis to be applied without adjustment to both isolated molecules and extended systems.
Where Pith is reading between the lines
- The approach may reduce the number of basis functions needed for large-scale simulations while preserving accuracy, lowering computational cost for systems with hundreds of atoms.
- Because the construction avoids plane-wave periodicity artifacts, it could be especially useful for surfaces, clusters, or defects where artificial image interactions would otherwise distort results.
- The same contraction principle might be tested in other orbital-based methods such as Hartree-Fock or post-DFT correlation techniques to check whether the kinetic-trace minimum remains effective.
Load-bearing premise
Minimizing the kinetic trace through TSW contraction yields bases whose accuracy and improvability are not compromised by truncation or contraction artifacts.
What would settle it
A calculation in which the new NAO basis fails to reach the same converged value as a large reference plane-wave basis for total energy or band gap of a chosen molecule or crystal would show the method does not deliver the claimed systematic improvability.
Figures
read the original abstract
To solve the Kohn-Sham equation within the framework of density functional theory, we develop a scheme to construct numerical atomic orbital (NAO) basis sets by contracting truncated spherical waves (TSWs). The contraction minimizes the trace of the kinetic operator in the residual space, generalizing the spillage minimizing scheme [M. Chen et al., J. Phys. Condens. Matter 22, 445501 (2010); P. Lin et al., Phys. Rev. B 103, 235131 (2021)]. In addition to the systematic improvability inherited from previous schemes, the use of TSW instead of plane waves as the expansion basis bridges reference states and NAOs more effectively, and eliminates spurious interactions between periodic images, thereby enabling better transferability through the inclusion of extensive reference states. Benchmarks demonstrate that the constructed NAO achieves satisfactory precision for various properties of both molecules and bulk systems, including total energy, bond length, atomization energy, lattice constant, cohesive energy, band gap, and energy-level alignment. By incorporating unoccupied states, the improved transferability in describing the conduction band is demonstrated to be effective and substantial.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a scheme for constructing numerical atomic orbital (NAO) basis sets by contracting truncated spherical waves (TSWs). The contraction is achieved by minimizing the trace of the kinetic operator in the residual space, generalizing prior spillage-minimization approaches. The TSW choice is motivated by eliminating periodic-image artifacts, which enables inclusion of more extensive reference states (including unoccupied orbitals) and is claimed to yield improved transferability. Benchmarks are reported for both molecular and periodic systems, asserting satisfactory accuracy in total energies, bond lengths, atomization and cohesive energies, lattice constants, band gaps, and energy-level alignments, with notable gains in conduction-band descriptions.
Significance. If the central construction and benchmarks hold, the work supplies a systematically improvable NAO basis with demonstrably better transferability than plane-wave contractions, particularly for unoccupied states. The independent kinetic-trace criterion and TSW truncation provide concrete technical advances that could reduce basis-set errors in large-scale DFT calculations for molecules and solids.
major comments (2)
- [Abstract] Abstract and results sections: the central claim that the constructed NAO 'achieves satisfactory precision' for total energy, bond length, atomization energy, lattice constant, cohesive energy, band gap, and energy-level alignment is not supported by any quantitative metrics, error bars, direct comparisons to reference plane-wave or other NAO results, or details on data selection/exclusion. This absence is load-bearing for the transferability and improvability assertions.
- [Method] Method section: while the kinetic-trace minimization is presented as generalizing the spillage schemes of Chen et al. and Lin et al., the manuscript does not supply the explicit functional form or demonstrate that the new objective is independent rather than reducing to a reparameterization of prior fitted quantities; this must be clarified to establish the claimed systematic improvability.
minor comments (1)
- [Abstract] The abstract would be strengthened by including at least one concrete numerical benchmark (e.g., mean absolute error in total energy or band gap relative to a reference) to illustrate the 'satisfactory precision' claim.
Simulated Author's Rebuttal
We thank the referee for the constructive comments and the opportunity to improve the manuscript. We address each major point below and will revise accordingly to strengthen the quantitative support and methodological clarity.
read point-by-point responses
-
Referee: [Abstract] Abstract and results sections: the central claim that the constructed NAO 'achieves satisfactory precision' for total energy, bond length, atomization energy, lattice constant, cohesive energy, band gap, and energy-level alignment is not supported by any quantitative metrics, error bars, direct comparisons to reference plane-wave or other NAO results, or details on data selection/exclusion. This absence is load-bearing for the transferability and improvability assertions.
Authors: We agree that explicit quantitative metrics, error bars, and direct comparisons are needed to substantiate the claims. In the revised manuscript, we will add tables in the results section reporting mean absolute deviations (MADs) and maximum errors for total energies, bond lengths, atomization/cohesive energies, lattice constants, band gaps, and energy-level alignments, benchmarked against high-cutoff plane-wave calculations and existing NAO sets (e.g., from Chen et al. and Lin et al.). We will include error bars from statistical sampling over the test sets, specify all molecular and periodic systems used (with selection criteria and exclusion details), and add figures comparing conduction-band descriptions with/without unoccupied reference states. These additions will directly support the transferability and systematic improvability assertions. revision: yes
-
Referee: [Method] Method section: while the kinetic-trace minimization is presented as generalizing the spillage schemes of Chen et al. and Lin et al., the manuscript does not supply the explicit functional form or demonstrate that the new objective is independent rather than reducing to a reparameterization of prior fitted quantities; this must be clarified to establish the claimed systematic improvability.
Authors: We will supply the explicit functional form in the revised method section. The objective minimizes Tr(P_res K P_res), where K is the kinetic-energy operator and P_res is the projector onto the residual space orthogonal to the contracted NAO basis (derived via the TSW expansion). This differs from spillage minimization, which targets the L2 norm of the wavefunction or density residual. We will add a derivation subsection proving independence: the kinetic-trace criterion weights higher-momentum components variationally, yielding better unoccupied-state accuracy without reducing to a reparameterized spillage functional (numerical tests will confirm distinct convergence behavior). This clarification will establish the systematic improvability. revision: yes
Circularity Check
No significant circularity: derivation introduces independent kinetic-trace objective and TSW basis
full rationale
The paper generalizes prior spillage-minimization schemes (Chen 2010, Lin 2021) by replacing the objective with minimization of the kinetic-operator trace in the residual space and by adopting truncated spherical waves instead of plane waves. This change is explicitly motivated and does not reduce the new NAO construction to a re-labeling or re-fitting of the cited inputs. Benchmarks on molecules and solids are presented as external validation rather than internal tautologies. Self-citation exists but is not load-bearing for the central claim of improved transferability via unoccupied states and elimination of periodic-image artifacts. No self-definitional loop, fitted-input-as-prediction, or ansatz-smuggling is exhibited in the stated equations or method.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Minimizing the trace of the kinetic operator in the residual space produces systematically improvable NAO contractions that generalize spillage minimization.
Reference graph
Works this paper leans on
-
[1]
Generalized Spillage In the CGH [31] scheme, the spillage function S= X nk ⟨ψPW nk |[1− ˆPk]|ψPW nk ⟩ = X nk [1− ˆPk] ψPW nk 2 (8) is minimized, where{ ψPW nk }denotes the reference states obtained from high-quality plane-wave calculations,nis the band index andklabels different sampling points in the Brillouin zone. In particular, the projection operator...
-
[2]
Reference Systems To ensure the transferability of NAOs, reference sys- tems should cover potential polarization and bonding in real chemical environments. Blum et al. [3] and Chen et al. [31] used homonuclear dimers with various bond lengths. Lin et al. [40] included additional trimers in the generation of triple-zeta level NAOs. Their NAOs all show stab...
-
[3]
The angular momentum components included cover the range from 0 tol max
Converging truncated spherical wave parameters With a given kinetic energy cutoff for spherical waves, the maximum angular momentuml max and a uniform real-space cutoff radiusr c together define a unique set of NSWs. The angular momentum components included cover the range from 0 tol max. By increasingr c and lmax of NSWs, systematic converging behaviors ...
-
[4]
Molecules We first performed benchmarks on the built NAOs in 11 diatomic molecular systems including Br 2, CO, Cl 2, F2, I 2, LiH, Li 2, N 2, Na 2, O 2 and S 2. Among these molecules, there are single (LiH, Li 2, Na 2, F 2, Cl 2, Br 2, I2), double (O 2, S 2) and triple (CO, N 2) bonds. In the following, if not specifically mentioned, the molecules O 2 and...
-
[5]
Bulks We also benchmark our NAOs in 26 bulk systems, in- cluding the classical covalent and ionic crystals, insula- tors and semiconductors, on the relative energies, struc- tural (lattice constant), thermodynamic (cohesive ener- gies), mechanical (bulk moduli) properties, and proper- ties that are indirectly related to electron transport, in- cluding ban...
- [6]
- [7]
-
[8]
V. Blum, R. Gehrke, F. Hanke, P. Havu, V. Havu, X. Ren, K. Reuter, and M. Scheffler, Computer Physics Communications180, 2175 (2009)
work page 2009
-
[9]
C. Fonseca Guerra, J. Snijders, G. t. Te Velde, and E. J. Baerends, Theoretical Chemistry Accounts99, 391 (1998)
work page 1998
-
[10]
E. Artacho, D. S´ anchez-Portal, P. Ordej´ on, A. Gar- cia, and J. M. Soler, physica status solidi (b)215, 809 (1999)
work page 1999
- [11]
-
[12]
J. Hutter, M. Iannuzzi, F. Schiffmann, and J. Vande- 18 □ H N □ P H|PH|P N 0 5 10E − Ef (eV) (a) PW pVTZ □ H N □ P H|PH|P N 0 5 10E − Ef (eV) (b) PW pVTZ (PW) □ H N □ P H|PH|P N 0 5 10E − Ef (eV) (c) PW pVTZ-2V □ H N □ P H|PH|P N 0 5 10E − Ef (eV) (d) PW pVTZ-2V (PW) FIG. A.1. Na body-centered-cubic band structures calculated with PW, pVTZ, and pVTZ-2V ba...
work page 2014
- [13]
- [14]
-
[15]
E. Apra, E. J. Bylaska, W. A. De Jong, N. Govind, K. Kowalski, T. P. Straatsma, M. Valiev, H. J. van Dam, Y. Alexeev, J. Anchell,et al., The Journal of chemical physics152(2020)
work page 2020
- [16]
- [17]
- [18]
-
[19]
M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scal- mani, V. Barone, G. A. Petersson, H. Nakatsuji, X. Li, M. Caricato, A. V. Marenich, J. Bloino, B. G. Janesko, R. Gomperts, B. Mennucci, H. P. Hratchian, J. V. Or- tiz, A. F. Izmaylov, J. L. Sonnenberg, D. Williams- Young, F. Ding, F. Lipparini, F. Egidi, J. G...
work page 2016
-
[20]
Y. Shao, Z. Gan, E. Epifanovsky, A. T. Gilbert, M. Wor- mit, J. Kussmann, A. W. Lange, A. Behn, J. Deng, X. Feng,et al., Molecular Physics113, 184 (2015)
work page 2015
- [21]
-
[22]
G. M. Barca, C. Bertoni, L. Carrington, D. Datta, N. De Silva, J. E. Deustua, D. G. Fedorov, J. R. Gour, A. O. Gunina, E. Guidez,et al., The Journal of chemical physics152(2020). 19 TABLE B.1. Occupation weighted energy levels differences (η-metrics) of NAOs with NSW and PW basis calculated on 11 molecule systems (in eV) Molecule pVDZ pVTZ − pVTZ pVQZ = N...
work page 2020
-
[23]
M. F. Guest, I. J. Bush, H. J. Van Dam, P. Sherwood, J. M. Thomas, J. H. Van Lenthe, R. W. Havenith, and J. Kendrick, Molecular Physics103, 719 (2005)
work page 2005
-
[24]
R. Ahlrichs, M. B¨ ar, M. H¨ aser, H. Horn, and C. K¨ olmel, Chemical Physics Letters162, 165 (1989)
work page 1989
-
[25]
J. M. Turney, A. C. Simmonett, R. M. Parrish, E. G. Hohenstein, F. A. Evangelista, J. T. Fermann, B. J. TABLE B.2. Delta values with respect to the PW calculated by NAOs for 26 bulk systems (in meV/atom) Bulk pVDZ pVTZ − pVTZ pVQZ = AlAs 1.261 0.551 0.027 0.718 AlP 1.240 0.521 0.172 0.863 AlSb 0.848 0.597 0.180 0.325 AlN 3.459 0.873 0.125 1.722 BN 0.978 0...
work page 2012
-
[26]
T. L. Beck, Reviews of Modern Physics72, 1041 (2000)
work page 2000
-
[27]
P. D. Haynes and M. C. Payne, Computer Physics Com- munications102, 17 (1997)
work page 1997
-
[28]
E. Hern´ andez, M. Gillan, and C. Goringe, Physical Re- 20 TABLE B.3. Occupation weighted band structure differences with PW basis (η-metrics) calculated on 26 bulk systems (in eV) Bulk pVDZ pVTZ − pVTZ pVQZ = AlAs 0.017 0.008 0.004 0.006 AlN 0.065 0.040 0.022 0.024 AlP 0.020 0.013 0.008 0.008 AlSb 0.007 0.003 0.002 0.002 BN 0.004 0.002 0.001 0.002 BP 0.0...
work page 1997
- [29]
-
[30]
L. Genovese, A. Neelov, S. Goedecker, T. Deutsch, S. A. Ghasemi, A. Willand, D. Caliste, O. Zilberberg, M. Rayson, A. Bergman,et al., The Journal of chemical physics129(2008)
work page 2008
-
[31]
O. F. Sankey and D. J. Niklewski, Physical Review B 40, 3979 (1989), publisher: American Physical Society
work page 1989
-
[32]
D. Porezag, T. Frauenheim, T. K¨ ohler, G. Seifert, and R. Kaschner, Physical Review B51, 12947 (1995), pub- lisher: American Physical Society
work page 1995
-
[33]
A. P. Horsfield, Physical Review B56, 6594 (1997), pub- lisher: American Physical Society
work page 1997
-
[34]
J. Junquera, ´O. Paz, D. S´ anchez-Portal, and E. Arta- cho, Physical Review B64, 235111 (2001), publisher: American Physical Society
work page 2001
-
[35]
Ozaki, Physical Review B67, 155108 (2003), pub- lisher: American Physical Society
T. Ozaki, Physical Review B67, 155108 (2003), pub- lisher: American Physical Society
work page 2003
-
[36]
M. Chen, G.-C. Guo, and L. He, Journal of Physics: Condensed Matter22, 445501 (2010)
work page 2010
-
[37]
J. M. Soler, E. Artacho, J. D. Gale, A. Garc´ ıa, J. Jun- quera, P. Ordej´ on, and D. S´ anchez-Portal, Journal of Physics: Condensed Matter14, 2745 (2002)
work page 2002
- [38]
-
[39]
T. Ozaki and H. Kino, Physical Review B69, 195113 (2004), publisher: American Physical Society
work page 2004
-
[40]
Dunning, Thom H., The Journal of Chemical Physics 90, 1007 (1989)
J. Dunning, Thom H., The Journal of Chemical Physics 90, 1007 (1989)
work page 1989
-
[41]
I. Y. Zhang, X. Ren, P. Rinke, V. Blum, and M. Schef- fler, New Journal of Physics15, 123033 (2013), pub- lisher: IOP Publishing
work page 2013
-
[42]
D. S´ anchez-Portal, E. Artacho, and J. M. Soler, Solid State Communications95, 685 (1995)
work page 1995
-
[43]
D. S´ anchez-Portal, E. Artacho, and J. M. Soler, Journal of Physics: Condensed Matter8, 3859 (1996)
work page 1996
-
[44]
S. D. Kenny, A. P. Horsfield, and H. Fujitani, Physical Review B62, 4899 (2000), publisher: American Physi- cal Society
work page 2000
-
[45]
P. Lin, X. Ren, and L. He, Physical Review B103, 235131 (2021), publisher: American Physical Society
work page 2021
-
[46]
A. T. T. Pham, A. Vora-ud, C. Chananonnawathorn, P. Muthitamongkol, S. Ruamruk, S. Limwichean, M. Horprathum, R. Syariati, S. Park, F. Ishii,et al., Vacuum , 114879 (2025)
work page 2025
-
[47]
S. Soni, K. Deshmukh, J. K. Saluja, R. Kumari, M. A. Pateria, and M. L. Verma, Journal of Materials Science: Materials in Electronics36, 1 (2025)
work page 2025
- [48]
-
[49]
F. V. Calderan, K. F. Andriani, P. Felicio-Sousa, G. A. Pinheiro, J. L. Da Silva, and M. G. Quiles, ACS Omega (2026)
work page 2026
-
[50]
I. Mandzhieva, F. Theiss, X. He, A. Ortmeier, A. Koirala, S. J. McBride, S. J. DeVience, M. S. Rosen, V. Blum, and T. Theis, Journal of Chemical Informa- tion and Modeling65, 7554 (2025)
work page 2025
-
[51]
O. Perevozchikov, V. Perevoztchikov, A. Khan, and A. Fedoseyev, inJournal of Physics: Conference Series, Vol. 3145 (IOP Publishing, 2025) p. 012006
work page 2025
-
[52]
H. Liu, W. He, H. Xu, X. Wang, and Q. An, Current Applied Physics54, 75 (2023)
work page 2023
-
[53]
C. L. Box, W. G. Stark, and R. J. Maurer, Electronic Structure5, 035005 (2023)
work page 2023
-
[54]
S. Kumar and R. Moudgil, Physical Chemistry Chemi- cal Physics24, 21283 (2022)
work page 2022
-
[55]
A. S. Nair, L. Foppa, and M. Scheffler, Scientific Data 12, 1518 (2025)
work page 2025
-
[56]
B. P. Klein, S. J. Hall, and R. J. Maurer, Journal of Physics: Condensed Matter33, 154005 (2021)
work page 2021
- [57]
- [58]
-
[59]
R. K. Pandey, Journal of Electronic Materials53, 6303 (2024)
work page 2024
-
[60]
W. G. Stark, C. L. Box, M. Sachs, N. Hertl, and R. J. Maurer, Physical Review B112, 035422 (2025)
work page 2025
-
[61]
A. Tellez-Mora, X. He, E. Bousquet, L. Wirtz, and A. H. Romero, npj Computational Materials10, 20 (2024)
work page 2024
-
[62]
A. K. Ibrahim and A. A. Al-Jobory, Journal of Molec- ular Modeling30, 409 (2024)
work page 2024
-
[63]
R. Latypov, S. Sozykin, and V. Beskachko, Journal of Surface Investigation: X-ray, Synchrotron and Neutron Techniques18, 63 (2024)
work page 2024
-
[64]
M. Biliroglu, M. T¨ ure, A. Ghita, M. Kotyrov, X. Qin, D. Seyitliyev, N. Phonthiptokun, M. Abdelsamei, J. Chai, R. Su,et al., Nature , 1 (2025)
work page 2025
-
[65]
T. J. Frankcombe, Physical Chemistry Chemical Physics27, 12041 (2025). 21 TABLE B.4. Occupation weighted band structure differences with PW basis (η 10-metrics) calculated on 24 bulk systems (in eV, Fermi levels are up-shifted for 10 eV) Bulk pVDZ pVTZ − pVTZ pVTZ-2V pVTZ-3V pVQZ = AlAs 0.228 0.068 0.052 0.047 0.034 0.042 AlN 0.106 0.050 0.023 0.020 0.027...
work page 2025
-
[66]
Z. Li, J. Huang, X. Ren, J. Li, R. Xiao, and H. Li, Journal of Power Sources640, 236632 (2025)
work page 2025
-
[67]
G. Ru, W. Qi, K. Xue, M. Wang, and X. Liu, Nanoscale 17, 7497 (2025)
work page 2025
-
[68]
Y. Liu, X. Liu, and M. Chen, Journal of Nuclear Mate- rials545, 152733 (2021)
work page 2021
-
[69]
Y. Liu, X. Ding, M. Chen, and S. Xu, Physical Chem- istry Chemical Physics24, 15511 (2022)
work page 2022
-
[70]
J. Luo, X. Zhu, X. Lian, Y. Zheng, R. Thottathil, W. Chen, S. Liu, A. Ariando, and J. Hu, 2D Materi- als12, 015022 (2024)
work page 2024
-
[71]
M. Sun, B. Jin, X. Yang, and S. Xu, Nature Communi- cations16, 3600 (2025)
work page 2025
-
[72]
F.-M. Jing, Z.-X. Shen, G.-Q. Qin, W.-K. Zhang, T. Lin, R. Cai, Z.-Z. Zhang, G. Cao, L. He, X.-X. Song,et al., Physical Review Applied23, 044053 (2025)
work page 2025
- [73]
- [74]
- [75]
-
[76]
Y. Tatetsu and Y. Gohda, in2018 IEEE International Magnetics Conference (INTERMAG)(IEEE, 2018) pp. 1–1
work page 2018
-
[77]
D. Bennett, M. Pizzochero, J. Junquera, and E. Kaxi- ras, Physical Review B111, 125123 (2025)
work page 2025
-
[78]
J. D. Talman, Journal of computational physics29, 35 (1978)
work page 1978
-
[79]
O. A. Sharafeddin, H. F. Bowen, D. J. Kouri, and D. K. Hoffman, Journal of Computational Physics100, 294 (1992)
work page 1992
-
[80]
M. Toyoda and T. Ozaki, Computer Physics Communi- cations181, 277 (2010)
work page 2010
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.