计算化学公社

标题: g16 包含碘原子的有机分子计算ESP电荷报错 no radius [打印本页]
作者
Author: chemcomp    时间: 2025-10-15 13:57
标题: g16 包含碘原子的有机分子计算ESP电荷报错 no radius
本帖最后由 chemcomp 于 2025-10-15 15:15 编辑

各位老师好，我在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.
复制代码









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

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



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3