如何使用CP2K中的SIRIUS模块计算铁胞，并输出molden文件

牧生

近几日看了http://bbs.keinsci.com/thread-50094-1-1.html，社长提到了“某些磁性的d族金属表面体系SCF难收敛，典型例子就是诸如Ni(111)表面、Fe(001)表面这种，而且改各种SCF收敛相关设置也不太好解决。碰到这种情况可以用CP2K的SIRIUS模块基于平面波计算”。于是我开始了相关的学习，重新编译了一遍CP2K-2024.3，并安装了SIRIUS，也测试tests文件夹中的Fe.inp，飞快且顺利的跑完了，产生一个out文件，但我需要产生molden文件，以期望后续使用Multiwfn做分析。

第一个操作的问题：
在这里 http://sobereva.com/651 提到

1. 在CP2K输入文件的&DFT字段中加入以下内容就可以在SCF收敛后在当前目录下产生.molden为后缀的文件
   
2. 
   
3. &PRINT
   
4.   &MO_MOLDEN
   
5.     NDIGITS 9
   
6.   &END MO_MOLDEN
   
7. &END PRINT

复制代码

于是我在Fe.inp的文本中加入相应的字段，并且注释NGRIDK这一行，变成
……
  &DFT
    &XC
      &XC_FUNCTIONAL
        &LDA_C_PZ
        &END LDA_C_PZ
        &LDA_X
        &END LDA_X
      &END XC_FUNCTIONAL
    &END XC
  &PRINT
    &MO_MOLDEN
      NDIGITS 9
    &END MO_MOLDEN
  &END PRINT
  &END DFT
……


计算完成后，也只得到了Fe.out，并没有得到molden。我猜是因为换用了SIRIUS，就不存在“在SCF收敛后在当前目录下产生.molden为后缀的文件”这种情况了。那么在用了SIRIUS这种情况下，我该如何得到molden文件呢？


第二个操作的问题：
tests文件夹中的Fe.inp，应该就只是计算了单点。如果我要用结构优化，且就只想用CP2K，不用QE或者vasp，如http://bbs.keinsci.com/thread-41115-1-1.html中描述的情况，那么在inp文件中，我该选择何种理论方法，该在哪个位置加什么关键词呢？


第三个操作的问题：
我把tests文件夹中的Fe.inp载入Multiwfn，扩胞为2X2X2，然后用CP2K去跑这个超胞。把原来Fe.inp中cell的大小换成这个超胞cell信息，用的是      
A [bohr] 20.94416766     0.00000000     0.00000000     
B [bohr]  0.00000000    20.94416766     0.00000000
C [bohr]  0.00000000     0.00000000     20.94416766
把超胞的每个原子替换原来的那两个原子。然后倒是可以正常跑起来，见附件中的Fe.out，很多这种提示，有什么问题吗
Warning in function "get_rho_up_dn" at /home/jing/software/cp2k-2024.3/tools/toolchain/build/SIRIUS-7.5.2/src/density/density.hpp:854
Interstitial charge density has negative values
most negatve value : -0.0066859



Fe.inp (2.4 KB, 下载次数 Times of downloads: 20)

2024-12-3 12:22 上传 Uploaded

点击下载
Click to download

Fe.out (45.59 KB, 下载次数 Times of downloads: 4)

2024-12-3 12:22 上传 Uploaded

点击下载
Click to download

顺带在此实名赞扬一下google，我发帖后一分钟，在google搜“CP2K+MOLDEN+SIRIUS”，都能搜到本贴了。

qczcb2

molden文件是包含GTF和基函数信息的，sirius只有平面波，是不是原理上就不可行呢?收敛后应该已经有准确的电荷密度了，不知道能不能用于默认方法初猜

牧生

qczcb2 发表于 2024-12-3 19:02
molden文件是包含GTF和基函数信息的，sirius只有平面波，是不是原理上就不可行呢?收敛后应该已经有准确的电 ...

那这样的话，总得输出个什么文件用于可视化分析啊。最好是Multiwfn可以读取的

sobereva

molden文件本来就不能记录平面波基函数描述的轨道，而且Multiwfn目前也不支持平面波基函数

loneliness

请问cp2k的SIRIUS模块比起QE的优缺点是啥呢

牧生

loneliness 发表于 2025-3-2 12:06
请问cp2k的SIRIUS模块比起QE的优缺点是啥呢

我试过在CP2K中使用SIRIUS模块来计算我这个体系，但速度很慢。最终我放弃了

在http://bbs.keinsci.com/thread-50094-1-1.html中

1. 值得一提的是，CP2K里有个专门做纯平面波计算的SIRIUS模块，可以得到和Quantum ESPRESSO极为相近的结果。如果你就是有特殊原因非要做平面波计算也可以用这个模块。

复制代码

我不会QE，而且没有Multiwfn这样好用的傻瓜工具，我也没打算去学QE，无法解答。