Improved algorithm for geometry optimisation using damped molecular dynamics

Research output: Contribution to journalArticlepeer-review


We present a new algorithm for optimising the configuration of a system of particles using damped molecular dynamics. This new algorithm has a much faster rate of convergence to the ground-state structure than other MD-based schemes. We also show how it is possible to use this algorithm to select dynamically the optimum time step or to precondition the masses for the modified dynamics, so as to produce additional significant efficiency gains especially as convergence is approached. The performance of the new algorithm is illustrated with several model systems and also by relaxing a defect in bulk silicon. The defect calculation was also relaxed using a standard BFGS algorithm and the performance was found to be comparable. (C) 2003 Elsevier Science B.V. All rights reserved.

Original languageEnglish
Pages (from-to)130-146
Number of pages16
JournalJournal of Computational Physics
Issue number1
Publication statusPublished - 10 Oct 2003


  • 65B99 numerical analysis
  • acceleration of convergence
  • 70-08 mechanics of particles and systems
  • computational methods

Cite this