计算化学公社

标题: 三重态结构优化一直不收敛问题 [打印本页]

作者
Author: liujinwang 时间: 2025-3-11 10:01
标题: 三重态结构优化一直不收敛问题
由于第一次对这种多原子结构进行结构优化，所以根据文献选择了PBE1PBE泛函，6-31G*(d,p)/SDD基组，但是一开始跑的时候一直不收敛，后面根据sob社长关于SCF不收敛问题选择了受力点最小的一个结构进行优化，并加入了gdiis关键词，也尝试过加入(gdiis,maxstep=3,notrust)关键词，但是也一直不收敛。后面又尝试过用b3lyp/6-311G(d,p)进行小基组做一个初猜的结构优化，但是小基组也不收敛，一直震荡，所以发帖想问一下现在这个还要继续等他跑完么，以及是否有其他更好的方法解决这个问题。

作者
Author: linjunyi 时间: 2025-3-11 10:12
试试calcfc，算力够的话试试calcall

作者
Author: sobereva 时间: 2025-3-11 11:53
没有叫6-31G*(d,p)的东西
加上nosymm，便于监控哪里震荡
用calcfc（recalc=10或15更好）
当前溶剂极性非常低，去掉溶剂模型有可能收敛容易一点，耗时更低而且也不会影响结果合理性
另外，当前体系当成三重态算意义不明

作者
Author: liujinwang 时间: 2025-3-11 15:28

sobereva 发表于 2025-3-11 11:53
没有叫6-31G*(d,p)的东西
加上nosymm，便于监控哪里震荡
用calcfc（recalc=10或15更好）

好的，谢谢社长，当前体系算三重态想优化得到磷光发射能

作者
Author: liujinwang 时间: 2025-3-18 10:35

sobereva 发表于 2025-3-11 11:53
没有叫6-31G*(d,p)的东西
加上nosymm，便于监控哪里震荡
用calcfc（recalc=10或15更好）

果然一下就收敛了，但是发现输出的单点能HOMO-LUMO轨道只有一半，是不是加入POP=full就可以了。

作者
Author: sobereva 时间: 2025-3-19 10:05

liujinwang 发表于 2025-3-18 10:35
果然一下就收敛了，但是发现输出的单点能HOMO-LUMO轨道只有一半，是不是加入POP=full就可以了。

“只有一半”语义不明
搞清楚什么叫限制性闭壳层和非限制性开壳层，这是最基本的量化概念。跟pop=full没半毛关系

欢迎光临计算化学公社 (http://bbs.keinsci.com/)