Quantum-classical dynamics of scattering processes in adiabatic and diabatic representations

We demonstrate the workability of a TDDVR based [J. Chem. Phys. 118, 5302 (2003)], novel quantum-classical approach, for simulating scattering processes on a quasi-Jahn-Teller model [J. Chem. Phys. 105, 9141 (1996)] surface. The formulation introduces a set of DVR grid points defined by the Hermite part of the basis set in each dimension and allows the movement of grid points around the central trajectory. With enough trajectories (grid points), the method converges to the exact quantum formulation whereas with only one grid point, we recover the conventional molecular dynamics approach. The time-dependent Schrodinger equation and classical equations of motion are solved self-consistently and electronic transitions are allowed anywhere in the configuration space among any number of coupled states. Quantum-classical calculations are performed on diabatic surfaces (two and three) to reveal the effects of symmetry on inelastic and reactive state-to-state transition probabilities, along with calculations on an adiabatic surface with ordinary Born-Oppenheimer approximation. Excellent agreement between TDDVR and DVR results is obtained in both the representations.     

Semi‐classical formulation of time‐dependent discrete variable representation method

In this article, we apply a novel time‐dependent discrete variable representation (TDDVR) method proposed by Barkakaty and Adhikari to investigate tunneling through an Eckart barrier. This semi‐classical method is theoretically rigorous and straightforward to implement. Among the TDDVR formulations, this report presents the first derivation of a rigorous form of quantum force (QF) for the present perspective. The validity of this semi‐classical approach is demanded based on the excellent agreement of the tunneling probability with the corresponding quantum results. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

A “classical” trajectory driven nuclear dynamics by a parallelized quantum‐classical approach to a realistic model Hamiltonian of benzene radical cation

We explore the workability of a parallelized algorithm of time‐dependent discrete variable representation (TDDVR) methodology formulated by involving “classical” trajectories on each DOF of a multi‐mode multi‐state Hamiltonian to reproduce the population dynamics, photoabsorption spectra and nuclear dynamics of the benzene radical cation. To perform such dynamics, we have used a realistic model Hamiltonian consists of five lowest electronic states (X²E_1g, B²E_2g, C²A_2u, D²E_1u, and E²B_2u) which are interconnected through several conical intersections with nine vibrational modes. The calculated nuclear dynamics and photoabsorption spectra with the advent of our parallelized TDDVR approach show excellent agreement with the results obtained by multiconfiguration time‐dependent Hartree method and experimental findings, respectively. The major focus of this article is to demonstrate how the “classical” trajectories for the different modes and the “classical” energy functional for those modes on each surface can enlight the time‐dependent feature of nuclear density and its' nodal structure. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Multi-state multi-mode nuclear dynamics on three isomers of C6H4F2+ using parallelized TDDVR approach

We have performed molecular dynamics on the three isomers of the difluorobenzene radical cation (C(6)H(4)F(2)(+)) after excitation from the ground state to a specific higher electronically excited state by using our recently implemented parallelized Time-Dependent Discrete Variable Representation (TDDVR) methodology. A five-state eleven-mode realistic model Hamiltonian for o-C(6)H(4)F(2)(+) and two separate five-state ten-mode Hamiltonians for m- and p-isomer of the same radical cation are considered, where those five electronic states are interconnected through several conical intersections in the vicinity of the Franck-Condon (FC) region and thus the dynamics for each case become complex. The photoelectron, mass analyzed threshold ionization spectra and population profiles obtained by using our TDDVR approach show reasonably good agreement with the results obtained by multiconfiguration time dependent Hartree (MCTDH) method. It is worthwhile to mention that the parallelized TDDVR algorithm reduces the computation time by more than an order of magnitude compared to its serial analog and, therefore, such approach appears to be a good compromise between accuracy and speed for a large molecular system.     

Time‐dependent schrödinger equation with Markovian outgoing wave boundary conditions: Applications to quantum tunneling dynamics and photoionization

Markovian outgoing wave boundary conditions are introduced as an approximate method to reduce the size of the computational grid for time integration of the time‐dependent Schrödinger equation. The ratio and polynomial methods developed as open boundary conditions are applied to the wave function at the boundaries of the computational grid. This computational method is used to study the wave packet dynamics for a metastable well, a double well, and strong‐field ionization of a model atom. Accurate results demonstrate that this method can significantly reduce the number of grid points required in a dynamical calculation for quantum dynamical problems. © 2012 Wiley Periodicals, Inc.     

An approach for generating trajectory-based dynamics which conserves the canonical distribution in the phase space formulation of quantum mechanics. I. Theories

We have reformulated and generalized our recent work [J. Liu and W. H. Miller, J. Chem. Phys. 126, 234110 (2007)] into an approach for generating a family of trajectory-based dynamics methods in the phase space formulation of quantum mechanics. The approach (equilibrium Liouville dynamics) is in the spirit of Liouville's theorem in classical mechanics. The trajectory-based dynamics is able to conserve the quantum canonical distribution for the thermal equilibrium system and approaches classical dynamics in the classical (? → 0), high temperature (β → 0), and harmonic limits. Equilibrium Liouville dynamics provides the framework for the development of novel theoretical∕computational tools for studying quantum dynamical effects in large∕complex molecular systems.     

The effect of phonon modes on the H2(v, j)/D2(v, j)-Cu(1nn) scattering processes

We include the effect of the phonon modes originating from the three layers of Cu(1nn) surface atoms on the dynamics of incoming molecular [H(2)(v, j)/D(2)(v, j)] degrees of freedom (DOFs) through a mean-field approach, where the surface temperature is incorporated into the effective potential by considering Bose-Einstein probability (BEP) factor for the initial state distribution of the surface modes calculated within harmonic approximation. Such time and temperature dependent effective Hamiltonian is further subdivided assuming a weak coupling between the two sets of molecular DOFs, namely, (x, y, z, Z) and (X, Y), respectively, in particular, to reduce the computational cost and the corresponding coupled quantum dynamical equations of motion have been formulated in terms of Time Dependent Discrete Variable Representation (TDDVR) approach. We demonstrate the workability of TDDVR method to investigate the scattering of H(2)(v, j) on Cu(1nn) surface by calculating the reaction probabilities and scattering cross-sections. Calculated results show that the phonon modes affect (a) the state-to-state transition probabilities of the scattered H(2) molecule substantially but chemisorption and physisorption processes negligibly and (b) the reaction probability of the incoming D(2) molecule noticeably.     

Quantized Hamilton Dynamics

The paper describes the quantized Hamilton dynamics (QHD) approach that extends classical Hamiltonian dynamics and captures quantum effects, such as zero point energy, tunneling, decoherence, branching, and state-specific dynamics. The approximations are made by closures of the hierarchy of Heisenberg equations for quantum observables with the higher order observables decomposed into products of the lower order ones. The technique is applied to the vibrational energy exchange in a water molecule, the tunneling escape from a metastable state, the double-slit interference, the population transfer, dephasing and vibrational coherence transfer in a two-level system coupled to a phonon, and the scattering of a light particle off a surface phonon, where QHD is coupled to quantum mechanics in the Schrödinger representation. Generation of thermal ensembles in the extended space of QHD variables is discussed. QHD reduces to classical mechanics at the first order, closely resembles classical mechanics at the higher orders, and requires little computational effort, providing an efficient tool for treatment of the quantum effects in large systems.     

Time-dependent probability of quantum tunneling in terms of the quasisemiclassical method

In view of the rapid progress in experiments of the tunneling dynamics in the time domain, we develop a quasisemiclassical method that is aimed at a study of the proton-transfer dynamics in a large system such as tropolone and its interesting derivatives, to which not only full quantum mechanics, but even a standard semiclassical theory is never easy to apply. In our very tractable method for multidimensional systems, the tunneling paths are generated in terms of the generalized classical mechanics, but the quantum phases arising from the action integral, the Maslov index, and the semicalssical amplitude factor as well in the semiclassical kernels are entirely neglected. This approach is called the quasisemiclassical method. One of the technical issues involved in the general semiclassical scheme is how to locate points from which a tunneling path emanates. Hence the studies of such tunneling points and the quasisemiclassical method should be examined collectively. We test several ways of determining the tunneling point, including those already proposed in the literature and a newly proposed one. It is shown numerically that the quasisemiclassical method with an appropriate choice of tunneling points reproduces the full quantum mechanical tunneling probability reasonably well. This case study indicates that the present conventional approach is promising to the study of large systems. The role of tunneling points in the initial process of tunneling is also discussed.     

Tunnelling Dynamics of Kicked Systems: A Route to Tunnelling Control

The effects of discontinuously time-varying perturbations on the dynamics of a particle moving in harmonic, symmetric double well and symmetric triple well potentials, are investigated both classically and quantum mechanically. The quantum dynamics is followed using the time-dependent Fourier grid Hamiltonian (TDFGH) method while the classical dynamics is analyzed within the framework of classical Hamiltonian mechanics. Depending on the spatial symmetry of the perturbation and the characteristic features of the reversal time , different types of ‘phase space’ structures are observed in each of the potentials. For symmetric double and triple well potentials, quantum dynamics reveals that complete destruction of tunnelling (CDT) can be achieved in the presence of a time-dependent spatially asymmetric perturbing field that is continuous in time. Any discontinuity in time-variation of the perturbation may induce over the barrier transition. The relevance of these results in the context of (i) tunnelling control and (ii) quantum computing with 3-state or 2-state quantum registers is briefly discussed.     

Whether proton transition to the triphosphate tail of ATP occurs at protein kinase environment: A Car‐Parrinello ab initio molecular dynamics study

Protonation state of the triphosphate tail of ATP (adenosine triphosphate) in protein environment is a fundamental issue, which has significant impact on the mechanism investigation of biochemical processes with ATP involved. Proton transition from surroundings (water molecule coordinating to magnesium, H_W; amino group of Lys, H_L) to the ATP tail in the catalytic core of protein kinase found recently disproved the commonly accepted deprotonation state of ATP tail. In this account, Car‐Parrinello ab initio molecular dynamics (CP‐AIMD) method has been employed to examine whether the proton transition occurs. To provide a comparison basis for the dynamics simulations, static quantum mechanics (QM), and combined quantum mechanics and molecular mechanics (QM/MM) calculations have also been carried out. Consistent results have been obtained that complete transition of hydrogen from the surroundings to the triphosphate tail of ATP is not allowed. The most dominant conformations correspond to the ones with H_W bonding to O(W) and H‐bonding to O(ATP), [O(W)‐H_W···O(ATP)], H_L bonding to N(Lys) and H‐bonding to O(ATP), [N(Lys)‐H_L···O(ATP)]. Metastable structures with one hydrogen atom bonding with two heavy atoms (hydrogen acceptors) were also located by our dynamic simulations. This bonding mode can satisfy the hungering for hydrogen of the two heavy atoms simultaneously. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Elucidating collision induced dissociation products and reaction mechanisms of protonated uracil by coupling chemical dynamics simulations with tandem mass spectrometry experiments

In this study we have coupled mixed quantum‐classical (quantum mechanics/molecular mechanics) direct chemical dynamics simulations with electrospray ionization/tandem mass spectrometry experiments in order to achieve a deeper understanding of the fragmentation mechanisms occurring during the collision induced dissociation of gaseous protonated uracil. Using this approach, we were able to successfully characterize the fragmentation pathways corresponding to ammonia loss (m/z 96), water loss (m/z 95) and cyanic or isocyanic acid loss (m/z 70). Furthermore, we also performed experiments with isotopic labeling completing the fragmentation picture. Remarkably, fragmentation mechanisms obtained from chemical dynamics simulations are consistent with those deduced from isotopic labeling. Copyright © 2015 John Wiley & Sons, Ltd.     

A quantum-classical approach to the photoabsorption spectrum of pyrazine

We have used the time-dependent discrete variable representation (TDDVR) method to simulate the photoabsorption spectrum of pyrazine. The time-dependent molecular dynamics of pyrazine after excitation to the S2 electronic state is considered as a benchmark to investigate the S2 absorption spectrum. We have carried out the dynamics on a basic four-mode model of pyrazine with the inclusion of five major modes as well as the rest of the vibrational modes as bath modes. Investigations reveal the effect of bath modes such as energy and population transfer from the subsystem to the bath. Calculated results demonstrate excellent agreement with traditional quantum-mechanical findings during the entire propagation and converge to the exact quantum results when enough gridpoints are used. It appears that TDDVR, as a numerical quantum dynamics methodology, is a good compromise between accuracy and speed.     

Density dynamics in some quantum systems

The quantum domain behavior of classical nonintegrable systems is well‐understood by the implementation of quantum fluid dynamics and quantum theory of motion. These approaches properly explain the quantum analogs of the classical Kolmogorov–Arnold–Moser type transitions from regular to chaotic domain in different anharmonic oscillators. Field‐induced tunneling and chaotic ionization in Rydberg atoms are also analyzed with the help of these theories. Quantum fluid density functional theory may be used to understand different time‐dependent processes like ion‐atom/molecule collisions, atom‐field interactions, and so forth. Regioselectivity as well as confined atomic/molecular systems and their reactivity dynamics have also been explained. © 2013 Wiley Periodicals, Inc.     

Solution of the “Classical” Schrödinger equation for a driven symmetric triple well: A comparison with its classical counterpart

The behavior of a driven symmetric triple well potential has been studied by developing an algorithm where the well‐established Bohmian mechanics and time‐dependent Fourier Grid Hamiltonian method are incorporated and the quantum theory of motion (QTM) phase space structures of the particle are constructed, both in “nonclassical” and “classical” limits. Comparison of QTM phase space structures with their classical analogues shows both similarity as well as dissimilarities. The temporal nature and the spatial symmetry of applied perturbation play crucial roles in having similar phase space structures. © 2016 Wiley Periodicals, Inc.     

Exploring the importance of quantum effects in nucleation: the archetypical Ne(n) case

The effect of quantum mechanics (QM) on the details of the nucleation process is explored employing Ne clusters as test cases due to their semi-quantal nature. In particular, we investigate the impact of quantum mechanics on both condensation and dissociation rates in the framework of the microcanonical ensemble. Using both classical trajectories and two semi-quantal approaches (zero point averaged dynamics, ZPAD, and Gaussian-based time dependent Hartree, G-TDH) to model cluster and collision dynamics, we simulate the dissociation and monomer capture for Ne(8) as a function of the cluster internal energy, impact parameter and collision speed. The results for the capture probability P(s)(b) as a function of the impact parameter suggest that classical trajectories always underestimate capture probabilities with respect to ZPAD, albeit at most by 15%-20% in the cases we studied. They also do so in some important situations when using G-TDH. More interestingly, dissociation rates k(diss) are grossly overestimated by classical mechanics, at least by one order of magnitude. We interpret both behaviours as mainly due to the reduced amount of kinetic energy available to a quantum cluster for a chosen total internal energy. We also find that the decrease in monomer dissociation energy due to zero point energy effects plays a key role in defining dissociation rates. In fact, semi-quantal and classical results for k(diss) seem to follow a common "corresponding states" behaviour when the proper definition of internal and dissociation energies are used in a transition state model estimation of the evaporation rate constants.