This item is only available for download by members of the University of Illinois community. Students, faculty, and staff at the U of I may log in with your NetID and password to view the item. If you are trying to access an Illinois-restricted dissertation or thesis, you can request a copy through your library's Inter-Library Loan office or purchase a copy directly from ProQuest.
Permalink
https://hdl.handle.net/2142/81958
Description
Title
Longer Time Steps for Molecular Dynamics
Author(s)
Izaguirre Paz, Jesus Antonio
Issue Date
1999
Doctoral Committee Chair(s)
Skeel, Robert D.
Department of Study
Computer Science
Discipline
Computer Science
Degree Granting Institution
University of Illinois at Urbana-Champaign
Degree Name
Ph.D.
Degree Level
Dissertation
Keyword(s)
Applied Mechanics
Language
eng
Abstract
"In molecular dynamics one solves Newton's equations of motion. This is a challenging problem because the number of particles is large (thousands of atoms) but more importantly because stability severely limits the length of time steps relative to the total length needed for simulations---time steps are in the order of femtoseconds (10-15 seconds) whereas simulations of a few microseconds (10-6 seconds) up to one second are most desired. This situation also holds for multiple time stepping (MTS) integrators which evaluate different parts of the force at different frequencies. This work is an attempt to construct MTS integrators which are stable for long time steps. This has been done by applying a seminal idea that appeared in [GASSS98]: modify the potential by defining it at time-averaged positions Ax , where the time averaging Ax takes into account the high frequency vibrational motion. The resulting family of numerical integrators is called the mollified impulse method or MOLLY. Several averagings suggested, but not tested, in [GASSS98] are evaluated in this dissertation, along with more successful new ones. Finally, an efficient implementation of these averagings is devised. One of them is included in a production molecular dynamics code (NAMD 2). It is shown that the limitations in current MOLLY methods stem from a deficient representation of the eigenspace of the fast part of the forces. An alternative, less rigorous, approach based on ""Langevin stabilization"" is explored. A mild damping and random noise is added to the system to remove instabilities without significantly affecting its dynamics."
Use this login method if you
don't
have an
@illinois.edu
email address.
(Oops, I do have one)
IDEALS migrated to a new platform on June 23, 2022. If you created
your account prior to this date, you will have to reset your password
using the forgot-password link below.
