咨询嵌入簇计算用的点电荷数值如何获取？

zhanfei

大家好
  我相咨询两个嵌入簇的问题。
一个为PC区域电荷取值的问题
文献(J. Chem. Phys. 138, 204101 (2013); doi: 10.1063/1.4804607)中有一个NaK3[Fe（II）(Cl)6]-4的体系嵌入簇的例子，此例子为Fe为高自旋。
文献中QC+PC使用的点电荷设置为Fe2.042，Cl -0.804 Na 1.050 K1.052 ；文献中描述点电荷设置原则是在保持QC+PC为电中性原则下，尽量接近拟合出来的CHELP电荷。
我试着用了几个不同大小的团簇，计算结果和文献报道的数值差距都挺大。
第一行为计算使用的电荷和自旋多重度以及基组。第二列为CHELP电荷电荷，最后三列为xyz坐标。设置如下!B3LYP def2-TZVP def2/J Conv TightSCF LargePrint NOPOP

-4 5  def2SVP

Fe     1.401148     0.000000     0.000000     0.000000

Cl    -0.897247    -0.669060    -2.385030     0.405100

Cl    -0.910117     2.392740    -0.500070     0.569970

Cl    -0.893076    -0.376900     0.589490     2.410510

Cl    -0.897333     0.669060     2.385030    -0.405100

Cl   -0.910250    -2.392740     0.500070    -0.569970

Cl    -0.893126     0.376900    -0.589490    -2.410510

-4 5 def2TZVP

Fe     2.247229     0.000000     0.000000     0.000000

Cl    -1.036828    -0.669060    -2.385030     0.405100

Cl   -1.054594     2.392740    -0.500070     0.569970

Cl    -1.031831    -0.376900     0.589490     2.410510

Cl    -1.036978     0.669060     2.385030    -0.405100

Cl    -1.054887    -2.392740     0.500070    -0.569970

Cl    -1.032110     0.376900    -0.589490    -2.410510

+4  5 def2TZVP

Fe    -0.659699     0.000000     0.000000     0.000000

Cl    -0.498829    -0.669060    -2.385030     0.405100

Cl    -0.497521     2.392740    -0.500070     0.569970

Cl    -0.496772    -0.376900     0.589490     2.410510

Cl    -0.498249     0.669060     2.385030    -0.405100

Cl    -0.497101    -2.392740     0.500070    -0.569970

Cl    -0.496162     0.376900    -0.589490    -2.410510

Na     0.935856     2.080930     1.917830     1.996050

Na     0.935626    -2.080930    -1.917830    -1.996050

K     0.962047    -2.080930     3.159750     1.136220

K     0.962387    -3.113240    -1.009050     2.212420

K     0.962254     1.048620    -2.250970     3.072260

K     0.962266     2.080930    -3.159750    -1.136220

K     0.961947     3.113240     1.009050    -2.212420

K     0.961950    -1.048620     2.250970    -3.072260

请问应该如何同时得到相对合理的Fe Cl Na K CHELP电荷数值？例如最后一个计算再加入Na K 后体系Fe电荷数值明显有问题。

第二个问题为PC如何选择的问题

一般文献报道嵌入簇PC区域如何选择的时候不太详细，这个区域形状有什么要求吗？比如以某个原子为球心选择原子，或者以晶胞为单位选择多少晶胞？

感谢

此致