qunMaThesis.pdf (1.31 MB)
Novel Multiscale Algorithms for Molecular Dynamics
thesis
posted on 2003-07-03, 00:00 authored by Qun MaIn post-genomic computational biology and bioinformatics, long simulations of the dynamics of molecular systems, particularly biological molecules such as proteins and DNA, require advances in time stepping computational methods. The most severe problem of these algorithms is instability. The objective of this dissertation is to present original work in constructing multiscale multiple time stepping (MTS) algorithms for molecular dynamics (MD) that allow large time steps. First, through nonlinear stability analysis and numerical experiments, we reveal that MTS integrators such as Impulse suffer nonlinear overheating when $Delta t = T/3$ or possibly $Delta t = T/4$ when constant--energy MD simulations are attempted, where $Delta t$ is the longest step size and $T$ is the shortest period of the modes in the system. Second, we present Targeted MOLLY (TM), a new multiscale integrator for MD simulations. TM combines an efficient implementation of B-spline MOLLY exploiting analytical Hessians of energies and a self--consistent dissipative leapfrog integrator. Results show that TM allows very large time steps for slow forces (and thus multiscale) for the numerically challenging flexible TIP3P water systems (Jorgensen, {it et al.} J. Chem. Phys., vol 79, pp 926--935, 1983) while still computing the dynamical and structural properties accurately. Finally, we show yet another new MOLLY integrator, the Backward Euler (BE) MOLLY in which hydrogen bond forces can easily be included in the averaging and thus stability might be further improved.
History
Date Created
2003-07-03Date Modified
2018-10-05Defense Date
2003-06-26Research Director(s)
Jesus A IzaguirreCommittee Members
Menelaos I Karavelas Ryan K Roeder Gregory Madey J Daniel GezelterDegree
- Doctor of Philosophy
Degree Level
- Doctoral Dissertation
Language
- English
Alternate Identifier
etd-07032003-150352Publisher
University of Notre DameProgram Name
- Computer Science and Engineering
Usage metrics
Categories
No categories selectedKeywords
Licence
Exports
RefWorks
BibTeX
Ref. manager
Endnote
DataCite
NLM
DC