24小时热门版块排行榜

返回列表

lalasiji

捐助贵宾 (初入文坛)

应助: 0 (幼儿园)
金币: 154
帖子: 4
在线: 6.6小时
虫号: 4144147
注册: 2015-10-14
专业: 理论和计算化学

[求助] 用gaussian计算原子ESP或RESP电荷，如何指定原子半径？

计算一些含多个Fe Co Cu Cr Mo W等金属MOF的ESP电荷，不知道如何指定原子半径？
原子不同或者使用的基组（比如用SDD或者SBKJC基组）不同的话，原子半径是否也不相同？
另外，如果想计算它们的RESP charge，是否直接用Gaussian的输出文件，
antechamber -i MOF.log -fi gout -o MOF.mol2 -fo mol2 -c resp -nc ${total chage} -s 2
就可以了？我已经试过antechamber，但得到的结果比较奇怪，不知道是哪儿的问题？
我的gussian(gjf)输入文件是：
%Mem=4GB
%nprocshared=12
%chk=test.chk
#P B3LYP/Gen Integral=Ultrafine Pseudo=Read SCF(Conver=8,XQC,MaxConventionalCycles=64,MaxCycle=128) NoSymm Test
5D 7F Pop=(chelpg,ReadRadii,Regular) IOp(6/33=2) GFinput GFprint

电荷自旋多重度
原子坐标
... ...
P C N O 0
6-31G(d)
****
Fe Cr W 0
CEP-31G
****
Fe Cr W 0
CEP-31G

Fe  1.8
Cr  1.8
W  1.8
另外，我也没有重复出别人报导同一分子,都采用SBKJC+(df)的电子能量:我的结果是-4000 au，他们的是-1000au
不知道我的gjf文件有无明显问题？
非常感谢！

回复此楼