Nonadiabatic coupling and diabatic electronic population dynamics on 11A2 and 11B1 states of ozone molecule

In this work, we examine nonadiabatic population dynamics for 1¹B₁ and 1¹A₂ states of ozone molecule (O₃). In O₃, two lowest singlet excited states, ¹A₂ and ¹B₁, can be coupled. Thus, population transfer between them occurs through the seam involving these two states. At any point of the seam (conical intersection), the Born-Oppenheimer approximation breaks down, and it is necessary to investigate nonadiabatic dynamics. We consider a linear vibronic coupling Hamiltonian model and evaluate vibronic coupling constant, diabatic frequencies for three modes of O₃, bilinear and quadratic coupling constants for diabatic potentials, displacements, and Huang-Rhys coupling constants using ab initio calculations. The electronic structure calculations have been performed at the multireference configuration interaction and complete active space with second-order perturbation theory with a full-valence complete active space self-consistent field methods and augmented Dunning's standard correlation-consistent-polarized quadruple zeta basis set to determine ab initio potential energy surfaces for the ground state and first two excited states of O₃, respectively. We have chosen active space comprising 18 electrons distributed over 12 active orbitals. Our calculations predict the linear vibronic coupling constant 0.123 eV. We have obtained the population on the 1¹B₁ and 1¹A₂ excited electronic states for the first 500 fs after photoexcitation.     

Adiabatic versus diabatic descriptions of the lowest Rydberg and valence 1Σ+ states of HCl

In this contribution we first report new ab initio self-consistent field configuration interaction calculations of the first excited adiabatic potential of (1)Σ(+) symmetry, the 2(1)Σ(+) or B(1)Σ(+) state, which presents two minima and can thus be seen as made up of the Rydberg E(1)Σ(+) and the valence V(1)Σ(+) states. Based on the computed 2(1)Σ(+) potential, we devised a theoretical procedure to compute the vibronic structure in order to try to explain the energy levels observed in the region above 76 254.4 cm(-1) which display an irregular vibrational structure, indicative of spectral perturbations. We try to find out which representation of the electronic states, the diabatic or the adiabatic one, is best suited to replicate the lowest observed vibronic levels of the E and V states. To this end, we deduce, from the 2(1)Σ(+) potential and its complementary adiabatic potential, two diabatic potentials. We then carry out a coupled equation treatment based on these diabatic potentials. The results of this treatment indicate that, in the present case, the adiabatic representation is better than the diabatic one to describe the observed vibronic levels. This is due, as expected, to the existence of a strong electrostatic interaction between the two diabatic potentials.     

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)).     

Proton-coupled electron transfer versus hydrogen atom transfer: generation of charge-localized diabatic states

The distinction between proton-coupled electron transfer (PCET) and hydrogen atom transfer (HAT) mechanisms is important for the characterization of many chemical and biological processes. PCET and HAT mechanisms can be differentiated in terms of electronically nonadiabatic and adiabatic proton transfer, respectively. In this paper, quantitative diagnostics to evaluate the degree of electron-proton nonadiabaticity are presented. Moreover, the connection between the degree of electron-proton nonadiabaticity and the physical characteristics distinguishing PCET from HAT, namely, the extent of electronic charge redistribution, is clarified. In addition, a rigorous diabatization scheme for transforming the adiabatic electronic states into charge-localized diabatic states for PCET reactions is presented. These diabatic states are constructed to ensure that the first-order nonadiabatic couplings with respect to the one-dimensional transferring hydrogen coordinate vanish exactly. Application of these approaches to the phenoxyl-phenol and benzyl-toluene systems characterizes the former as PCET and the latter as HAT. The diabatic states generated for the phenoxyl-phenol system possess physically meaningful, localized electronic charge distributions that are relatively invariant along the hydrogen coordinate. These diabatic electronic states can be combined with the associated proton vibrational states to generate the reactant and product electron-proton vibronic states that form the basis of nonadiabatic PCET theories. Furthermore, these vibronic states and the corresponding vibronic couplings may be used to calculate rate constants and kinetic isotope effects of PCET reactions.     

Extremely narrow peaks in predissociation of sodium dimer due to rovibronic coupling

In sodium dimer the 2 (3)Pi(g), 3 (3)Pi(g), and 4 (3)Sigma(g) (+) electronic states are coupled; the coupling of the two (3)Pi(g) states is due to vibrational motion while the nonadiabatic interaction between the (3)Sigma(g) (+) and the (3)Pi(g) states-in particular, the 3 (3)Pi(g) state-is mediated by rotational interaction. The resulting vibronic problem is studied in some detail. The bound vibrational states of the 3 (3)Pi(g) and 4 (3)Pi(g) (+) states lie in the dissociation continuum of the 2 (3)Pi(g) state and become resonances due to the prevailing nonadiabatic coupling. The resonances are calculated using the complex scaling method and the available ab initio adiabatic potential energy curves. It is demonstrated that the resonances associated with rotational nonadiabatic coupling are narrower by several orders of magnitude than those that emerge from the vibrational nonadiabatic coupling. The predissociation cross section is computed and compared with experiment.     

Solute-solvent intermolecular vibronic coupling as manifested by the molecular near-field effect in resonance hyper-Raman scattering

Vibronic coupling within the excited electronic manifold of the solute all-trans-β-carotene through the vibrational motions of the solvent cyclohexane is shown to manifest as the "molecular near-field effect," in which the solvent hyper-Raman bands are subject to marked intensity enhancements under the presence of all-trans-β-carotene. The resonance hyper-Raman excitation profiles of the enhanced solvent bands exhibit similar peaks to those of the solute bands in the wavenumber region of 21,700-25,000 cm(-1) (10,850-12,500 cm(-1) in the hyper-Raman exciting wavenumber), where the solute all-trans-β-carotene shows a strong absorption assigned to the 1A(g) → 1B(u) transition. This fact indicates that the solvent hyper-Raman bands gain their intensities through resonances with the electronic states of the solute. The observed excitation profiles are quantitatively analyzed and are successfully accounted for by an extended vibronic theory of resonance hyper-Raman scattering that incorporates the vibronic coupling within the excited electronic manifold of all-trans-β-carotene through the vibrational motions of cyclohexane. It is shown that the major resonance arises from the B-term (vibronic) coupling between the first excited vibrational level (v = 1) of the 1B(u) state and the ground vibrational level (v = 0) of a nearby A(g) state through ungerade vibrational modes of both the solute and the solvent molecules. The inversion symmetry of the solute all-trans-β-carotene is preserved, suggesting the weak perturbative nature of the solute-solvent interaction in the molecular near-field effect. The present study introduces a new concept, "intermolecular vibronic coupling," which may provide an experimentally accessible∕theoretically tractable model for understanding weak solute-solvent interactions in liquid.     

Allene and pentatetraene cations as models for intramolecular charge transfer: vibronic coupling Hamiltonian and conical intersections

We consider the vibronic coupling effects involving cationic states with degenerate components that can be represented as charge localized at either end of the short cumulene molecules allene and pentatetraene. Our aim is to simulate dynamically the charge transfer process when one component is artificially depopulated. We model the Jahn-Teller vibronic interaction within these states as well as their pseudo-Jahn-Teller coupling with some neighboring states. For the manifold of these states, we have calculated cross sections of the ab initio adiabatic potential energy surfaces along all nuclear degrees of freedom, including points at large distances from the equilibrium to increase the physical significance of our model. Ab initio calculations for the cationic states of allene and pentatetraene were based on the fourth-order M?ller-Plesset method and the outer valence Green's function method. In some cases we had to go beyond this method and use the more involved third-order algebraic diagrammatic construction method to include intersections with satellite states. The parameters for a five-state, all-mode diabatic vibronic coupling model Hamiltonian were least-square fitted to these potentials. The coupling parameters for the diabatic model Hamiltonian are such that, in comparison to allene, an enhanced preference for indirect charge transfer is predicted for pentatetraene.     

Analysis of pseudo jahn–teller distortion based on natural bond orbital theory: Case study for silicene

Ground state (GS) instability of nondegenerate molecules in high symmetric structures is understood through Pseudo Jahn–Teller mixing of the electronic states through the vibronic coupling. The general approach involves setting up of a Pseudo Jahn–Teller (PJT) problem wherein one or more symmetry allowed excited states couple to the GS to create vibrational instability along a normal mode. This faces two major complications namely (1) estimating the adiabatic potential energy surfaces for the excited states which are often difficult to describe in case the excited states have charge-transfer or multi-excitonic (ME) character and (2) finding out how many such excited states (all satisfying the symmetry requirements for vibronic coupling) of increasing energies need to be coupled with the GS for a particular PJT problem. An analogous alternative approach presented here for the well-known case of symmetry breaking of planar (D_6h) hexasilabenzene (Si₆H₆) to the buckled (D_3d) structure involves identifying the second-order donor–acceptor, hyperconjugative interactions (E²_{i → j}) that stabilize the distorted structure. Following the recent work of Nori-Shargh and Weinhold, one observes that the orbitals involved in the vibronic coupling between the S₀/S_n states and those for the donor (filled)–acceptor (empty) interactions are identical. In fact, deletion of any particular pair of E²_{i → j} interaction creates vibrational instability in the buckled structure and as a corollary, deleting it for the planar structure removes its instability. The one-to-one correlation between the natural bond orbital theory and PJT theory assists in an intuitive identification of the relevant (few) excited states from a manifold of computed ones that cause symmetry breaking by vibronic coupling. © 2019 Wiley Periodicals, Inc.     

Quantifying the breakdown of the Born-Oppenheimer approximation in surface chemistry

The Born-Oppenheimer Approximation (BOA) forms the basis for calculating electronically adiabatic potential energy surfaces, thus providing the framework for developing a molecular level understanding of a variety of important chemical problems. For surface chemistry at metal surfaces, it is now clear that for some processes electronically nonadiabatic effects can be important, even dominant; however, the magnitude of BOA breakdown may vary widely from one chemical system to another. In this paper we show that molecular-beam surface scattering experiments can be used to derive quantitative information about the magnitude of BOA breakdown. A state-to-state rate model is used to interpret the pre-exponential factor of the well-known Arrhenius surface temperature dependence of the electronically nonadiabatic vibrational excitation. We also show that reference to a "thermal limit" provides a quick and simple rule of thumb for quantifying BOA breakdown. We demonstrate this approach by comparing electronically nonadiabatic vibrational inelasticity for NO(ν = 0 → 1) to NO(ν = 15 →ν'? 15) and show that the electronically nonadiabatic coupling strengths are of a similar magnitude. We compare experiments for NO and HCl scattering from Au(111) and derive the quantitative relative magnitude for the electronically nonadiabatic influences in each system. The electronically nonadiabatic influences are 300-400 times larger for NO than for HCl, for incidence energies near 0.9 eV.     

A diabatic representation including both valence nonadiabatic interactions and spin-orbit effects for reaction dynamics

A diabatic representation is convenient in the study of electronically nonadiabatic chemical reactions because the diabatic energies and couplings are smooth functions of the nuclear coordinates and the couplings are scalar quantities. A method called the fourfold way was devised in our group to generate diabatic representations for spin-free electronic states. One drawback of diabatic states computed from the spin-free Hamiltonian, called a valence diabatic representation, for systems in which spin-orbit coupling cannot be ignored is that the couplings between the states are not zero in asymptotic regions, leading to difficulties in the calculation of reaction probabilities and other properties by semiclassical dynamics methods. Here we report an extension of the fourfold way to construct diabatic representations suitable for spin-coupled systems. In this article we formulate the method for the case of even-electron systems that yield pairs of fragments with doublet spin multiplicity. For this type of system, we introduce the further simplification of calculating the triplet diabatic energies in terms of the singlet diabatic energies via Slater's rules and assuming constant ratios of Coulomb to exchange integrals. Furthermore, the valence diabatic couplings in the triplet manifold are taken equal to the singlet ones. An important feature of the method is the introduction of scaling functions, as they allow one to deal with multibond reactions without having to include high-energy diabatic states. The global transformation matrix to the new diabatic representation, called the spin-valence diabatic representation, is constructed as the product of channel-specific transformation matrices, each one taken as the product of an asymptotic transformation matrix and a scaling function that depends on ratios of the spin-orbit splitting and the valence splittings. Thus the underlying basis functions are recoupled into suitable diabatic basis functions in a manner that provides a multibond generalization of the switch between Hund's cases in diatomic spectroscopy. The spin-orbit matrix elements in this representation are taken equal to their atomic values times a scaling function that depends on the internuclear distances. The spin-valence diabatic potential energy matrix is suitable for semiclassical dynamics simulations. Diagonalization of this matrix produces the spin-coupled adiabatic energies. For the sake of illustration, diabatic potential energy matrices are constructed along bond-fission coordinates for the HBr and the BrCH(2)Cl molecules. Comparison of the spin-coupled adiabatic energies obtained from the spin-valence diabatics with those obtained by ab initio calculations with geometry-dependent spin-orbit matrix elements shows that the new method is sufficiently accurate for practical purposes. The method formulated here should be most useful for systems with a large number of atoms, especially heavy atoms, and/or a large number of spin-coupled electronic states.     

Ab initio study of the VUV-induced multistate photodynamics of formaldehyde

Although formaldehyde, H?CO, has been extensively studied there are still several issues not-well understood, specially regarding its dynamics in the VUV energy range, mainly due to the amount of nonadiabatic effects governing its dynamics. Most of the theoretical work on this molecule has focused on vertical excitation energies of Rydberg and valence states. In contrast to photodissociation processes involving the lowest-lying electronic states below 4.0 eV, there is little known about the photodynamics of the high-lying electronic states of formaldehyde (7-10 eV). One question of particular interest is why the (π, π*) electronic state is invisible experimentally even though it corresponds to a strongly dipole-allowed transition. In this work we present a coupled multisurface 2D photodynamics study of formaldehyde along the CO stretching and the symmetric HCH bending motion, using a quantum time-dependent approach. Potential energy curves along all the vibrational normal modes of formaldehyde have been computed using equation-of-motion coupled cluster including single and double excitations with a quadruply augmented basis set. In the case of the CO stretching coordinate, state-averaged complete active space self-consistent field followed by multireference configuration interaction was used for large values of this coordinate. 2D (for the CO stretching coordinate and the HCH angle) and 3D (including the out-of-plane distortion) potential energy surfaces have been computed for several Rydberg and valence states. Several conical intersections (crossings between potential energy surfaces of the same multiplicity) have been characterized and analyzed and a 2D 5 × 5 diabatic model Hamiltonian has been constructed. Based on this Hamiltonian, electronic absorption spectra, adiabatic and diabatic electronic populations and vibrational densities have been obtained and analyzed. The experimental VUV absorption spectrum in the 7-10 eV energy range is well reproduced, including the vibrational structure and the high irregularity in the regime of strong interaction between the (π, π*) electronic state and neighboring Rydberg states.     

On the excess-energy dependence of radiationless decay rate constants

The results of calculations of the dependence of the radiationless rate constant on the excess of excitation energy within the two-electronic states model under the weak coupling and statistical limits are presented. It is assumed that the exact molecular states for a given electronic configuration are global in character containing equal contributions from all degenerated vibrational levels at a given excitation energy due to intramolecular vibrational relaxation (IVR). The results of calculations indicate an important role of the low-frequency vibrational modes, the potential energy surfaces of which cross between the two electronic states involved into the radiationless process. The sharp increase of the rate constant is predicted for the excitation energy below the diabatic crossing point, followed by saturation at higher energies. The calculated rate constants for the T₁→S₀ intersystem crossing in pyrazine and benzene are in good agreement with experimental observations. Some comments concerning the “channel-three” phenomenon in benzene are presented.     

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.     

Quantum fluctuation of electronic wave-packet dynamics coupled with classical nuclear motions

An ab initio electronic wave-packet dynamics coupled with the simultaneous classical dynamics of nuclear motions in a molecule is studied. We first survey the dynamical equations of motion for the individual components. Reflecting the nonadiabatic dynamics that electrons can respond to nuclear motions only with a finite speed, the equations of motion for nuclei include a force arising from the kinematic (nuclear momentum) coupling from electron cloud. To materialize these quantum effects in the actual ab initio calculations, we study practical implementation of relevant electronic matrix elements that are related to the derivatives with respect to the nuclear coordinates. Applications of the present scheme are performed in terms of the configuration state functions (CSF) using the canonical molecular orbitals as basis functions without transformation to particular diabatic basis. In the CSF representation, the nonadiabatic interaction due to the kinematic coupling is anticipated to be rather small, and instead it should be well taken into account through the off-diagonal elements of the electronic Hamiltonian matrix. Therefore it is expected that the nonadiabatic dynamics based on this CSF basis neglecting the kinematic coupling may work. To verify this anticipation and to quantify the actual effects of the kinematic coupling, we compare the dynamics with and without the kinematic-coupling terms using the same CSF set. Applications up to the fifth electronically excited states in a nonadiabatic collision between H(2) and B(+) shows that the overall behaviors of these two calculations are surprisingly similar to each other in an average sense except for a fast fluctuation reflecting the electronic time scale. However, at the same time, qualitative differences in the collision events are sometimes observed. Therefore it turns out after all that the kinematic-coupling terms cannot be neglected in the CSF-basis representation. The present applications also demonstrate that the nonadiabatic electronic wave-packet dynamics within ab initio quantum chemical calculation is feasible.     

Intramolecular vibrational redistribution following optical excitation

We study the time-dependence of vibrational redistribution following the excitation of a zero-order harmonic state. Two models are considered: (i) the (unphysical) case of a unique anharmonic coupling matrix element gives rise to very strong coherence effects, resulting in a very slow redistribution, at the time scale of the recurrence time 2 π?, characteristic of the manifold; (ii) on the contrary, a random distribution of the coupling matrix elements results into a redistribution rate; Δ_a = 2 π? ΔV²_a — where ΔV²_a is the variance of the distribution — with no influence of its average value. This last model is extended to the case of an initial excitation of mixed states which result from a selective coupling of a radiative state with some vibronic states of a non-radiative manifold. The experimental results available in the literature are reviewed; redistribution following an optical excitation of large molecules, with a vibrational excess energy as high as 5000 to 10000 cm^?1, is generally a slow process (≤ 10⁷ s^?1), in contrast with vibrational relaxation in the condensed phases.     

The potential energy surface of the jahn-teller-distorted 2E' ground state of copper trimer

It is shown that the configuration interaction calculations of Walch and Laskowski for the ground ²E' state of Cu₃ are consistent with truncating the potential energy at second order in the distance from the D_3h conical intersection, which leads to a diabatic representation (the “quadratic coupling model”) containing three parameters. Fitting these parameters to three of the vibronic spacings of Rohlfing and Valentini leads to a reasonable explanation of their spectrum. The resulting potential energy surfaces show a Jahn-Teller stabilization energy of 221 cm⁻¹ and a pseudorotation barrier of 95 cm⁻¹.     

X-ray absorption and resonant Auger spectroscopy of O2 in the vicinity of the O 1s-->sigma* resonance: experiment and theory

We report on an experimental and theoretical investigation of x-ray absorption and resonant Auger electron spectra of gas phase O(2) recorded in the vicinity of the O 1s-->sigma(*) excitation region. Our investigation shows that core excitation takes place in a region with multiple crossings of potential energy curves of the excited states. We find a complete breakdown of the diabatic picture for this part of the x-ray absorption spectrum, which allows us to assign an hitherto unexplained fine structure in this spectral region. The experimental Auger data reveal an extended vibrational progression, for the outermost singly ionized X (2)Pi(g) final state, which exhibits strong changes in spectral shape within a short range of photon energy detuning (0 eV>Omega>-0.7 eV). To explain the experimental resonant Auger electron spectra, we use a mixed adiabatic/diabatic picture selecting crossing points according to the strength of the electronic coupling. Reasonable agreement is found between experiment and theory even though the nonadiabatic couplings are neglected. The resonant Auger electron scattering, which is essentially due to decay from dissociative core-excited states, is accompanied by strong lifetime-vibrational and intermediate electronic state interferences as well as an interference with the direct photoionization channel. The overall agreement between the experimental Auger spectra and the calculated spectra supports the mixed diabatic/adiabatic picture.     

Stark shifts of the visible absorption systems of the NO−2 ion in NaNO2

Stark effect measurements on the lowest triplet and singlet transitions of the nitrite ion show that the change of dipole moment on excitation is very small (ca. 0.30 D) and of opposite sign for the two states. The dipole moment change in the triplet manifold is also found to distinctly depend on vibrational (ν₂) excitation. This is explained as a result of vibronic coupling among B₁ triplet states.     

The Optical Spectrum of Au2+

The electronic structure of the Au₂⁺ cation is essential for understanding its catalytic activity. We present the optical spectrum of mass-selected Au₂⁺ measured via photodissociation spectroscopy. Two vibrationally resolved band systems are observed in the 290–450 nm range (at ca. 440 and ca. 325 nm), which both exhibit rather irregular structure indicative of strong vibronic and spin-orbit coupling. The experimental spectra are compared to high-level quantum-chemical calculations at the CASSCF-MRCI level including spin-orbit coupling. The results demonstrate that the understanding of the electronic structure of this simple, seemingly H₂⁺-like diatomic molecular ion strictly requires multireference and relativistic treatment including spin-orbit effects. The calculations reveal that multiple electronic states contribute to each respective band system. It is shown that popular DFT methods completely fail to describe the complex vibronic pattern of this fundamental diatomic cation.