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.

Biological and heterogenous catalysts for the electrochemical CO2 Reduction Reaction (CO2RR) often exhibit a high degree of electronic delocalization that serves to minimize overpotential and maximize selectivity over the hydrogen evolution reaction (HER). Here, we report a molecular iron(II) system that captures this design concept in a homogeneous setting through the use of a redox non-innocent terpyridine-based pentapyridine ligand (tpyPY2Me). As a result of strong metal-ligand exchange coupling between the Fe(II) center and ligand, [Fe(tpyPY2Me)] 2+ exhibits redox behavior at potentials 640 mV more positive than the isostructural [Zn(tpyPY2Me)] 2+ analog containing the redox-inactive Zn(II) ion. This shift in redox potential is attributed to the requirement for both an open-shell metal ion and a redox non-innocent ligand. The metalligand cooperativity in [Fe(tpyPY2Me)] 2+ drives the electrochemical reduction of CO2 to CO at low overpotentials with high selectivity for CO2RR (> 90%) and turnover frequencies of 100,000 s -1 with no degradation over 20 h. The decrease in the thermodynamic barrier engendered by this coupling also enables homogeneous CO2 reduction catalysis in water without compromising selectivity or rates. Synthesis of the two-electron reduction product, [Fe(tpyPY2Me)] 0 , and characterization by X-ray crystallography, Mössbauer spectroscopy, X-ray absorption spectroscopy (XAS), variable temperature NMR, and density functional theory (DFT) calculations, support assignment of an open-shell singlet electronic structure that maintains a formal Fe(II) oxidation state with a doubly-reduced ligand system. This work provides a starting point for the design of systems that exploit metal-ligand cooperativity for electrocatalysis where the electrochemical potential of redox non-innocent ligands can be tuned through secondary metal-dependent interactions.

Second order Møller-Plesset theory provides a remarkably simple form for the electron correlation energy with many desirable properties, e.g. it is size-consistent, free of self-interaction error, and scales with the fifth power of system size. However, MP2 exhibits well-known shortcomings including an incomplete description of dispersion interactions and sizable failures for transition metal chemistry. Herein, we first explore multiple physically justified forms of single-parameter regularization and then demonstrate that with appropriate parameter choice, regularized MP2 with Hartree-Fock reference orbitals yields high and transferable accuracy across a wide variety of noncovalent interactions (S22, S66, XB40, A24, and L7 test sets) and (mostly closedshell) transition metal thermochemistry (metal-carbonyl dissociations and a subset of MOR41). We find that, especially for systems with interacting π systems relevant to dispersion interactions and dative bonding, regularization serves to damp overestimated pair-wise additive contributions to the first-order amplitudes that affect correlation energy and charge-density. The optimal parameter values for the noncovalent and transition metal sets are 1.1 and 0.4 for two regularizers, κ and σ 2 , respectively. These two regularizers slightly degrade the accuracy of conventional MP2 for some small-molecule test sets which are well-known to be sensitive to charge-density distribution (radical stabilization energies, barrier heights, dipole moments, and polarizabilities), most of which have relatively large gaps. Due to the relatively straightforward implementations of nuclear gradient and other properties, we recommend κ-MP2 with κ = 1.1 as a more accurate alternative to conventional MP2 and other related variants. Our results suggest that appropriately regularized MP2 models represent promising forms for the nonlocal correlation part of double hybrid density functionals, at no additional cost over conventional MP2.

Given the piecewise approach to modeling intermolecular interactions for force fields, they can be difficult to parameterize since they are fit to data like total energies that only indirectly connect to their separable functional forms. Furthermore, by neglecting certain types of molecular interactions such as charge penetration and charge transfer, most classical force fields must rely on, but do not always demonstrate, how cancellation of errors occurs among the remaining molecular interactions accounted for such as exchange repulsion, electrostatics, and polarization. In this work we present the first generation of the (many-body) MB-UCB force field that explicitly accounts for the decomposed molecular interactions commensurate with a variational energy decomposition analysis, including charge transfer, with force field design choices 1 arXiv:1905.07816v3 [physics.chem-ph] 27 Jul 2019 that reduce the computational expense of the MB-UCB potential while remaining accurate. We optimize parameters using only single water molecule and water cluster data up through pentamers, with no fitting to condensed phase data, and we demonstrate that high accuracy is maintained when the force field is subsequently validated against conformational energies of larger water cluster data sets, radial distribution functions of the liquid phase, and the temperature dependence of thermodynamic and transport water properties. We conclude that MB-UCB is comparable in performance to MB-Pol, but is less expensive and more transferable by eliminating the need to represent short-ranged interactions through large parameter fits to high order polynomials.

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.