Equilibrium thermodynamics and thermodynamic processes in nonlinear systems

Abstract

A general method for the calculation of thermodynamic quasi-equilibrium processes by molecular dynamics (MD) simulation in canonical ensemble is developed. The method is suitable for classical systems with arbitrary interaction potentials. Though this MD method does not allow to calculate directly partition function and entropy, it is possible to calculate necessary partial derivatives which enter into the expressions for the full derivatives dT/dV and d β/dV for adiabitic and isobaric processes. Namely the solutions of these ordinary differential equations define the corresponding thermodynamic processes. The adiabatic process for the 1D Toda lattice is analyzed in details. The usage of the Toda potential allows to perform all analytical calculus up to accurate answers and to compare numerical and analytical results. Exact analytical expressions for the thermodynamics of 1D lattices with few types of nearest neighbor interactions are obtained as a necessary interim solutions. MD-simulation of quasi equilibrium processes in canonical ensemble demands the achievement of thermodynamic equilibrium, thus the thermalization kinetics is briefly discussed.