Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

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.

Dehydrative cyclocondensation mechanisms of hydrogen thioperoxide and of alkanesulfenic acids

Structural features of hydrogen thioperoxide (oxadisulfane, H-S-O-H) and of alkanesulfenic acids (R-S-O-H; R = CH3, CH2CH3, CH2CH2CH3, CH(CH3)2, C(CH3)3, CF3, CCl3) and the mechanisms of their dehydrative cyclocondensation to the respective sulfinothioic acid (H-(S═O)-S-H) and alkyl alkanethiosulfinates (R-(S═O)-S-R) have been studied using coupled cluster theory with single and double and perturbative triple excitations [CCSD(T)] and quadratic configuration interaction with single and double and perturbative triple excitations [QCISD(T)] with the cc-pVDZ basis set and also using second-order Møller-Plesset perturbation theory (MP2) and the hybrid density functionals B3LYP, B3PW91, and PBE1PBE with the 6-311+G(d,p) basis set. The concerted cyclodehydration mechanisms include cyclic five-center transition states with relatively long distance sulfur-sulfur bonding interactions. Attractive and repulsive nonbonding interactions are predicted in the sulfenic acids, transition states, and thiosulfinates. In the alkyl alkanethiosulfinates attractive cyclic C-H----O═S nonbonding interactions are predicted. CCSD(T) and QCISD(T) predict similar values for the relative energies and CCSD(T) predicts the barrier to the cyclocondensation of H-S-O-H to sulfinothioic acid (H-(S═O)-S-H) to be 41.8 kcal/mol, and barriers in the range of 37.5 to 39.6 kcal/mol are predicted for the condensation of alkanesulfenic acids to alkyl alkanethiosulfinates.

Advances in methods and algorithms in a modern quantum chemistry program package

Advances in theory and algorithms for electronic structure calculations must be incorporated into program packages to enable them to become routinely used by the broader chemical community. This work reviews advances made over the past five years or so that constitute the major improvements contained in a new release of the Q-Chem quantum chemistry package, together with illustrative timings and applications. Specific developments discussed include fast methods for density functional theory calculations, linear scaling evaluation of energies, NMR chemical shifts and electric properties, fast auxiliary basis function methods for correlated energies and gradients, equation-of-motion coupled cluster methods for ground and excited states, geminal wavefunctions, embedding methods and techniques for exploring potential energy surfaces.