Andrew J. Haslam

Developing optimal Wertheim-like models of water for use in Statistical Associating Fluid Theory (SAFT) and related approaches

The statistical associating fluid theory (SAFT) is now well established as an approach for the description of the thermodynamics and phase equilibria of a wide variety of fluid systems. Numerous SAFT studies of the fluid phase equilibria of pure water and aqueous mixtures have been made with the various incarnations of the theory, yet there is no consensus on what the ‘optimal’ values of the intermolecular parameters are for water, or what the association scheme should be (two-, three-, or four-site models). We show that the conventional use of vapour-pressure and saturated-liquid-density data on their own leads to a degeneracy in the values of the parameters, particularly in the relative values of the dispersion and hydrogen-bonding energies. A discretized energy-parameter space is examined and a long valley of ‘optimal’ parameter sets for the vapour–liquid equilibria is found for the various models, ranging from low dispersion (high hydrogen-bonding) to high dispersion (low hydrogen-bonding) energies. Other thermodynamic information such as the heat of vaporization does not allow one to discriminate between the values of the parameters or the association scheme. An examination of the degree of association (hydrogen bonding) allows such a discrimination to be made: the four-site model is found to provide the best overall description of the thermodynamics, fluid phase equilibria, and degree of association. The set of parameters obtained in this way also provides the best description of the fluid-phase equilibria for a mixture of water and methanol. This degeneracy of parameters could also be important when models of water are refined to vapour–liquid equilibria in simulation studies of aqueous systems, and for other systems of associating molecules.

The phase behavior of a binary mixture of rodlike and disclike mesogens: Monte Carlo simulation, theory, and experiment

The phase behavior of a binary mixture of rodlike and disclike hard molecules is studied using Monte Carlo NVT (constant number of particles N, volume V, and temperature T) computer simulation. The rods are modeled as hard spherocylinders of aspect ratio LHSC/DHSC=5, and the discs as hard cut spheres of aspect ratio LCS/DCS=0.12. The diameter ratio DCS/DHSC=3.62 is chosen such that the molecular volumes of the two particles are equal. The starting configuration in the simulations is a mixed isotropic state. The phase diagram is mapped by changing the overall density of the system. At low densities stabilization of the isotropic phase relative to the ordered states is seen on mixing, and at high densities nematic–columnar and smectic A–columnar phase coexistence is observed. Biaxiality in the nematic phase is not seen. The phase diagram of the mixture is also calculated using the second virial theory of Onsager for nematic ordering, together with the scaling of Parsons and Lee to take into account the high...

A generalisation of the Onsager trial function approach: Describing nematic liquid crystals with an algebraic equation of state

The molecular theory of Onsager [L. Onsager, Ann. N.Y. Acad. Sci. 51, 627 (1949)] for liquid crystals is developed and extended to describe ordering transitions in systems of generic cylindrically symmetrical molecules. A number of new analytical results are discussed for particles characterised by a general form of the excluded-volume interaction. Our description makes use of the Onsager trial function (OTF) to represent the orientational distribution and degree of anisotropy. Algebraic expressions for the thermodynamic properties, which provide a particularly tractable description of the isotropic–nematic equilibria, are also presented. The degree of orientational order can be represented by a simple cubic equation in the molecular parameters (molecular diameter and aspect ratio) and thermodynamic variables (temperature and number density). Onsagers theory was originally developed at the level of the second virial coefficient; here the Parsons–Lee decoupling approximation is used to describe the higher-body contributions in a straightforward manner. The adequacy of treating the scaled Onsager (Parsons–Lee) free-energy functional within the OTF formalism to describe anisotropic states is illustrated by examining systems of hard spherocylinders. An excellent representation of the equation of state of the isotropic and nematic phases and the ordering transition is demonstrated for molecules of moderate aspect ratio (L/D = 5). Algebraic equations of state of the type developed here are suitable for practical engineering applications involving anisotropic fluids particularly in the case of multicomponent systems; our general analytical results for the averages of orientational functions will turn out to be useful in the development of a description of molecules with more realistic attractive and Maier–Saupe interactions.

The A in SAFT: developing the contribution of association to the Helmholtz free energy within a Wertheim TPT1 treatment of generic Mie fluids

An accurate representation of molecular association is a vital ingredient of advanced equations of state (EOSs), providing a description of thermodynamic properties of complex fluids where hydrogen bonding plays an important role. The combination of the first-order thermodynamic perturbation theory (TPT1) of Wertheim for associating systems with an accurate description of the structural and thermodynamic properties of the monomer fluid forms the basis of the statistical associating fluid theory (SAFT) family of EOSs. The contribution of association to the free energy in SAFT and related EOSs is very sensitive to the nature of intermolecular potential used to describe the monomers and, crucially, to the accuracy of the representation of the thermodynamic and structural properties. Here we develop an accurate description of the association contribution for use within the recently developed SAFT-VR Mie framework for chain molecules formed from segments interacting through a Mie potential [T. Lafitte, A. Apostolakou, C. Avendaño, A, Galindo, C. S. Adjiman, E. A. Müller, and G. Jackson, J. Chem. Phys. 139, 154504 (2013)]. As the Mie interaction represents a soft-core potential model, a method similar to that adopted for the Lennard-Jones potential [E. A. Müller and K. E. Gubbins, Ind. Eng. Chem. Res. 34, 3662 (1995)] is employed to describe the association contribution to the Helmholtz free energy. The radial distribution function (RDF) of the Mie fluid (which is required for the evaluation of the integral at the heart of the association term) is determined for a broad range of thermodynamic conditions (temperatures and densities) using the reference hyper-netted chain (RHNC) integral-equation theory. The numerical data for the association kernel of Mie fluids with different association geometries are then correlated for a range of thermodynamic states to obtain a general expression for the association contribution which can be applied for varying values of the Mie repulsive exponent. The resulting SAFT-VR Mie EOS allows for a much improved description of the vapour-liquid equilibria and single-phase properties of associating fluids such as water, methanol, ammonia, hydrogen sulphide, and their mixtures. A comparison is also made between the theoretical predictions of the degree of association for water and the extent of hydrogen bonding obtained from molecular simulations of the SPC/E and TIP4P/2005 atomistic models.

Pseudo hard-sphere potential for use in continuous molecular-dynamics simulation of spherical and chain molecules

We present a continuous pseudo-hard-sphere potential based on a cut-and-shifted Mie (generalized Lennard-Jones) potential with exponents (50, 49). Using this potential one can mimic the volumetric, structural, and dynamic properties of the discontinuous hard-sphere potential over the whole fluid range. The continuous pseudo potential has the advantage that it may be incorporated directly into off-the-shelf molecular-dynamics code, allowing the user to capitalise on existing hardware and software advances. Simulation results for the compressibility factor of the fluid and solid phases of our pseudo hard spheres are presented and compared both to the Carnahan-Starling equation of state of the fluid and published data, the differences being indistinguishable within simulation uncertainty. The specific form of the potential is employed to simulate flexible chains formed from these pseudo hard spheres at contact (pearl-necklace model) for m(c) = 4, 5, 7, 8, 16, 20, 100, 201, and 500 monomer segments. The compressibility factor of the chains per unit of monomer, m(c), approaches a limiting value at reasonably small values, m(c) < 50, as predicted by Wertheims first order thermodynamic perturbation theory. Simulation results are also presented for highly asymmetric mixtures of pseudo hard spheres, with diameter ratios of 3:1, 5:1, 20:1 over the whole composition range.

Modelling of the thermodynamic and solvation properties of electrolyte solutions with the statistical associating fluid theory for potentials of variable range

An improved formulation of the extension of the statistical associating fluid theory for potentials of variable range to electrolytes (SAFT-VRE) is presented, incorporating a representation for the dielectric constant of the solution that takes into account the temperature, density and composition of the solvent. The proposed approach provides an excellent correlation of the dielectric-constant data available for a number of solvents including water, representative alcohols and carbon dioxide, and it is shown that the methodology can be used to treat mixed-solvent electrolyte solutions. Models for strong electrolytes of the metal-halide family are considered here. The salts are treated as fully dissociated and ion-specific interaction parameters are presented. Vapour pressure, density, and mean ionic activity coefficient data are used to determine the ion–ion and solvent–ion parameters, and mixed-salt electrolyte solutions (brines) are then treated predictively. We find that the resulting intermolecular potential models follow physical trends in terms of energies and ion sizes with a close relationship observed with well-established ionic diameters. A good description is obtained for the densities, mean ionic activity coefficients, and vapour pressures of the electrolyte solutions studied. The theory is also seen to provide excellent predictions of the osmotic coefficient and of the depression of the freezing temperature, and provides a qualitative estimate of the solvation free energy. The vapour pressure of aqueous brines is predicted accurately, as is the density of these solutions, although not at the highest pressures considered. Calculations for the vapour–liquid and liquid–liquid equilibria of salts in water+methanol and water+n-butan-1-ol are presented. In addition, it is shown that the salting-out of carbon dioxide in sodium chloride solutions is captured well using a predictive model.

Advances in generalised van der Waals approaches for the isotropic–nematic fluid phase equilibria of thermotropic liquid crystals–an algebraic equation of state for attractive anisotropic particles with the Onsager trial function

In this contribution we provide a review and reformulation of perturbation theories (generalised van der Waals approaches) for the description of the fluid phase behaviour and orientational ordering transitions of thermotropic nematic liquid crystals. Free-energy functionals of the basic Onsager form are used as the platform for the development of a general formulation that reduces to the specific forms of the various theories that have found common use. A novel closed analytical description of the thermodynamic properties and degree of nematic order is then obtained by employing the Onsager trial function to represent the orientational distribution function in terms of a single parameter. The latter essentially constitutes an algebraic equation of state for the nematic phase appropriate for use in engineering applications. The description of ordering transitions with scaled-Onsager theories and suitable trial functions has already been illustrated by its application to systems of hard spherocylinders (HSCs), indicating that the approach provides an excellent representation of the orientational order of the hard-core (athermal) system [Mol. Phys. 106, 649 (2008)]. Here, the hard-body model is extended to account for attractive interactions (treated at the van der Waals level) of a general isotropic/anisotropic form (e.g., Lennard-Jonesium (LJ), square-well (SW), Maier–Saupe (MS), etc.). The adequacy of our generalised van der Waals–Onsager theory is exemplified by an analysis of the vapour–liquid, liquid–nematic, and vapour–nematic phase equilibria for hard spherocylinders with attractive square-well interactions (HSC-SW). The effect of the potential range and molecular aspect ratio on the vapour–liquid–nematic equilibria and orientational ordering transitions is examined to investigate the van der Waals limit (corresponding states) of the ordering phase behaviour. In the case of systems with an aspect ratio of ∼5 the corresponding-states limit is reached when the range is about 16 times the molecular diameter. For progressively longer molecules with an attractive range that follows their long dimension, the fluid–nematic equilibrium is enhanced to the point that the vapour–liquid boundary becomes metastable relative to fluid–nematic equilibria. In the case of molecules of very large aspect ratio (∼50), an additional region of nematic–nematic coexistence is exhibited by the system.

Modelling the effect of methanol, glycol inhibitors and electrolytes on the equilibrium stability of hydrates with the SAFT-VR approach

In this work we integrate the statistical associating fluid theory for fluids interacting through potentials of variable range (SAFT-VR) into a traditional van der Waals and Platteeuw framework for modelling clathrate hydrates. We incorporate a new water–guest cell potential for the hydrate phase that can be related to the potential adopted in the familiar SAFT-VR equation of state for modelling fluids. We show how the ability of this equation of state to treat a wide range of complex fluids increases the scope of hydrate modelling to incorporate, in a single framework, the presence of various inhibitors (alcohols, glycols) or brines – or, indeed, any fluid for which a model is available (for use within SAFT-VR) or can be conveniently obtained. Agreement with experimental results is good throughout and, in many cases, excellent.

Subtleties in the calculation of the pressure and pressure tensor of anisotropic particles from volume-perturbation methods and the apparent asymmetry of the compressive and expansive contributions

An efficient and versatile method to calculate the components of the pressure tensor for hard-body fluids of generic shape from the perspective of molecular simulation is presented. After due consideration of all the possible repulsive contributions exerted by molecules upon their surroundings during an anisotropic system expansion, it is observed that such a volume change can, for non-spherical molecules, give rise to configurations where overlaps occur. This feature of anisotropic molecules has to be taken into account rigorously as it can lead to discrepancies in the calculation of tensorial contributions to the pressure. Using the condition of detailed balance as a basis, a perturbation method developed for spherical molecules has been extended so that it is applicable to non-spherical and non-convex molecules. From a series of ‘ghost’ anisotropic volume perturbations the residual contribution to the components of the pressure tensor may be accurately calculated. Comparisons are made with prior methods and, where relevant, results are evaluated against existing data. For inhomogeneous systems this method provides a particularly convenient route to the calculation of the interfacial tension (surface free energy) from molecular simulations.

Development of intermolecular potential models for electrolyte solutions using an electrolyte SAFT-VR Mie equation of state

ABSTRACT We present a theoretical framework and parameterisation of intermolecular potentials for aqueous electrolyte solutions using the statistical associating fluid theory based on the Mie interaction potential (SAFT-VR Mie), coupled with the primitive, non-restricted mean-spherical approximation (MSA) for electrolytes. In common with other SAFT approaches, water is modelled as a spherical molecule with four off-centre association sites to represent the hydrogen-bonding interactions; the repulsive and dispersive interactions between the molecular cores are represented with a potential of the Mie (generalised Lennard-Jones) form. The ionic species are modelled as fully dissociated, and each ion is treated as spherical: Coulombic ion–ion interactions are included at the centre of a Mie core; the ion–water interactions are also modelled with a Mie potential without an explicit treatment of ion–dipole interaction. A Born contribution to the Helmholtz free energy of the system is included to account for the process of charging the ions in the aqueous dielectric medium. The parameterisation of the ion potential models is simplified by representing the ion–ion dispersive interaction energies with a modified version of the London theory for the unlike attractions. By combining the Shannon estimates of the size of the ionic species with the Born cavity size reported by Rashin and Honig, the parameterisation of the model is reduced to the determination of a single ion–solvent attractive interaction parameter. The resulting SAFT-VRE Mie parameter sets allow one to accurately reproduce the densities, vapour pressures, and osmotic coefficients for a broad variety of aqueous electrolyte solutions; the activity coefficients of the ions, which are not used in the parameterisation of the models, are also found to be in good agreement with the experimental data. The models are shown to be reliable beyond the molality range considered during parameter estimation. The inclusion of the Born free-energy contribution, together with appropriate estimates for the size of the ionic cavity, allows for accurate predictions of the Gibbs free energy of solvation of the ionic species considered. The solubility limits are also predicted for a number of salts; in cases where reliable reference data are available the predictions are in good agreement with experiment.