氮化碳周期性体系在高斯中优化问题求助

liusci · 发表于 Post on 2023-10-23 09:03:47

本帖最后由 liusci 于 2023-10-23 09:05 编辑

如图这个是C3N5的结构图，

我想算某个分子在上面的吸附能，由于我对CP2K使用还不熟，所以用高斯进行计算。

我搭建了模型进行优化，边缘C原子固定并加氢，结果优化完不是平面的结构，不知道哪里出了问题，请各位老师指点一下。

这个是我的优化条件：# opt freq scrf=(solvent=water) def2tzvp em=gd3bj pbe1pbe

谢谢各位老师！

Ahey · 发表于 Post on 2023-10-23 09:15:56

感觉还是这个簇模型太小，此外你这个优化用的基组太大了，可以放小点节约耗时

zjxitcc · 发表于 Post on 2023-10-23 13:46:18

这些帖子问题都是一样的，可归为一类。应立即改用擅长周期性计算的程序（例如QE, VASP, CP2K等）。我都快成复读机了:(
http://bbs.keinsci.com/thread-40548-1-1.html
http://bbs.keinsci.com/thread-40559-1-1.html

KSeGaSn · 发表于 Post on 2023-10-23 17:34:03

把平面的周期性结构改成用部分结构的团簇模型时优化后出现结构扭曲非常常见。另外你周期性结构的图片里N=N连接的长度很小甚至可能周期性结构里也不是完全平面的（N=N间垂直于纸面了）。还是好好用周期性结构的方法计算，否则这样算出来恐怕也没有任何价值（比如催化等方面应用的话对周期性结构的模拟是标配）

sobereva · 发表于 Post on 2023-10-24 04:25:56

鼓捣用Gaussian优化这种体系花的时间，用Multiwfn早就创建好输入文件用CP2K跑完了
使用Multiwfn非常便利地创建CP2K程序的输入文件
http://sobereva.com/587（http://bbs.keinsci.com/thread-21668-1-1.html）

随便给你举个例子，是北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）里我讲单层石墨炔变胞优化的创建输入文件的过程

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

[新手求助] 氮化碳周期性体系在高斯中优化问题求助