An explicit numerical scheme for Milne's phase-amplitude equations

TL;DR

This paper introduces an explicit numerical method for Milne's phase-amplitude equations that improves stability by solving a linear third-order equation instead of the nonlinear amplitude equation, reducing error sensitivity.

Contribution

The paper presents a novel explicit numerical scheme based on a linear third-order equation, enhancing stability and accuracy in solving Milne's phase-amplitude equations.

Findings

The method reduces error sensitivity compared to previous approaches.

It avoids instability caused by rapidly varying amplitude components.

The linear equation approach facilitates more reliable numerical solutions.

Abstract

We propose an explicit numerical method to solve Milne's phase-amplitude equations. Previously proposed methods solve directly Milne's nonlinear equation for the amplitude. For that reason, they exhibit high sensitivity to errors and are prone to instability through the growth of a spurious, rapidly varying component of the amplitude. This makes the systematic use of these methods difficult. On the contrary, the present method is based on solving a linear third-order equation which is equivalent to the nonlinear amplitude equation. This linear equation was derived by Kiyokawa, who used it to obtain analytical results on Coulomb wavefunctions [Kiyokawa, AIP Advances, 2015]. The present method uses this linear equation for numerical computation, thus resolving the problem of the growth of a rapidly varying component.