A statistical associating fluid theory for electrolyte solutions (SAFT-VRE)

Toward Modeling the Structure of Electrolytes at Charged Mineral Interfaces Using Classical Density Functional Theory

The organization of water molecules and ions between charged mineral surfaces determines the stability of colloidal suspensions and the strength of phase-separated particulate gels. In this article, we assemble a density functional that measures the free energy due to the interaction of water molecules and ions in electric double layers. The model accounts for the finite size of the particles using fundamental measure theory, hydrogen-bonding between water molecules using Wertheim's statistical association theory, long-range dispersion interactions using Barker and Henderson's high-temperature expansion, electrostatic correlations using a functionalized mean-spherical approximation, and Coulomb forces through the Poisson equation. These contributions are shown to produce highly correlated structures, aptly rendering the layering of counterions and co-ions at highly charged surfaces and permitting the solvation of ions and surfaces to be measured by a combination of short-range associations and long-ranged attractions. The model is tested in a planar geometry near soft, charged surfaces to reproduce the structure of water near graphene and mica. For mica surfaces, explicitly representing the density of the outer oxygen layer of the exposed silica tetrahedra allows water molecules to hydrogen-bond to the solid. When electrostatic interactions are included, water molecules assume a hybrid character, being accounted for implicitly in the dielectric constant but explicitly otherwise. The disjoining pressure between approaching like-charged surfaces is calculated, demonstrating the model's ability to probe pressure oscillations that arise during the expulsion of ions and water layers from the interfacial gap and predict strong interattractive stresses that form at narrow interfacial spacing when the surface charge is overscreened. This interattractive stress arises not due to in-plane correlations under strong electrostatic coupling but due to the out-of-plane structuring of associating ions and water molecules.

The Contribution of the Ion-Ion and Ion-Solvent Interactions in a Molecular Thermodynamic Treatment of Electrolyte Solutions

Developing molecular equations of state to treat electrolyte solutions is challenging due to the long-range nature of the Coulombic interactions. Seminal approaches commonly used are the mean spherical approximation (MSA) and the Debye-Hückel (DH) theory to account for ion-ion interactions and, often, the Born theory of solvation for ion-solvent interactions. We investigate the accuracy of the MSA and DH approaches using each to calculate the contribution of the ion-ion interactions to the chemical potential of NaCl in water, comparing these with newly computer-generated simulation data; the ion-ion contribution is isolated by selecting an appropriate primitive model with a Lennard-Jones force field to describe the solvent. A study of mixtures with different concentrations and ionic strengths reveals that the calculations from both MSA and DH theories are of similar accuracy, with the MSA approach resulting in marginally better agreement with the simulation data. We also demonstrate that the Born theory provides a good qualitative description of the contribution of the ion-solvent interactions; we employ an explicitly polar water model in these simulations. Quantitative agreement up to moderate salt concentrations and across the relevant range of temperature is achieved by adjusting the Born radius using simulation data of the free energy of solvation. We compute the radial and orientational distribution functions of the systems, thereby providing further insight on the differences observed between the theory and simulation. We thus provide rigorous benchmarks for use of the MSA, DH, and Born theories as perturbation approaches, which will be of value for improving existing models of electrolyte solutions, especially in the context of equations of state.

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.

Predictive models for the phase behaviour and solution properties of weak electrolytes: nitric, sulphuric, and carbonic acids

The distribution of ionic species in electrolyte systems is important in many fields of science and engineering, ranging from the study of degradation mechanisms to the design of systems for electrochemical energy storage. Often, other phenomena closely related to ionic speciation, such as ion pairing, clustering and hydrogen bonding, which are difficult to investigate experimentally, are also of interest. Here, we develop an accurate molecular approach, accounting for reactions as well as association and ion pairing, to deliver a predictive framework that helps validate experiment and guides future modelling of speciation phenomena of weak electrolytes. We extend the SAFT-VRE Mie equation of state [D. K. Eriksen et al., Mol. Phys., 2016, 114, 2724-2749] to study aqueous solutions of nitric, sulphuric, and carbonic acids, considering complete and partially dissociated models. In order to incorporate the dissociation equilibria, correlations to experimental data for the relevant thermodynamic equilibrium constants of the dissociation reactions are taken from the literature and are imposed as a boundary condition in the calculations. The models for water, the hydronium ion, and carbon dioxide are treated as transferable and are taken from our previous work. We present new molecular models for nitric acid, and the nitrate, bisulfate, sulfate, and bicarbonate anions. The resulting framework is used to predict a range of phase behaviour and solution properties of the aqueous acids over wide ranges of concentration and temperature, including the degree of dissociation, as well as the activity coefficients of the ionic species, and the activity of water and osmotic coefficient, density, and vapour pressure of the solutions. The SAFT-VRE Mie models obtained in this manner provide a means of elucidating the mechanisms of association and ion pairing in the systems studied, complementing the experimental observations reported in the literature.

Thermodynamic properties of confined square-well fluids with multiple associating sites

In this work, a molecular simulation study of confined hard-spheres particles with square-well (SW) attractive interactions with two and four associating SW sites based on the first-order perturbation form of Wertheim's theory is presented. An extended version of the Gibbs ensemble technique for inhomogeneous fluids [A. Z. Panagiotopoulos, Mol. Phys. 62, 701 (1987)] is used to predict the adsorption density profiles for associating fluids confined between opposite parallel walls. The fluid is confined in four kinds of walls: hard-wall, SW wall, Lennard-Jones (LJ) 12-6 wall potential, and LJ 10-4 wall potential. We analyze the behavior of the confined system for several supercritical temperatures as a function of variation of molecular parameters: potential range λ, bulk densities ρ_b^*, pore width H, cutoff range interaction r_c^*, and range of the potential and depth of the particle-wall (λ_w, ε_w^*). Additionally, we include predictions for liquid-vapor coexistence of bulk associative particles and how their critical properties are modified by the presence of associative sites in the molecule. The molecular simulation data presented in this work are of prime importance to the development of theoretical approaches for inhomogeneous fluids as classical density functional theory. The simulation results presented here are resourceful for predicting adsorption isotherms of real associating fluids such as water.

On the use of molecular-based thermodynamic models to assess the performance of solvents for CO2 capture processes: monoethanolamine solutions

Predictive models play an important role in the design of post-combustion processes for the capture of carbon dioxide (CO₂) emitted from power plants. A rate-based absorber model is presented to investigate the reactive capture of CO₂ using aqueous monoethanolamine (MEA) as a solvent, integrating a predictive molecular-based equation of state: SAFT-VR SW (Statistical Associating Fluid Theory-Variable Range, Square Well). A distinctive physical approach is adopted to model the chemical equilibria inherent in the process. This eliminates the need to consider reaction products explicitly and greatly reduces the amount of experimental data required to model the absorber compared to the more commonly employed chemical approaches. The predictive capabilities of the absorber model are analyzed for profiles from 10 pilot plant runs by considering two scenarios: (i) no pilot-plant data are used in the model development; (ii) only a limited set of pilot-plant data are used. Within the first scenario, the mass fraction of CO₂ in the clean gas is underestimated in all but one of the cases, indicating that a best-case performance of the solvent can be obtained with this predictive approach. Within the second scenario a single parameter is estimated based on data from a single pilot plant run to correct for the dramatic changes in the diffusivity of CO₂ in the reactive solvent. This parameter is found to be transferable for a broad range of operating conditions. A sensitivity analysis is then conducted, and the liquid viscosity and diffusivity are found to be key properties for the prediction of the composition profiles. The temperature and composition profiles are sensitive to thermodynamic properties that correspond to major sources of heat generation or dissipation. The proposed modelling framework can be used as an early assessment of solvents to aid in narrowing the search space, and can help in determining target solvents for experiments and more detailed modelling.

Molecular thermodynamic modeling of ionic liquids using the perturbation-based linear Yukawa isotherm regularity

In this paper, we present the results of an extensive study on a novel approach to the molecular modeling of pure ionic liquids (ILs) that incorporates the perturbed thermodynamic linear Yukawa isotherm regularity (LYIR), which is derived based on an effective nearest neighboring pair attractive interaction of the Yukawa potential. The LYIR was used to model the densities of ILs up to high pressures (35 MPa) and in the temperature range 293.15 to 393.15 K. To use the LYIR for ILs, a simple molecular model was proposed to describe their molecular structure, in which they were considered as a liquid consisting of the ion pairs moving together in the fluid, and each ion pair was assumed to be a one-center spherical united atom. The ILs under consideration contained one of the IL cations [C2mim](+), [C4mim](+), [C7mim](+), [C8mim](+), [C3mpy](+), [C3mpip](+), [C3mpyr](+) or [C4mpyr](+), and one of the IL anions [BF4](-), [C(CN)3](-), [CF3SO4](-) or [NTf2](-). The reliability and physical significance of the parameters as well as the proposed molecular model were tested by calculating the densities of pure imidazolium-, pyridinium-, piperidinium- and pyrrolidimium-based ILs. The results showed that the LYIR can be used to predict and reproduce the density of ILs in good agreement with the experimental data. In addition, the LYIR enabled us to determine the physical quantities, such as an effective Yukawa screening length, λ eff, the product of the effective energy well depth and the effective coordination number, (ε eff/k)z eff, the contribution of the non-reference thermal pressure and also the influence of the anionic and cationic structure on the λ eff parameter. The standard deviation of the IL densities predicted in this work is lower than those calculated by the one other important equation of state reported in the literature.

Force-Field Parameters from the SAFT-γ Equation of State for Use in Coarse-Grained Molecular Simulations

A description of fluid systems with molecular-based algebraic equations of state (EoSs) and by direct molecular simulation is common practice in chemical engineering and the physical sciences, but the two approaches are rarely closely coupled. The key for an integrated representation is through a well-defined force field and Hamiltonian at the molecular level. In developing coarse-grained intermolecular potential functions for the fluid state, one typically starts with a detailed, bottom-up quantum-mechanical or atomic-level description and then integrates out the unwanted degrees of freedom using a variety of techniques; an iterative heuristic simulation procedure is then used to refine the parameters of the model. By contrast, with a top-down technique, one can use an accurate EoS to link the macroscopic properties of the fluid and the force-field parameters. We discuss the latest developments in a top-down representation of fluids, with a particular focus on a group-contribution formulation of the statistical associating fluid theory (SAFT-γ). The accurate SAFT-γ EoS is used to estimate the parameters of the Mie force field, which can then be used with confidence in direct molecular simulations to obtain thermodynamic, structural, interfacial, and dynamical properties that are otherwise inaccessible from the EoS. This is exemplified for several prototypical fluids and mixtures, including carbon dioxide, hydrocarbons, perfluorohydrocarbons, and aqueous surfactants.

SAFT-γ force field for the simulation of molecular fluids: 2. Coarse-grained models of greenhouse gases, refrigerants, and long alkanes

In the first paper of this series [C. Avendaño, T. Lafitte, A. Galindo, C. S. Adjiman, G. Jackson, and E. A. Müller, J. Phys. Chem. B2011, 115, 11154] we introduced the SAFT-γ force field for molecular simulation of fluids. In our approach, a molecular-based equation of state (EoS) is used to obtain coarse-grained (CG) intermolecular potentials that can then be employed in molecular simulation over a wide range of thermodynamic conditions of the fluid. The macroscopic experimental data for the vapor-liquid equilibria (saturated liquid density and vapor pressure) of a given system are represented with the SAFT-VR Mie EoS and used to estimate effective intermolecular parameters that provide a good description of the thermodynamic properties by exploring a wide parameter space for models based on the Mie (generalized Lennard-Jones) potential. This methodology was first used to develop a simple single-segment CG Mie model of carbon dioxide (CO2) which allows for a reliable representation of the fluid-phase equilibria (for which the model was parametrized), as well as an accurate prediction of other properties such as the enthalpy of vaporization, interfacial tension, supercritical density, and second-derivative thermodynamic properties (thermal expansivity, isothermal compressibility, heat capacity, Joule-Thomson coefficient, and speed of sound). In our current paper, the methodology is further applied and extended to develop effective SAFT-γ CG Mie force fields for some important greenhouse gases including carbon tetrafluoride (CF4) and sulfur hexafluoride (SF6), modeled as simple spherical molecules, and for long linear alkanes including n-decane (n-C10H22) and n-eicosane (n-C20H42), modeled as homonuclear chains of spherical Mie segments. We also apply the SAFT-γ methodology to obtain a CG homonuclear two-segment Mie intermolecular potential for the more challenging polar and asymmetric compound 2,3,3,3-tetrafluoro-1-propene (HFO-1234yf), a novel replacement refrigerant with promising properties. The description of the fluid-phase behavior and the prediction of the other thermophysical properties obtained by molecular simulation using our SAFT-γ CG Mie force fields are found to be of comparable quality (and sometimes superior) to that obtained using the more sophisticated all-atom (AA) and united-atom (UA) models commonly employed in the field. We should emphasize that though the focus of our current work is on simple homonuclear models, the SAFT-γ methodology is based on a group contribution methodology which is naturally suited to the development of more sophisticated heteronuclear models.

Anomalous columnar order of charged colloidal platelets

Monte Carlo computer simulations are carried out for a model system of like-charged colloidal platelets in the isothermal-isobaric ensemble (NpT). The aim is to elucidate the role of electrostatic interactions on the structure of synthetic clay systems at high particle densities. Short-range repulsions between particles are described by a suitable hard-core model representing a discotic particle. This potential is supplemented with an electrostatic potential based on a Yukawa model for the screened Coulombic potential between infinitely thin disklike macro-ions. The particle aspect-ratio and electrostatic parameters were chosen to mimic an aqueous dispersion of thin, like-charged, rigid colloidal platelets at finite salt concentration. An examination of the fluid phase diagram reveals a marked shift in the isotropic-nematic transition compared to the hard cut-sphere reference system. Several statistical functions, such as the pair correlation function for the center-of-mass coordinates and structure factor, are obtained to characterize the structural organization of the platelets phases. At low salinity and high osmotic pressure we observe anomalous hexagonal columnar structures characterized by interpenetrating columns with a typical intercolumnar distance corresponding to about half of that of a regular columnar phase. Increasing the ionic strength leads to the formation of glassy, disordered structures consisting of compact clusters of platelets stacked into finite-sized columns. These so-called "nematic columnar" structures have been recently observed in systems of charge-stabilized gibbsite platelets. Our findings are corroborated by an analysis of the static structure factor from a simple density functional theory.

Development of an equation of state for electrolyte solutions by combining the statistical associating fluid theory and the mean spherical approximation for the nonprimitive model

A statistical associating fluid theory to model electrolyte fluids that explicitly accounts for solvent molecules by modeling water as a dipolar square-well associating fluid is presented. Specifically the statistical associating fluid theory for potentials of variable range (SAFT-VR) is combined with integral equation theory and the generalized mean spherical approximation using the nonprimitive model to describe the long-range ion-ion, ion-dipole, and dipole-dipole interactions. Isothermal-isobaric ensemble Monte Carlo simulations have been performed in order to test the new theoretical approach. In particular, simulations are performed for different ion concentrations and different ratios of the cation, anion, and solvent segment diameters. Predictions for the thermodynamic properties from the new equation of state are compared with the computer simulation data. Additionally, results from a combination of the SAFT-VR approach with Debye-Huckel theory and the primitive model are also presented and compared to those obtained with the nonprimitive model to illustrate the advantages of the new statistical associating fluid theory for potentials of variable range plus dipole and electrolytes (SAFT-VR+DE) approach. The results show that the proposed equation of state provides a good description of the PVT properties of electrolyte fluids with different sizes of ions and solvent.

Liquid Phase Behavior of Perfluoroalkylalkane Surfactants

We have performed a combined experimental and theoretical study of the thermodynamic properties of semifluorinated alkanes. In particular, the liquid density of perfluorohexylhexane (F6H6) and perfluorohexyloctane (F6H8) has been measured as a function of temperature from 273.15 to 353.15 K and at four temperatures as a function of pressure up to 600 bar. The results were interpreted using the SAFT-VR equation of state. The perfluoroalkylalkanes were modeled as heterosegmented diblock chains using parameters for the alkyl and perfluoroalkyl segments developed in earlier work. Through this simple approach, we are able to predict the thermodynamic behavior of the perfluoroalkylalkanes studied without fitting to any experimental data for the systems being studied.

Predicting the Phase Behavior of Nitrogen + n-Alkanes for Enhanced Oil Recovery from the SAFT-VR Approach: Examining the Effect of the Quadrupole Moment

The phase behavior of nitrogen + n-alkanes is studied within the framework of the statistical associating fluid theory for potentials of variable range (SAFT-VR). The effect of the quadrupole moment of nitrogen on the phase behavior is considered through an extension of the SAFT-VR equation that includes an additional contribution to the Helmholtz free energy due to quadrupolar interactions. A significant improvement in the description of the phase diagram of the binary mixtures of nitrogen with different n-alkanes is obtained with the new approach when compared to predictions from the original SAFT-VR EOS (i.e., without the quadrupolar term). The experimental value for the quadrupole moment of nitrogen is used in the new equation; thus, no additional parameters are employed. Given the nonideal nature of the binary mixtures, a binary interaction parameter is needed to describe the full-phase diagram and high-pressure critical lines of these systems; however, this can be fitted to a single system and successfully used to predict the phase behavior of other binary mixtures without additional fitting. Furthermore, only a single, transferable, cross-energy parameter is required when the quadrupolar term is considered, whereas a cross-range parameter is also needed with the original SAFT-VR approach. The inclusion of the quadrupolar term in the equation of state therefore reduces the need to use effective parameters by explicitly including at the molecular level interactions due to the quadrupole moment.

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.