Search

The increasing availability of GPUs for scientific computing has prompted interest in accelerating quantum chemical calculations through their use. The complexity of integral kernels for high angular momentum basis functions however often limits the utility of GPU implementations with large basis sets or for metal containing systems. In this work, we report implementation of $f$ function support in the GPU-accelerated TeraChem software package through the development of efficient kernels for the evaluation of Hamiltonian integrals. The high efficiency of the resulting code is demonstrated through density functional theory (DFT) calculations on increasingly large organic molecules and transition metal complexes, as well as coupled cluster singles and doubles (CCSD) calculations on water clusters. Preliminary investigations into Ni(I) catalysis with DFT and the photochemistry of MnH(CH$_3$) with complete active space self-consistent field (CASSCF) are also carried out. Overall, our GPU-accelerated software appears to be well-suited for fast simulation of large transition metal containing systems, as well as organic molecules., Comment: Added additional data

The fate of thymine upon excitation by ultraviolet radiation has been the subject of intense debate over the past three decades. Today, it is widely believed that its ultrafast excited state decay stems from a radiationless transition from the bright ${\pi}{\pi}^*$ state to a dark $n{\pi}^*$ state. However, conflicting theoretical predictions have made the experimental data difficult to interpret. Here we simulate the ultrafast dynamics in thymine at the highest level of theory to date, performing wavepacket dynamics with a new coupled cluster method. Our simulation confirms an ultrafast ${\pi}{\pi}^*$ to $n{\pi}^*$ transition (${\tau} = 41 \pm 14$ fs). Furthermore, the predicted oxygen-edge X-ray absorption spectra agree quantitatively with the experimental results. Our simulation also predicts an as-yet uncharacterized photochemical pathway: a ${\pi}{\sigma}^*$ channel that leads to hydrogen dissociation at one of the two N-H bonds in thymine. Similar behavior has been identified in other heteroaromatic compounds, including adenine, and several authors have speculated that a similar pathway may exist in thymine. However, this was never confirmed theoretically or experimentally. This prediction calls for renewed efforts to experimentally identify or exclude the presence of this channel., Comment: 42 pages, 23 figures

Simulations of photochemical reaction dynamics have been a challenge to the theoretical chemistry community for some time. In an effort to determine the predictive character of current approaches, we predict the results of an upcoming ultrafast diffraction experiment on the photodynamics of cyclobutanone after excitation to the lowest lying Rydberg state (S$_2$). A picosecond of nonadiabatic dynamics is described with ab initio multiple spawning. We use both time dependent density functional theory and equation-of-motion coupled cluster for the underlying electronic structure theory. We find that the lifetime of the S$_2$ state is more than a picosecond (with both TDDFT and EOM-CCSD). The predicted UED spectrum exhibits numerous structural features, but weak time dependence over the course of the simulations.

This study implements the full multicomponent third-order (MP3) and fourth-order (MP4) many-body perturbation theory methods for the first time. Previous multicomponent studies have only implemented a subset of the full contributions and the present implementation is the first multicomponent many-body method to include any connected triples contribution to the electron-proton correlation energy. The multicomponent MP3 method is shown to be comparable in accuracy to the multicomponent coupled-cluster doubles method for the calculation of proton affinities, while the multicomponent MP4 method is of similar accuracy as the multicomponent coupled-cluster singles and doubles (CCSD) method. From the results in this study, it is hypothesized that the relative accuracy of multicomponent methods is more similar to their single-component counterparts than previously assumed and that the previous explanation for the better performance of multicomponent CCSD compared to other ab initio multicomponent methods may be incorrect. It is demonstrated that for multicomponent MP4, the fourth-order triple-excitation contributions can be split into electron-electron and electron-proton contributions and the electron-electron contributions ignored with very little loss of accuracy of protonic properties., Comment: Minor revisions and removal of multicomponent CCSD material

Multicomponent methods seek to treat select nuclei, typically protons, fully quantum mechanically and equivalent to the electrons of a chemical system. In such methods, it is well known that due to the neglect of electron-proton correlation, a Hartree-Fock (HF) description of the electron-proton interaction catastrophically fails leading to qualitatively incorrect protonic properties. In single-component quantum chemistry, the qualitative failure of HF is normally indicative of the need for multireference methods such as complete active space self-consistent field (CASSCF). While a multicomponent CASSCF method was implemented nearly twenty years ago, it is only able to perform calculations with very small active spaces (~105 multicomponent configurations). Therefore, in order to extend the realm of applicability of the multicomponent CASSCF method, this study derives and implements a new two-step multicomponent CASSCF method that uses multicomponent heat-bath configuration interaction for the configuration interaction step, enabling calculations with very large active spaces (up to 16 electrons in 48 orbitals). We find that large electronic active spaces are needed to obtain qualitatively accurate protonic densities for the HCN and FHF- molecules. Additionally, the multicomponent CASSCF method implemented here should have further applications for double-well protonic potentials and systems that are inherently electronically multireference., Comment: 31 pages, Minor Revisions