计算化学公社

标题: 求助寡聚物的单-三线态能级差与双自由基指数的计算 [打印本页]

作者
Author: vverawei 时间: 2023-11-9 10:18
标题: 求助寡聚物的单-三线态能级差与双自由基指数的计算
各位老师好，我目前在计算寡聚物的单线态-三线态能级差和双自由基指数y0值，想请问各位老师，我计算这两个的泛函用的不一样可以吗？我计算ΔEST的时候用的是UPBE0/6-31G**，计算y0值的时候使用的UHF/6-31G** guess=mix，如果计算y值也使用UPBE0/6-31G** guess=mix 的话计算结果会差很大。还有就是想请问各位老师关于三线态电子自旋密度分布与y值的关系，y值越接近于1则双自由基特征越明显，是否在三线态自旋密度分布上表示为密度逐渐朝两端分布呢？

下面这个图是（D-A）n, n=5的时候我计算的自旋密度分布，但此时如果使用UHF/6-31g(d,p) guess=mix SCF=conver=4 来计算的y值却为0.99，而使用UPBE0/6-31g(d,p)的话y值计算结果为0.44，较为符合，但是我们实验上对它的预测是希望它能向两端分布更合理，图2是随着链长的增长的三线态自旋密度分布变化的趋势，想请教一下各位老师n=5时出现的一个异常的沿链分布是因为什么原因呢？

作者
Author: zjxitcc 时间: 2023-11-9 10:30
本帖最后由 zjxitcc 于 2023-11-9 10:33 编辑

SCF=conver=4结果不可信，基于此结果的分析不可信。计算不是随便算。另外，此问题不属于VASP分区。

作者
Author: vverawei 时间: 2023-11-9 10:43

zjxitcc 发表于 2023-11-9 10:30
SCF=conver=4结果不可信，基于此结果的分析不可信。计算不是随便算。另外，此问题不属于VASP分区。

好的谢谢老师，我改了分区。然后关于SCF=conver=4，我是因为用SCF=conver=6和SCF(novaracc,noincfock) int=acc2e=12都无法收敛以后才用的这个，我目前还在用SCF=vshift=300去算，请问老师关于y值和ΔEST的计算用不同泛函是可以的吗？

作者
Author: zjxitcc 时间: 2023-11-9 11:03

vverawei 发表于 2023-11-9 10:43
好的谢谢老师，我改了分区。然后关于SCF=conver=4，我是因为用SCF=conver=6和SCF(novaracc,noincfock) in ...

这个问题也不属于Quantum ESPRESSO分区。
SCF=conver=6通常只能作为中间计算结果使用，供收敛到正常精度。SCF不收敛是个人计算技术问题，必可收敛，而且双自由基体系一定要检验波函数稳定性，这两个你都没做到，因此还不能讨论结果如何如何。

y值和ΔEST的计算可以用不同泛函，不过HF一般不称作泛函。如果你用UHF/6-31G**，我建议你看看y1和y2，不要只盯着y0，万一你的体系有明显四自由基特征，不能误称双自由基特征。

作者
Author: vverawei 时间: 2023-11-9 11:35

zjxitcc 发表于 2023-11-9 11:03
这个问题也不属于Quantum ESPRESSO分区。
SCF=conver=6通常只能作为中间计算结果使用，供收敛到正常精度 ...

好的谢谢老师，我再去检查一下

作者
Author: sobereva 时间: 2023-11-9 17:44
量子化学问题别发到第一性原理版块，这次给你移动了，以后删帖处理

如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “单三线态能级与双自由基指数的计算 ” 改了，以后务必注意

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