高斯使用TPSS0-D3报错

baidie · 发表于 Post on 2022-8-2 14:38:03

#p TPSSTPSS/gen IOp(3/76=0750002500) empiricaldispersion=gd3bj pseudo=read scrf=(smd,solvent=Methanol)在使用TPSS0-D3（含25%HF交换能）此方法的时候，高斯16C输入命令如上
会出现如下报错信息：
Two-electron integral symmetry is turned on.
R6DS8: Unable to choose the S8 parameter, IExCor= 2523 IXCFnc= 0 ScaHFX= 0.250000 IDFTD=4

请问友友们我是哪里写错了呢？该怎么改正呢？

喵星大佬 · 发表于 Post on 2022-8-2 15:38:57

本帖最后由喵星大佬于 2022-8-2 15:40 编辑

因为你没有提供TPSS0的D3参数，虽然你用IOp3/76指定了HF的比例，但是自定义泛函肯定是不会有内置的D3参数的，所以也是要给的

根据文献Phys. Chem. Chem. Phys. 2017, 19, 32184的SI里给的参数(实际来源是J. Comput. Chem. 2011, 32, 1456)
应该用下面的IOp来指定

IOp(3/174=1000000,3/175=1257600,3/177=376800,3/178=4586500) em=gd3bj

复制代码

baidie · 发表于 Post on 2022-8-2 16:12:51

喵星大佬发表于 2022-8-2 15:38
因为你没有提供TPSS0的D3参数，虽然你用IOp3/76指定了HF的比例，但是自定义泛函肯定是不会有内置的D3参数的 ...

谢谢解答那请问25%HF交换能的表示是直接在 IOp(3/174=1000000,3/175=1257600,3/177=376800,3/178=4586500) 后面加上3/76=0750002500吗？除此之外 TPSS0方法在高斯中输入命令是应该写成TPSS0还是TPSSTPSS呢？

喵星大佬 · 发表于 Post on 2022-8-2 16:52:51

本帖最后由喵星大佬于 2022-8-2 17:06 编辑

http://bbs.keinsci.com/thread-3997-1-1.html

#p tpsstpss IOp(3/76=1000002500,3/77=0750007500) 6-31G* IOp(3/174=1000000,3/175=1257600,3/177=376800,3/178=4586500) em=gd3bj

复制代码

TPSS0在高斯里并没有内置，直接写显然报错

baidie · 发表于 Post on 2022-8-2 19:16:56

谢谢您！
还有个问题是为什么只能用em=gd3bj，换成em=gd3就会报错呢？

喵星大佬 · 发表于 Post on 2022-8-2 20:17:45

baidie 发表于 2022-8-2 19:16
谢谢您！
还有个问题是为什么只能用em=gd3bj，换成em=gd3就会报错呢？

因为gd3是0阻尼形式，需要定义的参数跟BJ阻尼形式不一样

baidie · 发表于 Post on 2022-8-5 08:52:29

如果是用D3的话，那请问输入文件的哪些地方需要修改呢？

sobereva · 发表于 Post on 2022-8-5 09:12:52

baidie 发表于 2022-8-5 08:52
如果是用D3的话，那请问输入文件的哪些地方需要修改呢？

除非用的是G09 D.01，否则对TPSS0没任何理由用零阻尼的D3

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 高斯使用TPSS0-D3报错

评分 Rate