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.     

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.     

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.     

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.     

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

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

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.     

Three-dimensional potential energy surface of the Ar-OH(2Pi i) complex

Pure rotational transitions in the ground state for Ar-OH and Ar-OD [Y. Ohshima et al., J. Chem. Phys. 95, 7001 (1991) and Y. Endo et al., Faraday Discuss. 97, 341 (1994)], those in the excited states of the OH vibration, nu(s)=1 and 2, observed by Fourier-transform microwave spectroscopy in the present study, rotation-vibration transitions observed by infrared-ultraviolet double-resonance spectroscopy [K. M. Beck et al., Chem. Phys. Lett. 162, 203 (1989) and R. T. Bonn et al., J. Chem. Phys. 112, 4942 (2000)], and the P-level structure observed by stimulated emission pumping spectroscopy [M. T. Berry et al., Chem. Phys. Lett. 178, 301 (1991)] have been simultaneously analyzed to determine the potential energy surface of Ar-OH in the ground state. A Schrodinger equation, considering all the freedom of motions for an atom-diatom system in the Jacobi coordinate, R, theta, and r, was numerically solved to obtain energies of the rovibrational energy levels using the discrete variable representation method. A three-dimensional potential energy surface is determined by a least-squares fitting. In the analysis the potential parameters, obtained by ab initio calculations at the RCCSD(T) level of theory with a set of basis functions of aug-cc-pVTZ and midbond functions, are used as initial values. The determined intermolecular potential energy surface and its dependence on the OH monomer bond length are compared with those of an isovalent radical complex, Ar-SH.     

Phase coexistence in polydisperse multi-Yukawa hard-sphere fluid: high temperature approximation

High temperature approximation (HTA) is used to describe the phase behavior of polydisperse multi-Yukawa hard-sphere fluid mixtures. It is demonstrated that in the frames of the HTA the model belongs to the class of "truncatable free energy models," i.e., the models with thermodynamical properties (Helmholtz free energy, chemical potential, and pressure) defined by the finite number of generalized moments. Using this property we were able to calculate the complete phase diagram (i.e., cloud and shadow curves as well as binodals) and size distribution functions of the coexisting phases of several different models of polydisperse fluids. In particular, we consider polydisperse one-Yukawa hard-sphere mixture with factorizable Yukawa coefficients and polydisperse Lennard-Jones (LJ) mixture with interaction energy parameter and/or size polydispersity. To validate the accuracy of the HTA we compare theoretical results with previously published results of more advanced mean spherical approximation (MSA) for the one-Yukawa model and with the Monte Carlo (MC) computer simulation results of [Wilding et al. J. Chem. Phys. 121, 6887 (2004); Phys. Rev. Lett. 95, 155701 (2005)] for the LJ model. We find that overall predictions of the HTA are in reasonable agreement with predictions of the MSA and MC, with the accuracy range from semiquantitative (for the phase diagram) to quantitative (for the size distribution functions).     

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.     

Theoretical scheme for the shear viscosity of Lennard-Jones fluids

We use the shear viscosity expression from the Enskog theory of dense gases in a perturbative scheme for the Lennard-Jones (LJ) fluid. This perturbative scheme is formulated by combining the analytic rational function approximation method of Bravo Yuste and Santos [Phys. Rev. A 43, 5418 (1991)] for the radial distribution function of hard-sphere fluids and the well known Mansoori-Canfield/Rasaiah-Stell perturbation theory to determine an effective diameter for the LJ fluid. The scheme is reliable on a wide range of temperatures and densities, and is very accurate around the critical point. Using this information, we build an accurate empirical formula for the shear viscosity in the liquid phase, which fits the recent data [K. Meier et al., J. Chem. Phys. 121, 3671 (2004)] in the whole simulation range.     

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.     

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.     

Optimized theory for simple and molecular fluids

An optimized closure approximation for both simple and molecular fluids is presented. A smooth interpolation between Perkus-Yevick and hypernetted chain closures is optimized by minimizing the free energy self-consistently with respect to the interpolation parameter(s). The molecular version is derived from a refinement of the method for simple fluids. In doing so, a method is proposed which appropriately couples an optimized closure with the variant of the diagrammatically proper integral equation recently introduced by this laboratory [K. M. Dyer et al., J. Chem. Phys. 123, 204512 (2005)]. The simplicity of the expressions involved in this proposed theory has allowed the authors to obtain an analytic expression for the approximate excess chemical potential. This is shown to be an efficient tool to estimate, from first principles, the numerical value of the interpolation parameters defining the aforementioned closure. As a preliminary test, representative models for simple fluids and homonuclear diatomic Lennard-Jones fluids were analyzed, obtaining site-site correlation functions in excellent agreement with simulation data.     

Rovibrational states of the H2O-H2 complex: an ab initio calculation

All bound rovibrational levels of the H(2)O-H(2) dimer are calculated for total angular momentum J = 0-5 on two recent intermolecular potential surfaces reported by Valiron et al. [J. Chem. Phys. 129, 134306 (2008)] and Hodges et al. [J. Chem. Phys. 120, 710 (2004)] obtained through ab initio calculations. The method used handles correctly the large amplitude internal motions in this complex; it involves a discrete variable representation of the intermolecular distance coordinate R and a basis of coupled free rotor wave functions for the hindered internal rotations and the overall rotation of the dimer. The basis is adapted to the permutation symmetry associated with the para/ortho (p/o) nature of both H(2)O and H(2) as well as to inversion symmetry. Dimers containing oH(2) are more strongly bound than dimers with pH(2), as expected, with dissociation energies D(0) of 33.57, 36.63, 53.60, and 59.04 cm(-1)for pH(2)O-pH(2), oH(2)O-pH(2), pH(2)O-oH(2), and oH(2)O-oH(2), respectively, on the potential of Valiron et al. that corresponds to a binding energy D(e) of 235.14 cm(-1). Rovibrational wave functions are computed as well and the nature of the bound states in the four different dimer species is discussed. Converged rovibrational levels on both potentials agree well with the high-resolution spectrum reported by Weida and Nesbitt [J. Chem. Phys. 110, 156 (1999)]; the hindered internal rotor model that was used to interpret this spectrum is qualitatively correct.     

Molecular nitrogen-N2 properties: the intermolecular potential and the equation of state

Quantum mechanical (QM) high precision calculations were used to determine N(2)-N(2) intermolecular interaction potential. Using QM numerical data the anisotropic potential energy surface was obtained for all orientations of the pair of the nitrogen molecules in the rotation invariant form. The new N(2)-N(2) potential is in reasonably good agreement with the scaled potential obtained by van der Avoird et al. using the results of Hartree-Fock calculations [J. Chem. Phys. 84, 1629 (1986)]. The molecular dynamics (MD) of the N(2) molecules has been used to determine nitrogen equation of state. The classical motion of N(2) molecules was integrated in rigid rotor approximation, i.e., it accounted only translational and rotational degrees of freedom. Fincham [Mol. Simul. 11, 79 (1993)] algorithm was shown to be superior in terms of precision and energy stability to other algorithms, including Singer [Mol. Phys. 33, 1757 (1977)], fifth order predictor-corrector, or Runge-Kutta, and was therefore used in the MD modeling of the nitrogen pressure [S. Krukowski and P. Strak, J. Chem. Phys. 124, 134501 (2006)]. Nitrogen equation of state at pressures up to 30 GPa (300 kbars) and temperatures from the room temperature to 2000 K was obtained using MD simulation results. Results of MD simulations are in very good agreement (the error below 1%) with the experimental data on nitrogen equation of state at pressures below 1 GPa (10 kbars) for temperatures below 1800 K [R. T. Jacobsen et al., J. Phys. Chem. Ref. Data 15, 735 (1986)]. For higher temperatures, the deviation is slightly larger, about 2.5% which still is a very good agreement. The slightly larger difference may be attributed to the vibrational motion not accounted explicitly by rigid rotor approximation, which may be especially important at high temperatures. These results allow to obtain reliable equation of state of nitrogen for pressures up to 30 GPa (300 kbars), i.e., close to molecular nitrogen stability limit, determined by Nellis et al. [Phys. Rev. Lett. 53, 1661 (1984)].     

Connection between the observable and centroid structural properties of a quantum fluid: application to liquid para-hydrogen

It is shown that the discrepancy between path integral Monte Carlo [M. Zoppi et al., Phys. Rev. B 65, 092204 (2002)] and path integral centroid molecular dynamics [F. J. Bermejo et al., Phys. Rev. Lett. 84, 5359 (2000)] calculations of the static structure factor of liquid para-hydrogen can be explained based on a deconvolution equation connecting centroid and physical radial distribution functions. An explicit expression for the kernel of the deconvolution equation has been obtained using functional derivative techniques. In the superposition approximation, this kernel is given by the functional derivative of the effective potential with respect to the pairwise classical potential. Results of path integral Monte Carlo calculations for the radial distribution function and the static structure factor of liquid para-hydrogen are presented.     

Closed-shell ring coupled cluster doubles theory with range separation applied on weak intermolecular interactions

We explore different variants of the random phase approximation to the correlation energy derived from closed-shell ring-diagram approximations to coupled cluster doubles theory. We implement these variants in range-separated density-functional theory, i.e., by combining the long-range random phase approximations with short-range density-functional approximations. We perform tests on the rare-gas dimers He(2), Ne(2), and Ar(2), and on the weakly interacting molecular complexes of the S22 set of Jurec?ka et al. [P. Jurec?ka, J. S?poner, J. C?erny?, and P. Hobza, Phys. Chem. Chem. Phys. 8, 1985 (2006)]. The two best variants correspond to the ones originally proposed by Szabo and Ostlund [A. Szabo and N. S. Ostlund, J. Chem. Phys. 67, 4351 (1977)]. With range separation, they reach mean absolute errors on the equilibrium interaction energies of the S22 set of about 0.4 kcal/mol, corresponding to mean absolute percentage errors of about 4%, with the aug-cc-pVDZ basis set.     

Local response dispersion method. II. Generalized multicenter interactions

Recently introduced local response dispersion method [T. Sato and H. Nakai, J. Chem. Phys. 131, 224104 (2009)], which is a first-principles alternative to empirical dispersion corrections in density functional theory, is implemented with generalized multicenter interactions involving both atomic and atomic pair polarizabilities. The generalization improves the asymptote of intermolecular interactions, reducing the mean absolute percentage error from about 30% to 6% in the molecular C(6) coefficients of more than 1000 dimers, compared to experimental values. The method is also applied to calculations of potential energy curves of molecules in the S22 database [P. Jure?ka et al., Phys. Chem. Chem. Phys. 8, 1985 (2006)]. The calculated potential energy curves are in a good agreement with reliable benchmarks recently published by Molnar et al. [J. Chem. Phys. 131, 065102 (2009)]. These improvements are achieved at the price of increasing complexity in the implementation, but without losing the computational efficiency of the previous two-center (atom-atom) formulation. A set of different truncations of two-center and three- or four-center interactions is shown to be optimal in the cost-performance balance.     

Modeling the phase behavior of H2S+n-alkane binary mixtures using the SAFT-VR+D approach

A statistical associating fluid theory for potential of variable range has been recently developed to model dipolar fluids (SAFT-VR+D) [Zhao and McCabe, J. Chem. Phys. 2006, 125, 104504]. The SAFT-VR+D equation explicitly accounts for dipolar interactions and their effect on the thermodynamics and structure of a fluid by using the generalized mean spherical approximation (GMSA) to describe a reference fluid of dipolar square-well segments. In this work, we apply the SAFT-VR+D approach to real mixtures of dipolar fluids. In particular, we examine the high-pressure phase diagram of hydrogen sulfide+n-alkane binary mixtures. Hydrogen sulfide is modeled as an associating spherical molecule with four off-center sites to mimic hydrogen bonding and an embedded dipole moment (micro) to describe the polarity of H2S. The n-alkane molecules are modeled as spherical segments tangentially bonded together to form chains of length m, as in the original SAFT-VR approach. By using simple Lorentz-Berthelot combining rules, the theoretical predictions from the SAFT-VR+D equation are found to be in excellent overall agreement with experimental data. In particular, the theory is able to accurately describe the different types of phase behavior observed for these mixtures as the molecular weight of the alkane is varied: type III phase behavior, according to the scheme of classification by Scott and Konynenburg, for the H2S+methane system, type IIA (with the presence of azeotropy) for the H2S+ethane and+propane mixtures; and type I phase behavior for mixtures of H2S and longer n-alkanes up to n-decane. The theory is also able to predict in a qualitative manner the solubility of hydrogen sulfide in heavy n-alkanes.