Al液诱导凝固过程的分子动力学模拟

采用分子动力学(molecular dynamics,MD)方法模拟了嵌入不同半径(0.37～2.4 nm)固态Al晶核的Al液诱导凝固过程.结果表明,体系临界温度和嵌入晶核半径的倒数呈线性关系,符合吉布斯-汤姆森(Gibbs-Thomson,G-T)效应,得到的G-T系数Γ为1.4×10~(-7)K·m,体材料熔点T_m~(bulk)为(985.36±11.25)K.借助Γ计算出固液界面能为(140.35±9.05) mJ/m~2,与通过毛细波动法获得的149 mJ/m~2比较接近,再次证实了Turnbull的实验值(93 mJ/m~2)较低.外推得到诱导凝固失效半径为0.91 nm,此时体系温度已达到临界温度极限值.在该极限值下,体系可以自发形核,孕育时间呈现随机性,微观结构表现为不稳定的层错互相交织.在有效晶核尺寸范畴内,孕育时间随嵌入晶核半径的增大而增加,长大速度却随半径的增大而降低,微观结构是稳定的lamellar结构.     

分子动力学模拟固态Al的自扩散性质

采用EAM势,用分子动力学模拟方法研究了固态Al的自扩散,用2种方法给出了Al在固态下较长的温度区域内的自扩散系数以及自扩散系数随温度的变化规律,并从微观机制上进行了分析和解释,结果较好地符合Arrhenius关系。     

氩粒子系统的分子动力学模拟

分子动力学方法是一种计算机模拟实验的方法,这种方法不仅可以得到原子的运动轨迹,还可以观察到原子运动过程中的各种微观细节。通过对256个氩原子系统进行分子动力学模拟得出了氩原子系统的能量演化过程、粒子的运动轨迹、系统的径向分布函数等一些有意义的结果。     

模型溶液的分子动力学模拟及扩散系数计算

通过对模型溶液的分子动力学模拟，确定了这些溶液的径向分布函数，自扩散系数Ｄ１、Ｄ２和互扩散系数Ｄ１２。结果表明，用Ｅｉｎｓｔｅｉｎ法和Ｇｒｅｅｎ－Ｋｕｂｏ法得到的扩散系数在数值上是一致的；溶液互扩散系数Ｄ１２与自扩散系数Ｄ１和Ｄ２满足关系式Ｄ１２＝ｘ１Ｄ１＋ｘ２Ｄ２。     

水合镁离子的分子动力学模拟研究

基于水合离子[M(H2O)m]n+的概念,利用Moldy分子动力学模拟并结合分子内相互作用势MCY来研究水中的镁离子.该模拟系统包括一个二价的镁离子和216个水分子,在温度为330K下,计算了系统内不同原子对之间的径向分布函数(RDF),模拟结果和实验值基本一致.而且从其径向分布函数可看出,镁离子的出现并没有使得镁离子周围的第一水化层水分子和第二水化层水分子之间的距离有明显的减少.     

磁场对水内能作用的分子动力学模拟

将水系统视为正则系统，采用分子动力学模型方法计算了磁场条件下水的内能，比热和径向分布函数，结果表明，磁场可以影响径向分布函数g(r)的分布，使水的结构发生变化，从而导致水的内能，比热发生改变，水的内能，比热，径向分布函数伴峰的高度与磁场强度的关系均呈多极值特征，且在常温下磁感应强度为0．25T时，磁场作用最明显。     

受限Lennard-Jones流体自扩散系数的分子动力学模拟

基于平衡态分子动力学(EMD)方法,建立了受限空间中的Lennard-Jones(LJ)流体自扩散模型.采用径向分布函数对LJ流体微观结构进行了表征,模拟了LJ流体在纳米尺度受限空间中的自扩散系数,并将其与相应的自由空间内LJ流体自扩散系数进行了比较,同时从分子水平分析了温度、密度和受限尺度对自扩散系数的影响.研究结果...     

甲醇/水混合物的分子动力学模拟

用分子动力学模拟研究甲醇／水混合物的结构性质和动力学性质．模拟中甲醇分子采用点点刚性模型,水分子采用TIP4P模型,计算了不同浓度甲醇／水混合物甲醇分子的O-O和H-H径向分布函数,配住数以及自扩散系数．从径向分布函数,配位数和自扩散系数发现,甲醇分子间的结构随着甲醇浓度的增加而加强．     

液态Cu等温凝固的分子动力学模拟

使用Tight-binding势函数,对液态Cu在等温凝固过程中的结构变化进行了分子动力学模拟(MD模拟)计算,得到体系在不同温度下的双体分布函数和配位数分布等静态结构信息,对等温凝固过程中FCC短程有序结构可能发生的变化以及由此导致的H-A键型变化进行了分析,并结合键对分析方法计算了不同弛豫时间下典型短程有序结构的分布. 计算表明,在Cu凝固结晶相变过程中1551键应是先向1541键转化,初始三维结构的形成可能主要依赖于Cu原子在两个方向上的扩散和弛豫.     

电流变液的分子动力学模拟

运用分子动力学方法研究了电流变液的在外电场下的液－固相变过程,原子间的相互作用采用简单偶极势来描写,计算结果表明：对于一个给定的外电场,当系统温度逐步降低的时间,电流变液中的颗粒首先会沿电场方向排列成链状结构,然后逐步转变成块状体四方结构。研究结果与作者前期用统计力学方法所得结论相一致,也与实验结果比较吻合。     

分子尺度下的水损沥青静态裂缝自愈合

为探寻水分对沥青自愈合的影响机理,采用分子动力学方法从分子角度对比分析常规沥青、水损沥青的自愈合过程,并通过宏观试验进行验证。首先,结合沥青体系能量变化,将愈合状态分为3个阶段：分子迁移低密度阶段、裂纹愈合及修复阶段、分子自由运动阶段;然后,运用径向分布函数和均方位移函数分析2种沥青在裂纹宽度、含水率、沥青组分的扩散机理和聚集状态,解释沥青的自愈合及水对沥青自愈合的影响机理;最后,进行宏观试验验证。结果表明：在裂缝自愈合过程中,水分将阻碍沥青分子间的作用力,削弱沥青体系能量,导致沥青扩散力降低,且含水量越多,扩散力下降越明显;裂缝减弱了沥青各组分的关联性,破坏了沥青质胶团和胶质过渡相的稳定性,且裂缝处的分子数目随着裂缝尺寸的增加而减少,促使沥青自愈合能力降低。采用分子动力学从分子角度解释沥青愈合机理的可靠性,有助于了解沥青分子在水损过程中结构形态和扩散行为的变化特征,为未来水损沥青自愈合的研究提供参考。     

温度和外加电场对水分子团簇影响的分子动力学模拟

应用分子动力学（MD）模拟方法研究了温度从273～373K，电场强度从0-2×10^10V／m时水分子团簇的热力学性质、微观结构和动力学性质．研究发现，内能与温度成正比与电场强度成反比，通过分析径向分布函数，我们得知，温度升高使水的氢键作用减小，外加电场增加使氢键作用增强，使水的有序程度加强．水的自扩散系数随温度升高而增大，随电场强度的增加而减小，且在低温时，有电场时水的自扩散系数是无电场时的十分之一．     

分子动力学模拟烷烃链条数对其在氧修饰石墨烯纳米孔中迁移的影响

为了分析高分子在纳米孔中迁移问题,更好地为纳米器件的设计提供理论支持,采用分子动力学方法模拟了一系列链条数不同的烷烃分子在氧修饰石墨烯纳米孔中的钻孔,研究了链条数N对钻孔条数的影响,给出了均方位移MSD和扩散系数D随N变化的规律.研究结果表明:穿孔条数随着链条数的增加上升明显,但当链条数多到一定程度后,穿孔烷烃数较为缓慢地下降.     

LiCl熔盐分子动力学模拟的Delaunay单纯形分析

利用Delaunay单纯形理论,对以LiCl为代表的熔盐系的动力学模拟构型进行了众多参数的实现及其统计.重点分析了熔化过程中Delaunay单纯形体积、面积、四面体系数及其Kirie单纯形随温度变化的趋势.借助逾渗理论和染色方法,熔化的过程可以认为是大面积Delaunay单纯形出现、成串直至逾渗的过程,而液体流动性的根源则是大体积Delaunay单纯形的出现和增加.同时,四面体系数的分布进一步证实Kirie单纯形的变化趋势.     

经激光束处理的金属Fe的分子动力学研究

采用分子动力学和EAM势函数,建立了金属Fe固液两相模型,实现了液态金属Fe在两种不同冷却速率下的凝固过程,用以模拟激光束对金属Fe的处理。通过原子势能、径向分布函数、共同近邻分析和原子构型演化等计算方法,研究了体系演化过程和终态结构特性,得到了在冷却速率为2 K/ps和5 K/ps时固液分界面稳定传播速度分别为1.74 m/(s·K)和0.89m/(s·K);在终态时由于无序原子的存在,原固液分界面处存在能量起伏,同时在原液态部分发现了高能态的原子团簇,为研究经激光束处理后的金属Fe的力学性能提供了微观结构基础。     

瓦斯流动的分子动力模拟中势函数的研究

势能函数作为分子动力模拟的重要问题,它的选取及求解成为这一研究的热点．探讨了分子动力模拟中所涉及的一系列势能函数,选取了所研究的煤层瓦斯流动分子动力模拟所需要的势函数,并根据临界常数法求取了此函数所涉及到的势能参数。最后计算出一个甲烷分子与截断半径处另一分子之间的L-J势能,为后期的模拟研究奠定了基础。     

超薄膜润滑的分子动力学模拟

超薄膜（膜厚趋于分子量级）的摩擦特性与宏观流体膜有很大的不同,与超薄膜的微观结构有密切的关系。本文应用可以同时模拟超薄膜宏观和微观特性的分子动力学模拟（ＭＤＳ）方法,研究了超薄膜的微观结构与摩擦学特性间的关系,发现了平行于壁面的层状类固态结构;固液作用强度及膜厚大小对类固态结构有着明显的影响;超薄膜的极限膜厚由类固态与液相的比例决定。这种微观结构的相变改变了超薄膜的摩擦学特性。模拟中还发现了剪切诱导的微观构型。     

Molecular dynamics simulation of kink in 〈100〉 edge dislocation in body centred cubic iron

Using the molecular dynamics method,we have constructed two kink models corresponding to the 〈100〉{010} and 〈100〉{011} edge dislocations (EDs) in body centred cubic (bcc) Fe. It is found that the geometric structure of a kink depends on the type of edge dislocation and the structural energies of the atoms sites in the dislocation core region. The formation energies,migration energies and widths of the kinks in different types of EDs are calculated. The results show that formation and migration of the kink in the 〈100〉{010} edge dislocation are difficult. The 〈100〉{011} edge dislocation moves primarily through kink nucleation,rather than kink migration.