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

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

Quantum wavepacket ab initio molecular dynamics: an approach for computing dynamically averaged vibrational spectra including critical nuclear quantum effects

We have introduced a computational methodology to study vibrational spectroscopy in clusters inclusive of critical nuclear quantum effects. This approach is based on the recently developed quantum wavepacket ab initio molecular dynamics method that combines quantum wavepacket dynamics with ab initio molecular dynamics. The computational efficiency of the dynamical procedure is drastically improved (by several orders of magnitude) through the utilization of wavelet-based techniques combined with the previously introduced time-dependent deterministic sampling procedure measure to achieve stable, picosecond length, quantum-classical dynamics of electrons and nuclei in clusters. The dynamical information is employed to construct a novel cumulative flux/velocity correlation function, where the wavepacket flux from the quantized particle is combined with classical nuclear velocities to obtain the vibrational density of states. The approach is demonstrated by computing the vibrational density of states of [Cl-H-Cl]-, inclusive of critical quantum nuclear effects, and our results are in good agreement with experiment. A general hierarchical procedure is also provided, based on electronic structure harmonic frequencies, classical ab initio molecular dynamics, computation of nuclear quantum-mechanical eigenstates, and employing quantum wavepacket ab initio dynamics to understand vibrational spectroscopy in hydrogen-bonded clusters that display large degrees of anharmonicities.     

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 wave packet ab initio molecular dynamics: an approach to study quantum dynamics in large systems

A methodology to efficiently conduct simultaneous dynamics of electrons and nuclei is presented. The approach involves quantum wave packet dynamics using an accurate banded, sparse and Toeplitz representation for the discrete free propagator, in conjunction with ab initio molecular dynamics treatment of the electronic and classical nuclear degree of freedom. The latter may be achieved either by using atom-centered density-matrix propagation or by using Born-Oppenheimer dynamics. The two components of the methodology, namely, quantum dynamics and ab initio molecular dynamics, are harnessed together using a time-dependent self-consistent field-like coupling procedure. The quantum wave packet dynamics is made computationally robust by using adaptive grids to achieve optimized sampling. One notable feature of the approach is that important quantum dynamical effects including zero-point effects, tunneling, as well as over-barrier reflections are treated accurately. The electronic degrees of freedom are simultaneously handled at accurate levels of density functional theory, including hybrid or gradient corrected approximations. Benchmark calculations are provided for proton transfer systems and the dynamics results are compared with exact calculations to determine the accuracy of the approach.     

Exploring dynamical electron theory beyond the Born-Oppenheimer framework: from chemical reactivity to non-adiabatically coupled electronic and nuclear wavepackets on-the-fly under laser field

Chemical theory and its application to dynamical electrons in molecules under intense electromagnetic fields is explored, in which we take an explicit account of nuclear nonadiabatic (kinematic) interactions along with simultaneous coupling with intense optical interactions. All the electronic wavefunctions studied here are necessarily time-dependent, and thereby beyond stationary state quantum chemistry based on the Born-Oppenheimer framework. As a general and tractable alternative framework with which to track the electronic and nuclear simultaneous dynamics, we propose an on-the-fly method to calculate the electron and nuclear wavepackets coupled along the branching non-Born-Oppenheimer paths, through which their bifurcations, strong quantum entanglement between nuclear electronic motions, and coherence and decoherence among the phases associated with them are properly represented. Some illustrative numerical examples are also reported, which are aimed at our final goals; real time tracking of nonadiabatic electronic states, chemical dynamics in densely degenerate electronic states coupled with nuclear motions and manipulation and/or creation of new electronic states in terms of intense lasers, and so on. Other examples are also presented as to how the electron wavepacket dynamics can be used to analyze chemical reactions, shedding a new light on some typical and conventional chemical reactions such as proton transfer followed by tautomerization.     

Nonadiabatic chemical dynamics in an intense laser field: electronic wave packet coupled with classical nuclear motions

Dynamics of molecules in an intense laser field is studied in terms of the quantum electronic wave packet coupled with classical nuclear motions. The equations of motion are derived taking a proper account of molecular interactions with the vector potential of a classical electromagnetic field, along with the nonadiabatic interaction due to the breakdown of the Born-Oppenheimer approximation. With the aid of electronic structure calculations, the present method enables us to track, in an ab initio manner, the dynamics of polyatomic molecules in an intense field. Preliminary calculations are carried out for the vibrational state of LiF and a collision of Li+F under an intense laser pulse, which are limited to the domain of no ionization.     

Quantum wavepacket ab initio molecular dynamics: generalizations using an extended Lagrangian treatment of diabatic states coupled through multireference electronic structure

We present a generalization to our previously developed quantum wavepacket ab initio molecular dynamics (QWAIMD) method by using multiple diabatic electronic reduced single particle density matrices, propagated within an extended Lagrangian paradigm. The Slater determinantal wavefunctions associated with the density matrices utilized may be orthogonal or nonorthogonal with respect to each other. This generalization directly results from an analysis of the variance in electronic structure with quantum nuclear degrees of freedom. The diabatic electronic states are treated here as classical parametric variables and propagated simultaneously along with the quantum wavepacket and classical nuclei. Each electronic density matrix is constrained to be N-representable. Consequently two sets of new methods are derived: extended Lagrangian-QWAIMD (xLag-QWAIMD) and diabatic extended Lagrangian-QWAIMD (DxLag-QWAIMD). In both cases, the instantaneous potential energy surface for the quantum nuclear degrees of freedom is constructed from the diabatic states using an on-the-fly nonorthogonal multireference formalism. By introducing generalized grid-based electronic basis functions, we eliminate the basis set dependence on the quantum nucleus. Subsequent reuse of the two-electron integrals during the on-the-fly potential energy surface computation stage yields a substantial reduction in computational costs. Specifically, both xLag-QWAIMD and DxLag-QWAIMD turn out to be about two orders of magnitude faster than our previously developed time-dependent deterministic sampling implementation of QWAIMD. Energy conservation properties, accuracy of the associated potential surfaces, and vibrational properties are analyzed for a family of hydrogen bonded systems.     

On-the-fly localization of electronic orbitals in Car-Parrinello molecular dynamics

The ab initio molecular-dynamics formalism of Car and Parrinello is extended to preserve the locality of the orbitals. The supplementary term in the Lagrangian does not affect the nuclear dynamics, but ensures "on the fly" localization of the electronic orbitals within a periodic supercell in the Gamma-point approximation. The relationship between the resulting equations of motion and the formation of a gauge-invariant Lagrangian combined with a gauge-fixing procedure is briefly discussed. The equations of motion can be used to generate a very stable and easy to implement numerical integration algorithm. It is demonstrated that this algorithm can be used to compute the trajectory of the maximally localized orbitals, known as Wannier orbitals, in ab initio molecular dynamics with only a modest increase in the overall computer time. In the present paper, the new method is implemented within the generalized gradient approximation to Kohn-Sham density-functional theory employing plane wave basis sets and atomic pseudopotentials. In the course of the presentation, we briefly discuss how the present approach can be combined with localized basis sets to design fast linear scaling ab initio molecular-dynamics methods.     

First principles semiclassical calculations of vibrational eigenfunctions

Vibrational eigenfunctions are calculated on-the-fly using semiclassical methods in conjunction with ab initio density functional theory classical trajectories. Various semiclassical approximations based on the time-dependent representation of the eigenfunctions are tested on an analytical potential describing the chemisorption of CO on Cu(100). Then, first principles semiclassical vibrational eigenfunctions are calculated for the CO(2) molecule and its accuracy evaluated. The multiple coherent states initial value representations semiclassical method recently developed by us has shown with only six ab initio trajectories to evaluate eigenvalues and eigenfunctions at the accuracy level of thousands trajectory semiclassical initial value representation simulations.     

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.     

Efficient evaluation of the accuracy of molecular quantum dynamics on an approximate analytical or interpolated ab initio potential energy surface

Ab initio electronic structure methods have reached a satisfactory accuracy for the calculation of static properties, but remain too expensive for quantum dynamical calculations. Recently, an efficient semiclassical method was proposed to evaluate the accuracy of quantum dynamics on an approximate potential without having to perform the expensive quantum dynamics on the accurate potential. Here, this method is applied for the first time to evaluate the accuracy of quantum dynamics on an approximate analytical or interpolated potential in comparison to the quantum dynamics on an accurate potential obtained by an ab initio electronic structure method. Specifically, the vibrational dynamics of H₂ on a Morse potential is compared with that on the full CI potential, and the photodissociation dynamics of CO₂ on a LEPS potential with that on the excited ¹Π surface computed at the EOM‐CCSD/aug‐cc‐pVDZ level of theory. Finally, the effect of discretization of a potential energy surface on the quantum dynamics is evaluated. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110:2426–2435, 2010     

Herman-Kluk semiclassical dynamics in action-angle representation: new approaches to mapping quantum degrees of freedom

A general approach to mapping a discrete quantum mechanical problem by a continuous Hamiltonian is presented. The method is based on the representation of the quantum number by a continuous action variable that extends from -infinity to infinity. The projection of this Hilbert space onto the set of integer quantum numbers reduces the Hamiltonian to a discrete matrix of interest. The theory allows the application of the semiclassical methods to discrete quantum mechanical problems and, in particular, to problems where quantum Hamiltonians are coupled to continuous degrees of freedom. The Herman Kluk semiclassical propagation is used to calculate the nonadiabatic dynamics for a model avoided crossing system. The results demonstrate several advantages of the new theory compared to the existing mapping approaches.     

Ab initio quantum mechanical charge field (QMCF) molecular dynamics: a QM/MM – MD procedure for accurate simulations of ions and complexes

A new formalism for quantum mechanical / molecular mechanical (QM/MM) dynamics of chemical species in solution has been developed, which does not require the construction of any other potential functions except those for solvent–solvent interactions, maintains all the advantages of large simulation boxes and ensures the accuracy of ab initio quantum mechanics for all forces acting in the chemically most relevant region. Interactions between solute and more distant solvent molecules are incorporated by a dynamically adjusted force field corresponding to the actual molecular configuration of the simulated system and charges derived from the electron distribution in the solvate. The new formalism has been tested with some examples of hydrated ions, for which accurate conventional ab initio QM/MM simulations have been previously performed, and the comparison shows equivalence and in some aspects superiority of the new method. As this simulation procedure does not require any tedious construction of two-and three-body interaction potentials inherent to conventional QM/MM approaches, it opens the straightforward access to ab initio molecular dynamics simulations of any kind of solutes, such as metal complexes and other composite species in solution.     

Classical dynamics of 3D Hydrogen molecular ion in intense laser fields

The classical trajectory method is used to study the dynamics of 3D Hydrogen molecular ion interacting with intense laser fields. In the 3D classical model, a three-body Hamiltonian with one-dimensional nuclear motion restricted to the direction of the laser field is considered. The motion of electron and nucleus is described by the classical Hamiltonian canonical equations. The probabilities of ionization, dissociation and Coulomb explosion as functions of time are calculated and the average distances from electron to the mass-center for various laser parameters are implemented by symplectic method. The dynamics of in two-color laser fields are also investigated. We compare our results with the corresponding quantum-mechanical calculations and find they produce similar qualitative features in many cases.     

Multistage ab initio quantum wavepacket dynamics for electronic structure and dynamics in open systems: momentum representation, coupled electron-nuclear dynamics, and external fields

We recently proposed a multistage ab initio wavepacket dynamics (MS-AIWD) treatment for the study of delocalized electronic systems as well as electron transport through donor-bridge-acceptor systems such as those found in molecular-wire/electrode networks. In this method, the full donor-bridge-acceptor open system is treated through a rigorous partitioning scheme that utilizes judiciously placed offsetting absorbing and emitting boundary conditions. In this manner, the electronic coupling between the bridge molecule and surrounding electrodes is accounted. Here, we extend MS-AIWD to include the dynamics of open-electronic systems in conjunction with (a) simultaneous treatment of nuclear dynamics and (b) external electromagnetic fields. This generalization is benchmarked through an analysis of wavepackets propagated on a potential modeled on an Al(27) - C(7) - Al(27) nanowire. The wavepacket results are inspected in the momentum representation and the dependence of momentum of the wavepacket as well as its transmission probabilities on the magnitude of external bias are analyzed.     

Semiclassical many-mode floquet theory

It is shown that the single-mode Floquet formalism of Shirley can be extended to a generalized many-mode Floquet theory, yielding a practical non-perturbative technique for the semiclassical treatment of the interaction of a quantum system several monochromatic oscillating fields. The theory is illustrated by a detailed study of the population dynamics of a three-level system driven by two monochromatic radiation fields.     

Matching-pursuit split-operator Fourier-transform simulations of excited-state intramolecular proton transfer in 2-(2'-hydroxyphenyl)-oxazole

The excited-state intramolecular proton-transfer dynamics associated with the keto-enolic tautomerization reaction in 2-(2(')-hydroxyphenyl)-oxazole is simulated according to a numerically exact quantum-dynamics propagation method and a full-dimensional excited-state potential energy surface, based on an ab initio reaction surface Hamiltonian. The reported simulations involve the propagation of 35-dimensional wave packets according to the recently developed matching-pursuit/split-operator-Fourier-transform (MP/SOFT) method by Wu and Batista. The underlying propagation scheme recursively applies the time-evolution operator as defined by the Trotter expansion to second order accuracy in dynamically adaptive coherent-state expansions. Computations of time-dependent survival amplitudes, photoabsorption cross sections, and time-dependent reactant(product) populations are compared to the corresponding calculations based on semiclassical approaches, including the Herman-Kluk semiclassical initial value representation method. The reported results demonstrate the capabilities of the MP/SOFT method as a valuble computational tool to study ultrafast reaction dynamics in polyatomic systems as well as to validate semiclassical simulations of complex (nonintegrable) quantum dynamics in multidimensional model systems.     

Direct assessment of quantum nuclear effects on hydrogen bond strength by constrained-centroid ab initio path integral molecular dynamics

The impact of quantum nuclear effects on hydrogen (H-) bond strength has been inferred in earlier work from bond lengths obtained from path integral molecular dynamics (PIMD) simulations. To obtain a direct quantitative assessment of such effects, we use constrained-centroid PIMD simulations to calculate the free energy changes upon breaking the H-bonds in dimers of HF and water. Comparing ab initio simulations performed using PIMD and classical nucleus molecular dynamics (MD), we find smaller dissociation free energies with the PIMD method. Specifically, at 50 K, the H-bond in (HF)(2) is about 30% weaker when quantum nuclear effects are included, while that in (H(2)O)(2) is about 15% weaker. In a complementary set of simulations, we compare unconstrained PIMD and classical nucleus MD simulations to assess the influence of quantum nuclei on the structures of these systems. We find increased heavy atom distances, indicating weakening of the H-bond consistent with that observed by direct calculation of the free energies of dissociation.