AB2型聚合物流体的表面结构性质

在密度泛函理论(DFT)框架下, 应用改进的基本度量理论(MFMT)表达硬球作用对自由能泛函的贡献, 根据统计力学理论结合加权密度近似(WDA)表达聚合作用对自由能泛函的贡献, 建立了描述AB2型聚合物流体的化学势, 得到了聚合物流体在硬球颗粒表面的密度分布表达式, 计算了聚合物流体在硬球颗粒表面附近的密度分布, 并探讨了体积分数、聚合程度和硬球颗粒尺度对体系密度分布的影响. 此外, 通过体系密度分布, 进一步分析了体积分数、聚合程度和硬球颗粒尺度与剩余吸附的关系.     

平行硬壁间荷电硬球流体的密度泛函理论研究

基于经典流体的密度泛函理论并结合改进的基本度量理论, 研究了受限于平行硬壁间的荷电硬球流体的平衡密度分布. 通过对比有、无外电场时分子数密度分布的相应变化情况, 分析了荷电硬球分子间的库仑排斥势能与场致势能间的竞争作用, 研究了电场对于受限荷电硬球流体系统聚集态结构的影响.     

约束条件下的硬球流体

利用密度泛函理论和分子动力学方法,对处于两平行硬墙之间的硬球流体的密度分布进行了计算通过比较两种方法的结果,发现在墙之间距离较大时,Rosenfeld密度泛函理论的结果与分子动力学模拟的结果符合很好;当两堵墙间的距离很小时,这两个结果之间存在明显的不一致.另外,还研究了约束条件下密度分布的结构.     

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

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

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

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

羟基化氧化锌催化甲醇裂解反应机理的理论研究

采用密度泛函理论的B3LYP方法在两种基组水平上研究了羟基化氧化锌催化裂解甲醇反应机理, 优化了反应过程中各反应物、中间体、过渡态和产物几何构型, 并且应用了自然轨道理论(NBO)和分子中的原子理论(AIM)分析了这些物质的成键特征和轨道间的相互作用. 研究中分别采用了羟基化一个ZnO和羟基化两个ZnO作为催化剂的两种简化模型. 结果表明两种不同基组研究所得的反应历程及相应能量变化趋势一致; 羟基化两个ZnO催化剂活化裂解甲醇反应更容易发生.     

氢键流体中Janus粒子的过量熵

Janus粒子在结构和性质上的特殊性使其在不同溶剂中呈现出丰富的聚集态结构和相态结构,这与其在溶液中的过量熵密切相关.本文以Janus粒子与氢键流体的稀溶液体系为例,利用经典流体的密度泛函理论研究了氢键流体中Janus粒子的过量熵.首先在极稀条件下给出Janus粒子外氢键流体中的局域密度分布,由此得到Janus粒子与氢键流体的二体分布函数,进而计算了不同情况下的过量熵.在此基础上,阐明了Janus粒子与氢键流体间的相互作用、氢键流体的体相密度、氢键强度和氢键官能度等因素对过量熵的影响.本研究旨在揭示这些因素对过量熵的调控机制,在定量水平上明确Janus粒子在溶液中组装的驱动力,从而为深入研究其聚集态结构提供可借鉴的理论线索.     

镁锡水滑石中的超分子作用

通过构建镁锡水滑石(Mg3Sn-LDHs-yH2O)周期性计算模型,基于密度泛函理论,用CASTEP程序模块,进行模型的几何全优化.对各体系的结构参数、Mulliken电荷布居、态密度(DOS)、逐级水合能等进行分析,探究客体CO320-和H2O的分布形态及其与主体层板的超分子作用.结果表明,[Mg6Sn2(OH)16]4+层间插入客体阴离子CO320-和水分子后,主客体间存在着较强的超分子作用力,主要包括静电和氢键作用,且氢键作用在水合过程中起主导作用,总体上层板-水(L-W)型/层板-阴离子(L-A)型氢键强度要强于阴离子-水(A-W)型/水-水(W-W)型氢键.随着层间水分子的增多,层间距先增大后又稍降低.当y=0、1时,客体所在的平面与主体层板平行,且与两层板的距离基本相等;当y=2、3时,客体以偏向某一层板的形式存在.与层间H2O相比,层间阴离子CO230-对体系态密度的影响更显著,层板与CO320-的总体作用力大于与H2O的总体作用力.随着水分子数的增加,Mg3Sn-LDHs-yH2O体系的逐级水合能绝对值逐渐降低,说明Mg3Sn-LDHs的水合程度不会无限增加,而是具有饱和量.     

缔合流体在纳米缝隙中的相行为研究

结合改进的基础度量理论(modified fundamental measure theory, MFMT)和密度泛函理论(Density functional theory, DFT), 研究了二缔合Lennard-Jones(LJ)流体在纳米缝隙中的相行为. 根据平衡时两相温度、化学势及巨势相等的原则, 计算了二缔合LJ流体在纳米缝隙中的相平衡. 阐明了缔合能量和缝隙宽度对平衡时分子和单体密度分布、吸附-脱吸等温线及汽液共存等热力学性质的影响.     

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

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

Software Package: An Advanced Theoretical Tool for Inhomogeneous Fluids(Atif)

In spite of the impending flattening of Moore’s law,the complexity and size of the systems we are interested in keep on increasing.This challenges the computer simulation tools due to the expensive computational cost.Fortunately,advanced theoretical methods can be considered as alternatives to accurately and efficiently capture the structural and thermodynamic properties of complex inhomogeneous fluids.In the last decades,classical density functional theory(cDFT)has proven to be a sophisticated,robust,and efficient approach for studying complex inhomogeneous fluids.In this work,we present a pedagogical introduction to a broadly accessible open-source density functional theory software package named"an advanced theoretical tool for inhomogeneous fluids"(Atif)and of the underlying theory.To demonstrate Atif,we take three cases as examples using a typical laptop computer:(i)electric double-layer of asymmetric electrolytes;(ii)adsorptions of sequencedefined semiflexible polyelectrolytes on an oppositely charged surface;and(iii)interactions between surfaces mediated by polyelectrolytes.We believe that this pedagogical introduction will lower the barrier to entry to the use of Atif by experimental as well as theoretical groups.A companion website,which provides all of the relevant sources including codes and examples,is attached.     

甲醇钠引发的环氧乙烷开环聚合反应过程

运用密度泛函理论(DFT)的Dmol3方法, 计算了甲醇钠引发的环氧乙烷开环聚合的反应过程. 并运用前线轨道理论对该聚合反应的各步反应历程进行了分析. 计算结果表明, 链引发为无能垒的放热反应, 放出的能量达到92.560 kJ·mol-1, 而链增长过程则需越过100.951 kJ·mol-1的反应能垒, 链增长物种与环氧乙烷的前线轨道相对称, 可以使开环聚合反应继续进行下去. 当向反应体系中加入草酸、磷酸等质子酸时, 会立即发生链终止反应. 此外, 还对链增长过渡态的合理性进行了确认, 绘出了相应的反应势能曲线.     

Computational approaches to design a molecular imprinted polymer for high selective extraction of 3,4-methylenedioxymethamphetamine from plasma

In this work, a molecular imprinted polymer (MIP) as a novel selective sorbent for extraction of 3,4-methylenedioxymethamphetamine (MDMA) from plasma samples was prepared. For selecting a more suitable monomer and polymerization solvent a methodology based on density functional theory calculations was developed. This computational design is based on the comparison of stabilization energies of the prepolymerization adducts between the template and different functional monomers. The effect of polymerization solvent was studied using of polarizable continuum model (PCM). The computational results revealed that the best suitable monomer and polymerization solvent for preparation of MIP is methacrylic acid (MAA) and chloroform, respectively. Also, another MIP with methacrylic acid (MAA) as monomer in acetonitrile was prepared to evaluate the validity of polarizable continuum model for selection of polymerization solvent. The performance of each polymer was evaluated by using Langmuir-Freundlich (LF) isotherm. As it is expected, the best results were obtained for the MIP which was prepared in chloroform. This MIP was used as a selective sorbent in solid-phase extraction coupled with high performance liquid chromatography-ultraviolet detector (MISPE-HPLC-UV) for rapid screening of MDMA in human plasma. For the proposed MISPE-HPLC-UV method, the linearity between responses (peak areas) and concentrations was found over the range of 3.6-11500 ng mL(-1) with a linear regression coefficient of 0.998. The limit of detection (LOD) and quantification (LOQ) in plasma were 1.0 and 3.3 ng mL(-1), respectively. The %RSD (n=5) data for five plasma samples containing 15, 25, 50, and 100 ng mL(-1) of MDMA were 1.02, 1.12, 2.05, 2.54, respectively.     

Polymerization Mechanism of α-Linear Olefin

The density functional theory on the level of B3LYP/6-31G was empolyed to study the chain growth mechanism in polymerization process of α-linear olefin in TiCl₃/AlEt₂Cl catalytic system to synthesize drag reduction agent. Full parameter optimization without symmetryrestrictions for reactants, products, the possible transition states, and intermediates wascalculated. Vibration frequency was analyzed for all of stagnation points on the potential energy surface at the same theoretical level. The internal reaction coordinate was calculated from the transition states to reactants and products respectively. The results showed as flloes:(i) Coordination compounds were formed on the optimum configuration of TiCl₃/AlEt₂Cl.(ii) The transition states were formed. The energy di?erence between transition states and the coordination compounds was 40.687 kJ/mol. (iii) Double bond opened and Ti-C(4) bond fractured, and the polymerization was completed. The calculation results also showedthat the chain growth mechanism did not essentially change with the increase of carbon atom number of α-linear olefin. From the relationship between polymerization activation energy and carbon atom number of the α-linear olefin, it can be seen that the α-linear olefin monomers with 6-10 carbon atoms had low activation energy and wide range. It was optimum to synthesize drag reduction agent by polymerization.     

硅烷类外给电子体对丙烯-丁烯序贯聚合的调控作用

外给电子体(ED)作为负载型Ziegler-Natta催化剂的一个重要组分, 在影响α-烯烃的催化活性及聚合物的立构规整性方面发挥着重要作用. 本文研究了4种不同结构及电子密度的硅烷类外给电子体[二甲基二甲氧基硅烷(D1)、 二丁基二甲氧基硅烷(D2)、 二苯基二甲氧基硅烷(D3)及二环戊基二甲氧基硅烷(D4)]对丙烯均聚及丙烯(一段)-丁烯(二段)序贯聚合的影响. 结果表明, ED对烯烃聚合的催化活性、 活性中心数及活性中心定向能力都具有显著的影响. 密度泛函理论(DFT)模拟计算表明, 随着ED的空间位阻和电子密度增加, ED在MgCl₂表面的吸附能降低, 吸附稳定性降低; ED的空间位阻和电子密度增加有利于提高丙烯聚合活性中心的定向能力, 当n(D4)/n(Ti)=20时, 合成的聚丙烯(PP)中的等规聚丙烯(iPP)组分含量达到92.8%. 当n(ED)/n(Ti)=15时, 丙烯聚合的聚合速率常数达到最大值; 具有更大空间位阻和电子密度的ED使得丙烯-丁烯序贯聚合的活性中心具有更强的定向能力, ED对丁烯(二段)聚合活性及聚丙烯/聚丁烯合金（PBA）中等规聚丁烯(iPB)组分的熔点影响更显著.     

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.     

Density of methanol-carbon dioxide mixtures at three temperatures: Comparison with vapor-liquid equilibria measurements and results obtained from chromatography

The density of methanol-carbon dioxide mixtures has been measured at 40, 50, and 60°C and the data presented as constant density curves in plots of pressure against composition. Mixtures containing from 1 to 23.5% methanol were covered. The density information was compared with liquid-vapor equilibria data from the literature and with results from chromatographic separations using both sub- and supercritical conditions. The results indicate that phase transitions do not usually occur when the concentration of the modifier is increased while holding temperature and pressure constant. Similarly, no phase transition occurs when the temperature is increased from 28 to 80°C, at a constant (high) pressure and modifier concentration. Binary fluids with many desirable compositions cannot be prepared at low densities: densities below 0.4–0.5 g/cm³ cannot be obtained from methanol-carbon dioxide mixtures because the fluids separate into two phases. This means that density programming with many binary fluid combinations will be of secondary importance compared with composition programming.     

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

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