We propose a novel theoretical and practical alternative to the Maupertuis functional in the field of molecular dynamics: the Tonelli functional. Our aim is to adapt this technique to the study of rare events where the initial and the final states of a system are known, and we look for transition paths. We reach it with a rigorous mathematical development of the functional and an efficient numerical algorithm. We couple the Tonelli functional with an exact finite dimensional reduction, and we prove error estimates for its implementation. This is not far from a multiscale approach; indeed the reduction will help the study of the high frequencies of the systems: it can be seen as a magnifying glass for atomic trajectories. We test our techniques first on simple models in order to show the details and the main features of these new tools: an harmonic oscillator, a one dimensional double-well potential, the well-tested Mueller potential, and a short oscillatory trajectory of a 4-atom cluster are studied under different points of view. Then we pass to a more demanding test: the isomerization of a Lennard-Jones cluster of 38 interacting atoms.

Tonelli Principle: Finite Reduction and Fixed Energy Molecular Dynamics Trajectories

CARDIN, FRANCO
2009

Abstract

We propose a novel theoretical and practical alternative to the Maupertuis functional in the field of molecular dynamics: the Tonelli functional. Our aim is to adapt this technique to the study of rare events where the initial and the final states of a system are known, and we look for transition paths. We reach it with a rigorous mathematical development of the functional and an efficient numerical algorithm. We couple the Tonelli functional with an exact finite dimensional reduction, and we prove error estimates for its implementation. This is not far from a multiscale approach; indeed the reduction will help the study of the high frequencies of the systems: it can be seen as a magnifying glass for atomic trajectories. We test our techniques first on simple models in order to show the details and the main features of these new tools: an harmonic oscillator, a one dimensional double-well potential, the well-tested Mueller potential, and a short oscillatory trajectory of a 4-atom cluster are studied under different points of view. Then we pass to a more demanding test: the isomerization of a Lennard-Jones cluster of 38 interacting atoms.
File in questo prodotto:
Non ci sono file associati a questo prodotto.
Pubblicazioni consigliate

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11577/2376935
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus ND
  • ???jsp.display-item.citation.isi??? 2
social impact