First-order mean-spherical approximation for interfacial phenomena: a unified method from bulk-phase equilibria study

The recently proposed first-order mean-spherical approximation (FMSA) [Y. Tang, J. Chem. Phys. 121, 10605 (2004)] for inhomogeneous fluids is extended to the study of interfacial phenomena. Computation is performed for the Lennard-Jones fluid, in which all phase equilibria properties and direct correlation function for density-functional theory are developed consistently and systematically from FMSA. Three functional methods, including fundamental measure theory for the repulsive force, local-density approximation, and square-gradient approximation, are applied in this interfacial investigation. Comparisons with the latest computer simulation data indicate that FMSA is satisfactory in predicting surface tension, density profile, as well as relevant phase equilibria. Furthermore, this work strongly suggests that FMSA is very capable of unifying homogeneous and inhomogeneous fluids, as well as those behaviors outside and inside the critical region within one framework.     

First-order mean spherical approximation for inhomogeneous fluids

The first-order mean-spherical approximation (FMSA) [Y. Tang, J. Chem. Phys., 118, 4140 (2003)] is extended to the studies of inhomogeneous fluids by combining with Rosenfeld's perturbative method [Y. Rosenfeld, J. Chem. Phys. 98, 8126 (1993)]. In the extension, the key input-direct correlation function of FMSA-is applied to constructing the free energy density functional. Preserving its high fidelity at the bulk limit, the FMSA shows satisfactory performance for Yukawa fluids near hard and attractive walls. The results are better than or comparable to several other theories reported before for the geometry. The FMSA is found, in particular, more satisfactory than the traditional mean-field theory for predicting density profiles around hard walls. The FMSA is also compared with the full MSA for inhomogeneous fluids, showing no appreciable differences. The inhomogeneous FMSA goes successfully through the self-consistency test for reproducing the radial distribution function of the bulk Yukawa fluid. As far as the computation is concerned, the FMSA can be executed much faster than any nonmean-field theories, and the speed is virtually identical to that of the mean-field theory.     

Prediction of phase behavior of nanoconfined Lennard-Jones fluids with density functional theory based on the first-order mean spherical approximation

The recently proposed first-order mean spherical approximation (FMSA) [Y. Tang, J. Chem. Phys. 121, 10605 (2004)] for inhomogeneous fluids is extended to study the phase behavior of nanoconfined Lennard-Jones fluids, which is consistent with the phase equilibria calculation of the corresponding bulk fluid. With a combination of fundamental measure theory, FMSA provides Helmholtz free energy and direct correlation function to formulate density functional theory, which implementation is as easy as the mean-field theory. Following previous success in predicting density profiles inside slit pores, this work is focused specially on the vapor-liquid equilibrium of the Lennard-Jones fluids inside these pores. It is found that outside the critical region FMSA predicts well the equilibrium diagram of slit pores with the sizes of 5.0, 7.5, and 10 molecular diameters by comparing with available computer simulation data. As a quantitative method, FMSA can be treated as an extension from its bulk calculation, while the mean-field theory is only qualitative, as its bulk version.     

由L-J流体的FMSA状态方程计算流体的PVT性质

运用Tang等提出的Lennard-Jones (L-J)流体两参数的一阶平均球形近似(FMSA)状态方程, 计算了流体的汽液共存相图和饱和蒸汽压曲线, 以及非饱和区的PVT性质, 并与文献数据进行比较. L-J参数由Tr＜0.95的汽液相共存数据回归得到. 计算结果表明, 对于分子较接近球形的流体, 除临界点附近外, 该方程可以在较大的温度和压力范围内计算真实流体的PVT性质, 结果满意. 对于球形分子, 该方程的精确度随分子尺寸的变大基本保持稳定. 该方程不适用于强极性物质. 在高密度区, 该方程的计算结果明显优于P-R方程. 对于分子偏离球形较远的流体, 该方程的适用性变差, 此时要考虑分子形状的影响, 可采用三参数的FMSA状态方程进行计算.     

Thermodynamic and structural properties of mixed colloids represented by a hard-core two-Yukawa mixture model fluid: Monte Carlo simulations and an analytical theory

The interaction between colloidal particles is well represented by a hard-core two-Yukawa potential. In order to assess the accuracy of theoretical predictions for the thermodynamic and structural properties of mixed colloids, standard Monte Carlo simulations are carried out for the hard-core two-Yukawa mixtures. In the simulations, one range parameter in the two-Yukawa potential is taken as 1.8 or 2.8647, and another is taken as 4, 8, or 13.5485. Both attractive and repulsive dominant cases of the potential outside the hard core are considered. The effects of temperature, density, composition, size and energy parameter ratios on internal energy, compressibility factor, and radial distribution function are investigated extensively. Theoretical calculations are performed in the framework of analytical solution for the Ornstein-Zernike equation with the first-order mean spherical approximation (FMSA). Our analysis shows that the FMSA is very accurate for the prediction of the compressibility factor of the hard-core two-Yukawa mixtures at all conditions studied. The FMSA generally predicts accurate internal energy, but overestimates the internal energy of the systems at lower temperatures. Furthermore, we found that a simplified exponential version of the FMSA predicts fairly good radial distribution function at contact for the mixed two-Yukawa fluids. The comparison of the theoretical compressibility factor with that from the Monte Carlo simulations suggests that the FMSA can be used to investigate the fluid-fluid equilibria of hard-core two-Yukawa mixtures.     

Direct correlation function for complex square barrier-square well potentials in the first-order mean spherical approximation

The direct correlation function of the complex discrete potential model fluids is obtained as a linear combination of the first-order mean spherical approximation (FMSA) solution for the simple square well model that has been reported recently [Hlushak et al., J. Chem. Phys. 130, 234511 (2009)]. The theory is employed to evaluate the structure and thermodynamics of complex fluids based on the square well-barrier and square well-barrier-well discrete potential models. Obtained results are compared with theoretical predictions of the hybrid mean spherical approximation, already reported in the literature [Guillen-Escamilla et al., J. Phys.: Condens. Matter 19, 086224 (2007)], and with computer simulation data of this study. The compressibility route to thermodynamics is then used to check whether the FMSA theory is able to predict multiple fluid-fluid transitions for the square barrier-well model fluids.     

Theoretical study of Sutherland fluids with long-range, short-range, and highly short-range potential parameters

On the basis of the first-order mean spherical approximation (FMSA) theory the behavior of Sutherland fluids with a number of parameters (gamma=3.1-36) is investigated. The investigation includes its modification by the simplified exponent approximation, renormalization group (RG) transformation, and density functional theory (DFT). For long-range parameters, the original FMSA is found sufficiently good to describe the global phase behavior, including inside the critical region. For short-range parameters, the modified FMSA by the simplified exponent approximation outside the critical region and RG transform inside the critical region are applied. For extremely short-range forces, the success is achieved by its combination with the DFT. This work gives a general sense about the capability of a theory for different ranges of potential, as well as for different thermodynamic states.     

Prediction of global vapor-liquid equilibria for mixtures containing polar and associating components with improved renormalization group theory

Our recently improved renormalization group (RG) theory is further reformulated within the context of density functional theory. To improve the theory for polar and associating fluids, an explicit and complete expression of the theory is derived in which the density fluctuation is expanded up to the third-order term instead of the original second-order term. A new predictive equation of state based on the first-order mean spherical approximation statistical associating fluid theory (FMSA-SAFT) and the newly improved RG theory is proposed for systems containing polar and associating fluids. The calculated results for both pure fluids and mixtures are in good agreement with experimental data both inside and outside the critical region. This work demonstrates that the RG theory incorporated with the solution of FMSA is a promising route for accurately describing the global phase behavior of complex fluids and mixtures.     

高阶FMSA展开研究缔合流体界面张力

结合一阶平均球近似(First-order mean-spherical approximation, FMSA)与重整化群(Renormalization group, RG)变换计算了流体全局性相行为. 应用FMSA理论解析得到的径向分布函数(Radial distribution function, RDF)和直接相关函数(Direct correction function, DCF)建立密度泛函方法, 并在其展开项中考虑了高阶微扰项作用, 即考虑了主体流体密度不一致性, 避免原有方法在计算密度分布时存在难以收敛、误差大等问题. 将高阶展开扩展应用到缔合流体, 计算表明, 和分子模拟数据相比, 界面密度分布和界面张力较之原有的密度泛函方法均有了明显改善.     

Structure and thermodynamics of hard-core Yukawa fluids: thermodynamic perturbation approaches

The thermodynamic perturbation theories, which are based on the power series of a coupling constant (λ-expansion), have been proposed for studying the structural and thermodynamic properties of a hard-core Yukawa (HCY) fluid: one (A1-approximation) is the perturbation theory based on the hard-sphere repulsion as a reference system. The other (A2-approximation) is the perturbation theory based on the reference system which incorporates both the repulsive and short-range attractive interactions. The first-order mean-spherical approximation (FMSA) provided by Tang and Lu [J. Chem. Phys. 99, 9828 (1993)] has been employed for investigating the thermodynamic properties of a HCY fluid using the alternative method via the direct correlation function. The calculated results show that (i) the A1 and A2 approximations are in excellent agreements with previous computer simulation results in the literature and compare with the semi-empirical works of Shukla including the higher-order free energy terms, (ii) the A1 and A2 approximations are better than the FMSA and the mean-spherical approximation, (iii) the A2-approximation compares with the A1-approximation, even though the perturbation effect of an A2-approximation is much smaller than that of an A1-approximation, and that (iv) the FMSA study is particularly of advantage in providing the structure and thermodynamics in a simple and analytic manner.     

Computer-assisted liquid chromatography for automated qualitative and quantitative analysis of toxic drugs

A computer-assisted liquid chromatographic system (MCASYST) is evaluated for automatic qualitative and quantitative analysis of toxic drugs in poisoned human fluids. The system performance is confirmed by the actual analyses of several cases of poisoning. The results indicate that the accuracy of identification by retention prediction and UV spectral search is very high, and quantitation by the memorized calibration curves can be performed very conveniently. However, it is difficult at present to identify all peaks that appear in chromatograms for human urine and serum because the drugs may be changed to their metabolites in those fluids. Metabolite identification will be the next step to improve the performance of this system.     

Theoretical prediction of multiple fluid-fluid transitions in monocomponent fluids

The authors use the analytical equation of state obtained by the discrete perturbation theory [A. L. Benavides and A. Gil-Villegas, Mol. Phys. 97, 1225 (1999)] to study the phase diagram of fluids with discrete spherical potentials formed by a repulsive square-shoulder plus an attractive square-well interaction (SS+SW). This interaction is characterized by the usual energy and size parameters plus three dimensionless parameters: two of them measuring the widths of the SS and the SW and the third the relative height of the SS. The matter of interest is that, for certain values of the interaction parameters, the SS+SW systems exhibit more than one first-order fluid-fluid transition. The evidence that several real substances (such as water, phosphorus, carbon, and silica, among others) exhibit an extra liquid-liquid transition has drawn interest into the study of interactions responsible for this behavior. The simple SS+SW fluid is one of the systems that, in spite of being spherically symmetric, shows multiple fluid-fluid transitions. In this work the authors investigate systematically the effect on the phase diagram of varying the interaction parameters. The use of an analytical free-energy equation gives a clear thermodynamic picture of the emergence of different types of critical points, throwing new light on the phase behavior of these fluids and thus clarifying previous results obtained by other techniques. The interplay of attractive and repulsive forces with several scale lengths produces very rich phase diagrams, including cases with three critical points. The region of the interaction-parameter space where multiple critical points appear is mapped for various families of interactions.     

Perturbed chain-statistical associating fluid theory extended to dipolar and quadrupolar molecular fluids

The perturbed chain statistical associating fluid theory (PC-SAFT) is extended to polar molecular fluids, namely dipolar and quadrupolar fluids. The extension is based on the perturbation theory for polar fluids by Stell and co-workers. Appropriate expressions are proposed for dipole-dipole, quadrupole-quadrupole, and dipole-quadrupole interactions. Furthermore, induced dipole interactions are calculated explicitly in the model. The new polar PC-SAFT model is relatively complex; for this purpose, a truncated polar PC-SAFT model is proposed using only the leading term in the polynomial expansion for polar interactions. The new model is used for the calculation of thermodynamic properties of various quadrupolar pure fluids. In all cases, the agreement between experimental data and model predictions is very good.     

Time evolution of density fluctuation in supercritical region. I. Non-hydrogen-bonded fluids studied by dynamic light scattering

The time evolution of the density fluctuation of molecules inhomogeneously dispersing in a mesoscopic volume is investigated by dynamic light scattering in several fluids in supercritical states. This study is the first time-domain investigation to compare the dynamics of density fluctuation among several fluids. The samples used are non-hydrogen-bonded fluids in the supercritical states: CHF(3), C(2)H(4), CO(2), and xenon. These four molecules have different properties but are of similar size. Under these conditions, the relationship between dynamic and static density inhomogeneities is studied by measuring the time correlation function of the density fluctuation. In all cases, this function is characterized by a single exponential function, decaying within a few microseconds. While the correlation times in the four fluids show noncoincidence, those values agree well with each other when scaled to a dimensionless parameter. From the results of this scaling based on the Kawasaki theory and Landau-Placzek theory, the relation between dynamics and static structures is analyzed, and the following four insights are obtained: (i) viscosity is the main contributor to the time evolution of density fluctuation; (ii) the principle of corresponding state is observed by the use of time-domain data; (iii) the Kawasaki theory and the Landau-Placzek theory are confirmed to be applicable to polar, nonpolar, and nondipolar fluids that have no hydrogen bonding, at temperatures relatively far from critical temperature; and (iv) the density fluctuation correlation length and the value of density fluctuation are estimated from the time-domain data and agree with the values from other experiments and calculations.     

Chromatography with supercritical fluids

Molecular interactions in the gas phase may be utilized for enhancing the volatility of substances in gas chromatography. A study of the effect of pressure on equilibria in GLC revealed that at higher pressures these interactions are quite appreciable, especially near the critical point. Measurements of plate heights show that with suitable adaptation of columns to the specific properties of high-pressure, dense fluids, the separating efficiency can be high and the speed of analysis fast. Experiments conducted with n-pentane and isopropanol as mobile fluids under supercritical conditions demonstrate that the volatility of heavy compounds may be enhanced by a factor of as much as 10⁴. With either a liquid or a solid adsorbent as a stationary phase, chromatography with supercritical fluids proves to be very suitable for the rapid analysis of heavy compounds. The effect of various operational parameters on the separation is discussed. Certain features of the technique arc illustrated with several examples.bl]     

Local thermodynamic derivation of Young's equation

A derivation of Young's equation based on the energy balance near the contact line is presented. Our proposal is rigorous and avoids the errors identified in the usual local derivation. It is valid under very general conditions (for any geometry, in a gravitational field and for compressive fluids). Deviations of the contact angle from Young's equation are discussed in several cases: surfaces of high curvature and line tension. Finally, the relationship between surface tensions and surface energies comes as an additional, natural result. Our derivation also provides a new physical insight into the equilibrium of forces acting near the contact line. Its local character makes the recourse to integral analysis unnecessary, which results in a great simplification when compared to other general treatments.     

Extending the GC-SAFT-VR approach to associating functional groups: Alcohols, aldehydes, amines and carboxylic acids

The statistical associating fluid theory is a widely used molecular-based equation of state that has been successfully applied to study a broad range of fluid systems. It provides a framework in which the effects of molecular shape and interactions on the thermodynamics and phase behavior of fluids can be separated and quantified. In the original approach, molecules were modeled as chains composed of identical segments; the heterogeneity of molecules in terms of structure and functional groups was described implicitly through effective parameters. To overcome this limitation, in recent works [Peng et al. Fluid Phase Equilib. 227(2), 131 (2009); Ind. Eng. Chem. Res. 49(3), 1378 (2010)] the GC-SAFT-VR approach has been developed to extend the theory to model chains composed of segments of different size and/or energy of interaction and enable the development of a group-contribution approach within the SAFT-VR framework in which molecular heterogeneity and connectivity is explicitly accounted for. The parameters for several key functional groups (CH₃, CH₂, CH, CH₂CH, CO, C₆H₅, esters, ethers, cis-alkenes and trans-alkenes groups) were determined by fitting to experimental vapor pressure and saturated liquid density data for a number of small molecules containing the functional groups of interest and transferability of the parameters tested by comparing the theoretical predictions with experimental data for pure fluids not included in the fitting process and binary mixtures of both simple fluids and the VLE and LLE of small molecules in polymer systems. In this work, we further extend the applicability of the GC-SAFT-VR approach through the study of the vapor-liquid phase behavior of associating systems, such as linear and branched alcohols, primary and secondary amines, aldehydes, and carboxylic acids, and their mixtures. In the study of these new molecules several new functional groups (OH (linear and branched), HCO, NH₂, NH and COOH) are defined and their molecular parameters characterized. The transferability of the parameters is again tested by comparing the theoretical predictions with experimental data for pure fluids and binary mixtures not included in the fitting process. The GC-SAFT-VR approach is found to predict the phase behavior of the systems studied in most cases in good agreement with experimental data and accurately captures the effects of changes in structure and molecular composition on phase behavior.     

Role of the first coordination shell in determining the equilibrium structure and dynamics of simple liquids

The traditional view that the physical properties of a simple liquid are determined primarily by its repulsive forces was recently challenged by Berthier and Tarjus, who showed that in some cases ignoring the attractions leads to large errors in the dynamics [L. Berthier and G. Tarjus, Phys. Rev. Lett. 103, 170601 (2009); J. Chem. Phys. 134, 214503 (2011)]. We present simulations of the standard Lennard-Jones liquid at several condensed-fluid state points, including a fairly low density state and a very high density state, as well as simulations of the Kob-Andersen binary Lennard-Jones mixture. By varying the range of the forces via a shifted-forces cutoff, results for the thermodynamics, dynamics, and structure show that the determining factor for getting the correct statics and dynamics is not whether or not the attractive forces per se are included in the simulations. What matters is whether or not interactions are included from all particles within the first coordination shell - the attractive forces can thus be ignored, but only at extremely high densities. The recognition of the importance of a local shell in condensed fluids goes back to van der Waals; our results confirm this idea and thereby the basic picture of the old hole and cell theories for simple condensed fluids.     

Shape factors in equations of state. Part II. Repulsion phenomena in multicomponent chain fluids

An equation of state for the multicomponent fluid phase of nonattracting rigid particles of arbitrary shape is presented. The equation is a generalization of a previously presented equation of state for pure fluids of rigid particles; the approach describes the volumetric properties of a pure fluid in terms of a shape factor, zeta, which can be back calculated by scaling the volumetric properties of pure fluids to that of a hard sphere. The performance of the proposed equation is tested against mixtures of chain fluids immersed in a "monomeric" solvent of hard spheres of equal and different sizes. Extensive new Monte Carlo simulation data are presented for 19 binary mixtures of hard homonuclear tangent freely-jointed hard sphere chains (pearl-necklace) of various lengths (three to five segments), with spheres of several size ratios and at various compositions. The performance of the proposed equation is compared to the hard-sphere SAFT approach and found to be of comparable accuracy. The equation proposed is further tested for mixtures of spheres with spherocylinders. In all cases, the equation proved to be accurate and simple to use.     

Pair structure of the hard-sphere Yukawa fluid: an improved analytic method versus simulations, Rogers-Young scheme, and experiment

We present a comprehensive study of the equilibrium pair structure in fluids of nonoverlapping spheres interacting by a repulsive Yukawa-like pair potential, with special focus on suspensions of charged colloidal particles. The accuracy of several integral equation schemes for the static structure factor, S(q), and radial distribution function, g(r), is investigated in comparison to computer simulation results and static light scattering data on charge-stabilized silica spheres. In particular, we show that an improved version of the so-called penetrating-background corrected rescaled mean spherical approximation (PB-RMSA) by Snook and Hayter [Langmuir 8, 2880 (1992)], referred to as the modified PB-RMSA (MPB-RMSA), gives pair structure functions which are in general in very good agreement with Monte Carlo simulations and results from the accurate but nonanalytical and therefore computationally more expensive Rogers-Young integral equation scheme. The MPB-RMSA preserves the analytic simplicity of the standard rescaled mean spherical (RMSA) solution. The combination of high accuracy and fast evaluation makes the MPB-RMSA ideally suited for extensive parameter scans and experimental data evaluation, and for providing the static input to dynamic theories. We discuss the results of extensive parameter scans probing the concentration scaling of the pair structure of strongly correlated Yukawa particles, and we determine the liquid-solid coexistence line using the Hansen-Verlet freezing rule.