Generalized equation of state for square-well potentials of variable range

The physics of empty liquids: from patchy particles to water

Empty liquids represent a wide class of materials whose constituents arrange in a random network through reversible bonds. Many key insights on the physical properties of empty liquids have originated almost independently from the study of colloidal patchy particles on one side, and a large body of theoretical and experimental research on water on the other side. Patchy particles represent a family of coarse-grained potentials that allows for a precise control of both the geometric and the energetic aspects of bonding, while water has arguably the most complex phase diagram of any pure substance, and a puzzling amorphous phase behavior. It was only recently that the exchange of ideas from both fields has made it possible to solve long-standing problems and shed new light on the behavior of empty liquids. Here we highlight the connections between patchy particles and water, focusing on the modelling principles that make an empty liquid behave like water, including the factors that control the appearance of thermodynamic and dynamic anomalies, the possibility of liquid-liquid phase transitions, and the crystallization of open crystalline structures.

Extended law of corresponding states: square-well oblates

The vapour-liquid coexistence collapse in the reduced temperature,T_r=T/T_c, reduced density,ρ_r=ρ/ρ_c, plane is known as a principle of corresponding states, and Noro and Frenkel have extended it for pair potentials of variable range. Here, we provide a theoretical basis supporting this extension, and show that it can also be applied to short-range pair potentials where both repulsive and attractive parts can be anisotropic. We observe that the binodals of oblate hard ellipsoids for a given aspect ratio (κ= 1/3) with varying short-range square-well interactions collapse into a single master curve in theΔB2*-ρ_rplane, whereΔB2*=(B2(T)-B2(Tc))/v0,B₂is the second virial coefficient, andv₀is the volume of the hard body. This finding is confirmed by both REMC simulation and second virial perturbation theory for varying square-well shells, mimicking uniform, equator, and pole attractions. Our simulation results reveal that the extended law of corresponding states is not related to the local structure of the fluid.

Generalized equation of state for fluids: From molecular liquids to colloidal dispersions

In this work, a new parameterization for the Statistical Association Fluid Theory for potentials of Variable Range (SAFT-VR) is coupled to the discrete potential theory to represent the thermodynamic properties of several fluids, ranging from molecular liquids to colloidal-like dispersions. In this way, this version of the SAFT-VR approach can be straightforwardly applied to any kind of either simple or complex fluid. In particular, two interaction potentials, namely, the Lennard-Jones and the hard-core attractive Yukawa potentials, are discretized to study the vapor-liquid equilibrium properties of both molecular and complex liquids, respectively. Our results are assessed with Monte Carlo computer simulations and available and accurate theoretical results based on the self-consistent Ornstein-Zernike approximation.

A Thermodynamic Approach for the Prediction of Oiling Out Boundaries from Solubility Data

Many pharmaceutical molecules, fine chemicals, and proteins exhibit liquid–liquid phase separation (LLPS, also known as oiling out) during solution crystallization. LLPS is of significant concern in crystallization process development, as oiling out can compromise the effectiveness of a crystallization and can lead to operational problems. A comprehensive methodology that allows a process scientist/engineer to characterize the various phase boundaries relevant to oiling out is currently lacking. In this work, we present a modeling framework useful in predicting the binodal, spinodal, and gelation boundaries starting from the solubility data of a solute that is prone to oiling out. We collate the necessary theoretical concepts from the literature and describe a unified approach to model the phase equilibria of solute–solvent systems from first principles. The modeling effort is validated using experimental data reported in the literature for various solute–solvent systems. The predictive methods presented in this work can be easily implemented and help a process engineer establish the design space for a crystallization process that is affected by liquid–liquid phase separation.

Vapor-liquid equilibrium and equation of state of two-dimensional fluids from a discrete perturbation theory

The interest in the description of the properties of fluids of restricted dimensionality is growing for theoretical and practical reasons. In this work, we have firstly developed an analytical expression for the Helmholtz free energy of the two-dimensional square-well fluid in the Barker-Henderson framework. This equation of state is based on an approximate analytical radial distribution function for d-dimensional hard-sphere fluids (1 ≤ d ≤ 3) and is validated against existing and new simulation results. The so-obtained equation of state is implemented in a discrete perturbation theory able to account for general potential shapes. The prototypical Lennard-Jones and Yukawa fluids are tested in its two-dimensional version against available and new simulation data with semiquantitative agreement.

The constant force continuous molecular dynamics for potentials with multiple discontinuities

In this work, we present an extension of the constant force approach [P. Orea and G. Odriozola, J. Chem. Phys. 138, 214105 (2013)] to the case of potentials with multiple discontinuities. To illustrate the method, we selected the square well potential of range λ=1.5 that exhibits two discontinuities. Square well single phase properties, vapor-liquid phase diagram, and surface tension were calculated and compared with available simulation data. Besides, we analyzed the internal energies of a square well plus a square shoulder potential having three discontinuities. For both potentials, a good agreement has been found when compared with results of other simulation techniques (discontinuous molecular dynamics and Monte Carlo methods). This extension can be easily implemented to more general and efficient continuous molecular dynamics packages (HOOMD, GROMACS, NAMD, etc.).

Critical asymmetry in renormalization group theory for fluids

The renormalization-group (RG) approaches for fluids are employed to investigate critical asymmetry of vapour-liquid equilibrium (VLE) of fluids. Three different approaches based on RG theory for fluids are reviewed and compared. RG approaches are applied to various fluid systems: hard-core square-well fluids of variable ranges, hard-core Yukawa fluids, and square-well dimer fluids and modelling VLE of n-alkane molecules. Phase diagrams of simple model fluids and alkanes described by RG approaches are analyzed to assess the capability of describing the VLE critical asymmetry which is suggested in complete scaling theory. Results of thermodynamic properties obtained by RG theory for fluids agree with the simulation and experimental data. Coexistence diameters, which are smaller than the critical densities, are found in the RG descriptions of critical asymmetries of several fluids. Our calculation and analysis show that the approach coupling local free energy with White's RG iteration which aims to incorporate density fluctuations into free energy is not adequate for VLE critical asymmetry due to the inadequate order parameter and the local free energy functional used in the partition function.

Coupling-parameter expansion in thermodynamic perturbation theory

An approach to the coupling-parameter expansion in the liquid state theory of simple fluids is presented by combining the ideas of thermodynamic perturbation theory and integral equation theories. This hybrid scheme avoids the problems of the latter in the two phase region. A method to compute the perturbation series to any arbitrary order is developed and applied to square well fluids. Apart from the Helmholtz free energy, the method also gives the radial distribution function and the direct correlation function of the perturbed system. The theory is applied for square well fluids of variable ranges and compared with simulation data. While the convergence of perturbation series and the overall performance of the theory is good, improvements are needed for potentials with shorter ranges. Possible directions for further developments in the coupling-parameter expansion are indicated.

Statistical thermodynamics of fluids with both dipole and quadrupole moments

New Gibbs ensemble simulation data for a polar fluid modeled by a square-well potential plus dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole interactions are presented. This simulation data is used in order to assess the applicability of the multipolar square-well perturbation theory [A. L. Benavides, Y. Guevara, and F. del Río, Physica A 202, 420 (1994)] to systems where more than one term in the multipole expansion is relevant. It is found that this theory is able to reproduce qualitatively well the vapor-liquid phase diagram for different multipolar moment strengths, corresponding to typical values of real molecules, except in the critical region. Hence, this theory is used to model the behavior of substances with multiple chemical bonds such as carbon monoxide and nitrous oxide and we found that with a suitable choice of the values of the intermolecular parameters, the vapor-liquid equilibrium of these species is adequately estimated.

Molecular dynamics in the isothermal-isobaric ensemble: the requirement of a "shell" molecule. III. Discontinuous potentials

Based on the approach of Gruhn and Monson [Phys. Rev. E 63, 061106 (2001)], we present a new method for deriving the collisions dynamics for particles that interact via discontinuous potentials. By invoking the conservation of the extended Hamiltonian, we generate molecular dynamics (MD) algorithms for simulating the hard-sphere and square-well fluids within the isothermal-isobaric (NpT) ensemble. Consistent with the recent rigorous reformulation of the NpT ensemble partition function, the equations of motion impose a constant external pressure via the introduction of a shell particle of known mass [M. J. Uline and D. S. Corti, J. Chem. Phys. 123, 164101 (2005); 123, 164102 (2005)], which serves to define uniquely the volume of the system. The particles are also connected to a temperature reservoir through the use of a chain of Nose-Hoover thermostats, the properties of which are not affected by a hard-sphere or square-well collision. By using the Liouville operator formalism and the Trotter expansion theorem to integrate the equations of motion, the update of the thermostat variables can be decoupled from the update of the positions of the particles and the momentum changes upon a collision. Hence, once the appropriate collision dynamics for the isobaric-isenthalpic (NpH) equations of motion is known, the adaptation of the algorithm to the NpT ensemble is straightforward. Results of MD simulations for the pure component square-well fluid are presented and serve to validate our algorithm. Finally, since the mass of the shell particle is known, the system itself, and not a piston of arbitrary mass, controls the time scales for internal pressure and volume fluctuations. We therefore consider the influence of the shell particle algorithm on the dynamics of the square-well fluid.

Thermodynamic properties of short-range square well fluid

The interfacial properties of short-range square well fluid with lambda=1.15, 1.25, and 1.375 were determined by using single canonical Monte Carlo simulations. Simulations were carried out in the vapor-liquid region. The coexistence curves of these models were calculated and compared to those previously reported in the literature and good agreement was found among them. We found that the surface tension curves for any potential model of short range form a single master curve when we plot gamma* vs TT(c). It is demonstrated that the critical reduced second virial coefficient B(2)* as a function of interaction range or T(c)* is not constant.

Statistical Associating Fluid Theory Coupled with Restrictive Primitive Model Extended to Bivalent Ions. SAFT2: 1. Single Salt + Water Solutions

Statistical associating fluid theory coupled with the restricted primitive model is extended to multivalent ions by relaxing the range of the square-well width parameter, which leads to a new dispersion term approximation and calls for a new set of salt and ion parameters. This new approximation, referred to as SAFT2, requires a single set of parameters derived from the salt (mean ionic) activity coefficients and liquid densities of single-salt solutions for five cations (Li(+), Na(+), K(+), Ca(2+), Mg(2+)), six anions (Cl(-), Br(-), I(-), NO(3)(-), SO(4)(-2), HCO(3)(-)), and 24 salts. These parameters, in turn, are shown to predict the osmotic coefficients for single salt + water solutions.

A general perturbation approach for equation of state development: Applications to simple fluids, ab initio potentials, and fullerenes

A new perturbation scheme based on the Barker-Henderson perturbation theory [J. Chem. Phys. 47, 4714 (1967)] is proposed to predict the thermodynamic properties of spherical molecules. Accurate predictions of second virial coefficients and vapor-liquid coexistence properties are obtained for a large variety of potential functions (square well, Yukawa, Sutherland, Lennard-Jones, Buckingham, Girifalco). New Gibbs ensemble Monte Carlo simulations of the generalized exp-m Buckingham potential are reported. An extension of the perturbation approach to mixtures is proposed, and excellent predictions of vapor-liquid equilibria are obtained for Lennard-Jones mixtures. The perturbation scheme can be applied to complex potential functions fitted to ab initio data to predict the properties of real molecules such as neon. The new approach can also be used as an auxiliary tool in molecular simulation studies, to efficiently optimize an intermolecular potential on macroscopic properties or match force fields based on different potential functions.

Vapor-liquid equilibrium and critical behavior of the square-well fluid of variable range: A theoretical study

The vapor-liquid phase behavior and the critical behavior of the square-well (SW) fluid are investigated as a function of the interaction range, lambdain [1.25, 3], by means of the self-consistent Ornstein-Zernike approximation (SCOZA) and analytical equations of state based on a perturbation theory [A. L. Benavides and F. del Rio, Mol. Phys. 68, 983 (1989); A. Gil-Villegas, F. del Rio, and A. L. Benavides, Fluid Phase Equilib. 119, 97 (1996)]. For this purpose the SCOZA, which has been restricted up to now to a few model systems, has been generalized to hard-core systems with arbitrary interaction potentials requiring a fully numerical solution of an integro-partial differential equation. Both approaches, in general, describe well the liquid-vapor phase diagram of the square-well fluid when compared with simulation data. SCOZA yields very precise predictions for the coexistence curves in the case of long ranged SW interaction (lambda>1.5), and the perturbation theory is able to predict the binodal curves and the saturated pressures, for all interaction ranges considered if one stays away from the critical region. In all cases, the SCOZA gives very good predictions for the critical temperatures and the critical pressures, while the perturbation theory approach tends to slightly overestimate these quantities. Furthermore, we propose analytical expressions for the critical temperatures and pressures as a function of the square-well range.

Predicting Mixture Phase Equilibria and Critical Behavior Using the SAFT-VRX Approach

The SAFT-VRX equation of state combines the SAFT-VR equation with a crossover function that smoothly transforms the classical equation into a nonanalytical form close to the critical point. By a combinination of the accuracy of the SAFT-VR approach away from the critical region with the asymptotic scaling behavior seen at the critical point of real fluids, the SAFT-VRX equation can accurately describe the global fluid phase diagram. In previous work, we demonstrated that the SAFT-VRX equation very accurately describes the pvT and phase behavior of both nonassociating and associating pure fluids, with a minimum of fitting to experimental data. Here, we present a generalized SAFT-VRX equation of state for binary mixtures that is found to accurately predict the vapor-liquid equilibrium and pvT behavior of the systems studied. In particular, we examine binary mixtures of n-alkanes and carbon dioxide + n-alkanes. The SAFT-VRX equation accurately describes not only the gas-liquid critical locus for these systems but also the vapor-liquid equilibrium phase diagrams and thermal properties in single-phase regions.