应力诱发NiAl单晶马氏体相变的分子动力学模拟

利用嵌入原子势（EAM）,对NiAl单晶在外应力作用下的动态拉伸过程进行了分子动力学模拟.应力-应变曲线分析以及原子构型分析表明外应力诱发NiAl合金发生了马氏体相变,原子结构由B2相转变为L1₀相.通过研究原子构型的演化过程,发现马氏体相变是通过多个{110}孪晶面的扩展和湮灭作用来完成的.同时探讨了马氏体相变的微观机理. 关键词： 马氏体相变 NiAl 分子动力学模拟     

铁的冲击相变机制的分子动力学研究

 用分子动力学方法模拟计算了在冲击波加载条件下,单晶铁中的结构相变（由体心立方结构α相到六角密排结构ε相）,相互作用势采用铁的嵌入式原子势（EAM）,单晶铁样品的尺寸为28.7 nm×22.9 nm×22.9 nm,总原子数为1.28×10⁶个。通过推动一个运动活塞对静止靶的作用来产生冲击压缩,加载方向沿单晶铁的[100]晶向。通过对原子位置的追踪,揭示了铁的冲击相变机制,计算结果表明相变机制包括两步：首先是在{011}面上的原子受到沿〈100〉晶向的压缩,使{011}面转化成正六角形密排面;然后是在{011}面上原子沿〈0-11〉晶向的滑移,完成由bcc结构到hcp结构的相变。同时发现滑移面只出现在与冲击波加载方向平行的(011)和(0-11)面上。     

铁冲击相变的分子动力学研究

用分子动力学方法模拟了单晶铁(Fe)在一定初始温度下冲击相变(α相→ε相)的微观过程，结果显示温度会导致冲击相变压力阈值降低.基于此微观过程，对加卸载波系的传播规律进行了相应计算和分析，结果表明在卸载过程中逆相变波(ε相→α相)相对于波前以当地纵波声速传播，而相对波后以亚声速传播，这可由卸载压力-密度曲线给出相应解释；计算了不同初态的卸载压力-密度状态曲线，并给出了逆相变带的分布，其分布规律显示了卸载过程逆相变的滞后现象. 关键词： 分子动力学 多体势 冲击波 相变     

热峰作用下单斜ZrO2相变过程的分子动力学模拟

ZrO2陶瓷耐高温、耐腐蚀、抗辐照性能强,是极具前景的反应堆惰性基质燃料和锕系元素固化材料.本文联合使用热峰模型和分子动力学方法,模拟了核辐射环境下ZrO2的相变过程:基于热峰模型,从快速重离子注入后能量沉积和传导的多物理过程出发,建立热扩散方程,求得ZrO2晶格温度时空演变特性;然后运用分子动力学方法模拟了该热峰作用...     

考虑非傅立叶效应的固液相变分子动力学模拟

采用耦合双温度模型的分子动力学方法对飞秒激光照射金箔的固液相变过程进行了模拟研究,利用序参数法对固液原子进行判定从而确定了金箔发生相变时的固液界面位置和温度,对基于傅立叶定律的抛物线模型和考虑非傅立叶效应的双曲线模型模拟得到的结果进行对比分析,在此基础上采用耦合双曲线模型的分子动力学方法研究了激光能量密度和脉冲宽度对金箔相变过程的影响.结果表明,当激光作用于金箔时,金箔上表面首先熔化,固液界面随时间不断向金箔底部移动,并且在相同条件下,双曲线模型下的金箔熔化深度和固液界面温度均大于抛物线模型的结果.当考虑非傅里叶效应时,激光能量密度越大,固液界面温度越高,金箔熔化时间越短;激光脉冲宽度越小,固液界面温度越大,金箔熔化速度越快.     

铝纳米颗粒的热物性及相变行为的分子动力学模拟

采用分子动力学方法模拟了纳米金属铝在粒径为0.8-3.2 nm 时的熔点、密度和声子热导率的变化, 研究了粒径为1.6 nm的铝纳米颗粒的密度、比热和声子热导率随温度的变化. 采用原子嵌入势较好地模拟了纳米金属铝的热物性及相变行为, 根据能量-温度曲线和比热容-温度曲线对铝纳米颗粒的相变温度进行了研究, 并利用表面能理论、尺寸效应理论对铝纳米颗粒熔点的变化进行了分析. 随着纳米粒径的不断增大, 铝纳米颗粒的熔点呈递增状态, 当粒径在2.2-3.2 nm时, 熔点的增幅减缓, 但仍处于递增趋势. 随着纳米粒径的增大, 铝纳米颗粒的密度呈单调递减, 热导率则呈线性单调递增, 且热导率的变化情况符合声子理论. 随着温度的升高, 粒径为1.6 nm的铝纳米颗粒的密度、热导率均减小. 该模拟从微观原子角度对纳米材料的热物性进行了研究, 对设计基于铝纳米颗粒的相变材料具有指导意义.     

高密度氦相变的分子动力学研究

采用分子动力学方法结合对关联函数分析计算了0–1000 K范围内氦的固–液相变曲线, 与实验数据的对比显示, 在0–500 K之间与实验数据符合很好, 500 K以上还没有相应的实验数据. 另外, 计算了钛金属中不同尺寸氦泡的压强, 并与高密度氦的固–液相变曲线进行了对比. 结果显示, 在低温条件下, 随着温度的降低, 钛晶体中可能会出现固态氦泡; 在300 K以上不会存在固态氦泡.     

多级孔碳负载聚乙二醇相变复合材料的分子动力学模拟研究

多级孔碳(Hierarchical porous carbon,HPC)具有大的比表面积(1345 m~2/g)和高的孔体积(2.69 cm~3/g),对聚乙二醇(Polyethylene glycol,PEG)的负载量高达92.5% (wt),复合材料展现出优异的储能特性,因此采用分子动力学模拟的方法对其开展机理研究。复合后PEG@HPC熔点为324.5 K,相互作用能和径向分布函数结果表明,孔径越小,基材对芯材的作用力越强,有利于防止芯材的泄漏,保持复合材料的结构稳定性。振动态密度和重叠能表明,复合后声子振动匹配增强,有益于能量的传递和材料热导率的提升。     

高温高压下NaCl热力学特性的分子动力学模拟

欲了解固体物理和地球内部特征就要弄清固体在高温高压的热力学特性，而碱金属卤化物的研究有利于在实验室难以接近的高温高压条件下建立金属氧化物的相转变模型。文中应用分子动力学方法模拟NaCl离子晶体的高压熔化、高温高压下热力学参数，如热容、等压体膨胀系数、等温体模量、Gruneisen参数等。     

水导热系数的分子动力学模拟

本文选用平衡分子动力学模拟方法,应用不同的势能模型对强极性分子水的导热系数进行了模拟计算.结果表明,用分子动力学模拟方法计算水的导热系数时,模拟结果对模型的选取极为敏感.在目前应用较为广泛的几种模型中以TIP4P模型较为适用,所得的导热系数与实验值较为接近.     

扩散系数的分子动力学模拟

本文选用平衡分子动力学模拟方法,应用不同的势能模型对强极性分子水的扩散系数进行了模拟计算．结果表明,用分子动力学模拟方法计算水的扩散系数时,模拟结果对模型的选取极为敏感。在目前应用较为广泛的几种模型中以 SPCE模型较为适用,所得的扩散系数与实验值较为接近。     

多核环境下的分子动力学模拟

本文在多核环境下,使用OpenMP实现了经典分子动力学模拟程序的并行;同时对分子动力学模拟进行了两项主要的优化:分子排序及运用SIMD指令运算.在4核下获得了4.13倍的计算性能提升,将经典分子动力学模拟的模拟规模提高至4000分子×10~7模拟总步数.     

纤维素热解过程的分子动力学模拟

本文利用 PM3 方法对生物质主要组分纤维素进行了结构优化并得到了一系列的结构参数.对不同力场下,聚合度为 9 的纤维素单链热分解进行了分子动力学方法模拟研究,得到不同力场下的模拟过程中参数.通过与相关生物质热解实验结果对比,利用 Amber 力场模拟得到的结果与实验值吻合较好.基于 Amber 力场并结合量子力学对纤维素单元热解过程进行了研究,模拟得到纤维素单元分子链在加热过程中的主要分解温度范围、断键顺序以及一次热解的基团,并对一次产物进行了分析.     

分子动力学模拟粒子在表面上的俘获过程

采用分子动力学方法研究了粒子在金属衬底上的俘获过程。发现粒子垂直入射时,俘获概率随入射能量的增加而减小,而45°入射时,俘获概率表现出非单调行为,这些结果与实验观察相一致。同时还研究了俘获概率对入射角的依赖。     

高密度氦状态方程的分子动力学模拟研究

 用分子动力学模拟和与固氦实验等温线比较的方法,确定了高密度氦指数6势函数中劲度系数α的优选值为12.7。然后我们用这个新的α值计算了fcc固氦的状态方程及其径向分布函数,发现当把ρ限定为1.60 g/cm³时,其径向分布函数的第二个峰将在2 000 K到3 040 K区间消失,表明此时固氦的长程有序性降低,或者说发生了固-液相变过程。     

铂薄膜导热特性的分子动力学模拟

采用非平衡态分子动力学(NEMD)方法模拟分析了纳米铂(Pt)薄膜的导热性能与脉冲激光作用下的温度响应特性.结果表明,100～500 nm铂薄膜的法向导热系数比体材料值低很多,而且低于其面向导热系数;微米铂薄膜的温度响应时间在纳秒量级;在脉冲加热的初始阶段,有一快速非傅立叶热波沿铂薄膜的厚度方向传递.     

分子动力学模拟铜薄膜的热导率

采用分子动力学(MD)方法模拟铜薄膜的热导率,给出了厚度在100~400nm、温度在100~600K范围内铜薄膜热导率对尺寸及温度的依赖关系.     

分形理论在分子动力学模拟中的应用

本文基于分形理论提出了一个假说,认为实际流体分子的无规运动可以用分数布朗函数作为概率密度函数来描写,而其分数维数可根据分子运动的图像确定。本文以流体Ar作为对象进行了分子动力学模拟,根据分子动力学模拟的结果提取了运动的分数维数,构造了描述分子无规运动的分数布朗函数,并对所提出的假说进行了验证。     

金刚石/硅(001)异质界面的分子动力学模拟研究

采用分子动力学方法模拟研究了未重构的金刚石/硅(001)面相接触时界面层原子的弛豫过程及所形成的异质界面的结构特征.硅碳二元系统中原子间的相互作用采用Tersoff多体经验势描述.弛豫前沿[110]与[110]方向界面碳硅原子数之比均为3∶2.界面碳硅原子总数之比为9∶4.弛豫后金刚石与硅界面处晶格匹配方式改变为[110]方向基本上以3∶2关系对准,而[110]方向大致以1∶1关系对准.相应地,界面碳硅原子总数之比接近3∶2.界面下方部分第二层硅原子在弛豫过程中向上迁移至界面是引起这种变化的原因,同时该层其他原子及其底下一到两个原子层厚度的区域在[001]方向上出现一定程度的无序化转变倾向.金刚石/硅异质界面处的硅碳原子发生强烈键合,形成平均键长为0.189nm的硅碳键.研究证实,晶格匹配主要呈现界面及其附近硅原子迎合界面碳原子排列的特点. 关键词： 金刚石 硅 异质界面 分子动力学