小橙子cc 发表于 2026-3-24 15:33
感谢老师回复！非常不好意思，由于新手入门，我对您回复的仍有一些疑问：
1、我阅读了您举得例子，您的体 ...

1. 我在那个帖子里展示的计算步骤，每一步都有nosymm，有的步骤有stable=opt，都是写好的。你说你要加上nosymm stable=opt，我不理解。如果我没写这两个关键词，你说加上，这倒可以理解。既然我写了，你说要加上，那我就不知道你是什么意思了。比如一个关键词写两遍肯定是错误的。

2. 对于non-singlet计算，你写B3LYP，Gaussian默认使用UB3LYP。对于单重态计算，如果你一开始就想进行开壳层单重态计算，那么需要显式写明UB3LYP guess=mix nosymm；如果你写B3LYP，这会导致一开始做的是闭壳层计算，实际计算与用户期望不符。

3. 有机化学知识；计算经验。

wzkchem5 发表于 2026-3-20 16:56
除楼上说的之外，利用Hund规则、共振式等，可以排除掉一部分自旋多重度。例如你这个体系，画共振式容易发现 ...

感谢老师回复！

zjxitcc 发表于 2026-3-20 12:53
思路是对的，但具体操作仍需小心。例如对于单重态的单点能计算，这类体系经常是开壳层单重态，不仅需要gues ...

感谢老师回复！非常不好意思，由于新手入门，我对您回复的仍有一些疑问：
1、我阅读了您举得例子，您的体系是开壳层单重态。请问当我不知道体系是单重态还是三重态时，算开壳层单重态我按照您所说的的方法进行计算；对于三重态，我按照您说的加上nosymm stable=opt计算，然后比较哪个是最低的能量吗？
2、当我写输入文件的计算方法时，我往往只写B3LYP/6-311G(d,p)，没有在输入文件中使用UB3LYP或是RB3LYP，这个会对计算结果产生影响吗？
3、同问题2，我没太理解您的这句话：“对于单重态的单点能计算，这类体系经常是开壳层单重态”我该怎么判断到底是开壳层单重态还是闭壳层单重态？

除楼上说的之外，利用Hund规则、共振式等，可以排除掉一部分自旋多重度。例如你这个体系，画共振式容易发现pi电子没有自由基，上面有两个sigma自由基，因此自旋多重度只可能是1或3，且如果是1那么必然是开壳层单重态。因此根本不需要试自旋多重度为5的情况

思路是对的，但具体操作仍需小心。例如对于单重态的单点能计算，这类体系经常是开壳层单重态，不仅需要guess=mix构建对称破缺初猜，还需要nosymm stable=opt检验波函数稳定性，确保收敛到稳定的波函数。如果用RKS闭壳层单重态随便算，很容易踩坑。对于non-singlet的计算，虽然不需要构建对称破缺初猜，但仍建议nosymm stable=opt检验波函数稳定性。

我在论坛上回答过一个三角烯体系的计算步骤（这个分子是真正无歧义的开壳层单重态体系之一），可以参考
http://bbs.keinsci.com/forum.php ... 339957&fromuid=2632