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

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

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

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

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

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.     

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

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

基于化学缔合统计理论的链状流体状态方程

基于化学缔合统计理论的链状流体状态方程(EOS)能够反映实际分子的形状、链节成链、缔合等具体信息,在实际流体热力学性质计算中有着广泛应用.一般的链状流体EOS仅考虑相邻链节间的相关性,我们则借助统计力学和计算机模拟结果在模型中纳入了相间链节间的相关性,获得的硬球链流体(HSCF)模型能够更好地预测模型流体的压缩因子和第二维里系数.以HSCF为参考,引入方阱色散微扰项获得了实际方阱链流体(SWCF)EOS;结合根据黏滞球模型导得的缔合项,进一步构建了缔合流体EOS.最近,我们根据微扰理论和积分方程方法又开发了一新的变阱宽方阱链流体(SWCF-VR)模型.SWCF和SWCF-VREOSs可很好地用于计算小分子、聚合物、离子液体等纯流体及混合物的相行为、热焓、表面张力、黏度等热力学及传递性质,显示了模型良好的工程应用价值.本文就本课题组多年来在自由空间范畴内基于化学缔合统计理论开发链状流体EOS及其实际应用作系统的总结.     

用密度梯度理论研究缔合流体的汽液成核性质

在密度梯度展开的基础上,将影响参数k 表达成温度的函数,建立了一个适用于均相和非均相缔合流体的状态方程。应用流体的蒸汽压和液相密度实验数据关联分子参数。在密度梯度理论的框架下,计算了水,重水,甲醇,乙醇,正丙醇,正丁醇,正戊醇和正己醇的成核速率并与实验数据进行了对比,计算结果令人满意。结果表明,密度梯度理论与密度泛函理论一样,可研究液核的结构和性质,但通过调整影响参数k, 可获得更为准确的成核速率。     

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

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

化学反应中的双稳性的临界现象及Landau相变理论

本文通过两个典型事例,即Langmuir形式的交换过程和Schlogle用以模拟一级相变的自催化反应,阐明化学反应系统的双稳性临界现象同样可以纳入Landau相变理论的模式。从而进一步揭示了双稳系统一级相变线临界点的二级(连续)相变特征,同时也揭示了作为远离平衡的化学反应双稳性临界现象的共性,以至于在更大范围内临界现象的共性。将相变的Landau理论,由原来的处于热平衡的凝聚物质推广到范围更宽广,表现更为丰富多彩的远离平衡系统的临界现象。     

柔性高分子/液晶体系的相平衡理论

在Flory-HugginS柔性高分子溶液理论和Lebwohl-Lasher液晶模型的基础上;导出了柔性高分子/液晶复合体系的混合自由能和化学位,并计算了相应的相图。计算结果表明,体系可出现各向同性相(Ⅰ)和各向同性相(Ⅰ),以及各向同性相(Ⅰ)和向列相(N)间的相平衡。并且发现,高分子链段与液晶分子间的各向同性相互作用对I/N的相平衡有重要影响。与实验结果的比较表明,本理论能完整、正确地描述高分子/液晶体系的相平衡。     

方阱二聚物相平衡性质的计算机模拟

在巨正则系综下, 采用Monte Carlo方法模拟了阱宽为1.5的方阱二聚物的临界性质, 利用完整标度结合histogram reweighting技术以及有限尺寸标度外推得到该流体的临界温度 , 临界密度 . 并采用Hyper-parallel tempering方法模拟了远离临界区域的气液相平衡性质, 大大提高了模拟效率. 所得临界性质较大幅度改进了原有模拟数据的准确度.     

含时密度矩阵重正化群的理论与应用

密度矩阵重正化群(DMRG)作为计算低维强关联体系强有力的方法为人熟知,在量子化学电子结构计算中得到广泛应用.最近几年,含时密度矩阵重正化群(TD-DMRG)的理论取得较快发展, TD-DMRG逐渐成为复杂体系量子动力学理论模拟的重要新兴方法之一.本文综述了基于矩阵乘积态(MPS)和矩阵乘积算符(MPO)的DMRG基本理论,并重点介绍了若干最常见的TD-DMRG时间演化算法,包括基于演化再压缩(P&C)的算法、基于含时变分原理(TDVP)的算法和时间步瞄准(TST)算法;还对利用TD-DMRG模拟有限温体系的纯化(Purification)算法和最小纠缠典型量子热态(METTS)算法进行了介绍.最后,对近年TD-DMRG在复杂体系量子动力学中的应用进行了总结.     

频域空间密度矩阵重正化群的研究进展

密度矩阵重正化群(DMRG)作为低维强关联体系中电子结构计算的强有力方法被广泛熟知,并被迅速地应用于量子化学,不仅在电子结构计算中发挥重要作用,同时也在近几年迅速地成为复杂体系量子动力学计算的重要方法.在DMRG框架中,衍生出了一系列计算动态响应性质的有效方法,并得到了广泛应用.本文简述了DMRG的基本理论,其矩阵乘积态(MPS)表示有效地扩展了该方法的应用范围.重点介绍了基于线性响应理论的动态DMRG,在频率空间求解系统在零温以及有限温度下响应性质的算法,并介绍其在电子关联问题和电子-声子关联问题中的应用,最后展望了该领域的未来发展方向.     

利用统计缔合流体理论状态方程预测混合气体水合物的平衡形成条件

利用统计缔合流体理论(SAFT)状态方程结合van der Waals-Platteeuw统计力学模型用于预测含有甲烷、乙烷、丙烷、乙烯、丙烯、H₂S, CO₂, N₂和H₂二元气体水合物的平衡形成条件. 对于气相和液相, 应用SAFT方程来描述. 在气相和液相相平衡的模拟过程中, SAFT方程考虑了硬球斥力、成链力、色散力以及缔合相互作用. 对于水合物相, 采用van der Waals-Platteuw模型来计算. 该方法的预测结果与实验数据吻合.     

Time-Dependent Density Matrix Renormalization Group Coupled with n-Mode Representation Potentials for the Excited State Radiationless Decay Rate: Formalism and Application to Azulene?

We propose a method for calculating the nonradiative decay rates for polyatomic molecules including anharmonic effects of the potential energy surface (PES) in the Franck-Condon region. The method combines the n-mode representation method to construct the ab initio PES and the nearly exact time-dependent density matrix renormalization group method (TD-DMRG) to simulate quantum dynamics. In addition, in the framework of TD-DMRG, we further develop an algorithm to calculate the final-state-resolved rate coefficient which is very useful to analyze the contribution from each vibrational mode to the transition process. We use this method to study the internal conversion (IC) process of azulene after taking into account the anharmonicity of the ground state PES. The results show that even for this semi-rigid molecule, the intramode anharmonicity enhances the IC rate significantly, and after considering the two-mode coupling effect, the rate increases even further. The reason is that the anharmonicity enables the C-H vibrations to receive electronic energy while C-H vibrations do not contribute on the harmonic PES as the Huang-Rhys factor is close to 0.     

A Density Matrix Renormalization Group Study of the Low-Lying Excited States of a Molybdenum Carbonyl-Nitrosyl Complex

A density matrix renormalization group-self consistent field (DMRG-SCF) study has been carried out to calculate the low-lying excited states of CpMo(CO)₂NO, a molybdenum complex containing NO and CO ligands. In order to automatically select an appropriate active space, a novel procedure employing the maximum single-orbital entropy for several states has been introduced and shown to be efficient and easy-to-implement when several electronic states are simultaneously considered. The analysis of the resulting natural transition orbitals and charge-transfer numbers shows that the lowest five excited electronic states are excitation into metal-NO antibonding orbitals, which offer the possibility for nitric oxide (NO) photorelease after excitation with visible light. Higher excited states are metal-centered excitations with contributions of metal-CO antibonding orbitals, which may serve as a gateway for carbon monoxide (CO) delivery. Time-dependent density functional theory calculations done for comparison, show that the state characters agree remarkably well with those from DMRG-SCF, while excitation energies are 0.4–1.0 eV red-shifted with respect to the DMRG-SCF ones.     

PSRK: A Group Contribution Equation of State Based on UNIFAC

A group contribution equation of state called PSRK (Predictive Soave-Redlich-Kwong) which is based on the Soave-Redlich-Kwong equation (Soave, 1972) has been developed. It uses the UNIFAC method to calculate the mixture parameter a and includes all already existing UNIFAC parameters. This concept makes use of recent developments by Michelsen (1990b) and has the main advantage, that vapor-uquid-equilibria (VLB) can be predicted for a large number of systems without introducing new model parameters that must be fitted to experimental VLB-data. The PSRK equation of state can be used for VLB-predictions over a much larger temperature and pressure range than the UNIFAC γ--approach and is easily extended to mixtures containing supercritical compounds. Additional PSRK parameters, which allow the calculation of gas/gas and gas/alkane phase equilibria, are given in this paper. In addition to those mixtures covered by UNIFAC, phase equilibrium calculations may also include gases like CH₄ C₂H₆, C₃H₆, c₄H₁₀, CO₂, N₂, H₂ and CO.     

Modelling the VLE properties using van der Waals-type equation of state

The vapour–liquid equilibrium (VLE) properties of polar and non-polar fluids have been modelled by the use of two modified van der Waals (vdW)-type equations of state (EOSs). In this article, a revised method is applied to the above-mentioned EOSs to improve the representation of VLE properties of different class of fluids. In this respect, the repulsion parameter b is considered to be temperature dependent and also a temperature-dependent revision factor α(T) is introduced to the liquid fugacity coefficient expression derived from traditional isothermal integration to reproduce the vapour pressure (P^s) of pure liquids. The present method is also extended to represent the VLE properties of binary mixtures containing noble gases, refrigerants and hydrocarbons. This method outperforms the original vdW-type EOSs in predicting the VLE and pressure-volume-temperature (PVT) properties of 22 pure substances and 7 binary mixtures.