A summary of the technical advances that are incorporated in the fourth major release of the Q-Chem quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and openshell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Møller-Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Some specific chemical examples include calculations on the strongly correlated Cr 2 dimer, exploring zeolitecatalysed ethane dehydrogenation, energy decomposition analysis of a charged ter-molecular complex arising from glycerol photoionisation, and natural transition orbitals for a Frenkel exciton state in a nine-unit model of a self-assembling nanotube.Keywords quantum chemistry, software, electronic structure theory, density functional theory, electron correlation, computational modelling, Q-Chem Disciplines Chemistry CommentsThis article is from Molecular Physics: An International Journal at the Interface Between Chemistry and Physics 113 (2015): 184, doi:10.1080/00268976.2014. RightsWorks produced by employees of the U.S. Government as part of their official duties are not copyrighted within the U.S. The content of this document is not copyrighted. Authors 185A summary of the technical advances that are incorporated in the fourth major release of the Q-CHEM quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and open-shell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Møller-Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Some specific chemical examples include calculations on the strongly corre...

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.

Previously, we reported an electron-water pseudopotential designed to be used in conjunction with a polarizable water model, in order to describe the hydrated electron [L. D. Jacobson et al., J. Chem. Phys. 130, 124115 (2009)]. Subsequently, we found this model to be inadequate for the aqueous electron in bulk water, and here we report a reparametrization of the model. Unlike the previous model, the current version is not fit directly to any observables; rather, we use an ab initio exchange-correlation potential, along with a repulsive potential that is fit to reproduce the density maximum of the excess electron's wave function within the static-exchange approximation. The new parametrization performs at least as well as the previous model, as compared to ab initio benchmarks for (H(2)O)(n) (-) clusters, and also predicts reasonable values for the diffusion coefficient, radius of gyration, and absorption maximum of the bulk species. The new model predicts a vertical electron binding energy of 3.7 eV in bulk water, which is 1.4 eV smaller than the value obtained using nonpolarizable models; the difference represents the solvent's electronic reorganization energy following electron detachment. We find that the electron's first solvation shell is quite loose, which may be responsible for the electron's large, positive entropy of hydration. Many-body polarization alters the electronic absorption line shape in a qualitative way, giving rise to a high-energy tail that is observed experimentally but is absent in previous simulations. In our model, this feature arises from spatially diffuse excited states that are bound only by electronic reorganization (i.e., solvent polarization) following electronic excitation.

We report a fragment-based electronic structure method, intended for the study of clusters and molecular liquids, that incorporates electronic polarization (induction) in a self-consistent fashion but treats intermolecular exchange and dispersion interactions perturbatively, as post-self-consistent field corrections, using a form of pairwise symmetry-adapted perturbation theory. The computational cost of the method scales quadratically as a function of the number of fragments (monomers), but could be made to scale linearly by exploiting distance-dependent thresholds. Extensive benchmark calculations are reported using the S22 database of high-level ab initio binding energies for dimers, and we find that average errors can be reduced to <1 kcal/mol with a suitable choice of basis set. Comparison to ab initio benchmarks for water clusters as large as (H 2 O) 20 demonstrates that the method recovers 90% of the binding energy in these systems, at a tiny fraction of the computational cost. As such, this approach represents a promising path toward accurate, systematically improvable, and parameter-free simulation of molecular liquids.

A method that we have recently introduced for rapid computation of intermolecular interaction energies is reformulated and subjected to further tests. The method employs monomer-based self-consistent field calculations with an electrostatic embedding designed to capture many-body polarization (the "XPol" procedure), augmented by pairwise symmetry-adapted perturbation theory (SAPT) to capture dispersion and exchange interactions along with any remaining induction effects. A rigorous derivation of the XPol+SAPT methodology is presented here, which demonstrates that the method is systematically improvable, and moreover introduces some additional intermolecular interactions as compared to the more heuristic derivation that was presented previously. Applications to various non-covalent complexes and clusters are presented, including geometry optimizations and one-dimensional potential energy scans. The performance of the XPol+SAPT methodology in its present form (based on second-order intermolecular perturbation theory and neglecting intramolecular electron correlation) is qualitatively acceptable across a wide variety of systems-and quantitatively quite good in certain cases-but the quality of the results is rather sensitive to the choice of one-particle basis set. Basis sets that work well for dispersion-bound systems offer less-than-optimal performance for clusters dominated by induction and electrostatic interactions, and vice versa. A compromise basis set is identified that affords good results for both induction and dispersion interactions, although this favorable performance ultimately relies on error cancellation, as in traditional low-order SAPT. Suggestions for future improvements to the methodology are discussed.