计算化学公社

标题: 求助, 金属酞菁结构优化,偶极矩,LUMO,HOMO计算时,泛函和基组的选择 [打印本页]

作者
Author: dyj_ 时间: 2025-10-29 11:26
标题: 求助, 金属酞菁结构优化,偶极矩,LUMO,HOMO计算时,泛函和基组的选择
本帖最后由 dyj_ 于 2025-10-29 11:27 编辑

请问下, 对于金属酞菁(中心金属一般为Fe,Co,Ni,Cu,Zn)结构优化,偶极矩,LUMO,HOMO计算时,泛函和基组应该如何选择?
尝试过 b3lyp/6-31++g(d,p) em=gd3bj, PBE1PBE/genecp em=gd3bj (6-311++g(d,p),金属用SDD), b3lyp/genecp em=gd3bj (6-311++g(d,p),金属用SDD), b3lyp/ma-TZVP em=gd3bj, b3lyp/ma-SVP em=gd3bj, PBE1PBE/ma-TZVP em=gd3bj, PBE1PBE/ma-SVP em=gd3bj,

数据来看, b3lyp/6-31++g(d,p) em=gd3bj和b3lyp/ma-TZVP em=gd3bj结果接近,并且与实验参数也接近, b3lyp/ma-TZVP em=gd3bj 较b3lyp/6-31++g(d,p) em=gd3bj时间增加了大约10倍.

已发表文章中有用b3lyp/6-31g(d,p) em=gd3bj和b3lyp/def2-SVP em=gd3bj, 没有明显说明是否添加弥散. 是否说后面计算可以直接全部用 b3lyp/6-31++g(d,p) em=gd3bj? 图是论坛中的,仅作示例.

作者
Author: 54687 时间: 2025-10-29 12:37
对于酞菁这类大共轭体系，结构优化时泛函用wB97XD可能更好（尽管wB97XD是给主族参数化的，但计算含过渡金属的体系结果并不是总是那么差），计算HOMO-LUMO gap则用B3LYP-D3(BJ)，具体参考
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）
对于不同任务而言，基组该如何选择，是否要加弥散函数，看下文
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）

作者
Author: dyj_ 时间: 2025-10-29 12:51
本帖最后由 dyj_ 于 2025-10-29 17:28 编辑

54687 发表于 2025-10-29 12:37
对于酞菁这类大共轭体系，结构优化时泛函用wB97XD可能更好（尽管wB97XD是给主族参数化的，但计算含过渡金属 ...

结构优化用的是wB97XD泛函,def2tzvp基组. 能量计算用B3LYP(d3bj), 卢老师有个回复说用PBE0我也试了.现在更想确定基组用6-31+g(d,p)是否合适.

作者
Author: sobereva 时间: 2025-10-30 08:54
除了算偶极矩，你这体系根本不需要加弥散。仔细看
谈谈弥散函数和“月份”基组
http://sobereva.com/119

6-31G**即便加了弥散，算偶极矩仍然是很差的选择。适合用什么下文说了
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）

作者
Author: wzkchem5 时间: 2025-10-30 09:05
你确定要算偶极矩，不是跃迁偶极矩？
这个分子除了巯基以外是完全对称的，只要巯基的取向是对称的，偶极矩就精确为0，不需要算

作者
Author: dyj_ 时间: 2025-10-30 09:26

wzkchem5 发表于 2025-10-30 09:05
你确定要算偶极矩，不是跃迁偶极矩？
这个分子除了巯基以外是完全对称的，只要巯基的取向是对称的，偶极矩 ...

抱歉,没有说清楚. 我想咨询的是金属酞菁衍生物这一大类的分子, 图片是从论坛里找的一个金属酞菁分子. 金属酞菁衍生物需要算偶极矩的情况是不对称的分子.

作者
Author: dyj_ 时间: 2025-10-30 09:31

sobereva 发表于 2025-10-30 08:54
除了算偶极矩，你这体系根本不需要加弥散。仔细看
谈谈弥散函数和“月份”基组
http://sobereva.com/119
...

谢谢社长解惑,请问下,这样是否可行: 结构优化使用wb97xd/def2-svp(或者wb97x-d3/def2-svp), 单点能: 不算偶极矩使用b3lyp/def2-tzvp, 算偶极矩使用b3lyp/ma-tzvp?

作者
Author: sobereva 时间: 2025-10-31 06:02

dyj_ 发表于 2025-10-30 09:31
谢谢社长解惑,请问下,这样是否可行: 结构优化使用wb97xd/def2-svp(或者wb97x-d3/def2-svp), 单点能: 不算 ...

可以
但计算偶极矩的话，def2-TZVPD是比ma-TZVP精度更好的选择

作者
Author: dyj_ 时间: 2025-11-5 10:42

sobereva 发表于 2025-10-31 06:02
可以
但计算偶极矩的话，def2-TZVPD是比ma-TZVP精度更好的选择

好的,感谢社长.

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