计算化学公社

标题: 求助，间断进行的吸附计算与持续进行的计算结果是否一样？ [打印本页]

作者
Author: gch123 时间: 2019-11-18 16:58
标题: 求助，间断进行的吸附计算与持续进行的计算结果是否一样？
请教一下各位老师，在进行吸附过程计算时遇到突发状况导致计算停止（如停电等），把计算停止时的优化结构作为初始结果继续进行计算，这样的计算结果与一次性计算完成的计算结果会有差别吗？

作者
Author: liyuanhe211 时间: 2019-11-18 17:12
MS不清楚，对分子体系，高斯的优化路径会有一定差别，但通常会收敛至同一个极小点，结果通常没有区别（但也有例外）。路径不同的具体原因我不特别清楚，猜测或许与优化算法会利用之前n步的位置和梯度有关，或许还有数值问题（另外如果提供精确的二阶导数，也只有计算的第一步有二阶导数）。

作者
Author: gch123 时间: 2019-11-18 20:06

liyuanhe211 发表于 2019-11-18 17:12
MS不清楚，对分子体系，高斯的优化路径会有一定差别，但通常会收敛至同一个极小点，结果通常没有区别（但也 ...

意思是有很小的概率会出现不一样的结果吗？

作者
Author: 卡开发发 时间: 2019-11-18 21:00

gch123 发表于 2019-11-18 20:06
意思是有很小的概率会出现不一样的结果吗？

最终构型一般无差别（无差别指的是落在数值误差之内），即便最终不太一样也能够通过调整初始结构来避免。除非几种稳定结构的构型特别相近，这种情况其实一般的构型优化也就不是特别好解决了。

同样使用QN-BFGS和GDIIS的情况，这种误差来源应该主要是：
1、Hessian初始获取的方式，如果完全重新算，Hessian一般是做n步速降使得梯度足够小然后再用BFGS更新或者是采用价力场进行初猜等。
2、GDIIS过程中，DIIS部分利用了前m步的梯度和Hessian信息来获取步长，重新算的情况GDIIS也是重新更新。

作者
Author: gch123 时间: 2019-11-18 22:21

卡开发发发表于 2019-11-18 21:00
最终构型一般无差别（无差别指的是落在数值误差之内），即便最终不太一样也能够通过调整初始结构来避免。 ...

谢谢老师

有点明白了

作者
Author: gch123 时间: 2019-11-19 08:46

卡开发发发表于 2019-11-18 21:00
最终构型一般无差别（无差别指的是落在数值误差之内），即便最终不太一样也能够通过调整初始结构来避免。 ...

老师，请问一下在结构优化过程中，能不能先使用DMol3进行优化得到相对稳定的结构之后再使用castep进行优化？这样应该能减少优化的时间吧

作者
Author: 卡开发发 时间: 2019-11-19 12:31

gch123 发表于 2019-11-19 08:46
老师，请问一下在结构优化过程中，能不能先使用DMol3进行优化得到相对稳定的结构之后再使用castep进行优 ...

不一定。也可能因为赝势或基组差异，导致结构朝着你初始的结构相反的方向跑。

欢迎光临计算化学公社 (http://ccc.keinsci.com/)