Yukawa流体的相平衡和界面张力研究

应用二阶微扰理论, Duh-Mier-Y-Teran状态方程和在平均球近似(mean spherical approximation, MSA)的基础上获得的直接相关函数, 建立了适用于均匀流体和非均匀流体的状态方程. 结合此状态方程, 重整化群理论(renormalization group theory, RG)和密度泛函理论(density functional theory, DFT), 分别研究了Yukawa流体的相平衡和界面张力. 结果与分子模拟数据吻合良好.     

用Yukawa状态方程和重整化群理论计算链状流体的汽液相平衡

结合重整化(Renormalization Group，RG)理论和适用于链状流体的Yukawa状态方程，研究了10个非极性链状流体和6个缔合流体的汽液相平衡和临界性质．16个体系的临界温度Tc，临界压力Pc和临界密度ρc的计算值和实验值之间的平均相对偏差(ARD％)分别为0．44％，1．31％和1．01％．所得结果与实验数据基本吻合．     

基于格子的链状流体分子热力学模型

格子模型是高分子溶液理论最常采用的流体模型,在此基础上建立的分子热力学模型在流体混合物的热力学性质、相平衡行为的计算等方面有广泛的应用．最近,我们对基于格子的混合亥氏函数模型重新进行了审视,采用统计力学理论推导与计算机模拟相结合的现代分子热力学研究方法建立了新的分子热力学模型,可以反映高分子链枝化、配位数、链刚性、共聚物链组成、氢键、压力等的影响,对小分子系统和高分子溶液的热力学性质和相平衡关系的计算不仅与Monte Carlo模拟结果吻合,比现有其他理论预测效果更好,对高分子溶液、离子液体混合物等实际系统的相平衡计算也取得令人满意的结果,显示出模型优越的工程应用价值．本文对上述工作进行了系统总结．     

多元气-液-固复杂体系相平衡

石油流体中含有气相、液相及可能遇到的固相包括水合物、石蜡和沥青质等,涉及多元气-液-固复杂体系的相平衡问题.为防止这些沉积物堵塞造成安全隐患,需要确定水合物、石蜡、沥青质沉积起始条件以及沉积量.本文针对化学热力学理论在含水合物、石蜡和沥青质的多元-多相平衡研究中的应用进行了综述.水合物相平衡模型较为成熟,主要有两类,其一为基于等温吸附理论的van der Waals-Platteeuw型热力学模型;其二为基于双过程水合物生成机理的Chen-Guo水合物热力学模型.石蜡沉积一般采用活度系数法、状态方程法及多固相模型描述.沥青质絮凝、沉积则可采用溶解度参数模型、状态方程法、胶体模型和标度理论模型进行计算.同时对多元气-液-固复杂体系的相平衡研究发展方向进行了展望.     

跨接VTSRK方程计算烃类的pVT性质

烃类pVT性质的精细表征对能源动力、化工等领域应用有重要价值,临界区热力性质描述是难点之一.本文建立了烷烃(C1-C20)的跨接比容平移Soave-Redlich-Kwong(SRK)(跨接VTSRK)状态方程,在SRK状态方程的基础上引入了比容平移和跨接方法,以改善饱和液相密度和近临界区域热力学性质的计算精度,方程参数被表达为物质临界参数和偏心因子的函数.比较结果表明,跨接方程对烷烃(C1-C20)饱和蒸气压、饱和气相密度、饱和液相密度的计算平均偏差分别为1.01%、1.83%和0.93%,显著优于原方程,单相区和近临界区的pVT性质计算精度也比原状态方程有较大改善.进一步将方程推广到环烷烃(环丙烷、环戊烷和环己烷)和苯、甲苯的计算,也获得了较好效果,验证了方程的预测能力.     

跨接VTSRK方程计算烃类的pVT性质

烃类pVT性质的精细表征对能源动力、化工等领域应用有重要价值,临界区热力性质描述是难点之一.本文建立了烷烃（C1-C20）的跨接比容平移Soave-Redlich-Kwong（SRK）（跨接VTSRK）状态方程,在SRK状态方程的基础上引入了比容平移和跨接方法,以改善饱和液相密度和近临界区域热力学性质的计算精度,方程参数被表达为物质临界参数和偏心因子的函数. 比较结果表明,跨接方程对烷烃（C1-C20）饱和蒸气压、饱和气相密度、饱和液相密度的计算平均偏差分别为1.01%、1.83%和0.93%,显著优于原方程,单相区和近临界区的pVT性质计算精度也比原状态方程有较大改善. 进一步将方程推广到环烷烃（环丙烷、环戊烷和环己烷）和苯、甲苯的计算,也获得了较好效果,验证了方程的预测能力.     

A_aD_d型氢键流体的统计理论(V):球形微腔中的相平衡

利用密度泛函理论并结合改进的基本度量理论研究了球形微腔中AaDd型氢键流体的相态结构.首先,根据氢键流体在球腔中的吸附-脱附等温线以及相应的巨势等温线获得不同条件下氢键流体的相图.在此基础上,重点讨论了氢键作用、球腔尺寸以及腔壁与流体之间的相互作用等因素对氢键流体相平衡特征的影响.结果表明,流体层化转变和毛细凝聚的临界温度、临界密度和临界相区域等相态特征与这些因素密切相关.研究结果可为进一步揭示几何约束下氢键流体的相平衡及聚集态结构提供可能的理论线索.     

醛酮类化合物的分子力场参数推导及热力学性质计算

针对醛酮类分子, 基于第一原理的量子化学计算结果开发了准确、可迁移的全原子力场. 利用所得到的力场, 采用分子动力学和蒙特卡罗方法对醛酮类分子的气相分子结构、振动频率、构象能和凝聚相液体性质(密度和蒸发焓), 以及流体的多种热力学性质包括气液相平衡和临界性质, 与涨落相关的等压热容, 传递性质如剪切粘度等进行了计算和预测. 计算结果表明, 该力场可准确地反映气相分子的结构、振动频率、构象能和凝聚相液体密度和蒸发焓等性质, 并准确地预测一系列醛酮类分子的多种热力学性质. 所开发的醛酮类分子力场函数形式简单, 并具有良好的可迁移性和准确性, 在应用范围和计算精度上都有明显的改善和提高.     

用重整化群理论研究胶体模型体系的相行为

用yukawa势能函数表达胶体颗粒之间的吸引作用。用Duh-Mier-Y-Teran状态方 程表达液相Helmholtz自由能。用一阶微扰理论、固体硬球径向分布函数解析式和 改进的胞腔模型建立固相状态方程，结合建立的状态方程和重整化群理论。研究了 胶体模型体系的液-液相平衡和液-固相平衡。研究表明，颗粒之间色散作用量程参 数的变化对胶休到本世纪末茶杯 系的相行为有特殊需要影响。所得结果与分子模 拟数据吻合良好。     

AaDd型氢键流体的统计理论（V）：球形微腔中的相平衡

利用密度泛函理论并结合改进的基本度量理论研究了球形微腔中AaDd型氢键流体的相态结构．首先,根据氢键流体在球腔中的吸附一脱附等温线以及相应的巨势等温线获得不同条件下氢键流体的相图．在此基础上,重点讨论了氢键作用、球腔尺寸以及腔壁与流体之间的相互作用等因素对氢键流体相平衡特征的影响．结果表明,流体层化转变和毛细凝聚的临界温度、临界密度和临界相区域等相态特征与这些因素密切相关．研究结果可为进一步揭示几何约束下氢键流体的相平衡及聚集态结构提供可能的理论线索．     

Thermodynamics for fluid mixtures near to and far from the vapor–liquid critical point

To describe thermodynamic properties of fluid mixtures near to and far from the vapor–liquid critical point, we need a method where a classical equation-of-state is augmented with a correction based on renormalization group (RG) theory. The advantage of the method described here is that, subject to well-defined assumptions, it can be applied not only to binary mixtures but to mixtures containing any number of components. While our method is based on White’s recursion procedure, our extension to mixtures is based on the isomorphism assumption and on an approximation suggested by Kiselev. To illustrate, calculations are presented for the critical loci of some alkane mixtures containing one discrete component and one pseudo-component whose composition is characterized by a continuous distribution of molecular weight. While critical loci calculated with the RG correction are similar to those calculated by the classical equation-of-state alone, inclusion of the RG correction provides better agreement with experiment.     

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 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.     

An improved renormalization group theory for real fluids

On the basis of White's theory, an improved renormalization group (RG) theory is developed for chain bonding fluids inside the critical region. Outside the critical region, the statistical associating fluid theory based on the first-order mean sphere approximation [Fluid Phase Equilibria 171, 27 (2000)] is adopted and all the microscopic parameters are taken directly from its earlier application of real fluids. Inside the critical region, the RG transformation for long-range density fluctuation is derived in the k space, which illustrates explicitly the contributions from the mean-field term, the local density fluctuation, and the nonlocal density fluctuation. The RG theory is applied to describe physical behavior of ten n alkanes (C1-C10) both near to and far from the critical point. With no additional parameters for chain bonding fluids, good results are obtained for critical specific heat and phase coexistence curves and the resulting critical exponents are in good agreement with the reported nonclassic values.     

Wetting study of imidazolium ionic liquids

In this work, we present a systematic contact angles study of a series of 1-alkyl, 3-methyl-imidazolium ionic liquids (ILs) on well-defined polar and nonpolar monolayer surfaces supported on Si wafers. The advancing and receding contact angles of ILs were used to determine the surface energy of the monolayer surfaces using Neumann's equation-of-state and Zisman's critical surface tension approaches. In parallel, the contact angles of conventional probe fluids (molecular liquids) including water, formamide, methylene iodide, ethylene glycol, and hexadecane were determined on the same surfaces. The results obtained showed a great deal of similarity in wetting behavior of ionic vs molecular probe fluids: the contact angles of both sets of liquids followed the same patterns in accord with the surface tension of the fluid. A good agreement was found between the surface energy determined by different sets of liquids.     

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.     

Global fluid phase equilibria and critical phenomena of selected mixtures using the crossover soft-SAFT equation

We present here the extension of the crossover soft-statistical associating fluid theory (soft-SAFT) equation of state to mixtures, as well as some illustrative applications of the methodology to mixtures of particular scientific and technological interest. The procedure is based on White's work (White, J. A. Fluid Phase Equilib. 1992, 75, 53) from the renormalization group theory, as for the pure fluids, with the isomorphism assumption applied to the mixtures. The equation is applied to three groups of mixtures: selected mixtures of n-alkanes, the CO2/n-alkane homologous series, and the CO2/1-alkanol homologous series. The crossover equation is first applied to the pure components of the mixtures, CO2 and the 1-alkanol family, while an available correlation is used for the molecular parameters of the n-alkane series (Llovell et al. J. Chem. Phys 2004, 121, 10715). A set of transferable molecular parameters is provided for the 1-alkanols series; these are accurate for the whole range of thermodynamic conditions. The crossover soft-SAFT equation is able to accurately describe these compounds near to and far from the critical point. The theory is then used to represent the phase behavior and the critical phenomena of the selected mixtures. We use binary interaction parameters xi and eta for dissimilar mixtures. These parameters are fitted at some particular conditions (one subcritical temperature or binary critical data) and used to predict the behavior of the mixture at different conditions (other subcritical conditions and/or critical conditions). The equation is able to capture the continuous change in the critical behavior of the CO2/n-alkane and the CO2/1-alkanol homologous series as the chain length of the second compound increases. Excellent agreement with experimental data is obtained, even in the most nonideal cases. The new equation is proved to be a powerful tool to study the global phase behavior of complex systems, as well as other thermodynamic properties of very challenging mixtures.     

Cubic equation-of-state correlation of the solubility of some anti-inflammatory drugs in supercritical carbon dioxide

The accurate experimental determination of solid drug solubility in supercritical fluids (SCFs) and its correlation are crucially important to the development of supercritical technologies for the pharmaceutical industry. In this work, the solubilities of flurbiprofen, ketoprofen, naproxen and ibuprofen in supercritical carbon dioxide (scCO₂) were correlated using the Peng–Robinson (PR), Soave–Redlich–Kwong (SRK) and Patel–Teja–Valderrama (PTV) equations-of-state with van der Waals (vdW), Panagiotopoulos–Reid (mrPR) and Mukhopadhyay–Rao (MR) mixing rules. Several modeling and correlation computational programs were developed in Mathematica^®, a powerful symbolic computational language. Correlations were compared and discussed on the basis of the employed equation-of-state and mixing/combining rules. The importance of the estimation methods used for the determination of critical and other drug physical properties (namely drug's sublimation pressure) was also underlined in the discussion of the correlation results quality. A potential helpful procedure (based on simultaneous correlation of one of the critical or physicochemical properties) is suggested, in order to satisfactorily choose property estimation methods to use.     

Phase Behavior of Polymer Blends

Major emphasis is placed on the phase behavior of miscible polymer blends. To understand the complex phase behavior of blends, a refined version of an equation-of-state theory is discussed. This theory makes the simultaneous occurrence of upper critical solution temperature and lower critical solution temperature in blends of high molar mass polymers conceivable. The kinetics of isothermal phase dissolution as it emanates from different experimental routes is discussed in terms of CahN′s linearized theory of phase separation. The rate of phase dissolution varies as a function of quench depth, which indicates the rate is directed by both the thermodynamic driving force and the mobility.