高斯计算报错，在计算跃迁电场时想要在坐标后同时输入两个原子半径，应该如何书写

mayyLYF

因为想要用momap计算速率，最近在学习有关包含金属原子的跃迁电场计算，在计算过程中输入文件为：
%nprocshared=16
%chk=Cu6-lanl2dz-6-31gd-nacme.chk
%mem=2000MW
#p td PBE1PBE/genecp prop=(fitcharge,field) iop(6/22=-4,6/29=1,6/30=0,6/17=2) scrf=(solvent=Toluene,pcm) nosymm pop=(CHELPG,ReadRadii)

PABBUT

0 1
Cu                -0.56079500    0.28173600    0.06384400

... ...
Br                 6.15374700    0.07728900   -0.21110900

C N H O Br 0
6-31g*
****
Cu 0
lanl2dz
****

Cu 0
lanl2dz

Cu 1.4


结果文件报错信息：
  37     1    1.45
   38     6    1.50
   39     1    1.45
   40     6    1.50
   41     1    1.45
   42     6    1.50
   43    35    0.00
GetVDW:  no radius for atom  43 atomic number  35.
Error termination via Lnk1e in /soft/install/gauss/g16c01/g16/l602.exe at Thu Jan 12 10:18:10 2023.
Job cpu time:       0 days  2 hours  8 minutes 20.6 seconds.
Elapsed time:       0 days  0 hours  8 minutes  5.0 seconds.
File lengths (MBytes):  RWF=   1695 Int=      0 D2E=      0 Chk=     16 Scr=      1

其中43号原子是Br原子，我这里在输入文件中只加入了金属的半径，根据这个报错信息应该也需要Br的原子半径，想求助下各位大佬有没有知道应该怎么在加上金属半径的同时加上Br原子的半径呢？

wzkchem5

这种从跃迁电场反推非绝热耦合矩阵元的做法等价于用Hellmann-Feynman定理来做计算，非常不准，而且因为高斯16现在已经支持了更准确的做法（用NAC关键字计算含有Pople项的非绝热耦合矩阵元），这种做法已经过时了。
关于为什么基于Hellmann-Feynman定理计算的非绝热耦合矩阵元不准，参见https://pubs.acs.org/doi/10.1021/acs.accounts.1c00312?ref=PDF
关于如何用Gaussian 16计算含有Pople项的非绝热耦合矩阵元，参见http://gaussian.com/td/

mayyLYF

wzkchem5 发表于 2023-1-12 22:17
这种从跃迁电场反推非绝热耦合矩阵元的做法等价于用Hellmann-Feynman定理来做计算，非常不准，而且因为高斯 ...

好的，非常感谢，我了解一下相关文献