A CASSCF and CASPT2 study on the excited states of s-trans-formaldazine

The low-lying excited states of s-trans-formaldazine (H2CN-NCH2) have been investigated using the complete active space self-consistent field (CASSCF) and the multiconfigurational second-order perturbation (CASPT2) methods. The vertical excitation energies have been calculated at the state-average CASSCF and multistate CASPT2 levels employing the cc-pVTZ basis set. The photodissociation mechanisms starting from the S1 state have been determined. The lowest energy points along the seams of surface intersections have been located in both the Franck-Condon region and the N-N dissociation pathway in the S1 state. Once the system populates the S1 state, in the viewpoint of energy, the radiationless decay via S1/S0(3) conical intersection followed by the N-N bond fission in the ground-state is more favorable in comparison with the N-N dissociation process in the S1 state. A three-surface crossing region (S1/T1/T2), where the S1, T1, and T2 states intersect, was also found. However, the intersystem crossing process via S1/T1/T2 is not energetically competitive with the internal conversion via S1/S0(3).     

Effects of water re-location and cavity trimming on the CASPT2//CASSCF/AMBER excitation energy of Rhodopsin

We have used quantum-mechanics/molecular-mechanics computations based on ab initio multiconfigurational perturbation theory to determine and rationalize the effect of the re-location of one crystallographic water molecule on the vertical excitation energy of the visual pigment rhodopsin. It is found that the re-location of one water molecule to the opposite side of the 11-cis retinal chromophore leads to a large 0.7–0.8 Å contraction in the chromophore—counterion salt-bridge distance. In spite of this structural effect, the change in excitation energy is found to be limited (< 1.5 kcal mol^?1). Through an analysis of different rhodopsin models in terms of “components” (isolated chromophore, isolated chromophore—counterion ion-pair and models deprived of the counterion charges) we show that the limited change of the excitation energy can be related to a displacement of the retinal chromophore to a different spot of the protein cavity.     

CASSCF/CASPT2 and TD-DFT study of valence and Rydberg electronic transitions in fluorene, carbazole, dibenzofuran, and dibenzothiophene

A combination of multireference CASSCF/CASPT2 and time-dependent DFT (TD-B3P86) theoretical treatments was employed to test their predictions against recently proposed assignments of the vacuum-UV spectra of fluorene and its three heteroanalogues-dibenzofuran, carbazole, and dibenzothiophene-up to the ionization threshold. For the low-lying transitions, the theoretically based assignments are generally not problematic because of the well-resolved bands, although, even in this region, the two methods yield some opposing predictions. Further on toward the vacuum region, the assignments prove increasingly challenging because of predicted crowding of transitions, many of which exhibit significant intensity. Some of the transitions in this region and beyond-toward the ionization thresholds-are thus necessarily assigned only tentatively. Overall, the two methods are frequently found to complement each other well, and equivalent transitions usually appear as bracketed from the high- (CASPT2) and low- (TD-B3P86) energy sides.     

Low-lying electronic states of HNCS and its ions: a CASSCF/CASPT2 study

Complete active space self-consistent-field (CASSCF) and multiconfigurational second-order perturbation theory (CASPT2) calculations in conjunction with the ANO-L basis set were performed to investigate systematically the low-lying electronic states of HNCS and its ions in C _s symmetry. Our highly accurate calculation indicated that theoretically determined geometric parameters and harmonic vibrational frequencies for the ground-state X ¹A′ are in good agreement with observed experimental data. The geometry of triplet HNCS is clearly favored C ₁ symmetry, and the relative energy is predicted to be 3.000 eV (69.2 kcal/mol). The vertical transition energies for the selected excited states of HNCS were calculated at CASSCF/CASPT2/ANO-L level of theory based on CASSCF optimized geometry. Except for a few linear states of X ²Π (1²A′, 1²A″), 1⁴Σ⁻ (1⁴A″), and 1²Σ⁺ (3²A′) states of HNCS⁺, our results confirmed that the majority of excited states are twisted trans-bend structures. The existence of bound excited anion states has been found for the first time in HNCS⁻. A more elaborate examination of ionization potential of HNCS (AIP, VIP) than previous reports has been presented.     

Probing the rhodopsin cavity with reduced retinal models at the CASPT2//CASSCF/AMBER level of theory

We show that the ab initio CASPT2//CASSCF strategy previously used to investigate the ground and excited states of the chromophore of the vision receptor rhodopsin (Rh) in vacuo can be successfully implemented in a QM/MM scheme allowing for CASPT2//CASSCF/AMBER geometry optimization and excited state property evaluation in proteins. Two receptor models (Rh-1 and Rh-2) incorporating different reduced chromophores are investigated. It is shown that Rh-2 features a chromophore equilibrium structure with the correct helicity and a lambdamax that is only 52 nm blue-shifted from the observed value. This result should open the way to a qualitatively correct ab initio QM/MM modeling of the early excited state transient species involved in the vision process.     

CASSCF and CASPT2 study on O- and Cl-loss predissociation mechanisms of OClO (A(2)A2)

For studying O- and Cl-loss predissociation mechanisms of OClO (A(2)A(2)), we calculated O- and Cl-loss dissociation potential energy curves (adiabatic minimum-energy dissociation paths) of several low-lying doublet and quartet states at the CASPT2 level and located the MECPs (minimum energy crossing points) for many pairs of the potential energy surfaces (PESs) at the CASPT2 and CASSCF levels. On the basis of our calculation results (including the spin-orbit couplings at the MECPs), we predict three processes for O-loss predissociation of A(2)A(2) and four processes for Cl-loss predissociation of A(2)A(2). The most favorable process for O-loss predissociation is OClO (A(2)A(2)) → A(2)A(2)/1 (2)B(2) MECP → 1 (2)B(2) (1 (2)A') → O ((3)P(g)) + ClO (X(2)Π) (the first O-loss limit), and the needed energy for this process from X(2)B(1) is 2.92 eV. The most favorable process for Cl-loss predissociation is OClO (A(2)A(2)) → A(2)A(2)/1 (2)B(2) MECP → TS1 (1 (2)B(2)) → 1 (2)B(2)/1 (2)A(1) MECP → Cl ((2)P(u)) + O(2) (X(3)Σ(g)(-)) (the first limit), and the needed energy is 3.08 eV. In the previously suggested mechanisms (processes), the A(2)A(2) state was considered to go to the important 1 (2)B(2) state via 1 (2)A(1) (A(2)A(2) → 1 (2)A(1) → 1 (2)B(2)). In the present study we have found that the A(2)A(2) state can directly go to 1 (2)B(2) (at the located A(2)A(2)/1 (2)B(2) MECP the CASPT2 energy (relative to X(2)B(1)) and CASSCF spin-orbit coupling are 2.92 eV and 61.3 cm(-1), respectively). We have compared our processes (A(2)A(2) → 1 (2)B(2) → ...) with the processes (A(2)A(2) → 1 (2)A(1) → 1 (2)B(2) → ...) suggested in the previous MRCI studies and rewritten by us using our calculation results. Energetically the MRCI process for O-loss predissociation (to the first limit) is only slightly (0.13 eV) more favorable than our process, and the MRCI processes for Cl-loss predissociation (to the first and second limits) need the same energies as our processes. By considering the probabilities of radiationless transitions, the MRCI processes are less favorable than our processes since the MRCI processes proceed via more PES/PES crossings (more MECPs). The experimental facts concerning the photodissociation are explained.     

A CASSCF/CASPT2 study on the low-lying excited states of HSiCN, HSiNC and their ions

Equilibrium geometries of low-lying electronic states of cyanosilylene (HSiCN), isocyanosilylene (HSiNC), and their ions have been investigated using the complete active space self-consistent field (CASSCF) approach. The harmonic vibrational frequencies on the optimized geometries were calculated using the multiconfiguration linear response (MCLR) method. Taking the further correlation effects into account, the complete active space perturbation theory of second-order (CASPT2) was carried out for the energetic correction. The CASPT2 calculations have been performed to obtain the vertical excitation energies of selected low-lying excited states of HSiCN and HSiNC. Computed results show that the singlet-triplet splittings are calculated to be 0.99 and 1.30 eV for HSiCN and HSiNC, respectively. The vertical excitation energies of the lowest singlet and triplet excited states in HSiCN are lower than those in HSiNC. The first vertical ionization energy of HSiCN (10.04 eV) is higher than that of HSiNC (9.97 eV). The ground-state adiabatic electron affinities are found to be rather high, and the value of HSiCN (1.85 eV) higher than that of HSiNC (1.52 eV). The existences of dipole-bound excited negative ion states have been discovered within HSiCN and HSiNC.     

Properties of the emitting state of the green fluorescent protein resolved at the CASPT2//CASSCF/CHARMM level

A "brute-force" ab initio CASPT2//CASSCF/CHARMM computational approach is employed to investigate the properties of the emitting state of the wild-type green fluorescence protein. The results indicate that the emitting moiety corresponds to a slightly perturbed H2O- - -chromophore complex. Thus, the protein matrix seems to be designed in such a way to mimic an environment that is more similar to gas-phase than water solution.     

Probing highly efficient photoisomerization of a bridged azobenzene by a combination of CASPT2//CASSCF calculation with semiclassical dynamics simulation

Mechanism of phototriggered isomerization of azobenzene and its derivatives is of broad interest. In this paper, the S(0) and S(1) potential energy surfaces of the ethylene-bridged azobenzene (1) that was recently reported to have highly efficient photoisomerization were determined by ab initio electronic structure calculations at different levels and further investigated by a semiclassical dynamics simulation. Unlike azobenzene, the cis isomer of 1 was found to be more stable than the trans isomer, consistent with the experimental observation. The thermal isomerization between cis and trans isomers proceeds via an inversion mechanism with a high barrier. Interestingly, only one minimum-energy conical intersection was determined between the S(0) and S(1) states (CI) for both cis → trans and trans → cis photoisomerization processes and confirmed to act as the S(1) → S(0) decay funnel. The S(1) state lifetime is ～30 fs for the trans isomer, while that for the cis isomer is much longer, due to a redistribution of the initial excitation energies. The S(1) relaxation dynamics investigated here provides a good account for the higher efficiency observed experimentally for the trans → cis photoisomerization than the reverse process. Once the system decays to the S(0) state via CI, formation of the trans product occurs as the downhill motion on the S(0) surface, while formation of the cis isomer needs to overcome small barriers on the pathways of the azo-moiety isomerization and rotation of the phenyl ring. These features support the larger experimental quantum yield for the cis → trans photoisomerization than the trans → cis process.     

Origin, nature, and fate of the fluorescent state of the green fluorescent protein chromophore at the CASPT2//CASSCF resolution

Ab initio CASPT2//CASSCF relaxation path computations are employed to determine the intrinsic (e.g., in vacuo) mechanism underlying the rise and decay of the luminescence of the anionic form of the green fluorescent protein (GFP) fluorophore. Production and decay of the fluorescent state occur via a two-mode reaction coordinate. Relaxation along the first (totally symmetric) mode leads to production of the fluorescent state that corresponds to a planar species. The second (out-of-plane) mode controls the fluorescent state decay and mainly corresponds to a barrierless twisting of the fluorophore phenyl moiety. While a "space-saving" hula-twist conical intersection decay channel is found to lie only 5 kcal mol(-1) above the fluorescent state, the direct involvement of a hula-twist deformation in the decay is not supported by our data. The above results indicate that the ultrafast fluorescence decay observed for the GFP chromophore in solution is likely to have an intrinsic origin. The possible effects of the GFP protein cavity on the fluorescence lifetime of the investigated chromophore model are discussed.     

Energetics of cytosine singlet excited-state decay paths--a difficult case for CASSCF and CASPT2

Three deactivation paths for singlet excited cytosine are calculated at the CASPT2//CASSCF (complete active space second-order perturbation//complete active space self-consistent field) level of theory, using extended active spaces that allow for a reliable characterization of the paths and their energies. The lowest energy path, with a barrier of approximately 0.1 eV, corresponds to torsion of the C5-C6 bond, and the decay takes place at a conical intersection analogous to the one found for ethylene and its derivatives. There is a further path with a low energy barrier of approximately 0.2 eV associated with the (n(N),pi*) state which could also be populated with a low energy excitation. The path associated with a conical intersection between the ground and (n(O),pi*) states is significantly higher in energy (> 1 eV). The presence of minima on the potential energy surface for the (n,pi*) states that could contribute to the biexponential decay found in the gas phase was investigated, but could not be established unequivocally.     

A CASPT2//CASSCF study of the quartet excited state \tilde{a}^{4}A^{\prime\prime} of the HCNN radical

Complete active space self-consistent field and second-order multiconfigurational perturbation theory methods have been performed to investigate the quartet excited state ${\tilde{a}}^{4}{A^{\prime\prime}}$ potential energy surface of HCNN radical. Two located minima with respective cis and trans structures could easily dissociate to CH $({\tilde{a}}^{4}\Sigma^{-})$ and $N_{2} ({\tilde{X}}^{1}\Sigma_{\rm g}^{+})$ products with similar barrier of about 16.0 kcal/mol. In addition, four minimum energy crossing points on a surface of intersection between ${\tilde{a}}^{4}A^{\prime\prime}$ and X ( $X={\tilde{X}}^{2}A^{\prime\prime}$ and ${\tilde{A}}^{2}A^{\prime}$ ) states are located near to the minima. However, the intersystem crossing ${\tilde{a}}^{4}A^{\prime\prime} \rightarrow X$ is weak due to the vanishingly small spin–orbit interactions. It further indicates that the direct dissociation on the ${\tilde{a}}^{4}{A^{\prime\prime}}$ state is more favored. This information combined with the comparison with isoelectronic HCCO provides an indirect support to the recent experimental proposal of photodissociation mechanism of HCNN.     

Efficient photochemical merocyanine-to-spiropyran ring closure mechanism through an extended conical intersection seam. A model CASSCF/CASPT2 study

A mechanism of the thermal and photochemical bleaching of merocyanine to spiropyran is proposed on the basis of CASSCF/CASPT2 calculations on the 6-(2-propenyliden)cyclohexadienone model system. Our results suggest that this photochemical transformation takes place in two steps. First, the initially pumped 1(pi-pi) S2 undergoes radiationless decay to 1(n-pi) S1 via an extended S2/S1 conical intersection seam that runs approximately parallel to the trans-to-cis isomerization coordinate, a few kilocalories per mole higher in energy. Thus, S2 --> S1 internal conversion is possible at all values of the S2 trans-to-cis reaction coordinate. Second, on the S1 potential energy surface, there is a barrierless ring closure reaction path from the S1 cis minimum that leads to a peaked S1/S0 conical intersection where the deactivation to the ground state takes place. The inertia of the moving nuclei then drives the system toward the ground-state minimum of the 2H-chromene product. Thus, the extended seam topology of the S2/S1 conical intersection and the coordinate of the branching space of the S1/S0 conical intersection are essential to explain the efficiency and high speed of this reaction.     

Electronic spectra of uranyl chloride complexes in acetone: a CASSCF/CASPT2 investigation

A theoretical study is presented of the electronic spectra of the complexes UO(2)Cl(2)ac(4), UO(2)Cl(2)ac(3), [UO(2)Cl(3)ac(2)](-) and [UO(2)Cl(3)ac](-) (ac = acetone) using perturbation theory based on a complete-active-space type wavefunction (CASSCF/CASPT2). Both scalar relativistic effects and spin-orbit coupling were included in the calculations. The calculated excitation energies and oscillator strength values have been compared to the experimental absorption spectrum for uranyl chloride complexes in acetone solution, for chloride-to-uranyl ratios between two and three. The main purpose of this work was to investigate the origin of the remarkable intensity increase observed in the lower part of the experimental absorption spectra, upon addition of chloride to uranyl complexes in acetone. The calculated excitation energies for the different complexes are similar and closely correspond to the experimental data. However, in none of the theoretical spectra, the high intensities observed in the lower part of the experimental spectrum are reproduced.     

A CASSCF/CASPT2 study on the low‐lying electronic states of the CH3SS and its cation

Complete active space self‐consistent field (CASSCF) and multiconfiguration second‐order perturbation theory (CASPT2) calculations with contracted ANO‐RCC basis set were performed for low‐lying electronic states of CH₃SS and its cation in C_s symmetry. For the ground state X²A″ of CH₃SS, the calculated S‐S stretching mode is in good agreement with experimental reports. The electron transitions of CH₃SS⁺, X¹A′ → 1¹A″, X¹A′ → 2¹A′, and X¹A′ → 2¹A″, are predicted at 1.055, 3.247, and 3.841 eV. Moreover, the calculated adiabatic and vertical ionization potential and adiabatic affinity are in reasonable agreement with the experimental data. The CASSCF/CASPT2 potential energy curves (PECs) were calculated for S₂‐loss dissociation from the X²A″, 1²A′, and 2²A″ states. The electronic states of the CH₃ radical and S₂ molecule as the dissociation products were carefully determined by checking energies and geometries of the asymptote products. The S₂‐loss PEC for CH₃SS indicate that S₂‐loss dissociation occurs from the X²A″ state leading to CH₃ (1²A″) + S₂ (X³Σ), the 1²A′ and 2²A″ leading to CH₃ (1²A″) + S₂ (¹Δ_g). © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012.     

A CASSCF/CASPT2 insight into excited‐state intramolecular proton transfer of four imidazole derivatives

Excited‐state intramolecular proton transfer (ESIPT) of four imidazole derivatives, 2‐(2′‐hydroxyphenyl)imidazole (HPI), 2‐(2′‐hydroxyphenyl)benzimidazole (HPBI), 2‐(2′‐hydroxyphenyl)‐1H‐phenanthro[9,10‐d]imidazole (HPPI) and 2‐(2′‐hydroxyphenyl)‐1‐phenyl‐1H‐phenanthro[9,10‐d]imidazole (HPPPI), were studied by the sophisticated CASSCF/CASPT2 methodology. The state‐averaged SA‐CASSCF method was used to optimize their geometry structures of S₀ and S₁ electronic states, and the CASPT2 calculations were used for the calibration of all the single‐point energies, including the absorption and emission spectra. A reasonable agreement is found between the theoretical predictions and the available experimental spectral data. The forward ESIPT barriers of four target compounds gradually decrease with the increase of molecular size. On the basis of the present calculations, it is a plausible speculation that the larger the size, the faster is the ESIPT rate, and eventually, HPPPI molecule can undergo a completely barrierless ESIPT to the more stable S₁ keto form. Additionally, taking HPI as a representative example, the radiationless decays connecting the S₀ and S₁/S₀ conical intersection structures were also studied by constructing a linearly interpolated internal coordinate (LIIC) reaction path. The qualitative analysis shows that the LIIC barrier of HPI in the keto form is remarkably lower than that of its enol‐form, indicating that the former has a big advantage over the latter in the nonradiative process. © 2015 Wiley Periodicals, Inc.     

Dibenzo-p-dioxin. An ab initio CASSCF/CASPT2 study of the pi-pi* and n-pi* valence excited states

The pi-pi* and n-pi* valence excited states of dibenzo-p-dioxin (DD) were studied via the complete active space SCF and multiconfigurational second-order perturbation theory employing the cc-pVDZ basis set and the full pi-electron active spaces of 16 electrons in 14 active orbitals. The geometry and harmonic vibrational wavenumbers of the ground state correlate well with the experimental and other theoretical data. In particular, significant improvements over previously reported theoretical results are observed for the excitation energies. All of the pi-pi* excited states exhibit planar D(2h)minima. Thus no evidence was found for a C(2v) butterfly-like relaxation, although the wavenumbers of the b(3u) butterfly flapping mode proved exceedingly low in both the ground S(0)((1)A(g)) and the lowest dipole allowed excited S(1)((1)B(2u)) state. The calculations of oscillator strengths established the 2(1)B(2u) <-- 1(1)A(g) and 2(1)B(1u) <-- 1(1)A(g) transitions as by far the most intense, whereas the only allowed of the n-pi* transitions ((1)B(3u)) should possess only a modest intensity. Studies into dependence of the oscillator strengths on the extent of the butterfly-like folding showed that the electronic spectrum is more consistent with a folded equilibrium geometry assumed by DD in solution.     

Hydrogen transfer vs proton transfer in 7-hydroxy-quinoline.(NH3)3: a CASSCF/CASPT2 study

Multiconfigurational CASSCF and CASPT2 calculations were performed to investigate the enol --> keto tautomerization in the lowest singlet excited state of the 7-hydroxyquinoline.(NH3)3 cluster. Two different reaction mechanisms were explored. The first one corresponds to that proposed previously by Tanner et al. (Science 2003, 302, 1736) on the basis of experimental observations and CASSCF optimizations under Cs-symmetry constraints. This mechanism comprises four consecutive steps and involves nonadiabatic transitions between the valence 1pipi* state and a pisigma* Rydberg-type state, resulting in hydrogen-atom transfer. Single-point CASPT2 calculations corroborate that for Cs-symmetry pathways hydrogen-atom transfer is clearly preferred over proton transfer. The second mechanism, predicted by CASSCF optimizations without constraints, implies proton transfer along a pathway on the 1pipi* surface in which one or more ammonia molecules depart significantly from the molecular plane defined by the hydroxyquinoline ring. The results suggest that both mechanisms may be competitive with proton transfer being somewhat favorable over hydrogen-atom transfer.