FLUCTUATIONS AND THERMODYNAMIC RESPONSE FUNCTIONS IN A LENNARD-JONES SOLID

Thermodynamic response functions of a nearest-neighbor Lennard-Jones solid-heat capacity, thermal-expansion coefficient, compressibility, and elastic constants-are calculated directly from fluctuations using molecular-dynamics simulations. The algorithm used is the earlier Parrinello-Rahman molecular dynamics modified to take into account symmetry and rotation invariance of the system under investigation. The convergence is very fast and results are in good agreement with existing Monte Carlo and molecular-dynamics results.