Adiabatic and Quasidiabatic States in a Gauge Theoretical Framework

Compensation States Approach in the Hybrid Diabatization Scheme: Extension to Multidimensional Data and Properties

The diabatization of reactive systems for more than just a couple of states is a very demanding problem and generally requires advanced diabatization techniques. Especially for dissociative processes, the drastic changes in the adiabatic wave functions often would require large diabatic state bases, which quickly become impractical. Recently, we addressed this problem by the compensation states approach developed in the context of our hybrid diabatization scheme. This scheme utilizes wave function as well as energy data in combination with a diabatic potential model. In regions where the initial diabatic state basis becomes insufficient for an appropriate representation of the adiabatic states, new model states are generated. The new model states compensate for the state space not spanned by the initial diabatic basis. Such a compensation state is obtained by projecting the initial diabatic state space out of the adiabatic wave function. This yields a very efficient basis representation of the electronic Hamiltonian. The present work presents two new aspects. First, it is shown how other operators like the spin-orbit operator in the framework of the Effective Relativistic Coupling by Asymptotic Representation (ERCAR) can be evaluated in this compact model state space without losing the correct wave function information and accuracy. Second, the extension of the approach to multidimensional potential energy surface models is presented for methyl iodide including the C-I dissociation coordinate and the angular H3C-I bending coordinates.

Vibronic coupling and ultrafast relaxation dynamics in the first five excited singlet electronic states of bithiophene

The vibronic structure and nuclear dynamics in the first five excited singlet electronic states of bithiophene (2T) are investigated here. Specifically, considerations are given to comprehend the first two structureless and broad electronic absorption bands and the role of nonadiabatic coupling in the excited state relaxation mechanism of 2T in the gas phase. Associated potential energy surfaces (PESs) are established by constructing a model vibronic coupling Hamiltonian using 18 vibrational degrees of freedom and extensive ab initio electronic structure calculations. The topographies of these PESs are critically examined, and multiple conical intersections are established. The nuclear dynamics calculations are performed by propagating wave packets on the coupled electronic manifold. The present theoretical results are in good agreement with the experimental observations. It is found that strong nonadiabatic coupling between the S1-S4 and S1-S5 states along totally symmetric modes is predominantly responsible for the structureless and broad first absorption band, and overlapping S2, S3, S4, and S5 states form the second absorption band. Photorelaxation from the highly excited S5 to the lowest S1 state takes place through a cascade of diabatic population transfers among the S1-S4-S5 electronic manifold within the first ∼100 fs. Totally symmetric C=C stretching, C-S stretching, C-H wagging, ring puckering, and inter-ring bending modes collectively drive such relaxation dynamics.

Time-Dependent Resonant Inelastic X-ray Scattering of Pyrazine at the Nitrogen K-Edge: A Quantum Dynamics Approach

We calculate resonant inelastic X-ray scattering spectra of pyrazine at the nitrogen K-edge in the time domain including wavepacket dynamics in both the valence and core-excited state manifolds. Upon resonant excitation, we observe ultrafast non-adiabatic population transfer between core-excited states within the core-hole lifetime, leading to molecular symmetry distortions. Importantly, our time-domain approach inherently contains the ability to manipulate the dynamics of this process by detuning the excitation energy, which effectively shortens the scattering duration. We also explore the impact of pulsed incident X-ray radiation, which provides a foundation for state-of-the-art time-resolved experiments with coherent pulsed light sources.

The effective relativistic coupling by asymptotic representation approach for molecules with multiple relativistic atoms

The Effective Relativistic Coupling by Asymptotic Representation (ERCAR) approach is a method to generate fully coupled diabatic potential energy surfaces (PESs) including relativistic effects, especially spin-orbit coupling. The spin-orbit coupling of a full molecule is determined only by the atomic states of selected relativistically treated atoms. The full molecular coupling effect is obtained by a diabatization with respect to asymptotic states, resulting in the correct geometry dependence of the spin-orbit effect. The ERCAR approach has been developed over the last decade and initially only for molecules with a single relativistic atom. This work presents its extension to molecules with more than a single relativistic atom using the iodine molecule as a proof-of-principle example. The theory for the general multiple atomic ERCAR approach is given. In this case, the diabatic basis is defined at the asymptote where all relativistic atoms are separated from the remaining molecular fragment. The effective spin-orbit operator is then a sum of spin-orbit operators acting on isolated relativistic atoms. PESs for the iodine molecule are developed within the new approach and it is shown that the resulting fine structure states are in good agreement with spin-orbit ab initio calculations.

A Local Diabatisation Method for Two-State Adiabatic Conical Intersections

A methodology to locally characterize conical intersections (CIs) between two adiabatic electronic states for which no nonadiabatic coupling (NAC) vectors are available is presented. Based on the Hessian and gradient at the CI, the branching space coordinates are identified. The potential energy surface around the CI in the branching space is expressed in the diabatic representation, from which the NAC vectors can be calculated in a wave-function-free, energy-based approach. To demonstrate the universality of the developed methodology, the minimum-energy CI (MECI) between the first (S1) and second (S2) singlet excited states of formamide is investigated at the state-averaged complete active space self-consistent field (SA-CASSCF) and extended multistate complete active space second-order perturbation theory (XMS-CASPT2) levels of theory. In addition, the asymmetrical MECI between the ground state (S0) and S1 of cyclopropanone is evaluated using SA-CASSCF, as well as (ME)CIs between the S1 and S2 states of benzene using SA-CASSCF and time-dependent density functional theory (TDDFT). Finally, a CI between the S1 and S2 excited states of thiophene was analyzed using TDDFT.

Nonadiabatic Field on Quantum Phase Space: A Century after Ehrenfest

Nonadiabatic transition dynamics lies at the core of many electron/hole transfer, photoactivated, and vacuum field-coupled processes. About a century after Ehrenfest proposed "Phasenraum" and the Ehrenfest theorem, we report a conceptually novel trajectory-based nonadiabatic dynamics approach, nonadiabatic field (NAF), based on a generalized exact coordinate-momentum phase space formulation of quantum mechanics. It does not employ the conventional Born-Oppenheimer or Ehrenfest trajectory in the nonadiabatic coupling region. Instead, in NAF the equations of motion of the independent trajectory involve a nonadiabatic nuclear force term in addition to an adiabatic nuclear force term of a single electronic state. A few benchmark tests for gas phase and condensed phase systems indicate that NAF offers a practical tool to capture the correct correlation of electronic and nuclear dynamics for processes where the states remain coupled all the time as well as for the asymptotic region where the coupling of electronic states vanishes.

Quantum Dynamics of Photodissociation: Recent Advances and Challenges

Recent advances in constructing accurate potential energy surfaces and nonadiabatic couplings from high-level ab initio data have revealed detailed potential landscapes in not only the ground electronic state but also excited ones. They enabled quantitatively accurate characterization of photoexcited reactive systems using quantum mechanical methods. In this Perspective, we survey the recent progress in quantum mechanical studies of adiabatic and nonadiabatic photodissociation dynamics, focusing on initial state control and product energy disposal. These new insights helped to understand quantum effects in small prototypical systems, and the results serve as benchmarks for developing more approximate theoretical methods.

Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals

Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.

Cooperative molecular structure in polaritonic and dark states

The interaction of quantum light with matter is known to give rise to mixed light-matter states. An ensemble of identical molecules is discussed. The resulting hybrid light-matter states exhibit complex structure even if only a single vibrational coordinate per molecule is considered. Starting from the uniform situation where all molecules possess the same value of this coordinate, polaritons and dark states follow like in atoms, but are functions of this coordinate. It is proven that any point on a resulting polariton energy curve is a (local) minimum or maximum for distorting molecules perpendicular to this curve. It is shown how to explicitly compute the impact of distortion solely based on the data of a free molecule. The structure of the dark states and their behavior upon distortion is analyzed as well. Useful techniques are introduced and general results on, for example, minimum energy path, symmetry breaking and restoration, are obtained. The developed strategy is transferred to include several or even many nuclear degrees of freedom per molecule and it is demonstrated that the interplay of several vibrational degrees of freedom in a single molecule of the ensemble is expected to lead to qualitatively different physics. General consequences are discussed.

Importance of dynamical electron correlation in diabatic couplings of electron-exchange processes

We demonstrate the importance of the dynamical electron correlation effect in diabatic couplings of electron-exchange processes in molecular aggregates. To perform a multireference perturbation theory with large active space of molecular aggregates, an efficient low-rank approximation is applied to the complete active space self-consistent field reference functions. It is known that kinetic rates of electron-exchange processes, such as singlet fission, triplet-triplet annihilation, and triplet exciton transfer, are not sufficiently explained by the direct term of the diabatic couplings but efficiently mediated by the low-lying charge transfer states if the two molecules are in close proximity. It is presented in this paper, however, that regardless of the distance of the molecules, the direct term is considerably underestimated by up to three orders of magnitude without the dynamical electron correlation, i.e., the diabatic states expressed in the active space are not adequate to quantitatively reproduce the electron-exchange processes.

Incorporating Berry Force Effects into the Fewest Switches Surface-Hopping Algorithm: Intersystem Crossing and the Case of Electronic Degeneracy

We present a preliminary surface-hopping approach for modeling intersystem crossing (ISC) dynamics between four electronic states: one singlet and one (triply degenerate) triplet. In order to incorporate all Berry force effects, the algorithm requires that, when moving along an adiabatic surface associated with the triplet manifold, one must also keep track of a quasi-diabatic index (akin to a "m_s" quantum number) for each trajectory. For a simple model problem, we find that a great deal of new physics can be captured by our algorithm, setting the stage for larger, more realistic (or perhaps even ab initio) simulations in the future.

Diabatic States of Molecules

Quantitative simulations of electronically nonadiabatic molecular processes require both accurate dynamics algorithms and accurate electronic structure information. Direct semiclassical nonadiabatic dynamics is expensive due to the high cost of electronic structure calculations, and hence it is limited to small systems, limited ensemble averaging, ultrafast processes, and/or electronic structure methods that are only semiquantitatively accurate. The cost of dynamics calculations can be made manageable if analytic fits are made to the electronic structure data, and such fits are most conveniently carried out in a diabatic representation because the surfaces are smooth and the couplings between states are smooth scalar functions. Diabatic representations, unlike the adiabatic ones produced by most electronic structure methods, are not unique, and finding suitable diabatic representations often involves time-consuming nonsystematic diabatization steps. The biggest drawback of using diabatic bases is that it can require large amounts of effort to perform a globally consistent diabatization, and one of our goals has been to develop methods to do this efficiently and automatically. In this Feature Article, we introduce the mathematical framework of diabatic representations, and we discuss diabatization methods, including adiabatic-to-diabatic transformations and recent progress toward the goal of automatization.

Extended Mulliken-Hush Method with Applications to the Theoretical Study of Electron Transfer

A novel adiabatic-to-diabatic (ATD) transformation strategy, namely, the extended Mulliken-Hush (XMH) method, is proposed to evaluate diabatic properties including electronic couplings, potential energy surfaces, and their crossings. The XMH method is developed by adopting our recently proposed ATD transformation formula of a general vectorial physical observable, in which a useful ATD transformation is further determined by using an auxiliary dipole between localized frontier orbitals as a simple approximation of the diabatic transition dipole. The XMH method is simple and practical that provides a flexible way to construct diabatic states. To some extent, it can be regarded as an extension of the generalized Mulliken-Hush (GMH) method since the latter takes a stronger approximation, in which the diabatic transition dipole is assumed to be vanishing. Test calculations on the HeH₂⁺ system show that the electronic couplings predicted by the XMH method are closer to the ones calculated by the valence bond block-diagonalization approach than the GMH ones since the XMH method takes into account both the magnitude and direction of the diabatic transition dipole, which is consistent with the properties of this molecule. In the study of electron transfer in the two kinds of donor-bridge-acceptor systems, the XMH method maintains the simplicity of the GMH method and gives reasonable results even when the latter fails, wherein the diabatic transition dipole is nearly perpendicular to the difference of the initial and final adiabatic dipoles. More importantly, the XMH method can be easily combined with high-level electronic structure methods, in which the properties of the ground and excited states may be more accurately calculated, and hence, one may expect that further development of the XMH method would result in a general computational model for studying electron transfer reactions.

Why the CC Stretch in HCC Is So Anharmonic

The marked anharmonicity of the CC stretching mode in the ethynyl radical (C₂H) is investigated in terms of a vibronic coupling model. It is demonstrated that the large difference between the harmonic frequency and the fundamental level (about 5-10 times that for triple-bond stretches in the related species HCCH, HCN, HNC and CN) can be attributed to the well-known vibronic interaction between the X̃²Σ⁺ and Ã²Π states of CCH. Although the mode has σ symmetry and it is the perturbations of π symmetry that mix the two electronic states, a combination of large intrinsic coupling strength, modest energy gap, and-most importantly-the strong tuning of the gap energy by the CC stretch mode leads to a profound vibronic influence on parts of the anharmonic force field that sample the CC stretch. Finally, calculations of the force field for the X̃²Σ⁺ state with different flavors of coupled-cluster (CC) theory provide insight and underscore an intrinsic advantage of equation-of-motion CC (EOM-CC) methods.

Vacuum Ultraviolet Excited State Dynamics of the Smallest Ketone: Acetone

We combined tunable vacuum-ultraviolet time-resolved photoelectron spectroscopy (VUV-TRPES) with high-level quantum dynamics simulations to disentangle multistate Rydberg-valence dynamics in acetone. A femtosecond 8.09 eV pump pulse was tuned to the sharp origin of the A₁(n3d_yz) band. The ensuing dynamics were tracked with a femtosecond 6.18 eV probe pulse, permitting TRPES of multiple excited Rydberg and valence states. Quantum dynamics simulations reveal coherent multistate Rydberg-valence dynamics, precluding simple kinetic modeling of the TRPES spectrum. Unambiguous assignment of all involved Rydberg states was enabled via the simulation of their photoelectron spectra. The A₁(ππ*) state, although strongly participating, is likely undetectable with probe photon energies ≤8 eV and a key intermediate, the A₂(nπ*) state, is detected here for the first time. Our dynamics modeling rationalizes the temporal behavior of all photoelectron transients, allowing us to propose a mechanism for VUV-excited dynamics in acetone which confers a key role to the A₂(nπ*) state.

Electronic Structure Methods for the Description of Nonadiabatic Effects and Conical Intersections

Nonadiabatic effects are ubiquitous in photophysics and photochemistry, and therefore, many theoretical developments have been made to properly describe them. Conical intersections are central in nonadiabatic processes, as they promote efficient and ultrafast nonadiabatic transitions between electronic states. A proper theoretical description requires developments in electronic structure and specifically in methods that describe conical intersections between states and nonadiabatic coupling terms. This review focuses on the electronic structure aspects of nonadiabatic processes. We discuss the requirements of electronic structure methods to describe conical intersections and nonadiabatic couplings, how the most common excited state methods perform in describing these effects, and what the recent developments are in expanding the methodology and implementing nonadiabatic couplings.

How important are the residual nonadiabatic couplings for an accurate simulation of nonadiabatic quantum dynamics in a quasidiabatic representation?

Diabatization of the molecular Hamiltonian is a standard approach to remove the singularities of nonadiabatic couplings at conical intersections of adiabatic potential energy surfaces. In general, it is impossible to eliminate the nonadiabatic couplings entirely-the resulting "quasidiabatic" states are still coupled by smaller but nonvanishing residual nonadiabatic couplings, which are typically neglected. Here, we propose a general method for assessing the validity of this potentially drastic approximation by comparing quantum dynamics simulated either with or without the residual couplings. To make the numerical errors negligible to the errors due to neglecting the residual couplings, we use the highly accurate and general eighth-order composition of the implicit midpoint method. The usefulness of the proposed method is demonstrated on nonadiabatic simulations in the cubic Jahn-Teller model of nitrogen trioxide and in the induced Renner-Teller model of hydrogen cyanide. We find that, depending on the system, initial state, and employed quasidiabatization scheme, neglecting the residual couplings can result in wrong dynamics. In contrast, simulations with the exact quasidiabatic Hamiltonian, which contains the residual couplings, always yield accurate results.

Non-adiabatic quantum dynamics without potential energy surfaces based on second-quantized electrons: Application within the framework of the MCTDH method

A first principles quantum formalism to describe the non-adiabatic dynamics of electrons and nuclei based on a second quantization representation (SQR) of the electronic motion combined with the usual representation of the nuclear coordinates is introduced. This procedure circumvents the introduction of potential energy surfaces and non-adiabatic couplings, providing an alternative to the Born-Oppenheimer approximation. An important feature of the molecular Hamiltonian in the mixed first quantized representation for the nuclei and the SQR representation for the electrons is that all degrees of freedom, nuclear positions and electronic occupations, are distinguishable. This makes the approach compatible with various tensor decomposition Ansätze for the propagation of the nuclear-electronic wavefunction. Here, we describe the application of this formalism within the multi-configuration time-dependent Hartree framework and its multilayer generalization, corresponding to Tucker and hierarchical Tucker tensor decompositions of the wavefunction, respectively. The approach is applied to the calculation of the photodissociation cross section of the HeH⁺ molecule under extreme ultraviolet irradiation, which features non-adiabatic effects and quantum interferences between the two possible fragmentation channels, He + H⁺ and He⁺ + H. These calculations are compared with the usual description based on ab initio potential energy surfaces and non-adiabatic coupling matrix elements, which fully agree. The proof-of-principle calculations serve to illustrate the advantages and drawbacks of this formalism, which are discussed in detail, as well as possible ways to overcome them. We close with an outlook of possible application domains where the formalism might outperform the usual approach, for example, in situations that combine a strong static correlation of the electrons with non-adiabatic electronic-nuclear effects.

Compact Bases for Vibronic Coupling in Spectral Simulations: The Photoelectron Spectrum of Cyclopentoxide in the Full 39 Internal Modes

We report an algorithm to automatically generate compact multimode vibrational bases for the Köppel-Domcke-Cederbaum (KDC) vibronic coupling wave function used in spectral simulations of moderate-sized molecules. As a full quantum method, the size of the vibronic expansion grows exponentially with respect to the number of vibrational modes, necessitating compact bases for moderate-sized systems. The problem of generating such a basis consists of two parts: one is the choice of vibrational normal modes, and the other is the number of phonons allowed in each mode. A previously developed final-state-biased technique addresses the former part, and this work focuses on the latter part: proposing an algorithm for generating an optimal phonon distribution. By virtue of this phonon distribution, compact and affordable bases can be automatically generated for systems with on the order of 15 atoms. Our algorithm is applied to determine the nonadiabatic photoelectron spectrum of cyclopentoxide in the full 39 internal modes.

The generality of the GUGA MRCI approach in COLUMBUS for treating complex quantum chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

Analytic non-adiabatic derivative coupling terms for spin-orbit MRCI wavefunctions. I. Formalism

Analytic gradients of electronic eigenvalues require one calculation per nuclear geometry, compared to at least 3n + 1 calculations for finite difference methods, where n is the number of nuclei. Analytic nonadiabatic derivative coupling terms (DCTs), which are calculated in a similar fashion, are used to remove nondiagonal contributions to the kinetic energy operator, leading to more accurate nuclear dynamics calculations than those that employ the Born-Oppenheimer approximation, i.e., that assume off-diagonal contributions are zero. The current methods and underpinnings for calculating both of these quantities, gradients and DCTs, for the State-Averaged MultiReference Configuration Interaction with Singles and Doubles (MRCI-SD) wavefunctions in COLUMBUS are reviewed. Before this work, these methods were not available for wavefunctions of a relativistic MRCI-SD Hamiltonian. Calculation of these terms is critical in successfully modeling the dynamics of systems that depend on transitions between potential energy surfaces split by the spin-orbit operator, such as diode-pumped alkali lasers. A formalism for calculating the transition density matrices and analytic derivative coupling terms for such systems is presented.

A generalized vibronic-coupling Hamiltonian for molecules without symmetry: Application to the photoisomerization of benzopyran

We present a model for the lowest two potential energy surfaces (PESs) that describe the photoinduced ring-opening reaction of benzopyran taken as a model compound to study the photochromic ring-opening reaction of indolinobenzospiropyran and its evolution toward its open-chain analog. The PESs are expressed in terms of three effective rectilinear coordinates. One corresponds to the direction between the equilibrium geometry in the electronic ground state, referred to as the Franck-Condon geometry, and the minimum of conical intersection (CI), while the other two span the two-dimensional branching space at the CI. The model correctly reproduces the topography of the PESs. The ab initio calculations are performed with the extended multiconfiguration quasidegenerate perturbation theory at second order method. We demonstrate that accounting for electron dynamic correlation drastically changes the global energy landscape since some zwitterionic states become strongly stabilized. Quantum dynamics calculations using this PES model produce an absorption spectrum that matches the experimental one to a good accuracy.

Up to a Sign. The Insidious Effects of Energetically Inaccessible Conical Intersections on Unimolecular Reactions

It is now well established that conical intersections play an essential role in nonadiabatic radiationless decay where their double-cone topography causes them to act as efficient funnels channeling wave packets from the upper to the lower adiabatic state. Until recently, little attention was paid to the effect of conical intersections on dynamics on the lower state, particularly when the total energy involved is significantly below that of the conical intersection seam. This energetic deficiency is routinely used as a sufficient condition to exclude consideration of excited states in ground state dynamics. In this account, we show that, this energy criterion notwithstanding, energy inaccessible conical intersections can and do exert significant influence on lower state dynamics. The origin of this influence is the geometric phase, a signature property of conical intersections, which is the fact that the real-valued electronic wave function changes sign when transported along a loop containing a conical intersection, making the wave function double-valued. This geometric phase is permitted by an often neglected property of the real-valued adiabatic electronic wave function; namely, it is determined only up to an overall sign. Noting that in order to change sign a normalized, continuous function must go through zero, for loops of ever decreasing radii, demonstrating the need for an electronic degeneracy (intersection) to accompany the geometric phase. Since the total wave function must be single-valued a compensating geometry dependent phase needs to be included in the total electronic-nuclear wave function. This Account focuses on how this consequence of the geometric phase can modify nuclear dynamics energetically restricted to the lower state, including tunneling dynamics, in directly measurable ways, including significantly altering tunneling lifetimes, thus confounding the relation between measured lifetimes and barrier heights and widths, and/or completely changing product rotational distributions. Some progress has been made in understanding the origin of this effect. It has emerged that for a system where the lower adiabatic potential energy surface exhibits a topography comprised of two saddle points separated by a high energy conical intersection, the effect of the geometric phase can be quite significant. In this case topologically distinct paths through the two adiabatic saddle points may lead to interference. This was pointed out by Mead and Truhlar almost 50 years ago and denoted the Molecular Aharonov-Bohm effect. Still, the difficulty in anticipating a significant geometric phase effect in tunneling dynamics due to energetically inaccessible conical intersections leads to the attribute insidious that appears in the title of this Account. Since any theory is only as relevant as the prevalence of the systems it describes, we include in this Account examples of real systems where these effects can be observed. The accuracy of the reviewed calculations is high since we use fully quantum mechanical dynamics and construct the geometric phase using an accurate diabatic state fit of high quality ab initio data, energies, energy gradients, and interstate couplings. It remains for future work to establish the prevalence of this phenomenon and its deleterious effects on the conventional wisdom discussed in this work.

Charge Migration in Propiolic Acid: A Full Quantum Dynamical Study

Ionization of molecules very often populates several cationic states launching pure electron dynamics that appear as ultrafast migration of the hole charge throughout the system. A crucial question in the emerging field of attochemistry is whether these pure electronic coherences last long enough to allow for their efficient observation and eventual manipulation with ultrashort laser pulses. We report a full-dimensional quantum calculation of concerted electron-nuclear dynamics initiated by outer-valence ionization of propiolic acid molecule, showing that the charge will oscillate between the carbon triple bond and the carbonyl oxygen for more than 10 fs before getting trapped by the nuclear motion. This time is enough for the charge migration to be observed and controlled. We argue that the molecule is very suitable for experimental studies.

Ab initio adiabatic and quasidiabatic potential energy surfaces of H(+) + CO system: A study of the ground and the first three excited electronic states

The global ground and first three excited electronic state adiabatic as well as the corresponding quasidiabatic potential energy surfaces is reported as a function of nuclear geometries in the Jacobi coordinates (R→,r→,γ) using Dunning's cc-pVTZ basis set at the internally contracted multi-reference (single and double) configuration interaction level of accuracy. Nonadiabatic couplings, arising out of relative motion of proton and the vibrational motion of CO, are also reported in terms of coupling potentials. The quasidiabatic potential energy surfaces and the coupling potentials have been obtained using the ab initio procedure [Simah et al., J. Chem. Phys. 111, 4523 (1999)] for the purpose of dynamics studies.