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.     

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.     

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.     

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.     

First-order mean spherical approximation for attractive, repulsive, and multi-Yukawa potentials

The first-order mean spherical approximation (FMSA) theory proposed by Tang et al. [Fluid Phase Equilib., 134, 21(1997)] is applied for studying several typical Yukawa fluids, including attractive, repulsive, and multi-Yukawa cases. The FMSA study is particularly advantageous in providing thermodynamics and structure information in a simple, analytical, and consistent manner. Comparisons with the latest reported computer simulation data for compressibility factor, internal energy, and radial distribution function show that FMSA performs very well and the performance is very close to the full MSA and to several other theories, developed individually for the above-mentioned cases or properties. The present study provides solid evidence to support FMSA applications to more complex fluids.     

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

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

受限于纳米狭缝中的四缔合Lennard-Lones流体的相平衡和界面张力研究

分别用改进的基础测量理论和平均球近似理论表达短程作用和长程作用对四缔合Lennard-Jones流体的过剩自由能的贡献. 在密度函泛理论的框架下, 研究了平均密度等温线, 密度分布, 未缔合分子在平衡汽相和液相中的分布, 相平衡以及平衡时的界面张力等热力学性质. 分析了缔合能量, 流体-固体作用和孔宽对受限于纳米狭缝中的四缔合Lennard-Jones流体相行为的影响.     

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.     

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

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

The QCHB model of fluids and their mixtures

The essentials of the QCHB (quasi-chemical hydrogen-bonding) equation-of-state model are presented along with some applications for calculations of phase equilibria and interfacial properties of fluids and their mixtures. This is a model applicable to non-polar systems as well as to highly non-ideal systems with strong specific interactions, to systems of small molecules as well as to macromolecules, including polydisperse polymers, glasses, and gels, to liquids as well as to vapours including supercritical systems, to homogeneous as well as to inhomogeneous systems. A quasi-thermodynamic approach of inhomogeneous systems is used for modeling the fluid–fluid interface. Consistent expressions for the interfacial tension and interfacial profiles for various properties are presented. A satisfactory agreement is obtained between experimental and calculated surface tensions. Extension of the approach to mixtures is examined along with the associated problems for the numerical calculations of the interfacial profiles. A new equation is derived for the chemical potentials in the interfacial region, which facilitates very much the calculation of the composition profiles across the interface. The relation of the model with the COSMO-RS approach is also discussed.     

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.     

Criticality of a liquid-vapor interface from an inhomogeneous integral equation theory

A microscopic theory is developed to study the liquid-vapor interfacial properties of simple fluids with ab initio treatment of the inhomogeneous two-body correlation functions, without any interpolation. It consists of the inhomogeneous Ornstein-Zernike equation coupled with the Duh-Henderson-Verlet closure and the Lovett-Mou-Buff-Wertheim equation. For the liquid-vapor interface of the Lennard-Jones fluid, we obtained the density profile and the surface tension, as well as their critical behaviour. In particular, we identified non-classical critical exponents. The theory accurately predicts the phase diagram and the interfacial properties in a very good agreement with simulations. We also showed that the method leads to true capillary-wave asymptotics in the macroscopic limit.     

Phase equilibria and plate-fluid interfacial tensions for associating hard sphere fluids confined in slit pores

The excess Helmholtz free energy functional for associating hard sphere fluid is formulated by using a modified fundamental measure theory [Y. X. Yu and J. Z. Wu, J. Chem. Phys. 117, 10156 (2002)]. Within the framework of density functional theory, the thermodynamic properties including phase equilibria for both molecules and monomers, equilibrium plate-fluid interfacial tensions and isotherms of excess adsorption, average molecule density, average monomer density, and plate-fluid interfacial tension for four-site associating hard sphere fluids confined in slit pores are investigated. The phase equilibria inside the hard slit pores and attractive slit pores are determined according to the requirement that temperature, chemical potential, and grand potential in coexistence phases should be equal and the plate-fluid interfacial tensions at equilibrium states are predicted consequently. The influences of association energy, fluid-solid interaction, and pore width on phase equilibria and equilibrium plate-fluid interfacial tensions are discussed.     

纯流体在临界区内外的表面张力研究

An equation of state(EOS)applicable for both the uniform and non-uniform fluids was established by using thedensity-gradient expansion,in which the influence parameter к[p(r),T] was obtained by the use of direct correlationfunction.The density functional theory(DFT)provides a framework under which both the phase equilibria and in-terfacial properties can be investigated within a single set of molecular parameters.The phase equilibria inside thecritical region can be improved by the renormalization group theory(RGT).However,the correction of interracialproperties by DFT and RGT is computationally difficult.In the present work,the density gradient theory(DGT)inwhich к[p(r),T] is treated as a constant is used to combine with the RGT for interfacial properties inside the criticalregion.     

Investigation of excess adsorption, solvation force, and plate-fluid interfacial tension for Lennard-Jones fluid confined in slit pores

The excess Helmholtz free energy functional is formulated in terms of a modified fundamental measure theory [Y. X. Yu and J. Z. Wu, J. Chem. Phys. 117, 10156 (2002)] for a short ranged repulsion and a first-order mean-spherical approximation theory [Y. P. Tang, J. Chem. Phys. 118, 4140 (2003)] for a long ranged attraction. Within the framework of the density functional theory, the density profile, excess adsorption, solvation force, and plate-fluid interfacial tension of a Lennard-Jones fluid confined in slit pores are predicted, and the results agree well with the simulation data. The phase equilibria inside the slit pores are determined according to the requirement that temperature, chemical potential, and grand potential in coexistence phases should be equal, and the plate-fluid interfacial tensions at equilibrium states are predicted consequently.     

Interfacial colloidal sedimentation equilibrium. I. Intensity based confocal microscopy

This paper reports confocal microscopy measurements of inhomogeneous colloidal sedimentation equilibrium profiles near planar wall surfaces for conditions when colloid dimensions are comparable to the characteristic gravitational length scale. The intensity based confocal method developed in this work enables real-space measurements of one-dimensional density profiles of Brownian colloids without identifying many single colloid centers in large imaging volumes. Measured sedimentation equilibrium profiles for single-phase interfacial fluids and for coexisting inhomogeneous fluid and solid phases are in agreement with a perturbation theory and Monte Carlo simulations within the local density approximation. Monte Carlo simulated colloid scale density profiles display some minor differences with confocal images in terms of microstructural transitions involving the onset of interfacial crystallization and the precise elevation of the fluid-solid interface. These discrepancies are attributed to polydispersity unaccounted for in the analyses, sensitivity of the perturbation theory to the effective hard sphere size, and the influence of ensemble, system size, and box shape in Monte Carlo simulations involving anisotropic/inhomogeneous solids. Successful demonstration of intensity based confocal microscopy provides a basis for future measurements of three-dimensional colloidal interactions, dynamics, and structure near surfaces.     

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.     

Structural properties of a model system with effective interparticle interaction potential applicable in modeling of complex fluids

We report grand canonical ensemble Monte Carlo (MC) simulation and theoretical studies of the structural properties of a model system described by an effective interparticle interaction potential, which incorporates basic interaction terms used in modeling of various complex fluids composed of mesoscopic particles dispersed in a solvent bath. The MC results for the bulk radial distribution function are employed to test the validity of the hard-sphere bridge function in combination with a modified hypernetted chain approximation (MHNC) in closing the Ornstein-Zernike (OZ) integral equation, while the MC data for the density profiles in different inhomogeneous environments are used to assess the validity of the third-order+second-order perturbation density functional theory (DFT). We found satisfactory agreement between the results predicted by the pure theories and simulation data, which classifies the proposed theoretical approaches as convenient tools for the investigation of complex fluids. The present investigation indicates that the bridge function approximation and density functional approximation, which are traditionally used for the study of neutral atomic fluids, also perform well for complex fluids only on condition that the underlying effective potentials include a highly repulsive core as an ingredient.     

Phase behavior of ionic fluids in slitlike pores: a density functional approach for the restricted primitive model

We present a density functional theory of nonuniform ionic fluids. This theory is based on the application of the electrostatic contribution to the free energy functional arising from mean spherical approximation for a bulk restricted primitive model and from the energy route bulk equation of state. In order to employ this functional we define a reference fluid and additional averaged densities, according to the approach introduced by Gillespie, Nonner and Eisenberg [J. Phys.: Condens. Matter 14, 12129 (2002)]. In the case of bulk systems the proposed theory reduces to the mean spherical approximation equation of state, arising from the energy route and thus it predicts the first-order phase transition. We use this theory to investigate the effects of confinement on the liquid-vapor equilibria. Two cases are considered, namely an electrolyte confined to the pore with uncharged walls and with charged walls. The dependence of the capillary evaporation diagrams on the pore width and on the electrostatic potential is determined.