计算化学公社

标题: 有机分子晶体结构优化请教 [打印本页]

作者
Author: xjtuglr 时间: 2024-7-14 15:43
标题: 有机分子晶体结构优化请教
各位老师，我仿照卢老师上课计算无水草酸晶体太赫兹光谱例子对另外一个分子晶体进行计算，只含有C H O三种元素，优化结构和晶胞，在PBE-D3（BJ）/pob-DZVP-rev2级别下计算，计算的cutoff目前设置的400，开启了OT，为了避免有虚频率在严格的阈值基础上将EPS-scf设小到1e-8，maxforce设小到1e-5，计算了第一轮优化以后都没有完全收敛，经过400步优化以后max gradient结果大了一个数量级没有收敛，其他指标均达到，第二次续算400步以后还是没有收敛，并且5个参数最后一步都没有显示YES，尝试缩小了置信半径还是没有成功收敛，请问这种情况有什么好的解决办法吗？

同时一开始结构优化的时候尝试使用过DZVP-MOLOPT-SR-GTH结合PBEsol，采用不严格的阈值计算过，成功收敛，针对这种分子晶体计算光谱的时候一般在是用赝势基组较好还是全电子基组较好？
分子晶体的cif文件是从https://www.crystallography.net/cod/2000967.html 这里下载的
谢谢各位！

作者
Author: sobereva 时间: 2024-7-15 02:11
没输入输出文件没法准确判断
可以考虑改用CG优化算法
只要没有GTH赝势下对CUTOFF要求极高的元素，全电子基组、GTH赝势基组、NLCC赝势基组都可以用。在确保不出虚频的条件下，全电子基组和NLCC赝势基组对CUTOFF的要求更低。

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