计算化学公社

标题: 求助金属络合物激发态能量计算相关问题 [打印本页]

作者
Author: Sally-L 时间: 2025-6-10 21:34
标题: 求助金属络合物激发态能量计算相关问题
本帖最后由 Sally-L 于 2025-6-10 21:35 编辑

各位老师好，最近用高斯计算三重态络合物（Cu2+与有机物络合）的能量，优化和频率计算配体部分用PBE0-D3(BJ)/6-311G(d)，金属部分用PBE0-D3(BJ)/SDD，单点能用PBE0-D3(BJ)/def2-TZVPP，输入文件如下所示，电荷和自旋多重度分别为0 3，但在Gaussian View中打开却变为 0 2，自旋没有triplet的选项，请问这是什么原因，是我的分子结构有问题吗?
%nprocshared=40
%mem=36GB
%chk=Cu-CBBP_triplet_gas_opt.chk
# opt freq PBE1PBE/genecp em=GD3BJ

Cu-CBBP_triplet_gas_opt

0 3
C                2.30937800 0.54721400 -0.17153100
C                3.78209800 0.54741100 -0.28275900
C                4.47470400 1.75901300 -0.34365400
C                5.85649700 1.75446300 -0.45750700
C                6.55833900 0.54532800 -0.48583800
C                5.85846800 -0.66333800 -0.42539200
C                4.47241800 -0.66359500 -0.33364100
O                1.64240000 -0.53818100 -0.12080500
O                1.64234700 1.63191900 -0.12149600
C                8.04845300 0.60018700 -0.66372600
C                8.92286700 -0.43973000 -0.04115100
C                10.16043200 -0.70014600 -0.63880900
C                11.02855700 -1.62697500 -0.07747000
C                10.67824000 -2.28573700 1.10055400
C                9.45759900 -2.01604600 1.71298700
C                8.57757600 -1.10247100 1.14043200
O                8.53723500 1.50224100 -1.32419500
H                3.91953700 2.69002300 -0.30809600
H                6.41137300 2.68358100 -0.53099100
H                6.39355100 -1.60575600 -0.47283600
H                3.91582000 -1.59381600 -0.30024000
H                10.42514800 -0.16306900 -1.54359900
H                11.98225300 -1.83320800 -0.55369600
H                11.36011000 -3.00564200 1.54357500
H                9.18926500 -2.51653800 2.63835800
H                7.63259400 -0.88945000 1.62930500
C                -2.30944700 0.54706900 0.17153100
C                -3.78212800 0.54735600 0.28271600
C                -4.47459500 1.75902100 0.34383300
C                -5.85639000 1.75457900 0.45763800
C                -6.55834500 0.54549700 0.48569000
C                -5.85858000 -0.66322100 0.42502200
C                -4.47252700 -0.66360100 0.33331400
O                -1.64247700 -0.53833400 0.12060200
O                -1.64227100 1.63172900 0.12164000
C                -8.04844800 0.60043000 0.66347900
C                -8.92284800 -0.43963700 0.04111200
C             -10.16025700 -0.70019900 0.63902700
C             -11.02838000 -1.62714400 0.07787400
C             -10.67822200 -2.28587300 -1.10021500
C                -9.45773800 -2.01603800 -1.71289900
C                -8.57771000 -1.10235100 -1.14053100
O                -8.53721200 1.50226600 1.32426400
H                -3.91933700 2.68998600 0.30850100
H                -6.41117900 2.68373600 0.53128500
H                -6.39375300 -1.60559900 0.47225400
H                -3.91603300 -1.59388200 0.29974100
H             -10.42486300 -0.16314000 1.54385900
H             -11.98195300 -1.83349000 0.55429600
H             -11.36009200 -3.00586600 -1.54309100
H                -9.18953000 -2.51650700 -2.63831900
H                -7.63284600 -0.88922900 -1.62958900
Cu             -0.00001000 0.54710200 -0.00004500

C H O 0
6-311G(d)
****
Cu 0
SDD
****

Cu 0
SDD

(, 下载次数 Times of downloads: 9)

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-6-10 21:49
如果关键词啊电荷和自旋多重度啊都是自己写的，别拿GaussView这个简陋的Calculation Setup创建输入文件说事……

鉴于俩有机羧酸根阴离子都是偶数电子，但唯一的Cu2+（[Ar]3d9）是奇数电子，整个体系的自旋多重度只能说2，4，6等偶数。此外这个模型不是很合理，铜和羧酸根的真实配位结构可能是类似醋酸铜[Cu2(OAc)4(H2O)2]的双核结构。

作者
Author: Sally-L 时间: 2025-6-10 22:08

Uus/pMeC6H4-/キ发表于 2025-6-10 21:49
如果关键词啊电荷和自旋多重度啊都是自己写的，别拿GaussView这个简陋的Calculation Setup创建输入文件说事 ...

好的谢谢提醒，确实也有看到铜离子和羧酸的双核结构，在不确定配位结构的时候，是需要将每种可能的结构都算一遍吗？还是说有什么别的方法可以确定络合物结构

作者
Author: wal 时间: 2025-6-10 22:10
由于奇数电子体系自旋多重度是偶数，gview判断出了这一点，贴心地关掉了triplet选项，以免你交上去立即报The combination of multiplicity xx and yy electrons is impossible的错

作者
Author: Sally-L 时间: 2025-6-10 22:16

wal 发表于 2025-6-10 22:10
由于奇数电子体系自旋多重度是偶数，gview判断出了这一点，贴心地关掉了triplet选项，以免你交上去立即报Th ...

好的谢谢，像这种算络合物三重态能量的话，请问自旋多重度应该怎么确定呀

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-6-10 22:23

Sally-L 发表于 2025-6-10 22:16
好的谢谢，像这种算络合物三重态能量的话，请问自旋多重度应该怎么确定呀

怎么可能算奇数电子体系“络合物三重态能量”这种根本都不存在的东西（愣）推荐阅读此文的Q5。“不确定配位结构”又是啥情况，做了实验没表征么，如果一切都要从计算来说明（包括没溶液相的结构得做阴阳离子加溶剂团簇构型搜索，没单晶结构得做晶体结构预测）那工作量可不小。

作者
Author: Sally-L 时间: 2025-6-10 23:53
本帖最后由 Sally-L 于 2025-6-11 00:06 编辑

Uus/pMeC6H4-/キ发表于 2025-6-10 22:23
怎么可能算奇数电子体系“络合物三重态能量”这种根本都不存在的东西（愣）推荐阅读此文的Q5。“不确定配 ...

谢谢分享，确实还没做实验，络合物的电荷除了是0还会有其他情况吗，我发现电荷好像会影响自旋

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-6-11 07:27

Sally-L 发表于 2025-6-10 23:53
谢谢分享，确实还没做实验，络合物的电荷除了是0还会有其他情况吗，我发现电荷好像会影响自旋

这好像是初中化学级别的结论……计算程序能从结构里原子的元素符号（以及赝势使用情况）知道总的核电荷数，但在此基础上要进一步知道总的电子数当然需要净电荷信息，看http://bbs.keinsci.com/thread-266-1-1.html的介绍。络合物整体是什么电荷、铜是几价离子还得看实际情况，不能为了得到结果就随便设嘛。

作者
Author: Sally-L 时间: 2025-6-11 09:42

Uus/pMeC6H4-/キ发表于 2025-6-11 07:27
这好像是初中化学级别的结论……计算程序能从结构里原子的元素符号（以及赝势使用情况）知道总的核电荷数 ...

好的谢谢

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