计算化学公社

标题: CP2K在计算石墨烯结构的单点能的报错求助 [打印本页]

作者
Author:
learn-GAUSSIAN    时间: 2023-5-16 21:13
标题: CP2K在计算石墨烯结构的单点能的报错求助
本帖最后由 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 ***
这是不收敛的问题,但是我觉得设的参数应该问题不大,我把文件已上传,还请各位不吝赐教!

作者
Author:
wzkchem5    时间: 2023-5-16 21:31
先确认几个问题:
(1)你确定你要算裸露的质子?注意但凡体系里有溶剂,质子必然不是裸露的,例如水溶液里质子应该按H3O+甚至H9O4+来算
(2)如果你确定你的质子是裸露的,那么你确定石墨烯不会把质子还原成氢原子?裸露的质子的氧化性是很强的,石墨烯的pi电子又比较容易给出来,因此石墨烯到质子可能会发生单电子转移
作者
Author:
sobereva    时间: 2023-5-16 21:34
Z方向尺寸太小

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

作者
Author:
learn-GAUSSIAN    时间: 2023-5-16 21:34
wzkchem5 发表于 2023-5-16 21:31
先确认几个问题:
(1)你确定你要算裸露的质子?注意但凡体系里有溶剂,质子必然不是裸露的,例如水溶液 ...

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

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




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3