Reply to Comment on "Efficient implementation of the superposition of atomic potentials initial guess for electronic structure calculations in Gaussian basis sets"
Abstract
We reply to the Comment by Surjuse, Deng, Asadchev, and Valeev [arXiv:2603.16989] on our article "Efficient implementation of the superposition of atomic potentials initial guess for electronic structure calculations in Gaussian basis sets" [J. Chem. Phys. 152, 144105 (2020)]. We do not disagree with their work, but point out that they appear to have misinterpreted a statement in our original abstract: our suggestion to evaluate the matrix elements of the superposition of atomic potentials (SAP) using two-electron integral code was made on grounds of portability rather than as a recommendation against a dedicated one-electron implementation. Such an implementation is in fact already available in existing codes - in Dirac, for example, where a Gaussian SAP potential is formally indistinguishable from a Gaussian finite-nucleus model at the level of the nuclear attraction integral routine. We further note that both routes reduce to the same three-Gaussian overlap kernel through the Obara-Saika Laplace transform, and that SAP itself enjoys a rigorous theoretical justification through Theophilou's theorem on spherical densities and its extensions by Nagy.
Related
- Training of Neural Network Machine Learning Potentials for Activation Barrier Calculations in Predicting Growth Modes of Deposited Metals
- First-Principles Calculations of Structural, Electronic, Optical, Vibrational, Dielectric, and Elastic Properties of Tetragonal Mg2TiO4 and Mg2TiO4:Mn4+
- Exchange-Correlation Potentials and Energies from Inverse Generalized Kohn-Sham Calculations
- DeepH-pack: a general-purpose neural network package for deep-learning electronic structure calculations
- Biquaternion Algebra with Bilinear Multiplication: A General, Elegant, and Computationally Advantageous Framework for Relativistic Electronic Structure Calculations on CPUs and GPUs
BibTeX
@misc{lehtola2026reply,
title = {{Reply to Comment on "Efficient implementation of the superposition of atomic potentials initial guess for electronic structure calculations in Gaussian basis sets"}},
author = {Susi Lehtola and Lucas Visscher and Eberhard Engel},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.25160},
eprint = {2609.25160},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.25160},
}
RIS
TY - UNPB
TI - Reply to Comment on "Efficient implementation of the superposition of atomic potentials initial guess for electronic structure calculations in Gaussian basis sets"
AU - Susi Lehtola
AU - Lucas Visscher
AU - Eberhard Engel
PY - 2026
DA - 2026/09/21
PB - arXiv
DO - 10.48550/arxiv.2609.25160
UR - https://arxiv.org/abs/2609.25160
AB - We reply to the Comment by Surjuse, Deng, Asadchev, and Valeev [arXiv:2603.16989] on our article "Efficient implementation of the superposition of atomic potentials initial guess for electronic structure calculations in Gaussian basis sets" [J. Chem. Phys. 152, 144105 (2020)]. We do not disagree with their work, but point out that they appear to have misinterpreted a statement in our original abstract: our suggestion to evaluate the matrix elements of the superposition of atomic potentials (SAP) using two-electron integral code was made on grounds of portability rather than as a recommendation against a dedicated one-electron implementation. Such an implementation is in fact already available in existing codes - in Dirac, for example, where a Gaussian SAP potential is formally indistinguishable from a Gaussian finite-nucleus model at the level of the nuclear attraction integral routine. We further note that both routes reduce to the same three-Gaussian overlap kernel through the Obara-Saika Laplace transform, and that SAP itself enjoys a rigorous theoretical justification through Theophilou's theorem on spherical densities and its extensions by Nagy.
ER -