计算化学公社

标题: 求助:使用cp2k对γ-Al2O3结构优化后和文献优化结果不同 [打印本页]

作者
Author: 卡卡西永恒 时间: 2023-9-12 20:48
标题: 求助:使用cp2k对γ-Al2O3结构优化后和文献优化结果不同
各位老师好，我使用MS软件按照图1中的数据建立了γ-Al2O3的晶胞，接着在multiwfn里分别设置理论方法PBEsol,设置k点3,2,2，设置任务优化原子和坐标，其他保持默认不变，生成cp2k的输入文件进行计算，最后收敛，但是优化后的结果和别人文献中的结果有差别，想问各位老师这种情况正常吗，我优化的结果是否是正确的
文献1中优化结果为：a=5.59 Å, b=8.42 Å, c=8.09 Å ， β=90.58°
文献2中优化结果为：a=5.585 Å，b=8.406 Å，c=8.079 Å，β=90.555°
我的优化结果为：a=5.54622 Å，b=8.36104 Å，c=8.07013 Å，β=90.514°

作者
Author: sobereva 时间: 2023-9-12 21:55
基本没相差多少
差的那一点是基组、CUTOFF、收敛限之类因素

作者
Author: 卡卡西永恒 时间: 2023-9-12 23:27

sobereva 发表于 2023-9-12 21:55
基本没相差多少
差的那一点是基组、CUTOFF、收敛限之类因素

谢谢sob老师，那我这个优化是不是也没有问题，可以继续进行后续切面等操作

作者
Author: sobereva 时间: 2023-9-12 23:35

卡卡西永恒发表于 2023-9-12 23:27
谢谢sob老师，那我这个优化是不是也没有问题，可以继续进行后续切面等操作

没有明显问题
如果你较真的话，也可以把基组提升一点（TZ档次）、CUTOFF提升一点、收敛限设严一点，看看是否相符程度更高

作者
Author: 卡卡西永恒 时间: 2023-9-13 18:06

sobereva 发表于 2023-9-12 23:35
没有明显问题
如果你较真的话，也可以把基组提升一点（TZ档次）、CUTOFF提升一点、收敛限设严一点，看看 ...

sob老师，我把原来的基组改为了TZVP-MOLOPT-GTH，运算的时候就报错了，这是什么原因呢

作者
Author: sobereva 时间: 2023-9-13 21:56

卡卡西永恒发表于 2023-9-13 18:06
sob老师，我把原来的基组改为了TZVP-MOLOPT-GTH，运算的时候就报错了，这是什么原因呢

提示写明了，此基组在基组文件里对此元素没有定义。
BASIS_MOLOPT_UZH里的基组种类丰富得多

作者
Author: 卡卡西永恒 时间: 2023-9-14 16:05

sobereva 发表于 2023-9-13 21:56
提示写明了，此基组在基组文件里对此元素没有定义。
BASIS_MOLOPT_UZH里的基组种类丰富得多

您是指下面图中的这些基组吗，我看Multiwfn里好像没有这些基组，是要手动修改输入文件吗，请问sob老师有没有关于修改基组的推荐资料呀，对这还不太了解，我手动改了一下报错了

作者
Author: sobereva 时间: 2023-9-14 19:00

卡卡西永恒发表于 2023-9-14 16:05
您是指下面图中的这些基组吗，我看Multiwfn里好像没有这些基组，是要手动修改输入文件吗，请问sob老师有 ...

自己改输入文件里的基组/赝势文件名、基组/赝势名就完了

这都是CP2K最基本的使用常识知识，在北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）里面有非常全面详细的讲解

作者
Author: 卡卡西永恒 时间: 2023-9-15 10:10

sobereva 发表于 2023-9-14 19:00
自己改输入文件里的基组/赝势文件名、基组/赝势名就完了

这都是CP2K最基本的使用常识知识，在北京科音 ...

好的，谢谢sob老师

作者
Author: 辣辣辣 时间: 2025-7-6 14:26
请问这是参考的那一篇文献呀

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