计算化学公社

标题: 阴离子体系加入弥散函数，使用def2-TZVPD基组报错了怎么解决 [打印本页]

作者
Author: sshanshan 时间: 2024-10-31 09:52
标题: 阴离子体系加入弥散函数，使用def2-TZVPD基组报错了怎么解决
已在附件中上传gjf文件，求大佬指点，我是一名计算小白

作者
Author: zjxitcc 时间: 2024-10-31 10:12
本帖最后由 zjxitcc 于 2024-10-31 10:19 编辑

（1）Gaussian没有内置def2TZVPD基组数据，需要用户采用自定义基组格式使用，基组数据可以在https://www.basissetexchange.org找到。如果不熟悉自定义基组，可以看
《详解Gaussian中混合基组、自定义基组和赝势基组的输入》http://sobereva.com/60
（2）以后说报错的同时，要给出报错信息截图（至少15行内容）。
（3）不建议对整个体系所有原子都添加弥散函数，这样耗时较高，SCF收敛和结构优化收敛难度也显著增大。应该先考察一下结构特点，如果体系的负电荷集中在几个原子上，那只需给那几个原子使用弥散函数。
（4）当前结构看上去不是很合理，建议先使用半经验方法（如xTB程序的GFN2-xTB方法）优化一下。不合理的结构一上来就用这么大的基组会浪费大量机时。
（5）仔细阅读此文《谈谈“计算时是否需要加DFT-D3色散校正？”》http://sobereva.com/413

作者
Author: sshanshan 时间: 2024-10-31 16:59

zjxitcc 发表于 2024-10-31 10:12
（1）Gaussian没有内置def2TZVPD基组数据，需要用户采用自定义基组格式使用，基组数据可以在https://www.ba ...

感谢大佬，我去学习一下

作者
Author: JCenter 时间: 2025-5-7 20:57

zjxitcc 发表于 2024-10-31 10:12
（1）Gaussian没有内置def2TZVPD基组数据，需要用户采用自定义基组格式使用，基组数据可以在https://www.ba ...

老师，您好。针对您说的（3），计算单点和使用波函数计算原子电荷时，也是这么考虑吗，还是说需要对整体使用弥散更好。

作者
Author: cokie 时间: 2025-5-7 22:26

JCenter 发表于 2025-5-7 20:57
老师，您好。针对您说的（3），计算单点和使用波函数计算原子电荷时，也是这么考虑吗，还是说需要对整体 ...

• 不加弥散结果一定没法用，必须加弥散：
计算偶极矩、多极矩
计算极化率、超极化率
计算里德堡激发态
计算阴离子体系能量、电子亲和能

谈谈弥散函数和“月份”基组

作者
Author: zjxitcc 时间: 2025-5-7 22:30

JCenter 发表于 2025-5-7 20:57
老师，您好。针对您说的（3），计算单点和使用波函数计算原子电荷时，也是这么考虑吗，还是说需要对整体 ...

仍然成立。

作者
Author: JCenter 时间: 2025-5-9 22:14

cokie 发表于 2025-5-7 22:26
谈谈弥散函数和“月份”基组

你没仔细审题哦，我问的是，整体原子加弥散和单独对类似于羧基氧等集中负电荷的原子加弥散在计算单点能和波函数的区别。

作者
Author: JCenter 时间: 2025-5-9 22:14

zjxitcc 发表于 2025-5-7 22:30
仍然成立。

感谢老师

作者
Author: cokie 时间: 2025-5-9 22:21

JCenter 发表于 2025-5-9 22:14
你没仔细审题哦，我问的是，整体原子加弥散和单独对类似于羧基氧等集中负电荷的原子加弥散在计算单点能和 ...

确实，看花眼了，所答非所问。

计算原子电荷时，局部弥散和整体弥散的区别还是有的。主要体现在和阴离子基团有波函数重叠的原子。我之前计算过含有苯甲酸阴离子基团材料的原子电荷，其中采用 Def2-TZVP/ma-TZVP 和 ma-TZVP 所得到的结果，羧基碳有不到 0.05 左右的差值。在这种级别下，无伤大雅。但 6-31G**/6-31+G** 和 6-31+G** ，羧基碳有将近 0.1 左右的区别，这个还是不能忽略的。

作者
Author: JCenter 时间: 2025-5-9 22:37
是滴我采用的就是def2tzvp/ma-tzvp

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