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.

Low-Order Scaling G0W0 by Pair Atomic Density Fitting

Low-order Scaling $G_0W_0$ by Pair Atomic Density Fitting

An O(N^3) implementation of Hedin's GW approximation for molecules

Scalable Molecular GW Calculations: Valence and Core Spectra

Two-component $GW$ calculations: Cubic scaling implementation and comparison of vertex corrected and partially self-consistent $GW$ variants

Solving multi-pole challenges in the GW100 benchmark enables precise low-scaling GW calculations

Linear scaling density functional calculations with Gaussian orbitals

Extrapolation of G0W0 energy levels from small basis sets for elements from H to Cl

Advancing non-atom-centered basis methods for more accurate interaction energies: benchmarks and large-scale applications

A Density-Based Basis-Set Incompleteness Correction for GW Methods

All-electron periodic G 0 W 0 implementation with numerical atomic orbital basis functions: Algorithm and benchmarks

Resolution-of-identity approach to Hartree–Fock, hybrid density functionals, RPA, MP2 and GW with numeric atom-centered orbital basis functions

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

Small Basis Set Density Functional Theory Method for Cost-efficient, Large-scale Condensed Matter Simulations

Automated Workflow for Accurate High-Throughput GW Calculations

AB-G0W0: A practical G0W0 method without frequency integration based on an auxiliary boson expansion

Accelerating Time-Dependent Density Functional Theory and GW Calculations for Molecules and Nanoclusters with Symmetry Adapted Interpolative Separable Density Fitting

Time-Dependent Density Functional Theory Excited-State Forces with the Gaussian and Augmented Plane Wave Method

Projected Augmented Waves (PAW) inspired mixed basis sets for small molecules

Toward Pair Atomic Density Fitting for Correlation Energies with Benchmark Accuracy

Realizing Effective Cubic-Scaling Coulomb Hole Plus Screened Exchange Approximation in Periodic Systems Via Interpolative Separable Density Fitting with a Plane-Wave Basis Set.