如何在Gaussian中优化结构后得到指定电子态后计算NMR？

lukas

各位好，我想请教一下，如果优化得到了对称破缺态，现在希望计算这个态的分子中芳香环nics(1)zz。仅以其out结果中的几何结构,在这样的关键词#pb3lyp/6-31+G* NMR geom=connectivity guess=mix 下计算，但发现最后NMR任务的fchk文件的轨道并不是对称破缺态，而是闭壳层单重态。请问gaussian中怎么计算得到这样的对称破缺态的NMR结果？

sobereva

lukas 发表于 2024-12-2 22:25
老师好，请问哪一种对称破缺态是合理的呢？比如计算得到很多对称破缺态，有很多种自旋分布在不同的原子上 ...

能量最低的

顺带提一句，Gaussian对开壳层做NICS、结合Multiwfn做ICSS考察芳香性是可以的，参考我下面这篇，里面通过磁属性对16碳环不同自旋态对比了芳香性差异
Yang Wu, Xiufen Yan, Zeyu Liu*, Tian Lu,* Double Aromaticity of Franck-Condon Excited States of Cyclo[16]carbon, ChemRxiv (2023) DOI: 10.26434/chemrxiv-2023-mtx00

lukas

sobereva 发表于 2024-11-25 23:03
Gaussian不支持CASSCF下算NMR
做CASSCF必须根据实际情况选择合适的活性空间，不是光写个CAS(2,2)就了事的
...

老师好，请问哪一种对称破缺态是合理的呢？比如计算得到很多对称破缺态，有很多种自旋分布在不同的原子上的对称破缺态，请问是能量最低的态最合理还是和三重态自旋分布是一致的态最合理呢？

zjxitcc

lukas 发表于 2024-11-25 22:04
谢谢老师。但我不清楚是说高斯中计算nmr在dft产生的轨道下不严谨还是说高斯的nmr计算方法不严谨，如果使c ...

不能，十万八千里。建议使用MOKIT联用多个量化程序自动进程CASSCF NMR计算，见https://jeanwsr.gitlab.io/mokit- ... ml#54-icss-and-nics

sobereva

Gaussian不支持CASSCF下算NMR
做CASSCF必须根据实际情况选择合适的活性空间，不是光写个CAS(2,2)就了事的

先确保得到了对称破缺态，然后guess=read读取其波函数当计算NMR任务的初猜。guess=mix没那么万能，仔细看下文
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82

lukas

wzkchem5 发表于 2024-11-25 13:19
高斯应该算不了unrestricted体系的核磁，就算能算，结果也不可靠，一是unrestricted体系的核磁有一些额外项 ...

谢谢老师。但我不清楚是说高斯中计算nmr在dft产生的轨道下不严谨还是说高斯的nmr计算方法不严谨，如果使casscf 是不是先hf得到初猜，然后 #P CAS(2,2)/6-31G* guess=read nmr ，是否能实现您说的算法。

wzkchem5

高斯应该算不了unrestricted体系的核磁，就算能算，结果也不可靠，一是unrestricted体系的核磁有一些额外项的贡献，二是对称破缺态不满足时间反演不变性。
这种体系建议用CASSCF算核磁，有可能的话用CASPT2更好，不过我不知道有没有程序支持