最近,尝试计算twin-graphene的热导率,为了简化,直接先使用石墨烯的tersoff势能进行EMD下热导率的计算,测试发现NVT和NPT下势能表现正常,但是从系统NVT或NPT切换到NVE开始热导率计算时,观测到温度发散,最终材料解体。
针对这个问题,修改了很多计算参数。最终发现原因如下:
猜测可能来源于势能本身的特征,应用于twin graphene这种内部有局域双层的结构时,可能存在部分原子距离急剧变小的情况,势能没有更好的描述,导致原子在分动模拟中,如果步长太大,会有反常的持续加速的结果(类似于原子靠的很近的势能表现),表现就是温度持续上升,除非有温度控制的系综,解决方法有以下:
1.就是使用很小的步长来准确描述势能急剧变化情况下原子的运动,经测试0.1fs步长,温度能够稳定。但是EMD用这个步长计算太耗时间。
2.当然可以选择专门针对相应体系的其他更好的势能模型,不过没有合适势能的情况下,经测试发现可以加入LJ势能,来避免原子太近或势能不合理的地方。但是tersoff和lj无法直接hybrid,一对原子要么采用tersoff,要么采用LJ描述,只能使用一种,因此,需要对晶胞内的原子类型做区分,并修改原来的tersoff势能,将局域双层结构中的上下层原子区分,关掉相应的tersoff部分,用LJ替代,经过测试发现修正后的tersoff势能能够很好的描述twin-graphene这种结构。
3.使用fix NVE/limit也能够回复正常,这个是系综直接控制原子每步偏离不能超过某个临界值,只是对热输运计算的影响需要进一步思考。