面心立方Ce同构相变的分子动力学模拟

摘要

基于嵌入原子法,本文给出了一个金属Ce原子间的相互作用势.利用该势分别计算了γ-Ce和α-Ce的晶格常数、结合能、弹性常数,计算结果与实验或第一性原理研究中得出的数值符合得较好.给出了两相Ce中如点缺陷形成能、表面能、层错能以及孪晶能等晶体缺陷形成能.通过分析两相Ce的声子谱,得出了不同温度下两相的晶格振动熵差,其中在室温条件下约为0.67kB/atom.还利用分子动力学模拟得出了该相变的等温线,并且利用径向分布函数分析了相变前后两相的晶体结构,确认了该相变为面心立方同构相变,即Ce的α-γ相变.由此表明,本文的嵌入原子法势,不仅可以分别合理地描述γ-Ce和α-Ce,还可以反映γ-Ce和α-Ce两相之间的相变.