lammps案例：石墨烯热导率模拟计算（EMD方法)

大家好，我是小马老师。

在前面的文章中，介绍了非平衡态下石墨烯的热导率模拟方法，本文介绍第二种热导率模拟方法：使用平衡态分子动力学（EMD）计算热导率。

本文仍然以石墨烯热导率计算为例，以供大家对比参考。

在平衡态下计算热导率，主要计算公式为Green-Kubo。

用到的主要命令为compute heat/flux。

用法为：

compute myFlux all heat/flux myKE myPE myStress

其中，myKE为原子动能，myPE为原子势能，myStress为原子应力。

在使用compute heat/flux命令前，必须提前计算出这三个量的值。

下面给出石墨烯EMD热导模拟代码，代码已经注释。data文件可自己建模，也可加微信sunnyfirst888联系获取。

热导率具体计算方法在集训营会有详细介绍，如有需要可微信联系。

案例代码

本文作者小马老师正式推出一对一咨询辅导服务，根据课题方向不同详细讲解对应的in文件编写方法、模拟关键技术、数据后处理，经一对一咨询辅导后能够独立编写出in文件。

公众号:

320科技工作室