Thermodynamic properties of triangle-well fluids in two dimensions: MC and MD simulations.

With the aim of providing complementary data of the thermodynamics properties of the triangular well potential, the vapor/liquid phase diagrams for such potential with different interaction ranges were calculated in two dimensions by Monte Carlo and molecular dynamics simulations; also, the vapor/liquid interfacial tension was calculated. As reported for other interaction potentials, it was observed that the reduction of the dimensionality makes the phase diagram to shrink. Finally, with the aid of reported data for the same potential in three dimensions, it was observed that this potential does not follow the principle of corresponding states.

Constant-force approach to discontinuous potentials.

Aiming to approach the thermodynamical properties of hard-core systems by standard molecular dynamics simulation, we propose setting a repulsive constant-force for overlapping particles. That is, the discontinuity of the pair potential is replaced by a linear function with a large negative slope. Hence, the core-core repulsion, usually modeled with a power function of distance, yields a large force as soon as the cores slightly overlap. This leads to a quasi-hardcore behavior. The idea is tested for a triangle potential of short range. The results obtained by replica exchange molecular dynamics for several repulsive forces are contrasted with the ones obtained for the discontinuous potential and by means of replica exchange Monte Carlo. We found remarkable agreements for the vapor-liquid coexistence densities as well as for the surface tension.

Common behavior of the critical properties of the 2D and 3D square-well fluids.

We have analyzed the behavior of the critical properties and second virial coefficient of the square well fluids in two (2D) and three dimensions (3D) as a function of the interaction range. In both systems, the critical density shows an oscillating-like behavior as the interaction range increases. The second virial coefficient evaluated at the critical temperature as a function of the interaction range shows a general behavior for both cases, and quite surprisingly, there is a minimum of this parameter, for the 2D and 3D fluids, located approximately at the same interaction range. These findings are discussed in terms of the structure of the fluids, via the analysis of the radial distribution function evaluated at the critical point.

Structural properties of fluids interacting via piece-wise constant potentials with a hard core.

The structural properties of fluids whose molecules interact via potentials with a hard core plus two piece-wise constant sections of different widths and heights are presented. These follow from the more general development previously introduced for potentials with a hard core plus n piece-wise constant sections [A. Santos, S. B. Yuste, and M. Lopez de Haro, Condens. Matter Phys. 15, 23602 (2012)] in which use was made of a semi-analytic rational-function approximation method. The results of illustrative cases comprising eight different combinations of wells and shoulders are compared both with simulation data and with those that follow from the numerical solution of the Percus-Yevick and hypernetted-chain integral equations. It is found that the rational-function approximation generally predicts a more accurate radial distribution function than the Percus-Yevick theory and is comparable or even superior to the hypernetted-chain theory. This superiority over both integral equation theories is lost, however, at high densities, especially as the widths of the wells and/or the barriers increase.

Interfacial and coexistence properties of soft spheres with a short-range attractive Yukawa fluid: molecular dynamics simulations.

Molecular dynamics simulations have been carried out to obtain the interfacial and coexistence properties of soft-sphere attractive Yukawa (SAY) fluids with short attraction range, κ = 10, 9, 8, 7, 6, and 5. All our simulation results are new. These data are also compared with the recently reported results in the literature of hard-core attractive Yukawa (HAY) fluids. We show that the interfacial and coexistence properties of both potentials are different. For the surveyed systems, here we show that all coexistence curves collapse into a master curve when we rescale with their respective critical points and the surface tension curves form a single master curve when we plot γ* vs. T/T(c).

Liquid/vapor coexistence and surface tension of the Sutherland fluid with a variable range of interaction: computer simulation and perturbation theory studies.

The liquid-vapor phase diagram and surface tension for the Sutherland fluids with a variable range of interaction have been determined by canonical Monte Carlo simulations and compared with the augmented van der Waals analytic theory. The theory, based on a short-range Yukawa reference, performs well for medium-range models but its accuracy deteriorates with the shortening range due to deteriorating accuracy of the underlined analytic mean spherical approximation solution for the reference Yukawa potential. The simulation results are also analyzed from the point of the extended principle of corresponding states and it is shown that the surface tension and also, to a certain degree of accuracy, the vapor-liquid coexistence curve satisfy the principle.

Modeling of colloid agglomeration inhibition inside a slitlike pore.

An extended particle agglomeration control model and Monte Carlo simulation methodology were used to describe the behavior of the colloidal aggregation process in presence of inhibitor as a confined fluid. In this work results about the mean cluster size, Z, with respect to different variables, temperature, inhibitor concentration and pore size, are presented and showed that Z strongly depends on the slitlike pore size. In addition, a temperature interval where the heating of the system leads to the particle clustering was found.

Thermodynamic properties of short-range attractive Yukawa fluid: simulation and theory.

Coexistence properties of the hard-core attractive Yukawa potential with inverse-range parameter kappa=9, 10, 12, and 15 are calculated by applying canonical Monte Carlo simulation. As previously shown for longer ranges, we show that also for the ranges considered here the coexistence curves scaled by the critical density and temperature obey the law of corresponding states, and that a linear relationship between the critical density and the reciprocal of the critical temperature holds. The simulation results are compared to the predictions of the self-consistent Ornstein-Zernike approximation, and a good agreement is found for both the critical points and the coexistence curves, although some slight discrepancies are present.

Phase diagrams of model C60 and C70 fullerenes from short-range attractive potentials.

We report a computer-simulation study of six model fluids interacting through short-range attractive potentials in order to calculate the vapor-liquid (VL) diagrams using canonical Monte Carlo simulation. It is found that the binodal curves of these systems correctly reproduce those reported in the literature for C(60) and C(70) Girifalco potentials. Besides, we found that all coexistence curves collapse into a master curve when we rescale with their respective critical points.

On the corresponding states law of the Yukawa fluid.

We have analyzed the currently available simulation results as well as performed some additional Monte Carlo simulation for the hard-core attractive Yukawa fluid in order to study its corresponding state behavior. We show that the values of reduced surface tension map onto the master curve and a universal equation of state can be obtained in the wide range of the attractive Yukawa tail length after a certain rescaling of the number density. Some comparisons with other nonconformal potentials are presented and discussed.

Phase diagram and surface tension of the hard-core attractive Yukawa model of variable range: Monte Carlo simulations.

The liquid-vapor phase diagram and surface tension for hard-core Yukawa potential with 4

Oscillatory surface tension due to finite-size effects.

The simulation results of surface tension at the liquid-vapor interface are presented for fluids interacting with Lennard Jones and square-well potentials. From the simulation of liquids we have reported [M. González-Melchor et al., J. Chem. Phys. 122, 4503 (2005)] that the components of pressure tensor in parallelepiped boxes are not the same when periodic boundary conditions and small transversal areas are used. This fact creates an artificial oscillatory stress anisotropy in the system with even negative values. By doing direct simulations of interfaces we show in this work that surface tension has also an oscillatory decay at small surface areas; this behavior is opposite to the monotonic decay reported previously for the Lennard Jones fluid. It is shown that for small surface areas, the surface tension of the square-well potential artificially takes negative values and even increases with temperature. The calculated surface tension using a direct simulation of interfaces might have two contributions: one from finite-size effects of interfacial areas due to box geometry and another from the interface. Thus, it is difficult to evaluate the true surface tension of an interface when small surface areas are used. Care has to be taken to use the direct simulation method of interfaces to evaluate the predicted surface tension as a function of interfacial area from capillary-wave theory. The oscillations of surface tension decay faster at temperatures close to the critical point. It is also discussed that a surface area does not show any important effect on coexisting densities, making this method reliable to calculate bulk coexisting properties using small systems.

Stress anisotropy induced by periodic boundary conditions.

Finite size effects due to periodic boundary conditions are investigated using computer simulations in the canonical ensemble. We study liquids with densities corresponding to typical liquid coexistence densities, and temperatures between the triple and critical points. The components of the pressure tensor are computed in order to analyze the finite size effects arising from the size and geometry of the simulation box. Two different box geometries are considered: cubic and parallelepiped. As expected the pressure tensor is isotropic in cubic boxes, but it becomes anisotropic for small noncubic boxes. We argue this is the origin of the anomalous behavior observed recently in the computation of the surface tension of liquid-vapor interfaces. Otherwise, we find that the bulk pressure is sensitive to the box geometry when small simulation boxes are considered. These observations are general and independent of the model liquid considered. We report results for liquids interacting through short range forces, square well and Lennard-Jones, and also long range Coulombic interactions. The effect that small surface areas have on the surface tension is discussed, and some preliminary results at the liquid vapor-interface for the square well potential are given.

Liquid-vapor interface of square-well fluids of variable interaction range.

Properties of the liquid-vapor interface of square-well fluids with ranges of interaction lambda=1.5, 2.0, and 3.0 are obtained by Monte Carlo simulations and from square-gradient theories that combine the Carnahan-Starling equation of state for hard spheres with the second and third virial coefficients. The predicted surface tensions show good agreement with the simulation results for lambda=2 and for lambda=3 in a temperature range reasonably close to the critical point, 0.8

Surface tension of associating fluids by Monte Carlo simulations.

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.