MP法计算石墨烯纳米带热导率所遇见的几个问题

vayne

我最近在用MP法计算石墨烯纳米带热导率的过程中，遇见了一些问题，自己尝试了一些方法然而并没能很好地解决。因此，想在此请教论坛里的各位前辈，希望老师们能不吝赐教。
Question1    计算得到的温度分布曲线两边温度响应明显不一样，如图1（a）和图1（b）。
Question2    在某些算例中，输出的温度数据tmp.file里的冷源附近数层 温度为0，如图2。
Question3    在某些算例中，设定的温度为600K，但输出的温度数据却一直在升高，如图3。

我第一步是令模型在NVT系综下驰豫至温度稳定，命令在这 “fix       NVT all nvt temp $T $T 0.1 drag 0.05”
我一开始认为是NVT下驰豫时间不够，从原来的0.25ns换成2ns，但是输出的温度反而更紊乱；
也试过更改 交换原子的频率（步数）；更改 模型分的层数；更改边界条件；等等
都没有解决上述三个问题！！！

请问有老师遇见过这些情况吗？希望能指点一下我。谢谢！

zhaoyp123

请问解决了吗