Abstract
A new family of P-stable two-step methods with minimal phase lag is developed for the numerical integration of the special second-order initial value problem. An application to the one-dimensional Schrödinger equation indicates that these new methods are generally more accurate than previously developed finite difference methods.