计算化学公社

标题: QCHEM用约束密度泛函CDFT遇到问题 [打印本页]
作者
Author: xyz0101    时间: 2023-6-27 21:36
标题: QCHEM用约束密度泛函CDFT遇到问题
在qchem中用CDFT来规定片段的电荷量，可以用一下命令来添加
$cdft
2
1 1 25
-1 26 38
$end
但这种方法只设置了不同片段的带电量，从而引起自旋多重度的变化该怎么设置呢？使用手册里没有很好说明。

作者
Author: myzpku    时间: 2023-6-28 14:05
看这页的例子5.24 https://manual.q-chem.com/latest/Ch5.S11.SS4.html
作者
Author: xyz0101    时间: 2023-6-28 14:17
我看到这个例子了
$cdft
2
1 1 3 s
-1 17 19 s
$end
我不清楚这个s设置的是什么含义
作者
Author: xyz0101    时间: 2023-6-28 14:19
myzpku 发表于 2023-6-28 14:05
看这页的例子5.24 https://manual.q-chem.com/latest/Ch5.S11.SS4.html

我更想了解的是他之后5.25这个例子是怎么来写cdft部分的？
作者
Author: myzpku    时间: 2023-6-28 17:10
xyz0101 发表于 2023-6-28 14:17
我看到这个例子了
$cdft
2

表示1-3号原子的spin (rho_alpha - rho_beta) 减去17-19号原子的spin等于2，即该三重态体系的自旋集中在第一个Cu center的情况。你可以跑一下这个sample然后看输出的“CDFT Becke Populations”进行验算：
  1.                   CDFT Becke Populations
  2. ----------------------------------------------------
  3.             Atom      Excess Electrons     Population (a.u.)    Net Spin
  4.        1      Cu        2.380846       31.380846          1.035867
  5.        2       O       -0.471315        7.528685          0.399140
  6.        3       O       -0.388478        7.611522          0.395532
  7.        4       C        0.033295        6.033295          0.034481
  8.        5       N       -0.439179        6.560821         -0.001311
  9.        6       N       -0.409184        6.590816         -0.016496
  10.        7       N       -0.348841        6.651159         -0.036626
  11.        8       H       -0.151261        0.848739         -0.001412
  12.        9       H       -0.152936        0.847064         -0.002978
  13.       10       H       -0.133281        0.866719         -0.003413
  14.       11       H       -0.152378        0.847622         -0.004083
  15.       12       H       -0.163553        0.836447         -0.001505
  16.       13       H       -0.152522        0.847478         -0.002061
  17.       14       H       -0.152939        0.847061         -0.002921
  18.       15       H       -0.162067        0.837933         -0.001498
  19.       16       H       -0.142934        0.857066         -0.003377
  20.       17      Cu        2.444358       31.444358          0.433466
  21.       18       O       -0.470643        7.529357         -0.264487
  22.       19       O       -0.384459        7.615541         -0.338440
  23.       20       C        0.032065        6.032065         -0.009103
  24.       21       N       -0.453973        6.546027          0.129939
  25.       22       N       -0.420439        6.579561          0.126258
  26.       23       N       -0.357597        6.642403          0.101261
  27.       24       H       -0.153066        0.846934          0.003889
  28.       25       H       -0.146049        0.853951          0.005190
  29.       26       H       -0.165527        0.834473          0.002712
  30.       27       H       -0.165050        0.834950          0.003652
  31.       28       H       -0.153038        0.846962          0.004610
  32.       29       H       -0.159174        0.840826          0.003617
  33.       30       H       -0.153214        0.846786          0.002830
  34.       31       H       -0.152314        0.847686          0.003004
  35.       32       H       -0.133289        0.866711          0.004293
复制代码
作者
Author: myzpku    时间: 2023-6-28 17:12
xyz0101 发表于 2023-6-28 14:19
我更想了解的是他之后5.25这个例子是怎么来写cdft部分的？

5.25第一个constraint是表示第一个fragment (F-)有一个excess electron，第二个constraint表示第一个fragment的spin为0



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