大共轭平面自旋极化单重态计算求助

阿社社社 · 发表于 Post on 2023-4-25 20:57:34

各位老师好，请教一个共轭平面自旋极化单重态能量计算，结构如下图，我看了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去优化得到最后结构。感觉逻辑上没问题，我正在尝试，老师们觉得这个方法如何？

wzkchem5 · 发表于 Post on 2023-4-25 23:53:58

本帖最后由 wzkchem5 于 2023-4-25 17:46 编辑

这个体系通过画最优共振式，以及利用化学常识（alpha-氨基碳自由基和酚氧自由基的氧化还原电位）可知基态几乎不可能是自旋极化单重态。不是所有存在自旋极化单重态共振式的体系，都存在自旋极化单重态的稳定基态波函数。建议用TDDFT，以闭壳层单重态为参考态，把自旋极化单重态当作激发态来算，或者用MOM等方法做DeltaSCF计算也行

sobereva · 发表于 Post on 2023-4-26 06:55:03

那叫基于片段初猜计算，看下文，根本用不着你自己分成两块去计算
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 大共轭平面自旋极化单重态计算求助

浏览过的版块