Independent trajectory implementation of the semiclassical Liouville method: application to multidimensional reaction dynamics

We describe an independent trajectory implementation of semiclassical Liouville method for simulating quantum processes using classical trajectories. In this approach, a single ensemble of trajectories describes all semiclassical density matrix elements of a coupled electronic state problem, with the ensemble evolving classically under a single reference Hamiltonian chosen on the basis of physical grounds. In this paper, we introduce an additional uncoupled trajectory approximation, allowing the members of the ensemble to evolve independently of one another and eliminating the major computational costs of our previous coupled trajectory implementation. The accuracy of the method is demonstrated for model one-dimensional problems. In addition, the approach is applied to the chemical reaction dynamics of a collinear triatomic system, yielding excellent agreement with exact calculations. This method allows molecular dynamics involving coupled electronic surfaces to be modeled with essentially the same effort as classical molecular dynamics and ensemble averaging.     

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.     

A semiclassical initial-value representation for quantum propagator and boltzmann operator

Starting from the position-momentum integral representation, we apply the correction operator method to the derivation of a uniform semiclassical approximation for the quantum propagator and then extend it to approximate the Boltzmann operator. In this approach, the involved classical dynamics is determined by the method itself instead of given beforehand. For the approximate Boltzmann operator, the corresponding classical dynamics is governed by a complex Hamiltonian, which can be described as a pair of real Hamiltonian systems. It is demonstrated that the semiclassical Boltzmann operator is exact for linear systems. A quantum propagator in the complex time is thus proposed and preliminary numerical results show that it is a reasonable approximation for calculating thermal correlation functions of general systems. © 2018 Wiley Periodicals, Inc.     

A centroid molecular dynamics study of liquid para-hydrogen and ortho-deuterium

Centroid molecular dynamics (CMD) is applied to the study of collective and single-particle dynamics in liquid para-hydrogen at two state points and liquid ortho-deuterium at one state point. The CMD results are compared with the results of classical molecular dynamics, quantum mode coupling theory, a maximum entropy analytic continuation approach, pair-product forward- backward semiclassical dynamics, and available experimental results. The self-diffusion constants are in excellent agreement with the experimental measurements for all systems studied. Furthermore, it is shown that the method is able to adequately describe both the single-particle and collective dynamics of quantum liquids.     

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.     

Nonadiabatic electron wavepacket dynamics of molecules in an intense optical field: an ab initio electronic state study

A theory of quantum electron wavepacket dynamics that nonadiabatically couples with classical nuclear motions in intense optical fields is studied. The formalism is intended to track the laser-driven electron wavepackets in terms of the linear combination of configuration-state functions generated with ab initio molecular orbitals. Beginning with the total quantum Hamiltonian for electrons and nuclei in the vector potential of classical electromagnetic field, we reduce the Hamiltonian into a mixed quantum-classical representation by replacing the quantum nuclear momentum operators with the classical counterparts. This framework gives equations of motion for electron wavepackets in an intense laser field through the time dependent variational principle. On the other hand, a generalization of the Newtonian equations provides a matrix form of forces acting on the nuclei for nonadiabatic dynamics. A mean-field approximation to the force matrix reduces this higher order formalism to the semiclassical Ehrenfest theory in intense optical fields. To bring these theories into a practical quantum chemical package for general molecules, we have implemented the relevant ab initio algorithms in it. Some numerical results in the level of the semiclassical Ehrenfest-type theory with explicit use of the nuclear kinematic (derivative) coupling and the velocity form for the optical interaction are presented.     

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.     

Are azobenzenophanes rotation-restricted?

We simulated the photoisomerization dynamics of an azobenzenophane with a semiclassical surface hopping approach and a semiempirical reparametrized quantum mechanics/molecular mechanics Hamiltonian. Only one of the two azobenzene chromophores in the molecule is taken into account quantum mechanically: the other one is treated by molecular mechanics. Both n-->pi* and pi-->pi* excitations are considered. Our results show that the photoisomerization reaction mainly involves the rotation around the N=N double bond. The excited state relaxation features are in qualitative agreement with experimental time-resolved fluorescence results.     

Spectroscopic interpretation: the high vibrations of CDBrClF

We extract the dynamics implicit in an algebraic fitted model Hamiltonian for the deuterium chromophore's vibrational motion in the molecule CDBrClF. The original model has four degrees of freedom, three positions and one representing interbond couplings. A conserved polyad allows in a semiclassical approach the reduction to three degrees of freedom. For most quantum states we can identify the underlying motion that when quantized gives the said state. Most of the classifications, identifications, and assignments are done by visual inspection of the already available wave function semiclassically transformed from the number representation to a representation on the reduced dimension toroidal configuration space corresponding to the classical action and angle variables. The concentration of the wave function density to lower dimensional subsets centered on idealized simple lower dimensional organizing structures and the behavior of the phase along such organizing centers already reveals the atomic motion. Extremely little computational work is needed.     

A semiclassical study of the thermal conductivity of low temperature liquids

The conventional classical energy current auto-correlation function has been extended into a quantum mechanical version and then approximated by the linearized semiclassical initial value representation approach. Comparison of the thermal conductivity to simulation results shows that about 15% quantum correction to the classical molecular dynamics results for liquid neon are quantitatively predicted. For liquid para-hydrogen the quantum effects are sufficiently large that the linearized semiclassical approach is only 20% accurate, while for both liquid He(4) and He(3) the thermal conductivity disagrees by a factor of 2, although exchange effects appear to play a minor role.     

Including quantum effects in the dynamics of complex (i.e., large) molecular systems

The development in the 1950s and 1960s of crossed molecular beam methods for studying chemical reactions at the single-collision molecular level stimulated the need and desire for theoretical methods to describe these and other dynamical processes in molecular systems. Chemical dynamics theory has made great strides in the ensuing decades, so that methods are now available for treating the quantum dynamics of small molecular systems essentially completely. For the large molecular systems that are of so much interest nowadays (e.g., chemical reactions in solution, in clusters, in nanostructures, in biological systems, etc.), however, the only generally available theoretical approach is classical molecular dynamics (MD) simulations. Much effort is currently being devoted to the development of approaches for describing the quantum dynamics of these complex systems. This paper reviews some of these approaches, especially the use of semiclassical approximations for adding quantum effects to classical MD simulations, also showing some new versions that should make these semiclassical approaches even more practical and accurate.     

Mixed quantum-classical simulations of charge transport in organic materials: numerical benchmark of the Su-Schrieffer-Heeger model

The electron-phonon coupling is critical in determining the intrinsic charge carrier and exciton transport properties in organic materials. In this study, we consider a Su-Schrieffer-Heeger (SSH) model for molecular crystals, and perform numerical benchmark studies for different strategies of simulating the mixed quantum-classical dynamics. These methods, which differ in the selection of initial conditions and the representation used to solve the time evolution of the quantum carriers, are shown to yield similar equilibrium diffusion properties. A hybrid approach combining molecular dynamics simulations of nuclear motion and quantum-chemical calculations of the electronic Hamiltonian at each geometric configuration appears as an attractive strategy to model charge dynamics in large size systems "on the fly," yet it relies on the assumption that the quantum carriers do not impact the nuclear dynamics. We find that such an approximation systematically results in overestimated charge-carrier mobilities, with the associated error being negligible when the room-temperature mobility exceeds ～4.8 cm(2)∕Vs (～0.14 cm(2)/Vs) in one-dimensional (two-dimensional) crystals.     

An accurate and simple quantum model for liquid water

The path-integral molecular dynamics and centroid molecular dynamics methods have been applied to investigate the behavior of liquid water at ambient conditions starting from a recently developed simple point charge/flexible (SPC/Fw) model. Several quantum structural, thermodynamic, and dynamical properties have been computed and compared to the corresponding classical values, as well as to the available experimental data. The path-integral molecular dynamics simulations show that the inclusion of quantum effects results in a less structured liquid with a reduced amount of hydrogen bonding in comparison to its classical analog. The nuclear quantization also leads to a smaller dielectric constant and a larger diffusion coefficient relative to the corresponding classical values. Collective and single molecule time correlation functions show a faster decay than their classical counterparts. Good agreement with the experimental measurements in the low-frequency region is obtained for the quantum infrared spectrum, which also shows a higher intensity and a redshift relative to its classical analog. A modification of the original parametrization of the SPC/Fw model is suggested and tested in order to construct an accurate quantum model, called q-SPC/Fw, for liquid water. The quantum results for several thermodynamic and dynamical properties computed with the new model are shown to be in a significantly better agreement with the experimental data. Finally, a force-matching approach was applied to the q-SPC/Fw model to derive an effective quantum force field for liquid water in which the effects due to the nuclear quantization are explicitly distinguished from those due to the underlying molecular interactions. Thermodynamic and dynamical properties computed using standard classical simulations with this effective quantum potential are found in excellent agreement with those obtained from significantly more computationally demanding full centroid molecular dynamics simulations. The present results suggest that the inclusion of nuclear quantum effects into an empirical model for water enhances the ability of such model to faithfully represent experimental data, presumably through an increased ability of the model itself to capture realistic physical effects.     

Real time quantum propagation on a Monte Carlo trajectory guided grids of coupled coherent states: 26D simulation of pyrazine absorption spectrum

In this work we apply the coupled coherent states technique of quantum molecular dynamics to simulation of the absorption spectrum of pyrazine. All 24 vibrational modes are taken into account. The nonadiabatic coupling obetween the S(1) and S(2) electronic states is treated by a mapping approach that adds two extra degrees of freedom to the effective vibronic Hamiltonian. The results are in a good agreement with experiment and with previous calculations by quantum multiconfigurational time dependent Hartree and semiclassical Herman-Kluk methods.     

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.     

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.     

First ultraviolet absorption band of methane: an ab initio study

Quantum mechanical calculations of the cross sections for photodissociation of CH4 and CD4 in the 1t2-->3s band are presented. The potential energy surfaces for the three states correlating with the 1 1T2 state at tetrahedral geometries are calculated. The elements of the (3x3) matrix representing the electronic Hamiltonian in the diabatic basis are expanded in powers of nuclear coordinates, up to the second order. The expansion coefficients are based on accurate multireference configuration interaction calculations. The electronically nonadiabatic dynamics is treated with the multiconfiguration time-dependent Hartree approach. All nine internal degrees of methane are included in the quantum dynamics simulations. The calculated cross section agrees well with experiment. Semiclassical calculations using the reflection principle suggest that the peaks in the spectrum correspond to the three adiabatic electronic states correlating with the 1 1T2 state at Td geometries. However, the non-Born-Oppenheimer terms in the Hamiltonian have a strong effect on the positions of the peaks in the absorption spectrum. The results of semiclassical calculations, which neglect these terms, are therefore quite different from the accurate quantum results and experiment.     

Optimized Monte Carlo sampling in forward-backward semiclassical dynamics

Forward-backward semiclassical dynamics (FBSD) provides a rigorous and powerful methodology for calculating time correlation functions in condensed phase systems characterized by substantial quantum mechanical effects associated with zero-point motion, quantum dispersion, or identical particle exchange symmetries. The efficiency of these simulations arises from the use of classical trajectories to capture all dynamical information. However, full quantization of the density operator makes these calculations rather expensive compared to fully classical molecular dynamics simulations. This article discusses the convergence properties of various correlation functions and introduces an optimal Monte Carlo sampling scheme that leads to a significant reduction of statistical error. A simple and efficient procedure for normalizing the FBSD results is also discussed. Illustrative examples on model systems are presented.