计算化学公社

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

作者
Author: mayyLYF 时间: 2023-1-12 21:52
标题: 高斯计算报错，在计算跃迁电场时想要在坐标后同时输入两个原子半径，应该如何书写
因为想要用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原子的半径呢？

作者
Author: wzkchem5 时间: 2023-1-12 22:17
这种从跃迁电场反推非绝热耦合矩阵元的做法等价于用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/

作者
Author: mayyLYF 时间: 2023-1-13 16:50

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

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

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