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.

Comparison of the performance of exact-exchange-based density functional methods

How to describe nondynamic electron correlation is still a major challenge to density functional theory (DFT). Recent models designed particularly for this problem, such as Becke'05 (B05) and Perdew-Staroverov-Tao-Scuseria (PSTS) functionals employ the exact-exchange density, the efficient calculation of which is technically quite challenging. We have recently implemented self-consistently the B05 functional based on an efficient resolution-identity (RI) technique. In this study, we report a self-consistent RI implementation of the PSTS functional. In contrast to its original implementation, our version brings no limitation on the choice of the basis set. We have also implemented the Mori-Sanchez-Cohen-Yang-2 (MCY2) functional, another recent DFT method that includes full exact exchange. The performance of PSTS, B05, and MCY2 is validated on thermochemistry, reaction barriers, and dissociation energy curves, with an emphasis on nondynamic correlation effects in the discussion. All three methods perform rather well in general, B05 and MCY2 being on average somewhat better than PSTS. We include also results with other functionals that represent various aspects of the development in this field in recent years, including B3LYP, M06-HF, M06-2X, ωB97X, and TPSSh. The performance of the heavy-parameterized functionals M06-2X and ωB97X is on average better than that of B05, MCY2, and PSTS for standard thermodynamic properties and reactions, while the latter functionals do better in hydrogen abstraction reactions and dissociation processes. In particular, B05 is found to be the only functional that yields qualitatively correct dissociation curves for two-center symmetric radicals like He(2)(+). Finally, we compare the performance of all these functionals on a strongly correlated exemplary case system, the NO dimer. Only PSTS, B05, and MCY2 describe the system qualitatively correctly. Overall, this new type of functionals show good promise of overcoming some of the difficulties DFT encounters for systems with strong nondynamic correlation.

Reactions within fluorobenzene-ammonia heterocluster ions: experiment and theory

Reactions occurring within gas phase fluorobenenze-ammonia heterocluster cations (FC(6)H(5)-(NH(3))(n=1-4)) have been studied through the use of a triple quadrupole mass spectrometer as well as employing density functional theory (DFT). Collision induced dissociation (CID) experiments were conducted in which mass selected cluster ions are accelerated into a cell containing argon gas and the resulting products then subsequently mass analyzed. Two dominate reaction channels are observed. The first is simple evaporative loss of neutral ammonia from the cluster ion. The second involves a substitution reaction occurring within the cluster ion to form the aniline cation, C(6)H(5)NH(2)(+), where the reactivity was found to vary as a function of cluster size. DFT calculations have been performed to both help analyze the structure and the reactivity of these cluster ions. Pronounced differences in activation energies were found that provide an explanation for the observed variation of reactivity as a function of cluster size. An ad hoc model based upon the Arrhenius equation was developed to fit both the experimental collision energy dependence of the reaction and the observed lowering of the reaction barrier to aniline formation as a function of cluster size.

Reactions within p-difluorobenzene/methanol heterocluster ions: a detailed experimental and theoretical investigation

The reactivity of p-difluorobenzene/methanol cluster ions has been investigated by using triple quadrupole mass spectrometry and DFT calculations. The present study was performed in light of a recent investigation of p-difluorobenzene/methanol (P = F-C(6)H(4)-F and M = CH(3)OH) heterocluster ions where the solvent-catalyzed formation of p-fluoroanisole (A = CH(3)O-C(6)H(4)-F) was observed in P(M)(2)(+) clusters and not in PM(+) clusters. The results of our mass selected cluster ion study and theoretical calculations confirm that a single extra molecule of methanol can lower the reaction activation energy barrier in agreement with previous work for smaller clusters (PM(+) and P(M)(2)(+)). However, we also observe that P(M)(3)(+) and P(M)(4)(+) clusters undergo evaporative loss of neutral methanol to establish the P(M)(2)(+) cluster before reacting. P(M)(n>4)(+) clusters are capable of reacting through multiple pathways, in some cases generating a 1,4-dimethoxybenzene (B = CH(3)O-C(6)H(4)-OCH(3)) product via two separate substitution reactions within the same cluster ion. DFT calculations were employed to model the structures of the parent cluster ions, and transition state calculations were used to evaluate the activation energy for the p-fluoroanisole-forming substitution reaction. The calculations suggest that the reaction proceeds through a transition state containing a six-member hydrogen-bonded ring involving a reacting methanol and a second methanol that significantly lowers the activation energy.

Combined QM/MM calculations of active-site vibrations in binding process of P450cam to putidaredoxin

Combined QM/MM calculations of the active-site of cytochrome P450cam have been performed before and after the binding of P450cam to putidaredoxin. The calculations were carried out for both a 5-coordinated and a 6-coordinated active-site of cytochrome P450cam, with either a water molecule or a carbon monoxide molecule as a 6th distal ligand. An experimentally observed increase in the Fe-S stretching frequency that occurs after cytochrome P450cam binds to putidaredoxin, has been reproduced in our study. Experimentally observed changes in the Fe-C and C-O vibration frequencies that occur after binding of both proteins, have also been reproduced in our study. The computed increase of the Fe-S and Fe-C stretching frequencies is correlated with a corresponding decrease of the Fe-S and Fe-C interatomic distances. According to our calculations, for the active-site with carbon monoxide in the triplet electronic state, the binding process increases the spin densities on the iron and sulfur atoms, which changes the Fe-C and C-O stretching frequencies in opposite directions, in agreement with experimental data.

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.

Lennard-Jones parameters for the combined QM/MM method using the B3LYP/6-31G*/AMBER potential

A combined DFT quantum mechanical and AMBER molecular mechanical potential (QM/MM) is presented for use in molecular modeling and molecular simulations of large biological systems. In our approach we evaluate Lennard-Jones parameters describing the interaction between the quantum mechanical (QM) part of a system, which is described at the B3LYP/6-31+G* level of theory, and the molecular mechanical (MM) part of the system, described by the AMBER force field. The Lennard-Jones parameters for this potential are obtained by calculating hydrogen bond energies and hydrogen bond geometries for a large set of bimolecular systems, in which one hydrogen bond monomer is described quantum mechanically and the other is treated molecular mechanically. We have investigated more than 100 different bimolecular systems, finding very good agreement between hydrogen bond energies and geometries obtained from the combined QM/MM calculations and results obtained at the QM level of theory, especially with respect to geometry. Therefore, based on the Lennard-Jones parameters obtained in our study, we anticipate that the B3LYP/6-31+G*/AMBER potential will be a precise tool to explore intermolecular interactions inside a protein environment.