Conical Intersection And Potential Energy Surface Features Of A Model Retinal Chromophore: Comparison Of Eom-cc And Multireference Methods

This work investigates the performance of equation-of-motion coupled-cluster (EOM-CC) methods for describing the changes in the potential energy surfaces of the penta-2,4-dieniminium cation, a reduced model of the retinal chromophore of visual pigments, due to dynamical electron correlation effects. The ground-state wave function of this model includes charge-transfer and diradical configurations whose weights vary along different displacements and are rapidly changing at the conical intersection between the ground and the first excited states, making the shape of the potential energy surface sensitive to a balanced description of nondynamical and dynamical correlation. Recently, variational (MRCISD) and perturbative (MRPT2) approaches for including dynamical correlation in CASSCF-based calculations were tested along three representative ground state paths. Here, we use the same three paths to compare the performance of single-reference EOM-CC methods against MRCISD and MRCISD+Q We find that the spin-flip variant of EOM-CCSD with perturbative inclusion of triple excitations (dT or IT) produces potential energy profiles of the two lowest electronic states in quantitative agreement with MRCISD+Q (our highest-quality reference method). The nonparallelity errors and differences in vertical energy differences of the two surfaces along these scans are less than 1.4 kcal/mol (EOM-SF-CCSD(dT) versus MRCISD+Q). For comparison, the largest error of MRCISD versus MRCISD+Q is 1.7 kcal/mol. Our results show that the EOM-CC methods provide an alternative to multireference approaches and may be used to study photochemical systems like the one used in this work

Combined Self-consistent-field And Spin-flip Tamm-dancoff Density Functional Approach To Potential Energy Surfaces For Photochemistry

We present a new approach to calculating potential energy surfaces for photochemical reactions by combining self-consistent-field calculations for single-reference ground and excited states with symmetry-corrected spin-flip Tamm-Dancoff approximation calculations for multireference electronic states. The method is illustrated by an application with the M05-2X exchange-correlation functional to cis-trans isomerization of the penta-2,4-dieniminium cation, which is a model (with three conjugated double bonds) of the protonated Schiff base of retinal. We find good agreement with multireference configuration interaction-plus-quadruples (MRCISD+Q) wave function calculations along three key paths in the strong-interaction region of the ground and first excited singlet states

Origin Of Fluorescence In 11-cis Locked Bovine Rhodopsin

The excited state lifetime of bovine rhodopsin (Rh) increases from ca. 100 fs to 85 ps when the C11=C12 bond of its chromophore is locked by a cyclopentene moiety (Rh5). To explain such an increase, we employed ab initio multiconfigurational quantum chemistry to construct computer models of Rh and Rh5 and to investigate the shape of their excited state potential energy surfaces in a comparative way. Our results show that the observed Rh5 fluorescence (lambda(f)(max) = 620 nm) is due to a previously unreported locally excited intermediate whose lifetime is controlled by a small energy barrier. The analysis of the properties and decay path of such an intermediate provides useful information for engineering rhodopsin variants with augmented fluorescence efficiencies

Assessment Of Density Functional Theory For Describing The Correlation Effects On The Ground And Excited State Potential Energy Surfaces Of A Retinal Chromophore Model

In the quest for a cost-effective level of theory able to describe a large portion of the ground and excited potential energy surfaces of large chromophores, promising approaches are rooted in various approximations to the exact density functional theory (DFT). In the present work, we investigate how generalized Kohn-Sham DFT (GKS-DFT), time-dependent DFT (TDDFT), and spin-restricted ensemble-DFT (REKS) methods perform along three important paths characterizing a model retinal chromophore (the penta-2,4-dieniminium cation) in a region of near-degeneracy (close to a conical intersection) with respect to reference high-level multiconfigurational wave function methods. If GKS-DFT correctly describes the closed-shell charge transfer state, only TDDFT and REKS approaches give access to the open-shell diradical, one which sometimes corresponds to the electronic ground state. It is demonstrated that the main drawback of the usual DFT-based methods lies in the absence of interactions between the charge transfer and the diradicaloid configurations. Hence, we test a new computational scheme based on the State-averaged REKS (SA-REKS) approach, which explicitly includes these interactions into account. The State-Interaction SA-REKS (SI-SA-REKS) method significantly improves on the REKS and the SA-REKS results for the target system. The similarities and differences between DFT and wave function-based approaches are analyzed according to (1) the active space dimensions of the wave function-based methods and (2) the relative electronegativities of the allyl and protonated Schiff base moieties

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

Probing Vibrationally Mediated Ultrafast Excited-state Reaction Dynamics With Multireference (caspt2) Trajectories

Excited-state trajectories computed at the complete active space second-order perturbation theory (CASPT2) reveal how vibrational excitation controls the molecular approach to the intersection space that drives the photodissociation of a prototypical halogenated methyl radical, namely CF2I. Translating the Franck-Condon structure along the ground-state CASPT2 vibrational modes in this system followed by propagating the displaced structures in the first excited doublet state simulates specific vibrational excitations and vibrationally mediated dynamics, respectively. Three distinct situations are encountered: the trajectories (i) converge to an energetically flat segment of the intersection space, (ii) locate a segment of the intersection space, and (iii) access a region where the intersection space degeneracy is lifted to form a ridge of avoided crossings. The computational protocol documented herein can be used as a tool to design control strategies based on selective excitation of vibrational modes, including adaptive feedback schemes using coherent light sources

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

Vacuum ultraviolet photoionization cross section of the hydroxyl radical

The absolute photoionization spectrum of the hydroxyl (OH) radical from 12.513 to 14.213 eV was measured by multiplexed photoionization mass spectrometry with time-resolved radical kinetics. Tunable vacuum ultraviolet (VUV) synchrotron radiation was generated at the Advanced Light Source. OH radicals were generated from the reaction of O(^1D) + H_2O in a flow reactor in He at 8 Torr. The initial O(^1D) concentration, where the atom was formed by pulsed laser photolysis of ozone, was determined from the measured depletion of a known concentration of ozone. Concentrations of OH and O(^3P) were obtained by fitting observed time traces with a kinetics model constructed with literature rate coefficients. The absolute cross section of OH was determined to be σ(13.436 eV) = 3.2 ± 1.0 Mb and σ(14.193 eV) = 4.7 ± 1.6 Mb relative to the known cross section for O(^3P) at 14.193 eV. The absolute photoionization spectrum was obtained by recording a spectrum at a resolution of 8 meV (50 meV steps) and scaling to the single-energy cross sections. We computed the absolute VUV photoionization spectrum of OH and O(^3P) using equation-of-motion coupled-cluster Dyson orbitals and a Coulomb photoelectron wave function and found good agreement with the observed absolute photoionization spectra

Toward An Understanding Of The Retinal Chromophore In Rhodopsin Mimics

Recently, a rhodopsin protein mimic was constructed by combining mutants of the cellular retinoic acid binding protein II (CRABPII) with an all-trans retinal chromophore. Here, we present a combine computational quantum mechanics/molecular mechanics (QM/MM) and experimental ultrafast kinetic study of CRABPII. We employ the QM/MM models to study the absorption (lambda(a)(max)), fluorescence (lambda(f)(max)), and reactivity of a CRABPII triple mutant incorporating the all-trans protonated chromophore (PSB-KLE-CRABPII). We also study the spectroscopy of the same mutant incorporating the unprotonated chromophore and of another double mutant incorporating the neutral unbound retinal molecule held inside the pocket. Finally, for PSB-KLE-CRABPII, stationary fluorescence spectroscopy and ultrafast transient absorption spectroscopy resolved two different evolving excited state populations which were computationally assigned to distinct locally excited and charge-transfer species. This last species is shown to evolve along reaction paths describing a facile isomerization of the biologically relevant 11-cis and 13-cis double bonds. This work represents a first exploratory attempt to model and study these artificial protein systems. It also indicates directions for improving the QM/MM models so that they could be more effectively used to assist the bottom-up design of genetically encodable probes and actuators employing the retinal chromophore

Vacuum ultraviolet photoionization cross section of the hydroxyl radical

The absolute photoionization spectrum of the hydroxyl (OH) radical from 12.513 to 14.213 eV was measured by multiplexed photoionization mass spectrometry with time-resolved radical kinetics. Tunable vacuum ultraviolet (VUV) synchrotron radiation was generated at the Advanced Light Source. OH radicals were generated from the reaction of O(^1D) + H_2O in a flow reactor in He at 8 Torr. The initial O(^1D) concentration, where the atom was formed by pulsed laser photolysis of ozone, was determined from the measured depletion of a known concentration of ozone. Concentrations of OH and O(^3P) were obtained by fitting observed time traces with a kinetics model constructed with literature rate coefficients. The absolute cross section of OH was determined to be σ(13.436 eV) = 3.2 ± 1.0 Mb and σ(14.193 eV) = 4.7 ± 1.6 Mb relative to the known cross section for O(^3P) at 14.193 eV. The absolute photoionization spectrum was obtained by recording a spectrum at a resolution of 8 meV (50 meV steps) and scaling to the single-energy cross sections. We computed the absolute VUV photoionization spectrum of OH and O(^3P) using equation-of-motion coupled-cluster Dyson orbitals and a Coulomb photoelectron wave function and found good agreement with the observed absolute photoionization spectra