Identification of unavoided crossings in nonadiabatic photoexcited dynamics involving multiple electronic states in polyatomic conjugated molecules

Radiationless transitions between electronic excited states in polyatomic molecules take place through unavoided crossings of the potential energy surfaces with substantial non-adiabatic coupling between the respective adiabatic states. While the extent in time of these couplings are large enough, these transitions can be reasonably well simulated through quantum transitions using trajectory surface hopping-like methods. In addition, complex molecular systems may have multiple "trivial" unavoided crossings between noninteracting states. In these cases, the non-adiabatic couplings are described as sharp peaks strongly localized in time. Therefore, their modeling is commonly subjected to the identification of regions close to the particular instantaneous nuclear configurations for which the energy surfaces actually cross each other. Here, we present a novel procedure to identify and treat these regions of unavoided crossings between non-interacting states using the so-called Min-Cost algorithm. The method differentiates between unavoided crossings between interacting states (simulated by quantum hops), and trivial unavoided crossings between non-interacting states (detected by tracking the states in time with Min-Cost procedure). We discuss its implementation within our recently developed non-adiabatic excited state molecular dynamics framework. Fragments of two- and four-ring linear polyphenylene ethynylene chromophore units at various separations have been used as a representative molecular system to test the algorithm. Our results enable us to distinguish and analyze the main features of these different types of radiationless transitions the molecular system undertakes during internal conversion.     

Non-adiabatic molecular dynamics with quantum solvent effects

Three novel approaches extending quantum-classical non-adiabatic (NA) molecular dynamics (MD) to include quantum effects of solvent environments are described. In a standard NA-MD the solute subsystem is treated quantum mechanically, while the larger solvent part of a system is treated classically. The three novel approaches presented here are based on the Bohmian formulation of quantum mechanics, the stochastic Schrödinger equation for the evolution of open quantum systems and the quantized Hamilton dynamics generalization of classical mechanics. The approaches extend the standard NA-MD to incorporate the following quantum effects of the solvent. (1) Branching, i.e. the ability of solvent quantum wave packets to split and follow asymptotically diverging trajectories correlated with different quantum states of the solute. (2) Decoherence, i.e. loss of quantum interference within the solute subsystem induced by the diverging solvent trajectories. (3) Zero point energy that contributes to NA coupling and must be preserved during the energy exchange between solvent and solute degrees of freedom. The Bohmian quantum-classical mechanics, stochastic mean-field and quantized mean-field approximations incorporate the quantum solvent effects into the standard quantum-classical NA-MD in a straightforward and efficient way that can be easily applied to quantum dynamics of condensed phase chemical systems.     

Shells and supershells in metal clusters

Concept of energy shells in finite bound quantum systems and their semi-classical interpretation are reviewed. Level densities are associated with classical periodic orbits of a particle in the mean-field. Orbits responsible for the bunching of levels, or the shells, are those with short trajectory length. The inverse of the trajectory length determines the energy spacing of the shells. When two short orbits have almost equal but slightly different lengths, the corresponding quantum level density shows a global beating pattern enveloping individual shell oscillations. This pattern is called “supershell” and may be observed in metal clusters with more than a thousand valence electrons.     

Quantum effects in ab-initio calculations of rate constants for chemical reactions occuring in the condensed phase

An overview of recent advances in the development of methods designed to calculate rate constants for chemical reactions obeying mass action kinetic equations in condensed phases is presented. A general framework addressing mixed quantum-classical systems is elaborated that enables quantum features such as tunneling effects, zero-point vibrations, dynamic quantum coherence, and non-adiabatic effects to be calculated. An efficient Monte Carlo sampling method for performing ab-initio calculations of rate constants and isotope effects in chemical processes in condensed phases is outlined, and the connection of isotope effects to reaction mechanism is explored     

Perspective: Quantum or classical coherence?

Some coherence effects in chemical dynamics are described correctly by classical mechanics, while others only appear in a quantum treatment--and when these are observed experimentally it is not always immediately obvious whether their origin is classical or quantum. Semiclassical theory provides a systematic way of adding quantum coherence to classical molecular dynamics and thus provides a useful way to distinguish between classical and quantum coherence. Several examples are discussed which illustrate both cases. Particularly interesting is the situation with electronically non-adiabatic processes, where sometimes whether the coherence effects are classical or quantum depends on what specific aspects of the process are observed.     

Non-adiabatic quantum molecular dynamics: basic formalism and case study

A general formalism is presented that treats selfconsistently and simultaneously classical atomic motion and quantum electronic excitations in dynamical processes of atomic many-body systems (non-adiabatic quantum molecular dynamics). On the basis of time-dependent density functional theory, coupled highly non-linear equations of motion are derived for arbitrary basis sets for the time-dependent Kohn-Sham orbitals. Possible approximations to make the approach practical for large atomic cluster systems are discussed. As a first application of the still exact equations of motion, non-adiabatic effects in the scattering of H⁺+H, as a case study, are investigated.     

Reduced classical trajectory equations for electronically non-adiabatic transition on photochemical pericyclic reactions

It may be difficult to use the classical trajectory equations (CTE) for the estimation of electronically non-adiabatic transition probabilities in photochemical pericyclic reactions because of many nuclear degrees of freedom. In order to avoid this difficulty, the CTE were reformulated in terms of the reaction path coordinates, and the reduced CTE were derived, in which the system was restricted to move one-dimensionally along the postulated reaction path. As an application, the non-adiabatic decay from the lowest excited state to the ground state was investigated for the conrotatory and disrotatory processes of the photochemical electrocyclic reaction of 1,3-cis-butadiene to form cyclobutene.     

Communication: Partial linearized density matrix dynamics for dissipative, non-adiabatic quantum evolution

An approach for treating dissipative, non-adiabatic quantum dynamics in general model systems at finite temperature based on linearizing the density matrix evolution in the forward-backward path difference for the environment degrees of freedom is presented. We demonstrate that the approach can capture both short time coherent quantum dynamics and long time thermal equilibration in an application to excitation energy transfer in a model photosynthetic light harvesting complex. Results are also presented for some nonadiabatic scattering models which indicate that, even though the method is based on a "mean trajectory" like scheme, it can accurately capture electronic population branching through multiple avoided crossing regions and that the approach offers a robust and reliable way to treat quantum dynamical phenomena in a wide range of condensed phase applications.     

Time-dependent Fourier grid Hamiltonian method for modelling real-time quantum dynamics: Theoretical models and applications

A local grid method for modelling real-time quantum dynamical events is formulated. The formulation is straightforward for 1-D systems. For more than one dimension, appeal has to be made to mean-field approximation of the appropriate kind. The simplest mean-field model results in time-dependent Hartree-Fourier grid method. The relationship of the proposed method with some other methods available in the literature is examined. A few examples of numerical applications dealing with (i) the dynamics of dissociation and ionization processes in diatoms and atoms respectively and (ii) tunnelling dynamics in the intramolecular H-atom transfer phenomenon are presented.     

面向有机光电材料的电子激发态结构与过程的计算方法

共轭聚合物与有机分子材料中的电子激发结构与过程决定了材料的光电功能：根据Kasha规则,低能级激发态的排序决定能否发光;最低激发态至基态的辐射跃迁与无辐射跃迁之间的竞争决定了发光效率,后者主要由非绝热耦合（声子作用）决定;电荷激发态载体的传输由电子分布与振动耦合或杂质和无序的散射弛豫过程决定．本文针对有机功能材料的发光性能,介绍两种理论方法的研究进展,即可用于计算共轭聚合物激发态结构的量子化学密度矩阵重整化群方法和计算发光效率的多模耦合无辐射跃迁速率方法．这些方法被应用于有机功能材料的性能预测和分子设计中．     

Dynamical mean-field theory from a quantum chemical perspective

We investigate the dynamical mean-field theory (DMFT) from a quantum chemical perspective. Dynamical mean-field theory offers a formalism to extend quantum chemical methods for finite systems to infinite periodic problems within a local correlation approximation. In addition, quantum chemical techniques can be used to construct new ab initio Hamiltonians and impurity solvers for DMFT. Here, we explore some ways in which these things may be achieved. First, we present an informal overview of dynamical mean-field theory to connect to quantum chemical language. Next, we describe an implementation of dynamical mean-field theory where we start from an ab initio Hartree-Fock Hamiltonian that avoids double counting issues present in many applications of DMFT. We then explore the use of the configuration interaction hierarchy in DMFT as an approximate solver for the impurity problem. We also investigate some numerical issues of convergence within DMFT. Our studies are carried out in the context of the cubic hydrogen model, a simple but challenging test for correlation methods. Finally, we finish with some conclusions for future directions.     

Implementation of ab initio multiple spawning in the Molpro quantum chemistry package

The ab initio multiple spawning (AIMS) method has been developed to solve the electronic and nuclear Schrodinger equations simultaneously for application to photochemical reaction dynamics. We discuss some details of the implementation of AIMS in the Molpro program package. A few aspects of the implementation are highlighted, including a new multiple timescale integrator and a scheme for solving the coupled-perturbed multiconfiguration self-consistent field (CP-MCSCF) equations in the context of ab initio molecular dynamics. The implementation is very efficient and we demonstrate calculations on the photoisomerization of ethylene using more than 5000 trajectory basis functions. We have included the capability for hybrid quantum mechanics/molecular mechanics (QM/MM) simulations within AIMS, and we investigate the role of an argon solvent in the photoisomerization of ethylene. Somewhat surprisingly, the surrounding argon has little effect on the timescale of non-adiabatic quenching in ethylene.     

Towards disentangling coupled electronic-vibrational dynamics in ultrafast non-adiabatic processes

Femtosecond time-resolved photoelectron spectroscopy is emerging as a new technique for investigating polyatomic excited state dynamics. Due to the sensitivity of photoelectron spectroscopy to both electronic configurations and vibrational dynamics, it is well suited to the study of non-adiabatic processes such as internal conversion, which often occur on sub-picosecond time scales. We discuss the technical requirements for such experiments, including lasers systems, energy- and angle-resolved photoelectron spectrometers and new detectors for coincidence experiments. We present a few examples of these methods applied to problems in diatomic wavepacket dynamics and ultrafast non-adiabatic processes in polyatomic molecules.     

Model calculations of harpooning elementary processes in the system K + Br2

The dynamics of non-adiabatic elementary processes in the system K + Br₂ (exchange reaction, electron transfer, dissociation) are investigated using the quasiclassical trajectory method. Based on a estimation of the transition probability from the covalent to the ionic potential surface, the treatment is restricted to a calculation of the motion (1D, 2D, and 3D) on the ionic potential which is represented by a model function. Differential and total cross sections are determined, briefly discussed and compared with molecular-beam data.     

Atoms in molecules: beyond Born–Oppenheimer paradigm

This contribution presents the first atoms in molecules study that goes beyond the Born–Oppenheimer paradigm employing the newly developed two-component quantum theory of atoms in molecules (TC-QTAIM). The LiH, LiD, and LiT systems containing quantum instead of clamped hydrogen nuclei are used as typical examples. The computational analysis that is done on non-adiabatic wavefunctions derived from the fully variational multicomponent molecular orbital approach (FV-MC-MO) results in hydrogen atomic basins without any clamped nucleus. The topological analysis of the Γ-field, the field that replaces the usual one-electron density used in the orthodox topological analysis, reveals delicate differences among the considered systems. The calculation of basin properties also demonstrates that the TC-QTAIM differentiates among atomic basins containing isotopes. Since the nuclear dynamics is contained intrinsically in non-adiabatic wavefunctions, the nuclear contribution to both topological analysis and basin properties naturally emerges from the TC-QTAIM analysis resolving the long-standing obstacle of consistent incorporation of nuclear dynamics within the context of the orthodox QTAIM. Also, a similar analysis is done on non-adiabatic wavefunctions describing excited instead of ground nuclear vibrations of the considered systems demonstrating the fact that TC-QTAIM is capable of being employed for both ground and excited nuclear vibrational states.     

Adiabatic and non-adiabatic quantum dynamics calculation of O(1D) + D2 → OD + D reaction

Adiabatic (1A' or 1A' state) and non-adiabatic (2A'/1A' states) quantum dynamics calculations have been carried out for the title reaction (O((1)D) + D(2) → OD + D) to obtain the initial state-specified (v(i) = 0, j(i) = 0) integral cross section and rate constant using the potential energy surfaces of Dobbyn and Knowles. A total of 50 partial wave contributions have been calculated using the Chebyshev wave packet method with full Coriolis coupling to achieve convergence up to the collision energy of 0.28 eV. The total integral cross section and rate constant are in excellent agreement with experimental as well as quasi-classical trajectory results. Contributions from the adiabatic pathway of the 1A' state and the non-adiabatic pathway of the 2A'/1A' states, increase significantly with the collision energy. Compared to the O((1)D) + H(2) system, the kinetic isotope effect (k(D)/k(H)) is found to be nearly temperature independent above 100 K and its value of 0.77 ± 0.01 shows excellent agreement with the experimental result of 0.81.     

Generalization of classical mechanics for nuclear motions on nonadiabatically coupled potential energy surfaces in chemical reactions

Classical trajectory study of nuclear motion on the Born-Oppenheimer potential energy surfaces is now one of the standard methods of chemical dynamics. In particular, this approach is inevitable in the studies of large molecular systems. However, as soon as more than a single potential energy surface is involved due to nonadiabatic coupling, such a naive application of classical mechanics loses its theoretical foundation. This is a classic and fundamental issue in the foundation of chemistry. To cope with this problem, we propose a generalization of classical mechanics that provides a path even in cases where multiple potential energy surfaces are involved in a single event and the Born-Oppenheimer approximation breaks down. This generalization is made by diagonalization of the matrix representation of nuclear forces in nonadiabatic dynamics, which is derived from a mixed quantum-classical representation of the electron-nucleus entangled Hamiltonian [Takatsuka, K. J. Chem. Phys. 2006, 124, 064111]. A manifestation of quantum fluctuation on a classical subsystem that directly contacts with a quantum subsystem is discussed. We also show that the Hamiltonian thus represented gives a theoretical foundation to examine the validity of the so-called semiclassical Ehrenfest theory (or mean-field theory) for electron quantum wavepacket dynamics, and indeed, it is pointed out that the electronic Hamiltonian to be used in this theory should be slightly modified.     

Massively parallel spin–orbit configuration interaction

We describe a new approach to incorporating quantum effects into chemical reaction rate theory using quantum trajectories. Our development is based on the entangled trajectory molecular dynamics method for simulating quantum processes using trajectory integration and ensemble averaging. By making dynamical approximations similar to those underlying classical transition state theory, quantum corrections are incorporated analytically into the quantum rate expression. We focus on a simple model of quantum decay in a metastable system and consider the deep tunneling limit where the classical rate vanishes and the process is entirely quantum mechanical. We compare our approximate estimate with the well-known WKB tunneling rate and find qualitative agreement.     

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

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