Inom molekyldynamik bör modeller av metalliska system i allmänhet ha systemets temperatur som en beroende variabel. I synnerhet bör termen i systemets Hamiltonian som representerar potentiell energi utöver det interpartikulära avståndet även vara beroende av temperatur. Detta temperaturberoende gör i allmänhet Hamiltonianen icke-separabel. Konventionella explicita numeriska metoder som är symplektiska då de används på system med separabel Hamiltonian är i allmänhet inte symplektiska då de används i system med icke-separabel Hamiltonian. På grund av detta eftersöks en integrator som behåller symplekticitet då den används i system med en icke-separabel Hamiltonian. En samling integratorer som är symplektiska även då de används på system med icke-separabel Hamiltonian visas prestera lika bra eller bättre än de konventionella Velocity Verlet- och fjärde ordningens Runge Kutta-integratorerna, med nackdelen att de undersökta integratorerna uppvisar numerisk instabilitet då de tillämpas på system där de interpartikulära krafterna beror exponentiellt på inverterade interpartikulära avstånd. Till författarnas kännedom är denna studie den första tillämpningen av de undersökta integratorerna inom molekyldynamik. Resultaten i denna studie ger en fingervisning om att de undersökta integratorerna inte ger en övergripande lösning på problemet att integrera rörelseekvationerna hos ett system med icke-separabel Hamiltonian. Emellertid behövs vidare undersökningar som använder den undersökta samlingen av integratorer i andra system i molekyldynamik än de som undersöks i denna studie för att ge en mer definitiv slutsats.

In molecular dynamics, mathematical models of metallic systems should in general have the temperature of the system as a dependent variable. In particular, the potential energy term of the Hamiltonian function of the interaction model should be dependent on temperature in addition to interparticular distances. This puts the Hamiltonian function on a form which is generally non-separable. Conventional explicit numerical methods which are symplectic when used to integrate the equations of motion of systems with separable Hamiltonians are not in general symplectic when used to integrate the equations of motion of systems with a non-separable Hamiltonian. Hence, an integrator which sustains symplecticity when used in a system with non-separable Hamiltonian is sought. A family of explicit integrators which are symplectic when integrating systems with a non-separable Hamiltonian are shown to exhibit similar or superior performance to the Velocity Verlet and fourth-order Runge-Kutta schemes, albeit with the drawback of numerical instability when used on a system where forces depend exponentially on the inverted interparticular distances. To the knowledge of the authors, this study is the first time this family of integrators is applied in the context of molecular dynamics. The results of this study provide a first indication that a comprehensive solution to the problem of integrating the equations of motion of a system with a non-separable Hamiltonian explicitly and symplectically is not provided by the considered family of integrators. However, further investigations into using this family of integrators in other molecular dynamics systems than those investigated here are needed to provide a more definitive conclusion.