ABSTRACT
Integral equations of uniform fluids have been considered unable to predict any characteristic feature of the fluid-solid phase transition, including the shoulder that arises in the second peak of the fluid-phase radial distribution function, RDF, of hard-core systems obtained by computer simulations, at fluid densities very close to the structural two-step phase transition. This reasoning is based on the results of traditional integral approximations, like Percus-Yevick, PY, which does not show such a shoulder in hard-core systems, neither in two nor three dimensions. In this work, we present results of three Ansätze, based on the PY theory, that were proposed to remedy the lack of PY analytical solutions in two dimensions. This comparative study shows that one of those Ansätze does develop a shoulder in the second peak of the RDF at densities very close to the phase transition, qualitatively describing this feature. Since the shoulder grows into a peak at still higher densities, this integral equation approach predicts the appearance of an orientational order characteristic of the hexatic phase in a continuous fluid-hexatic phase transition.
ABSTRACT
Canonical Monte Carlo (NVT-MC) simulations were performed to obtain surface tension and coexistence densities at the liquid-vapor interface of one-site associating Lennard-Jones and hard-core Yukawa fluids, as functions of association strength and temperature. The method to obtain the components of the pressure tensor from NVT-MC simulations was validated by comparing the equation of state of the associative hard sphere system with that coming from isothermal-isobaric Monte Carlo simulations. Surface tension of the associative Lennard-Jones fluid determined from NVT-MC is compared with previously reported results obtained by molecular dynamics simulations of a pseudomixture model of monomers and dimers. A good agreement was found between both methods. Values of surface tension of associative hard-core Yukawa fluids are presented here for the first time.