likjvau 发表于 2022-3-3 02:26
那老师，请问我帖子描述的计算思路有什么错误么，我在找到了基态自旋后，激发态计算直接固定了这个自旋4 ...

那叫<S**2>，不叫<s**>，那个**是乘方的意思，所以S**2就是S平方算符，<S**2>是S平方算符的期望。
振子强度是近似等于0还是精确等于0？如果只是近似等于0，那是正常的，稀土的f-f跃迁本来就是禁阻跃迁，振子强度接近0（但不严格等于0）

zjxitcc 发表于 2022-3-3 09:17
singlet、triplet、50-50只能应用于闭壳层单重态，当你自旋多重度不是1时，这三个关键词全失效，算出什么 ...

那老师，请问我帖子描述的计算思路有什么错误么，我在找到了基态自旋后，激发态计算直接固定了这个自旋4，但是计算出所有值振子强度全部为0，设置了开壳层计算，<s**>也全部为3.75左右

likjvau 发表于 2022-3-3 08:13
老师我还想问一下，激发态计算时，设置的自旋多重度一定要跟基态自旋保持一样么，我设置自旋为3，计算出 ...

singlet、triplet、50-50只能应用于闭壳层单重态，当你自旋多重度不是1时，这三个关键词全失效，算出什么不可控（哪些态能量低，就会被算出来）

wzkchem5 发表于 2022-3-3 04:44
你截图已经说了，不是把stable=opt去掉就万事大吉了，应该先只加stable=opt不加td算一遍，得到稳定波函数以 ...

老师我还想问一下，激发态计算时，设置的自旋多重度一定要跟基态自旋保持一样么，我设置自旋为3，计算出的所有激发能对应的<s**>都在2.00左右，设置自旋为4，计算出的所有激发能的<s**>都在3.75左右,不太清楚为什么计算什么态不是根据td=singlet、triplet、50-50决定的么

wzkchem5 发表于 2022-3-3 04:44
你截图已经说了，不是把stable=opt去掉就万事大吉了，应该先只加stable=opt不加td算一遍，得到稳定波函数以 ...

好的，感谢老师指导，我先算下试试

你截图已经说了，不是把stable=opt去掉就万事大吉了，应该先只加stable=opt不加td算一遍，得到稳定波函数以后，用guess=read加td但不加stable=opt，读取stable=opt不加td的波函数做SCF。不是说把stable=opt去掉、只写td就行了的。
这个体系大核赝势优化出来的结构是否靠谱我拿不准，尤其不知道你配合物的结构，建议你自己查一下类似的配合物别人用大核赝势优化的结构和实验晶体结构的符合情况如何，如果符合情况可以接受，就可以用。