Semiclassical Liouville method for the simulation of electronic transitions: single ensemble formulation

In this paper, we describe a single ensemble implementation of the semiclassical Liouville method for simulating quantum processes using classical trajectories. In this approach, one ensemble of trajectories supports the evolution of all semiclassical density matrix elements, rather than employing a distinct ensemble for each. The ensemble evolves classically under a single reference Hamiltonian, which is chosen based on physical grounds; for electronic relaxation of an initially excited state, the initially populated upper surface Hamiltonian is the natural choice. Classical trajectories evolving on the reference potential then represent the time-dependent upper state population density and also the electronic coherence and the ground state density created by electronic transition. The error made in the classical motion of the trajectories for these latter distributions is compensated for by incorporating the difference between the correct and reference Liouville propagators into the calculation of the coefficients of the individual trajectories. This approach gives very accurate results for a number of model problems and cases describing ultrafast electronic relaxation dynamics.     

Interference and quantization in semiclassical response functions

Application of the Herman-Kluk semiclassical propagator to the calculation of spectroscopic response functions for anharmonic oscillators has demonstrated the quantitative accuracy of these approximate dynamics. In this approach, spectroscopic response functions are expressed as multiple phase-space integrals over pairs of classical trajectories and their associated stability matrices. Here we analyze the Herman-Kluk semiclassical approximation to a linear response function and determine the origin of the capacity of this method to reproduce quantum effects in a response function from classical dynamical information. Our analysis identifies those classical trajectories that contribute most significantly to the response function on different time scales. This finding motivates a procedure for computing the linear response function in which the interference between pairs of classical trajectories is treated approximately, resulting in an integral over a single average trajectory, as in a purely classical calculation.     

Simulation of environmental effects on coherent quantum dynamics in many-body systems

In this paper we describe an application of the trajectory-based semiclassical Liouville method for modeling coherent molecular dynamics on multiple electronic surfaces to the treatment of the evolution and decay of quantum electronic coherence in many-body systems. We consider a model representing the coherent evolution of quantum wave packets on two excited electronic surfaces of a diatomic molecule in the gas phase and in rare gas solvent environments, ranging from small clusters to a cryogenic solid. For the gas phase system, the semiclassical trajectory method is shown to reproduce the evolution of the electronic-nuclear coherence nearly quantitatively. The dynamics of decoherence are then investigated for the solvated systems using the semiclassical approach. It is found that, although solvation in general leads to more rapid and extensive loss of quantum coherence, the details of the coupled system-bath dynamics are important, and in some cases the environment can preserve or even enhance quantum coherence beyond that seen in the isolated system.     

Measuring nonadiabaticity of molecular quantum dynamics with quantum fidelity and with its efficient semiclassical approximation

We propose to measure nonadiabaticity of molecular quantum dynamics rigorously with the quantum fidelity between the Born-Oppenheimer and fully nonadiabatic dynamics. It is shown that this measure of nonadiabaticity applies in situations where other criteria, such as the energy gap criterion or the extent of population transfer, fail. We further propose to estimate this quantum fidelity efficiently with a generalization of the dephasing representation to multiple surfaces. Two variants of the multiple-surface dephasing representation (MSDR) are introduced, in which the nuclei are propagated either with the fewest-switches surface hopping or with the locally mean field dynamics (LMFD). The LMFD can be interpreted as the Ehrenfest dynamics of an ensemble of nuclear trajectories, and has been used previously in the nonadiabatic semiclassical initial value representation. In addition to propagating an ensemble of classical trajectories, the MSDR requires evaluating nonadiabatic couplings and solving the Schro?dinger (or more generally, the quantum Liouville-von Neumann) equation for a single discrete degree of freedom. The MSDR can be also used in the diabatic basis to measure the importance of the diabatic couplings. The method is tested on three model problems introduced by Tully and on a two-surface model of dissociation of NaI.     

利用基于直接动力学的轨线面跳跃方法处理非绝热过程

势能面交叉引起的非绝热过程广泛存在于光化学和光物理中。对这一过程进行描述是理论化学的重要挑战之一。非绝热过程涉及原子核与电子之间的耦合运动,因此量子化学的基本假设之一"玻恩-奥本海默"近似被打破,所以对其进行描述需要发展新的动力学理论方法。在这些方法中,Tully发展的最少轨线面跳跃方法凭借易于程序化、便于计算等优点已经发展成为处理非绝热问题的主要动力学方法之一。其中原子核以经典的方式在单一势能面上进行演化,电子以量子的方式沿着同一轨线进行演化。在整个演化过程中,非绝热跃迁通过轨线在不同势能面间的跃迁来描述,其中跳跃发生的几率与电子的演化有关。如果将该方法与从头算直接动力学相结合,可以在全原子水平上研究实际分子体系的非绝热动力学,给出其激发态寿命、非绝热动力学中分子的主要运动方式、反应通道以及分支比等重要信息。本文旨在讨论最少面跳跃直接动力学方法研究非绝热问题的一些进展,包括动力学基本理论,特别关注将最少面跳跃方法和直接动力学结合的数值实现细节,同时讨论该方法在研究实际体系当中的一些应用,并对轨线面跳跃方法下一步发展的一些方向进行合理的展望。     

Nonstatistical effects in the dissociation of ethyl radical: Finding order in chaos

How does one identify order in complex dynamical systems? A Born-Oppenheimer molecular dynamics simulation of the dissociation of ethyl radical, C(2)H(5), produces an ensemble of classical trajectories which are decomposed in the time-frequency domain using wavelets. A time-dependent scalar metric, the normalized instantaneous orbital complexity, is constructed and shown to correlate not only to the more conventional Lyapunov exponents but also to the dissociation time for an individual trajectory. The analysis of the ensemble of trajectories confirms that the long-lived trajectories are associated with a low degree of ergodicity. While the analysis of molecular dissociation dynamics is the narrow focus of the present work, the method is more general for discovery and identification of ordered regimes within large sets of chaotic data.     

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.     

Algorithms for sampling a quantum microcanonical ensemble of harmonic oscillators at potential minima and conical intersections

Algorithms are presented for sampling quantum microcanonical ensembles for a potential energy minimum and for the conical intersection at the minimum energy crossing point of two coupled electronic states. These ensembles may be used to initialize trajectories for chemical dynamics simulations. The unimolecular dynamics of a microcanonical ensemble about a potential energy minimum may be compared with the dynamics predicted by quantum Rice-Ramsperger-Kassel-Marcus (RRKM) theory. If the dynamics is non-RRKM, it will be of particular interest to determine which states have particularly long lifetimes. Initializing a microcanonical ensemble for the electronically excited state at a conical intersection is a model for electronic nonadiabatic dynamics. The trajectory surface-hopping approach may be used to study the ensuing chemical dynamics. A strength of the model is that zero-point energy conditions are included for the initial nonadiabatic dynamics at the conical intersection.     

Determination of molecular vibrational state energies using the ab initio semiclassical initial value representation: application to formaldehyde

We have demonstrated the use of ab initio molecular dynamics (AIMD) trajectories to compute the vibrational energy levels of molecular systems in the context of the semiclassical initial value representation (SC-IVR). A relatively low level of electronic structure theory (HF/3-21G) was used in this proof-of-principle study. Formaldehyde was used as a test case for the determination of accurate excited vibrational states. The AIMD-SC-IVR vibrational energies have been compared to those from curvilinear and rectilinear vibrational self-consistent field/vibrational configuration interaction with perturbation selected interactions-second-order perturbation theory (VSCF/VCIPSI-PT2) and correlation-corrected vibrational self-consistent field (cc-VSCF) methods. The survival amplitudes were obtained from selecting different reference wavefunctions using only a single set of molecular dynamics trajectories. We conclude that our approach is a further step in making the SC-IVR method a practical tool for first-principles quantum dynamics simulations.     

Semiclassical calculation of the vibrational echo

The infrared echo measurement probes the time scales of the molecular motions that couple to a vibrational transition. Computation of the echo observable within rigorous quantum mechanics is problematic for systems with many degrees of freedom, motivating the development of semiclassical approximations to the nonlinear optical response. We present a semiclassical approximation to the echo observable, based on the Herman-Kluk propagator. This calculation requires averaging over a quantity generated by two pairs of classical trajectories and associated stability matrices, connected by a pair of phase-space jumps. Quantum, classical, and semiclassical echo calculations are compared for a thermal ensemble of noninteracting anharmonic oscillators. The semiclassical approach uses input from classical mechanics to reproduce the significant features of a complete, quantum mechanical calculation of the nonlinear response.     

Generation of initial trajectories for transition path sampling of chemical reactions with ab initio molecular dynamics

Transition path sampling is an innovative method for focusing a molecular dynamics simulation on a reactive event. Although transition path sampling methods can generate an ensemble of reactive trajectories, an initial reactive trajectory must be generated by some other means. In this paper, the authors have evaluated three methods for generating initial reactive trajectories for transition path sampling with ab initio molecular dynamics. The authors have tested each of these methods on a set of chemical reactions involving the breaking and making of covalent bonds: the 1,2-hydrogen elimination in the borane-ammonia adduct, a tautomerization, and the Claisen rearrangement. The first method is to initiate trajectories from the potential energy transition state, which was effective for all reactions in the test set. Assigning atomic velocities found using normal mode analysis greatly improved the success of this method. The second method uses a high temperature molecular dynamics simulation and then iteratively reduces the total energy of the simulation until a low temperature reactive trajectory is found. This was effective in generating a low temperature trajectory from an initial trajectory run at 3000 K of the tautomerization reaction, although it failed for the other two. The third uses an orbital based bias potential to find a reactive trajectory and uses this trajectory to initiate an unbiased trajectory. The authors found that a highest occupied molecular orbital-lowest unoccupied molecular orbital bias could be used to find a reactive trajectory for the Claisen rearrangement, although it failed for the other two reactions. These techniques will help make it practical to use transition path sampling to study chemical reaction mechanisms that involve bond breaking and forming.     

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.     

Quantum dynamics of inelastic scattering with a moving grid

The time‐dependent discrete variable representation (TDDVR) of a wave function with grid points defined by the Hermite part of the Gauss–Hermite (G‐H) basis set introduces quantum corrections to classical mechanics. The grid points in this method follow classical trajectory and the approach converges to the exact quantum formulation with sufficient trajectories (TDDVR points) but just with a single grid point; only classical mechanics performs the dynamics. This newly formulated approach (developed for handling time‐dependent molecular quantum dynamics) has been explored to calculate vibrational transitions in the inelastic scattering processes. Traditional quantum mechanical results exhibit an excellent agreement with TDDVR profiles during the entire propagation when enough grid points are included in the quantum‐classical dynamics. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Quantum tunneling effect in entanglement dynamics

Quantum tunneling effect in entanglement dynamics between two coupled particles with separable Gaussian initial state is investigated using entangled trajectory molecular dynamics method in terms of the reduced‐density linear entropy. It has been presented through showing distinguish contribution of single trajectory to linear entropy between classical trajectory and entangled trajectory with same initial state. We find that quantum tunneling effect makes single trajectory's contribution remarkably decrease under quantum dynamics compared to classical dynamics. The nonlocality of quantum entanglement is presented, and the energy transfer between two coupled particles through quantum correlations and classical ones is also discussed in the end. © 2015 Wiley Periodicals, Inc.     

Computing the density of paths in complex systems

Trajectories of duration tau joining two points q(0) and q(1) in the configuration space of a classical system satisfy Hamilton's principle: they are stationary points of the classical action. The second variation (fluctuations) of the action around the stationary point signals whether the latter is or not a minimum and delivers the density in trajectory space around the points q(0) and q(1). This concept of density of paths is of great importance in semiclassical quantum theory, where it weights the contribution to the propagator from the single classical trajectories. In this paper, two algorithms based on the concepts of molecular dynamics simulation are introduced for computing the density of paths, also called van Vleck [Proc. Natl. Acad. Sci. U.S.A. 14, 178 (1928)] determinant. Examples for realistic systems are presented, together with a suggestion about possible applications in the field of rare events in physics and chemistry.     

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.     

Semiclassical nonadiabatic dynamics based on quantum trajectories for the O(3P,1D) + H2 system

The O(3P,1D) + H2 --> OH + H reaction is studied using trajectory dynamics within the approximate quantum potential approach. Calculations of the wave-packet reaction probabilities are performed for four coupled electronic states for total angular momentum J = 0 using a mixed coordinate/polar representation of the wave function. Semiclassical dynamics is based on a single set of trajectories evolving on an effective potential-energy surface and in the presence of the approximate quantum potential. Population functions associated with each trajectory are computed for each electronic state. The effective surface is a linear combination of the electronic states with the contributions of individual components defined by their time-dependent average populations. The wave-packet reaction probabilities are in good agreement with the quantum-mechanical results. Intersystem crossing is found to have negligible effect on reaction probabilities summed over final electronic states.     

Comparisons of classical chemical dynamics simulations of the unimolecular decomposition of classical and quantum microcanonical ensembles

Previous studies have shown that classical trajectory simulations often give accurate results for short-time intramolecular and unimolecular dynamics, particularly for initial non-random energy distributions. To obtain such agreement between experiment and simulation, the appropriate distributions must be sampled to choose initial coordinates and momenta for the ensemble of trajectories. If a molecule's classical phase space is sampled randomly, its initial decomposition will give the classical anharmonic microcanonical (RRKM) unimolecular rate constant for its decomposition. For the work presented here, classical trajectory simulations of the unimolecular decomposition of quantum and classical microcanonical ensembles, at the same fixed total energy, are compared. In contrast to the classical microcanonical ensemble, the quantum microcanonical ensemble does not sample the phase space randomly. The simulations were performed for CH(4), C(2)H(5), and Cl(-)---CH(3)Br using both analytic potential energy surfaces and direct dynamics methods. Previous studies identified intrinsic RRKM dynamics for CH(4) and C(2)H(5), but intrinsic non-RRKM dynamics for Cl(-)---CH(3)Br. Rate constants calculated from trajectories obtained by the time propagation of the classical and quantum microcanonical ensembles are compared with the corresponding harmonic RRKM estimates to obtain anharmonic corrections to the RRKM rate constants. The relevance and accuracy of the classical trajectory simulation of the quantum microcanonical ensemble, for obtaining the quantum anharmonic RRKM rate constant, is discussed.     

Initial value representation for the SU(n) semiclassical propagator

The semiclassical propagator in the representation of SU(n) coherent states is characterized by isolated classical trajectories subjected to boundary conditions in a doubled phase space. In this paper, we recast this expression in terms of an integral over a set of initial-valued trajectories. These trajectories are monitored by a filter that collects only the appropriate contributions to the semiclassical approximation. This framework is suitable for the study of bosonic dynamics in n modes with fixed total number of particles. We exemplify the method for a Bose-Einstein condensate trapped in a triple-well potential, providing a detailed discussion on the accuracy and efficiency of the procedure.