QCHEM用约束密度泛函CDFT遇到问题

xyz0101

在qchem中用CDFT来规定片段的电荷量，可以用一下命令来添加 $cdft 2 1 1 25 -1 26 38 $end 但这种方法只设置了不同片段的带电量，从而引起自旋多重度的变化该怎么设置呢？使用手册里没有很好说明。

myzpku

看这页的例子5.24 https://manual.q-chem.com/latest/Ch5.S11.SS4.html

xyz0101

我看到这个例子了
$cdft
2
1 1 3 s
-1 17 19 s
$end
我不清楚这个s设置的是什么含义

xyz0101

myzpku 发表于 2023-6-28 14:05
看这页的例子5.24 https://manual.q-chem.com/latest/Ch5.S11.SS4.html

我更想了解的是他之后5.25这个例子是怎么来写cdft部分的？

myzpku

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

复制代码

myzpku

xyz0101 发表于 2023-6-28 14:19
我更想了解的是他之后5.25这个例子是怎么来写cdft部分的？

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