RESUMO
Operator-difference multilayer schemes for solving the time-dependent Schrödinger equation up to sixth order of accuracy in the time step are presented. Reduced schemes for solving a set of coupled time-dependent Schrödinger equations with respect to the hyper-radial variable are devised using expansion of a wave packet over the set of appropriate basis angular functions. Further discretization of the resulting problem is realized by means of the finite-element method. The convergence of the expansion with respect to the number of basis functions and the efficiency of the numerical schemes are demonstrated in the exactly solvable model of an electric-field-driven two-dimensional oscillator (or a charged particle in a constant uniform magnetic field), in which we explicitly observed an effect of the periodical focusing and defocusing of the probability density flux.