{"id":{"repo_id":"uiuc","oai_identifier":"oai:www.ideals.illinois.edu:2142/81958"},"canonical_url":"https://search.dev.ndltd.org/etd/uiuc/oai:www.ideals.illinois.edu:2142/81958","repository":{"repo_id":"uiuc","name":"University of Illinois - Urbana-Champaign","base_url":"https://www.ideals.illinois.edu/oai-pmh"},"display":{"title":"Longer Time Steps for Molecular Dynamics","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.\"","abstract_html":""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."","abstract_has_math":false,"creators":["Izaguirre Paz, Jesus Antonio"],"institution":"University of Illinois at Urbana-Champaign","degree_name":"Ph.D.","degree_level":"Dissertation","degree_discipline":"Computer Science","degree_department":null,"school":null,"contributors":["Skeel, Robert D."],"advisors":[],"committee_chairs":[],"committee_members":[],"year":2015,"date_issued":"2015-09-25T20:21:11Z","date_published":"2015-09-25T20:21:11Z","updated_at":"2026-07-22T22:26:17Z","subjects":["Applied Mechanics"],"languages":["eng"],"rights":[],"rights_urls":[],"identifier_entries":[{"key":"dc:identifier","label":"Identifier","values":["(MiAaPQ)AAI9953052"],"render_values":[{"text":"(MiAaPQ)AAI9953052","href":null,"code":true}]}]},"links":{"outbound_url":"http://hdl.handle.net/2142/81958","outbound_label":"Handle","outbound_source":"dc:identifier"},"metadata_groups":[{"id":"people","label":"People","entries":[{"key":"dc:contributor","label":"Contributor","values":["Skeel, Robert D."]},{"key":"dc:creator","label":"Author","values":["Izaguirre Paz, Jesus Antonio"]}]},{"id":"academic_context","label":"Academic Context","entries":[{"key":"dc:date","label":"Dc Date","values":["2015-09-25T20:21:11Z","10000-01-01","1999"]},{"key":"dc:type","label":"Dc Type","values":["text"]},{"key":"thesis:degree_discipline","label":"Discipline","values":["Computer Science"]},{"key":"thesis:degree_level","label":"Degree Level","values":["Dissertation"]},{"key":"thesis:degree_name","label":"Degree Name","values":["Ph.D."]},{"key":"thesis:institution_name","label":"Thesis Institution Name","values":["University of Illinois at Urbana-Champaign"]}]},{"id":"subjects_keywords","label":"Subjects and Keywords","entries":[{"key":"dc:subject","label":"Dc Subject","values":["Applied Mechanics"]}]},{"id":"language_rights","label":"Language and Rights","entries":[{"key":"dc:language","label":"Dc Language","values":["eng"]}]},{"id":"identifiers","label":"Identifiers","entries":[{"key":"dc:identifier","label":"Identifier","values":["http://hdl.handle.net/2142/81958","(MiAaPQ)AAI9953052"]}]},{"id":"additional","label":"Additional Metadata","entries":[{"key":"dc:description","label":"Description","values":["\"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.\"","Made available in DSpace on 2015-09-25T20:21:11Z (GMT). No. of bitstreams: 2 license.txt: 4848 bytes, checksum: 96035ab3f5e1c23cc7138a224ce498bd (MD5) 9953052.pdf: 4792026 bytes, checksum: 49b1461ab53ee521b8597eab3b6119eb (MD5) Previous issue date: 1999","Embargo set by: Seth Robbins for item 83239 Lift date: Forever Reason: Restricted to the U of I community idenfinitely during batch ingest of legacy ETDs","Restricted to the U of I community idenfinitely during batch ingest of legacy ETDs","U of I Only","115 p.","Thesis (Ph.D.)--University of Illinois at Urbana-Champaign, 1999."]},{"key":"dc:title","label":"Title","values":["Longer Time Steps for Molecular Dynamics"]}]}],"canonical_facts":{"dc:contributor":["Skeel, Robert D."],"dc:creator":["Izaguirre Paz, Jesus Antonio"],"dc:date":["2015-09-25T20:21:11Z","10000-01-01","1999"],"dc:description":["\"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.\"","Made available in DSpace on 2015-09-25T20:21:11Z (GMT). No. of bitstreams: 2 license.txt: 4848 bytes, checksum: 96035ab3f5e1c23cc7138a224ce498bd (MD5) 9953052.pdf: 4792026 bytes, checksum: 49b1461ab53ee521b8597eab3b6119eb (MD5) Previous issue date: 1999","Embargo set by: Seth Robbins for item 83239 Lift date: Forever Reason: Restricted to the U of I community idenfinitely during batch ingest of legacy ETDs","Restricted to the U of I community idenfinitely during batch ingest of legacy ETDs","U of I Only","115 p.","Thesis (Ph.D.)--University of Illinois at Urbana-Champaign, 1999."],"dc:identifier":["http://hdl.handle.net/2142/81958","(MiAaPQ)AAI9953052"],"dc:language":["eng"],"dc:subject":["Applied Mechanics"],"dc:title":["Longer Time Steps for Molecular Dynamics"],"dc:type":["text"],"thesis:degree_discipline":["Computer Science"],"thesis:degree_level":["Dissertation"],"thesis:degree_name":["Ph.D."],"thesis:institution_name":["University of Illinois at Urbana-Champaign"]},"updated_at":"2026-07-22T22:26:17Z"}