Bipolar Reaction Path Hamiltonian Approach for Reactive Scattering Problems

Moving Boundary Truncated Grid Method: Application to the Time Evolution of Distribution Functions in Phase Space

The moving boundary truncated grid (TG) method, previously developed to integrate the time-dependent Schrödinger equation and the imaginary time Schrödinger equation, is extended to the time evolution of distribution functions in phase space. A variable number of phase space grid points in the Eulerian representation are used to integrate the equation of motion for the distribution function, and the boundaries of the TG are adaptively determined as the distribution function evolves in time. Appropriate grid points are activated and deactivated for propagation of the distribution function, and no advance information concerning the dynamics in phase space is required. The TG method is used to integrate the equations of motion for phase space distribution functions, including the Klein-Kramers, Wigner-Moyal, and modified Caldeira-Leggett equations. Even though the initial distribution function is nonnegative, the solutions to the Wigner-Moyal and modified Caldeira-Leggett equations may develop negative basins in phase space originating from interference effects. Trajectory-based methods for propagation of the distribution function do not permit the formation of negative regions. However, the TG method can correctly capture the negative basins. Comparisons between the computational results obtained from the full grid and TG calculations demonstrate that the TG method not only significantly reduces the computational effort but also permits accurate propagation of various distribution functions in phase space.

Complex quantum Hamilton-Jacobi equation with Bohmian trajectories: application to the photodissociation dynamics of NOCl

The complex quantum Hamilton-Jacobi equation-Bohmian trajectories (CQHJE-BT) method is introduced as a synthetic trajectory method for integrating the complex quantum Hamilton-Jacobi equation for the complex action function by propagating an ensemble of real-valued correlated Bohmian trajectories. Substituting the wave function expressed in exponential form in terms of the complex action into the time-dependent Schrödinger equation yields the complex quantum Hamilton-Jacobi equation. We transform this equation into the arbitrary Lagrangian-Eulerian version with the grid velocity matching the flow velocity of the probability fluid. The resulting equation describing the rate of change in the complex action transported along Bohmian trajectories is simultaneously integrated with the guidance equation for Bohmian trajectories, and the time-dependent wave function is readily synthesized. The spatial derivatives of the complex action required for the integration scheme are obtained by solving one moving least squares matrix equation. In addition, the method is applied to the photodissociation of NOCl. The photodissociation dynamics of NOCl can be accurately described by propagating a small ensemble of trajectories. This study demonstrates that the CQHJE-BT method combines the considerable advantages of both the real and the complex quantum trajectory methods previously developed for wave packet dynamics.