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.     

On the vibronic coupling approximation: a generally applicable approach for determining fully quadratic quasidiabatic coupled electronic state Hamiltonians

In this report we introduce an iterative procedure for constructing a quasidiabatic Hamiltonian representing N(state)-coupled electronic states in the vicinity of an arbitrary point in N(int)-dimensional nuclear coordinate space. The Hamiltonian, which is designed to compute vibronic spectra employing the multimode vibronic coupling approximation, includes all linear terms which are determined exactly using analytic gradient techniques. In addition, all [N(state)][N(int)] quadratic terms, where [n]=n(n+1)/2, are determined from energy gradient and derivative coupling information obtained from reliable multireference configuration interaction wave functions. The use of energy gradient and derivative coupling information enables the large number of second order parameters to be determined employing ab initio data computed at a limited number of points (N(int) being minimal) and assures a maximal degree of quasidiabaticity. Numerical examples are given in which quasidiabatic Hamiltonians centered around three points on the C(3)H(3)N(2) potential energy surface (the minimum energy point on the ground state surface and the minimum energy points on the two- and three-state seams of conical intersection) were computed and compared. A method to modify the conical intersection based Hamiltonians to better describe the region of the ground state minimum is introduced, yielding improved agreement with ab initio results, particularly in the case of the Hamiltonian defined at the two-state minimum energy crossing.     

On the characterization of three state conical intersections using a group homomorphism approach: mapping the full N-5 dimensional seam space

A method for characterizing the degeneracy preserving seam space in the vicinity of a three state conical intersection is introduced. Second order degenerate perturbation theory is used to construct an approximately diabatic Hamiltonian whose eigenenergies and eigenstates accurately describe the vicinity of the three state conical intersection in its full dimensionality. The perturbative analysis enables the large number, 6(N(int)(N(int)+1)2), of unique second order parameters needed to construct this accurate Hamiltonian to be determined from ab initio data at a limited number of nuclear configurations, with (N(int)+10) being minimal. Using the minimum energy three state conical intersection of the pyrazolyl radical (N(int) = 18), the potential of this approach is illustrated. A Hamiltonian comprised of the ten characteristic (linear) parameters and over 1440 second order parameters is constructed and used to determine the locus of the conical intersection seam as well as to describe the 18 dimensional space in the vicinity of that point of intersection. Our results demonstrate the ability of this methodology to quantitatively reproduce the ab initio potential energy surfaces near a three state conical intersection.     

Theoretical analysis of photoinduced H-atom elimination in thiophenol

The photoinduced hydrogen elimination reaction in thiophenol via the conical intersections of the dissociative (1)πσ? excited state with the bound (1)ππ? excited state and the electronic ground state has been investigated with ab initio electronic-structure calculations and time-dependent quantum wave-packet calculations. A screening of the coupling constants of the symmetry-allowed coupling modes at the (1)ππ?-(1)πσ? and (1)πσ?-S(0) conical intersection shows that the SH torsional mode is by far the most important coupling mode at both conical intersections. A model including three intersecting potential-energy surfaces (S(0), (1)ππ?, (1)πσ?) and two nuclear degrees of freedom (SH stretch and SH torsion) has been constructed on the basis of ab initio complete-active-space self-consistent field and multireference second-order perturbation theory calculations. The nonadiabatic quantum wave-packet dynamics initiated by optical excitation of the (1)ππ? and (1)πσ? states has been explored for this three-state two-coordinate model. The photodissociation dynamics is characterized in terms of snapshots of time-dependent wave packets, time-dependent electronic population probabilities, and the branching ratio of the (2)σ/(2)π electronic states of the thiophenoxyl radical. The dependence of the timescale of the photodissociation process and the branching ratio on the initial excitation of the SH stretching and SH torsional vibrations has been analyzed. It is shown that the node structure, which is imposed on the nuclear wave packets by the initial vibrational preparation as well as by the transitions through the conical intersections, has a profound effect on the photodissociation dynamics. The effect of additional weak coupling modes of CC twist (ν(16a)) and ring-distortion (ν(16b)) character has been investigated with three-dimensional and four-dimensional time-dependent wave-packet calculations, and has been found to be minor.     

On the characterization of three state conical intersections: a quasianalytic theory using a group homomorphism approach

In this work, degenerate perturbation theory through second order is used to characterize the vicinity of a three state conical intersection. This report extends our recent demonstration that it is possible to describe the branching space (in which the degeneracy is lifted linearly) and seam space (in which the degeneracy is preserved) in the vicinity of a two state conical intersection using second order perturbation theory. The general analysis developed here is based on a group homomorphism approach. Second order perturbation theory, in conjunction with high quality ab initio electronic structure data, produces an approximately diabatic Hamiltonian whose eigenenergies and eigenstates can accurately describe the three adiabatic potential energy surfaces, the interstate derivative couplings, and the branching and seam spaces in their full dimensionality. The application of this approach to the minimum energy three state conical intersection of the pyrazolyl radical demonstrates the potential of this method. A Hamiltonian comprised of the ten characteristic (linear) parameters and over 300 second order parameters is constructed to describe the branching space associated with a point of conical intersection. The second order parameters are determined using data at only 30 points. In the vicinity of the conical intersection the energy and derivative couplings are well reproduced and the singularity in the derivative coupling is analyzed.     

A two-dimensional wavepacket study of the nonadiabatic dynamics of CH2BrCl

The nonadiabatic photodissociation dynamics of CH2BrCl into CH2Br + Cl or CH2Cl + Br is studied using two-dimensional wavepacket propagations on ab initio multiconfigurational MS-CASPT2 potential energy surfaces. Using a three-state diabatic model, we investigate the electronic states responsible for the two competing fragmentation channels and how the conical intersection present between the two lowest excited states affects the dissociation rate. Within this model, we find that the Br/Cl branching ratio depends on the irradiation wavelength. Predominant C-Br fragmentation occurs for wavelengths longer than 200 nm, while nonadiabatic C-Cl dissociation with a constant branching ratio of 0.4 is predicted upon absorption of photons in the range of 170-180 nm. Additionally, we observe complete nonadiabatic population transfer in less than 100 fs, that is, before the wavepacket can reach the conical intersection. As a consequence, there is no three-body CH2 + Br + Cl dissociation.     

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.     

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.     

Intramolecular electron transfer in bis(methylene) adamantyl radical cation: a case study of diabatic trapping

Our characterization of the potential energy surface for electron transfer (ET) in the bis(methylene)adamantane (BMA) model radical cation shows that the surface topology is prone to diabatic trapping (competition between ET and upward hops to the excited state). The general conditions for this phenomenon have been derived. The surface is centered around a conical intersection, and diabatic trapping occurs because one of the branching space coordinates (coordinates that lift the degeneracy at first order) corresponds to a vector of small length. For BMA, this coordinate is an antisymmetric breathing mode of the rigid carbon framework. Other modes (including methylene torsions and pyramidalizations) may lift the degeneracy at second-order but do not affect the energy gap at the intersection region effectively. The resulting topology is similar to that of an (n - 1) dimensional seam (where n is the number of nuclear degrees of freedom of the molecule) that cannot be avoided along the reaction coordinate, thus favoring recrossing to the upper surface. This analysis is extended by ab initio semiclassical dynamics using an Ehrenfest and a trajectory surface hopping algorithm implemented at the CASSCF level. Examination of the trajectories shows that there is no single mode that controls the diabatic trap, in agreement with the condition that there is no predominant degeneracy-lifting coordinate. Thus the reactivity depends on a combination of small effects, where presumably higher-order effects come into play. This should be the general behavior of dynamics at a diabatic trapping situation.     

A three-state model for the photophysics of adenine

An ab initio theoretical study at the CASPT2 level is reported on minimum energy reaction paths, state minima, transition states, reaction barriers, and conical intersections on the potential energy hypersurfaces of two tautomers of adenine: 9H- and 7H-adenine. The obtained results led to a complete interpretation of the photophysics of adenine and derivatives, both under jet-cooled conditions and in solution, within a three-state model. The ultrafast subpicosecond fluorescence decay measured in adenine is attributed to the low-lying conical intersection (gs/pipi* La)(CI), reached from the initially populated 1(pipi* La) state along a path which is found to be barrierless only in 9H-adenine, while for the 7H tautomer the presence of an intermediate plateau corresponding to an NH2-twisted conformation may explain the absence of ultrafast decay in 7-substituted compounds. A secondary picosecond decay is assigned to a path involving switches towards two other states, 1(pipi* Lb) and 1(npi*), ultimately leading to another conical intersection with the ground state, (gs/npi*), with a perpendicular disposition of the amino group. The topology of the hypersurfaces and the state properties explain the absence of secondary decay in 9-substituted adenines in water in terms of the higher position of the 1(npi*) state and also that the 1(pipi* Lb) state of 7H-adenine is responsible for the observed fluorescence in water. A detailed discussion comparing recent experimental and theoretical findings is given. As for other nucleobases, the predominant role of a pipi*-type state in the ultrafast deactivation of adenine is confirmed.     

On the representation of coupled adiabatic potential energy surfaces using quasi-diabatic Hamiltonians: a distributed origins expansion approach

In two previous papers we have introduced a method to generate coupled quasi-diabatic Hamiltonians (H(d)) that are capable of representing adiabatic energies, energy gradients, and derivative couplings over a wide range of geometries including seams of conical intersection. In this work, two new synergistic features are introduced. Firstly, the functional form of H(d) is generalized. Rather than requiring there to be a low energy point of high symmetry to serve as the unique origin, functions centered on points distributed in nuclear coordinate space are used in the polynomials that comprise the matrix elements in H(d). The use of functions with distributed origins, allows reproduction of the ab initio data with lower order expansions, and offers the possibility of describing multichannel dissociation. The fitting algorithm is combined with a three-step procedure in which the domain of H(d) is extended from a core set of nuclear configurations to a region of nuclear coordinate space appropriate for nuclear dynamics, with a prescribed accuracy. This significant extension of the domain of definition compared to our original work, which is facilitated by the distributed origin approach, is achieved largely through the use of surface hopping trajectories. The 1,2(1)A states of NH(3), which provide an archetypical example of nonadiabatic dynamics, are used to demonstrate the utility of this approach. The representation describes 21 points on the 1(1)A-2(1)A seam of conical intersection and their local topography flawlessly and on the entire domain, the electronic structure data is represented to an accuracy of 77.00 (46.90) cm(-1), as measured by the root mean square (mean unsigned) error for energies lower than 50 000 cm(-1). This error is a factor of 10 lower than that of the most accurate representation of high quality ab initio data, on a comparable domain, previously reported for this system.     

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.     

Relationship between photoisomerization path and intersection space in a retinal chromophore model

A low-lying segment of the intersection space (IS) between the excited-state and the ground-state energy surfaces of a retinal chromophore model has been mapped using ab initio CASSCF computations. Analysis of the structural relationship between the computed IS cross-section and the excited state Z --> E isomerization path shows that these are remarkably close both in energy and in structure. Indeed, the IS segment and the Z --> E path remain roughly parallel and merge only when the double bond reaches a 70 degree twisting. This finding supports the idea that, in certain chromophores, a more extended segment of IS, and not a single conical intersection, contributes to the decay to the ground state.     

Cyclic-N3. II. Significant geometric phase effects in the vibrational spectra

An accurate theoretical prediction of the vibrational spectra for a pure nitrogen ring (cyclic-N(3)) molecule is obtained up to the energy of the (2)A(2)/(2)B(1) conical intersection. A coupled-channel approach using the hyperspherical coordinates and the recently published ab initio potential energy surface [D. Babikov, P. Zhang, and K. Morokuma, J. Chem. Phys. 121, 6743 (2004)] is employed. Two independent sets of calculations are reported: In the first set, the standard Born-Oppenheimer approximation is used and the geometric phase effects are totally neglected. In the second set, the generalized Born-Oppenhimer approximation is used and the geometric phase effects due to the D(3h) conical intersection are accurately treated. All vibrational states are analyzed and assigned in terms of the normal vibration mode quantum numbers. The magnitude of the geometric phase effect is determined for each state. One important finding is an unusually large magnitude of the geometric phase effects in the cyclic-N(3): it is approximately 100 cm(-1) for the low-lying vibrational states and exceeds 600 cm(-1) for several upper states. On average, this is almost two orders of magnitude larger than in the previously reported studies. This unique example suggests a favorable path to experimental validation.     

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.     

Energetics and kinetics of the reaction of HOCO with hydrogen atoms

The potential energy surface for the reaction of HOCO radicals with hydrogen atoms has been explored using the CCSD(T)/aug-cc-pVQZ ab initio method. Results show that the reaction occurs via a formic acid (HOC(O)H) intermediate, and produces two types of products: H(2)O+CO and H(2)+CO(2). Reaction enthalpies (0 K) are obtained as -102.0 kcalmol for the H(2)+CO(2) products, and -92.7 kcalmol for H(2)O+CO. Along the reaction pathways, there exists a nearly late transition state for each product channel. However, the transition states locate noticeably below the reactant asymptote. Direct ab initio dynamics calculations are also carried out for studying the kinetics of the H+HOCO reaction. At room temperature, the rate coefficient is predicted to be 1.07x10(-10)cm(3) molec(-1) s(-1) with a negligible activation energy E(a)=0.06 kcalmol, and the branching ratios are estimated to be 0.87 for H(2)+CO(2), and 0.13 for H(2)O+CO. In contrast, the product branching ratios have a strong T dependence. The branching ratio for H(2)O+CO could increase to 0.72 at T=1000 K.     

Conical intersection of the ground and first excited states of water: energies and reduced density matrices from the anti-Hermitian contracted Schrödinger equation

A conical intersection between the ground and first-excited states of water is computed through the direct calculation of two-electron reduced density matrices (2-RDMs) from solutions of the anti-Hermitian contracted Schr?dinger equation (ACSE). This study is an extension of a previous study in which the ACSE was used to compute the energies around a conical intersection in the triplet excited states of methylene [Snyder, J. W., Jr.; Rothman, A. E.; Foley, J. J.; Mazziotti, D. A. J. Chem. Phys. 2010, 132, 154109]. We compute absolute energies of the 1(1)A' and 2(1)A' states of water (H(2)O) and the location of the conical intersection. The ACSE energies are compared to those from ab initio wave function methods. To treat multireference correlation, we seed the ACSE with an initial 2-RDM from a multiconfiguration self-consistent field (MCSCF) calculation. Unlike the situation for methylene, the two states in the vicinity of the conical intersection of water both have the same spatial symmetry. Hence, the study demonstrates the ability of the ACSE to resolve states of the same spatial symmetry that are nearly degenerate in energy. The 2-RDMs from the ACSE nearly satisfy necessary N-representability conditions. Comparison of the results from double-ζ and augmented double-ζ basis sets demonstrates the importance of augmented (or diffuse) functions for determining the location of the conical intersection.     

Ab initio nonadiabatic quantum dynamics of cyclohexadiene/hexatriene ultrafast photoisomerization

Reaction mechanisms of the ultrafast photoisomerization between cyclohexadiene and hexatriene have been elucidated by the quantum dynamics on the ab initio potential energy surfaces calculated by multireference configuration interaction method. In addition to the quantum wave-packet dynamics along the two-dimensional reaction coordinates, the semiclassical analyses have also been carried out to correctly estimate the nonadiabatic transition probabilities around conical intersections in the full-dimensional space. The reaction time durations of radiationless decays in the wave-packet dynamics are found to be generally consistent with the femtosecond time-resolution experimental observations. The nonadiabatic transition probabilities among the ground (S0), first (S1), and second (S2) excited states have been estimated by using the semiclassical Zhu-Nakamura formula considering the full-dimensional wave-packet density distributions in the vicinity of conical intersections under the harmonic normal mode approximation. The cyclohexadiene (CHD) ring-opening process proceeds descending on the S1(1 1B) potential after the photoexcitation. The major part of the wave-packet decays from S1(1 1B) to S1(2 1A) by the first seam line crossing along the C2-symmetry-breaking directions. The experimentally observed ultrafast S1-S0 decay can be explained by the dynamics through the S1-S0 conical intersection along the direction toward the five-membered ring. The CHD: hexatriene (HT) branching ratio is estimated to be approximately 5:5, which is in accordance with the experiment in solution. This branching ratio is found to be mainly governed by the location of the five-membered ring S1-S0 conical intersection along the ground state potential ridge between CHD and HT.     

改进的半经典动力学模拟二苯乙烯光致顺反异构化反应

发展了一种改进的半经典动力学模拟方法, 并将其程序化用于气相二苯乙烯光致顺反异构化反应的机理研究.新的方法不仅采用e 指数模型改进了原有Zhu-Nakamura 理论中计算电子非绝热跃迁几率的计算方法, 而且将约束哈密顿方法用于限制性分子动力学模拟过程中. 计算结果表明, 采用此方法得到的统计平均的量子产率及反应机理与以前的实验与理论结果吻合较好, 从而可以应用于全量子动力学方法无法进行的大分子体系的动力学研究.     

Competitive decay at two- and three-state conical intersections in excited-state intramolecular proton transfer

We demonstrate the existence of a simultaneous degeneracy (not required by symmetry) of three electronic states in malonaldehyde. This is one of the first reports of such a triple degeneracy involving S0, S1, and S2 in a molecule with a closed-shell ground state. We further report on a two-state S2/S1 conical intersection which is higher in energy than the three-state intersection, but closer to the Franck-Condon point. First-principles quantum dynamics calculations of the photochemistry after excitation to S2 show that there is a competition between these intersections, with more than half of the population decaying to S1 through the higher energy S2/S1 intersection. Surprisingly, much of the population which makes it to the triple degeneracy point is not funneled directly to S0, but rather remains trapped on S1. We attribute this to the large dimensionality of the branching plane at a three-state intersection (the degeneracy is lifted along at least five distinct molecular displacements).