计算自旋极化单重态划分片段求助

adreamy

我想计算的分子是具有自由基阴离子的苝酰亚胺，构建这样的分子式是否合理（如下图）？按照谈谈片段组合波函数与自旋极化单重态 - 思想家公社的门口：量子化学·分子模拟·二次元 (sobereva.com)帖子中利用guess=mix stable 以及nosymm关键词计算自旋极化单重态时，计算结束提示波函数是稳定的，但是用Multiwfn计算自旋密度提示Summing up positive value and multiply differential element：2.15E-5,得不到计算自旋极化单重态。计算时按照中性体系带自由基计算的，是否需要对整个分子设置上带两个负电荷？另外，如果按照划分片段计算自旋极化单重态，应该如何划分片段？能设置三个片段马？如每个带自由基的氧原子为一个片段，剩下的为另一个片段。 应该是我分子式弄错了，打扰大家了，我不知道怎么删除帖子，请求sob老师看到的时候帮忙删下贴子，谢谢。

adreamy

应该是我画错分子式了，打扰大家了。但是我不知道怎么删帖。如果有知道的可以留言，我把帖子删除了

zjxitcc

你说“自由基阴离子”，那你这阴离子符号"-"画在哪里？自由基符号"·"画在哪里？你这看着是个中性分子。合理自由基阴离子示例：(O2)·-

有机体系不需要使用片段组合波函数作为初猜，如果你nosymm guess=mix stable发现波函数是稳定的、且<S**2>=0，那说明在当前计算级别下，它真的就是闭壳层的，不用再试了。

DFT泛函种类无数。对于有极强双自由基特征、四自由基特征的分子，你用大多数泛函算，总会算出<S**2>大于0；对于很弱双自由基特征的分子，你用大多数泛函算，总会算出<S**2>=0；处于中间阶段的分子，显然有的泛函会给出闭壳层的结果，有的泛函会给出开壳层单重态的结果。

adreamy

zjxitcc 发表于 2023-8-26 14:36
你说“自由基阴离子”，那你这阴离子符号"-"画在哪里？自由基符号"·"画在哪里？你这看着是个中性分子。合 ...

谢谢老师，我之前把分子式画错了，打扰老师了，我在自己试试。