计算化学公社

标题: 怎样设置cp2k输入文件优化得到D6h对称性的石墨烯 [打印本页]

作者
Author: mmxll 时间: 2024-9-15 13:28
标题: 怎样设置cp2k输入文件优化得到D6h对称性的石墨烯
老师们好，想要用cp2k优化得出D6h对称性的石墨烯，但是优化得出的是D2h，想要问一下如何在输入文件中设置才可以。下面是我的输入文件

作者
Author: sobereva 时间: 2024-9-16 00:49
没法设置
想要保持对称性尽可能高，应当用Gaussian等支持利用点群、原子中心积分格点的程序算。而CP2K基于立方数值格点计算，数值计算误差注定导致体系的对称性没法非常理想，而且优化还是基于笛卡尔坐标的，对这种孤立体系优化收敛也慢。

而且还有各种硬伤。哪能用DZVP这么烂的基组，起码也该用DZVP-MOLOPT-SR-GTH。貌似你要算弱相互作用问题，若有特殊原因非要用全电子基组也应该pob-TZVP-rev2。对碳原子考虑相对论效应作甚，&RELATIVISTIC字段毫无意义，而且你这基组根本不是给相对论计算设计的。当前盒子X、Y方向尺寸对于MT当非周期性计算根本不够，违背MT的基本要求，而Z方向又严重浪费。
缺乏CP2K和理论常识就尽量用Multiwfn直接产生输入文件，减少出现硬伤的可能，而且不懂的不要乱改
使用Multiwfn非常便利地创建CP2K程序的输入文件
http://sobereva.com/587（http://bbs.keinsci.com/thread-21668-1-1.html）
并且强烈建议通过北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）系统性学一遍基本知识免得像这样瞎算。

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