Quasiclassical trajectory studies using 3D spline interpolation ofabinitiosurfaces

摘要

The accuracy of spline fitting procedures has been investigated for Morse potentials, Lennardhyphen;Jones potentials, a collinear Dndash;Clndash;H potential surface, and for three dimensional DIM (Hendash;H+2) and a valencehyphen;bond Dndash;Clndash;H surface. The adequacy of 3D spline surfaces in quasiclassical trajectory studies has been examined for both the Hendash;H+2and D+HCl systems. It is found that while one dimensional spline fits are very accurate, this accuracy decreases substantially for a twohyphen;dimensional spline fit. There is an additional, but smaller, decrease in accuracy for threehyphen;dimensional splines. The surface gradients are found to be less accurately fitted than the corresponding function values. In general, (15times;15times;15) cubic spline fits are found to lack the necessary accuracy to produce a pointhyphen;byhyphen;point match of a quasiclassical trajectory to that obtained on the original analytic surface. However, total reaction cross sections, energy partitioning distributions, and spatial scattering distributions computed on spline surfaces are found to be in good accord with those obtained from the full analytic surface. This suggests that in spite of their limited accuracy, spline fits may provide a very useful tool for interpolation ofabinitiosurfaces. The computer time requirements for use of 3D cubic splines in quasiclassical trajectory calculations are found to be within range of present computational facilities.