Curvature fluctuations and Lyapunov exponent at Melting

Abstract

We calculate the maximal Lyapunov exponent in constant-energy molecular dynamics simulations at the melting transition for finite clusters of 6 to 13 particles (model rare-gas and metallic systems) as well as for bulk rare-gas solid. For clusters, the Lyapunov exponent generally varies linearly with the total energy, but the slope changes sharply at the melting transition. In the bulk system, melting corresponds to a jump in the Lyapunov exponent, and this corresponds to a singularity in the variance of the curvature of the potential energy surface. In these systems there are two mechanisms of chaos -- local instability and parametric instability. We calculate the contribution of the parametric instability towards the chaoticity of these systems using a recently proposed formalism. The contribution of parametric instability is a continuous function of energy in small clusters but not in the bulk where the melting corresponds to a decrease in this quantity. This implies that the melting in small clusters does not lead to enhanced local instability.