计算化学公社

标题: 大共轭平面自旋极化单重态计算求助 [打印本页]

作者
Author:
阿社社社    时间: 2023-4-25 20:57
标题: 大共轭平面自旋极化单重态计算求助
各位老师好,请教一个共轭平面自旋极化单重态能量计算,结构如下图,我看了sob老师帖子说可以分段计算波函数,最后在整体计算。我是把这个分子从中间得双键那断开,考虑到是双键,我个人是为是两各片段各带三个单电子,自旋多重度分别为4和-4。但是计算结果提示No NMR shielding tensors so no spin-rotation constants. The archive entry for this job was punched. 这部得到的chk文件在用于计算整体,一秒钟就算完了,也导不出波函数文件。是不是这种结构使用分段计算的设置有不同?
后来我看到一篇文献计算方法(也是一个共轭平面),他是先用限制b3lyp优化结构(这步得到的是闭壳),然后用stable=opt计算波函数稳定性,得到一个能量较低的解(开壳),再用ub3lyp guess=read去优化得到最后结构。感觉逻辑上没问题,我正在尝试,老师们觉得这个方法如何?





作者
Author:
wzkchem5    时间: 2023-4-25 23:53
本帖最后由 wzkchem5 于 2023-4-25 17:46 编辑

这个体系通过画最优共振式,以及利用化学常识(alpha-氨基碳自由基和酚氧自由基的氧化还原电位)可知基态几乎不可能是自旋极化单重态。不是所有存在自旋极化单重态共振式的体系,都存在自旋极化单重态的稳定基态波函数。建议用TDDFT,以闭壳层单重态为参考态,把自旋极化单重态当作激发态来算,或者用MOM等方法做DeltaSCF计算也行
作者
Author:
sobereva    时间: 2023-4-26 06:55
那叫基于片段初猜计算,看下文,根本用不着你自己分成两块去计算
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82

No NMR shielding tensors so no spin-rotation constants. The archive entry for this job was punched根本不是报错




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3