你的第一个计算文件在第6步的时候梯度飞了，随后振子强度突然走低，感觉是B3LYP低估了某电荷转移态能量，随后follow歪了，最后跟着ghost态走到了一个本来不应该有的圆锥交叉

用你的初始结构换cam-B3LYP-d3bj/6-31g*，8步收敛，看起来一切正常。 SNMI-S1-opt-3.zip (256.45 KB, 下载次数 Times of downloads: 1)

2025-5-10 23:36 上传 Uploaded

点击下载
Click to download

S1振子强度挺高的，跟原先疑似走歪的近0振子强度大相径庭。S3这个电荷转移态看着有点可疑，可能是follow它最后follow到圆锥交叉里去了


这也是我面对陌生体系算激发任务不习惯B3LYP or PBE0的原因。。一不小心follow到ghost态可能就不知道给你带哪去了

wal 发表于 2025-5-11 12:28
建议cam-B3LYP-D3BJ起手，这也是我的习惯。如果发现感兴趣的态是LE激发，高估激发能太厉害，此时用已收敛 ...

好滴感谢w老师，学到了

czjtry 发表于 2025-5-11 12:11
谢谢w老师，狠狠学习辽，我是实验组自学的计算，类似的这种有机小分子我一直是使用b3lyp/6-31G**，那激发 ...

wB97XD可以当准万能泛函使，基态和激发态的优化都可以用。
泛函的选择可以看《简谈量子化学计算中DFT泛函的选择》（http://sobereva.com/272）

我也是实验组的哈哈哈，我和wal是科音的高级量子化学班的同学

czjtry 发表于 2025-5-11 12:11
谢谢w老师，狠狠学习辽，我是实验组自学的计算，类似的这种有机小分子我一直是使用b3lyp/6-31G**，那激发 ...

建议cam-B3LYP-D3BJ起手，这也是我的习惯。如果发现感兴趣的态是LE激发，高估激发能太厉害，此时用已收敛结构切换PBE0 B3LYP之类的重新算也不会如从头开始算那样耗时很久。
也可以考虑wb97xd，这个泛函跟cam-B3LYP算电荷转移效果差不多，但我个人不是很喜欢这个泛函，体感算柔性体系容易印刷虚频。

PS:我也是实验组自学计算哈哈

Stardust0831 发表于 2025-5-11 01:05
我也来交作业啦，wal佬用的是CAM-B3LYP，这是我用wB97XD计算出来的结果，与他的结果相近。

谢谢Stardust老师！ღ( ´･ᴗ･` )

wal 发表于 2025-5-10 23:37
你的第一个计算文件在第6步的时候梯度飞了，随后振子强度突然走低，感觉是B3LYP低估了某电荷转移态能量，随 ...

谢谢w老师，狠狠学习辽，我是实验组自学的计算，类似的这种有机小分子我一直是使用b3lyp/6-31G**，那激发任务您一般推荐是使用哪个呢？

我也来交作业啦，wal佬用的是CAM-B3LYP，这是我用wB97XD计算出来的结果，与他的结果相近。

Freeman 发表于 2025-5-10 17:02
输出文件显示，优化到最后，激发能只有0.0几电子伏特了，说明出现了S0/S1势能面交叉。这是几何优化振荡和SC ...

我和朋友讨论了一下。觉得不应该只看到激发能小，还应该注意到S1的振子强度非常之小，有ghost态的嫌疑。

Freeman 发表于 2025-5-10 17:02
输出文件显示，优化到最后，激发能只有0.0几电子伏特了，说明出现了S0/S1势能面交叉。这是几何优化振荡和SC ...

谢谢free老师，这个情况我还没见过，感谢free老师解惑，我也去学习一下CAS试试能不能解决这个问题

sobereva 发表于 2025-5-10 17:56
有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，导致信息零零碎碎， ...

不好意思sob老师，潜水至今第一次发，忘记规矩了，已更正，谨以为戒orz

wal 发表于 2025-5-10 21:08
这有点抽象，看起来不像是我熟悉的那种容易走圆锥交叉的分子

我也觉得。有可能是结构不合理而导致的问题，我这边正在尝试算他这个。有结果了会发。

这有点抽象，看起来不像是我熟悉的那种容易走圆锥交叉的分子

czjtry 发表于 2025-5-10 16:10
输出文件压缩

有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，导致信息零零碎碎，这点在置顶的新社员必读贴里明确说了。

输出文件显示，优化到最后，激发能只有0.0几电子伏特了，说明出现了S0/S1势能面交叉。这是几何优化振荡和SCF不收敛的主要原因。此时不能用TDDFT了。可以尝试CAS或MRSF-TDDFT，但是它们的使用细节颇有讲究，而且不能用高斯做。