计算化学公社

标题: 如何计算有机分子晶体内分子间的相互作用？ [打印本页]
作者
Author: Rosefinch    时间: 2023-10-17 14:08
标题: 如何计算有机分子晶体内分子间的相互作用？
请问一下大家，我想用CP2K来计算一个有机分子晶体内分子间的弱相互作用能大小，可以直接用BSSE来计算嘛？或者有什么更加高效的方法或软件？
作者
Author: sobereva    时间: 2023-10-18 03:24
先搞清楚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更高效的。

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

好的好的，谢谢sob老师，我想用DZVP-MOLOPT-GTH计算，请问我可以一次设置4个片段嘛？
作者
Author: Rosefinch    时间: 2023-10-18 16:48
本帖最后由 Rosefinch 于 2023-10-18 17:40 编辑
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
复制代码


作者
Author: Rosefinch    时间: 2023-10-19 10:51
不好意思，我写的是不全的，但是我按照我写的这个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文件显示错误，好在看最后的结果似乎没有什么问题。



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