New approaches to potential energy minimization and molecular dynamics algorithms

Research output: Contribution to journalArticle

Abstract

We describe two new algorithms for macromolecular simulations: a truncated Newton method for potential energy minimization and an implicit integration scheme for molecular dynamics (MD). The truncated Newton algorithm is specifically adapted for large-scale potential energy functions. It uses analytic second derivatives and exploits the separability structure of the Hessian into bonded and nonbonded terms. The method is rapidly convergent (with a quadratic convergence rate) and allows variations for avoiding analytic computation of the nonbonded Hessian terms. The MD algorithm combines the implicit Euler scheme for integration with the Langevin dynamics formulation. The implicit scheme permits a wide range of time steps without loss of numerical stability. In turn, it requires that a nonlinear system be solved at every step. We accomplish this task by formulating a related minimization problem-not to be confused with minimization of the potential energy-that can be solved rapidly with the truncated Newton method. Additionally, the MD scheme permits the introduction of a "cutoff" frequency (ωc) which, in particular, can be used to mimic the quantum-mechanical discrimination among activity of the various vibrational modes.

Original languageEnglish (US)
Pages (from-to)251-260
Number of pages10
JournalComputers and Chemistry
Volume15
Issue number3
DOIs
StatePublished - 1991

Fingerprint

Dynamic Algorithms
Energy Minimization
Molecular Dynamics Simulation
Truncated Newton Method
Potential energy
Molecular Dynamics
Molecular dynamics
Implicit Scheme
Newton-Raphson method
Potential energy functions
Langevin Dynamics
Euler Scheme
Quadratic Convergence
Convergence of numerical methods
Cutoff frequency
Numerical Stability
Separability
Term
Second derivative
Potential Function

ASJC Scopus subject areas

  • Biotechnology
  • Chemical Engineering(all)
  • Applied Microbiology and Biotechnology

Cite this

New approaches to potential energy minimization and molecular dynamics algorithms. / Schlick, Tamar.

In: Computers and Chemistry, Vol. 15, No. 3, 1991, p. 251-260.

Research output: Contribution to journalArticle

@article{61b5a15f2dc7417ea90be5601550fcc8,
title = "New approaches to potential energy minimization and molecular dynamics algorithms",
abstract = "We describe two new algorithms for macromolecular simulations: a truncated Newton method for potential energy minimization and an implicit integration scheme for molecular dynamics (MD). The truncated Newton algorithm is specifically adapted for large-scale potential energy functions. It uses analytic second derivatives and exploits the separability structure of the Hessian into bonded and nonbonded terms. The method is rapidly convergent (with a quadratic convergence rate) and allows variations for avoiding analytic computation of the nonbonded Hessian terms. The MD algorithm combines the implicit Euler scheme for integration with the Langevin dynamics formulation. The implicit scheme permits a wide range of time steps without loss of numerical stability. In turn, it requires that a nonlinear system be solved at every step. We accomplish this task by formulating a related minimization problem-not to be confused with minimization of the potential energy-that can be solved rapidly with the truncated Newton method. Additionally, the MD scheme permits the introduction of a {"}cutoff{"} frequency (ωc) which, in particular, can be used to mimic the quantum-mechanical discrimination among activity of the various vibrational modes.",
author = "Tamar Schlick",
year = "1991",
doi = "10.1016/0097-8485(91)80014-D",
language = "English (US)",
volume = "15",
pages = "251--260",
journal = "Computational Biology and Chemistry",
issn = "1476-9271",
publisher = "Elsevier Limited",
number = "3",

}

TY - JOUR

T1 - New approaches to potential energy minimization and molecular dynamics algorithms

AU - Schlick, Tamar

PY - 1991

Y1 - 1991

N2 - We describe two new algorithms for macromolecular simulations: a truncated Newton method for potential energy minimization and an implicit integration scheme for molecular dynamics (MD). The truncated Newton algorithm is specifically adapted for large-scale potential energy functions. It uses analytic second derivatives and exploits the separability structure of the Hessian into bonded and nonbonded terms. The method is rapidly convergent (with a quadratic convergence rate) and allows variations for avoiding analytic computation of the nonbonded Hessian terms. The MD algorithm combines the implicit Euler scheme for integration with the Langevin dynamics formulation. The implicit scheme permits a wide range of time steps without loss of numerical stability. In turn, it requires that a nonlinear system be solved at every step. We accomplish this task by formulating a related minimization problem-not to be confused with minimization of the potential energy-that can be solved rapidly with the truncated Newton method. Additionally, the MD scheme permits the introduction of a "cutoff" frequency (ωc) which, in particular, can be used to mimic the quantum-mechanical discrimination among activity of the various vibrational modes.

AB - We describe two new algorithms for macromolecular simulations: a truncated Newton method for potential energy minimization and an implicit integration scheme for molecular dynamics (MD). The truncated Newton algorithm is specifically adapted for large-scale potential energy functions. It uses analytic second derivatives and exploits the separability structure of the Hessian into bonded and nonbonded terms. The method is rapidly convergent (with a quadratic convergence rate) and allows variations for avoiding analytic computation of the nonbonded Hessian terms. The MD algorithm combines the implicit Euler scheme for integration with the Langevin dynamics formulation. The implicit scheme permits a wide range of time steps without loss of numerical stability. In turn, it requires that a nonlinear system be solved at every step. We accomplish this task by formulating a related minimization problem-not to be confused with minimization of the potential energy-that can be solved rapidly with the truncated Newton method. Additionally, the MD scheme permits the introduction of a "cutoff" frequency (ωc) which, in particular, can be used to mimic the quantum-mechanical discrimination among activity of the various vibrational modes.

UR - http://www.scopus.com/inward/record.url?scp=0011682977&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0011682977&partnerID=8YFLogxK

U2 - 10.1016/0097-8485(91)80014-D

DO - 10.1016/0097-8485(91)80014-D

M3 - Article

VL - 15

SP - 251

EP - 260

JO - Computational Biology and Chemistry

JF - Computational Biology and Chemistry

SN - 1476-9271

IS - 3

ER -