sobereva 发表于 2024-7-30 22:01
可靠性有可能遭审稿人质疑。可以拿一个体系结合0.5 fs跑跑和之前的对比。就算不报错，过大的步长也可能导 ...

谢谢老师耐心解答

gehan 发表于 2024-7-30 10:19
好的，谢谢老师的建议，我下来尝试一下。老师因为我已经算了很多这种在1100K下的带氢体系，timestep都用 ...

可靠性有可能遭审稿人质疑。可以拿一个体系结合0.5 fs跑跑和之前的对比。就算不报错，过大的步长也可能导致动力学行为明显不对。

你说的那种改法解决了问题大概率是巧合，不是原理上解决的办法。

sobereva 发表于 2024-7-30 06:14
1000多K跑带氢的AIMD用1fs不是很稳妥，特别是对于有可能成键断键的情况，尝试减小步长到0.5 fs或更小

好的，谢谢老师的建议，我下来尝试一下。老师因为我已经算了很多这种在1100K下的带氢体系，timestep都用默认的1fs，不过之前计算的大部分都没有什么明显的报错，那这样的结果还可靠吗？还有一个问题就是，我自己改了restart文件中的截断能(改大一些)和CSVR热浴中timecon(改小了)。这样，有些时候上述报错就会消失，可以正常续算。请问老师这样的做法可行吗，得到的结果暂未出现明显的不合理现象。

1000多K跑带氢的AIMD用1fs不是很稳妥，特别是对于有可能成键断键的情况，尝试减小步长到0.5 fs或更小