多形性相变的分子动力学模拟

本文利用分子动力学方法研究了KCl晶体在ρ＝ρ0时的温度相变。面心立方（fcc）和体心立方（bcc）两种结构的径向分布函数随温度的变化的情况说明，在高压下，发生着bcc结构相fcc结构的转变，bcc结构是不稳定的。为了选取合理的势参数，利用了分子动力学程序在T＝0时的性质，计算了NaCl和KCl晶体的零温状态方程，研究了它们在压力作用下发生的多形性相变。计算表明，NaCl和KCl晶体将分别于18.8GPa和5.9GPa发生从fcc到bcc的多形性相变，这些值相当接近实验结果。着说明本文选用的势参数势有一定精度的。