Symmetrical Windowing for Quantum States in Quasi-Classical Trajectory Simulations

A microscopically reversible approach toward computing reaction probabilities via classical trajectory simulation has been developed that bins trajectories symmetrically on the basis of their initial and final classical actions. The symmetrical quasi-classical (SQC) approach involves defining a classical action window function centered at integer quantum values of the action, choosing a width parameter that is less than unit quantum width, and applying the window function to both initial reactant and final product vibrational states. Calculations were performed using flat histogram windows and Gaussian windows over a range of width parameters. Use of the Wigner distribution function was also investigated as a possible choice. It was demonstrated for collinear H + H2 reactive scattering on the BKMP2 potential energy surface that reaction probabilities computed via the SQC methodology using a Gaussian window function of 1/2 unit width produces good agreement with quantum mechanical results over the 0.4-0.6 eV energy range relevant to the ground vibrational state to the ground vibrational state reactive transition.