This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange–correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear–electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an “open teamware” model and an increasingly modular design.

An energy decomposition analysis (EDA) separates a calculated interaction energy into as many interpretable contributions as possible; for instance, permanent and induced electrostatics, Pauli repulsions, dispersion and charge transfer. The challenge is to construct satisfactory definitions of all terms in the chemically relevant regime where fragment densities overlap, rendering unique definitions impossible. Towards this goal, we present an improved EDA for Kohn-Sham density functional theory (DFT) with properties that have previously not been simultaneously attained. Building on the absolutely localized molecular orbital (ALMO)-EDA, this second generation ALMO-EDA is variational and employs valid antisymmetric electronic wavefunctions to produce all five contributions listed above. These contributions moreover all have non-trivial complete basis set limits. We apply the EDA to the water dimer, the T-shaped and parallel-displaced benzene dimer, the p-biphthalate dimer "anti-electrostatic" hydrogen bonding complex, the biologically relevant binding of adenine and thymine in stacked and hydrogen-bonded configurations, the triply hydrogen-bonded guanine-cytosine complex, the interaction of Cl(-) with s-triazine and with the 1,3-dimethyl imidazolium cation, which is relevant to the study of ionic liquids, and the water-formaldehyde-vinyl alcohol ter-molecular radical cationic complex formed in the dissociative photoionization of glycerol.

In energy decomposition analysis of Kohn-Sham density functional theory calculations, the so-called frozen (or pre-polarization) interaction energy contains contributions from permanent electrostatics, dispersion, and Pauli repulsion. The standard classical approach to separate them suffers from several well-known limitations. We introduce an alternative scheme that employs valid antisymmetric electronic wavefunctions throughout and is based on the identification of individual fragment contributions to the initial supersystem wavefunction as determined by an energetic optimality criterion. The density deformations identified with individual fragments upon formation of the initial supersystem wavefunction are analyzed along with the distance dependence of the new and classical terms for test cases that include the neon dimer, ammonia borane, water-Na(+), water-Cl(-), and the naphthalene dimer.

AMOEBA is a molecular mechanics force field that addresses some of the shortcomings of a fixed partial charge model, by including permanent atomic point multipoles through quadrupoles, as well as many-body polarization through the use of point inducible dipoles. In this work, we investigate how well AMOEBA formulates its non-bonded interactions, and how it implicitly incorporates quantum mechanical effects such as charge penetration (CP) and charge transfer (CT), for water-water and water-ion interactions. We find that AMOEBA's total interaction energies, as a function of distance and over angular scans for the water dimer and for a range of water-monovalent cations, agree well with an advanced density functional theory (DFT) model, whereas the water-halides and water-divalent cations show significant disagreement with the DFT result, especially in the compressed region when the two fragments overlap. We use a second-generation energy decomposition analysis (EDA) scheme based on absolutely localized molecular orbitals (ALMOs) to show that in the best cases AMOEBA relies on cancellation of errors by softening of the van der Waals (vdW) wall to balance permanent electrostatics that are too unfavorable, thereby compensating for the missing CP effect. CT, as another important stabilizing effect not explicitly taken into account in AMOEBA, is also found to be incorporated by the softened vdW interaction. For the water-halides and water-divalent cations, this compensatory approach is not as well executed by AMOEBA over all distances and angles, wherein permanent electrostatics remains too unfavorable and polarization is overdamped in the former while overestimated in the latter. We conclude that the DFT-based EDA approach can help refine a next-generation AMOEBA model that either realizes a better cancellation of errors for problematic cases like those illustrated here, or serves to guide the parametrization of explicit functional forms for short-range contributions from CP and/or CT.

Energy decomposition analysis (EDA) of electronic structure calculations has facilitated quantitative understanding of diverse intermolecular interactions. Nevertheless, such analyses are usually performed at a single geometry and thus decompose a "single-point" interaction energy. As a result, the influence of the physically meaningful EDA components on the molecular structure and other properties are not directly obtained. To address this gap, the absolutely localized molecular orbital (ALMO)-EDA is reformulated in an adiabatic picture, where the frozen, polarization, and charge transfer energy contributions are defined as energy differences between the stationary points on different potential energy surfaces (PESs), which are accessed by geometry optimizations at the frozen, polarized and fully relaxed levels of density functional theory (DFT). Other molecular properties such as vibrational frequencies can thus be obtained at the stationary points on each PES. We apply the adiabatic ALMO-EDA to different configurations of the water dimer, the water-Cl and water-Mg/Ca complexes, metallocenes (Fe, Ni, Cu, Zn), and the ammonia-borane complex. This method appears to be very useful for unraveling how physical effects such as polarization and charge transfer modulate changes in molecular properties induced by intermolecular interactions. As an example of the insight obtained, we find that a linear hydrogen bond geometry for the water dimer is preferred even without the presence of polarization and charge transfer, while the red shift in the OH stretch frequency is primarily a charge transfer effect; by contrast, a near-linear geometry for the water-chloride hydrogen bond is achieved only when charge transfer is allowed.

Charge-transfer (CT) is an important binding force in the formation of intermolecular complexes, and there have been a variety of theoretical models proposed to quantify this effect. These approaches, which typically rely on a definition of a "CT-free" state based on a partition of the system, sometimes yield significantly different results for a given intermolecular complex. Two widely used definitions of the "CT-free" state, the absolutely localized molecular orbitals (ALMO) method (where only on-fragment orbital mixings are permitted) and the constrained density functional theory (CDFT) approach (where fragment electron populations are fixed), are carefully examined in this work. Natural bond orbital (NBO) and the regularized symmetry-adapted perturbation theory (SAPT) are also briefly considered. Results for the ALMO and CDFT definitions of CT are compared on a broad range of model systems, including hydrogen-bonding systems, borane complexes, metal-carbonyl complexes, and complexes formed by water and metal cations. For most of these systems, CDFT yields a much smaller equilibrium CT energy compared to that given by the ALMO-based definition. This is mainly because the CDFT population constraint does not fully inhibit CT, which means that the CDFT "CT-free" state is in fact CT-contaminated. Examples of this contamination include (i) matching forward and backward donation (e.g., formic acid dimer) and (ii) unidirectional CT without changing fragment populations. The magnitude of the latter effect is quantified in systems such as the water dimer by employing a 3-space density constraint in addition to the orbital constraint. Furthermore, by means of the adiabatic EDA, it is shown that several observable effects of CT, such as the "pyramidalization" of the planar BH molecule upon the complexation with Lewis bases, already appear on the "CT-free" CDFT surface. These results reveal the essential distinctions between the ALMO and CDFT definitions of CT and suggest that the former is more consistent with accepted understanding of the role of CT in intermolecular binding.

Quantum chemistry in the form of density functional theory (DFT) calculations is a powerful numerical experiment for predicting intermolecular interaction energies. However, no chemical insight is gained in this way beyond predictions of observables. Energy decomposition analysis (EDA) can quantitatively bridge this gap by providing values for the chemical drivers of the interactions, such as permanent electrostatics, Pauli repulsion, dispersion, and charge transfer. These energetic contributions are identified by performing DFT calculations with constraints that disable components of the interaction. This review describes the second-generation version of the absolutely localized molecular orbital EDA (ALMO-EDA-II). The effects of different physical contributions on changes in observables such as structure and vibrational frequencies upon complex formation are characterized via the adiabatic EDA. Example applications include red- versus blue-shifting hydrogen bonds; the bonding and frequency shifts of CO, N2, and BF bound to a [Ru(II)(NH3)5]2 + moiety; and the nature of the strongly bound complexes between pyridine and the benzene and naphthalene radical cations. Additionally, the use of ALMO-EDA-II to benchmark and guide the development of advanced force fields for molecular simulation is illustrated with the recent, very promising, MB-UCB potential. Expected final online publication date for the Annual Review of Physical Chemistry, Volume 72 is April 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

We have designed an energy decomposition analysis (EDA) to gain a deeper understanding of single chemical bonds, that is, those in which the interacting fragments are doublet open-shell systems but the supersystem is closed-shell. The method is a spin-pure extension of the absolutely localized molecular orbital (ALMO) EDA to the one-pair perfect pairing energy (equivalently to an active space of two electrons in two orbitals). The total interaction energy is broken up into four terms: frozen interactions, spin-coupling, polarization, and charge-transfer. A variety of single bonds are analyzed and, in addition, we use this method to show how solvation changes the nature of bonds, producing different results in the gas-phase and with explicit solvent molecules.