g16 包含碘原子的有机分子计算ESP电荷报错 no radius

chemcomp

各位老师好，我在Gaussian中，使用如下 gjf 文件，对一个包含碘原子的有机分子进行优化和ESP电荷计算，得到了下文所示的报错信息。

请各位老师帮忙看看是什么问题，应当如何修改，谢谢！

gjf 文件：

--Link1--
%nprocshared=16
%chk=ligand
%mem=16GB
#HF/gen pseudo=read opt SCF=tight Test Pop=(MK,ReadRadii) iop(6/33=2,6/42=6,6/50=1) nosymm

remark line goes here

0   1
    C    1.2620000000       -0.0460000000        1.0430000000     
    H    1.3824000000        0.8136000000        1.7031000000
    xxxx
    I    3.0060000000        0.0110000000       -0.2860000000     

C H 0
6-31G*
****
I 0
LANL2DZ
****

I 0
LANL2DZ

I 1.99

out_resp.gesp




报错：

1. Merz-Kollman atomic radii used.
   
2. Read replacement radii for each type of atom:
   
3. Atom Element Radius
   
4.     1     6    1.50
   
5.     2     1    1.20
   
6.     3     1    1.20
   
7.     4     6    1.50
   
8.     5     1    1.20
   
9.     6     1    1.20
   
10.     7     6    1.50
    
11.     8     1    1.20
    
12.     9     1    1.20
    
13.    10     6    1.50
    
14.    11     6    1.50
    
15.    12     6    1.50
    
16.    13     6    1.50
    
17.    14     1    1.20
    
18.    15     6    1.50
    
19.    16     1    1.20
    
20.    17     6    1.50
    
21.    18     1    1.20
    
22.    19     1    1.20
    
23.    20     1    1.20
    
24.    21    53    0.00
    
25. <font color="#ff0000"> GetVDW:  no radius for atom  21 atomic number  53.</font>
    
26. Error termination via Lnk1e in /software/g16/l602.exe at xxx.
    

复制代码

sobereva

用Multiwfn算RESP电荷是最好的做法，没有之一
仔细看
RESP拟合静电势电荷的原理以及在Multiwfn中的计算
http://sobereva.com/441（http://bbs.keinsci.com/thread-10880-1-1.html）

直接用fch文件当输入文件。提示缺半径时直接用Multiwfn在屏幕上推荐的半径