请教一个ASE对接VASP进行结构优化的无法收敛的问题

sun111

在进行结构优化的时候一直无法收敛，体系是在氧化亚铜上吸附小分子，修改过ediff、kpts、maxmove、dt、dtmax参数都无法收敛，下面附上文件，希望大家能指导一下，谢谢

含光君

用BFGS优化器不要用FIRE

乐平

ASE 相当于是用来准备 VASP 的输入文件吗？ 然后再当做作业管理来提交计算？

含光君

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

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

乐平

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

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

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

含光君

乐平 发表于 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

卡开发发

含光君 发表于 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确实可以提升鲁棒性。

卡开发发

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

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