Complete active space analysis of a reaction pathway: Investigation of the oxygen–oxygen bond formation

Implementation of frozen density embedding in CP2K and OpenMolcas: CASSCF wavefunctions embedded in a Gaussian and plane wave DFT environment

Most chemical processes happen at a local scale where only a subset of molecular orbitals is directly involved and only a subset of covalent bonds may be rearranged. To model such reactions, Density Functional Theory (DFT) is often inadequate, and the use of computationally more expensive correlated wavefunction (WF) methods is required for accurate results. Mixed-resolution approaches backed by quantum embedding theory have been used extensively to approach this imbalance. Based on the frozen density embedding freeze-and-thaw algorithm, we describe an approach to embed complete active space self-consistent field simulations run in the OpenMolcas code in a DFT environment calculated in CP2K without requiring any external tools. This makes it possible to study a local, active part of a chemical system in a larger and relatively static environment with a computational cost balanced between the accuracy of a WF method and the efficiency of DFT, which we test on environment-subsystem pairs. Finally, we apply the implementation to an oxygen molecule leaving an aluminum (111) surface and a ruthenium(IV) oxide (110) surface.

Strategies to Obtain Reliable Energy Landscapes from Embedded Multireference Correlated Wavefunction Methods for Surface Reactions

Embedded correlated wavefunction (ECW) theory is a powerful tool for studying ground- and excited-state reaction mechanisms and associated energetics in heterogeneous catalysis. Several factors are important to obtaining reliable ECW energies, critically the construction of consistent active spaces (ASs) along reaction pathways when using a multireference correlated wavefunction (CW) method that relies on a subset of orbital spaces in the configuration interaction expansion to account for static electron correlation, e.g., complete AS self-consistent field theory, in addition to the adequate partitioning of the system into a cluster and environment, as well as the choice of a suitable basis set and number of states included in excited-state simulations. Here, we conducted a series of systematic studies to develop best-practice guidelines for ground- and excited-state ECW theory simulations, utilizing the decomposition of NH3 on Pd(111) as an example. We determine that ECW theory results are relatively insensitive to cluster size, the aug-cc-pVDZ basis set provides an adequate compromise between computational complexity and accuracy, and that a fixed-clean-surface approximation holds well for the derivation of the embedding potential. Additionally, we demonstrate that a merging approach, which involves generating ASs from the molecular fragments at each configuration, is preferable to a creeping approach, which utilizes ASs from adjacent structures as an initial guess, for the generation of consistent potential energy curves involving open-d-shell metal surfaces, and, finally, we show that it is essential to include bands of excited states in their entirety when simulating excited-state reaction pathways.

How does multi-reference computation change the catalysis chemistry? DFT and CASPT2 studies of the Cu-catalysed coupling reactions between aryl iodides and β-diketones

The molecular mechanism of a Cu-catalysed coupling reaction was theoretically studied using density functional theory (DFT) and the complete active space self-consistent field method followed by the second-order perturbation theory (CASSCF/CASPT2) to investigate the effects of the strong electron correlation of the Cu centre on the reaction profile. Both DFT and CASSCF/CASPT2 calculations showed that the catalytic cycle proceeds via an oxidative addition (OA) reaction, followed by a reductive elimination (RE) reaction, where OA is the rate-determining step. Although the DFT-calculated activation energies of the OA and RE steps are highly dependent on the choice of functionals, the CASSCF/CASPT2 results are less affected by the choice of DFT-optimised geometries. Therefore, with a careful assessment based on the CASSCF/CASPT2 single-point energy evaluation, an optimal choice of the DFT geometry is of good qualitative use for energetics at the CASPT2 level of theory. Based on the changes in the electron populations of the 3d orbitals during the OA and RE steps, the characteristic features of the DFT-calculated electronic structure were qualitatively consistent with those calculated using the CASSCF method. Further electronic structure analysis by the natural orbital occupancy of the CASSCF wavefunction showed that the ground state is almost single-reference in this system and the strong electron correlation effect of the Cu centre can be dealt with using the MP2 or CCSD method, too. However, the slightly smaller occupation numbers of the 3dπ orbital in the course of reactions suggested that the electron correlation effect of the Cu(III) centre appears through the interaction between the 3dπ orbital and the C-I antibonding σ* orbital in the OA step, and between the 3dπ orbital and the Cu-C antibonding σ* orbital in the RE step.