实验测不出磁矩吗？

个人建议啊，这几个的基组不大，自己算下就是了
特别ecp可以做的大一些

对于整体多重度是4的时候，三个CrIII原子有几种不同的电子排布方式，如果不说明白结果没有可比性。用2015那篇文章的结构算了一下：

Cr Cr(F) Cr 多重度	4 4 4	4 -4 4	2 2 2	4 2 -2	4 4 -4
B3LYP D3BJ def2-SVP SDD（赝势）	0	-0.68	68.73	45.84	-0.26
B3LYP D3BJ def2-TZVP SDD（赝势）	0	-0.67	68.84	45.94
B3LYP D3BJ def2-SVP	0	-0.72
PBE def2-SVP SDD（赝势）	0	-1.76	failed	failed
M06L def2-SVP SDD（赝势）	0	-1.62		failed
PBE0 D3BJ def2-SVP SDD（赝势）	0	-0.44		49.76
PWPB95 D3BJ def2-QZVPP	0	-0.43

可见几个泛函都能给出9个单电子反铁磁性排布的4重态比10重态稳定约1kcal/mol的结果，而有d电子成对的时候会比10重态高几十kcal/mol，题主可以上传一下输出文件或者自己检查自旋布居看看是否这个原因。

其中有些泛函在使用片段组合波函数方法仍无法得出Cr原子只有一个单电子的排布，比如M06L、PBE、M11L、MN15L，刚好他们不会出现B3LYP对应单电子排布时的“异常”情况。

另外2015那篇文章是因为Cr的10重态和4重态（9个单电子，反铁磁性排布）结果差不多，所以使用10重态。“This is done  for the Cr case, and it will be shown that the parallel and antiparallel spin configurations show almost identical results for the interaction energies and frequencies. This allows us to avoid the broken-symmetry DFT approach, where issues such as convergence problems and spin contamination usually appear.”

从结果可以看出这两个态能量差别这么小用哪个都没问题

阿锋001 发表于 2019-7-22 15:34
顺便问一下社长，我之前算 过渡金属体系也有遇到楼主的情况，有时两种自旋多重度下能量相差零点几 kcal/m ...

这两种自旋多重度可能都会同时出现，而且出现比率都不小。这种情况两个自旋多重度之间可能很容易互相自旋翻转。
对于后续反应的计算，可以两种自旋多重度都算算，有可能其中一个从势垒等方面看明显有优势

sobereva 发表于 2019-7-22 01:59
体系也不大，不如用明显高一档的双杂化泛函进行测试

从你已经做的测试来看，实际情况估计四重态和十重态 ...

顺便问一下社长，我之前算 过渡金属体系也有遇到楼主的情况，有时两种自旋多重度下能量相差零点几 kcal/mol，高精度下依然如此，那么（1.）是把低一些的作为基态，还是两者作为简并态都是基态  （2.）那后续的计算是以低的继续，还是两者都可以？

sobereva 发表于 2019-7-22 01:59
体系也不大，不如用明显高一档的双杂化泛函进行测试

从你已经做的测试来看，实际情况估计四重态和十重态 ...

谢谢社长的回复和建议

体系也不大，不如用明显高一档的双杂化泛函进行测试

从你已经做的测试来看，实际情况估计四重态和十重态能量本身就是差不多，现有的普通泛函很难给出绝对有说服力的结论。