vasp计算分子动力学，结构能量一直小范围升高怎么解决

helpcal

请教一下vasp算分子动力学(MD)为什么结构能量一直升高。文献里是震荡的。
这是按照和文献一样的参数设的。用的是Nose热浴，NVT系综。步长2fs，在400K算5000步（10ps）。测试的材料是单层In2Se3, 一个unit cell是5颗原子。文献是用的441cell和221k点，我测试过441cell的111和221k点以及551cell的111k点结果都一样。图里放的是和文献完全一样的输入（441cell和221k点）（由于集群时间限制还剩大概三四十步没跑完，所以最后一点点是空的，不过不影响，其他2个case图是一样的）

下面是相关参数：
#MD
POTIM = 2
TEBEG = 400
TEEND = 400
MDALGO = 2
SMASS  = 1.0


NSW    = 5000
IBRION =      0

ISIF = 2

温度震荡是正常的。但是能量显然有问题。一个是一直在增加，还有一个是完全没有文献中大致2eV的震荡，才增加0.2eV大概。有没有大佬知道这是什么问题？中英文资料和chatgpt都已经搜过但也没找到解决办法。

helpcal

查了下发现这个结果根本没到平衡，请问有人知道怎么做么？除了换热浴有可能，其他的是难道给一个 随机的初始结构么

helpcal

乐平 发表于 2023-6-26 11:46
你并不能完全保证输入文件一模一样，因为你的结构可能和文献有区别。比如你的初始结构也许优化得更精确 ...

好的，感谢。他那情况感觉是额外做了些什么没在文献里说出来。

乐平

helpcal 发表于 2023-6-26 08:20
右图是文献的结果。应当和其一样能量有一定范围的震动，但这用一样的输入文件得到完全不同的结果很奇怪。 ...

用一样的输入文件得到完全不同的结果很奇怪

你并不能完全保证输入文件一模一样，因为你的结构可能和文献有区别。比如你的初始结构也许优化得更精确一些？或者文献里还有某些细节没说？

另外，你的结果一点也不奇怪，反而比文献里的能量变化幅度更小，更稳定。

Weldingspock

helpcal 发表于 2023-6-26 08:20
右图是文献的结果。应当和其一样能量有一定范围的震动，但这用一样的输入文件得到完全不同的结果很奇怪。 ...

可以的，你再统计到每个原子的能量变化情况，1meV左右的话已经是非常稳定了

helpcal

乐平 发表于 2023-6-25 22:46
你可以看看左上图能量变化的纵坐标， -286.520 到 -286.505，变化范围才 0.15 eV 已经很小了。

你可以把 ...

右图是文献的结果。应当和其一样能量有一定范围的震动，但这用一样的输入文件得到完全不同的结果很奇怪。这是我非常疑惑的一点。

不过单看我自己的结果，这种情况是能判断结构是热力学稳定的吗，因为能量变化非常小？

乐平

你可以看看左上图能量变化的纵坐标， -286.520 到 -286.505，变化范围才 0.15 eV 已经很小了。

你可以把左上图能量变化的“纵坐标范围”和右上图左上图能量变化的“纵坐标范围”统一，你看看有多平，多直。