液态Au和Ag冷却过程的分子动力学模拟

采用反映原子多体相互作用的FS势模型,通过分子动力学方法研究了贵金属Au、Ag在温度2000~300 K之间的冷却过程,考察其结构变化特点,给出了冷却过程微观结构转变的重要信息.     

耐蚀合金Au3Cu高温冷却过程中能量及结构转变的分子动力学模拟

用分子动学模拟方法对液态Au3Cu冷却过程进行了研究,考察了不同冷却速度下Au3Cu结构变化特点,原子间相互作用势采用F-S多体势,结构分析采用键取向序和对分析技术。计算结果表明,冷却速度对液态Au3Cu能量及结构转变有重要影响,给出了不同冷却速度下液态Au3Cu结构转变的微观信息。     

金属Cu熔化结晶过程的分子动力学模拟

采用常温、常压分子动力学模拟技术,研究了在周期性边界条件下,由８６４个Ｃｕ原子构成的模型系统的熔化、结晶过程。原子间相互作用势采用ＥＡＭ势。模拟结果表明：在连续升温过程中,金属Ｃｕ在１５２０Ｋ熔化;以不同的冷速进行冷却,在较慢冷却条件下,液态Ｃｕ在１０１０Ｋ结晶;当冷速较快时,液态Ｃｕ形成非晶态。分析了升降温过程中熔体偶分布函数、原子体积、能量、ＭＳＤ随温度的变化特征。     

Cu3Au熔体快速冷却过程微观结构演化的分子动力学模拟

利用分子动力学模拟技术研究了金属间化合物Cu3Au熔体的双体分布函数g(r)在快速凝固条件下随温度的变化情况，结果表明，Cu3Au降温至700K时第二峰已发生劈裂，液态金属中已经产生了非晶态；用键对分析技术详细考察了Cu3Au中微观组团随温度的演化特点，液体中的键对数及多面体数与温度的关系都表明，Cu3Au在向非晶转变的过程中，的确发生了微观结构组态的变化，其中以液体中的缺陷多面体随温度变化最为剧烈。非晶不是过冷液态的“冻结”。     

Cu3Au熔体快速冷却过程微观结构演化的分子动力学模拟

利用分子动力学模拟技术研究了金属间化合物Cu3Au熔体的双体分布函数g(r)在快速凝固条件下随温度的变化情况,结果表明,Cu3Au降温至700K时第二峰已发生劈裂,液态金属中已经产生了非晶态;用键对分析技术详细考察了Cu3Au中微观组团随温度的演化特点,液体中的键对数及多面体数与温度的关系都表明,Cu3Au在向非晶转变的过程中,的确发生了微观结构组态的变化,其中以液体中的缺陷多面体随温度变化最为剧烈.非晶不是过冷液态的"冻结".     

低密度Ar熔化及结晶的分子动力学模拟

采用有位移力的ＬｅｎｎａｒｄＪｏｎｅｓ（１２－６）势对微正则系综下低密度Ａｒ系统（简约数密度为ρ＝０．８５）的一级相变过程进行了细致的分子动力学模拟,发现Ａｒ系统的熔化过程是原子的崩塌过程,结晶过程和理想的完整晶体不同,是一活化过程：形核长大过程随温度的降低进行,原来均匀分布在系统中的自由体积呈集中分布,由此系统达到了更稳定的结晶状态。     

熔体局部结构的分子动力学模拟研究

计算机技术的迅速发展使分子动力学模拟成为可能。这是一种直接方法用以更详细地了解液态的结构。文章介绍在分子动力学模拟产生的一系列平衡瞬态构型基础上的等近邻键序参数方法和使用这一方法对若干熔体局部结构的研究，包括若干碱金属卤化物熔盐，ＣａＦ２快离子态和熔融态，熔融态ＺｎＣｌ２和Ｒｂ２ＺｎＣｌ４，ＢＢＯ熔体和晶体生长母液等。     

液态氧化铝结构的分子动力学模拟

利用分子动力学模拟液态氧化铝从2663 K到过冷温度2223 K的结构变化,通过总双体分布函数得到的液态氧化铝的原子第一、第二近邻距离与Ansell研究小组的实验结果基本一致.通过配位数和偏角分布函数的分析,液态氧化铝的局域结构主要是由氧的三重铝原子配位和铝原子四重(AlO₄)和五重(AlO₅)氧原子配位组成.没有观察到可能因铝原子高、低配位数氧原子的团簇分离而出现的液-液相变.     

GaAs晶体在不同应变下生长过程的分子动力学模拟

GaAs晶体的高质量生长对于制造高性能高频微波电子器件和发光器件具有重要意义.本文通过分子动力学方法对GaAs晶体沿[110]晶向的诱导结晶进行模拟,并采用最大标准团簇分析、双体分布函数和可视化等方法研究应变对生长过程和缺陷形成的影响.结果表明,不同应变条件下GaAs晶体的结晶过程发生显著变化.在初始阶段,施加一定拉应变和较大的压应变后,体系的晶体生长速率发生降低,且应变越大,结晶速率越低.此外,随着晶体的生长,体系形成以{111}小平面为边界的锯齿形界面,生长平面与{111}小平面之间的夹角影响固液界面的形态,进而影响孪晶的形成.施加拉应变越大,此夹角越小,形成孪晶缺陷越多,结构越不规则.同时,体系中极大部分的位错与孪晶存在伴生关系,应变的施加可以抑制或促进位错的形核,合适的应变甚至可以使晶体无位错生长.本文从原子尺度上研究GaAs的微观结构演化,可为晶体生长理论提供理论指导.     

TiAl合金薄膜在冷却过程中结构变化的原子尺度计算研究

采用基于嵌入原子势模型的分子动力学方法模拟了熔融TiAl合金薄膜分别在急冷降温和连续降温两个冷却过程中微观结构的变化.通过应用原子平均能量、对分布函数和键对分析技术等对薄膜内原子局域结构进行的分析. 研究表明,在降温过程中,薄膜内局域结构出现分阶段变化的特点.在急冷降温过程中,薄膜结构的变化可以分为3个阶段,而在连续降温过程中,薄膜内局域结构的变化分为2个阶段. 关键词： 分子动力学 TiAl合金 薄膜 计算机模拟     

液态合金NiAl凝固过程中微观结构转变的分子动力学模拟

采用分子动力学模拟方法对液态NiAl凝固过程进行了研究 ,考察了不同冷却速度下液态NiAl结构变化特点 ,原子间相互作用势采用F S多体势 ,结构分析采用键取向序和对分析技术 .计算结果表明 ,冷却速度对液态NiAl结构转变有重要影响 ,在不同的冷却速度下 ,NiAl凝固过程出现了明显不同 ,冷速为 4× 10 13 和4× 10 12 K/s时 ,NiAl快速凝固为无序的非晶体结构 ;而在较慢的 8× 10 11K/s冷速下 ,NiAl凝固为晶态结构 .给出了不同冷却速度下液态NiAl结构转变的微观信息 .     

热镀铜层冷凝过程的分子动力学模拟

采用常温、常压分子动力学模拟方法和FS(Finnis Sinclair)势 ,研究了在周期性边界条件下由 5 0 0个原子构成的液态Cu模型系统的凝固过程 ,考察了不同降温速率下Cu的凝固行为 ,得到了不同温度、不同冷却速率下Cu的双体分布函数 ;采用HA键型指数法统计了各种小原子团在不同温度下所占比例 ,采用键取向序分析了体系降温全过程的局域取向对称性 ,得到原子体系微观结构组态变化的重要信息 ;最后 ,利用能量分析的方法对体系微观结构的变化进行了说明 ,给出了液态Cu冷凝过程中微观结构转变的重要信息 .     

Au纳米薄膜的熔化机理及其结构演变的分子动力学模拟

利用分子动力学模拟详细研究了不同厚度的Au纳米薄膜的熔化机理和结构演变. 模拟结果表明所有Au纳米薄膜的熔化行为分为两个阶段,即表面预熔和均相熔化. 只有最外层原子出现了预熔化行为, 其他内层原子在均相熔化之前始终保持稳定的固态,这与零维的Au纳米团簇和一维的Au纳米线的预熔化行为是不同的. 同时Au纳米薄膜的熔化温度随着薄膜厚度的增加而升高. 在预熔化过程中,在原子水平上发现了所有的Au纳米薄膜的f100g晶面向f111g晶面转变的表面重建过程. 对于最薄的L2纳米薄膜,当温度低于500 K 时表面应力不能诱导这样的表面重建. 然而一维的Au纳米线在更低温度下就能够观察到了由表面应力诱导的表面重建过程. 这主要是因为Au纳米线具有更高的比表面积所导致的. 另外研究结果还表明当模拟温度达到某一特定值时,由双原子层组成的Au纳米薄膜能够分裂成一维的纳米线.     

耐蚀合金Al3Ni和Ni3Al液态冷凝过程的比较研究

采用F S多体势对液态合金Al3Ni和Ni3Al在不同冷却速度下的微观结构及其转变机制进行了分子动力学模拟 ,得到了不同冷速下各温度的双体分布函数 ;采用HA键型指数法对其结构进行了分析 ,结果表明 :Al3Ni在两种冷速下均以非晶的形式出现 ,只是慢冷时体系的有序度略有升高 ;而Ni3Al的结构及能量转变受冷速影响较大 ,快冷时形成非晶 ,而慢冷时出现明显结晶 ;同样冷速下Al含量较少的Ni3Al体系的有序度高 ,更易形成晶体 ,晶体的形成过程中有能量突变 .     

量子分子动力学模拟液体钚的输运性质

采用第一原理分子动力学（QMD）方法模拟液体钚的输运性质.计算的粘性和扩散系数在较低温度时与文献有明显差异,在实验测量范围内,模拟结果与实验一致,温度升高时数值模拟结果趋于一致.利用QMD的模拟结果计算了应力自相关函数和速度自相关函数,结果表明：在温度较低时,液体钚呈现明显的强关联特性.对于具有强关联特性的液体,利用较短时间的QMD模拟结果,通过简单e指数拟合外推到t→∞得到的扩散系数和粘性具有较大偏差,这是造成本文模拟结果与文献结果出现差异的主要原因.通过增加QMD模拟时间步数,获得了更为准确的输运性质.     

CaF2熔化温度的分子动力学模拟

利用壳层模型分子动力学方法,研究了高温高压条件下CaF2的熔化温度,同时计算了温度为300K、压强上升到100GPa时CaF2 的状态方程.研究中考虑了分子动力学模拟的过热熔化,通过晶体的现代熔化理论,对CaF2 的分子动力学模拟熔化温度进行了修正, 获得了高温高压下CaF2的熔化温度.因此,常压下壳层模型分子动力学方法为研究物质熔化提供了一个很好的方法.     

V5+掺杂TiO2金红石结构的分子动力学

采用经典离子晶体相互作用势模型， 进行量子化修正，对TiO2金红石相替位式掺杂0.5mol％V5+进行分子动力学模拟，计算了常温常压下的键长、键角的分布函数，分析了晶体的微结构变化以及V5+对TiO6八面体的影响．模拟结果表明，掺杂V5+会导致Ti-O、O-O键长以及O-Ti-O键角的分布范围变宽，其中，键长的主峰位置未发生显著的移动， 而键角主峰向低角度方向移动． 由于V5+比Ti4+有较高的正电荷和较小的离子半径， 因此V5+向氧八面体间空隙扩散，且掺杂V5+迁移出原来的O6八面体，从而导致空隙邻近的TiO6八面体发生严重扭曲．在低掺杂情况下(0.5 mol％)虽然掺杂易造成结构畸变，但模拟体系整体仍维持金红石结构；所得模拟的键长变化、掺杂V5+迁移、晶相等结构特性与FTIR、Raman、XRD、ESR等实验结果相一致．     

基于接触数引导的迭代多组独立分子动力学模拟实现配体-蛋白质结合与解离过程

配体的结合与解离过程在蛋白质实现其生物学功能方面非常关键,因此对这些高度动态过程的研究变得非常重要. 尽管已有实验方法可以确定蛋白质-配体复合物的三维结构,但一般仅可获得静态图片. 随着计算机算力的快速提高以及算法的优化,分子动力学模拟在探索配体的结合与解离过程方面具有诸多优势. 然而,当系统变得足够大时,分子动力学模拟的时间和空间尺度成为了巨大的挑战. 本工作提出了一种研究配体-蛋白质结合与解离的增强采样工具,它基于配体和蛋白质之间形成的接触数来引导迭代多组独立分子动力学模拟. 在腺苷酸激酶的模拟结果中,观测到配体的结合和解离过程,而使用传统分子动力学模拟在同一时间尺度下则无法实现这一过程.