Journal of Chemical Physics, Vol.101, No.9, 7934-7938, 1994
Multibody Distribution Function Contributions to the Heat-Capacity for the Truncated Lennard-Jones Fluid
We have performed a series of standard NVT Monte Carlo simulations in which we have calculated the configurational heat capacity, C-v, of a truncated Lennard-Jones model fluid by the fluctuation theorem. As we found in a previous investigation, the heat capacity exhibited two extrema as a function of density on sufficiently low temperature isotherms. During the course of these runs, we have decomposed the heat capacity into the sum of quantities, C-v((2)), C-v((3)), and C-v((4)), which are averages over two, three, and four body distribution functions. We analyze combinations of the three contributions to gain insight into the nature of the extrema. In particular, we find that the sum C-v((2))+C-v((3)) has a maximum and a minimum at all temperatures studied. The magnitude and slope of C-v((2))+C-v((3)) compared with C-v((4)) determines whether or not C-v itself has maxima and minima.
Keywords:MONTE-CARLO SIMULATION