Semiclassical dynamics in Wigner phase space II: Nonadiabatic hybrid Wigner dynamics.

We present an approximate semiclassical (SC) framework for mixed quantized dynamics in Wigner phase space in a two-part series. In the first article, we introduced the Adiabatic Hybrid Wigner Dynamics (AHWD) method that allows for a few important "system" degrees of freedom to be quantized using high-level double Herman-Kluk SC theory while describing the rest (the "bath") using classical-limit linearized SC theory. In this second article, we extend our hybrid Wigner dynamics to nonadiabatic processes.

Semiclassical dynamics in Wigner phase space I: Adiabatic hybrid Wigner dynamics.

The Wigner phase space formulation of quantum mechanics is a complete framework for quantum dynamic calculations that elegantly highlights connections with classical dynamics. In this series of two articles, building upon previous efforts, we derive the full hierarchy of approximate semiclassical (SC) dynamic methods for adiabatic and non-adiabatic problems in Wigner phase space. In Paper I, focusing on adiabatic single surface processes, we derive the well-known double Herman-Kluk (DHK) approximation for real-time correlation functions in Wigner phase space and connect it to the linearized SC (LSC) approximation through a stationary phase approximation.

GPU-accelerated on-the-fly nonadiabatic semiclassical dynamics.

GPU-accelerated on-the-fly nonadiabatic dynamics is enabled by interfacing the linearized semiclassical dynamics approach with the TeraChem electronic structure program. We describe the computational workflow of the "PySCES" code interface, a Python code for semiclassical dynamics with on-the-fly electronic structure, including parallelization over multiple GPU nodes. We showcase the abilities of this code and present timings for two benchmark systems: fulvene solvated in acetonitrile and a charge transfer system in which a photoexcited zinc-phthalocyanine donor transfers charge to a fullerene acceptor through multiple electronic states on an ultrafast timescale.

Mean-Field Ring Polymer Rates Using a Population Dividing Surface.

Mean-field ring polymer molecular dynamics offers a computationally efficient method for the simulation of reaction rates in multilevel systems. Previous work has established that, to model a nonadiabatic state-to-state reaction accurately, it is necessary to ensure reactive trajectories form kinked ring polymer configurations at the dividing surface. Building on this idea, we introduce a population difference coordinate and a reactive flux expression modified to only include contributions from kinked configurations.

A Linearized Semiclassical Dynamics Study of the Multiquantum Vibrational Relaxation of NO Scattering from a Au(111) Surface.

The vibrational relaxation of NO molecules scattering from a Au(111) surface has served as the focus of efforts to understand nonadiabatic energy transfer at metal-molecule interfaces. Experimental measurements and previous theoretical efforts suggest that multiquantal NO vibrational energy relaxation occurs via electron-hole pair excitations in the metal. Here, using a linearized semiclassical approach, we accurately predict the vibrational relaxation of NO from the ν = 3 state for different incident translational energies.