如何计算有机分子晶体内分子间的相互作用？

Rosefinch

请问一下大家，我想用CP2K来计算一个有机分子晶体内分子间的弱相互作用能大小，可以直接用BSSE来计算嘛？或者有什么更加高效的方法或软件？

sobereva

先搞清楚BSSE是什么意思再说，别乱用BSSE这个词，仔细看
谈谈BSSE校正与Gaussian对它的处理
http://sobereva.com/46

当你用BSSE问题很小的基组，比如TZV2PX-MOLOPT-GTH，算弱相互作用能时不特意用counterpoise或gCP考虑BSSE问题都是完全可以接受的。
倘若用诸如DZVP-MOLOPT-GTH，BSSE没小到可以忽略，此时可以通过CP2K直接支持的counterpoise任务在计算弱相互作用能时考虑对BSSE的修正，Multiwfn直接就能产生此任务的输入文件
使用Multiwfn非常便利地创建CP2K程序的输入文件
http://sobereva.com/587（http://bbs.keinsci.com/thread-21668-1-1.html）

几乎没有比CP2K更高效的。

Rosefinch

sobereva 发表于 2023-10-18 03:24
先搞清楚BSSE是什么意思再说，别乱用BSSE这个词，仔细看
谈谈BSSE校正与Gaussian对它的处理
http://sober ...

好的好的，谢谢sob老师，我想用DZVP-MOLOPT-GTH计算，请问我可以一次设置4个片段嘛？

Rosefinch

sobereva 发表于 2023-10-18 03:24
先搞清楚BSSE是什么意思再说，别乱用BSSE这个词，仔细看
谈谈BSSE校正与Gaussian对它的处理
http://sober ...

我在cp2k官网看到BSSE关键词下的GLB_CONF {integer} {integer} .. {integer}和SUB_CONF {integer} {integer} .. {integer}关键词有很多个{integer}，是不是说可以对多个片段间的相互作用能进行BSSE校正？然后我根据别的帖子写了一下包含四个分子片段的inp文件的BSSE部分，可以麻烦您帮忙看一下是否正确嘛？十分感谢！

1. &BSSE
   
2.     &FRAGMENT
   
3.       LIST       1..94
   
4.     &END FRAGMENT
   
5.     &FRAGMENT
   
6.       LIST      95..188
   
7.     &END FRAGMENT
   
8.     &FRAGMENT
   
9.       LIST      189..282
   
10.     &END FRAGMENT
    
11.     &FRAGMENT
    
12.       LIST      283..376
    
13.     &END FRAGMENT
    
14.     &CONFIGURATION # real(A)+real(B)+real(C)+real(D)
    
15.       GLB_CONF 1 1 1 1
    
16.       SUB_CONF 1 1 1 1
    
17.       CHARGE  0
    
18.       MULTIPLICITY   1
    
19.     &END CONFIGURATION
    
20.     &CONFIGURATION # real(A)
    
21.       GLB_CONF 1 0 0 0
    
22.       SUB_CONF 1 0 0 0
    
23.       CHARGE  0
    
24.       MULTIPLICITY   1
    
25.     &END CONFIGURATION
    
26.     &CONFIGURATION # real(B)
    
27.       GLB_CONF 0 1 0 0
    
28.       SUB_CONF 0 1 0 0
    
29.       CHARGE  0
    
30.       MULTIPLICITY   1
    
31.     &END CONFIGURATION
    
32.     &CONFIGURATION # real(C)
    
33.       GLB_CONF 0 0 1 0
    
34.       SUB_CONF 0 0 1 0
    
35.       CHARGE  0
    
36.       MULTIPLICITY   1
    
37.     &END CONFIGURATION
    
38.     &CONFIGURATION # real(D)
    
39.       GLB_CONF 0 0 0 1
    
40.       SUB_CONF 0 0 0 1
    
41.       CHARGE  0
    
42.       MULTIPLICITY   1
    
43.     &END CONFIGURATION
    
44.     &CONFIGURATION # real(A)+ghost(B)+ghost(C)+ghost(D)
    
45.       GLB_CONF 1 1 1 1
    
46.       SUB_CONF 1 0 0 0
    
47.       CHARGE  0
    
48.       MULTIPLICITY   1
    
49.     &END CONFIGURATION
    
50.     &CONFIGURATION # ghost(A)+real(B)+ghost(C)+ghost(D)
    
51.       GLB_CONF 1 1 1 1
    
52.       SUB_CONF 0 1 0 0
    
53.       CHARGE  0
    
54.       MULTIPLICITY   1
    
55.     &END CONFIGURATION
    
56.     &CONFIGURATION # ghost(A)+ghost(B)+real(C)+ghost(D)
    
57.       GLB_CONF 1 1 1 1
    
58.       SUB_CONF 0 0 1 0
    
59.       CHARGE  0
    
60.       MULTIPLICITY   1
    
61.     &END CONFIGURATION
    
62.     &CONFIGURATION # ghost(A)+ghost(B)+ghost(C)+real(D)
    
63.       GLB_CONF 1 1 1 1
    
64.       SUB_CONF 0 0 0 1
    
65.       CHARGE  0
    
66.       MULTIPLICITY   1
    
67.     &END CONFIGURATION
    
68.   &END BSSE

复制代码

Rosefinch

不好意思，我写的是不全的，但是我按照我写的这个inp投任务去计算，cp2k会自动完整正确地计算各个片段(1-body；2-body；3-body；4-body)的能量。
不过在计算我没有定义的片段时，out文件显示“CHARGE =     0           MULTIPLICITY =     0”，但是我看官网对于BSSE那一部分中MULTIPLICITY关键词里描述说“Default is 1 (singlet) for an even number and 2 (doublet) for an odd number of electrons.” 不知是不是out文件显示错误，好在看最后的结果似乎没有什么问题。