ZnCl2熔盐的分子动力学模拟

近年来，若干作者根据几种粒子间函数，对ZnCl2熔盐结构做过分子动力学模拟［1－3］其出的Zn－Cl和Cl－Cl离子间的偏径向分布函数与中子衍射实测值符合较好，但Zn－Zn离子间距以及Zn－Zn间配位数计算值多偏高．鉴于Busing势函数在多价卤化物馆盐的分子动力学计算中应用效果较好[4],我们试用Busing势函数为基础对ZnCl2熔盐结构和能量做分子动力学计算．1研究方法计算所用粒子势为Busing势函数此处，Zi为离子的电行数（ZZn。＋＝2，Zcl－＝1），几；是离子有效半径，人为＊离子的“硬度”参数·据文献问，f二0．005071，尸zn。十二0．…     

小尺寸金属团簇熔化过程的分子动力学模拟

运用分子动力学方法模拟了小尺寸金属团簇的熔化过程, 原子之间的作用采用嵌入原子法(EAM)模型, 计算了均方根键长涨落δ随温度的变化, 以及升温过程中团簇热容的变化. 包含55、56个原子的面心立方(FCC)结构Au团簇的熔化过程是基本相同的. 而同样结构和数目Cu团簇的熔化过程却呈现出不同的趋势. Cu55、Cu56在模拟过程中都出现了FCC结构到二十面体结构的转变. 但由于表面多出了一个原子, Cu56的热容曲线比Cu55多了一个峰, 体系出现了预熔化现象. 这表明小尺寸团簇的固液转变的过程与团簇的原子类型、几何结构和原子数目密切相关.     

HTPB与Al不同晶面结合能和力学性能的分子动力学模拟

采用分子力学(MM)和分子动力学(MD)方法, 在250、300、350、400、450 K, 对固体推进剂端羟基聚丁二烯(HTPB)和铝晶胞不同晶面结构所组成的层模型在COMPASS力场下, 进行模拟计算, 求得结合能和静态力学性能(弹性系数、模量和泊松比). 模拟结果表明, 在400 K时HTPB与Al(011)面的结合能最大, 从综合力学性能优劣上看, 各个面从优到劣的排序为(011)>(221)>(001), HTPB与Al的结合能与力学性能具有对应关系, 结合能大的力学性能优异, 结合能小的力学性能较差.     

高分子/棒状纳米粒子复合物的分子动力学模拟

采用非平衡态分子动力学模拟研究了剪切场下棒状纳米粒子对高分子基体的结构、 动力学和流变性质的影响. 通过比较多种体积分数(0.8%~10%)的纳米复合物及纯熔体的模拟结果发现, 随着纳米粒子的增加, 高分子链的扩散和松弛逐渐受到抑制, 而链尺寸几乎保持不变. 从Weissenberg number(W_i)角度看, 在剪切流场下, 高分子链的结构性质(如归一化的均方回转半径、 回转张量和取向抑制参数)几乎与纳米粒子的体积分数无关, 而高分子链的Tumbling运动受到抑制. 研究还发现, 纳米复合物与纯熔体的剪切黏度曲线趋势基本一致, 即W_i=1将曲线分为平台区和剪切变稀区. 纳米棒的加入仅定量地改变了流体的剪切黏度.     

侧链液晶聚合物相转变点的分子动力学模拟

液晶基本上可分为３类：向列型液晶、近晶型液晶和胆甾型液晶［１］．液晶材料的相转变温度的预测对选择和设计液晶材料具有非常重要的作用．本文用分子动力学模拟法对聚［ω－（４′－甲氧基联苯－４－氧基）己基］甲基丙烯酸酯（ＰＭ６ＭＰＰ）单链分子进行了研究,模拟...     

单嘧磺隆晶体-活性构象转换的分子动力学模拟

应用分子动力学模拟方法对单嘧磺隆在水、正辛醇和正辛烷3种不同溶剂中的构象行为、单嘧磺隆与3种溶剂之间的相互作用能及氢键相互作用进行了计算研究. 计算结果表明, 在3种不同的溶剂中, 单嘧磺隆的优势构象不同; 其构象转换过程, 特别是转换成活性构象的过程主要发生在水溶液中; 与溶剂分子间的相互作用是分子构象行为的决定因素; 单嘧磺隆的脲桥部分可以和含氢键接受体的溶剂形成氢键, 分子间与分子内氢键的竞争可能是从晶体构象转换成活性构象的主要驱动力.     

基于粗粒化模型对有机溶剂的分子动力学模拟

在基于Boltzmann分布对四种基本构象进行Monte Carlo取样后, 通过与全原子模型的范德华势比较得到了Gay-Berne (GB)参数. 又在对用量化计算得到的分子体系的电势进行电荷、偶极矩和四极矩的拟合后, 得到了电多极展开势(EMP)参数. 利用得到的粗粒化参数, 基于粗粒化模型, 对CHCl₃及四氢呋喃(THF)两种有机溶剂进行了分子动力学模拟(MDS), 并将结果同全原子模拟进行了比较. 计算结果表明用粗粒化模型从整体上能重复全原子模型的模拟结果, 但在某些细节的计算与全原子模型有偏差, 其原因可能是目前工作仅考虑了单位点情况, 为此今后在对具有复杂结构的分子进行粗粒化模拟时还应考虑合理放置及增加相互作用位点.     

超临界水的分子动力学模拟

采用分子动力学(MD)模拟的方法对超临界条件下水的结构及扩散性质进行了研究.模拟结果表明超临界条件下水分子之间的氢键作用明显减弱,分子极性大大降低.扩散性质与常温下相比,其大小约上升了两个数量级.     

中介尺度Au纳米团簇熔化的分子动力学模拟

采用分子动力学模拟技术，研究了原子个数为16～8628的 Au纳米团簇的熔化过程.采用 Johnson的EAM (embedded atom method) 模型,模拟结果表明，金属纳米团簇存在一中介尺度区域.对Au纳米团簇而言，当原子个数N >456时，团簇的热力学性质与团簇尺寸呈线性关系，熔化首先从表面开始，逐步向中心区域推进，且满足Tmb－Tmc(N)＝aN(－1/3)的关系.另外，计算了中介区域的团簇的尺寸、熔化温度、表面能、熵、焓等热力学量以及均方根位移（RMSD）等动力学量，为研究纳米团簇提供定量数据.     

N-脒基脲二硝酰胺盐晶体形貌的分子动力学研究

采用分子动力学方法,在正侧(NVT)系综下研究了N-脒基脲二硝酰胺盐(FOX-12)在溶剂中的晶体形貌.通过构建溶剂分子层-晶面的界面吸附模型模拟其动力学平衡构型,计算溶剂与晶体表面间的结合能,进而对真空附着能进行修正并获得溶剂条件下的晶貌.使用自然冷却法在水和水/甲醇中培养FOX-12晶体并利用扫描电子显微镜进行了表征.结果表明,在真空条件下决定FOX-12晶貌的6个重要晶面为(110),(200),(201),(011),(002)和(111);FOX-12在水溶液条件下的主要晶面为(110)和(011),在水/甲醇溶液条件下的主要晶面为(200)和(011),预测的晶体形貌与实验结果相吻合.对水分子和FOX-12的(110)面间的径向分布函数进行了计算,分析了水分子和晶面间的分子间作用力.     

金属Cu液固转变及晶体生长的分子动力学模拟

采用分子动力学模拟研究了液态Ｃｕ在不同冷却速度下的凝固特点,模拟采用ＥＡＭ作用势,计算了不同温度,不同冷却速度下Ｃｕ的偶相关函数,结果表明ＥＡＭ作用势能很好地描述液态Ｃｕ的结构特征,当冷却速度较快时,液Ｃｕ形成非晶;当冷却速度较慢时,液Ｃｕ形成晶体,分析了不同冷却速度下体系的相变热力学及相变动力学过程,最后采用液固两层构型法,描述了Ｃｕ晶体的生长过程。     

表面上锚定聚乙烯链聚集的分子动力学模拟

利用分子动力学模拟方法研究了锚定于二维无限大表面上的聚乙烯分子链的聚集和有序化过程.聚乙烯链与表面的相互作用采用一个以平面法向距离为函数的势能表达式S(z).以均方回转半径垂直于表面的分量s²⊥作为指标,考察了模拟过程中表面势能函数的强度效应和模拟体系的温度效应.研究了锚定聚乙烯在表面上聚集的有序化过程,考察了锚定聚乙烯链的聚集成核的动力学过程,聚集体的形貌及分子量依赖性.发现锚定聚乙烯在300K时只在锚定点附近形成一个局部有序结构,然后以此为核进行生长.在模拟得到的有序程度不同的聚集体结构中,锚定点都位于聚集体链轴的端面.     

气体在无定型聚异戊二烯中扩散的分子动力学模拟

采用分子动力学模拟方法研究了多个温度下氧气、氮气及甲烷在无定型顺式1,4-聚异戊二烯中的扩散系数。在模拟过程中,使用COMPASS力场作为分子力场。应用COMPASS力场的势能函数,聚合物的密度及玻璃化转变温度的计算结果与实验值有较好吻合。在278-378 K的温度范围内,通过3或1.5 ns时长的正则系综动力学模拟,计算了不同温度下氧气、氮气及甲烷的扩散系数。结果表明,根据爱因斯坦关系式计算得到的扩散系数与实验结果比较接近。对气体扩散系数与温度的关系进一步研究,发现在278-378 K温度范围内,甲烷的扩散系数随温度变化的半对数曲线图是非线性的,而氧气和氮气的扩散系数随温度变化的半对数曲线图是线性的。本文研究结果有助于理解温度对气体扩散的影响机制,并为高温下气体在天然橡胶中扩散系数的测定及天然橡胶热氧老化建模分析提供依据。     

贵金属Au的液态结构分子动力学模拟

用分子动力学模拟方法在1573－200K的温度范围内对液态Au的微正则系综进行了模拟研究。模拟采用或嵌原子相互作用势对时间和空间的平均，得到了不同温度下Au的润分布函数及原子组态变化的重要信息，并利用键对分析技术对模拟结果作了深入讨论．     

高压下MgO熔化特性的分子动力学模拟

Molecular dynamics simulation was used to study the melting of MgO at high pressures. The melting temperature of MgO was accurately obtained at elevated temperature and high pressure after corrections based on the modern theory of melting. The calculated melting curve was compared with the available experimental data and other theoretical results at the pressure range of 0-135 GPa. The corrected melting temperature of MgO is in good agreement with the results from Lindemann melting equation and the twophase simulated results below 15 GPa.     

高分子玻璃化转变过程的分子动力学模拟本科教学实验

玻璃化转变是高分子专业教学中的重要内容,转变过程中自由体积的变化和分子链段的运动是理解玻璃化转变的难点。在计算机上使用分子模拟的方法得到三维的可视化图形,可以更直观地观察到自由体积和分子链的变化。提炼近年分子模拟技术在玻璃化转变研究中的最新科研成果,得到既有理论基础又适合本科教学的课程素材。分子模拟可以实现了仪器测试无法达到的超快速升降温,并与仪器测试得到的实验结果相对照,验证了普通实验很难验证的理论观点,从而拓展了本科实验的范围,达到了很好的教学效果。     

受限状态聚合物熔体的分子动力学模拟

用粗粒化的分子动力学(MD)模拟方法从分子层次研究了受限于粗糙壁内的聚合物熔体的动力学性质. 结果表明, 对于链长较短的受限聚合物熔体体系, 随着膜厚的增加, 体系内部高分子链的松弛时间逐渐减少; 然而对于链长较长的受限体系, 聚合物链的松弛时间随着膜厚的增加先减少后增加. 推测这种由于链长的变化所引起的动力学性质的差异源自受限熔体内聚合物链聚集状态的改变, 并且通过考察交叠参数对这种改变进行了分析. 结果表明, 在膜厚增加的过程中, 决定受限状态高分子长链松弛机理的因素逐渐从受限效应转变成为链间的缠结效应.     

Molecular Dynamics Simulation Studies on the Cooling Process of Polyvinyl Chloride

IntroductionRecently,molecular dynamics(MD)simulationhas been used extensively in the study of polymer sys-tems[1—11].The increased interest in computer simula-tion is assisted by the fast evolution of computer tech-nology,which makes molecular simulatio…     

用分子动力学模拟水合物储氢

用分子动力学(MD)模拟方法系统研究了结构Ⅰ型(SⅠ)和结构H型(SH)氢气水合物中氢气的占据情况并确定了氢气水合物的稳定结构: SⅠ水合物氢气分子数小胞中为2, 大胞中为3; SH水合物氢气分子数小胞中为2, 中胞中为2, 大胞中为11. 分析了稳定情况下水合物各胞腔内氢气分子之间的径向分布函数(RDF), 得出了氢气分子在各胞腔内的稳定位置. 由稳定位置得到了稳定结构下氢气水合物的储氢质量分数: SⅠ为5.085%, SH为6.467%. 与实验对比得出结论: SH水合物稳定结构下的储氢能力最强.     

聚乙烯链在碳纳米管侧壁吸附的动力学模拟研究

利用经典的分子动力学模拟方法对聚乙烯(PE)分子在两种不同类型的碳纳米管(CNT)中的吸附进行了研究. 计算了两者之间的扩散系数和相互作用能; 利用PE链自身的扭转角分布和取向参数对PE链构象进行了分析. 结果表明, PE链可以在CNT上很好的吸附, 且PE的构象和吸附位置主要与温度和CNT的半径有关, 与管的类型关系不大.