计算化学公社

标题: 高斯中使用双杂DSD-PBEP86优化报错 [打印本页]

作者
Author: Luckynan 时间: 2022-11-18 15:21
标题: 高斯中使用双杂DSD-PBEP86优化报错
本帖最后由 Luckynan 于 2022-11-18 15:21 编辑

在使用高斯16优化结构时，使用双杂化泛函DSD-PBEP86计算报错。关键词#popt DSDPBEP86/6-31g(d) DENSITY，想请教一下，该如何修改？是什么原因？输入文件及输出文件如下。

作者
Author: wzkchem5 时间: 2022-11-18 15:54
先不论报错原因，这个计算即使正常结束结果也是不能用的，因为双杂化泛函至少要搭配3-zeta基组。双杂化泛函搭配2-zeta基组，结果还不如普通泛函搭配3-zeta基组靠谱。

作者
Author: niobium 时间: 2022-11-18 15:55
结构优化没有必要用双杂化泛函

作者
Author: wjc404 时间: 2022-11-18 16:26
看倒数第二行，硬盘使用超过了40GW(320GB)，可能是超过了限制然后进程被杀掉了。这种体系用高斯算双杂化泛函的梯度有些勉为其难，硬要算可以尝试用ORCA开RI算。

作者
Author: sobereva 时间: 2022-11-18 18:09
写DENSITY毫无意义

你这种体系用PBE0-D3(BJ)结合6-311G*优化得妥妥的，用双杂化结合6-31G*结果只会更烂，耗时还多两个数量级，而且rwf文件还会占极大的硬盘，临时文件目录写满了任务就会崩溃

作者
Author: Luckynan 时间: 2022-11-18 20:51

sobereva 发表于 2022-11-18 18:09
写DENSITY毫无意义

你这种体系用PBE0-D3(BJ)结合6-311G*优化得妥妥的，用双杂化结合6-31G*结果只会更烂 ...

好的，谢谢老师。还想问一下，我主要计算化合物的CT激发，使用cam-b3lyp和wb97xd均无法很好的符合实验值，PBE0-D3(BJ)的吸收值也会比实验值蓝移25nm左右，我有必要使用自定义泛函嘛？

作者
Author: wzkchem5 时间: 2022-11-18 21:22

Luckynan 发表于 2022-11-18 13:51
好的，谢谢老师。还想问一下，我主要计算化合物的CT激发，使用cam-b3lyp和wb97xd均无法很好的符合实验值 ...

25nm已经属于TDDFT误差范围以内了（当然，如果你算的是极紫外区的吸收，那还可以指望误差做得比25nm小），一般来说没必要也很难算得更准，即使更准也是碰巧。
一定要算得更准的话，第一必须做振动分辨光谱，第二必须显式考虑溶剂，第三必须做构象平均（包括溶剂分子的构象），第四激发能必须在DLPNO-STEOM-CCSD或更高级别下计算，TDDFT是不行的。
此外，如果你只研究吸收的话，结构优化和TDDFT的理论级别可以不统一，不要因为看到双杂化泛函算TDDFT准，就也用双杂化泛函做结构优化。

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