New optimized algorithms for molecular dynamics simulations
The method of molecular dynamics (MD) is a powerful tool for the prediction and investigation of various phenomena in physics, chemistry and biology. The development of efficient MD algorithms for integration of the equations of motion in classical and quantum many-body systems should therefore...
Збережено в:
Дата: | 2002 |
---|---|
Автори: | , , |
Формат: | Стаття |
Мова: | English |
Опубліковано: |
Інститут фізики конденсованих систем НАН України
2002
|
Назва видання: | Condensed Matter Physics |
Онлайн доступ: | http://dspace.nbuv.gov.ua/handle/123456789/120663 |
Теги: |
Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
|
Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
Цитувати: | New optimized algorithms for molecular dynamics simulations / I.P. Omelyan, I.M. Mryglod, R. Folk // Condensed Matter Physics. — 2002. — Т. 5, № 3(31). — С. 369-390. — Бібліогр.: 37 назв. — англ. |
Репозитарії
Digital Library of Periodicals of National Academy of Sciences of UkraineРезюме: | The method of molecular dynamics (MD) is a powerful tool for the prediction
and investigation of various phenomena in physics, chemistry and biology.
The development of efficient MD algorithms for integration of the equations
of motion in classical and quantum many-body systems should therefore
impact a lot of fields of fundamental research. In the present study it is
shown that most of the existing MD integrators are far from being ideal and
further significant improvement in the efficiency of the calculations can be
reached. As a result, we propose new optimized algorithms which allow
to reduce the numerical uncertainties to a minimum with the same overall
computational costs. The optimization is performed within the well recognized
decomposition approach and concerns the widely used symplectic
Verlet-, Forest-Ruth-, Suzuki- as well as force-gradient-based schemes. It
is concluded that the efficiency of the new algorithms can be achieved
better with respect to the original integrators in factors from 3 to 1000 for
orders from 2 to 12. This conclusion is confirmed in our MD simulations
of a Lennard-Jones fluid for a particular case of second- and fourth-order
integration schemes. |
---|