PSI4 is a free and open-source ab initio electronic structure program providing implementations of Hartree–Fock, density functional theory, many-body perturbation theory, configuration interaction, density cumulant theory, symmetry-adapted perturbation theory, and coupled-cluster theory. Most of the methods are quite efficient, thanks to density fitting and multi-core parallelism. The program is a hybrid of C++ and Python, and calculations may be run with very simple text files or using the Python API, facilitating post-processing and complex workflows; method developers also have access to most of PSI4’s core functionalities via Python. Job specification may be passed using The Molecular Sciences Software Institute (MolSSI) QCSCHEMA data format, facilitating interoperability. A rewrite of our top-level computation driver, and concomitant adoption of the MolSSI QCARCHIVE INFRASTRUCTURE project, makes the latest version of PSI4 well suited to distributed computation of large numbers of independent tasks. The project has fostered the development of independent software components that may be reused in other quantum chemistry programs.

A systematic examination of the computational expense and accuracy of Symmetry-Adapted Perturbation Theory (SAPT) for the prediction of non-covalent interaction energies is provided with respect to both method [SAPT0, DFT-SAPT, SAPT2, SAPT2+, SAPT2+(3), and SAPT2+3; with and without CCD dispersion for the last three] and basis set [Dunning cc-pVDZ through aug-cc-pV5Z wherever computationally tractable, including truncations of diffuse basis functions]. To improve accuracy for hydrogen-bonded systems, we also include two corrections based on exchange-scaling (sSAPT0) and the supermolecular MP2 interaction energy (δMP2). When considering the best error performance relative to computational effort, we recommend as the gold, silver, and bronze standard of SAPT: SAPT2+(3)δMP2/aug-cc-pVTZ, SAPT2+/aug-cc-pVDZ, and sSAPT0/jun-cc-pVDZ. Their respective mean absolute errors in interaction energy across the S22, HBC6, NBC10, and HSG databases are 0.15 (62.9), 0.30 (4.4), and 0.49 kcal mol⁻¹ (0.03 h for adenine·thymine complex).

A systematic study of techniques for treating noncovalent interactions within the computationally efficient density functional theory (DFT) framework is presented through comparison to benchmark-quality evaluations of binding strength compiled for molecular complexes of diverse size and nature. In particular, the efficacy of functionals deliberately crafted to encompass long-range forces, a posteriori DFT+dispersion corrections (DFT-D2 and DFT-D3), and exchange-hole dipole moment (XDM) theory is assessed against a large collection (469 energy points) of reference interaction energies at the CCSD(T) level of theory extrapolated to the estimated complete basis set limit. The established S22 [revised in J. Chem. Phys. 132, 144104 (2010)] and JSCH test sets of minimum-energy structures, as well as collections of dispersion-bound (NBC10) and hydrogen-bonded (HBC6) dissociation curves and a pairwise decomposition of a protein–ligand reaction site (HSG), comprise the chemical systems for this work. From evaluations of accuracy, consistency, and efficiency for PBE-D, BP86-D, B97-D, PBE0-D, B3LYP-D, B970-D, M05-2X, M06-2X, ωB97X-D, B2PLYP-D, XYG3, and B3LYP-XDM methodologies, it is concluded that distinct, often contrasting, groups of these elicit the best performance within the accessible double-ζ or robust triple-ζ basis set regimes and among hydrogen-bonded or dispersion-dominated complexes. For overall results, M05-2X, B97-D3, and B970-D2 yield superior values in conjunction with aug-cc-pVDZ, for a mean absolute deviation of 0.41 – 0.49 kcal/mol, and B3LYP-D3, B97-D3, ωB97X-D, and B2PLYP-D3 dominate with aug-cc-pVTZ, affording, together with XYG3/6-311+G(3df,2p), a mean absolute deviation of 0.33 – 0.38 kcal/mol.