Adams-Bashforth predictor method

One of the most common numerical methods used in molecular dynamics to solve Newton s equations of motions is the Velocity Verlet integrator. This is typically implemented as a second order method, and we find that it can become numerically tmstable dtuing the course of hyperthermal collision events, where the atom velocities are often far from equilibrium. As an alternative, we have implemented a fifth/ sixth order predictor-corrector scheme for our calculations. Specifically, the driver we chose utilizes the Adams-Bashforth predictor method together with the Adams-Moulton corrector method for approximating the solution to the equations of motion. 

In 30 the authors have developed trigonometrically fitted Adams-Bashforth-Moulton predictor-corrector (P-C) methods. It is the first time in the literature that these methods are applied for the efficient solution of the resonance problem of the Schrodinger equation. The new trigonometrically fitted P-C schemes are based on the well known Adams-Bashforth-Moulton methods. In particular, they are based on the fourth order Adams-Bashforth scheme (as predictor) and on the fifth order Adams-Moulton scheme (as corrector). More... 

In 26 the authors have developed a new trigonometrically-fitted predictor-corrector (P-C) scheme based on the Adams-Bashforth-Moulton P-C methods. In particular, the predictor is based on the fifth algebraic order Adams-Bashforth scheme and the corrector on the sixth algebraic order Adams-Moulton scheme. More specifically the new developed scheme integrates exactly any linear combination of the functions ... 

In 37 the authors have developed trigonometrically fitted Adams-Bashforth-Moulton predictor-corrector (P-C) methods. The new trigonometrically fitted P-C schemes are based on the well known Adams-Bashforth-Moulton... 

The Adam-Bashforth methods are frequently used as predictors and the Adam-Moulton methods are often used as correctors. The combination of the two formulas results in predictor-corrector schemes. 

In practice, implicit multistep methods are used to improve upon approximations obtained by explicit methods. This combination is the so-called predictor-corrector method. Predictor-corrector methods employ a single-step method, such as the Runge-Kutta of order 4, to generate the starting values to an explicit method, such as an Adams-Bashforth. Then the approximation from the explicit method is improved upon by use of an implicit method, such as an Adams-Moulton method. Also, there are variable step size algorithms associated with the predictor-corrector strategy in the literature [5,25]. 

Using the Runge-Kutta order four method for starting, the Adams-Bashforth method as the predictor, and the Adams-Moulton method as a corrector, approximate the solution to... 

The quasiclassical trajectory method was used to study this system, and the variable step size modified Bulirsch-Stoer algorithm was specially developed for recombination problems such as this one. Comparisons were made with the fourth order Adams-Bashforth-Moulton predictor-corrector algorithm, and the modified Bulirsch-Stoer method was always more efficient, with the relative efficiency of the Bulirsch-Stoer method increasing as the desired accuracy increased. We measure the accuracy by computing the rms relative difference between the initial coordinates and momenta and their back-integrated values. For example, for a rms relative difference of 0.01, the ratio of the CPU times for the two methods was 1.6, for a rms relative difference of 0.001 it was 2.0, and for a rms relative difference of 10 it was 3.3. Another advantage of the variable step size method is that the errors in individual trajectories are more similar, e.g. a test run of ten trajectories yielded rms errors which differed by a factor of 53 when using the modified Bulirsch-Stoer... 

The well-known predictor-corrector Adams-Bashforth method of algebraic order four. 

The implicit multistep methods add stability but require more computation to evaluate the implicit part. In addition, the error coefficient of the Adams-Moulton method of order k is smaller than that of the Adams Bashforth method of the same order. As a consequence, the implicit methods should give improved accuracy. In fact, the error coefficient for the imphcit fourth-order Adams Moulton method is 19/720, and for the explicit fourth-order Adams Bashforth method it is 251/720. The difference is thus about an order of magnitude. Pairs of exphcit and implicit multistep methods of the same order are therefore often used as predictor-corrector pairs. In this case, the explicit method is used to calculate the solution,, at v +i. Furthermore, the imphcit method (corrector) uses y + to calculate /(x +i,y +i), which replaces /(x +i,y +i). This allows the solution, y +i, to be improved using the implicit method. The combination of the Adams Bashforth and the Adams Moulton methods as predictor orrector pairs is implemented in some ODE solvers. The Matlab odel 13 solver is an example of a variable-order Adams Bashforth Moulton multistep solver. 

The Adams-Bashforth methods use information about prior points. In principle, one can form polynomials using forward points as well. Using the points Xfc+i, Xk,. .., Xfcj, to form a - - 1 polynomial generates a class of methods known as Adams-Moulton Methods. However, in these methods also calculation of yk+i requires the solution of fk+i implicitly. Implicit methods are discussed separately in a section which deals with stiff equations. One can also use a combination of an implicit method, such as an Adams-Moulton method, along with an explicit method, like an Adams-Bashforth method, to form an explicit method known as the Predictor-Corrector Method, which is discussed below. 

---