求助有机金属铜卤化物晶体多激发态发射计算 - 计算化学公社

各位老师好，我在用CP2K计算一个有机金属铜卤化物的晶体，具有0维的电子结构，电子云主要局域在金属铜卤部分，我想计算它的两个激发态（实验测到两个发光峰，经过实验确认应该是同一物质的两个不同激发态）。由于单个晶胞有600多个原子，因此，周期性结构下只做了基态结构优化，算了吸收。然后，想的用簇模型进行激发态结构优化的计算，抠出簇后（大概150个原子），用CP2K的cell方法做了基态的限制性结构优化，固定住周围的有机原子，对发光中心金属铜卤化物先进行基态优化，优化后计算了吸收，结果跟周期性结构下的计算结果能吻合，然后做限制性的激发态结构优化，要算两个态，高能级的态优化后的结果还不错，轨道，能量和振子强度等都没问题，但是第一激发态的结果不太对，算的振子强度基本为0，实际上算吸收的时候是个弱吸收，振子强度大概0.03，先后用了TDDFT和UKS的方法算第一激发态，振子强度都为0，然后TDDFT算的发射能偏高，跟高能级的激发态的结果（TDDFT算的）差不多，UKS算的能量、轨道能对得上。我看了TDDFT算的结果的每步结构优化后的激发能，应该没有发生势能面交叉，优化到其它态上去的情况，现在不知道该用较高的收敛限继续优化，还是换泛函和基组重新算，我的泛函用的是pbe0，用ADMM加速，然后基组用的DZVP，基组本来是想用混合基组的，但是算完吸收跟DZVP差距不太大，反而两个态的跃迁轨道变杂了，然后就都用DZVP了
用CP2K算完两个激发态后，由于第一激发态结果不太对，就用gaussian算了。gaussian算的时候用pbe0，6-31g(d,p)加def2tzvp基组，算完吸收原本CP2K的两个态1，8变成了11,47，给Cu用过SDD基组，结果差不多。对两个态做限制性激发态结构优化，收敛后，激发态能几乎没有降低，跟吸收差不多，然后看了各个态的能量变化，也没有交叉的地方，感觉没法用gaussian算这个，不知道为什么。

问题总结如下：
1、有机金属铜卤化物簇模型计算的第一激发态结构优化能量和振子强度不对，是否提高收敛限续算，还是换泛函和基组重新计算，或者是否有更具针对性的策略。
2、在相同泛函和基组下，gaussian的激发态数目远多于CP2K的cell方法，此外，进行限制性结构优化的激发能能量随着步数基本不降低，即发射能基本等于激发态能，是什么原因