CASPT2 Theory Research Articles

Multireference Perturbation Theory Combined with PCM and RISM Solvation Models: A Benchmark Study for Chemical Energetics.

The polarizable continuum model (PCM) has been one of the most widely used approaches to take into account the solvation effect in quantum chemical calculations. In this paper, we performed a series of benchmark calculations to assess the accuracy of the PCM scheme combined with the second-order complete-active-space perturbation theory (CASPT2) for molecular systems in polar solvents. For solute molecules with extensive conjugated π orbitals, exemplified by elongated conjugated arylcarbenes, we have incorporated the ab initio density matrix renormalization group algorithm into the PCM-CASPT2 method. In the previous work, we presented a combination of the DMRG-CASPT2 method with the reference interaction site model (RISM) theory for describing the solvation effect using the radial distribution function and compared its performance to the widely used density-functional approaches (PCM-TD-DFT). The work here allows us to further show a more thorough assessment of the RISM model compared to the PCM with an equal level of the wave function treatment, the (DMRG-)CASPT2 theory, toward a high-accuracy electronic structure calculations for solvated chemical systems. With the exception that the PCM models are not capable of properly describing the hydrogen bondings, accuracy of the PCM-CASPT2 model is in most cases quite comparable to the RISM counterpart.

Singlet Oxygen Oxidation of the Radical Cations of 8-Oxo-2'-deoxyguanosine and Its 9-Methyl Analogue: Dynamics, Potential Energy Surface, and Products Mediated by C5-O2 -Addition.

8-Oxo-2'-deoxyguanosine (OG) is the most common DNA lesion. Notably, OG becomes more susceptible to oxidative damage than the undamaged nucleoside, forming mutagenic products in vivo. Herein the reactions of singlet O2 with the radical cations of 8-oxo-2'-deoxyguanosine (OG.+ ) and 9-methyl-8-oxoguanine (9MOG.+ ) were investigated using ion-molecule scattering mass spectrometry, from which barrierless, exothermic O2 -addition products were detected for both reaction systems. Corroborated by static reaction potential energy surface constructed using multi-reference CASPT2 theory and molecular dynamics simulated in the presence of the reactants' kinetic and internal energies, the C5-terminal O2 -addition was pinpointed as the most probable reaction pathway. By elucidating the reaction mechanism, kinetics and dynamics, and reaction products and energetics, this work constitutes the first report unraveling the synergetic damage of OG by ionizing radiation and singlet O2 .

An alternative choice of the zeroth-order Hamiltonian in CASPT2 theory.

A zeroth-order Hamiltonian based on Koopmans matrices for complete active space second-order perturbation theory (CASPT2) is presented. This Hamiltonian involves three types of Fock matrices. The original CASPT2 Fock matrix is retained for excitation classes where the excitation does not change the number of electrons in the complete active space (CAS). For excitation classes involving a change in the number of electrons in the CAS, two alternative Fock matrices corresponding to either positive or negative ionization of the CAS are introduced. These are constructed such that they exactly reproduce the Koopmans matrices for a singly ionized CAS. Test calculations indicate that the method gives better excitation energies than CASPT2 without using empirical parameters, for example, the ionization potential-electron affinity shift, which is also designed to improve excitation energies. The method is also less prone to intruder states than conventional CASPT2. Moreover, the dissociation curve for the chromium dimer looks much more reasonable than the one obtained with conventional CASPT2.

Coordination of N2 ligands to lanthanum: the complexes La (N2)1\u20138

The structures and bonding properties of La(N2)x (x = 1–8) complexes were investigated by density functional theory (DFT) computations using the B3LYP exchange-correlation functional in conjunction with a quasi-relativistic pseudopotential for La. The quality of the DFT electronic structures was confirmed in selected cases by relativistic multireference calculations using CASPT2 theory. From the end-on and side-on dinitrogen coordination modes in general, the structures with end-on coordination were found to be the more stable. The first coordination sphere of the complexes is filled by eight and six N2 ligands in the end-on and side-on type species, respectively. The main bonding interaction is the donation of La 5d valence electrons to anti-bonding orbitals of N2 resulting in characteristic elongation of the NN bonds. These directional interactions determine the (from steric point of view in several cases less logic) equilibrium molecular structures. The charge transfer resulted in partial charges up to 1.5 e of the originally neutral components (La, N2) leading to electrostatic attractive interactions which compose the minor contribution in the bonding.

Kinetics of the Methanol Reaction with OH at Interstellar, Atmospheric, and Combustion Temperatures.

The OH radical is the most important radical in combustion and in the atmosphere, and methanol is a fuel and antifreeze additive, model biofuel, and trace atmospheric constituent. These reagents are also present in interstellar space. Here we calculate the rate constants, branching ratios, and kinetic isotope effects (KIEs) of the hydrogen abstraction reaction of methanol by OH radical in a broad temperature range of 30-2000 K, covering interstellar space, the atmosphere, and combustion by using the competitive canonical unified statistical (CCUS) model in both the low-pressure and high-pressure limits and, for comparison, the pre-equilibrium model. Coupled cluster CCSD(T)-F12a theory and multi-reference CASPT2 theory were used to carry out benchmark calculations of the stationary points on the potential energy surface to select the most appropriate density functional method for direct dynamics calculations of rate constants. We find a significant effect of the anharmonicity of high-frequency modes of transition states on the low-temperature rate constant, and we show how tunneling leads to an unusual negative temperature dependence of the rate constants in the range 200 K > T > 100 K. The calculations also demonstrate the importance of the extent of stabilization of the pre-reactive complex. The capture rate for the formation of the complex is the dominant dynamical bottleneck for T < 100 K, and it leads to weak temperature dependence of the rate below 100 K in the high-pressure-limit of the CCUS model. We also report the pressure dependence of branching ratios (which are hard to measure so theory is essential) and the KIEs, and we report an unusual nonmonotonic variation of the KIE in the high-pressure limit at low temperatures.

Analytical Derivative Coupling for Multistate CASPT2 Theory

The probability of nonradiative transitions in photochemical dynamics is determined by the derivative couplings, the couplings between different electronic states through the nuclear degrees of freedom. Efficient and accurate evaluation of the derivative couplings is, therefore, of central importance to realize reliable computer simulations of photochemical reactions. In this work, the derivative couplings for multistate multireference second-order perturbation theory (MS-CASPT2) and its "extended" variant (XMS-CASPT2) are studied, in which we present an algorithm for their analytical evaluation. The computational costs for evaluating the derivative couplings are essentially the same as those for calculating the nuclear energy gradients. The geometries and energies calculated with XMS-CASPT2 for small molecules at minimum energy conical intersections (MECIs) are in good agreement with those computed by multireference configuration interaction. As numerical examples, MECIs are optimized using XMS-CASPT2 for stilbene and a green fluorescent protein model chromophore (the 4-para-hydroxybenzylidene-1,2-dimethyl-imidazolin-5-one anion).

Optical Excitations in MnO and MnO:ZnO via Embedded CASPT2 Theory and Their Implications for Solar Energy Conversion

MnO:ZnO was recently proposed as a novel, potentially visible-light absorbing catalyst with many advantages for water splitting. However, little is known about its optical absorption properties and...

A theoretical study of the low-lying excited states of thieno[3,4-b]pyrazine

The low-lying electronic excited states of thieno[3,4-b]pyrazine have been studied using the multiconfigurational second-order perturbation CASPT2 theory with extended atomic natural orbital basis sets. The CASPT2 results allow for a full interpretation of the electronic absorption and emission spectra and provide valuable information for the rationalization of the experimental data. The nature, position, and intensity of the spectral bands have been analyzed in detail. A preliminary comparative study of the ground-state geometry of thieno[3,4-b]pyrazine has been performed at the coupled cluster single and doubles and density functional theory levels using a variety of correlation-consistent basis sets. Thieno[3,4-b]pyrazine exhibits a polyene-like structure in the ground state due to the bond localization in the pyrazine moiety. An aromatization of the pyrazine unit is predicted for the lowest-energy electronic excited states.

Electronic Excitation in a Saturated Chain: An MS-CASPT2 Treatment of the Anti Conformer of n-Tetrasilane

The singlet−singlet electronic spectrum of the anti conformer of n-tetrasilane has been studied using multiconfigurational wave functions (CASSCF), second-order perturbation theory (CASPT2), and its multi-state extension (MS-CASPT2), in conjunction with large ANO-type basis sets including Rydberg functions. The calculations include the 4s, 4p, and 3d members of the Rydberg series converging on the first ionization. Mixing of valence and Rydberg states observed in the CASSCF wave functions is not fully rectified by single-reference CASPT2 theory, whereas the MS-CASPT2 method separates the valence and Rydberg states effectively. At the MS-CASPT2 level, six valence excited states have been found below the lowest Rydberg transition, predicted at 7.40 eV. In terms of natural orbitals, they correspond to single electron promotions from the σ-symmetry HOMO to four σ* and two π* valence orbitals. Vertical dipole-allowed valence transition energies (oscillator strengths) are computed at 6.33 eV (1.12), 6.68 eV (0.01), and 6.96 eV (0.15), for excitation to 11Bu, 11Au, and 21Bu states, respectively. Three vertical dipole-forbidden valence transitions are predicted to be interleaved among them at 6.55 eV (21Ag), 6.87 eV (31Ag), and 7.10 eV (11Bg). The results are consistent with the available experimental data and are easily rationalized in terms of a simple Hückel model of σ delocalization.

A theoretical study of the low-lying excited states of ozone

A theoretical study has been performed on the five lowest excited states of the ozone molecule using multiconfigurational second-order perturbation theory (CASPT2). The predicted order of states is: 3A 2 ( T 0 = 1.15 eV), 3B 2 ( T 0 = 1.33 eV, 3B 1 ( T 0 = 1.33 eV), 1A 2 ( T 0 = 1.44 eV) and 1B 1 ( T 0 = 1.88 eV). Corresponding experimental data are: 1.18, 1.30, 1.45, 1.58, and 2.05 eV, respectively. Equilibrium geometries, harmonic frequencies of symmetric vibrations, and vertical excitation energies are also reported. The dissociation limit D e for the ground state of ozone is found to be 1.08 eV, in agreement with the experimental value (1.13 eV). The calculations make use of a modified Fock operator in the CASPT2 theory. Relative energies of states with a different number of open shells are substantially improved. The modified CASPT2 method was checked by calculating spectroscopic constants of the oxygen molecule.