Search papers, labs, and topics across Lattice.
This paper introduces a smoothed truncated Coulomb (sTC) potential that enhances the evaluation of electron repulsion integrals (ERIs) in periodic Hartree–Fock calculations using Gaussian basis sets. By employing a real-space Gaussian convolution to smooth truncation boundaries, the method allows for a dual-space algorithm that is efficient for both pseudopotential and all-electron calculations. Benchmark results demonstrate that sTC significantly improves thermodynamic-limit convergence across various material systems while closely matching traditional truncated Coulomb results, thereby broadening the applicability of HF methods in computational chemistry.
Smoothing the truncation boundary in Coulomb potentials can drastically enhance convergence in Hartree–Fock calculations without sacrificing computational efficiency.
Truncated Coulomb (TC) potentials reduce finite-size errors and accelerate thermodynamic-limit convergence in periodic Hartree--Fock (HF) calculations, but their use with Gaussian basis sets is complicated by the evaluation of electron repulsion integrals (ERIs), particularly for nonspherical truncation domains and all-electron calculations. We introduce the smoothed truncated Coulomb (sTC) potential as a systematically improvable approximation to a parent TC potential. A real-space Gaussian convolution smooths the sharp truncation boundary, and a single dimensionless parameter, $η$, controls the width of the smoothing window, which can be tightened to systematically approach the TC reference. The smoothing by sTC enables a dual-space algorithm for evaluating periodic Gaussian-basis ERIs that requires neither a large plane-wave basis nor new molecular integral kernels and is applicable to both pseudopotential and all-electron calculations, including the important Wigner--Seitz-cell truncation boundaries. Benchmarks spanning insulating, semiconducting, layered, metallic, and molecular-crystal systems show that sTC-based HF closely reproduces TC results in pseudopotential calculations and extends TC-quality calculations to all-electron settings. Across these systems, sTC substantially improves thermodynamic-limit convergence relative to the probe-charge Ewald method while remaining practical for all-electron calculations in which direct TC-based calculations are computationally challenging.
