Search

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.

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 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 correlated Cr2 dimer, exploring zeolite-catalysed 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.

We report an implementation of periodic boundary conditions for mixed quantum mechanics/molecular mechanics (QM/MM) simulations, in which atomic partial charges are used to represent periodic images of the QM region. These charges are incorporated into the Fock matrix in a manner that preserves the variational nature of the self-consistent field procedure, and their interactions with the MM charges are summed using the conventional Ewald technique. To ensure that the procedure is stable in arbitrary basis sets, the atomic charges are derived by least-squares fit to the electrostatic potential generated by the QM region. We formulate and implement analytic energy gradients for the QM/MM-Ewald method and demonstrate that stable molecular dynamics simulations are thereby obtained. As a proof-of-concept application, we perform QM/MM simulations of a hydrated electron in bulk liquid water at the level of Hartree-Fock theory plus empirical dispersion. These simulations demonstrate that the "cavity model" of the aqueous electron, in which the spin density of the anionic defect is localized within an excluded volume in the liquid, is stable at room temperature on a time scale of at least several picoseconds. These results validate cavity-forming pseudopotential models of e−(aq) that have previously been derived from static-exchange Hartree-Fock calculations, and cast doubt upon whether non-cavity-forming pseudopotentials are faithful to the underlying Hartree-Fock calculation from which they were obtained. [ABSTRACT FROM AUTHOR]

An implementation of Ewald summation for use in mixed quantum mechanics/molecular mechanics (QM/MM) calculations is presented, which builds upon previous work by others that was limited to semi-empirical electronic structure for the QM region. Unlike previous work, our implementation describes the wave function's periodic images using 'ChElPG' atomic charges, which are determined by fitting to the QM electrostatic potential evaluated on a real-space grid. This implementation is stable even for large Gaussian basis sets with diffuse exponents, and is thus appropriate when the QM region is described by a correlated wave function. Derivatives of the ChElPG charges with respect to the QM density matrix are a potentially serious bottleneck in this approach, so we introduce a ChElPG algorithm based on atom-centered Lebedev grids. The ChElPG charges thus obtained exhibit good rotational invariance even for sparse grids, enabling significant cost savings. Detailed analysis of the optimal choice of user-selected Ewald parameters, as well as timing breakdowns, is presented. [ABSTRACT FROM AUTHOR]