非平衡分子动力学法模拟计算SiC材料的热导率

为了探究更高效的碳化硅(SiC)材料热导率的模拟方法,应用逆非平衡分子动力学(rNEMD)法及传统非平衡分子动力学(NEMD)法对β晶型SiC(β-SiC)材料的热导率进行模拟计算和对比; 2种方法的模拟过程均先建立横截面尺度小而轴向尺度大的棒状模型,采用周期性边界条件、应用修正嵌入原子法(MEAM)势函数,先后进行正则系综(NVT)的弛豫和微正则系综(NVE)内的动态沿轴向生成温度梯度的过程,分别利用傅里叶定律模拟计算得到SiC材料的热导率。结果表明:2种方法的计算结果均出现热导率随生成温度梯度的材料轴向尺度增加而增大的有限尺度效应,应用倒数拟合的外推法可以计算模拟体系沿轴向为无穷大时的宏观体相β-SiC材料的热导率; r NEMD法具有较高的计算效率,更适合热导率的模拟计算。