请教RNEMD法求石墨烯纳米带热导率的问题。

vayne

在计算过程中，我得到的一维石墨烯纳米带温度分布如图1和图2所示。
出现以下两个问题：
1.明明温度设定的是400K，但温度分布峰值总高于40K
2.两边斜率不对称。我以为是弛豫步数的问题，但增加了弛豫步数，温度反而更高，斜率依然不对称。
请问有遇到过相似问题的前辈吗，求告知如何解决。
我的in文件在附件中给出，采用的边界条件是 s p p。

vayne

啊 问题出在哪里。