Локально-адаптивный выбор шага интегрирования в задачах молекулярной динамики

Работа посвящена численным методам решения задач молекулярной динамики. Дан обзор существующих в настоящее время методов и алгоритмов. Сформулировано и доказано необходимое условие сходимости произвольной разностной схемы решения дифференциальных уравнений движения Ньютона. Из условия ограниченности интенсивности схемного источника энергии выведен критерий выбора временного шага в задаче двух тел. На основе этого критерия построен оригинальный алгоритм численного интегрирования задачи N тел, cочетающий в себе все современные подходы к построению новых схем. На примере двух модельных задач показано, что предлагаемый алгоритм обеспечивает как контроль за флуктуациями энергии, так и значительное ускорение счета.