Devin A. Matthews

Calculation of Vibrational Transition Frequencies and Intensities in Water Dimer: Comparison of Different Vibrational Approaches

We have calculated frequencies and intensities of fundamental and overtone vibrational transitions in water and water dimer with use of different vibrational methods. We have compared results obtained with correlation-corrected vibrational self-consistent-field theory and vibrational second-order perturbation theory both using normal modes and finally with a harmonically coupled anharmonic oscillator local mode model including OH-stretching and HOH-bending local modes. The coupled cluster with singles, doubles, and perturbative triples ab initio method with augmented correlation-consistent triple-zeta Dunning and atomic natural orbital basis sets has been used to obtain the necessary potential energy and dipole moment surfaces. We identify the strengths and weaknesses of these different vibrational approaches and compare our results to the available experimental results.

Ab initio determination of the crystalline benzene lattice energy to sub-kilojoule/mole accuracy

Working out how to pack benzene in silico Many organic compounds crystallize in several different energetically similar packing arrangements, or polymorphs. This complicates processes such as drug formulation that rely on reproducible crystallization. Yang et al. have now achieved the long-standing goal of calculating a crystal packing arrangement from first principles to an accuracy that can distinguish polymorphs (see the Perspective by Price). They used benzene as a prototypical test case and applied quantum chemical methods that improve estimates of multibody interactions. The results bode well for future applications of theory to optimization of crystallization protocols. Science, this issue p. 640; see also p. 619 Theoretical calculations of molecular packing in crystals attain sufficient accuracy to distinguish polymorphs. [Also see Perspective by Price] Computation of lattice energies to an accuracy sufficient to distinguish polymorphs is a fundamental bottleneck in crystal structure prediction. For the lattice energy of the prototypical benzene crystal, we combined the quantum chemical advances of the last decade to attain sub-kilojoule per mole accuracy, an order-of-magnitude improvement in certainty over prior calculations that necessitates revision of the experimental extrapolation to 0 kelvin. Our computations reveal the nature of binding by improving on previously inaccessible or inaccurate multibody and many-electron contributions and provide revised estimates of the effects of temperature, vibrations, and relaxation. Our demonstration raises prospects for definitive first-principles resolution of competing polymorphs in molecular crystal structure prediction.

A massively parallel tensor contraction framework for coupled-cluster computations

Precise calculation of molecular electronic wavefunctions by methods such as coupled-cluster requires the computation of tensor contractions, the cost of which has polynomial computational scaling with respect to the system and basis set sizes. Each contraction may be executed via matrix multiplication on a properly ordered and structured tensor. However, data transpositions are often needed to reorder the tensors for each contraction. Writing and optimizing distributed-memory kernels for each transposition and contraction is tedious since the number of contractions scales combinatorially with the number of tensor indices. We present a distributed-memory numerical library (Cyclops Tensor Framework (CTF)) that automatically manages tensor blocking and redistribution to perform any user-specified contractions. CTF serves as the distributed-memory contraction engine in Aquarius, a new program designed for high-accuracy and massively-parallel quantum chemical computations. Aquarius implements a range of coupled-cluster and related methods such as CCSD and CCSDT by writing the equations on top of a C++ templated domain-specific language. This DSL calls CTF directly to manage the data and perform the contractions. Our CCSD and CCSDT implementations achieve high parallel scalability on the BlueGene/Q and Cray XC30 supercomputer architectures showing that accurate electronic structure calculations can be effectively carried out on top of general distributed-memory tensor primitives. We introduce Cyclops Tensor Framework (CTF), a distributed-memory library for tensor contractions.CTF is able to perform tensor decomposition, redistribution, and contraction at runtime.CTF enables the expression of massively-parallel coupled-cluster methods via a concise tensor contraction interface.The quantum chemistry software suite Aquarius employs CTF to execute two coupled-cluster methods: CCSD and CCSDT.The Aquarius CCSD and CCSDT codes scale well on BlueGene/Q and Cray XC30, comparing favorably to NWChem.

Calculated stretching overtone levels and Darling–Dennison resonances in water: a triumph of simple theoretical approaches

The coupled-cluster singles and doubles treatment with a perturbative treatment of triple excitations known as CCSD(T) has been used in conjunction with a hierarchy of atomic natural orbital basis sets to study stretching levels in water up to 2 eV above the zero-point level. Agreement with experiment obtained with perturbation theory augmented by a simple treatment of resonances suggested by Lehmann is quite remarkable in spite of the well-known and strong Darling–Dennison resonances in this system. With the largest basis sets, deviation from experiment is ca. 20 cm−1 at energies in the 15,000–16,000 cm−1 range.

Quantitative analysis of Fermi resonances by harmonic derivatives of perturbation theory corrections

Vibrational perturbation theory has proven to be a highly accurate and efficient method for extending the harmonic approximation in the treatment of polyatomic molecular vibrations. Unfortunately, accidental near-degeneracies of the harmonic vibrational levels can lead to resonance and a breakdown of the perturbation approximation. These resonances can be resolved by the diagonalization of a small effective Hamiltonian derived from either of the usual Rayleigh–Schrödinger or van Vleck perturbation theories. However, the proper choice of states for inclusion in the effective Hamiltonian is crucial to the accuracy of the results, and is not often clearly evident. It is proposed that the analytical partial derivatives of the anharmonic vibrational correction with respect to the various harmonic frequencies, called ‘Harmonic Derivatives’ in this work, can be used as a tool to quantitatively assess the existence and strength of first-order, or Fermi, resonances. These derivatives are shown to concisely and clearly reflect the quality of the perturbation approximation and the effect of its breakdown on the computed vibrational levels.

Stabilization of the Simplest Criegee Intermediate from the Reaction between Ozone and Ethylene: A High-Level Quantum Chemical and Kinetic Analysis of Ozonolysis

The fraction of the collisionally stabilized Criegee species CH2OO produced from the ozonolysis of ethylene is calculated using a two-dimensional (E, J)-grained master equation technique and semiclassical transition-state theory based on the potential energy surface obtained from high-accuracy quantum chemical calculations. Our calculated yield of 42 ± 6% for the stabilized CH2OO agrees well, within experimental error, with available (indirect) experimental results. Inclusion of angular momentum in the master equation is found to play an essential role in bringing the theoretical results into agreement with the experiment. Additionally, yields of HO and HO2 radical products are predicted to be 13 ± 6% and 17 ± 6%, respectively. In the kinetic simulation, the HO radical product is produced mostly from the stepwise decomposition mechanism of primary ozonide rather than from dissociation of hot CH2OO.

Non-orthogonal spin-adaptation of coupled cluster methods: A new implementation of methods including quadruple excitations

The theory of non-orthogonal spin-adaptation for closed-shell molecular systems is applied to coupled cluster methods with quadruple excitations (CCSDTQ). Calculations at this level of detail are of critical importance in describing the properties of molecular systems to an accuracy which can meet or exceed modern experimental techniques. Such calculations are of significant (and growing) importance in such fields as thermodynamics, kinetics, and atomic and molecular spectroscopies. With respect to the implementation of CCSDTQ and related methods, we show that there are significant advantages to non-orthogonal spin-adaption with respect to simplification and factorization of the working equations and to creating an efficient implementation. The resulting algorithm is implemented in the CFOUR program suite for CCSDT, CCSDTQ, and various approximate methods (CCSD(T), CC3, CCSDT-n, and CCSDT(Q)).

Spectral functions of the uniform electron gas via coupled-cluster theory and comparison to the G W and related approximations

We use ab initio coupled-cluster theory to compute the spectral function of the uniform electron gas at a Wigner-Seitz radius of r_s=4. The coupled-cluster approximations we employ go significantly beyond the diagrammatic content of state-of-the-art GW theory. We compare our calculations extensively to GW and GW-plus-cumulant theory, illustrating the strengths and weaknesses of these methods in capturing the quasiparticle and satellite features of the electron gas. Our accurate calculations further allow us to address the long-standing debate over the occupied bandwidth of metallic sodium. Our findings indicate that the future application of coupled-cluster theory to condensed phase material spectra is highly promising.

Assessment of the accuracy of coupled cluster perturbation theory for open-shell systems. II. Quadruples expansions.

The accuracy at which total energies of open-shell atoms and organic radicals may be calculated is assessed for selected coupled cluster perturbative triples expansions, all of which augment the coupled cluster singles and doubles (CCSD) energy by a non-iterative correction for the effect of triple excitations. Namely, the second- through sixth-order models of the recently proposed CCSD(T-n) triples series [J. J. Eriksen et al., J. Chem. Phys. 140, 064108 (2014)] are compared to the acclaimed CCSD(T) model for both unrestricted as well as restricted open-shell Hartree-Fock (UHF/ROHF) reference determinants. By comparing UHF- and ROHF-based statistical results for a test set of 18 modest-sized open-shell species with comparable RHF-based results, no behavioral differences are observed for the higher-order models of the CCSD(T-n) series in their correlated descriptions of closed- and open-shell species. In particular, we find that the convergence rate throughout the series towards the coupled cluster singles, doubles, and triples (CCSDT) solution is identical for the two cases. For the CCSD(T) model, on the other hand, not only its numerical consistency, but also its established, yet fortuitous cancellation of errors breaks down in the transition from closed- to open-shell systems. The higher-order CCSD(T-n) models (orders n > 3) thus offer a consistent and significant improvement in accuracy relative to CCSDT over the CCSD(T) model, equally for RHF, UHF, and ROHF reference determinants, albeit at an increased computational cost.

Accelerating the convergence of higher-order coupled cluster methods

The problem of the generally inferior convergence behavior of higher-order coupled cluster methods, such as CCSDT and CCSDTQ, compared to CCSD is analyzed in terms of Møller-Plesset perturbation theory. A new structure for the CCSDT and CCSDTQ equations (and various approximations of these) is proposed which reorders contributions between the various cluster amplitudes and emphasizes lower-order corrections to the energy at each iteration. Numerical testing of the proposed method compared to the widely used direct inversion in the iterative subspace convergence acceleration technique shows significant improvement in the rate of convergence and total time-to-solution, especially for methods including quadruple excitations.