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.     

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.     

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 properties of Lennard-Jones chain molecules: renormalization-group corrections to a modified statistical associating fluid theory

A modified version of the statistical associating fluid theory (SAFT), the so-called soft-SAFT equation of state (EOS), has been extended by a crossover treatment to take into account the long density fluctuations encountered when the critical region is approached. The procedure, based on White's work from the renormalization group theory [Fluid Phase Equilibria 75, 53 (1992); L. W. Salvino and J. A. White, J. Chem. Phys. 96, 4559 (1992)], is implemented in terms of recursion relations where the density fluctuations are successively incorporated. The crossover soft-SAFT equation provides the correct nonclassical critical exponents when approaching the critical point, and reduces to the original soft-SAFT equation far from the critical region. The accuracy of the global equation is tested by direct comparison with molecular simulation results of Lennard-Jones chains, obtaining very good agreement and clear improvements compared to the original soft-SAFT EOS. Excellent agreement with vapor-liquid equilibrium experimental data inside and outside the critical region for the n-alkane series is also obtained. We provide a set of transferable molecular parameters for this family, unique for the whole range of thermodynamic properties.     

Re-examining the properties of the aqueous vapor-liquid interface using dispersion corrected density functional theory

First-principles molecular dynamics simulations, in which the forces are computed from electronic structure calculations, have great potential to provide unique insight into structure, dynamics, electronic properties, and chemistry of interfacial systems that is not available from empirical force fields. The majority of current first-principles simulations are driven by forces derived from density functional theory with generalized gradient approximations to the exchange-correlation energy, which do not capture dispersion interactions. We have carried out first-principles molecular dynamics simulations of air-water interfaces employing a particular generalized gradient approximation to the exchange-correlation functional (BLYP), with and without empirical dispersion corrections. We assess the utility of the dispersion corrections by comparison of a variety of structural, dynamic, and thermodynamic properties of bulk and interfacial water with experimental data, as well as other first-principles and force field-based simulations.     

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.     

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.     

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.     

Wetting behavior of spherical nanoparticles at a vapor-liquid interface: a density functional theory study

The wetting behavior of spherical nanoparticles at a vapor-liquid interface is investigated by using density functional theory, and the line tension calculation method is modified by analyzing the total energy of the vapor-liquid-particle equilibrium. Compared with the direct measurement data from simulation, the results reveal that the thermodynamically consistent Young's equation for planar interfaces is still applicable for high curvature surfaces in predicting a wide range of contact angles. The effect of the line tension on the contact angle is further explored, showing that the contact angles given by the original and modified Young's equations are nearly the same within the region of 60° < θ < 120°. Whereas the effect is considerable when the contact angle deviates from the region. The wetting property of nanoparticles in terms of the fluid-particle interaction strength, particle size, and temperature is also discussed. It is found that, for a certain particle, a moderate fluid-particle interaction strength would keep the particle stable at the interface in a wide temperature range.     

Double-hybrid density functional theory for excited electronic states of molecules

Double-hybrid density functionals are based on a mixing of standard generalized gradient approximations (GGAs) for exchange and correlation with Hartree-Fock (HF) exchange and a perturbative second-order correlation part (PT2) that is obtained from the Kohn-Sham (GGA) orbitals and eigenvalues. This virtual orbital-dependent functional (dubbed B2PLYP) contains only two empirical parameters that describe the mixture of HF and GGA exchange (ax) and of the PT2 and GGA correlation (ac), respectively. Extensive testing has recently demonstrated the outstanding accuracy of this approach for various ground state problems in general chemistry applications. The method is extended here without any further empirical adjustments to electronically excited states in the framework of time-dependent density functional theory (TD-DFT) or the closely related Tamm-Dancoff approximation (TDA-DFT). In complete analogy to the ground state treatment, a scaled second-order perturbation correction to configuration interaction with singles (CIS(D)) wave functions developed some years ago by Head-Gordon et al. [Chem. Phys. Lett. 219, 21 (1994)] is computed on the basis of density functional data and added to the TD(A)-DFTGGA excitation energy. The method is implemented by applying the resolution of the identity approximation and the efficiency of the code is discussed. Extensive tests for a wide variety of molecules and excited states (of singlet, triplet, and doublet multiplicities) including electronic spectra are presented. In general, rather accurate excitation energies (deviations from reference data typically <0.2 eV) are obtained that are mostly better than those from standard functionals. Still, systematic errors are obtained for Rydberg (too low on average by about 0.3 eV) and charge-transfer transitions but due to the relatively large ax parameter (0.53), B2PLYP outperforms most other functionals in this respect. Compared to conventional HF-based CIS(D), the method is more robust in electronically complex situations due to the implicit account of static correlation effects by the GGA parts. The (D) correction often works in the right direction and compensates for the overestimation of the transition energy at the TD level due to the elevated fraction of HF exchange in the hybrid GGA part. Finally, the limitations of the method are discussed for challenging systems such as transition metal complexes, cyanine dyes, and multireference cases.     

Analysis of the validity of perturbation density functional theory: based on extensive simulation for simple fluid at supercritical and subcritical temperature under various external potentials

Because of the scarcity of available simulation data for confined hard-core attractive Yukawa model fluid, extensive Monte Carlo (MC) simulation research for this fluid under the influence of various external potentials were carried out. The present MC simulation results were employed to test a performance of the third-order perturbation density functional theory (DFT) based on a high order direct correlation function (DCF) [S. Zhou and E. Ruckenstein, Phys. Rev. E. 61, 2704 (2000)]. It was found that the present perturbation DFT formalism is soundly structured only if the imported second-order DCF is reliable. In this case, the accuracy of the results can be satisfactory or even very high for various types of external potentials. Further, the associated adjustable parameter can be universal, i.e., independent of the particular external field responsible for the generation of a nonuniform density profile. Dependence of both the maintenance of the reliability of the formalism and holding of the universality of the adjustable parameter on the accuracy of the imported bulk second-order DCF can be strengthened by the large difference between the external field investigated and that caused by a single hard wall used for specification of the adjustable parameter. In case the gaseous density in the subcritical region is below the coexistence density, an excellent performance of the present formalism is observed even for the mean spherical approximation's second-order DCF as an input. This advantageous property, combined with the fact that the present formalism needs only the second-order DCF of fluid at the coexistence state as an input, enables the present formalism to be a very good theoretical tool for the investigations of wetting and prewetting transitions.     

The role of fluid wall association on adsorption of chain molecules at functionalized surfaces: a density functional approach

We present a density functional theory to describe adsorption in systems where selected segments of chain molecules of fluids can bond (or associate) with functional groups attached to the surfaces. Association of active segments with the surface is modeled within the framework of the first-order thermodynamic perturbation theory. We discuss the influence of several parameters such as the density of surface active sites, the energy of association, the chain length, and the number of the active segment in the chain molecule on the structure of the fluid adjacent to the wall. The proposed model can be considered as a first step towards developing a density functional theory of molecular brushes chemically bonded to solid surfaces.     

Equation of state for square-well chain molecules with variable range II. Extension to mixtures

An equation of state (EOS) developed in our previous work for square-well chain molecules with variable range is further extended to the mixtures of non-associating fluids. The volumetric properties of binary mixtures for small molecules as well as polymer blends can well be predicted without using adjustable parameter. With one temperature-independent binary interaction parameter, satisfactory correlations for experimental vapor–liquid equilibria (VLE) data of binary normal fluid mixtures at low and elevated pressures are obtained. In addition, VLE of n-alkane mixtures and nitrogen + n-alkane mixtures at high pressures are well predicted using this EOS. The phase behavior calculations on polymer mixture solutions are also investigated using one-fluid mixing rule. The equilibrium pressure and solubility of gas in polymer are evaluated with a single adjustable parameter and good results are obtained. The calculated results for gas + polymer systems are compared with those from other equations of state.     

基于化学缔合统计理论的链状流体状态方程

基于化学缔合统计理论的链状流体状态方程(EOS)能够反映实际分子的形状、链节成链、缔合等具体信息,在实际流体热力学性质计算中有着广泛应用.一般的链状流体EOS仅考虑相邻链节间的相关性,我们则借助统计力学和计算机模拟结果在模型中纳入了相间链节间的相关性,获得的硬球链流体(HSCF)模型能够更好地预测模型流体的压缩因子和第二维里系数.以HSCF为参考,引入方阱色散微扰项获得了实际方阱链流体(SWCF)EOS;结合根据黏滞球模型导得的缔合项,进一步构建了缔合流体EOS.最近,我们根据微扰理论和积分方程方法又开发了一新的变阱宽方阱链流体(SWCF-VR)模型.SWCF和SWCF-VREOSs可很好地用于计算小分子、聚合物、离子液体等纯流体及混合物的相行为、热焓、表面张力、黏度等热力学及传递性质,显示了模型良好的工程应用价值.本文就本课题组多年来在自由空间范畴内基于化学缔合统计理论开发链状流体EOS及其实际应用作系统的总结.     

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.