计算化学公社

标题: cp2k做AIMD如何确保其模拟的结构能量在全局最小值而不是局部极小值呢？ [打印本页]

作者
Author: nusiew 时间: 2021-10-27 10:02
标题: cp2k做AIMD如何确保其模拟的结构能量在全局最小值而不是局部极小值呢？
本帖最后由 nusiew 于 2021-10-27 10:06 编辑

大家好，在用cp2k做AIMD的过程中，我发现建立的初始结构对最终输出的平衡结构有很大的影响，尤其是对于两相体系，对于均相体系，差异还不是特别明显。对于一个体系来说，其能量的局部极小点可能有很多个，但全局最小的点只有一个，我很好奇cp2k是如何能确定其模拟结果最终能落在全局最小点而不是局部极小点？相同原子组成，为什么初始结构会对最终模拟结果产生影响呢？

(, 下载次数 Times of downloads: 28)

作者
Author: sobereva 时间: 2021-10-27 12:32
你既然是做动力学，就跟极小点毫无关系，本来动力学也不可能恰好跑到极小点位置，顶多说是在什么或哪些极小点附近（或势阱里）运动。

根据遍历性假设，如果跑无穷长时间，势能面所有点都会被采样。而对于有限时间的模拟，显然得到的有限长度的轨迹会采样哪些部分和初始位置有直接关系。跑的时间越长、温度越高、势垒越低，得到的轨迹越可能采样势能面越广阔的区域。

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