An equation of state for electrolyte solutions by a combination of low-density expansion of non-primitive mean spherical approximation and statistical associating fluid theory

A two-parameter equation of state (EOS) for electrolyte solutions is developed. The equation is in terms of Helmholtz free energy and incorporated with our previous results of the low-density expansion of non-primitive mean spherical approximation (MSA). The concentration dependent dielectric constant is thus inherently included in the model. The statistical associating fluid theory (SAFT) is introduced to represent the association interactions, including the solvent–solvent and ion–solvent. The EOS is tested for 15 aqueous alkali halide solutions at ambient condition. The equation can represent simultaneously the mean ionic activity coefficients, the osmotic coefficients and densities in a good accuracy up to saturated concentration. The comparisons with EOSs published earlier in the literature are carried out. The limitations of the model are also discussed.     

Thermodynamic regularities for associating fluids from statistical associating fluid theory equation of state

In this work, we used a statistical associating fluid theory to analyze two important thermodynamic regularities for some associating fluids, including water, methanol and ethanol. The studied regularities included: (i) the common bulk modulus point on the isotherms of the reduced bulk modulus versus reduced density, (ii) near linearity of the reduced isothermal bulk modulus as a function of reduced pressure. In this work, we also reported the influence of the molecular size and interaction strength on the bulk modulus point.     

A study of a nonprimitive model for electrolyte solutions based on perturbation theory

The equation for the Helmholtz free energy for systems of small anisotropic molecules and ions is deduced by substituting the complete expression for various potential energies (including repulsive, dispersive, electrostatic, and induced energies) into the perturbation expansion. The equation is applied to pure water. The relative dielectric constant is set at unity. Based on the equal chemical potentials of equilibrated vapor and liquid phases, the molecular parameters of water are regressed from the densities of saturated vapor in the temperature range of 0 to 370°C. The ARD of regression is 1.16%. These parameters are used to predict the heat of vaporization and densities of saturated vapor and liquid phases of water in the same temperature range. The ARDs of prediction are 4.5% and 9.8%, respectively. The equation is used to correlate the osmotic coefficients of twelve 1:1 electrolyte solutions. The relative dielectric constant is set at unity. The parameters (Soft-sphere diameter and dispersive constant) of seven ions (Na⁺, K⁺, Rb⁺, Cs⁺, Cl^–, Br^–, and I^–) are obtained. The total average absolute deviation between calculated and experimental values of the osmotic coefficient is 0.041. The parameters of ions can keep constant in different systems.     

Two comments on the mean spherical approximation solution for the restricted primitive model of an electrolyte solution

Analytical radial distribution functions are proposed for the restricted primitive model of a 1:1 electrolyte solution which favourably compare with the Monte-Carlo results of Card and Valleau. Attention is also drawn to the electrical contribution to the activity coefficient in the mean spherical approximation.     

Equation of state for aqueous electrolyte systems based on the semirestricted non-primitive mean spherical approximation

The semirestricted non-primitive mean spherical approximation (npmsa) is used in combination with the PC-SAFT equation of state to model completely dissociating aqueous alkali halide systems. The salts are described using ion-specific parameters from tables and correlations. Upon analyzing aqueous electrolyte systems at infinite dilution of the salt it was concluded that for the arithmetic mean ion diameter of anion and cation, the semirestricted npmsa contribution gives no reliable results. Therefore, this parameter is adjusted in this work. The model was applied to aqueous alkali halide systems up to the solubility limit at T = 298.15 K. Mean ionic activity coefficients and osmotic coefficients were correlated with good results. The model was subsequently applied to temperatures up to T = 373.15 K and compared to liquid densities and to system pressures up to the solubility limit of the salts in water. The agreement between experimental data and the proposed equation of state is satisfactory for the liquid densities and excellent in case of the system pressures.     

Modeling of density of aqueous solutions of amino acids with the statistical associating fluid theory

The density of aqueous solutions of amino acids has been modeled with the statistical associating fluid theory (SAFT) equation of state. The modeling is accomplished by extending the previously developed new method to determine the SAFT parameters for amino acids. The modeled systems include α-alanine/H₂O, β-alanine/H₂O, proline/H₂O, l-asparagine/H₂O, l-glutamine/H₂O, l-histidine/H₂O, serine/H₂O, glycine/H₂O, alanine/H₂O/sucrose, dl-valine/H₂O/sucrose, arginine/H₂O/sucrose, serine/H₂O/ethylene glycol, and glycine/H₂O/ethylene glycol. The density of binary solutions of amino acids has been correlated or predicted with a high precision. And then the density of multicomponent aqueous solutions of amino acids has been modeled based on the modeling results of binary systems, and a high accuracy of density calculations has been obtained. Finally, the water activities of dl-valine/H₂O, glycine/H₂O, and proline/H₂O have been predicted without using binary interaction parameters, and good results have been obtained.     

On the mean spherical approximation for the Lennard–Jones fluid

Subtle details on the mean spherical approximation (MSA) theory for the Lennard–Jones potential [Fluid Phase Equilib. 134 (1997) 21] are presented. In order to enhance the appreciation of the theory, the accuracy of the mapping method and the contact approximation used in the theory are demonstrated by comparing with the exact results at certain extreme conditions. Technical derivation of internal energy of the MSA is also fully displayed. In addition, the typographic errors appeared in [Fluid Phase Equilib. 134 (1997) 21] are also corrected in this work.     

Microstructure and depletion forces in polymer-colloid mixtures from an interfacial statistical associating fluid theory

By using a classical density functional theory (interfacial statistical associating fluid theory), we investigate the structure and effective forces in nonadsorbing polymer-colloid mixtures. The theory is tested under a wide range of conditions and performs very well in comparison to simulation data. A comprehensive study is conducted characterizing the role of polymer concentration, particle/polymer-segment size ratio, and polymer chain length on the structure, polymer induced depletion forces, and the colloid-colloid osmotic second virial coefficient. The theory correctly captures a depletion layer on two different length scales, one on the order of the segment diameter (semidilute regime) and the other on the order of the polymer radius of gyration (dilute regime). The particle/polymer-segment size ratio is demonstrated to play a significant role on the polymer structure near the particle surface at low polymer concentrations, but this effect diminishes at higher polymer concentrations. Results for the polymer-mediated mean force between colloidal particles show that increasing the concentration of the polymer solution encourages particle-particle attraction, while decreasing the range of depletion attraction. At intermediate to high concentrations, depletion attraction can be coupled to a midrange repulsion, especially for colloids in solutions of short chains. Colloid-colloid second virial coefficient calculations indicate that the net repulsion between colloids at low polymer densities gives way to net attraction at higher densities, in agreement with available simulation data. Furthermore, the results indicate a higher tendency toward colloidal aggregation for larger colloids in solutions of longer chains.     

On a new approximation in the theory of electrolyte solutions

We compare the results of a new approximation for the interionic radial distribution function developed previously by Olivares and McQuarrie to those from the nonlinear Poisson—Boltzmann equation for the highly-charged system of the spherical protein bearing a charge as 20 in an aqueous electrolyte solution. In particular, we use both approximations to predict the experimental results for protein titration, which are the experimental data to which the results of the Poisson—Boltzmann equation had been directed.     

Semiclassical approach to model quantum fluids using the statistical associating fluid theory for systems with potentials of variable range

Thermodynamic properties of quantum fluids are described using an extended version of the statistical associating fluid theory for potentials of variable range (SAFT-VR) that takes into account quantum corrections to the Helmholtz free energy A, based on the Wentzel-Kramers-Brillouin approximation. We present the theoretical background of this approach (SAFT-VRQ), considering two different cases depending on the continuous or discontinuous nature of the particles pair interaction. For the case of continuous potentials, we demonstrate that the standard Wigner-Kirkwood theory for quantum fluids can be derived from the de Broglie-Bohm formalism for quantum mechanics that can be incorporated within the Barker and Henderson perturbation theory for liquids in a straightforward way. When the particles interact via a discontinuous pair potential, the SAFT-VR method can be combined with the perturbation theory developed by Singh and Sinha [J. Chem. Phys. 67, 3645 (1977); and ibid. 68, 562 (1978)]. We present an analytical expression for the first-order quantum perturbation term for a square-well potential, and the theory is applied to model thermodynamic properties of hydrogen, deuterium, neon, and helium-4. Vapor-liquid equilibrium, liquid and vapor densities, isochoric and isobaric heat capacities, Joule-Thomson coefficients and inversion curves are predicted accurately with respect to experimental data. We find that quantum corrections are important for the global behavior of properties of these fluids and not only for the low-temperature regime. Predictions obtained for hydrogen compare very favorably with respect to cubic equations of state.     

Toward making the mean spherical approximation of primitive model electrolytes analytic: an analytic approximation of the MSA screening parameter

The mean spherical approximation (MSA) for the primitive model of electrolytes provides reasonable estimates of thermodynamic quantities such as the excess chemical potential and screening length. It is especially widely used because of its explicit formulas so that numerically solving equations is minimized. As originally formulated, the MSA screening parameter Γ (akin to the reciprocal of the Debye screening length) does not have an explicit analytic formula; an equation for Γ must be solved numerically. Here, an analytic approximation for Γ is presented whose relative error is generally ?10(-5). If more accuracy is desired, one step of an iterative procedure (which also produces an explicit formula for Γ) is shown to give relative errors within machine precision in many cases. Even when ion diameter ratios are ～10 and ion valences are ～10, the relative error for the analytic approximation is still ?10(-3) and for the single iterative substitution it is ?10(-9).     

Modelling of phase equilibria for associating mixtures using an equation of state

In the present work, the group contribution with association equation of state (GCA-EoS) is extended to represent phase equilibria in mixtures containing acids, esters, and ketones, with water, alcohols, and any number of inert components. Association effects are represented by a group-contribution approach. Self- and cross-association between the associating groups present in these mixtures are considered. The GCA-EoS model is compared to the group-contribution method MHV2, which does not take into account explicitly association effects. The results obtained with the GCA-EoS model are, in general, more accurate when compared to the ones achieved by the MHV2 equation with less number of parameters. Model predictions are presented for binary self- and cross-associating mixtures.     

Electrolyte-UNIQUAC-NRF model for the correlation of the mean activity coefficient of electrolyte solutions

The new electrolyte-UNIQUAC-NRF excess Gibbs function is obtained for calculation of the activity coefficient of the binary electrolyte solutions. The excess Gibbs energy of the model consists of the Pitzer–Debye–Hückel equation, describing the long-range electrostatic contribution and the electrolyte-UNIQUAC-NRF model to account for the short-range contributions. With two adjustable parameters per electrolyte, the new model is applied to correlation of the mean activity coefficients of more than 130 binary aqueous electrolyte solutions at 25 °C. Also the binary parameters, obtaining from regression of mean activity data, are used for prediction of osmotic coefficient data for the same electrolytes. The results are compared with various local composition models such as the electrolyte-NRTL, electrolyte NRF-Wilson, electrolyte-NRTL-NRF, N-Wilson-NRF models. The comparison of the results with experiment demonstrates that the new model can correlate the experimental activity coefficient data and predict the osmotic coefficient data of binary electrolytes accurately.     

Solution of the mean spherical approximation for polydisperse multi-Yukawa hard-sphere fluid mixture using orthogonal polynomial expansions

The Blum-Hoye [J. Stat. Phys. 19 317 (1978)] solution of the mean spherical approximation for a multicomponent multi-Yukawa hard-sphere fluid is extended to a polydisperse multi-Yukawa hard-sphere fluid. Our extension is based on the application of the orthogonal polynomial expansion method of Lado [Phys. Rev. E 54, 4411 (1996)]. Closed form analytical expressions for the structural and thermodynamic properties of the model are presented. They are given in terms of the parameters that follow directly from the solution. By way of illustration the method of solution is applied to describe the thermodynamic properties of the one- and two-Yukawa versions of the model.     

Lanthanide salts solutions: representation of osmotic coefficients within the binding mean spherical approximation

Osmotic coefficients of aqueous solutions of lanthanide salts are described using the binding mean spherical approximation (BIMSA) model based on the Wertheim formalism for association. The lanthanide(III) cation and the co-ion are allowed to form a 1-1 ion pair. Hydration is taken into account by introducing concentration-dependent cation size and solution permittivity. An expression for the osmotic coefficient, derived within the BIMSA, is used to fit data for a wide variety of lanthanide pure salt aqueous solutions at 25 degrees C. A total of 38 lanthanide salts have been treated, including perchlorates, nitrates, and chlorides. For most solutions, good fits could be obtained up to high ionic strengths. The relevance of the fitted parameters has been discussed, and a comparison with literature values has been made (especially the association constants) when available.     

Investigation of the surface tension of methane and n-alkane mixtures by perturbed-chain statistical associating fluid theory combined with density-gradient theory

The perturbed-chain statistical associating fluid theory (PC-SAFT) and density-gradient theory are used to construct an equation of state to describe the phase behavior of binary methane–n-alkane mixtures. With the molecular parameters and influence parameters regressed from bulk properties and surface tensions of pure fluids, respectively as input, both the bulk and interfacial properties are investigated. The surface tension of the binary systems methane–propane, methane–pentane, methane–heptane and methane–decane are predicted, and the results are satisfactory compared with the experimental data. Our results show that PC-SAFT combined with density-gradient theory is able to describe the interfacial properties of binary methane–n-alkane mixtures in wide temperature and pressure ranges, and illustrate the influence of the equilibrium bulk properties and chain length of n-alkane molecule on the interfacial properties.     

Analytic solution of the mean spherical approximation for a multi-component plasma

The mean spherical approximation for a mixture of charged hard spheres in a uniform neutralizing background is solved analytically. The factor correlation functions and the excess thermodynamic properties are explicitly expressed through a single parameter, which can be obtained by solving an algebraic equation.