Ensembles of quantum trajectories are evolved to study time-dependent reaction dynamics in multidimensional systems with up to 25 vibrational modes. The equations of motion are formulated in curvilinear reaction path coordinates and all coupling terms are retained, including those involving curvature of the reaction path. The model potential is a Gaussian barrier along the translational coordinate coupled to $M$ vibrational modes. Spatial derivatives needed to propagate the trajectories are evaluated by least squares fitting in a contracted basis set. Stable propagation of the trajectory ensembles was carried out until complete bifurcation into reflected and reactive subensembles. The reaction probabilities were evaluated by Monte Carlo integration of the multidimensional smooth transmitted densities. Computational results, including trajectory plots and time-dependent reaction probabilities, are presented for $M=1$, 5, and 25 vibrational modes.

• Received 2 December 2005

DOI:https://doi.org/10.1103/PhysRevE.73.046701