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

查看数: 985 | 评论数: 7 | 收藏 Add to favorites 1

关灯 | 提示：支持键盘翻页<-左右->

发布时间: 2023-7-2 09:43

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

回复 Reply

卡开发发 发表于 Post on 2023-7-3 14:12:57

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

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

参与人数 Participants 2	eV +7	收起理由 Reason
sun111	+ 4	谢谢分享
含光君	+ 3	精品内容

卡开发发 发表于 Post on 2023-7-3 14:06:15

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

含光君 发表于 Post on 2023-7-2 20:07: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

参与人数 Participants 1	eV +5	收起理由 Reason
乐平	+ 5	谢谢分享

乐平 发表于 Post on 2023-7-2 19:10:59

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

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

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

含光君 发表于 Post on 2023-7-2 16:56:05

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

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

乐平 发表于 Post on 2023-7-2 15:35:26

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

含光君 发表于 Post on 2023-7-2 15:29:33

用BFGS优化器不要用FIRE

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register