cp2k中关于MULTIPLICITY的设置疑问

yiranfengbai · 发表于 Post on 2024-5-21 01:11:39

本帖最后由 yiranfengbai 于 2024-5-21 03:20 编辑

老师您好，我最近一直在用CP2K优化CoFe2O4晶胞，使用PBE泛函，2X2X2超晶胞，使用了DFT+U，U值取自参考文献 Co: 2 eV，Fe: 3.3 eV。查文献发现CoFe2O4是反铁磁，磁矩：八面体Co：2.68，Fe: 3.96，四面体Fe: -4.12。输入文件inp中设置：八面体Co：3，Fe: 4，四面体Fe: -4，Multiwfn给出MULTIPLICITY=1+3*8=25，结果优化和单点能一直不收敛，参考Sob老师的博文修改了CG、IRAC、FULL_KINETIC，还是不收敛。Sob老师提示 MULTIPLICITY=25有问题，没考虑到2X2X2超胞。所以修改为MULTIPLICITY=1+3*8*8=193，发现单点能可以收敛，但几何优化不收敛，最后测试发现把FULL_KINETIC修改为FULL_ALL后很快就成功收敛了。然而在群中，ScE老师提示 ”MULTIPLICITY=193太离谱了”，这让我突然就很迷茫了。ULTIPLICITY=自旋向上电子 - 自旋向下电子+1，我的设置应该没错吧？查了好几天资料还是没弄明白，所以想问问老师：MULTIPLICITY=193到底是错误的吗？如果是错误的应该怎么设置呢？
ScE老师还提示: Co 和 Fe的U值应该设为5或6eV，PLUS_U_METHOD的选择：MULLIKEN或MULLIKEN_CHARGES有很大影响。我之前设置：Co:2 Fe: 3.3 PLUS_U_METHOD MULLIKEN时cp2k.out中一直WARNING：DFT+U值贡献为负。而改为PLUS_U_METHOD MULLIKEN_CHARGES(后面发现FULL_ALL时，PLUS_U_METHOD MULLIKEN不再报错)时就没有出现这样的WARNING。我猜测ScE老师的意思是应该设置：Co:5（或6） Fe:5（或6） PLUS_U_METHOD MULLIKEN。可是之前Sob老师提示CP2K中U值通常都较低设为2~4eV，有些文献中较大U值都是套用QE中的设置。请问老师，PLUS_U_METHOD对U的选择影响大吗？还有，我按照Co:2 Fe: 3.3 PLUS_U_METHOD MULLIKEN设置计算band gap，cp2k中显示1.61eV，Multiwfn: cp/7/2/1,3/0 显示-8.599（文献中band gap范围：0.7~1.2eV），请问老师这是为什么啊？
已经上传了inp和cp2k.out文件，molden文件太大了，在集群上下了1个多小时还没下载完。

sobereva · 发表于 Post on 2024-5-21 02:06:34

没具体文件和具体命令没法说Multiwfn为什么判断成那样

先压缩再下载啊
计算化学研究者应当了解的文件压缩的常识
http://sobereva.com/672（http://bbs.keinsci.com/thread-37881-1-1.html）

PLUS_U_METHOD对U会产生一定影响，但不至于能影响到好几eV。CP2K一般不会用4eV以上的Ueff。弄清楚Ueff值是哪来的。与其用不当的还不如不用，或者咬牙用杂化泛函

yiranfengbai · 发表于 Post on 2024-5-21 03:48:10

sobereva 发表于 2024-5-21 02:06
没具体文件和具体命令没法说Multiwfn为什么判断成那样

先压缩再下载啊

老师，我已经上传百度网盘了
链接: https://pan.baidu.com/s/1AWHXau21cOLObsZTB6gVeQ 提取码: ga4h

sobereva · 发表于 Post on 2024-5-21 04:51:57

没上传给Multiwfn用的得到band gap的输入文件

yiranfengbai · 发表于 Post on 2024-5-21 10:21:18

sobereva 发表于 2024-5-21 04:51
没上传给Multiwfn用的得到band gap的输入文件

老师，我重新上传了CoFe2O4.gif和计算band gap的输入文件。下面是我在Multiwfn中设置输出文件的操作，产生的inp中 MULTIPLICITY 25，我改为193，同时还修改了MULLIKEN_CHARGES、CG、IRAC。为了输出band gap信息，我在&PRINT中加入&MO_CUBES
启动Multiwfn后依次输入：
CoFe2O4.gif
cp2k

-1//单点能计算
4//使用OT方法
-9//
9//设置磁矩
33-40//四面体位Fe原子序号（通过VESTA判断），自旋向下
-4//
41，43，46，47，52，54，55，56//八面体位Fe，自旋向上
4//
42，44，45，48，49，50，51，53//八面体位Co，自旋向上
3//
q//
8//打开DFT+U
6//输出磁矩信息
10//输出水平
2//
3//设置超胞
2，2，2//
5//设置截断能
500，60//
0//
-2//产生molden文件
0//产生输入文件

sobereva · 发表于 Post on 2024-5-22 06:06:23

yiranfengbai 发表于 2024-5-21 10:21
老师，我重新上传了CoFe2O4.gif和计算band gap的输入文件。下面是我在Multiwfn中设置输出文件的操作，产 ...

我说的是Multiwfn里用cp/7/2/1,3/0的时候给Multiwfn用的输入文件

yiranfengbai · 发表于 Post on 2024-5-22 10:36:00

老师，我已经把使用的所以文件都打包放在3L的网盘中了，不太清楚您需要哪个。另外，老师我的MULTIPLICITY 设置是否有错误呢

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

[CP2K] cp2k中关于MULTIPLICITY的设置疑问