The behavior of thermodynamic response functions and the thermodynamic scalar curvature in the supercritical region have been studied for a Lennard-Jones fluid based on a revised modified Benedict-Webb-Rubin equation of state. Response function extrema are sometimes used to estimate the Widom line, which is characterized by the maxima of the correlation lengths. We calculated the Widom line for the Lennard-Jones fluid without using any response function extrema. Since the volume of the correlation length is proportional to the Riemannian thermodynamic scalar curvature, the locus of the Widom line follows the slope of maximum curvature. We show that the slope of the Widom line follows the slope of the isobaric heat capacity maximum only in the close vicinity of the critical point and that, therefore, the use of response function extrema in this context is problematic. Furthermore, we constructed the vapor-liquid coexistence line for the Lennard-Jones fluid using the fact that the correlation length, and therefore the thermodynamic scalar curvature, must be equal in the two coexisting phases. We compared the resulting phase envelope with those from simulation data where multiple histogram reweighting was used and found striking agreement between the two methods.

We report on a molecular dynamics (MD) study of the connectivity of hydrogen bond networks in liquid water, focusing primarily on the microscopic distribution functions giving the weight fraction of molecules belonging to a ‘‘net’’ of M molecules (M=1,2,3,...). The MD data compare favorably—using no adjustable parameters—with predictions of random bond percolation theory. We also study the connectivity of those molecules with four intact hydrogen bonds, and compare the corresponding distribution functions with correlated-site percolation theory. Our analysis supports the proposal that when looking at the bond connectivity, water appears as a macroscopic space-filling network—as expected from continuum models of water. When looking at the correlated site percolation problem defined by the four-bonded molecules, water appears as a myriad of tiny ramified low-density patches, somewhat reminiscent of mixture theories and cluster models. In Appendix A, we find a strong correlation between the number of molecules within a sphere of radius rc around a given molecule and the total interaction energy of that molecule with its neighbors residing in that sphere; for most choices of rc, the energy becomes less negative when more molecules are in the sphere, in contrast to the behavior of a normal fluid. This result supports the finding of Geiger and Stanley that regions of high bond connectivity are correlated with regions of low density. In Appendix B we describe in detail how we adapt conventional percolation theory to the calculation of cluster size distribution functions for hydrogen bond networks in water.

No abstract

The thermodynamic curvature scalar R for the Lennard-Jones system is evaluated in phase space, including vapor, liquid, and solid state. We paid special attention to the investigation of R along vapor-liquid, liquid-solid, and vapor-solid equilibria. Because R is a measure of interaction strength, we traced out the line R=0 dividing the phase space into regions with effectively attractive (R<0) or repulsive (R>0) interactions. Furthermore, we analyzed the dependence of R on the strength of attraction applying a perturbation ansatz proposed by Weeks-Chandler-Anderson. Our results show clearly a transition from R>0 (for poorly repulsive interaction) to R<0 when loading attraction in the intermolecular potential.