关于激发态优化无法收敛的问题

liuxiao.shaw

各位大佬，我正在优化图片1中的分子，优化S1使用的方法是# opt td b3lyp/6-31g(d) nosymm em=gd3bj 电荷/自旋多重度分别为 4 1；优化T1使用的方法是# opt ub3lyp/6-31g(d) nosymm em=gd3bj 电荷/自旋多重度分别为 4 3。在优化过程中碰到几个问题：

1. UKS优化后stable=opt检测三线态的波函数稳定性显示收敛, <S**2>=2.600，算Tn能量显示为负 (换了几种泛函都出现了这种情况)；检测单线态波函数稳定性显示不收敛，stable=opt优化后波函数收敛，然后guess=read计算至triplet,nstates=15能量，结果显示的<S**2>不为整数，只有零点几。请问这种情况的Tn能量是准确的吗？

2. S1激发态优化过程中出现能量为负的情况，使用stable=opt优化波函数后用guess=read，继续优化S1结构；重新优化时前几轮S1为能量为正，但是比上面算出来的T1的能量还要低，优化几轮之后又出现激发态能量为负的情况，然后stable=opt优化波函数；如此往复几次之后，S1优化时能量仍然出现负值，而且结构优化不收敛，请问这种情况该如何解决？尝试过scf=xqc，没有效果。

请各位大佬帮帮忙！

xsc6

没太看懂你的目的是什么，一般来说如果你算的东西如果同时要用到S1和T1，建议是都用同一种方法算可能会比较好

wzkchem5

你这个分子确定不需要做显式溶剂化吗？

liuxiao.shaw

wzkchem5 发表于 2024-11-12 10:05
你这个分子确定不需要做显式溶剂化吗？

老师您好，目前没加溶剂，如果加水的话会好一些吗？只给胍基那边加显式水可以吗？

wzkchem5

liuxiao.shaw 发表于 2024-11-12 16:15
老师您好，目前没加溶剂，如果加水的话会好一些吗？只给胍基那边加显式水可以吗？

胍基和溶剂有强烈相互作用，如果实验加了溶剂，计算必须加溶剂。和胍基直接成氢键的溶剂显式描述，其他溶剂隐式描述

liuxiao.shaw

wzkchem5 发表于 2024-11-13 09:59
胍基和溶剂有强烈相互作用，如果实验加了溶剂，计算必须加溶剂。和胍基直接成氢键的溶剂显式描述，其他溶 ...

好的，谢谢老师，我先试一下