On the (Emultiply sign in circlee)-Jahn-Teller conical intersections in the 3p(E') and 3d(E") Rydberg electronic states of triatomic hydrogen

The static and dynamic aspects of the Jahn-Teller (JT) interactions in the 3p(E') and 3d(E") Rydberg electronic states of H3 are analyzed theoretically. The static aspects are discussed based on recent ab initio quantum chemistry results, and the dynamic aspects are examined in terms of the vibronic spectra and nonradiative decay behavior of these states. The adiabatic potential-energy surfaces of these degenerate electronic states are derived from extensive ab initio calculations. The calculated adiabatic potential-energy surfaces are diabatized following our earlier study on this system in its 2p(E') ground electronic state. The nuclear dynamics on the resulting conically intersecting manifold of electronic states is studied by a time-dependent wave-packet approach. Calculations are performed both for the uncoupled and coupled state situations in order to understand the importance of nonadiabatic interactions due to the JT conical intersections in these excited Rydberg electronic states.     

On the simulation of photoelectron spectra in molecules with conical intersections and spin-orbit coupling: the vibronic spectrum of CH3S

A method to simulate photoelectron spectra for states coupled by conical intersections and the spin-orbit interaction is reported. The algorithm is based on the multimode vibronic coupling model and treats the spin-orbit interaction in a nonperturbative manner. Since the algorithm is not dependent on molecular symmetry, the approach is generally applicable to accidental conical intersections as well as the symmetry required intersections found in Jahn-Teller molecules. The method is also computationally efficient using energy gradient and derivative coupling information to limit the number of nuclear configurations at which ab initio data are required. This approach is applied to simulate the negative ion photoelectron spectrum of the methylthio radical. The two-state Hamiltonian employed to describe this system was determined employing ab initio gradients and derivative couplings at only 17 nuclear configurations.     

The Jahn-Teller effect in the triply degenerate electronic state of methane radical cation

A quantum dynamics study is performed to examine the complex nuclear motion underlying the first photoelectron band of methane. The broad and highly overlapping structures of the latter are found to originate from transitions to the ground electronic state, X(2)T(2), of the methane radical cation. Ab initio calculations have also been carried out to establish the potential energy surfaces for the triply degenerate electronic manifold of CH(4)(+). A suitable diabatic vibronic Hamiltonian has been devised and the nonadiabatic effects due to Jahn-Teller conical intersections on the vibronic dynamics investigated in detail. The theoretical results show fair accord with experiment.     

Vibronic interactions in the photodetachment spectroscopy of phenide anion

Photodetachment spectroscopy of phenide anion C6H5- is theoretically studied with the aid of electronic structure calculations and quantum dynamical simulations of nuclear motion. The theoretical results are compared with the available experimental data. The vibronic structure of the first, second, and third photoelectron bands associated with the ground X 2A1 and low-lying excited A 2B1 and B 2A2 electronic states of the phenyl radical C6H5 is examined at length. While the X state of the radical is energetically well separated and its interaction is found to be rather weak with the rest, the A and B electronic states are found to be only approximately 0.57 eV apart in energy at the vertical configuration. Low-energy conical intersections between the latter two states are discovered and their impact on the nuclear dynamics underlying the second and third photoelectron bands is delineated. The nuclear dynamics in the X state solely proceeds through the adiabatic path and the theoretically calculated vibrational level structure of this state compares well with the experimental result. Two Condon active totally symmetric (a1) vibrational modes of ring deformation type form the most dominant progression in the first photoelectron band. The existing ambiguity in the assignment of these two vibrational modes is resolved here. The A-B conical intersections drive the nuclear dynamics via nonadiabatic paths, and as a result the second and third photoelectron bands overlap and particularly the third band due to the B state of C6H5 becomes highly diffused and structureless. Experimental photodetachment spectroscopy results are not available for these bands. However, the second band has been detected in electronic absorption spectroscopy measurements. The present theoretical results are compared with these absorption spectroscopy data to establish the nonadiabatic interactions between the A and B electronic states of C6H5.     

Vibronic coupling in the excited cationic states of ethylene: simulation of the photoelectron spectrum between 12 and 18 eV

The effect of vibronic coupling on structure and spectroscopy is investigated in the excited cationic states of ethylene. It is found from equation of motion coupled cluster singles and doubles method for ionization potential electronic structure calculations in a triple-zeta plus double polarization basis set that ethylene in its third (B (2)A(g)) and fourth (C (2)B(2u)) ionized states does not have a stable minimum-energy geometry. The potential-energy surfaces of these states are energetically distinct and well separated at the ground-state geometry of ethylene, but in a geometry optimization as the structure of the ion relaxes, these surfaces end up in conical intersections and finally in the stable equilibrium geometry of the second ionized state (A (2)B(3g)). The topology of the potential-energy surfaces can be clearly understood using a vibronic model Hamiltonian. Furthermore, by diagonalizing this model Hamiltonian, the photoelectron spectrum of ethylene corresponding to the second, third, and fourth ionized states (12-18 eV) is simulated. Spectra from vibronic simulations including up to quartic coupling constants and using various normal-mode basis sets are compared to those from vertical Franck-Condon simulations to understand the importance of vibronic coupling and nonadiabatic effects and to examine the influence of individual normal modes on the spectrum.     

Multimode Jahn-Teller and pseudo-Jahn-Teller interactions in the cyclopropane radical cation: complex vibronic spectra and nonradiative decay dynamics

The complex vibronic spectra and the nonradiative decay dynamics of the cyclopropane radical cation (CP+) are simulated theoretically with the aid of a time-dependent wave packet propagation approach using the multireference time-dependent Hartree scheme. The theoretical results are compared with the experimental photoelectron spectrum of cyclopropane. The ground and first excited electronic states of CP+ are of X2E' and A2E' type, respectively. Each of these degenerate electronic states undergoes Jahn-Teller (JT) splitting when the radical cation is distorted along the degenerate vibrational modes of e' symmetry. The JT split components of these two electronic states can also undergo pseudo-Jahn-Teller (PJT)-type crossings via the vibrational modes of e', a1' and a2' symmetries. These lead to the possibility of multiple multidimensional conical intersections and highly nonadiabatic nuclear motions in these coupled manifolds of electronic states. In a previous publication [J. Phys. Chem. A 2004, 108, 2256], we investigated the JT interactions alone in the X2E' ground electronic manifold of CP+. In the present work, the JT interactions in the A2E' electronic manifold are treated, and our previous work is extended by considering the coupling between the X2E' and A2E' electronic states of CP+. The nuclear dynamics in this coupled manifold of two JT split doubly degenerate electronic states is simulated by considering fourteen active and most relevant vibrational degrees of freedom. The vibronic level spectra and the ultrafast nonradiative decay of the excited cationic states are examined and are related to the highly complex entanglement of electronic and nuclear degrees of freedom in this prototypical molecular system.     

An unusually large nonadiabatic error in the BNB molecule

The vibronic coupling model of Ko?uppel, Domcke, and Cederbaum in one dimension is introduced as a means to estimate the effects of electronic nonadiabaticity on the vibrational energy levels of molecules that exhibit vibronic coupling. For the BNB molecule, the nonadiabatic contribution to the nominal fundamental vibrational energy of the antisymmetric stretching mode is approximately -80?cm(-1). The surprisingly large effect for this mode, which corresponds to an adiabatic potential that is essentially flat near the minimum due to the vibronic interaction, is contrasted with another model system that also exhibits a flat potential (precisely, a vanishing quadratic force constant) but has a significantly larger gap between interacting electronic states. For the latter case, the nonadiabatic contribution to the level energies is about two orders of magnitude smaller even though the effect on the potential is qualitatively identical. A simple analysis shows that significant nonadiabatic corrections to energy levels should occur only when the affected vibrational frequency is large enough to be of comparable magnitude to the energy gap involved in the coupling. The results provide evidence that nonadiabatic corrections should be given as much weight as issues such as high-level electron correlation, relativistic corrections, etc., in quantum chemical calculations of energy levels for radicals with close-lying and strongly coupled electronic states even in cases where conical intersections are not obviously involved. The same can be said for high-accuracy thermochemical studies, as the zero-point vibrational energy of the BNB example contains a nonadiabatic contribution of approximately -70?cm(-1) (-0.9?kJ mol(-1)).     

On the determination of partial differential cross sections for photodetachment and photoionization processes producing polyatomic molecules with electronic states coupled by conical intersections

A formalism is derived for the computation of partial differential cross sections for electron photodetachment and photoionization processes that leave the residual or target molecule in electronic states that are strongly coupled by conical intersections. Because the electronic states of the target are nonadiabatically coupled, the standard adiabatic states approach of solving the electronic Schro?dinger equation for the detached electron at fixed nuclear geometries and then vibrationally averaging must be fundamentally modified. We use a Lippmann-Schwinger equation based approach, which leads naturally to a partitioning of the transition amplitude into a Dyson orbital like part plus a scattering correction. The requisite Green's function is that developed in our previous paper for the direct determination of total integral cross sections. The method takes proper account of electron exchange, possible nonorthogonality of the orbital describing the detached electron, and nonadiabatic effects in the product molecule. The Green's function is constructed in an L(2) basis using complex scaling techniques. The accurate treatment of nonadiabatic effects in the residual molecule is accomplished using the multimode vibronic coupling model. For photodetachment, an approximate approach, which is less computationally demanding, is suggested.     

On the construction of quasidiabatic state representations of bound adiabatic state potential energy surfaces coupled by accidental conical intersections: incorporation of higher order terms

The quadratic vibronic coupling model is an important computational tool for simulating photoelectron spectra involving strongly coupled electronic states in polyatomic molecules. However, recent work has indicated the need for higher order terms, with most of the initial studies focusing on molecules with symmetry-required degeneracies. In this study we report an extension of our approach for constructing fully quadratic representations of bound electronic states coupled by conical intersections, which allows for the inclusion of higher order terms, demonstrated here employing a quartic expansion. Procedures are developed that eliminate unphysical behavior for large displacements, a problem likely to be an endemic to anharmonic expansions. Following work on representing dissociative electronic states, Lagrange multipliers are used to constrain the constructed representation to reproduce exactly the energy, energy gradients, and∕or derivative couplings at specific points, or nodes, in nuclear coordinate space. The approach is illustrated and systematically studied using the four lowest electronic states of triazolyl, (CH)(2)N(3).     

A proposal for the theoretical treatment of multi-electronic-state molecular dynamics: Hemiquantal dynamics with the whole dim basis (HWD). A test on the evolution of excited Ar3 cluster ions

A method which couples a model Hamiltonian to hemiquantal dynamics is proposed as a general theoretical treatment of multistate molecular processes. A test application to the Ar₃⁺ system is performed which involves four nuclear degrees of freedom and six quasidiabatic states. The method proves efficient for dealing with non-adiabatic transitions occurring at singularities like conical intersections and avoided crossings between potential surfaces.     

Photochemistry of pyrrole: time-dependent quantum wave-packet description of the dynamics at the 1pi sigma*-S0 conical intersections

The photoinduced hydrogen-elimination reaction in pyrrole via the conical intersections of the two (1)pi sigma(*) excited states with the electronic ground states [(1)B(1)(pi sigma(*))-S(0) and (1)A(2)(pi sigma(*))-S(0)] have been investigated by time-dependent quantum wave-packet calculations. Model potential-energy surfaces of reduced dimensionality have been constructed on the basis of accurate multireference ab initio electronic-structure calculations. For the (1)B(1)-S(0) conical intersection, the model includes the NH stretching coordinate as the tuning mode and the hydrogen out-of-plane bending coordinate as the coupling mode. For the (1)A(2)-S(0) conical intersection, the NH stretching coordinate and the screwing coordinate of the ring hydrogens are taken into account. The latter is the dominant coupling mode of this conical intersection. The electronic population-transfer processes at the conical intersections, the branching ratio between the dissociation channels, and their dependence on the initial preparation of the system have been investigated for pyrrole and deuterated pyrrole. It is shown that the excitation of the NH stretching mode strongly enhances the reaction rate, while the excitation of the coupling mode influences the branching ratio of different dissociation channels. The results suggest that laser control of the photodissociation of pyrrole via mode-specific vibrational excitation should be possible. The calculations provide insight into the microscopic details of ultrafast internal-conversion processes in pyrrole via hydrogen-detachment processes, which are aborted at the (1)pi sigma(*)-S(0) conical intersections. These mechanisms are of relevance for the photostability of the building blocks of life (e.g., the DNA bases).     

Theoretical study of excitations in furan: spectra and molecular dynamics

The excitation spectra and molecular dynamics of furan associated with its low-lying excited singlet states 1A2(3s), 1B2(V), 1A1(V'), and 1B1(3p) are investigated using an ab initio quantum-dynamical approach. The ab initio results of our previous work [J. Chem. Phys. 119, 737 (2003)] on the potential energy surfaces (PES) of these states indicate that they are vibronically coupled with each other and subject to conical intersections. This should give rise to complex nonadiabatic nuclear dynamics. In the present work the dynamical problem is treated using adequate vibronic coupling models accounting for up to four coupled PES and thirteen vibrational degrees of freedom. The calculations were performed using the multiconfiguration time-dependent Hartree method for wave-packet propagation. It is found that in the low-energy region the nuclear dynamics of furan is governed mainly by vibronic coupling of the 1A2(3s) and 1B2(V) states, involving also the 1A1(V') state. These interactions are responsible for the ultrafast internal conversion from the 1B2(V) state, characterized by a transfer of the electronic population to the 1A2(3s) state on a time scale of approximately 25 fs. The calculated photoabsorption spectrum of furan is in good qualitative agreement with experimental data. Some assignments of the measured spectrum are proposed.     

Nonequilibrium Fermi golden rule for electronic transitions through conical intersections

We consider photoinduced electronic transitions through conical intersections in large molecules. Starting from the linear vibronic model Hamiltonian and treating linear diabatic couplings within the second order cumulant expansion, we have developed a simple analytical expression for the time evolution of electronic populations at finite temperature. The derived expression can be seen as a nonequilibrium generalization of the Fermi golden rule due to a nonequilibrium character of the initial photoinduced nuclear distribution. All parameters in our model are obtained from electronic structure calculations followed by a diabatization procedure. The results of our model are found to agree well with those of quantum dynamics for a test set of systems: fulvene molecule, 2,6-bis(methylene) adamantyl cation, and its dimethyl derivative.     

Communication: multistate quantum dynamics of photodissociation of carbon dioxide between 120 nm and 160 nm

UV absorption cross section of CO(2) is studied using high level ab initio quantum chemistry for electrons and iterative quantum dynamics for nuclear motion on interacting global full dimensional potential energy surfaces. Six electronic states-1, 2, 3(1)A(') and 1, 2, 3(1)A(")-are considered. At linearity, they correspond to the ground electronic state X?(1)Σ(g) (+) and the optically forbidden but vibronically allowed valence states 1(1)Δ(u), 1(1)Σ(u) (-), and 1(1)Π(g). In the Franck-Condon region, these states interact via Renner-Teller and conical intersections and are simultaneously involved in an intricate network of non-adiabatic couplings. The absorption spectrum, calculated for many rotational states, reproduces the distinct two-band shape of the experimental spectrum measured at 190 K and the characteristic patterns of the diffuse structures in each band. Quantum dynamics unravel the relative importance of different vibronic mechanisms, while metastable resonance states, underlying the diffuse structures, provide dynamically based vibronic assignments of individual lines.     

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.     

Vibrationally mediated photodissociation of ethylene cation by reflectron multimass velocity map imaging

A new imaging technique, reflectron multimass velocity map ion imaging, is used to study the vibrationally mediated photodissociation dynamics in the ethylene cation. The cation ground electronic state is prepared in specific vibrational levels by two-photon resonant, three-photon ionization via vibronic bands of (pi, nf) Rydberg states in the vicinity of the ionization potential of ethylene, then photodissociated through the (B 2A(g)) excited state. We simultaneously record spatially resolved images of parent C2H4+ ions as well as photofragment C2H3+ and C2H2+ ions originating in dissociation from the vibronic excitations in two distinct bands, 7f 4(0)2 and 8f 0(0)0, at roughly the same total energy. By analyzing the images, we directly obtain the total translation energy distributions for the two dissociation channels and the branching between them. The results show that there exist differences for competitive dissociation pathways between H and H2 elimination from C2H4+ depending on the vibronic preparation used, i.e., on the vibrational excitation in the ground state of the cation prior to photodissociation. Our findings are discussed in terms of the possible influence of the torsional excitation on competition between direct dissociation, isomerization, and radiationless transitions through conical intersections among the numerous electronic states that participate in the dissociation.     

Conical intersections and semiclassical trajectories: comparison to accurate quantum dynamics and analyses of the trajectories

Semiclassical trajectory methods are tested for electronically nonadiabatic systems with conical intersections. Five triatomic model systems are presented, and each system features two electronic states that intersect via a seam of conical intersections (CIs). Fully converged, full-dimensional quantum mechanical scattering calculations are carried out for all five systems at energies that allow for electronic de-excitation via the seam of CIs. Several semiclassical trajectory methods are tested against the accurate quantum mechanical results. For four of the five model systems, the diabatic representation is the preferred (most accurate) representation for semiclassical trajectories, as correctly predicted by the Calaveras County criterion. Four surface hopping methods are tested and have overall relative errors of 40%-60%. The semiclassical Ehrenfest method has an overall error of 66%, and the self-consistent decay of mixing (SCDM) and coherent switches with decay of mixing (CSDM) methods are the most accurate methods overall with relative errors of approximately 32%. Furthermore, the CSDM method is less representation dependent than both the SCDM and the surface hopping methods, making it the preferred semiclassical trajectory method. Finally, the behavior of semiclassical trajectories near conical intersections is discussed.     

Remarkable impact of intermode couplings on multimode vibronic dynamics: the photoelectron spectrum of CH3F

Electronic and nuclear motions on intersecting potential energy surfaces are often intricately mixed and the spectrum can become very complex. Here we choose the strongly coupled Jahn-Teller system CH3F+ as a prototype example, and establish the importance of intermode coupling terms on multimode vibronic dynamics. The theoretical approach consists of a full second-order diabatic vibronic Hamiltonian, constructed from high-quality electronic structure calculations. Our results compare amazingly well with the experimental data. This highlights the success of the present theoretical approach in explaining the complex structure of vibronic spectra, ubiquitous in molecular systems.     

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.     

Electronic degeneracies in symmetric (Jahn-Teller) and nonsymmetric aliphatic radical cations: global topology of sigma-bonded molecules

The ground-state potential surfaces of five aliphatic radical cations are investigated using a spin-pairing model. It is shown that the ground-state surface of an n-atomic system supports several stationary points (minima and transition states, including second-order ones). In addition, there are numerous nuclear configurations at which the ground state is electronically degenerate. The electronic degeneracies due to interactions between atoms bound to the same atom are either 2-fold (conical intersections) or 3-fold degenerate but not of a higher dimension. Each 3-fold degeneracy is accompanied by an even number of conical intersections (four or two). A systematic procedure for locating all of these nuclear configurations (that are in fact 3n - 8 or 3n - 11 dimensional hypersurfaces) is described. The model allows for the qualitative determination of the structure and charge distribution of the system at all of the stationary points and electronic degeneracies. Quantum chemical calculations confirm the predictions of the model, which is used to direct and facilitate the calculations.