A Family of Multistep Methods with Zero Phase-Lag and Derivatives for the Numerical Integration of Oscillatory ODEs

TL;DR

This paper introduces a family of optimized 8-step methods with nullified phase-lag and its derivatives, significantly improving the numerical integration of oscillatory ODEs like the Schrödinger equation and N-body problems.

Contribution

The paper develops a new family of multistep methods with nullified phase-lag derivatives, enhancing efficiency and stability for oscillatory differential equations.

Findings

Methods with more phase-lag derivatives nullified have larger stability intervals.

Nullifying phase-lag derivatives improves accuracy in oscillatory problems.

The methods effectively integrate the Schrödinger equation and N-body problems.

Abstract

In this paper we develop a family of three 8-step methods, optimized for the numerical integration of oscillatory ordinary differential equations. We have nullified the phase-lag of the methods and the first r derivatives, where r=1,2,3. We show that with this new technique, the method gains efficiency with each derivative of the phase-lag nullified. This is the case for the integration of both the Schrodinger equation and the N-body problem. A local truncation error analysis is performed, which, for the case of the Schrodinger equation, also shows the connection of the error and the energy, revealing the importance of the zero phase-lag derivatives. Also the stability analysis shows that the methods with more derivatives vanished, have a bigger interval of periodicity.