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.     

Application of a renormalization-group treatment to the statistical associating fluid theory for potentials of variable range (SAFT-VR)

An accurate prediction of phase behavior at conditions far and close to criticality cannot be accomplished by mean-field based theories that do not incorporate long-range density fluctuations. A treatment based on renormalization-group (RG) theory as developed by White and co-workers has proven to be very successful in improving the predictions of the critical region with different equations of state. The basis of the method is an iterative procedure to account for contributions to the free energy of density fluctuations of increasing wavelengths. The RG method has been combined with a number of versions of the statistical associating fluid theory (SAFT), by implementing White's earliest ideas with the improvements of Prausnitz and co-workers. Typically, this treatment involves two adjustable parameters: a cutoff wavelength L for density fluctuations and an average gradient of the wavelet function Φ. In this work, the SAFT-VR (variable range) equation of state is extended with a similar crossover treatment which, however, follows closely the most recent improvements introduced by White. The interpretation of White's latter developments allows us to establish a straightforward method which enables Φ to be evaluated; only the cutoff wavelength L then needs to be adjusted. The approach used here begins with an initial free energy incorporating only contributions from short-wavelength fluctuations, which are treated locally. The contribution from long-wavelength fluctuations is incorporated through an iterative procedure based on attractive interactions which incorporate the structure of the fluid following the ideas of perturbation theories and using a mapping that allows integration of the radial distribution function. Good agreement close and far from the critical region is obtained using a unique fitted parameter L that can be easily related to the range of the potential. In this way the thermodynamic properties of a square-well (SW) fluid are given by the same number of independent intermolecular model parameters as in the classical equation. Far from the critical region the approach provides the correct limiting behavior reducing to the classical equation (SAFT-VR). In the critical region the β critical exponent is calculated and is found to take values close to the universal value. In SAFT-VR the free energy of an associating chain fluid is obtained following the thermodynamic perturbation theory of Wertheim from the knowledge of the free energy and radial distribution function of a reference monomer fluid. By determining L for SW fluids of varying well width a unique equation of state is obtained for chain and associating systems without further adjustment of critical parameters. We use computer simulation data of the phase behavior of chain and associating SW fluids to test the accuracy of the new equation.     

Thermodynamic properties of the Mie n-6 fluid: a comparison between statistical associating fluid theory of variable range approach and molecular dynamics results

Molecular dynamics (MD) simulations of direct and derivative thermodynamic properties of the Mie n-6 fluid (n=8, 10, and 12) have been performed for liquid to supercritical states. Using the results, an in depth test of the monomer-monomer interaction estimation of a recently derived statistical associating fluid theory of variable range (SAFT-VR) equation of state [Lafitte et al., J. Chem. Phys., 124, 024509 (2006)] has been carried out based on the Mie n-6 potential. For pure fluids, using an appropriate scaling, MD simulations show that density and isometric heat capacity are nearly independent of n, whereas sound velocity and thermal pressure coefficient tend to increase with n. In addition, the results show that predictions provided by the equation of state are consistent with those coming from MD and catch correctly the trends of each property with n except for the heat capacity. The comparison is next extended to binary mixtures with components differing only in the value of the n parameter and which demonstrate the reliability of the scheme (MX1b) used by Lafitte et al. to deal with this parameter in the SAFT-VR equation of state. In addition, a new empirical one-fluid approximation of the n parameter is proposed thanks to MD simulations, which very favorably compare with the one-fluid model on n previously proposed in the literature. The consistency of this approximation is addressed by making use of it in combination with the SAFT-VR Mie equation of state. It is shown that using such an approach, which is easier to handle than the MX1b one, yields slightly improved results compared to those of the MX1b.     

Predicting adsorption isotherms using a two-dimensional statistical associating fluid theory

A molecular thermodynamics approach is developed in order to describe the adsorption of fluids on solid surfaces. The new theory is based on the statistical associating fluid theory for potentials of variable range [A. Gil-Villegas et al., J. Chem. Phys. 106, 4168 (1997)] and uses a quasi-two-dimensional approximation to describe the properties of adsorbed fluids. The theory is tested against Gibbs ensemble Monte Carlo simulations and excellent agreement with the theoretical predictions is achieved. Additionally the authors use the new approach to describe the adsorption isotherms for nitrogen and methane on dry activated carbon.     

Modified interfacial statistical associating fluid theory: application to tethered polymer chains

Modified interfacial statistical associating fluid theory density functional theory is extended to tethered polymer chains in the absence or presence of free polymer chains. The structures of the "dry" and "wet" polymer brushes have been calculated and compared with simulation results available in the literature. The comparisons show that the theory accurately predicts the structure of the tethered polymer brush. The average brush heights calculated from the theory agree with well-established scaling theories for tethered polymers. However, these scaling theories cannot predict the detailed structure, accurately. The effects of the segment-segment interactions of the tethered polymer and the free polymer have been effectively captured by the theory.     

Modeling thermodynamic properties of natural gas mixtures using perturbed-chain statistical associating fluid theory

Perturbed-Chain Statistical Associating Fluid Theory (PC-SAFT) was employed for predicting thermodynamic properties of natural gas mixture. Thermodynamic properties like density, isobaric and isochoric heat capacity, enthalpy, entropy, and internal energy were calculated with the PC-SAFT. Results are validated against experimental data for natural gas and mixtures similar to natural gas. The validation show that the Average Absolute Deviation (AAD) for density is 1.10% for binary mixture and 1.08% for mixtures similar to natural gas. Also AAD value for enthalpy is 1.42%, for internal energy, 0.77, for entropy, 0.43, for isochoric heat capacity, 1.26%, and for isobaric heat capacity, 2.66%. Results show PC-SAFT to be able to predict all the thermodynamics properties of natural gas and mixtures similar to natural gas with high accuracy in a wide range of temperature and pressure.     

A study of square-well statistical associating fluid theory approximations

Six square-well (SW) statistical associating fluid theory (SAFT) models, fitted to the experimental saturated liquid volume and saturated vapor pressure for pure n-alkanes, are analyzed for predicting the coexisting densities, second virial coefficients, and binary phase equilibria. The models that result in low values of the segment energy and weak molecular weight dependence of the parameters are found to be more accurate for real fluids. The inclusion of the dimer structure in the SW chain term seems to produce no significant benefit for representing real substances.     

Modified interfacial statistical associating fluid theory: a perturbation density functional theory for inhomogeneous complex fluids

A density functional theory based on Wertheim's first order perturbation theory is developed for inhomogeneous complex fluids. The theory is derived along similar lines as interfacial statistical associating fluid theory [S. Tripathi and W. G. Chapman, J. Chem. Phys. 122, 094506 (2005)]. However, the derivation is more general and applies broadly to a range of systems, retaining the simplicity of a segment density based theory. Furthermore, the theory gives the exact density profile for ideal chains in an external field. The general avail of the theory has been demonstrated by applying the theory to lipids near surfaces, lipid bilayers, and copolymer thin films. The theoretical results show excellent agreement with the results from molecular simulations.     

Thermodynamic modeling of CO2 solubility in ionic liquid with heterosegmented statistical associating fluid theory

Heterosegmented statistical associating fluid theory is used to represent the CO₂ solubility in ionic liquids. As in our previous work, ionic liquid molecule is divided into several groups representing the alkyls, cation head, and anion. The cation of ionic liquid is modeled as a chain molecule that consists of one spherical segment representing the cation head and groups of segments of different types representing different substituents (alkyls). The anion of ionic liquid is modeled as a spherical segment of different type. To account for the electrostatic/polar interaction between the cation and anion, the spherical segments representing cation head and anion each have one association site, which can only cross associate. Carbon dioxide is modeled as a molecule with three association sites, two sites of type O and one site of type C, where sites of the same type do not associate with each other. The parameters of CO₂ are obtained from the fitting of the density and the saturation vapor pressure of CO₂. For the CO₂-ionic liquid systems, cross association between site of type C in CO₂ and another association site in anion is allowed to occur to account for the Lewis acid–base interaction. The parameters for cross association interactions and the binary interaction parameters used to adjust the dispersive interactions between unlike segments are obtained from the fitting of the available CO₂ solubility in ionic liquids. The model is found to well represent the CO₂ solubility in the imidazolium ionic liquids from 283 to 415 K and up to 200 bar.     

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.     

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.     

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.     

Molecular theory of thermal conductivity of the Lennard-Jones fluid

In this paper the thermal conductivity of the Lennard-Jones fluid is calculated by applying the combination of the density-fluctuation theory, the modified free volume theory of diffusion, and the generic van der Waals equation of state. A Monte Carlo simulation method is used to compute the equilibrium pair-correlation function necessary for computing the mean free volume and the coefficient in the potential-energy and virial contributions to the thermal conductivity. The theoretical results are compared with our own molecular dynamics simulation results and with those reported in the literature. They agree in good accuracy over wide ranges of density and temperature examined in molecular dynamics simulations. Thus the combined theory represents a molecular theory of thermal conductivity of the Lennard-Jones fluid and by extension simple fluids, which enables us to compute the nonequilibrium quantity by means of the Monte Carlo simulations for the equilibrium pair-correlation function.     

Thermodynamic properties and aggregate formation of surfactant-like molecules from theory and simulation

The goal of this work is twofold: to predict the phase equilibria behavior of simplified surfactant models and to predict the population of aggregates as a function of pressure. We compare Monte Carlo simulation results of these systems with predictions from a modified version of the statistical associating fluid theory (soft-SAFT). Surfactant-like molecules are modeled as Lennard-Jones chains of tangent segments with one or two association sites. We study the influence of the number and location of the association sites on the thermodynamic properties and fraction of nonbonded molecules in all cases. The influence of the chain length is also investigated for a particular location of the sites. Results are compared with NPT Monte Carlo simulations to test the accuracy of the theory, and to study the molecular configurations of the system. Soft-SAFT is able to quantitatively predict the MC PVT results, independently of the location of the association sites. The theory is also able to capture the qualitative trend of the population of aggregates with pressure. Quantitative agreement is only obtained for specific locations of the sites.     

Phase behavior of a fluid confined in slitlike pores with walls modified by preadsorbed chain molecules

We have studied the microscopic structure, thermodynamics of adsorption, and phase behavior of Lennard-Jones fluid in slitlike pores with walls modified due to preadsorption of chain molecules. The chain species are grafted at the walls by terminating segments. Our theoretical considerations are based on a density functional approach in the semigrand canonical ensemble. The applied constraint refers to the constant number of grafted chain molecules in the pore without restriction of the number of chains at each of the walls. We have observed capillary condensation of Lennard-Jones fluid combined with the change of the distribution of chains from nonsymmetric to symmetric with respect to the pore walls. The phase diagrams of the model are analyzed in detail, dependent on the pore width, length of chains, and grafted density.