Use of multigrids to reduce the cost of performing interpolative separable density fitting (2404.09373v1)
Abstract: In this article, we present an interpolative separable density fitting (ISDF) based algorithm to calculate exact exchange in periodic mean field calculations. In the past, decomposing the two-electron integrals into tensor hypercontraction (THC) form using ISDF was the most expensive step of the entire mean field calculation. Here we show that by using a multigrid-ISDF algorithm both the memory and the CPU cost of this step can be reduced. The CPU cost is brought down from cubic scaling to quadratic scaling with a low computational prefactor which reduces the cost by almost two orders of magnitude. Thus, in the new algorithm, the cost of performing ISDF is largely negligible compared to other steps. Along with the CPU cost, the memory cost of storing the factorized two-electron integrals is also reduced by a factor of up to 35. With the current algorithm, we can perform Hartree-Fock calculations on a Diamond supercell containing more than 17,000 basis functions and more than 1,500 electrons on a single node with no disk usage. For this calculation, the cost of constructing the exchange matrix is only a factor of four slower than the cost of diagonalizing the Fock matrix. Augmenting our approach with linear scaling algorithms can further speed up the calculations.
- Self-consistent equations including exchange and correlation effects. Physical review, 140(4A):A1133, 1965.
- A D Becke. No Title. J. Chem. Phys., 98:5648, 1993.
- Thirty years of density functional theory in computational chemistry: an overview and extensive assessment of 200 density functionals. Molecular physics, 115(19):2315–2372, 2017.
- Energy band gaps and lattice parameters evaluated with the Heyd-Scuseria-Ernzerhof screened hybrid functional. Journal of Chemical Physics, 123(17):1–8, 2005.
- Screened hybrid density functionals applied to solids. The Journal of chemical physics, 124(15), 2006.
- Excess electron states in reduced bulk anatase TiO2: Comparison of standard GGA, GGA+U, and hybrid DFT calculations. Journal of Chemical Physics, 129(15):154113, 2008.
- Accurate band gaps for semiconductors from density functional theory. Journal of Physical Chemistry Letters, 2(3):212–217, 2011.
- Stability of non-metal dopants to tune the photo-absorption of TiO2 at realistic temperatures and oxygen partial pressures: A hybrid DFT study. Scientific Reports, 9(1):20–26, 2019.
- Photocatalytic co2 reduction: A review of ab initio mechanism, kinetics, and multiscale modeling simulations. ACS catalysis, 10(24):14984–15007, 2020.
- Principles for a direct SCF approach to LICAO-MO ab-initio calculations. Journal of Computational Chemistry, 3(3):385–399, 1982.
- The continuous fast multipole method. Chemical Physics Letters, 230(1-2):8–16, 1994.
- Fast assembly of the coulomb matrix: A quantum chemical tree code. The Journal of chemical physics, 104(12):4685–4698, 1996.
- A fast multipole method for periodic systems with arbitrary unit cell geometries. Chemical Physics Letters, 283(1):61–68, 1998.
- Linear scaling computation of the Fock matrix. The Journal of Chemical Physics, 106(13):5526, 1997.
- Linear and sublinear scaling formation of Hartree-Fock-type exchange matrices. The Journal of Chemical Physics, 109(5):1663, 1998.
- Stefan Goedecker. Linear scaling electronic structure methods. Reviews of Modern Physics, 71(4):1085–1123, 1999.
- Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance. Journal of Chemical Theory and Computation, 16(6):3757–3785, 2020.
- Fast evaluation of the Coulomb potential for electron densities using multipole accelerated resolution of identity approximation. The Journal of Chemical Physics, 118:9136, 2003.
- Linear scaling density fitting. Journal of Chemical Physics, 125(19):194109, 2006.
- Fast Hartree-Fock theory using local density fitting approximations. Molecular Physics, 102(21-22):2311–2321, 2004.
- Hartree-Fock exchange computed using the atomic resolution of the identity approximation. Journal of Chemical Physics, 128:104106, 2008.
- Efficient implementation of the pair atomic resolution of the identity approximation for exact exchange for hybrid and range-separated density functionals. Journal of Chemical Theory and Computation, 11(2):518–527, 2015.
- Fast, accurate evaluation of exact exchange: The occ-RI-K algorithm. Journal of Chemical Physics, 143(2):024113, 2015.
- B. I. Dunlap. Robust variational fitting: Gaspar’s variational exchange can accurately be treated analytically. Journal of Molecular Structure: THEOCHEM, 501-502(August):221–228, 2000.
- B. I. Dunlap. Robust and variational fitting: Removing the four-center integrals from center stage in quantum chemistry. Journal of Molecular Structure: THEOCHEM, 529(1-3):37–40, 2000.
- Fast construction of the exchange operator in an atom-centred basis with concentric atomic density fitting. Molecular Physics, 115(17-18):2065–2076, 2017.
- Accurate localized resolution of identity approach for linear-scaling hybrid density functionals and for many-body perturbation theory. New Journal of Physics, 17(9):093020, sep 2015.
- Efficient all-electron hybrid density functionals for atomistic simulations beyond 10,000 atoms, 2024.
- NHF Beebe and J Linderberg. Simplifications in the generation and transformation of two-electron integrals in molecular calculations. International Journal of Quantum Chemistry, XII:683–705, 1977.
- Unbiased auxiliary basis sets for accurate two-electron integral approximations. Journal of Chemical Physics, 127(11):114107, 2007.
- Richard A. Friesner. Solution of self-consistent field electronic structure equations by a pseudospectral method. Chemical Physics Letters, 116(1):39–43, 1985.
- Efficient, approximate and parallel Hartree-Fock and hybrid DFT calculations. A ’chain-of-spheres’ algorithm for the Hartree-Fock exchange. Chemical Physics, 356(1-3):98–109, 2009.
- Communication: Tensor hypercontraction. III. Least-squares tensor hypercontraction for the determination of correlated wavefunctions. Journal of Chemical Physics, 137:221101, 2012.
- Tensor hypercontraction. II. Least-squares renormalization. Journal of Chemical Physics, 137(22):224106, 2012.
- Tensor hypercontraction density fitting. I. Quartic scaling second- and third-order Møller-Plesset perturbation theory. Journal of Chemical Physics, 137(4):1085, 2012.
- Compression of the electron repulsion integral tensor in tensor hypercontraction format with cubic scaling cost. Journal of Computational Physics, 302:329–335, 2015.
- Interpolative Separable Density Fitting Decomposition for Accelerating Hybrid Density Functional Calculations with Applications to Defects in Silicon. Journal of Chemical Theory and Computation, 13(11):5420–5431, 2017.
- Interpolative Separable Density Fitting through Centroidal Voronoi Tessellation with Applications to Hybrid Functional Electronic Structure Calculations. Journal of Chemical Theory and Computation, 14(3):1311–1320, 2018.
- Systematically Improvable Tensor Hypercontraction: Interpolative Separable Density-Fitting for Molecules Applied to Exact Exchange, Second- A nd Third-Order Møller-Plesset Perturbation Theory. Journal of Chemical Theory and Computation, 16(1):243–263, 2020.
- Fast Exchange with Gaussian Basis Set Using Robust Pseudospectral Method. Journal of Chemical Theory and Computation, 18(12):7306–7320, 2022.
- Even faster exact exchange for solids via tensor hypercontraction. arXiv preprint arXiv:2304.05505, 2023.
- Machine learning k-means clustering of interpolative separable density fitting algorithm for accurate and efficient cubic-scaling exact exchange plus random phase approximation within plane waves. Journal of Chemical Theory and Computation, 2024.
- D. R. Bowler and T. Miyazaki. O(N) methods in electronic structure calculations. Reports on Progress in Physics, 75(3):036503, 2012.
- Order- N implementation of exact exchange in extended insulating systems. Physical Review B - Condensed Matter and Materials Physics, 79:085102, 2009.
- Self-averaging stochastic kohn-sham density-functional theory. Physical Review Letters, 111:106402, 2013.
- Stochastic Optimally Tuned Range-Separated Hybrid Density Functional Theory. Journal of Physical Chemistry A, 120(19):3071–3078, 2016.
- Deterministic/fragmented-stoclhastic exchange for large-scale hybrid dft calculations. Journal of Chemical Theory and Computation, 19(24):9239–9247, 2023.
- Molecular dynamics simulation of liquid water: Hybrid density functionals. The Journal of Physical Chemistry B, 110(8):3685–3691, 2006.
- Lin Lin. Adaptively Compressed Exchange Operator. Journal of Chemical Theory and Computation, 12(5):2242–2249, 2016.
- Auxiliary Density Matrix Methods for Hartree-Fock Exchange Calculations. Journal of Chemical Theory and Computation, 6(8):2348–2364, 2010.
- Interpolative separable density fitting decomposition for accelerating hartree–fock exchange calculations within numerical atomic orbitals. The Journal of Physical Chemistry A, 124(27):5664–5674, 2020.
- The Gaussian and augmented-plane-wave density functional method for ab initio molecular dynamics simulations. Theoretical Chemistry Accounts, 103(2):124–140, 1999.
- Quickstep: Fast and accurate density functional calculations using a mixed Gaussian and plane waves approach. Computer Physics Communications, 167(2):103–128, 2005.
- An efficient real space multigrid qm/mm electrostatic coupling. Journal of Chemical Theory and Computation, 1(6):1176–1184, 2005.
- Thomas L Beck. 5 real-space and multigrid methods in computational chemistry. Reviews in Computational Chemistry, 26:223, 2009.
- Second-order møller–plesset perturbation theory in the condensed phase: An efficient and massively parallel gaussian and plane waves approach. Journal of chemical theory and computation, 8(11):4177–4188, 2012.
- CP2K: An electronic structure and molecular dynamics software package -Quickstep: Efficient and accurate electronic structure calculations. Journal of Chemical Physics, 152:194103, 2020.
- Attractive electron–electron interactions within robust local fitting approximations. Journal of Computational Chemistry, 34(17):1486–1496, 2013.
- On resolution-of-the-identity electron repulsion integral approximations and variational stability. Journal of Chemical Theory and Computation, 13(10):4897–4906, 2017.
- Finding Structure with Randomness: Probabilistic Algorithms for Constructing Approximate Matrix Decompositions. SIAM Review, 53(2):217–288, 2011.
- Randomized algorithms for the low-rank approximation of matrices. Proceedings of the National Academy of Sciences, 104(51):20167–20172, 2007.
- Devin A. Matthews. Improved grid optimization and fitting in least squares tensor hypercontraction. Journal of Chemical Theory and Computation, 16(3):1382–1385, 2020.
- Accurate molecular integrals and energies using combined plane wave and gaussian basis sets in molecular electronic structure theory. The Journal of chemical physics, 116(18):7795–7805, 2002.
- The Fourier transform Coulomb method: Efficient and accurate calculation of the Coulomb operator in a Gaussian basis. The Journal of Chemical Physics, 117(17):7827–7835, 11 2002.
- PySCF: the Python-based simulations of chemistry framework. Wiley Interdisciplinary Reviews: Computational Molecular Science, page (published online), 2017.
- Separable dual-space gaussian pseudopotentials. Phys. Rev. B, 54:1703–1710, Jul 1996.
- Relativistic separable dual-space gaussian pseudopotentials from h to rn. Phys. Rev. B, 58:3641–3662, Aug 1998.
- Juerg Hutter. New optimization of gth pseudopotentials for pbe, scan, pbe0 functionals. gth pseudopotentials for hartree-fock. nlcc pseudopotentials for pbe. https://github.com/juerghutter/gth, 2019.
- Correlation-consistent gaussian basis sets for solids made simple. Journal of Chemical Theory and Computation, 18(3):1595–1606, 2022. PMID: 35192359.