Three-state conical intersections in nucleic acid bases

The involvement of three-state conical intersections in the photophysics and radiationless decay processes of the nucleobases has been investigated using multireference configuration interaction methods. Three-state conical intersections have been located for the pyrimidine base, uracil, and the purine base, adenine. In uracil, a three-state degeneracy between the S(0), S(1), and S(2) states has been located at 6.2 eV above the ground-state minimum energy. This energy is 0.4 eV higher than vertical excitation to S(2) and at least 1.3 eV higher than the two-state conical intersections found previously. In adenine, two different three-state degeneracies between the S(1), S(2), and S(3) states have been located at energies close to the vertical excitation energies. The energetics of these three-state conical intersections suggest they can play a role in a radiationless decay pathway present in adenine. The existence of two different seams of three-state conical intersections indicates that these features are common and complicate the potential energy surfaces of adenine and possibly many other aromatic molecules.     

On the characterization of three-state conical intersections using a group homomorphism approach: the two-state degeneracy spaces

Second-order degenerate perturbation theory, in conjunction with the group homomorphism method for describing a similarity transformation, are used to characterize the subspace of two-state conical intersections contained in the branching space of a three-state conical intersection. It is shown by explicit calculation, using the lowest three-state conical intersection of (CH)3N2, that a second-order treatment yields highly accurate absolute energies, even at significant distances from the reference point of three-state intersection. The excellent agreement between the second order and ab initio results depends on the average energy component, which is computed using 5 first-order terms and 15 second-order terms. The second-order absolute energy change over the range rho = 0.0-0.3 au, where rho is the distance from the three-state conical intersection in the branching space coordinates, is approximately 6500 and 9500 cm(-1) for the E(1=2) and E(2=3) seams, respectively, with the maximum ab initio energy deviation from degeneracy of 200 cm(-1) occurring at rho = 0.3 au. The characteristic parameters gIJ and hIJ are also predicted to great accuracy, even at large rho, with the error growing to only 10-15% at rho = 0.3 au.     

Ab initio molecular dynamics of excited-state intramolecular proton transfer around a three-state conical intersection in malonaldehyde

Excited-state potential energy surface (PES) characterization is carried out at the CASSCF and MRSDCI levels, followed by ab initio dynamics simulation of excited-state intramolecular proton transfer (ESIPT) on the S2(pipi*) state in malonaldehyde. The proton-transfer transition state lies close to an S2/S1 conical intersection, leading to substantial coupling of proton transfer with electronic relaxation. Proton exchange proceeds freely on S2, but its duration is limited by competition with twisting out of the molecular plane. This rotamerization pathway leads to an intersection of the three lowest singlet states, providing the first detailed report of ab initio dynamics around a three-state intersection (3SI). There is a significant energy barrier to ESIPT on S1, and further pyramidalization of the twisted structure leads to the minimal energy S1/S0 intersection and energetic terminal point of excited-state dynamics. Kinetics and additional mechanistic details of these pathways are discussed. Significant depletion of the spectroscopic state and recovery of the ground state is seen within the first 250 fs after photoexcitation.     

Beyond two-state conical intersections. Three-state conical intersections in low symmetry molecules: the allyl radical

Using multireference configuration interaction expansions comprised of over 7 million configuration state functions, three-state conical intersections are reported for the closely spaced, spectroscopically observed (tilde)B(2A1), (tilde)C(2B1), and (tilde)D(2B2) states (in C(2v) symmetry) of the allyl radical. These conical intersections of states which were previously assigned as the 3,4,5(2)A states and are here reassigned as the 4,5,6(2)A states, are expected to be accessible using optical probes. This conclusion is obtained from the structure of the minimum energy point on the 4,5,6(2)A three-state conical intersection seam which is similar to the equilibrium structure of the ground (tilde)X(2A2) state and only 1.1 eV above the (tilde)D(2B2) state at its equilibrium geometry. The seam of three-state degeneracies joins two two-state seams of conical intersection, the 4,5(2)A and 5,6(2)A conical intersection seams. The energy of the minimum energy point on the 4,5(2)A two-state seam is only 0.15 eV above that of the (tilde)D(2B2) state at its equilibrium structure. Three-state intersections are also reported for the 3,4,5(2)A states.     

Three-state conical intersections in cytosine and pyrimidinone bases

Three-state conical intersections have been located and characterized for cytosine and its analog 5-methyl-2-pyrimidinone using multireference configuration-interaction ab initio methods. The potential energy surfaces for each base contain three different three-state intersections: two different S(0)-S(1)-S(2) intersections (gs/pi pi(*)/n(N)pi(*) and gs/pi pi(*)/n(O)pi(*)) and an S(1)-S(2)-S(3) intersection (pi pi(*)/n(N)pi(*)/n(O)pi(*)). Two-state seam paths from these intersections are shown to be connected to previously reported two-state conical intersections. Nonadiabatic coupling terms have been calculated, and the effects of the proximal third state on these quantities are detailed. In particular, it is shown that when one of these loops incorporates more than one seam point, there is a profound and predictable effect on the phase of the nonadiabatic coupling terms, and as such provides a diagnostic for the presence and location of additional seams. In addition, it is shown that each of the three three-state conical intersections located on cytosine and 5-methyl-2-pyrimidinone is qualitatively similar between the two bases in terms of energies and character, implying that, like with the stationary points and two-state conical intersections previously reported for these two bases, there is an underlying pattern of energy surfaces for 2-pyrimidinone bases, in general, and this pattern also includes three-state conical intersections.     

The curvature of the conical intersection seam: an approximate second-order analysis

We present a method for analyzing the curvature (second derivatives) of the conical intersection hyperline at an optimized critical point. Our method uses the projected Hessians of the degenerate states after elimination of the two branching space coordinates, and is equivalent to a frequency calculation on a single Born-Oppenheimer potential-energy surface. Based on the projected Hessians, we develop an equation for the energy as a function of a set of curvilinear coordinates where the degeneracy is preserved to second order (i.e., the conical intersection hyperline). The curvature of the potential-energy surface in these coordinates is the curvature of the conical intersection hyperline itself, and thus determines whether one has a minimum or saddle point on the hyperline. The equation used to classify optimized conical intersection points depends in a simple way on the first- and second-order degeneracy splittings calculated at these points. As an example, for fulvene, we show that the two optimized conical intersection points of C2v symmetry are saddle points on the intersection hyperline. Accordingly, there are further intersection points of lower energy, and one of C2 symmetry--presented here for the first time--is found to be the global minimum in the intersection space.     

Nuclear dynamics for a three-state Jahn-Teller model system

We report wavepacket dynamics on a model system with a three-state conical intersection. Quantum wavepacket dynamics using the multiconfigurational time-dependent Hartree method have been carried out for the T ? (e + t(2)) Jahn-Teller problem, using a Jahn-Teller vibronic model Hamiltonian. The effects of the magnitude of the coupling parameters and of the initial position of the wavepacket on the dynamics around the three-state conical intersection have been considered. It was found that the effect of the coupling strength is not dramatic for the population transfer in most cases, but the details of the dynamics and the involvement of the different modes are affected by it.     

环丁酮光化学反应中S1/T1/S0交叉点的重要作用

势能面间的交叉在光化学反应中起着重要的作用 ,是由激发态反应物到基态产物发生无辐射跃迁的机制 .在本文中 ,我们用 CASSCF和态平均 CASSCF方法分别对环丁酮光化学反应的势能剖面及 S1,T1和 S0三个势能面间交叉进行了研究 .结果发现 ,基态和三态产物的形成是通过 S1,T1和 S0三个势能面交叉于同一区域 (称为 S1/T1/S0交叉点 )这一有效途径完成的 .　　60年代末 ,实验 [1－ 9]发现环丁酮和其它烷基酮 ,如丙酮、环戊酮的光化学反应机理很不一致 .主要体现在 ,i)环丁酮 (n,π态 )的α解离发生在 S1态势能面上 ,而其它烷基酮 (n,π态 …     

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.     

The fluorescence mechanism of 5-methyl-2-pyrimidinone: an ab initio study of a fluorescent pyrimidine analog

The photophysically important potential energy surfaces of the fluorescent pyrimidine analog 5-methyl-2-pyrimidinone have been explored using multireference configuration-interaction ab initio methods at three levels of dynamical correlation, all of which support a fluorescence mechanism. At vertical excitation S1 (dark, n(N)pi*) and S2 (bright, pipi*) are almost degenerate at 4.4 eV, with S3 (dark, n(O)pi*) at 5.1 eV. The excited system can follow the S1-S2 seam of conical intersections, accessible from the Franck-Condon region, to its minimum and then evolve from this conical intersection on the S1 (pipi*) surface to a global minimum. At lower levels of correlation, the S1 surface shows two minima separated by a barrier of up to 0.18 eV. The secondary minimum found at the lower levels of correlation becomes the global minimum with higher correlation. The S1 population at this minimum can be trapped from accessing the lowest energy S0-S1 (pipi*/gs) conical intersection by an energy gap at least 0.3-0.4 eV higher than the S1 minimum. The calculated emission energy from this minimum is 2.80 eV. Gradient pathways connecting important S1 geometries are presented, as well as other excited state conical intersections.     

Ultrafast nonradiative decay of electronically excited States of malachite green: ab initio calculations

We have investigated the nonradiative deactivation process of malachite green in the singlet excited states, S(1) and S(2), by high-level ab initio quantum chemical calculations using the CASPT2//CASCF approach. The deactivation pathways connecting the Franck-Condon region and conical intersection regions are identified. The initial population in the S(1) state is on a flat surface and the relaxation involves a rotation of phenyl rings, which leads the molecule to reach the conical intersection between the S(1) and S(0) states, where it efficiently decays back to the ground state. There exists a small barrier connecting the Franck-Condon and conical intersection regions on the S(1) potential energy surface. The decay mechanism from the S(2) state also involves the twisting motion of phenyl rings. In contrast to the excitation to the S(1) state, the initial population is on a downhill ramp potential and the barrierless relaxation through the rotation of substituted phenyl rings is expected. During the course of relaxation, the molecule switches to the S(1) state at the conical intersection between S(2) and S(1), and then it decays back to the ground state through the intersection between S(1) and S(0). In relaxation from both S(1) and S(2), large distortion of phenyl rings is required for the ultrafast nonradiative decay to the ground state.     

A time-dependent wave packet study of the vibronic and spin-orbit interactions in the dynamics of Cl((2)P)+H(2)-->HCl(X (1)Sigma(g) (+))+H((2)S) reaction

We investigate the vibronic and spin-orbit (SO) coupling effects in the state-selected dynamics of the title reaction with the aid of a time-dependent wave packet approach. The ab initio potential energy surfaces of Capecchi and Werner [Science 296, 715 (2002)] have been employed for this purpose. Collinear approach of the Cl((2)P) atom to the H(2) molecule splits the degeneracy of the (2)P state and gives rise to (2)Sigma and (2)Pi electronic states. These two surfaces form a conical intersection at this geometry. These states transform as 1 (2)A('), 1 (2)A("), and 2 (2)A('), respectively, at the nonlinear configurations of the nuclei. In addition, the SO interaction due to Cl atom further splits these states into (2)Sigma(1/2), (2)Pi(3/2), and (2)Pi(1/2) components at the linear geometry. The ground-state reagent Cl((2)P(3/2))+H(2) correlates with (2)Sigma(1/2) and (2)Pi(3/2), where as the SO excited reagent Cl(*)((2)P(1/2))+H(2) correlates with (2)Pi(1/2) at the linear geometry. In order to elucidate the impact of the vibronic and SO coupling effects on the initial state-selected reactivity of these electronic states we carry out quantum scattering calculations based on a flux operator formalism and a time-dependent wave packet approach. In this work, total reaction probabilities and the time dependence of electronic population of the system by initiating the reaction on each of the above electronic states are presented. The role of conical intersection alone on the reaction dynamics is investigated with a coupled two-state model and for the total angular momentum J=0 (neglecting the electronic orbital angular momentum) both in a diabatic as well as in the adiabatic electronic representation. The SO interaction is then included and the dynamics is studied with a coupled three-state model comprising six diabatic surfaces for the total angular momentum J=0.5 neglecting the Coriolis Coupling terms of the Hamiltonian. Companion calculations are carried out for the uncoupled adiabatic and diabatic surfaces in order to explicitly reveal the impact of two different surface coupling mechanisms in the dynamics of this prototypical reaction.     

Intramolecular charge transfer in 4-aminobenzonitriles does not necessarily need the twist

In electron donor/acceptor species such as 4-(dimethylamino)benzonitrile (DMABN), the excitation to the S(2) state is followed by internal conversion to the locally excited (LE) state. Dual fluorescence then becomes possible from both the LE and the twisted intramolecular charge-transfer (TICT) states. A detailed mechanism for the ICT of DMABN and 4-aminobenzonitrile (ABN) is presented in this work. The two emitting S(1) species are adiabatically linked along the amino torsion reaction coordinate. However, the S(2)/S(1) CT-LE radiationless decay occurs via an extended conical intersection "seam" that runs almost parallel to this torsional coordinate. At the lowest energy point on this conical intersection seam, the amino group is untwisted; however, the seam is accessible for a large range of torsional angles. Thus, the S(1) LE-TICT equilibration and dual fluorescence will be controlled by (a) the S(1) torsional reaction path and (b) the position along the amino group twist coordinate where the S(2)/S(1) CT-LE radiationless decay occurs. For DMABN, population of LE and TICT can occur because the two species have similar stabilities. However, in ABN, the equilibrium lies in favor of LE, as a TICT state was found at much higher energy with a low reaction barrier toward LE. This explains why dual fluorescence cannot be observed in ABN. The S(1)-->S(0) deactivation channel accessible from the LE state was also studied.     

Three-state trajectory surface hopping studies of the photodissociation dynamics of formaldehyde on ab initio potential energy surfaces

Full-dimensional, three-state, surface hopping calculations of the photodissociation dynamics of formaldehyde are reported on ab initio potential energy surfaces (PESs) for electronic states S(1), T(1), and S(0). This is the first such study initiated on S(1) with ab initio-calculated spin-orbit couplings among the three states. We employ previous PESs for S(0) and T(1), and a new PES for S(1), which we describe here, as well as new spin-orbit couplings. The time-dependent electronic state populations and the branching ratio of radical products produced from S(0) and T(1) states and that of total radical products and molecular products at three total energies are calculated. Details of the surface hopping dynamics are described, and a novel pathway for isomerization on T(1) via S(0) is reported. Final translational energy distributions of H + HCO products from S(0) and T(1) are also reported as well as the translational energy distribution and final rovibrational distributions of H(2) products from the molecular channel. The present results are compared to previous trajectory calculations initiated from the global minimum of S(0). The roaming pathway leading to low rotational distribution of CO and high vibrational population of H(2) is observed in the present calculations.     

A simulation of the photoelectron spectrum of pyrazolide

Building on previous theoretical and spectroscopic studies of the pyrazolyl radical, a new three-state quasidiabatic Hamiltonian is reported which reproduces not only the equilibrium geometries and harmonic frequencies of the nominal X (2)A(2) state and low-lying A (2)B(1) excited state, but also the minimum energy points on the lowest two-state (X (2)A(2), A (2)B(1)) and three-state (X (2)A(2), A (2)B(1), B (2)B(2)) seams of conical intersection. The three-state Hamiltonian includes all terms through second order in both the diagonal and off-diagonal blocks. Its construction is accomplished in two steps. First, a nascent Hamiltonian, centered at the lowest energy two-state conical intersection, is determined using ab initio gradients and derivative couplings. Then, the nascent Hamiltonian is improved by optimizing selected contributions to the second-order coefficients to better reproduce relevant minima and harmonic frequencies. This Hamiltonian is then expressed in a basis tailored to describe the neutral states of interest under the multimode vibronic coupling approximation. The vibronic Hamiltonian is diagonalized to obtain negative ion photoelectron spectra for pyrazolide-h(3) and the completely deuterated analog pyrazolide-d(3). The resultant spectra, determined employing vibronic Hamiltonians as large as 500 million terms, compare favorably to recent theoretical and spectroscopic results for pyrazolyl-d(3) and to spectroscopic results for pyrazolyl-h(3), for which no reliable simulations had been available.     

Theoretical Study on Photoisomerization of 2-Methylfuran to 3-Methylfuran

1 INTRODUCTION 2-Methylfuran belongs to the basic heteroaromatic compounds relevant to many fields of modern che- mistry, ranging from the study of natural products and biologically active substances to the develop- ment of building blocks for organic synthesis and conducting polymers[1]. Since the photochemistry ofR-furan was gradually recognized in 1960s[2～7], lots of interest has been aroused. Herein we only study one branch of photoche- mistry of R-furan: the isomerization of 2-methy…     

Carbene formation in its lower singlet state from photoexcited 3H-diazirine or diazomethane. A combined CASPT2 and ab initio direct dynamics trajectory study

The potential energy surfaces of the ground and valence excited states of both 3H-diazirine and diazomethane have been studied computationally by mean of the CASSCF method in conjunction with the cc-pVTZ basis set. The energies of the critical points found on such surfaces have been recomputed at the CASPT2/cc-pVTZ level. Additionally, ab initio direct dynamic trajectory calculations have been carried out on the S(1) and S(2) surfaces, starting each trajectory run at the region dominated by the conformational molecular rearrangement of diazomethane. It is found that both isomers are interconnected along a C(s)() reaction coordinate on each potential surface. Radiationless deactivation of the corresponding S(1) state of each isomer occurs through the same point on the surface, an S(1)/S(0) conical intersection. Thereafter, the system has enough energy to surmount the barrier which leads to dissociation products (CH(2) + N(2)) on S(0) state. Therefore, photoexcitation to S(1) state of either diazirine of diazomethane produces methylene in its lower singlet state on a very short time scale (ca. 100 fs). Furthermore, both isomers can generate excited singlet carbene when they are excited onto the S(2) surface; in this case, they lose the activation energy passing through another common S(2)/S(1) conical intersection and then proceed to dissociation into carbene and N(2) on the S(1) surface. For the special case of methylene, it rapidly experiences deexcitation to S(0) state.     

Cyclooctatetraene computational photo- and thermal chemistry: a reactivity model for conjugated hydrocarbons

We use ab initio CASSCF and CASPT2 computations to construct the composite multistate relaxation path relevant to cycloocta-1,3,5,7-tetraene singlet photochemistry. The results show that an efficient population of the dark excited state (S(1)) takes place after ultrafast decay from the spectroscopic excited state (S(2)). A planar D(8)(h)-symmetric minimum represents the collecting point on S(1). Nonadiabatic transitions to S(0) appear to be controlled by two different tetraradical-type conical intersections, which are directly accessible from the S(1) minimum following specific excited-state reaction paths. The higher-energy conical intersection belongs to the same type of intersections previously documented in linear and cyclic conjugated hydrocarbons and features a triangular -(CH)(3)- kink. This point mediates both cis --> trans photoisomerization and cyclopropanation reactions. The lowest energy conical intersection has a boat-shaped structure. This intersection accounts for production of semibullvalene or for double-bond shifting. The mapping of both photochemical and thermal reaction paths (including also Cope rearrangements, valence isomerizations, ring inversions, and double-bond shifting) has allowed us to draw a comprehensive reactivity scheme for cyclooctatetraene, which rationalizes the experimental observations and documents the complex network of photochemical and thermal reaction path interconnections. The factors controlling the selection and accessibility of a number of conjugated hydrocarbon prototype conical intersections and ground-state relaxation channels are discussed.     

A simultaneous eigen valued expression to model the three coupled electronic states with triply and doubly degenerate seams together II

In part I a novel and yet complicated mathematical model to predict the presence of triply degenerate seam was presented and was illustrated its efficiency using some physical models. Here in this part, a discussion on the the model's ability to account for the triple degeneracy found in the potential energy surface is presented in addition to a discussion on the derivatives and its discontinuities and infiniteness at a specific point. Analytical evaluation of derivatives is possible in principle but in practice, numerical approximation to the exact derivatives is the only way out to solve this problem. Two examples (one being physical and the other being chemically realistic) are taken to the probe the above mentioned properties of the proposed simultaneous model. Integration of model function gives finite values everywhere irrespective of the point of degeneracy being present or not. This is just opposite to derivative evaluation.     

Reaction coordinate analysis of the S2-S1 internal conversion of phenylacetylene

The reaction coordinate of the S(2)-S(1) internal conversion (IC) of phenylacetylene (PA) was analyzed using the ab initio complete active space self-consistent field (CASSCF) method. In the first process after electronic excitation into S(2), the aromatic benzene ring is transformed into a nonaromatic quinoid structure. The ethynyl part (-C[triple bond]CH) takes an incomplete allenoid structure in which the CC bond elongates to an intermediate value between typical C[triple bond]C triple and C=C double bonds, but the bend angle of -CCH is 180 degrees . In the second process, PA takes a complete allenoid structure with an out-of-plane location of the beta-H atom (i.e., the H atom of the ethynyl part) and a further elongation of the CC bond so that PA is most stable in S(2) (S(2)-bent). The conical intersection between S(2) and S(1) (S(2)/S(1)-CIX) is located near the S(2)-bent geometry and is slightly unstable energetically. After transition at S(2)/S(1)-CIX, PA quickly loses both quinoid and allenoid structures and recovers the aromaticity of the benzene ring in S(1). Analysis of the dipole moment along the reaction coordinate shows that the weak electron-withdrawing group of the ethynyl part in S(0) suddenly changes into an electron-donating group in S(2) after the main transition of S(0)-S(2). The photoinduced change of the dipole moment is a driving force to the formation of a quinoid structure in S(2). Regarding the benefit of the reaction coordinate analysis of the multidimensional potential energy surfaces of PA, the present picture of the IC process is much more elaborate than our previous representation (Amatatsu, Y.; Hasebe, Y. J. Phys. Chem. A 2003 107, 11169-11173). Vibrational analyses along the reaction coordinate were also performed to support a time-resolved spectroscopic experiment on the S(2)-S(1) IC process of PA.