求助金属络合物激发态能量计算相关问题

Sally-L

各位老师好，最近用高斯计算三重态络合物（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

Uus/pMeC6H4-/キ

如果关键词啊电荷和自旋多重度啊都是自己写的，别拿GaussView这个简陋的Calculation Setup创建输入文件说事……

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

Sally-L

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

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

wal

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

Sally-L

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

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

Uus/pMeC6H4-/キ

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

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

Sally-L

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

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

Uus/pMeC6H4-/キ

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

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

Sally-L

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

好的谢谢