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.

Minimal Basis Set Hartree-Fock Corrected with Atom-Centered Potentials for Molecular Crystal Modeling and Crystal Structure Prediction

Dispersion corrected hartree-fock and density functional theory for organic crystal structure prediction

Crystal structure predictions for molecules with soft degrees of freedom using intermonomer force fields derived from first principles

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

Large-scale parameter estimation for Crystal Structure Prediction. Part 1: Dataset, Methodology, and Implementation

Powder-Diffraction-Based Structural Comparison for Crystal Structure Prediction without Prior Indexing

How Accurate Can Crystal Structure Predictions Be for High-Energy Molecular Crystals?

Predictive crystallography at scale: mapping, validating, and learning from 1,000 crystal energy landscapes

Reliable crystal structure predictions from first principles

A method for the automatic generation of a minimal basis set of structural templates for material phase-space exploration

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

CSPBench: a benchmark and critical evaluation of Crystal Structure Prediction

Reliable and Practical Computational Prediction of Molecular Crystal Polymorphs

Powder‐diffraction‐based structural comparison for crystal structure prediction without prior indexing

Enhancing crystal structure prediction by combining computational and experimental data via graph networks

Improved Description of the Structure of Molecular and Layered Crystals: Ab Initio DFT Calculations with van der Waals Corrections

Geometrical correction for the inter- and intramolecular basis set superposition error in periodic density functional theory calculations

Latent Conservative Objective Models for Data-Driven Crystal Structure Prediction

Assessment of the exchange-hole dipole moment dispersion correction for the energy ranking stage of the seventh crystal structure prediction blind test