计算化学公社

标题: 求助B3LYP*泛函中色散校正的输入方法 [打印本页]

作者
Author: Shejon 时间: 2022-6-30 15:57
标题: 求助B3LYP*泛函中色散校正的输入方法
本帖最后由 Shejon 于 2022-6-30 16:34 编辑

各位老师好，我在尝试优化一个Fe3+与没食子酸的单核配合物构型，查到相关文献有用B3LYP*泛函（HF15%的泛函），但是发现em=gd3会导致R6DS8: Unable to choose the S8 parameter错误，查阅论坛以后发现是使用 em=gd3时程序未存有相应泛函的D3色散校正参数导致的。
在Gaussian官网查阅发现G16里面是有相关参数定义的，原文：The table below gives the list of functionals in Gaussian 16 for which GD3 parameters are defined. For these functionals, dispersion is requested with EmpiricalDispersion=GD3，同时查到 IOps(3/174-176) 可以加入色散校正，但是本人不知道如何输入了，在此请了解的老师告诉我一下应该如何在这个任务中加入色散校正，谢谢。纯新手，第一次发帖多有冒犯还请各位多多指点，谢谢了。

使用版本Gaussian16 Revision C.01
任务信息
# opt genecp blyp iop(3/76=1000001500,3/77=0720008500,3/78=0810010000)

-1 6
O                -1.58300000 -0.26300000 -0.75500000
O                -6.51700000 -0.98200000 -3.98700000
C                -4.31400000 -0.52900000 -3.01300000
C                -5.71900000 -0.54700000 -2.99100000
O                0.58300000 -1.51400000 0.66500000
O                2.42700000 -6.14600000 3.89600000
C                1.71800000 -4.01100000 2.92200000
C                2.40500000 -5.23800000 2.90000000
O                7.08100000 -0.26300000 -0.75500000
O                2.14700000 -0.98200000 -3.98700000
C                4.35000000 -0.52900000 -3.01300000
C                2.94500000 -0.54700000 -2.99100000
O                0.87000000 -1.34800000 2.86500000
O                3.96000000 -5.26100000 -0.36700000
O                3.79800000 -6.75000000 1.76500000
C                1.01500000 -1.93000000 1.76500000
C                1.76700000 -3.23200000 1.76500000
C                2.46600000 -3.58000000 0.60700000
C                3.18400000 -4.78800000 0.63000000
C                3.12100000 -5.57800000 1.76500000
O                -7.72600000 -0.09700000 -1.85600000
C                -2.15900000 -0.09700000 -1.85600000
C                -3.66300000 -0.09700000 -1.85600000
C                -6.37200000 -0.09700000 -1.85600000
O                -1.58300000 0.06800000 -2.95600000
C                -4.31400000 0.33400000 -0.69800000
C                -5.71900000 0.35200000 -0.72000000
O                -6.51700000 0.78700000 0.27600000
O                -3.74900000 5.99000000 0.66500000
O                -1.90500000 1.35700000 3.89600000
C                -2.61400000 3.49200000 2.92200000
C                -1.92700000 2.26500000 2.90000000
O                0.93800000 -0.09700000 -1.85600000
C                6.50500000 -0.09700000 -1.85600000
C                5.00100000 -0.09700000 -1.85600000
C                2.29200000 -0.09700000 -1.85600000
O                7.08100000 0.06800000 -2.95600000
Fe                0.19900000 0.32300000 -0.04500000
C                4.35000000 0.33400000 -0.69800000
C                2.94500000 0.35200000 -0.72000000
O                2.14700000 0.78700000 0.27600000
O                -3.46200000 6.15500000 2.86500000
O                -0.37200000 2.24200000 -0.36700000
O                -0.53400000 0.75300000 1.76500000
C                -3.31700000 5.57400000 1.76500000
C                -2.56500000 4.27100000 1.76500000
C                -1.86600000 3.92300000 0.60700000
C                -1.14800000 2.71500000 0.63000000
C                -1.21100000 1.92500000 1.76500000
H                2.48800000 -2.95500000 -0.27300000
H                1.15900000 -3.70300000 3.79900000
H                -1.85700000 4.55100000 -0.28400000
H                -3.17500000 3.79200000 3.80600000
H                -3.77900000 0.66200000 0.18100000
H                -3.78100000 -0.86400000 -3.89700000
H                4.89000000 -0.86700000 -3.89600000
H                4.89000000 0.65500000 0.19300000
H                3.65000000 -7.22000000 2.60600000
H                1.88400000 -5.86900000 4.65500000
H                3.67400000 -4.95200000 -1.23800000
H                -6.04200000 0.88200000 1.11600000
H                -8.06200000 0.22700000 -1.00000000
H                -7.45600000 -0.91800000 -3.73000000
H                -3.96500000 6.96600000 2.62000000
H                8.03500000 0.02700000 -2.71200000
H                -2.43000000 1.60600000 4.68500000
H                2.62500000 -1.32700000 -4.77000000
H                -0.22319757 2.66727165 -1.21470388
H                2.68717434 1.10463281 1.00326961

C O H 0
6-31+G**
****
Fe 0
LANL2DZ
****

Fe 0
Lanl2dz

作者
Author: wzkchem5 时间: 2022-6-30 16:38
首先得查到B3LYP*的色散校正参数，只有从文献查到这个泛函的色散校正参数，才能谈怎么从输入文件输入的问题。如果查不到，建议改成其他类似泛函，比如TPSSh等

作者
Author: sobereva 时间: 2022-6-30 20:06
B3LYP*这种非主流泛函基本不可能有人给其专门搞D3参数。优化这个体系也没必要带D3，显然色散作用不会明显影响体系构型。
优化这种体系也没必要用这个，TPSSh、PBE0通常都是更好选择

作者
Author: Shejon 时间: 2022-7-1 08:52

wzkchem5 发表于 2022-6-30 16:38
首先得查到B3LYP*的色散校正参数，只有从文献查到这个泛函的色散校正参数，才能谈怎么从输入文件输入的问题 ...

感谢回答，明白了，我换一个泛函吧。

作者
Author: Shejon 时间: 2022-7-1 08:52

sobereva 发表于 2022-6-30 20:06
B3LYP*这种非主流泛函基本不可能有人给其专门搞D3参数。优化这个体系也没必要带D3，显然色散作用不会明显影 ...

谢谢sob老师，我换个泛函。

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