本帖最后由 wal 于 2025-5-22 16:33 编辑 〇66 发表于 2025-5-22 16:06 1-3。没看到输出文件,我这里先假设9999对应步数上限,502对应SCF不收敛,你自行检查一下有没有不符合的 9999报错明显与SCF不收敛无关,也谈不上CI,大概率是初猜太差了,达到步数上限还没到驻点,要么就是震荡了 你这里标红的就060和280看着可能是接近CI引起的,别的不像,应当作为几何优化不收敛问题去解决 100这个SCF不收敛可以看看是SCF震荡还是啥别的原因 4.f很低也可能是电荷分离态/对称禁阻跃迁等,做一下激发分析看看 |
本帖最后由 wal 于 2025-5-15 22:12 编辑 〇66 发表于 2025-5-15 21:46 越接近CI越可能出现SCF不收敛。但单单试图解决SCF不收敛很可能是徒劳的,因为这个点收敛了,再往前走几个点真的走进CI就彻底寄了。 如果你想算顺反异构的能量变化趋势(其中包含CI的情况下),就得用能描述CI的方法,比如ORCA的SF-TDDFT或GAMESS的MRSF-TDDFT,或者直接使用CASSCF进行结构优化。 |
wal 发表于 2025-3-21 14:54 谢谢你的解答,我还有个问题,我扫描不出来的点有很多个,但是有的激发能很低,f=0.0000,但是有的激发能就还好,可能会达到f=0.0500左右。如果说f=0.0000是进入了CI,f=0.0500左右的点又作何解释呢?DFT如果不行,这种情况还有别的更好的办法能算出来他的顺反异构过程的能量变化趋势吗?非常感谢解答 |
本帖最后由 wal 于 2025-3-21 15:00 编辑 〇66 发表于 2025-3-21 12:48 这是个像“多少根头发算一把头发”一样比较唯心的问题,没有标准答案 倒不如说scf怎么算都不收敛/出现波函数不稳定的情况组合上激发能很小/为负才是进入CI范围更好的提示(对于S0-Sn的CI来说) |
| 参与人数Participants 1 | eV +5 | 收起 理由Reason |
|---|---|---|
|
| + 5 | 赞! |
wal 发表于 2025-3-21 12:23 好的 谢谢您了,我已经删除并重新压缩上传了。一般激发能为多少的时候 就说明进入CI区域了呢 |
本帖最后由 wal 于 2025-3-21 12:25 编辑 〇66 发表于 2025-3-21 12:18 建议删掉 sob老师看见这么长的楼会揍你的 下次把输出文件压缩上传就好明显看到最后激发能只有0.09ev,走到CI区域了,普通的TD-DFT已无法胜任你的任务 |
1000185497.jpg (441.21 KB, 下载次数 Times of downloads: 37)
|
本帖最后由 〇66 于 2025-3-21 12:47 编辑 激发能不是激发态和基态的差值吗(可能我理解有误吧),这是我的分子S1态下醛基在转到60°的时候的输出文件,可以帮忙看看哪里是激发能吗?实在麻烦您了
pht0-1_060.zip
(387.6 KB, 下载次数 Times of downloads: 0)
|
本帖最后由 wal 于 2025-3-21 10:55 编辑 〇66 发表于 2025-3-21 10:39 不用单独算基态啊 你激发态计算没有打印激发能吗? CI在论坛里搜索"圆锥交叉" 简单来讲是同自旋的两个势能面相遇了 是个比较麻烦的玩意 遇见CI的时候普通DFT可能会失效 |
wal 发表于 2025-3-20 23:36 我还没有计算基态这个的转动情况,激发能还没有,等我算算基态的转动情况再看看,基态只计算了最稳定的构型。可以麻烦问一下什么是CI吗 |
|
本帖最后由 wal 于 2025-3-20 23:38 编辑 看看不收敛之前激发能是多少 如果激发能很低甚至是负的说明碰见CI了 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-1-25 18:28 , Processed in 0.287499 second(s), 26 queries , Gzip On.