计算化学公社

标题: 高斯使用TPSS0-D3报错 [打印本页]

作者
Author: baidie 时间: 2022-8-2 14:38
标题: 高斯使用TPSS0-D3报错
#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

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

作者
Author: 喵星大佬 时间: 2022-8-2 15:38
本帖最后由喵星大佬于 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

复制代码

作者
Author: baidie 时间: 2022-8-2 16:12

喵星大佬发表于 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呢？

作者
Author: 喵星大佬 时间: 2022-8-2 16:52
本帖最后由喵星大佬于 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在高斯里并没有内置，直接写显然报错

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

作者
Author: 喵星大佬 时间: 2022-8-2 20:17

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

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

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

作者
Author: sobereva 时间: 2022-8-5 09:12

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

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

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