Abstract:Oscillatory second order linear ordinary differential equations arise in many scientific calculations. Because the running times of standard solvers increase linearly with frequency when they are applied to such problems, a variety of specialized methods, most of them quite complicated, have been proposed. Here, we point out that one of the simplest approaches not only works, but yields a scheme for solving oscillatory second order linear ordinary differential equations which is significantly faster than current state-of-the-art techniques. Our method, which operates by constructing a slowly varying phase function representing a basis of solutions of the differential equation, runs in time independent of the frequency and can be applied to second order equations whose solutions are oscillatory in some regions and slowly varying in others. In the high-frequency regime, our algorithm discretizes the nonlinear Riccati equation satisfied by the derivative of the phase function via a Chebyshev spectral collocation method and applies the Newton-Kantorovich method to the resulting system of nonlinear algebraic equations. We prove that the iterates converge quadratically to a nonoscillatory solution of the Riccati equation. The quadratic convergence of the Newton-Kantorovich method and the simple form of the linearized equations ensure that this procedure is extremely efficient. Our algorithm then extends the slowly varying phase function calculated in the high-frequency regime throughout the solution domain by solving a certain third order linear ordinary differential equation related to the Riccati equation. We describe the results of numerical experiments showing that our algorithm is orders of magnitude faster than existing schemes, including the modified Magnus method [18], the current state-of-the-art approach [7] and the recently introduced ARDC method [1].

A generator of families of two-step numerical methods with free parameters and minimal phase-lag

A new family of two stage symmetric two-step methods with vanished phase-lag and its derivatives for the numerical integration of the Schrödinger equation

Family of symmetric linear six-step methods with vanished phase-lag and its derivatives and their application to the radial Schrödinger equation and related problems

Optimizing a Hybrid Two-Step Method for the Numerical Solution of the Schrödinger Equation and Related Problems with Respect to Phase-Lag

A Phase-Fitted Runge-Kutta-Nyström method for the Numerical Solution of Initial Value Problems with Oscillating Solutions

A New Methodology for the Development of Efficient Multistep Methods for First-Order IVPs with Oscillating Solutions

A family of iterative methods to solve nonlinear problems with applications in fractional differential equations

Uniformly accurate numerical schemes for highly oscillatory Klein-Gordon and nonlinear Schrödinger equations

Two-step Extended RKN Methods for Oscillatory Systems.

A numerical damped oscillator approach to constrained Schrödinger equations

A low-cost, two-step fourteenth-order phase-fitting approach to tackling problems in chemistry

A Parametric Six-Step Method for Second-Order IVPs with Oscillating Solutions

Efficient Multistep Algorithms for First-Order IVPs with Oscillating Solutions: II Implicit and Predictor–Corrector Algorithms

Numerical solution of highly oscillatory ordinary differential equations

Oscillation Criteria for Even-Order Nonlinear Dynamic Equations with Sublinear and Superlinear Neutral Terms on Time Scales

On the Simulations of Second-Order Oscillatory Problems with Applications to Physical Systems

Symmetric and arbitrarily high-order time-stepping methods for solving nonlinear Klein–Gordon equations

Functionally-fitted energy-preserving methods for solving oscillatory nonlinear Hamiltonian systems

Extended RKN-type methods with minimal dispersion error for perturbed oscillators

An accelerated frequency-independent solver for oscillatory differential equations

An Adaptive Spectral Method for Oscillatory Second-Order Linear ODEs with Frequency-Independent Cost