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

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

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

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

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

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

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

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

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

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

你可以把 ...

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

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

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

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