Initial excited‐state relaxation of locked retinal protonated schiff base chromophore. An insight from coupled cluster and multireference perturbation theory calculations

The initial S₁ excited‐state relaxation of retinal protonated Schiff base (RPSB) analog with central C11C12 double bond locked by eight‐membered ring (locked‐11.8) was investigated by means of multireference perturbation theory methods (XMCQDPT2, XMS‐CASPT2, MS‐CASPT2) as well as single‐reference coupled‐cluster CC2 method. The analysis of XMCQDPT2‐based geometries reveals rather weak coupling between in‐plane and out‐of‐plane structural evolution and minor energetical relaxation of three locked‐11.8 conformers. Therefore, a strong coupling between bonds length inversion and backbone out‐of‐plane deformation resulting in a very steep S₁ energy profile predicted by CASSCF/CASPT2 calculations is in clear contradiction with the reference XMCQDPT2 results. Even though CC2 method predicts good quality ground‐state structures, the excited‐state structures display more advanced torsional deformation leading to ca. 0.2 eV exaggerated energy relaxation and significantly red shifted (0.4–0.7 eV) emission maxima. According to our findings, the initial photoisomerization process in locked‐11.8, and possibly in other RPSB analogs, studied fully (both geometries and energies) by multireference perturbation theory may be somewhat slower than predicted by CASSCF/CASPT2 or CC2 methods. © 2018 Wiley Periodicals, Inc.     

Predictions of the geometries and fluorescence emission energies of oxyluciferins

The complete active space self-consistent field (CASSCF) method and multiconfigurational second-order perturbation theory (CASPT2) have been used to study the structures and spectra of oxyluciferins (OxyLH2). The ground and lowest-lying singlet excited states geometries have been optimized using CASSCF. CASPT2 has been used to predict relaxed emission energies. The focus is on the lowest-lying singlet excited states of the anionic keto and enol forms of OxyLH2(-1) at the optimized excited-state geometries. The planar keto and enol forms of OxyLH2(-1) are minima on both the S0 and the S1 potential energy surfaces. The twisted keto and enol forms of OxyLH2(-1) are transition states on the S0 and S1 potential energy surfaces. The S1 --> S0 fluorescence emission energies are in the range of 54.2-58.4 kcal/mol for the anionic planar keto forms of OxyLH2, and in the range of 55.7-63.2 kcal/mol for the anionic enol forms of OxyLH2. S0 and S1 potential energy surfaces and thus are not implicated in the emission spectra in the gas phase.     

乙基硫自由基及阴、阳离子低电子激发态从头算研究

采用CASSCF方法和6-311++(3df, 3pd)基组以及Cs对称性优化了乙基硫自由基和阳、阴离子3种分子的12个电子态的几何构型. 利用二级微扰方法(CASPT2)对这12个电子态做了单点能校正. 通过比较自由基与阴阳离子的能量, 得出了绝热电子亲和势和绝热电子电离能, 与实验结果在允许误差范围内基本一致.     

Conical Intersections and Low‐Lying Electronic States of Tetrafluoroethylene

The low‐lying electronic states of tetrafluoroethylene (C₂F₄) are characterized theoretically for the first time using equation‐of‐motion coupled cluster theory (EOM‐CCSD), and complete active space self‐consistent field (CASSCF) and second‐order perturbation theory (CASPT2). Computations are performed for vertical excitation energies, equilibrium geometries, minimum‐energy conical intersections, and potential energy curves along three geometric coordinates: 1) twisting of the F?C?C?F dihedral angle, 2) pyramidalization of the CF₂ group, and 3) migration of a fluorine atom resulting in an ethylidene‐like (CF₃CF) structure. The results suggest two relaxation pathways from the Rydberg‐3s excited electronic state to the ground state. These relaxation pathways are discussed in conjunction with the femtosecond photoionization spectroscopy results of Trushin et al. [ChemPhysChem­ 2004 , 5, 1389].     

Theoretical studies on structures and electronic spectra of linear carbon chains C2nH+ (n = 1−5)

The density functional theory (DFT) and the complete active space self‐consistent‐field (CASSCF) method have been used for full geometry optimization of carbon chains C_2nH⁺ (n = 1–5) in their ground states and selected excited states, respectively. Calculations show that C_2nH⁺ (n = 1–5) have stable linear structures with the ground state of X³Π for C₂H⁺ or X³Σ^? for other species. The excited‐state properties of C_2nH⁺ have been investigated by the multiconfigurational second‐order perturbation theory (CASPT2), and predicted vertical excitation energies show good agreement with the available experimental values. On the basis of our calculations, the unsolved observed bands in previous experiments have been interpreted. CASSCF/CASPT2 calculations also have been used to explore the vertical emission energy of selected low‐lying states in C_2nH⁺ (n = 1–5). Present results indicate that the predicted vertical excitation and emission energies of C_2nH⁺ have similar size dependences, and they gradually decrease as the chain size increases. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Conical intersections in thymine

The mechanisms which are responsible for the radiationless deactivation of the npi* and pipi* excited singlet states of thymine have been investigated with multireference ab initio methods (the complete-active-space self-consistent-field (CASSCF) method and second-order perturbation theory with respect to the CASSCF reference (CASPT2)) as well as with the CC2 (approximated singles and doubles coupled-cluster) method. The vertical excitation energies, the equilibrium geometries of the 1npi*and 1pipi* states, as well as their adiabatic excitation energies have been determined. Three conical intersections of the S1 and S0 energy surfaces have been located. The energy profiles of the excited states and the ground state have been calculated with the CASSCF method along straight-line reaction paths leading from the ground-state equilibrium geometry to the conical intersections. All three conical intersections are characterized by strongly out-of-plane distorted geometries. The lowest-energy conical intersection (CI1) arises from a crossing of the lowest 1pipi* state with the electronic ground state. It is found to be accessible in a barrierless manner from the minimum of the 1pipi* state, providing a direct and fast pathway for the quenching of the population of the lowest optically allowed excited states of thymine. This result explains the complete diffuseness of the absorption spectrum of thymine in supersonic jets. The lowest vibronic levels of the optically nearly dark 1npi* state are predicted to lie below CI1, explaining the experimental observation of a long-lived population of dark excited states in gas-phase thymine.     

Theoretical computation of low‐lying electronic states of HCNS: A CASPT2 study

Complete active space self‐consistent field (CASSCF) and complete active space second‐order perturbation theory (CASPT2) calculations in conjunction with the aug‐cc‐pVTZ basis set have been used to investigate the low‐lying electronic states of thiofulminic acid (HCNS), HCNS⁺, and HCNS^?. The result of geometry optimization using CASPT2/aug‐cc‐pVTZ shows that theoretically determined geometric parameters and harmonic vibrational frequencies for the HCNS ground state X¹Σ⁺(X¹A′) are in agreement with previous studies. The ionization energies, the electron affinity energies, the adiabatic excitation energies, and vertical excitation energies have been calculated and the corresponding cation and anion states are identified. By calculating adiabatic electron affinity, the states of HCNS^? have been identified to contain both π orbital states (X²A′ and 1²A″) and dipole‐bond states (1⁴A′ and 1⁴A″). © 2012 Wiley Periodicals, Inc.     

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.     

Molecular parameters of tetraatomic carbonyls X2CO and XYCO (X,Y = H,F, Cl) in the ground and lowest excited electronic states,part 1: A test of ab initio methods

Geometrical parameters of tetraatomic carbonyl molecules X₂CO and XYCO (X, Y = H, F, Cl) in the ground (S₀) and lowest excited singlet (S₁) and triplet (T₁) electronic states as well as values of barriers to inversion in S₁ and T₁ states and S₁ ← S₀ and T₁ ← S₀ adiabatic transition energies were systematically investigated by means of various quantum‐chemical techniques. The following methods were tested: HF, MP2, CIS, CISD, CCSD, EOM‐CCSD, CCSD(T), CR‐EOM‐CCSD(T), CASSCF, MR‐MP2, CASPT2, CASPT3, NEVPT2, MR‐CISD, and MR‐AQCC within cc‐pVTZ and cc‐pVQZ basis sets. The accuracy of quantum‐chemical methods was estimated in comparison with experimental data and rather accurate structures of excited electronic states were obtained. MP2 and CASPT2 methods appeared to be the most efficient and CCSD(T), CR‐EOM‐CCSD(T), and MR‐AQCC the most accurate. It was found that at equilibrium all the molecules under study are nonplanar in S₁ and T₁ electronic states with CO out‐of‐plane angle ranging from 34° (H₂CO, S₁) to 52° (F₂CO, T₁), and height of barrier to inversion varying from 300 (H₂CO, S₁) to 11,000 (F₂CO, T₁) cm^?1. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Assignment on the X,A, B,C, and D states of the C6H5Br+ cation based on high‐level calculations

Geometries, frequencies, and energies of the 1²B₁, 1²A₂, 1²B₂, 2²B₁, 2²B₂, and 1²A₁, of the C₆H₅Br⁺ ion were calculated by using CASSCF and CASPT2 methods in conjunction with an ANO‐RCC basis. The CASPT2//CASSCF adiabatic excitation energies and CASPT2 relative energies for the six states are in good agreement with experiment. The X, A, B, C, and D electronic states of the C₆H₅Br⁺ ion were assigned to be X²B₁, A²A₂, B²B₂, C²B₁, and D²B₂ based on the CASSCF and CASPT2 calculations. The assignment on the D state of the C₆H₅Br⁺ ion is different from the previously published works. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Spiro[4.4]nonatetraene and its Positive and Negative Radical Ions: Molectronic Structure Investigations

The electronic structure of spiro[4.4]nonatetraene 1 as well as that of its radical anion and cation were studied by different spectroscopies. The electron‐energy‐loss spectrum in the gas phase revealed the lowest triplet state at 2.98 eV and a group of three overlapping triplet states in the 4.5 – 5.0 eV range, as well as a number of valence and Rydberg singlet excited states. Electron‐impact excitation functions of pure vibrational and triplet states identified various states of the negative ion, in particular the ground state with an attachment energy of 0.8 eV, an excited state corresponding to a temporary electron attachment to the 2b₁ MO at an attachment energy of 2.7 eV, and a core excited state at 4.0 eV. Electronic‐absorption spectroscopy in cryogenic matrices revealed several states of the positive ion, in particular a richly structured first band at 1.27 eV, and the first electronic transition of the radical anion. Vibrations of the ground state of the cation were probed by IR spectroscopy in a cryogenic matrix. The results are discussed on the basis of density‐functional and CASSCF/CASPT2 quantum‐chemical calculations. In their various forms, the calculations successfully rationalized the triplet and the singlet (valence and Rydberg) excitation energies of the neutral molecule, the excitation energies of the radical cation, its IR spectrum, the vibrations excited in the first electronic absorption band, and the energies of the ground and the first excited states of the anion. The difference of the anion excitation energies in the gas and condensed phases was rationalized by a calculation of the Jahn‐Teller distortion of the anion ground state. Contrary to expectations based on a single‐configuration model for the electronic states of 1 , it is found that the gap between the first two excited states is different in the singlet and the triplet manifold. This finding can be traced to the different importance of configuration interaction in the two multiplicity manifolds.     

Is the choice of a standard zeroth‐order hamiltonian in CASPT2 ansatz optimal in calculations of excitation energies in protonated and unprotonated schiff bases of retinal?

To account for systematic error of CASPT2 method empirical modification of the zeroth‐order Hamiltonian with Ionization Potential‐Electron Affinity (IPEA) shift was introduced. The optimized IPEA value (0.25 a.u.), called standard IPEA (S‐IPEA), was recommended but due to its unsatisfactory performance in multiple metallic and organic compounds it has been questioned lately as a general parameter working properly for all molecules under CASPT2 study. As we are interested in Schiff bases of retinal, an important question emerging from this conflict of choice, to use or not to use S‐IPEA, is whether the introduction of the modified zeroth‐order Hamiltonian into CASPT2 ansatz does really improve their energetics. To achieve this goal, we assessed an impact of the IPEA shift value, in a range of 0–0.35 a.u., on vertical excitation energies to low‐lying singlet states of two protonated (RPSBs) and two unprotonated (RSBs) Schiff bases of retinal for which experimental data in gas phase are available. In addition, an effect of geometry, basis set, and active space on computed VEEs is also reported. We find, that for these systems, the choice of S‐IPEA significantly overestimates both S₀→S₁ and S₀→S₂ energies and the best theoretical estimate, in reference to the experimental data, is provided with either unmodified zeroth‐order Hamiltonian or small value of the IPEA shift in a range of 0.05–0.15 a.u., depending on active space and basis set size, equilibrium geometry, and character of the excited state. © 2018 Wiley Periodicals, Inc.     

Multiple Decay Mechanisms and 2D‐UV Spectroscopic Fingerprints of Singlet Excited Solvated Adenine‐Uracil Monophosphate

The decay channels of singlet excited adenine uracil monophosphate (ApU) in water are studied with CASPT2//CASSCF:MM potential energy calculations and simulation of the 2D‐UV spectroscopic fingerprints with the aim of elucidating the role of the different electronic states of the stacked conformer in the excited state dynamics. The adenine ¹L_a state can decay without a barrier to a conical intersection with the ground state. In contrast, the adenine ¹L_b and uracil S(U) states have minima that are separated from the intersections by sizeable barriers. Depending on the backbone conformation, the CT state can undergo inter‐base hydrogen transfer and decay to the ground state through a conical intersection, or it can yield a long‐lived minimum stabilized by a hydrogen bond between the two ribose rings. This suggests that the ¹L_b, S(U) and CT states of the stacked conformer may all contribute to the experimental lifetimes of 18 and 240 ps. We have also simulated the time evolution of the 2D‐UV spectra and provide the specific fingerprint of each species in a recommended probe window between 25 000 and 38 000 cm^?1 in which decongested, clearly distinguishable spectra can be obtained. This is expected to allow the mechanistic scenarios to be discerned in the near future with the help of the corresponding experiments. Our results reveal the complexity of the photophysics of the relatively small ApU system, and the potential of 2D‐UV spectroscopy to disentangle the photophysics of multichromophoric systems.     

CASSCF and CASPT2 studies on the structures, transition energies, and dipole moments of ground and excited states for azulene

The electronic structure of azulene molecule has been studied. We have obtained the optimized structures of ground and singlet excited states by using the complete active space self-consistent-field (CASSCF) method, and calculated vertical and 0-0 transition energies between the ground and excited states with second-order M?ller-Plesset perturbation theory (CASPT2). The CASPT2 calculations indicate that the bond-equalized C(2v) structure is more stable than the bond-alternating C(s) structure in the ground state. For a physical understanding of electronic structure change from C(2v) to C(s), we have performed the CASSCF calculations of Duschinsky matrix describing mixing of the b(2) vibrational mode between the ground (1A(1)) and the first excited (1B(2)) states based on the Kekule-crossing model. The CASPT2 0-0 transition energies are in fairly good agreement with experimental results within 0.1-0.3 eV. The CASSCF oscillator strengths between the ground and excited states are calculated and compared with experimental data. Furthermore, we have calculated the CASPT2 dipole moments of ground and excited states, which show good agreement with experimental values.     

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 CAS study on S‐loss and O‐loss dissociation mechanisms of the SO2+ ion in the C,D, and E states

In the present work, we mainly study dissociation of the C ²B₁, D²A₁, and E²B₂ states of the SO₂⁺ ion using the complete active‐space self‐consistent field (CASSCF) and multiconfiguration second‐order perturbation theory (CASPT2) methods. We first performed CASPT2 potential energy curve (PEC) calculations for S‐ and O‐loss dissociation from the X, A, B, C, D, and E primarily ionization states and many quartet states. For studying S‐loss predissociation of the C, D, and E states by the quartet states to the first, second, and third S‐loss dissociation limits, the CASSCF minimum energy crossing point (MECP) calculations for the doublet/quartet state pairs were performed, and then the CASPT2 energies and CASSCF spin‐orbit couplings were calculated at the MECPs. Our calculations predict eight S‐loss predissociation processes (via MECPs and transition states) for the C, D, and E states and the energetics for these processes are reported. This study indicates that the C and D states can adiabatically dissociate to the first O‐loss dissociation limit. Our calculations (PEC and MECP) predict a predissociation process for the E state to the first O‐loss limit. Our calculations also predict that the E²B₂ state could dissociate to the first S‐ and O‐loss limits via the A²B₂ ← E²B₂ transition. On the basis of the 13 predicted processes, we discussed the S‐ and O‐loss dissociation mechanisms of the C, D, and E states proposed in the previous experimental studies. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Theoretical study of Pu and Am tetracarbide molecules

The electronic structure and ground‐state molecular properties of Pu and Am tetracarbides have been investigated by relativistic multireference calculations using CASSCF/CASPT2 theory as well as by density functional theory in conjunction with relativistic pseudopotentials. The CASSCF/CASPT2 treatment has been extended by spin–orbit coupling effects for selected species using the CAS state‐interaction method. The five atoms can form various structural isomers, from which 12 ones have been identified in our study. The electronic ground state in both molecules corresponds to a planar fan‐type structure of C_2v symmetry, in which the actinide atom is connected to a bent C₄ moiety. The other structures are much higher in energy, the ones computed in this study appear between 250 and 1050 kJ/mol. The bonding characteristics in the most relevant structures have been analyzed on the basis of the valence molecular orbitals and natural bond orbital analysis. The most stable structures have been characterized by their spectroscopic (vibrational and electron) properties. © 2014 Wiley Periodicals, Inc.     

Photophysics of Schiff Bases: Theoretical Study of Salicylidene Methylamine

The proton‐transfer reaction in a model aromatic Schiff base, salicylidene methylamine (SMA), in the ground and in the lowest electronically‐excited singlet states, is theoretically analyzed with the aid of second‐order approximate coupled‐cluster model CC2, time‐dependent density functional theory (TD‐DFT) using the Becke, three‐parameter Lee–Yang–Parr (B3LYP) functional, and complete active space perturbation theory CASPT2 electronic structure methods. Computed vertical‐absorption spectra for the stable ground‐state isomers of SMA fully confirm the photochromism of SMA. The potential‐energy profiles of the ground and the lowest excited singlet state are calculated and four photophysically relevant isomeric forms of SMA; α, β, γ, and δ are discussed. The calculations indicate two S₁/S₀ conical intersections which provide non‐adiabatic gates for a radiationless decay to the ground state. The photophysical scheme which emerges from the theoretical study is related to recent experimental results obtained for SMA and its derivatives in the low‐temperature argon matrices (J. Grzegorzek, A. Filarowski, Z. Mielke, Phys. Chem. Chem. Phys. 2011 , 13, 16596–16605). Our results suggest that aromatic Schiff bases are potential candidates for optically driven molecular switches.     

Electric properties of the water molecule in 1A1, 1B1, and 3B1 electronic states: Refined CASSCF and CASPT2 calculations

Summary The dipole moments and dipole polarizabilities of the ¹A₁, ¹B₁, and ³B₁ electronic states of the water molecule have been calculated by using the CASSCF approach followed by the evaluation of the dynamic electron correlation contribution by the second-order perturbation scheme CASPT2. All calculations have been carried out in a specifically extended ANO basis set which accounts for the Rydberg character of the two excited states. In order to estimate the correctness and accuracy of the present data a scan over a variety of different active spaces for the CASSCF wave function has been made. The present results are superior to earlier CASSCF calculations, although their qualitative features remain essentially the same. The dipole moments in ¹B₁ and ³B₁ states are predicted to be about 0.49 a.u. and 0.33 a.u., respectively, and have the opposite orientation with respect to the ground state dipole moment. The dipole polarizability tensors of the excited states are characterized by high anisotropy and are dominated by the in-plane component perpendicular to the symmetry axis. All their components are found to be about an order of magnitude larger than those of the ground state polarizability tensor. The excitation energy dependence on the choice of the active orbital space in the CASSCF reference function is also considered and the analysis of the present data concludes in the concept of what is called the mutually compatible active spaces for the two states involved in excitation. All CASPT2 results are in good agreement with the results of recent calculations carried out in the framework of the open-shell coupled cluster formalism. This agreement confirms the high efficiency of the CASSCF/CASPT2 approach to the treatment of the electron correlation effects.