Abstract:A semi-empirical counterpoise-type correction for basis set superposition error (BSSE) in molecular systems is presented. An atom pair-wise potential corrects for the inter- and intra-molecular BSSE in supermolecular Hartree-Fock (HF) or density functional theory (DFT) calculations. This geometrical counterpoise (gCP) denoted scheme depends only on the molecular geometry, i.e., no input from the electronic wave-function is required and hence is applicable to molecules with ten thousands of atoms. The four necessary parameters have been determined by a fit to standard Boys and Bernadi counterpoise corrections for Hobza's S66×8 set of non-covalently bound complexes (528 data points). The method's target are small basis sets (e.g., minimal, split-valence, 6-31G*), but reliable results are also obtained for larger triple-ζ sets. The intermolecular BSSE is calculated by gCP within a typical error of 10%-30% that proves sufficient in many practical applications. The approach is suggested as a quantitative correction in production work and can also be routinely applied to estimate the magnitude of the BSSE beforehand. The applicability for biomolecules as the primary target is tested for the crambin protein, where gCP removes intramolecular BSSE effectively and yields conformational energies comparable to def2-TZVP basis results. Good mutual agreement is also found with Jensen's ACP(4) scheme, estimating the intramolecular BSSE in the phenylalanine-glycine-phenylalanine tripeptide, for which also a relaxed rotational energy profile is presented. A variety of minimal and double-ζ basis sets combined with gCP and the dispersion corrections DFT-D3 and DFT-NL are successfully benchmarked on the S22 and S66 sets of non-covalent interactions. Outstanding performance with a mean absolute deviation (MAD) of 0.51 kcal/mol (0.38 kcal/mol after D3-refit) is obtained at the gCP-corrected HF-D3/(minimal basis) level for the S66 benchmark. The gCP-corrected B3LYP-D3/6-31G* model chemistry yields MAD=0.68 kcal/mol, which represents a huge improvement over plain B3LYP/6-31G* (MAD=2.3 kcal/mol). Application of gCP-corrected B97-D3 and HF-D3 on a set of large protein-ligand complexes prove the robustness of the method. Analytical gCP gradients make optimizations of large systems feasible with small basis sets, as demonstrated for the inter-ring distances of 9-helicene and most of the complexes in Hobza's S22 test set. The method is implemented in a freely available FORTRAN program obtainable from the author's website.

Automatic generation of complementary auxiliary basis sets for explicitly correlated methods

Automatic Generation of Accurate and Cost-efficient Auxiliary Basis Sets

Automatic Generation of Auxiliary Basis Sets in Spherical Representation Using the Cholesky Decomposition

Basis-Set Limit CCSD(T) Energies for Large Molecules with Local Natural Orbitals and Reduced-Scaling Basis-Set Corrections

Auxiliary basis expansions for large-scale electronic structure calculations

Adaptive atomic basis sets

Systematically Improvable Optimized Atomic Basis Sets Forab Initiocalculations

Molecule-optimized Basis Sets and Hamiltonians for Accelerated Electronic Structure Calculations of Atoms and Molecules

The cc-pV5Z-F12 basis set: reaching the basis set limit in explicitly correlated calculations

Atomic basis functions for molecular electronic structure calculations

Correlation Consistent Gaussian Basis Sets for H, B–Ne with Dirac–Fock AREP Pseudopotentials: Applications in Quantum Monte Carlo Calculations

Large-scale ab initio simulations based on systematically improvable atomic basis

Approaching the Basis Set Limit for DFT Calculations Using an Environment-Adapted Minimal Basis with Perturbation Theory: Formulation, Proof of Concept, and a Pilot Implementation

Automatic Selection of Active Orbitals from Generalized Valence Bond Orbitals.

An explicitly correlated approach to basis set incompleteness in Full Configuration Interaction Quantum Monte Carlo

A geometrical correction for the inter- and intra-molecular basis set superposition error in Hartree-Fock and density functional theory calculations for large systems

Adaptive Algorithm for Electronic Structure Calculations Using Reduction of Gaussian Mixtures

Accelerating the Convergence of Auxiliary-Field Quantum Monte Carlo inSolids with Optimized Gaussian Basis Sets

Efficient and accurate local approximations to coupled-electron pair approaches: An attempt to revive the pair natural orbital method

An efficient method by combining different basis sets and SAPT levels