Accelerating convergence of molecular dynamics-based structural relaxation

Research output: Contribution to journalJournal articleResearchpeer-review

Abstract

We describe strategies to accelerate the terminal stage of molecular dynamics (MD)based relaxation algorithms, where a large fraction of the computational resources are used. First, we analyze the qualitative and quantitative behavior of the QuickMin family of MD relaxation algorithms and explore the influence of spectral properties and dimensionality of the molecular system on the algorithm efficiency. We test two algorithms, the MinMax and Lanczos, for spectral estimation from an MD trajectory, and use this to derive a practical scheme of time step adaptation in MD relaxation algorithms to improve efficiency. We also discuss the implementation aspects. Secondly, we explore the final state refinement acceleration by a combination with the conjugate gradient technique, where the key ingredient is an implicit corrector step. Finally, we test the feasibility of passive Hessian matrix accumulation from an MD trajectory, as another route for final phase acceleration. Our suggestions may be implemented within most MD quench implementations with a few, straightforward lines of code, thus maintaining the appealing simplicity of the MD quench algorithms. In this paper, we also bridge the conceptual gap between the MD quench algorithms inspired from physics and the mathematically rooted line search algorithms.
Original languageEnglish
JournalInternational Journal of Modern Physics C
Volume16
Issue number2
Pages (from-to)193-223
ISSN0129-1831
Publication statusPublished - 2005

Keywords

  • Lanczos algorithms
  • time step adaptation
  • spectral estimation
  • structural optimization
  • molecular dynamics relaxation

Fingerprint

Dive into the research topics of 'Accelerating convergence of molecular dynamics-based structural relaxation'. Together they form a unique fingerprint.

Cite this