如何限定过渡金属双原子体系的电子状态

pikachuupup

各位老师好，我有一个困惑，假如我有两个过渡金属原子组成的体系，每个金属原子的单电子都为3，那么计算时双金属体系的多重度就应该设置为7，即单电子总数为6。我想要得到每个金属上单电子数均为3的高自旋态能量（非对称性破缺态）。
但是在计算过程中，设置自旋多重度为7后，通过SPIN POPULATIONS发现SCF很容易收敛到1,5（即第一个原子的单电子数为1，第二个原子单电子数为5）等等。对于高自旋的多金属体系来说，其SCF解并不唯一且稳定的解也不一定唯一，那么应该如何限定才能得到想要的电子状态3,3呢？
我曾经尝试过Brokensym关键词限定每个位点的单电子数量，但是最终计算结果给出的是对称性破缺态的能量而并非High Spin态的能量，使用Flipspin好像也是同样的道理，有没有什么关键词或手段能够达到限定体系的电子状态的目的呢？

hebrewsnabla

具体是什么体系？

pikachuupup

hebrewsnabla 发表于 2023-7-20 18:14
具体是什么体系？

是两个Ir原子组成的体系，即Ir2，我想获得高自旋态下的能量

sobereva

“每个金属原子的单电子都为3，那么计算时双金属体系的多重度就应该设置为7” 逻辑不成立
俩金属原子之间又不是没成键作用

pikachuupup

sobereva 发表于 2023-7-20 19:12
“每个金属原子的单电子都为3，那么计算时双金属体系的多重度就应该设置为7” 逻辑不成立
俩金属原子之间 ...

是的，这个地方我说的不严谨，如果我的目标就是计算高自旋态3,3的能量，请问sob老师有实现限定电子状态的办法吗

sobereva

pikachuupup 发表于 2023-7-20 19:17
是的，这个地方我说的不严谨，如果我的目标就是计算高自旋态3,3的能量，请问sob老师有实现限定电子状态的 ...

随手用高斯算了下，直接就得到各个Ir自旋布居都是3的结果
%chk=C:\gtest\1.chk
#p tpssh sdd

Title Card Required

0 7
Ir 0. 0. 0.
Ir 0. 0. 1.8

hebrewsnabla

我随手用orca算了一下（2.0 angstrom，PBE0/def2-svp)，也是得到两个原子自旋布居都是3的结果。

pikachuupup

sobereva 发表于 2023-7-20 19:21
随手用高斯算了下，直接就得到各个Ir自旋布居都是3的结果
%chk=C:\gtest\1.chk
#p tpssh sdd

谢谢老师，那如果又想得到1,5这样的解，应该如何实现呢，对这个问题一直有疑惑

sobereva

pikachuupup 发表于 2023-7-20 19:37
谢谢老师，那如果又想得到1,5这样的解，应该如何实现呢，对这个问题一直有疑惑

谈谈约束性DFT(CDFT)
http://sobereva.com/271（http://bbs.keinsci.com/thread-535-1-1.html）

wzkchem5

用片段初猜方法；或者先收敛到1,5的态以后通过手动交换占据轨道和虚轨道，使得电子重新填充到正确的轨道上，从而得到你需要的波函数

pikachuupup

sobereva 发表于 2023-7-20 19:55
谈谈约束性DFT(CDFT)
http://sobereva.com/271（http://bbs.keinsci.com/thread-535-1-1.html）

谢谢sob老师，我去学习学习

pikachuupup

wzkchem5 发表于 2023-7-20 20:08
用片段初猜方法；或者先收敛到1,5的态以后通过手动交换占据轨道和虚轨道，使得电子重新填充到正确的轨道上 ...

谢谢您，我去学习一下