计算化学公社

标题: 请教一个ASE对接VASP进行结构优化的无法收敛的问题 [打印本页]

作者
Author: sun111 时间: 2023-7-2 09:43
标题: 请教一个ASE对接VASP进行结构优化的无法收敛的问题
在进行结构优化的时候一直无法收敛，体系是在氧化亚铜上吸附小分子，修改过ediff、kpts、maxmove、dt、dtmax参数都无法收敛，下面附上文件，希望大家能指导一下，谢谢 (, 下载次数 Times of downloads: 2) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0)

作者
Author: 含光君 时间: 2023-7-2 15:29
用BFGS优化器不要用FIRE

作者
Author: 乐平 时间: 2023-7-2 15:35
ASE 相当于是用来准备 VASP 的输入文件吗？然后再当做作业管理来提交计算？

作者
Author: 含光君 时间: 2023-7-2 16:56

乐平发表于 2023-7-2 15:35
ASE 相当于是用来准备 VASP 的输入文件吗？然后再当做作业管理来提交计算？

这应该是用ASE搭workflow算恒电势

作者
Author: 乐平 时间: 2023-7-2 19:10

含光君发表于 2023-7-2 16:56
这应该是用ASE搭workflow算恒电势

Workflow 就是设置参数，提交任务的流程吧？

相当于用 ASE 准备 VASP 输入文件 INCAR, 然后提交计算任务的工具？

作者
Author: 含光君 时间: 2023-7-2 20:07
本帖最后由含光君于 2023-7-2 20:10 编辑

乐平发表于 2023-7-2 19:10
Workflow 就是设置参数，提交任务的流程吧？

相当于用 ASE 准备 VASP 输入文件 INCAR, 然后提交计算任 ...

确实，他用的这个代码应该是参考段志遥老师的文章，https://pubs.acs.org/doi/full/10.1021/acs.jpcc.1c02998
通过不断修正体系的电荷量在恒电势条件下算opt
ASE的doc中做过benchmark，FIRE优化器收敛速度慢于其他优化器。https://wiki.fysik.dtu.dk/gpaw/d ... e/ase_optimize.html

作者
Author: 卡开发发 时间: 2023-7-3 14:06

含光君发表于 2023-7-2 20:07
确实，他用的这个代码应该是参考段志遥老师的文章，https://pubs.acs.org/doi/full/10.1021/acs.jpcc.1c0 ...

不太好说，其中BFGSLineSearch以及Scipy的CG这些算法需要做LineSearch，这些LineSearch的过程测试并未提到是否被记录，算上这些FIRE其实不见得太吃亏。虽然基于QN和CG的方案在严重偏离二次型函数的时候LineSearch确实可以提升鲁棒性。

作者
Author: 卡开发发 时间: 2023-7-3 14:12

含光君发表于 2023-7-2 15:29
用BFGS优化器不要用FIRE

LZ的问题应该可能先要好好测试下两个动能截断，ENCUT=400eV对波函数而言可能不一定是足够的，不只是这样prec其实也要做充分的设置，或是直接去设置ENAUG，这是因为引入隐式溶剂求PB方程对密度网格也会有要求。否则力噪声特别大的情况再好的算法也无济于事。

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