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.

Projected Hartree-Fock (PHF) theory has a long history in quantum chemistry. PHF is here understood as the variational determination of an N-electron broken symmetry Slater determinant that minimizes the energy of a projected state with the correct quantum numbers. The method was actively pursued for several decades but seems to have been abandoned. We here derive and implement a "variation after projection" PHF theory using techniques different from those previously employed in quantum chemistry. Our PHF methodology has modest mean-field computational cost, yields relatively simple expressions, can be applied to both collinear and non-collinear spin cases, and can be used in conjunction with deliberate symmetry breaking and restoration of other molecular symmetries like complex conjugation and point group. We present several benchmark applications to dissociation curves and singlet-triplet energy splittings, showing that the resulting PHF wavefunctions are of high quality multireference character. We also provide numerical evidence that in the thermodynamic limit, the energy in PHF is not lower than that of broken-symmetry HF, a simple consequence of the lack of size consistency and extensivity of PHF.

Water splitting by artificial catalysts is a critical process in the production of hydrogen gas as an alternative fuel. In this paper, we examine the essential role of theoretical calculations, with particular focus on density functional theory (DFT), in understanding the water-splitting reaction on these catalysts. First, we present an overview of DFT thermochemical calculations on water-splitting catalysts, addressing how these calculations are adapted to condensed phases and room temperature. We show how DFT-derived chemical descriptors of reactivity can be surprisingly good estimators for reactive trends in water-splitting catalysts. Using this concept, we recover trends for bulk catalysts using simple model complexes for at least the first-row transition-metal oxides. Then, using the CoPi cobalt oxide catalyst as a case study, we examine the usefulness of simulation for predicting the kinetics of water splitting. We demonstrate that the appropriate treatment of solvent effects is critical for computing accurate redox potentials with DFT, which, in turn, determine the rate-limiting steps and electrochemical overpotentials. Finally, we examine the ability of DFT to predict mechanism, using ruthenium complexes as a focal point for discussion. Our discussion is intended to provide an overview of the current strengths and weaknesses of the state-of-the-art DFT methodologies for condensed-phase molecular simulation involving transition metals and also to guide future experiments and computations toward the understanding and development of novel water-splitting catalysts.

Restricted open-shell Kohn-Sham (ROKS) theory provides a powerful computational tool for calculating singlet excited state energies and dynamics. However, the possibility of multiple solutions to the ROKS equations - with the associated difficulty of automatically selecting the physically meaningful solution - limits its usefulness for intensive applications such as long-time Born-Oppenheimer molecular dynamics. We present an implementation of ROKS for excited states which prescribes the physically correct solution from an overlap criterion and guarantees that this solution is stationary, allowing for straightforward evaluation of nuclear gradients. The method is used to benchmark ROKS for vertical excitation energies of small and large organic dyes and for the calculation of Stokes shifts. With common density functional approximations, ROKS vertical excitation energies, and Stokes shifts show similar accuracy to those from time-dependent density functional theory and Δ-self-consistent-field approaches. Advantages of the ROKS approach for excited state structure and molecular dynamics are discussed.