Numerical treatment of the time-dependent Schrodinger equation in rotating coordinates

The Schrodinger equation for an electron in the time-dependent field of two nuclei is solved numerically in a rotating coordinate system. Two space dimensions are treated by discretisation and the third by an analytical expansion in magnetic substates. Finite difference techniques together with the implicit Crank-Nicolson method for the time evolution are used to treat the resulting system of coupled equations. The method is applied to p+H scattering E lab =2 keV. Comparisons with experimental data and calculations using analytical expansions are made.