基于密度泛函理论研究二元排斥Yukawa流体的表面结构性质

基于自由能密度泛函理论(DFT)考察了二元排斥Yukawa (HCRY)流体在不同外场下的密度分布. 基于微扰理论, 体系的Helmholtz自由能泛函采用硬球排斥部分和长程色散部分贡献之和, 其中Kierlik和Rosinberg的加权密度近似(WDA)被用来计算硬球排斥部分, 而色散部分采用平均场理论(MFT)进行描述. 为了验证DFT计算结果的合理性, 研究中采用巨正则Monte Carlo(GCMC)模拟计算了在不同主体相密度、硬核直径和位能参数比的条件下二元HCRY混合流体的密度分布. 结果表明, 该DFT计算结果与GCMC模拟值吻合良好.     

牛血清蛋白在NaCl水溶液中渗透压的研究

将NaCl水溶液中带电蛋白质分子间的静电排斥作用用Yukawa位能函数描述,蛋白质分子间的色散作用用Lennard-Jones位能函数描述,代替了经典DLVO理论。基于McMillan-Mayer渗透压统计理论,将Duh和Mier-Y-Teran的Yukawa状态方程与Cotterman等的Lennard-Jones微扰公式相结合,建立了一个新的状态方程,并采用该状态方程研究了含有NaCl的牛血清蛋白(简称BSA)水溶液的渗透压,依据BSA在0.15mol/LNaCl水溶液中的回归参数可预测其在1～5mol/LNaCl水溶液中的渗透压。该状态方程采用了两个回归参数(硬球直径和LJ参数),具有一定的关联和预测效果,并与己有的其它理论模型进行了比较。     

硬核Yukawa流体中胶体粒子间的耗尽势

基于Roth、Evans和Dietrich有关耗尽势的密度泛函理论研究了硬核Yukawa(HCY)流体中胶体粒子间的耗尽势.在极稀溶液条件下,通过计算两个胶体粒子在不同条件下的耗尽势,分析了HCY流体的相关因素对耗尽势的影响.结果表明,胶体粒子与溶剂分子的尺寸比率、HCY流体分子间的相互作用、HCY流体分子的体相密度以及胶体粒子与流体分子之间的相互作用等因素均可对胶体粒子间耗尽势产生显著影响.研究结果可为实验上调控胶体粒子间的相互作用提供可能的线索.     

TBP体系活度系数的测定和微扰理论的研究

在前人工作的基础上, 改进了液上空间气相色谱测活度系数的方法, 实验测定了TBP. 稀释剂和水形成的多个二元系、三元系和四元系的活度系数和密度. 选用的稀释剂有n-C_6H_(14)、n-C_7H_(16)、n-C_8H_(18)、C_6H_6、cy-C_6H_(12)、CCl_4和CHCl_3. 在Pierotti理论的基础上, 采用新的硬球作用表达式和径向分布函数, 并计及分子间的排斥能、色散能、取向能和诱导能, 建立了简单的活度系数模型, 并用于TBP和稀释剂体系的计算. 从二元系回归得到的分子参数较好地预测了三元系的活度系数.     

渗透压法对牛血清蛋白在NaCl水溶液中相互作用的研究

测量了牛血清蛋白不同pH值和不同离子强度Nacl水溶液中的渗透压,计算了渗透第二维里系数,按照Prausnitz提出的分子热学模型计算了蛋白质之间的静电排斥能、色散吸引能和离子排斥体积产生的吸引势能,并对这三种势能与溶液pH值和离子强度的变化关系进行讨论。     

微扰理论非原始模型用于1:1价单一电解质水溶液

本文采用微扰理论非原始模型, 以带电硬球混合物为参考体系, 考虑粒子间各种作用能(色散、静电、偶极、四极、诱导), 首次取相对介电常数为1, 拟合了12处1:1价电解质水溶液的渗透系数, 获得了成功, 得到了7种1价离子Na^+,K^+, Rb^+, Cs^+, Cl^-, Br^-, I^-的微观参数(软球直径σ和色散能常数ε/k)。计算得到的电解质水溶液渗透系数的总平均绝对偏差是0.041。这些离子的微观参数在不同体系中维持不变。计算中未引入混合参数。     

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

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

Mg2NiQ6(Q=H,D, T)体系的热力学氢同位素效应

为研究Mg2NiQ4(Q=H, D, T)体系的热力学氢同位素效应, 基于量子力学第一原理, 采用密度泛函与赝势平面波相结合的方法, 计算了Mg2NiQ6(Q=H, D, T)体系的声子色散谱, 利用声子色散谱得到了热力学函数. 利用文献报道的氢同位素气体分子的热力学数据, 采用热力学方法分析了Mg2Ni吸氢形成Mg2NiQ6的同位素效应. 研究结果表明, Mg2Ni吸氢的同位素效应主要是原子相对振动的频率不同导致的. CaF2结构的Mg2NiQ6的氢同位素效应随温度升高由负同位素效应转变为正同位素效应.     

Mg2NiQ6（Q=H，D，T）体系的热力学氢同位素效应

为研究Mg2NiQ4(Q=H,D,T)体系的热力学氢同位素效应,基于量子力学第一原理,采用密度泛函与赝势平面波相结合的方法,计算了Mg2NiQ6(Q=H,D,T)体系的声子色散谱,利用声子色散谱得到了热力学函数.利用文献报道的氢同位素气体分子的热力学数据,采用热力学方法分析了Mg2Ni吸氢形成Mg2NiQ6的同位素效应.研究结果表明,Mg2Ni吸氢的同位素效应主要是原子相对振动的频率不同导致的.CaF2结构的Mg2NiQ6的氢同位素效应随温度升高由负同位素效应转变为正同位素效应.     

电子亲和势法判断负离子聚合单体的活性

提出了利用不同取代基烯类单体的电子亲和势来判断其负离子聚合反应活性的方法.采用密度泛函理论的B3LYP/6-31G(d)方法优化了不同取代基烯类单体几何构型,在B3LYP/6-311++G(3df,2p)水平上计算了其电子亲和势.通过电子亲和势计算值与文献报道实验数据比较,表明本文采用的计算方法是比较可靠的.结合不同取代基烯类单体的电子亲和势的计算结果,通过与Q-e关系及取代基常数σ数据进行比较表明,电子亲和势可以用来判断不同单体负离子聚合反应的活性高低.     

Systematic coarse-graining of potential energy landscapes and dynamics in liquids

Recent efforts have shown that the dynamic properties of a wide class of liquids can be mapped onto semi-universal scaling laws and constitutive relations that are motivated by thermodynamic analyses of much simpler models. In particular, it has been found that many systems exhibit dynamics whose behavior in state space closely follows that of soft-sphere particles interacting through an inverse power repulsion. In the present work, we show that a recently developed coarse-graining theory provides a natural way to understand how arbitrary liquids can be mapped onto effective soft-sphere models and hence how one might potentially be able to extract underlying dynamical scaling laws. The theory is based on the relative entropy, an information metric that quantifies how well a soft-sphere approximation to a liquid's multidimensional potential energy landscape performs. We show that optimization of the relative entropy not only enables one to extract effective soft-sphere potentials that suggest an inherent scaling of thermodynamic and dynamic properties in temperature-density space, but that also has rather interesting connections to excess entropy based theories of liquid dynamics. We apply the approach to a binary mixture of Lennard-Jones particles, and show that it gives effective soft-sphere scaling laws that well-describe the behavior of the diffusion constants. Our results suggest that the relative entropy formalism may be useful for "perturbative" type theories of dynamics, offering a general strategy for systematically connecting complex energy landscapes to simpler reference ones with better understood dynamic behavior.     

Thermodynamic properties of non-conformal soft-sphere fluids with effective hard-sphere diameters

In this work we study a set of soft-sphere systems characterised by a well-defined variation of their softness. These systems represent an extension of the repulsive Lennard-Jones potential widely used in statistical mechanics of fluids. This type of soft spheres is of interest because they represent quite accurately the effective intermolecular repulsion in fluid substances and also because they exhibit interesting properties. The thermodynamics of the soft-sphere fluids is obtained via an effective hard-sphere diameter approach that leads to a compact and accurate equation of state. The virial coefficients of soft spheres are shown to follow quite simple relationships that are incorporated into the equation of state. The approach followed exhibits the rescaling of the density that produces a unique equation for all systems and temperatures. The scaling is carried through to the level of the structure of the fluids.     

Correlation and Prediction of Osmotic Pressures for Aqueous Bovine Serum Albumin-NaCl Solutions Based on Two Yukawa Potentials

A new equation of state is proposed to correlate and predict the osmotic pressure data for aqueous bovine serum albumin (BSA) solutions with different NaCl concentrations and pH values with only one adjustable parameter. The Carnahan-Starling equation represents the contribution of the hard sphere repulsion to the osmotic pressure. The attractive dispersion and double-layer repulsion interactions are represented by two Yukawa potentials, respectively. The equation of state of Duh and Mier-Y-Teran for one Yukawa potential is expanded to two Yukawa potentials to describe the nonidealities of the charged BSA-aqueous NaCl solution, instead of the classical DLVO theory. The average relative deviation of correlation of the osmotic pressure in 0.15 M NaCl solution is 18%. The average relative deviation of prediction in 1-5 M NaCl solutions is 20.33%. A comparison with other models and the limitations of our model are discussed. Copyright 2001 Academic Press.     

Equation of state for expanded fluid mercury: variational theory with many-body interaction

A variational associating fluid theory is proposed to describe equations of state for expanded fluid mercury. The theory is based on the soft-sphere variational theory, incorporating an ab initio diatomic potential and an attractive many-body potential; the latter is evaluated with quantum chemical methods and expressed as a function of the local atomic coordination number and the nearest-neighbor distance. The resultant equation of state can reproduce the observed gas-liquid coexistence curve with good accuracy, without introducing phenomenological effective pair potentials. Various thermodynamic quantities such as pressure, isocloric thermal pressure coefficient, adiabatic sound velocity, and specific heat are calculated over a wide density-temperature range and compared with available experimental data.     

Electron interaction in molecules

An idea of electron interaction in molecule has been applied to the SCF MO calculations of theπ-electronic structure of some complex aromatic hydrocarbons and their derivatives. The theoretical results for singlet and triplet transition energies, first ionization potentials and bond lengths agree fairly well with the experimental data. A correlation equation between the valence state ionization potential and the one center electron repulsion integral has been proposed. It has been shown that the electron repulsion in molecule is considerably smaller than in free atom. The present calculation shows that we can treat sulphur as a normal heteroatom analogous to oxygen and nitrogen.     

H原子在Ru(1121)台阶面上的吸附位和吸附态

利用原子和表面簇合物相互作用的五参数Morse势方法(简称5-MP)研究了H-Ru(1121)台阶面吸附体系,并利用推广的LEPS方法研究了2H-Ru(1121)体系,探讨了表面吸附态之间的相互作用.研究结果表明,H原子在开放的Ru(1121)台阶面表面上的晶胞存在4种不等价的表面三重吸附态.还存在两类与表面有直接扩散通道的内层吸附态.理论分析表明,TDS实验出现的α态是与内层吸附态有强相互作用的强排斥表面吸附态,β态为弱排斥表面吸附态,而γ态为与内层吸附态无相互作用的表面吸附态.实验上将α态归属为四重吸附态的推测没有获得理论结果的支持.     

Study on osmotic pressures for aqueous lysozyme and alpha-chymotrypsin-electrolyte solutions with two Yukawa potentials

The equation of state (EOS) of Duh and Mier-y-Terán for one Yukawa potential is expanded to two Yukawa potentials to describe the nonidealities of the charged lysozyme and alpha-chymotrypsin solutions. Instead of the classical DLVO theory, the attractive dispersion and double-layer repulsion interactions are represented by two Yukawa potentials, respectively. For the aqueous lysozyme solutions, the only adjustable dispersion energy parameter epsilon/k is regressed and the average deviation is 1.76%. For the aqueous alpha-chymotrypsin solutions, two adjustable parameters (the molecular weight and dispersion energy parameter) are regressed and the average deviation is 7.62%. Some correlation and prediction results are discussed.     

Calculation of high pressure vapour—liquid equilibria for systems containing polar substances with the use of an augmented van der Waals equation of state

An augmented van der Waals equation of state based on a perturbation theory has been applied to the calculation of high pressure vapour—liquid equilibria for systems containing polar substances. The equation of state comprises four terms, which imply the contributions from repulsion, symmetric, non-polar asymmetric, and polar asymmetric interactions. The characteristic parameters of each pure substance have been determined by three methods with the use of vapour pressures and saturated liquid densities. Mixing models for the terms of the repulsion, symmetric, and non-polar asymmetric interactions are the same as used previously. Two types of mixing models based on a three-fluid model and/or a one-fluid model are developed for the polar asymmetric term. The polar asymmetric term has a large effect on the prediction of the vapour—liquid equilibrium. With the introduction of a binary interaction parameter, the equation is found to be useful in correlating the vapour—liquid equilibria for a system containing a polar substance except near a critical region.     

一维函数方法求解原子和分子薛定谔方程

对薛定谔方程的严格数值求解, 尤其是发展标准方法之外的、 包含新功能的解法, 一直是物理学研究的基本关注点. 本文介绍一种近些年发展的一维函数近似解方法, 该方法通过对波函数的不同坐标分量进行处理来求解原子和分子体系的薛定谔方程. 电子的试探波函数被离散化到实空间均匀格点上, 因此可以通过残差矢量校正的方法对其进行改进. 一维函数方法本身的特征决定其非常利于数值积分, 避免了很多由常规的多电子、 多中心势分子积分所带来的问题. 计算中, 最终能量是从严格的能量上限逐渐收敛所获得, 计算出的两电子薛定谔波函数呈现出常规单电子近似方法所含有的电子关联效应. 不同于密度泛函理论及Hartree-Fock的单电子解法, 本方法对电子-电子排斥能的多体效应的处理更加精确.