The structure and energetics of Cr(CO)6 and Cr(CO)5

The geometric structure of Cr(CO)6 is optimized at the modified coupled pair functional (MCPF), single and double excitation coupled-cluster (CCSD) and CCSD(T) levels of theory (including a perturbational estimate for connected triple excitations), and the force constants for the totally symmetric representation are determined. The geometry of Cr(CO)5 is partially optimized at the MCPF, CCSD, and CCSD(T) levels of theory. Comparison with experimental data shows that the CCSD(T) method gives the best results for the structures and force constants, and that remaining errors are probably due to deficiencies in the one-particle basis sets used for CO. The total binding energies of Cr(CO)6 and Cr(CO)5 are also determined at the MCPF, CCSD, and CCSD(T) levels of theory. The CCSD(T) method gives a much larger total binding energy than either the MCPF or CCSD methods. An analysis of the basis set superposition error (BSSE) at the MCPF level of treatment points out limitations in the one-particle basis used. Calculations using larger basis sets reduce the BSSE, but the total binding energy of Cr(CO)6 is still significantly smaller than the experimental value, although the first CO bond dissociation energy of Cr(CO)6 is well described. An investigation of 3s3p correlation reveals only a small effect. In the largest basis set, the total CO binding energy of Cr(CO)6 is estimated to be 140 kcal/mol at the CCSD(T) level of theory, or about 86 percent of the experimental value. The remaining discrepancy between the experimental and theoretical value is probably due to limitations in the one-particle basis, rather than limitations in the correlation treatment. In particular an additional d function and an f function on each C and O are needed to obtain quantitative results. This is underscored by the fact that even using a very large primitive set (1042 primitive functions contracted to 300 basis functions), the superposition error for the total binding energy of Cr(CO)6 is 22 kcal/mol at the MCPF level of treatment

Dynamical Insights into the Decomposition of 1,2-Dioxetane

Chemiluminescence in 1,2-dioxetane occurs through a thermally activated decomposition reaction into two formaldehyde molecules. Both ground-state and nonadiabatic dynamics (including singlet excited states) of the decomposition reaction have been simulated, starting from the first O-O bond-breaking transition structure. The ground-state dissociation occurs between t = 30 fs and t = 140 fs. The so-called entropic trap leads to frustrated dissociations, postponing the decomposition reaction. Specific geometrical conditions are necessary for the trajectories to escape from the entropic trap and for dissociation to be possible. The singlet excited states participate as well in the trapping of the molecule: dissociation including the nonadiabatic transitions to singlet excited states now occurs from t = 30 fs to t = 250 fs and later. Specific regions of the seam of the S0/S1 conical intersections that would "retain" the molecule for longer on the excited state have been identified

Symmetry breaking in O4(+): An application of the Brueckner coupled-cluster method

A recent calculation of the antisymmetric stretch frequency for the rectangular structure of quartet O4(+) using the singles and doubles quadratic configuration interaction method with a perturbational estimate of connected triple excitations (QCISD(T)) method gave a value of 3710 cm(exp -1). This anomalous frequency is shown to be a consequence of symmetry breaking effects, which occur even though the QCISD(T) solution derived from a delocalized SCF reference function lies energetically well below the two localized (symmetry-broken) solutions at the equilibrium geometry. The symmetry breaking is almost eliminated at the CCSD level of theory, but the small remaining symmetry breaking effects are magnified at the CCSD(T) level of theory so that the antisymmetric stretch frequency is still significantly in error. The use of the Brueckner coupled cluster method, however, leads to a symmetrical solution which is free of symmetry breaking effects, with an antisymmetric stretch frequency of 1322 cm(exp -1), in good agreement with our earlier calculations using the complete active space self consistent field/complete active space state interaction (CASSCF/CASSI) method

How Do Methyl Groups Enhance the Triplet Chemiexcitation Yield of Dioxetane?

Chemiluminescence is the emission of light as a result of a nonadiabatic chemical reaction. The present work is concerned with understanding the yield of chemiluminescence, in particular how it dramatically increases upon methylation of 1,2-dioxetane. Both ground-state and nonadiabatic dynamics (including singlet excited states) of the decomposition reaction of various methyl-substituted dioxetanes have been simulated. Methyl-substitution leads to a significant increase in the dissociation time scale. The rotation around the O-C-C-O dihedral angle is slowed; thus, the molecular system stays longer in the "entropic trap" region. A simple kinetic model is proposed to explain how this leads to a higher chemiluminescence yield. These results have important implications for the design of efficient chemiluminescent systems in medical, environmental, and industrial applications

Analytic Gradients for Complete Active Space Pair-Density Functional Theory

Analytic gradient routines are a desirable feature for quantum mechanical methods, allowing for efficient determination of equilibrium and transition state structures and several other molecular properties. In this work, we present analytical gradients for multiconfiguration pair-density functional theory (MC-PDFT) when used with a state-specific complete active space self-consistent field reference wave function. Our approach constructs a Lagrangian that is variational in all wave function parameters. We find that MC-PDFT locates equilibrium geometries for several small- to medium-sized organic molecules that are similar to those located by complete active space second-order perturbation theory but that are obtained with decreased computational cost

The Ultrafast Photoisomerizations Of Rhodopsin And Bathorhodopsin Are Modulated By Bond Length Alternation And Hoop Driven Electronic Effects

Rhodopsin (Rh) and bathorhodopsin (bathoRh) quantum-mechanics/molecular-mechanics models based on ab initio multiconfigurational wave functions are employed to look at the light induced pi-bond breaking and reconstitution occurring during the Rh -\u3e bathoRh and bathoRh -\u3e Rh isomerizations. More specifically, semiclassical trajectory computations are used to compare the excited (S-1) and ground (S-0) state dynamics characterizing the opposite steps of the Rh/bathoRh photochromic cycle during the first 200 fs following photoexcitation. We show that the information contained in these data provide an unprecedented insight into the sub-picosecond pi-bond reconstitution process which is at the basis of the reactivity of the protein embedded 11-cis and all-trans retinal chromophores. More specifically, the data point to the phase and amplitude of the skeletal bond length alternation stretching mode as the key factor switching the chromophore to a bonding state. It is also confirmed/found that the phase and amplitude of the hydrogen-out-of-plane mode controls the stereochemical outcome of the forward and reverse photoisomerizations

Towards an accurate molecular orbital theory for excited states : Ethene, butadiene, and hexatriene

A newly proposed quantum chemical approach for ab initio calculations of electronic spectra of molecular systems is applied to the molecules ethene, trans‐1,3‐butadiene, and trans‐trans‐1,3,5‐hexatriene. The method has the aim of being accurate to better than 0.5 eV for excitation energies and is expected to provide structural and physical data for the excited states with good reliability. The approach is based on the complete active space (CAS) SCF method, which gives a proper description of the major features in the electronic structure of the excited state, independent of its complexity, accounts for all near degeneracy effects, and includes full orbital relaxation. Remaining dynamic electron correlation effects are in a subsequent step added using second order perturbation theory with the CASSCF wave function as the reference state. The approach is here tested in a calculation of the valence and Rydberg excited singlet and triplet states of the title molecules, using extended atomic natural orbital (ANO) basis sets. The ethene calculations comprised the two valence states plus all singlet and triplet Rydberg states of 3s, 3p, and 3d character, with errors in computed excitation energies smaller than 0.13 eV in all cases except the V state, for which the vertical excitation energy was about 0.4 eV too large. The two lowest triplet states and nine singlet states were studied in butadiene. The largest error (0.37 eV) was found for the 2 1Bu state. The two lowest triplet and seven lowest singlet states in hexatriene had excitation energies in error with less than 0.17 [email protected] ; [email protected] ; [email protected]

Dynamic Electron Correlation Effects On The Ground State Potential Energy Surface Of A Retinal Chromophore Model

The ground state potential energy surface of the retinal chromophore of visual pigments (e.g., bovine rhodopsin) features a low-lying conical intersection surrounded by regions with variable charge transfer and diradical electronic structures. This implies that dynamic electron correlation may have a large effect on the shape of the force fields driving its reactivity. To investigate this effect, we focus on mapping the potential energy for three paths located along the ground state CASSCF potential energy surface of the penta-2,4-dieniminium cation taken as a minimal model of the retinal chromophore. The first path spans the bond length alternation coordinate and intercepts a conical intersection point. The other two are minimum energy paths along two distinct but kinetically competitive thermal isomerization coordinates. We show that the effect of introducing the missing dynamic electron correlation variationally (with MRCISD) and perturbatively (with the CASPT2, NEVPT2, and XMCQDPT2 methods) leads, invariably, to a stabilization of the regions with charge transfer character and to a significant reshaping of the reference CASSCF potential energy surface and suggesting a change in the dominating isomerization mechanism. The possible impact of such a correction on the photoisomerization of the retinal chromophore is discussed

Mapping The Excited State Potential Energy Surface Of A Retinal Chromophore Model With Multireference And Equation-of-motion Coupled-cluster Methods

The photoisomerization of the retinal chromophore of visual pigments proceeds along a complex reaction coordinate on a multidimensional surface that comprises a hydrogen-out-of-plane (HOOP) coordinate, a bond length alternation (BLA) coordinate, a single bond torsion and, finally, the reactive double bond torsion. These degrees of freedom are coupled with changes in the electronic structure of the chromophore and, therefore, the computational investigation of the photochemistry of such systems requires the use of a methodology capable of describing electronic structure changes along all those coordinates. Here, we employ the penta-2,4-dieniminium (PSB3) cation as a minimal model of the retinal chromophore of visual pigments and compare its excited state isomerization paths at the CASSCF and CASPT2 levels of theory. These paths connect the cis isomer and the trans isomer of PSB3 with two structurally and energetically distinct conical intersections (CIs) that belong to the same intersection space. MRCISD+Q energy profiles along these paths provide benchmark values against which other ab initio methods are validated. Accordingly, we compare the energy profiles of MRPT2 methods (CASPT2, QD-NEVPT2, and XMCQDPT2) and EOM-SF-CC methods (EOM-SF-CCSD and EOM-SF-CCSD(dT)) to the MRCISD+Q reference profiles. We find that the paths produced with CASSCF and CASPT2 are topologically and energetically different, partially due to the existence of a locally excited region on the CASPT2 excited state near the Franck-Condon point that is absent in CASSCF and that involves a single bond, rather than double bond, torsion. We also find that MRPT2 methods as well as EOM-SF-CCSD(dT) are capable of quantitatively describing the processes involved in the photoisomerization of systems like PSB3