{"id":{"repo_id":"uiuc","oai_identifier":"oai:www.ideals.illinois.edu:2142/81667"},"canonical_url":"https://search.dev.ndltd.org/etd/uiuc/oai:www.ideals.illinois.edu:2142/81667","repository":{"repo_id":"uiuc","name":"University of Illinois - Urbana-Champaign","base_url":"https://www.ideals.illinois.edu/oai-pmh"},"display":{"title":"Fast Polarizable Force Field Computation in Biomolecular Simulations","abstract":"Polarizable force fields are considered to be the single most significant development in the next-generation force fields used in biomolecular simulations. The self-consistent computation of induced atomic dipoles in a polarizable force field is expensive due to the cost of solving a large dense linear system at each timestep in molecular dynamics simulations. Methods are developed that reduce the cost of computing the electrostatic energy and force of a polarizable model from about 7.5 times the cost of computing those of a non-polarizable model to less than twice the cost. The reduction is achieved by an efficient implementation of the particle-mesh Ewald method, an accurate and robust predictor based on least squares fitting, and two non-stationary iterative methods whose fast convergence is empowered by a simple preconditioner. Furthermore, with these methods, we show that the self-consistent approach with a larger timestep is faster than the extended Lagrangian approach. The use of dipole moments from previous timesteps to calculate an accurate initial guess for iterative methods leads to an energy drift and compromises the volume-preserving property of the integration. Iterative methods with zero initial guess do not lead to perceptible energy drift if a reasonably strict convergence criterion for the iteration is imposed and the numerical integrator is volume-preserving. The approximate solution computed by an iterative method ruins the symplectic property of the integrator. To address this problem, a non-iterative method has been developed based on an approximation to the electrostatic potential energy and has been efficiently implemented. The method preserves the symplecticness of the integrator and is suitable for long time simulations. The research will help polarizable force fields modeling and computation to become a routine part of molecular dynamics simulations for biomolecular systems.","abstract_html":"Polarizable force fields are considered to be the single most significant development in the next-generation force fields used in biomolecular simulations. The self-consistent computation of induced atomic dipoles in a polarizable force field is expensive due to the cost of solving a large dense linear system at each timestep in molecular dynamics simulations. Methods are developed that reduce the cost of computing the electrostatic energy and force of a polarizable model from about 7.5 times the cost of computing those of a non-polarizable model to less than twice the cost. The reduction is achieved by an efficient implementation of the particle-mesh Ewald method, an accurate and robust predictor based on least squares fitting, and two non-stationary iterative methods whose fast convergence is empowered by a simple preconditioner. Furthermore, with these methods, we show that the self-consistent approach with a larger timestep is faster than the extended Lagrangian approach. The use of dipole moments from previous timesteps to calculate an accurate initial guess for iterative methods leads to an energy drift and compromises the volume-preserving property of the integration. Iterative methods with zero initial guess do not lead to perceptible energy drift if a reasonably strict convergence criterion for the iteration is imposed and the numerical integrator is volume-preserving. The approximate solution computed by an iterative method ruins the symplectic property of the integrator. To address this problem, a non-iterative method has been developed based on an approximation to the electrostatic potential energy and has been efficiently implemented. The method preserves the symplecticness of the integrator and is suitable for long time simulations. The research will help polarizable force fields modeling and computation to become a routine part of molecular dynamics simulations for biomolecular systems.","abstract_has_math":false,"creators":["Wang, Wei"],"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:19:56Z","date_published":"2015-09-25T20:19:56Z","updated_at":"2026-07-22T22:26:16Z","subjects":["Biophysics, General"],"languages":["eng"],"rights":[],"rights_urls":[],"identifier_entries":[{"key":"dc:identifier","label":"Identifier","values":["(MiAaPQ)AAI3182416"],"render_values":[{"text":"(MiAaPQ)AAI3182416","href":null,"code":true}]}]},"links":{"outbound_url":"http://hdl.handle.net/2142/81667","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":["Wang, Wei"]}]},{"id":"academic_context","label":"Academic Context","entries":[{"key":"dc:date","label":"Dc Date","values":["2015-09-25T20:19:56Z","10000-01-01","2005"]},{"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":["Biophysics, General"]}]},{"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/81667","(MiAaPQ)AAI3182416"]}]},{"id":"additional","label":"Additional Metadata","entries":[{"key":"dc:description","label":"Description","values":["Polarizable force fields are considered to be the single most significant development in the next-generation force fields used in biomolecular simulations. The self-consistent computation of induced atomic dipoles in a polarizable force field is expensive due to the cost of solving a large dense linear system at each timestep in molecular dynamics simulations. Methods are developed that reduce the cost of computing the electrostatic energy and force of a polarizable model from about 7.5 times the cost of computing those of a non-polarizable model to less than twice the cost. The reduction is achieved by an efficient implementation of the particle-mesh Ewald method, an accurate and robust predictor based on least squares fitting, and two non-stationary iterative methods whose fast convergence is empowered by a simple preconditioner. Furthermore, with these methods, we show that the self-consistent approach with a larger timestep is faster than the extended Lagrangian approach. The use of dipole moments from previous timesteps to calculate an accurate initial guess for iterative methods leads to an energy drift and compromises the volume-preserving property of the integration. Iterative methods with zero initial guess do not lead to perceptible energy drift if a reasonably strict convergence criterion for the iteration is imposed and the numerical integrator is volume-preserving. The approximate solution computed by an iterative method ruins the symplectic property of the integrator. To address this problem, a non-iterative method has been developed based on an approximation to the electrostatic potential energy and has been efficiently implemented. The method preserves the symplecticness of the integrator and is suitable for long time simulations. The research will help polarizable force fields modeling and computation to become a routine part of molecular dynamics simulations for biomolecular systems.","Made available in DSpace on 2015-09-25T20:19:56Z (GMT). No. of bitstreams: 2 license.txt: 4848 bytes, checksum: 96035ab3f5e1c23cc7138a224ce498bd (MD5) 3182416.pdf: 3487004 bytes, checksum: 4bbae74cb69d2e5d733fe7d52573bf18 (MD5) Previous issue date: 2005","Embargo set by: Seth Robbins for item 82948 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","93 p.","Thesis (Ph.D.)--University of Illinois at Urbana-Champaign, 2005."]},{"key":"dc:title","label":"Title","values":["Fast Polarizable Force Field Computation in Biomolecular Simulations"]}]}],"canonical_facts":{"dc:contributor":["Skeel, Robert D."],"dc:creator":["Wang, Wei"],"dc:date":["2015-09-25T20:19:56Z","10000-01-01","2005"],"dc:description":["Polarizable force fields are considered to be the single most significant development in the next-generation force fields used in biomolecular simulations. The self-consistent computation of induced atomic dipoles in a polarizable force field is expensive due to the cost of solving a large dense linear system at each timestep in molecular dynamics simulations. Methods are developed that reduce the cost of computing the electrostatic energy and force of a polarizable model from about 7.5 times the cost of computing those of a non-polarizable model to less than twice the cost. The reduction is achieved by an efficient implementation of the particle-mesh Ewald method, an accurate and robust predictor based on least squares fitting, and two non-stationary iterative methods whose fast convergence is empowered by a simple preconditioner. Furthermore, with these methods, we show that the self-consistent approach with a larger timestep is faster than the extended Lagrangian approach. The use of dipole moments from previous timesteps to calculate an accurate initial guess for iterative methods leads to an energy drift and compromises the volume-preserving property of the integration. Iterative methods with zero initial guess do not lead to perceptible energy drift if a reasonably strict convergence criterion for the iteration is imposed and the numerical integrator is volume-preserving. The approximate solution computed by an iterative method ruins the symplectic property of the integrator. To address this problem, a non-iterative method has been developed based on an approximation to the electrostatic potential energy and has been efficiently implemented. The method preserves the symplecticness of the integrator and is suitable for long time simulations. The research will help polarizable force fields modeling and computation to become a routine part of molecular dynamics simulations for biomolecular systems.","Made available in DSpace on 2015-09-25T20:19:56Z (GMT). No. of bitstreams: 2 license.txt: 4848 bytes, checksum: 96035ab3f5e1c23cc7138a224ce498bd (MD5) 3182416.pdf: 3487004 bytes, checksum: 4bbae74cb69d2e5d733fe7d52573bf18 (MD5) Previous issue date: 2005","Embargo set by: Seth Robbins for item 82948 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","93 p.","Thesis (Ph.D.)--University of Illinois at Urbana-Champaign, 2005."],"dc:identifier":["http://hdl.handle.net/2142/81667","(MiAaPQ)AAI3182416"],"dc:language":["eng"],"dc:subject":["Biophysics, General"],"dc:title":["Fast Polarizable Force Field Computation in Biomolecular Simulations"],"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:16Z"}