A simple methodology for analyzing association effects on response functions via Monte Carlo simulations.

A simple methodology was developed to analyze association effects on the thermodynamic response functions for a pure self-associated fluid via Monte Carlo simulations. The procedure essentially involves expressing the residual energy and volume of the fluid in terms of these properties for two hypothetical fluids consisting of monomers and associated molecules, respectively. This allows the thermodynamic response functions to be expressed in a perturbative form as a combination of the values for the property in the monomeric fluid and the contribution of association (the perturbative term). The proposed methodology was used to determine both contributions to the isobaric heat capacity and to the temperature and pressure derivatives of the volume for OPLS methanol along the 50 MPa isobar from 220 to 1500 K. Based on the results, both terms exert a substantial influence on the isobaric heat capacity; by contrast, the association term for the volumetric properties is negligible. These results are consistent with those of a previous work involving simulations with the same model under identical thermodynamic conditions but a different approach. They are also compared with others previously reported in context. Moreover, a comprehensive study of the different types of clusters present in the fluid was performed and the results were related to thermodynamic properties. A strong correlation between the heat capacity of the monomeric fluid and this structural analysis was found.