计算化学公社

标题: 咨询嵌入簇计算用的点电荷数值如何获取？ [打印本页]

作者
Author: zhanfei 时间: 2018-2-5 21:40
标题: 咨询嵌入簇计算用的点电荷数值如何获取？
本帖最后由 zhanfei 于 2018-2-5 21:41 编辑

大家好
我相咨询两个嵌入簇的问题。
一个为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区域如何选择的时候不太详细，这个区域形状有什么要求吗？比如以某个原子为球心选择原子，或者以晶胞为单位选择多少晶胞？

感谢

此致

欢迎光临计算化学公社 (http://bbs.keinsci.com/)