Abstract:We present a purely numerical approach in Cartesian grid, for efficient computation of Hartree-Fock (HF) exchange contribution in the HF and density functional theory models. This takes inspiration from a recently developed algorithm [Liu \emph{et al.}, J.~Chem.~Theor.~Comput.~ \textbf{13}, 2571 (2017)]. A key component is the accurate evaluation of electrostatic potential integral, which is the rate-determining step. This introduces the Fourier convolution theorem in conjunction with a range-separated Coulomb interaction kernel. The latter is efficiently mapped into real grid through a simple optimization procedure, giving rise to a constraint in the range-separated parameter. The overall process offers logarithmic scaling with respect to molecular size. It is then extended towards global hybrid functionals such as B3LYP, PBE0 and BHLYP within pseudopotential Kohn-Sham theory, through an LCAO-MO ansatz in Cartesian grid, developed earlier in our laboratory. For sake of comparison, a parallel semi-numerical approach has also been worked out that exploits the familiar Obara-Saika recursion algorithm. An excellent agreement between these two routes is demonstrated through total energy and orbital energy in a series of atoms and molecules (including 10 $\pi$-electron molecules), employing an LANL2DZ-type basis function. A critical analysis of these two algorithms reveals that the proposed numerical scheme could lead to very attractive and competitive scaling. The success of our approach also enables us for further development of optimally tuned range-separated hybrid and hyper functionals.

Elongation method for linear scaling

Linear Scaling Electronic Structure Calculations With Numerical Atomic Basis Set

A Nonempirical Scaling Correction Approach for Density Functional Methods Involving Substantial Amount of Hartree-Fock Exchange.

Linear Scaling Algorithm for Tight-Binding Molecular Dynamics Simulations.

Linear scaling DFT calculations for large Tungsten systems using an optimized local basis

Toward linear scaling auxiliary field quantum Monte Carlo with local natural orbitals

Efficient and Robust Ab Initio Self-Consistent Field Acceleration Algorithm Based on a Semiempirical Model Hamiltonian

Self-consistent field theory based molecular dynamics with linear system-size scaling

Optimization of Convergence Criteria for Fragmentation Methods

Comparison Study of Finite Element and Basis Set Methods for Finite Size Scaling

Linear-Scaling Subspace-Iteration Algorithm with Optimally Localized Nonorthogonal Wave Functions for Kohn-Sham Density Functional Theory

Parallel Implementation of Large-Scale Linear Scaling Density Functional Theory Calculations with Numerical Atomic Orbitals in HONPAS.

Linear Scaling Local Correlation Approach for Solving the Coupled Cluster Equations of Large Systems.

A Sub-Linear Scaling Algorithm for Computing the Electronic Structure of Materials

Linear scaling computation of the Fock matrix VII. Periodic Density Functional Theory at the $Γ$-point

Linear Scaling Solution of the Time-Dependent Self-Consistent-Field Equations

Efficient HF exchange evaluation through Fourier convolution in Cartesian grid for orbital-dependent density functionals

An evaluation of scaling factors for multiparameter scaling procedures based on DFT force fields

Efficient parallel linear scaling method to get the response density matrix in all-electron real-space density-functional perturbation theory

Linear-scaling ab-initio calculations for large and complex systems

Linear Scaling Calculation of Maximally Localized Wannier Functions with Atomic Basis Set.