CP2K在计算石墨烯结构的单点能的报错求助

learn-GAUSSIAN · 发表于 Post on 2023-5-16 21:13:25

本帖最后由 learn-GAUSSIAN 于 2023-5-16 21:13 编辑

本人是CP2K新手，想计算一个质子+石墨烯结构的单点能，建模后使用Multiwfn非常方便地构建了输入文件结构，构建了20*20*60的盒子，并且设置了XY周期性，采用B3LYP方法进行计算，质子位于石墨烯上方很远处(25埃)。
现在输出文件有三个问题：
一是：WARNING:: No connections detected for Hydrogen - Atom Nr: 153 !
应该说的是质子没有与之相连的部分
二是：*** WARNING in dbcsr_mm.F:294 :: Using a non-square number of MPI ranks ***
*** might lead to poor performance. Used ranks: 23 Suggested: 25 49 ***
应该说的是没有使用并行
三是： *** WARNING in qs_scf.F:598 :: SCF run NOT converged ***
这是不收敛的问题，但是我觉得设的参数应该问题不大，我把文件已上传，还请各位不吝赐教！

wzkchem5 · 发表于 Post on 2023-5-16 21:31:05

先确认几个问题：
（1）你确定你要算裸露的质子？注意但凡体系里有溶剂，质子必然不是裸露的，例如水溶液里质子应该按H3O+甚至H9O4+来算
（2）如果你确定你的质子是裸露的，那么你确定石墨烯不会把质子还原成氢原子？裸露的质子的氧化性是很强的，石墨烯的pi电子又比较容易给出来，因此石墨烯到质子可能会发生单电子转移

sobereva · 发表于 Post on 2023-5-16 21:34:30

Z方向尺寸太小

仔细看
解决CP2K中SCF不收敛的方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）
用SCREEN_ON_INITIAL_P T直接让程序产生初猜可还行

learn-GAUSSIAN · 发表于 Post on 2023-5-16 21:34:43

wzkchem5 发表于 2023-5-16 21:31
先确认几个问题：
（1）你确定你要算裸露的质子？注意但凡体系里有溶剂，质子必然不是裸露的，例如水溶液 ...

你好老师，非常感谢您的回复。我确定是质子。我现在只是想测试算一下体系中当质子位于石墨烯表面5埃时的能量，再算一下质子位于远处的能量，最后再将质子放在盒子中做一个能量。最后我想作个比较。在目前主要的问题在于我这个体系不收敛啊，我也不知道该怎么解决

wzkchem5 · 发表于 Post on 2023-5-16 23:44:04

learn-GAUSSIAN 发表于 2023-5-16 14:34
你好老师，非常感谢您的回复。我确定是质子。我现在只是想测试算一下体系中当质子位于石墨烯表面5埃时的 ...

那么石墨烯往质子的电子转移是热力学有利的吗？如果是，那么你要算的体系根本不是基态，而是激发态，用普通的SCF是算不了的

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

[CP2K] CP2K在计算石墨烯结构的单点能的报错求助