Automatic State Interaction with Large Localized Active Spaces for Multimetallic Systems

The localized active space self-consistent field method factorizes a complete active space wave function into an antisymmetrized product of localized active space wave function fragments. Correlation between fragments is then reintroduced through localized active space state interaction (LASSI), in which the Hamiltonian is diagonalized in a model space of LAS states. However, the optimal procedure for defining the LAS fragments and LASSI model space is unknown. We here present an automated framework to explore systematically convergent sets of model spaces, which we call LASSI[r, q]. This method requires the user to select only r, the number of electron hops from one fragment to another, and q, the number of fragment basis functions per Hilbert space, which converges to CASCI in the limit of r, q → ∞. Numerical tests of this method on the trimetal oxo-centered complexes [Fe(III)Al(III)Fe(II)(μ3-O)(HCOO)6] and [Fe(III)2Fe(II)(μ3-O)(HCOO)6] show efficient convergence to the CASCI limit with 4-10 orders of magnitude fewer states than CASCI.

Analytic Nuclear Gradients for Complete Active Space Linearized Pair-Density Functional Theory

Accurately modeling photochemical reactions is difficult due to the presence of conical intersections and locally avoided crossings, as well as the inherently multiconfigurational character of excited states. As such, one needs a multistate method that incorporates state interaction in order to accurately model the potential energy surface at all nuclear coordinates. The recently developed linearized pair-density functional theory (L-PDFT) is a multistate extension of multiconfiguration PDFT, and it has been shown to be a cost-effective post-MCSCF method (as compared to more traditional and expensive multireference many-body perturbation methods or multireference configuration interaction methods) that can accurately model potential energy surfaces in regions of strong nuclear-electronic coupling in addition to accurately predicting Franck-Condon vertical excitations. In this paper, we report the derivation of analytic gradients for L-PDFT and their implementation in the PySCF-forge software, and we illustrate the utility of these gradients for predicting ground- and excited-state equilibrium geometries and adiabatic excitation energies for formaldehyde, s-trans-butadiene, phenol, and cytosine.

State Preparation in Quantum Algorithms for Fragment-Based Quantum Chemistry

State preparation for quantum algorithms is crucial for achieving high accuracy in quantum chemistry and competing with classical algorithms. The localized active space-unitary coupled cluster (LAS-UCC) algorithm iteratively loads a fragment-based multireference wave function onto a quantum computer. In this study, we compare two state preparation methods, quantum phase estimation (QPE) and direct initialization (DI), for each fragment. We test the two state preparation methods on three systems, ranging from a model system, a set of interacting hydrogen molecules, to more realistic chemical problems, like the C-C double bond breaking in transbutadiene and the spin ladder in a bimetallic system. We analyze the impact of QPE parameters, such as the number of ancilla qubits and Trotter steps, on the prepared state. We find a trade-off between the methods, where DI requires fewer resources for smaller fragments, while QPE is more efficient for larger fragments. Our resource estimates highlight the benefits of system fragmentation in state preparation for subsequent quantum chemical calculations. These findings have broad applications for preparing multireference quantum chemical wave functions on quantum circuits that can be used for realistic chemical applications.

Minimum-Energy Conical Intersections by Compressed Multistate Pair-Density Functional Theory

Compressed multistate pair-density functional theory (CMS-PDFT) is a multistate version of multiconfiguration pair-density functional theory that can capture the correct topology of coupled potential energy surfaces (PESs) around conical intersections. In this work, we develop interstate coupling vectors (ISCs) for CMS-PDFT in the OpenMolcas and PySCF/mrh electronic structure packages. Yet, the main focus of this work is using ISCs to calculate minimum-energy conical intersections (MECIs) by CMS-PDFT. This is performed using the projected constrained optimization method in OpenMolcas, which uses ISCs to restrain the iterations to the conical intersection seam. We optimize the S1/S0 MECIs for ethylene, butadiene, and benzene and show that CMS-PDFT gives smooth PESs in the vicinities of the MECIs. Furthermore, the CMS-PDFT MECIs are in good agreement with the MECI calculated by the more expensive XMS-CASPT2 method.

The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

A Perspective on Sustainable Computational Chemistry Software Development and Integration

The power of quantum chemistry to predict the ground and excited state properties of complex chemical systems has driven the development of computational quantum chemistry software, integrating advances in theory, applied mathematics, and computer science. The emergence of new computational paradigms associated with exascale technologies also poses significant challenges that require a flexible forward strategy to take full advantage of existing and forthcoming computational resources. In this context, the sustainability and interoperability of computational chemistry software development are among the most pressing issues. In this perspective, we discuss software infrastructure needs and investments with an eye to fully utilize exascale resources and provide unique computational tools for next-generation science problems and scientific discoveries.

Understanding Antiferromagnetic and Ligand Field Effects on Spin Crossover in a Triple-Decker Dimeric Cr(II) Complex

Two possible explanations for the temperature dependence of spin-crossover (SCO) behavior in the dimeric triple-decker Cr(II) complex ([(η5-C5Me5)Cr(μ2:η5-P5)Cr(η5-C5Me5)]+) have been offered. One invokes variations in antiferromagnetic interactions between the two Cr(II) ions, whereas the other posits the development of a strong ligand-field effect favoring the low-spin ground state. We perform multireference electronic structure calculations based on the multiconfiguration pair-density functional theory to resolve these effects. We find quintet, triplet, and singlet electronic ground states, respectively, for the experimental geometries at high, intermediate, and low temperatures. The ground-state transition from quintet to triplet at an intermediate temperature derives from increased antiferromagnetic interactions between the two Cr(II) ions. By contrast, the ground-state transition from triplet to singlet at low temperature can be attributed to increased ligand-field effects, which dominate with continued variations in antiferromagnetic coupling. This study provides quantitative detail for the degree to which these two effects can act in concert for the observed SCO behavior in this complex and others subject to temperature-dependent variations in geometry.

Optical Properties of Neutral F Centers in Bulk MgO with Density Matrix Embedding

The optical spectra of neutral oxygen vacancies (F0 centers) in the bulk MgO lattice are investigated using density matrix embedding theory. The impurity Hamiltonian is solved with the complete active space self-consistent field and second-order n-electron valence state perturbation theory (NEVPT2-DMET) multireference methods. To estimate defect-localized vertical excitation energies at the nonembedding and thermodynamic limits, a double extrapolation scheme is employed. The extrapolated NEVPT2-DMET vertical excitation energy value of 5.24 eV agrees well with the experimental absorption maxima at 5.03 eV, whereas the excitation energy value of 2.89 eV at the relaxed triplet defect-localized state geometry overestimates the experimental emission at 2.4 eV by only nearly 0.5 eV, indicating the involvement of the triplet-singlet decay pathway.

Density Matrix Embedding Using Multiconfiguration Pair-Density Functional Theory

We present a quantum embedding method for ground and excited states of extended systems that uses multiconfiguration pair-density functional theory (MC-PDFT) with densities provided by periodic density matrix embedding theory (pDMET). We compute local excitations in oxygen mono- and divacancies on a magnesium oxide (100) surface and find absolute deviations within 0.05 eV between pDMET using the MC-PDFT, denoted as pDME-PDFT, and the more expensive, nonembedded MC-PDFT approach. We further use pDME-PDFT to calculate local excitations in larger supercells for the monovacancy defect, for which the use of nonembedded MC-PDFT is prohibitively costly.

Linearized Pair-Density Functional Theory

Multiconfiguration pair-density functional theory (MC-PDFT) is a post-SCF multireference method that has been successful at computing ground- and excited-state energies. However, MC-PDFT is a single-state method in which the final MC-PDFT energies do not come from diagonalization of a model-space Hamiltonian matrix, and this can lead to inaccurate topologies of potential energy surfaces near locally avoided crossings and conical intersections. Therefore, in order to perform physically correct ab initio molecular dynamics with electronically excited states or to treat Jahn-Teller instabilities, it is necessary to develop a PDFT method that recovers the correct topology throughout the entire nuclear configuration space. Here we construct an effective Hamiltonian operator, called the linearized PDFT (L-PDFT) Hamiltonian, by expanding the MC-PDFT energy expression to first order in a Taylor series of the wave function density. Diagonalization of the L-PDFT Hamiltonian gives the correct potential energy surface topology near conical intersections and locally avoided crossings for a variety of challenging cases including phenol, methylamine, and the spiro cation. Furthermore, L-PDFT outperforms MC-PDFT and previous multistate PDFT methods for predicting vertical excitations from a variety of representative organic chromophores.

Local Excitations of a Charged Nitrogen Vacancy in Diamond with Multireference Density Matrix Embedding Theory

We investigate the negatively charged nitrogen-vacancy center in diamond using periodic density matrix embedding theory (pDMET). To describe the strongly correlated excited states of this system, the complete active space self-consistent field (CASSCF) followed by n-electron valence state second-order perturbation theory (NEVPT2) was used as the impurity solver. Since the NEVPT2-DMET energies show a linear dependence on the inverse of the size of the embedding subspace, we performed an extrapolation of the excitation energies to the nonembedding limit using a linear regression. The extrapolated NEVPT2-DMET first triplet-triplet excitation energy is 2.31 eV and that for the optically inactive singlet-singlet transition is 1.02 eV, both in agreement with the experimentally observed vertical excitation energies of ∼2.18 eV and ∼1.26 eV, respectively. This is the first application of pDMET to a charged periodic system and the first investigation of the NV^- defect using NEVPT2 for periodic supercell models.

Localized Quantum Chemistry on Quantum Computers

Quantum chemistry calculations of large, strongly correlated systems are typically limited by the computation cost that scales exponentially with the size of the system. Quantum algorithms, designed specifically for quantum computers, can alleviate this, but the resources required are still too large for today's quantum devices. Here, we present a quantum algorithm that combines a localization of multireference wave functions of chemical systems with quantum phase estimation (QPE) and variational unitary coupled cluster singles and doubles (UCCSD) to compute their ground-state energy. Our algorithm, termed "local active space unitary coupled cluster" (LAS-UCC), scales linearly with the system size for certain geometries, providing a polynomial reduction in the total number of gates compared with QPE, while providing accuracy above that of the variational quantum eigensolver using the UCCSD ansatz and also above that of the classical local active space self-consistent field. The accuracy of LAS-UCC is demonstrated by dissociating (H₂)₂ into two H₂ molecules and by breaking the two double bonds in trans-butadiene, and resource estimates are provided for linear chains of up to 20 H₂ molecules.

Localized Active Space-State Interaction: a Multireference Method for Chemical Insight

Multireference electronic structure methods, like the complete active space (CAS) self-consistent field model, have long been used to characterize chemically interesting processes. Important work has been done in recent years to develop modifications having a lower computational cost than CAS, but typically these methods offer no more chemical insight than that from the CAS solution being approximated. In this paper, we present the localized active space-state interaction (LASSI) method that can be used not only to lower the intrinsic cost of the multireference calculation but also to improve interpretability. The localized active space (LAS) approach utilizes the local nature of the electron-electron correlation to express a composite wave function as an antisymmetrized product of unentangled wave functions in local active subspaces. LASSI then uses these LAS states as a basis from which to express complete molecular wave functions. This not only makes the molecular wave function more compact but also permits flexibility in choosing those states to be included in the basis. Such selective inclusion of states translates to the selective inclusion of specific types of interactions, thereby allowing a quantitative analysis of these interactions. We demonstrate the use of LASSI to study charge migration and spin-flip excitations in multireference organic molecules. We also compute the J coupling parameter for a bimetallic compound using various LAS bases to construct the Hamiltonian to provide insights into the coupling mechanism.

Large-Scale Benchmarking of Multireference Vertical-Excitation Calculations via Automated Active-Space Selection

We have calculated state-averaged complete-active-space self-consistent-field (SA-CASSCF), multiconfiguration pair-density functional theory (MC-PDFT), hybrid MC-PDFT (HMC-PDFT), and n-electron valence state second-order perturbation theory (NEVPT2) excitation energies with the approximate pair coefficient (APC) automated active-space selection scheme for the QUESTDB benchmark database of 542 vertical excitation energies. We eliminated poor active spaces (20–40% of calculations) by applying a threshold to the SA-CASSCF absolute error. With the remaining calculations, we find that NEVPT2 performance is significantly impacted by the size of the basis set the wave functions are converged in, regardless of the quality of their description, which is a problem absent in MC-PDFT. Additionally, we find that HMC-PDFT is a significant improvement over MC-PDFT with the translated PBE (tPBE) density functional and that it performs about as well as NEVPT2 and second-order coupled cluster on a set of 373 excitations in the QUESTDB database. We optimized the percentage of SA-CASSCF energy to include in HMC-PDFT when using the tPBE on-top functional, and we find the 25% value used in tPBE0 to be optimal. This work is by far the largest benchmarking of MC-PDFT and HMC-PDFT to date, and the data produced in this work are useful as a validation of HMC-PDFT and of the APC active-space selection scheme. We have made all the wave functions produced in this work (orbitals and CI vectors) available to the public and encourage the community to utilize this data as a tool in the development of further multireference model chemistries.

Periodic Density Matrix Embedding for CO Adsorption on the MgO(001) Surface

The adsorption of simple gas molecules to metal oxide surfaces is a primary step in many heterogeneous catalysis applications. Quantum chemical modeling of these reactions is a challenge in terms of both cost and accuracy, and quantum-embedding methods are promising, especially for localized chemical phenomena. In this work, we employ density matrix embedding theory (DMET) for periodic systems to calculate the adsorption energy of CO to the MgO(001) surface. Using coupled-cluster theory with single and double excitations and second-order Møller-Plesset perturbation theory as quantum chemical solvers, we perform calculations with embedding clusters up to 266 electrons in 306 orbitals, with the largest embedding models agreeing to within 1.2 kcal/mol of the non-embedding references. Moreover, we present a memory-efficient procedure of storing and manipulating electron repulsion integrals in the embedding space within the framework of periodic DMET.

Electronic structure of strongly correlated systems: recent developments in multiconfiguration pair-density functional theory and multiconfiguration nonclassical-energy functional theory

Strong electron correlation plays an important role in transition-metal and heavy-metal chemistry, magnetic molecules, bond breaking, biradicals, excited states, and many functional materials, but it provides a significant challenge for modern electronic structure theory. The treatment of strongly correlated systems usually requires a multireference method to adequately describe spin densities and near-degeneracy correlation. However, quantitative computation of dynamic correlation with multireference wave functions is often difficult or impractical. Multiconfiguration pair-density functional theory (MC-PDFT) provides a way to blend multiconfiguration wave function theory and density functional theory to quantitatively treat both near-degeneracy correlation and dynamic correlation in strongly correlated systems; it is more affordable than multireference perturbation theory, multireference configuration interaction, or multireference coupled cluster theory and more accurate for many properties than Kohn–Sham density functional theory. This perspective article provides a brief introduction to strongly correlated systems and previously reviewed progress on MC-PDFT followed by a discussion of several recent developments and applications of MC-PDFT and related methods, including localized-active-space MC-PDFT, generalized active-space MC-PDFT, density-matrix-renormalization-group MC-PDFT, hybrid MC-PDFT, multistate MC-PDFT, spin–orbit coupling, analytic gradients, and dipole moments. We also review the more recently introduced multiconfiguration nonclassical-energy functional theory (MC-NEFT), which is like MC-PDFT but allows for other ingredients in the nonclassical-energy functional. We discuss two new kinds of MC-NEFT methods, namely multiconfiguration density coherence functional theory and machine-learned functionals.

This feature article overviews recent work on active spaces, matrix product reference states, treatment of quasidegeneracy, hybrid theory, density-coherence functionals, machine-learned functionals, spin–orbit coupling, gradients, and dipole moments.

Excited States of Crystalline Point Defects with Multireference Density Matrix Embedding Theory

Accurate and affordable methods to characterize the electronic structure of solids are important for targeted materials design. Embedding-based methods provide an appealing balance in the trade-off between cost and accuracy─particularly when studying localized phenomena. Here, we use the density matrix embedding theory (DMET) algorithm to study the electronic excitations in solid-state defects with a restricted open-shell Hartree-Fock (ROHF) bath and multireference impurity solvers, specifically, complete active space self-consistent field (CASSCF) and n-electron valence state second-order perturbation theory (NEVPT2). We apply the method to investigate the electronic excitations in an oxygen vacancy (OV) on a MgO(100) surface and find absolute deviations within 0.05 eV between DMET using the CASSCF/NEVPT2 solver, denoted as CAS-DMET/NEVPT2-DMET, and the nonembedded CASSCF/NEVPT2 approach. Next, we establish the practicality of DMET by extending it to larger supercells for the OV defect and a neutral silicon vacancy in diamond where the use of nonembedded CASSCF/NEVPT2 is extremely expensive.

Localized Active Space Pair-Density Functional Theory

Accurate quantum chemical methods for the prediction of spin-state energy gaps for strongly correlated systems are computationally expensive and scale poorly with the size of the system. This makes calculations for many experimentally interesting molecules impractical even with abundant computational resources. Previous work has shown that the localized active space (LAS) self-consistent field (SCF) method can be an efficient way to obtain multiconfiguration SCF wave functions of comparable quality to the corresponding complete active space (CAS) ones. To obtain quantitative results, a post-SCF method is needed to estimate the complete correlation energy. One such method is multiconfiguration pair-density functional theory (PDFT), which calculates the energy based on the density and on-top pair density obtained from a multiconfiguration wave function. In this work, we introduce localized-active-space PDFT, which uses a LAS wave function for subsequent PDFT calculations. The method is tested by computing spin-state energies and gaps in conjugated organic molecules and a bimetallic compound and comparing to the corresponding CAS-PDFT values.

Multiconfiguration Density-Coherence Functional Theory

This paper presents a new theory called multiconfiguration density coherence functional theory (MC-DCFT). This theory provides a new route to define density functionals for multiconfiguration wave functions, in particular by using the one-particle density matrix in the coordinate representation. The theory is illustrated by calculating the dissociation curve of four heteronuclear and homonuclear diatomic molecules, namely, H₂, F₂, N₂, and HF, using density coherence functionals converted from PBE, BLYP, and PW91. By introducing two parameters in the converted density functionals, we are able to calculate bond dissociation energies of comparable accuracy as those calculated by multiconfiguration pair-density functional theory (MC-PDFT) and complete active space second-order perturbation theory (CASPT2). This demonstrates that it would be possible to build a successful multiconfiguration density functional theory based on density coherence.

Analytic gradients for multiconfiguration pair-density functional theory with density fitting: Development and application to geometry optimization in the ground and excited states

Density fitting reduces the computational cost of both energy and gradient calculations by avoiding the computation and manipulation of four-index electron repulsion integrals. With this algorithm, one can efficiently optimize the geometries of large systems with an accurate multireference treatment. Here, we present the derivation of multiconfiguration pair-density functional theory for energies and analytic gradients with density fitting. Six systems are studied, and the results are compared to those obtained with no approximation to the electron repulsion integrals and to the results obtained by complete active space second-order perturbation theory. With the new approach, there is an increase in the speed of computation with a negligible loss in accuracy. Smaller grid sizes have also been used to reduce the computational cost of multiconfiguration pair-density functional theory with little effect on the optimized geometries and gradient values.

A New Mixing of Nonlocal Exchange and Nonlocal Correlation with Multiconfiguration Pair-Density Functional Theory

We propose a hybrid multiconfiguration pair-density functional theory (HMC-PDFT) that is a weighted average of complete-active-space self-consistent-field (CASSCF) and multiconfiguration pair-density functional theory (MC-PDFT) energies with a semiempirical parameter to control the fraction of CASSCF energy. We also explore a more general two-parameter hybrid method with a scaled correlation energy that allows us to compare to the recently proposed λ-MC-PDFT method. We scan the parameter space for the scaled-correlation method using test sets consisting of electronic excitation energies and diatomic bond energies, and we find no significant improvement by introducing the scaling parameter. We find that unscaled HMC-PDFT offers significantly improved accuracy over both CASSCF and the original MC-PDFT for a wide range of systems, and we present as an example of this approach "tPBE0", the "translated" MC-PDFT generalization of the popular PBE0 hybrid Kohn-Sham density functional.

Analytic gradients for state-averaged multiconfiguration pair-density functional theory

Analytic gradients are important for efficient calculations of stationary points on potential energy surfaces, for interpreting spectroscopic observations, and for efficient direct dynamics simulations. For excited electronic states, as are involved in UV-Vis spectroscopy and photochemistry, analytic gradients are readily available and often affordable for calculations using a state-averaged complete active space self-consistent-field (SA-CASSCF) wave function. However, in most cases, a post-SA-CASSCF step is necessary for quantitative accuracy, and such calculations are often too expensive if carried out by perturbation theory or configuration interaction. In this work, we present the analytic gradients for multiconfiguration pair-density functional theory based on SA-CASSCF wave functions, which is a more affordable alternative. A test set of molecules has been studied with this method, and the stationary geometries and energetics are compared to values in the literature as obtained by other methods. Excited-state geometries computed with state-averaged pair-density functional theory have similar accuracy to those from complete active space perturbation theory at the second-order.

Spin-State Ordering in Metal-Based Compounds Using the Localized Active Space Self-Consistent Field Method

Quantitatively accurate calculations for spin-state ordering in transition-metal complexes typically demand a robust multiconfigurational treatment. The poor scaling of such methods with increasing size makes them impractical for large, strongly correlated systems. Density matrix embedding theory (DMET) is a fragmentation approach that can be used to specifically address this challenge. The single-determinantal bath framework of DMET is applicable in many situations, but it has been shown to perform poorly for molecules characterized by strong correlation when a multiconfigurational self-consistent field solver is used. To ameliorate this problem, the localized active space self-consistent field (LASSCF) method was recently described. In this work, LASSCF is applied to predict spin-state energetics in mono- and di-iron systems, and we show that the model offers an accuracy equivalent to that of CASSCF but at a substantially lower computational cost. Performance as a function of basis set and active space is also examined.

Multiconfigurational Self-Consistent Field Theory with Density Matrix Embedding: The Localized Active Space Self-Consistent Field Method

Density matrix embedding theory (DMET) is a fully quantum-mechanical embedding method which shows great promise as a method of defeating the inherent exponential cost scaling of multiconfigurational wave function-based calculations by breaking large systems into smaller, coupled subsystems. However, we recently [ Pham et al. J. Chem. Theory Comput. 2018 , 14 , 1960 .] encountered evidence that the approximate single-determinantal bath picture inherent to DMET is sometimes problematic when the complete active space self-consistent field (CASSCF) is used as a solver and the method is applied to realistic models of strongly correlated molecules. Here, we show this problem can be defeated by generalizing DMET to use a multiconfigurational wave function as a bath without sacrificing practically attractive features of DMET, such as a second-quantization form of the embedded subsystem Hamiltonian, by dividing the active space into unentangled active subspaces each localized to one fragment. We introduce the term localized active space (LAS) to refer to this kind of wave function. The LAS bath wave function can be obtained by the DMET algorithm itself in a self-consistent manner, and we refer to this approach, introduced here for the first time, as the localized active space self-consistent field (LASSCF) method. LASSCF exploits a modified DMET algorithm, but it is a variational wave function method; it does not require DMET's ambiguous error function minimization, and it reproduces full-molecule CASSCF in cases where comparable DMET calculations fail. Our results for test calculations on the nitrogen double-bond dissociation potential energy curves of several diazene molecules suggest that LASSCF can be an appropriate starting point for a perturbative treatment. Outside of the context of embedding, the LAS wave function is inherently an attractive alternative to a CAS wave function because of its favorable cost scaling, which is exponential only with respect to the size of individual fragment active subspaces, rather than the whole active space of the entire system.

Stochastic many-body perturbation theory for anharmonic molecular vibrations

A new quantum Monte Carlo (QMC) method for anharmonic vibrational zero-point energies and transition frequencies is developed, which combines the diagrammatic vibrational many-body perturbation theory based on the Dyson equation with Monte Carlo integration. The infinite sums of the diagrammatic and thus size-consistent first- and second-order anharmonic corrections to the energy and self-energy are expressed as sums of a few m- or 2m-dimensional integrals of wave functions and a potential energy surface (PES) (m is the vibrational degrees of freedom). Each of these integrals is computed as the integrand (including the value of the PES) divided by the value of a judiciously chosen weight function evaluated on demand at geometries distributed randomly but according to the weight function via the Metropolis algorithm. In this way, the method completely avoids cumbersome evaluation and storage of high-order force constants necessary in the original formulation of the vibrational perturbation theory; it furthermore allows even higher-order force constants essentially up to an infinite order to be taken into account in a scalable, memory-efficient algorithm. The diagrammatic contributions to the frequency-dependent self-energies that are stochastically evaluated at discrete frequencies can be reliably interpolated, allowing the self-consistent solutions to the Dyson equation to be obtained. This method, therefore, can compute directly and stochastically the transition frequencies of fundamentals and overtones as well as their relative intensities as pole strengths, without fixed-node errors that plague some QMC. It is shown that, for an identical PES, the new method reproduces the correct deterministic values of the energies and frequencies within a few cm(-1) and pole strengths within a few thousandths. With the values of a PES evaluated on the fly at random geometries, the new method captures a noticeably greater proportion of anharmonic effects.

Second-order many-body perturbation expansions of vibrational Dyson self-energies

Second-order many-body perturbation theories for anharmonic vibrational frequencies and zero-point energies of molecules are formulated, implemented, and tested. They solve the vibrational Dyson equation self-consistently by taking into account the frequency dependence of the Dyson self-energy in the diagonal approximation, which is expanded in a diagrammatic perturbation series up to second order. Three reference wave functions, all of which are diagrammatically size consistent, are considered: the harmonic approximation and diagrammatic vibrational self-consistent field (XVSCF) methods with and without the first-order Dyson geometry correction, i.e., XVSCF[n] and XVSCF(n), where n refers to the truncation rank of the Taylor-series potential energy surface. The corresponding second-order perturbation theories, XVH2(n), XVMP2[n], and XVMP2(n), are shown to be rigorously diagrammatically size consistent for both total energies and transition frequencies, yield accurate results (typically within a few cm(-1) at n = 4 for water and formaldehyde) for both quantities even in the presence of Fermi resonance, and have access to fundamentals, overtones, and combinations as well as their relative intensities as residues of the vibrational Green's functions. They are implemented into simple algorithms that require only force constants and frequencies of the reference methods (with no basis sets, quadrature, or matrix diagonalization at any stage of the calculation). The rules for enumerating and algebraically interpreting energy and self-energy diagrams are elucidated in detail.

Computational studies of the isomerization and hydration reactions of acetaldehyde oxide and methyl vinyl carbonyl oxide

Alkene ozonolysis is a major source of hydroxyl radical (*OH), the most important oxidant in the troposphere. Previous experimental and computational work suggests that for many alkenes the measured *OH yields should be attributed to the combined impact of both chemically activated and thermalized syn-alkyl Criegee intermediates (CIs), even though the thermalized CI should be susceptible to trapping by molecules such as water. We have used RRKM/master equation and variational transition state theory calculations to quantify the competition between unimolecular isomerization and bimolecular hydration reactions for the syn and anti acetaldehyde oxide formed in trans-2-butene ozonolysis and for the CIs formed in isoprene ozonolysis possessing syn-methyl groups. Statistical rate theory calculations were based on quantum chemical data provided by the B3LYP, QCISD, and multicoefficient G3 methods, and thermal rate constants were corrected for tunneling effects using the Eckart method. At tropospheric temperatures and pressures, all thermalized CIs with syn-methyl groups are predicted to undergo 1,4-hydrogen shifts from 2 to 8 orders of magnitude faster than they react with water monomer at its saturation number density. For thermalized anti acetaldehyde oxide, the rates of dioxirane formation and hydration should be comparable.