Thermodynamic response functions of fluids: a microscopic approach based on NpT Monte Carlo

On the basis of NpT Monte Carlo simulations, a detailed analysis on the microscopic origins of some specific features of thermodynamic response functions of fluids is performed. Specifically, the residual isobaric heat capacity C(p) (res), the isobaric thermal expansivity alpha(p), and the isothermal compressibility kappa(T) for Lennard-Jones methane and optimized potential for liquid simulations (OPLS) methanol have been determined via standard techniques. For the former, data along the liquid, gas, and supercritical regions are presented, while a wide temperature range at a single supercritical pressure is covered for the latter. They have been obtained by computing the various pairwise fluctuations contributing to each property. Attention is mainly focused on isothermal and isobaric maxima found for both C(p) (res) and alpha(p), which have been rationalized at a molecular level using qualitative arguments. It is encountered that maxima emerge as a natural consequence of the destruction of fluid structure as temperature is increased or as pressure is decreased. The results for Lennard-Jones methane reveal the competition of energetic and volumetric effects, while those for OPLS methanol evidence that hydrogen-bonding is dominant as energetic effects are concerned. Further discussion on previous results and alternative approaches using equations of state as well as on closely related topics such as "maxima and critical phenomena" is included.