Основні чисельні алгоритми методу молекулярної динаміки

Алгоритм Верле “leapfrog”

2.
3.

Алгоритм Верле у швидкісній формі “velocity Verlet”

2.

Цей алгоритм можна виразити також

Алгоритми Верле - порівняння

(a): Verlet
(b): Half-Step (Leapfrog)
(c ): Velocity Verlet

Алгоритм Беемана

2.

Більш складний ніж алгоритм Верле у швидкісній формі
Та вимагає

Алгоритми типу предиктор-коректор

Предиктор : звичайний розклад у ряд Тейлора

Коректор : визначає

Алгоритми типу предиктор-коректор. Алгоритм Рунге-Кутта

2.
3.
4.

Алгоритми типу предиктор-коректор. Алгоритми Гіра

c0=1/6, c1=5/6, c2=1, c3=1/3 - алгоритм Гіра

Стабільність алгоритмів при виборі часового кроку

Кружечки: Verlet
Квадратики: Gear 4th order
Трикутники: Gear

Стабільність алгоритмів при довгих симуляціях