计算化学公社

标题: 求助Gaussian优化NO分子输出文件有误 [打印本页]

作者
Author: 原子分子物理 时间: 2020-10-20 09:38
标题: 求助Gaussian优化NO分子输出文件有误
大家好，最近在做与NO分子有关的计算，但是在第一步优化的时候就出了问题，文件可以完整输出，但是里面的输出信息我感觉有点奇怪。优化用的泛函和基组是ub3lyp/def2svp，输出文件有问题部分如下：
Population analysis using the SCF density.

**********************************************************************

Orbital symmetries:
Alpha Orbitals:
   Occupied  (SG) (SG) (SG) (SG) (PI) (PI) (SG) (PI)
   Virtual (PI) (SG) (SG) (PI) (PI) (SG) (SG) (PI) (PI) (SG)
               (PI) (PI) (DLTA) (DLTA) (DLTA) (DLTA) (SG) (PI)
               (PI) (SG)
Beta  Orbitals:
   Occupied  (SG) (SG) (SG) (SG) (PI) (PI) (SG)
   Virtual (PI) (PI) (SG) (SG) (PI) (PI) (SG) (SG) (PI) (PI)
               (SG) (PI) (PI) (DLTA) (DLTA) (DLTA) (DLTA) (SG)
               (PI) (PI) (SG)
Unable to determine electronic state:  partially filled degenerate orbitals.
Alpha  occ. eigenvalues --  -19.27276 -14.46110  -1.23825  -0.67261  -0.53920
Alpha  occ. eigenvalues -- -0.51235  -0.48293  -0.21537
Alpha virt. eigenvalues -- -0.10422 0.33586 0.52321 0.65759 0.67142

我的问题是文中标红部分怎么解决。PS:我的目的是做CASPT2计算.
希望知道怎么解决的大神能告知一下，感谢

作者
Author: wzkchem5 时间: 2020-10-20 09:52
这是正常的，因为NO的基态是简并的，DFT不能描述简并基态，会把体系的点群对称性由C∞v破缺到C2v，所以程序自然不知道在C∞v下所属的不可约表示是什么。
换用(state-averaged) CASSCF、CASPT2等可以描述简并基态的方法后，这个问题就解决了

作者
Author: biogon 时间: 2020-10-20 09:54
本帖最后由 biogon 于 2020-10-20 10:01 编辑

NO基态简并了

作者
Author: 原子分子物理 时间: 2020-10-20 10:01

wzkchem5 发表于 2020-10-20 09:52
这是正常的，因为NO的基态是简并的，DFT不能描述简并基态，会把体系的点群对称性由C∞v破缺到C2v，所以程序 ...

是的。基态是X2pi,是简并的，“换用(state-averaged) CASSCF、CASPT2等可以描述简并基态的方法后，这个问题就解决了”这句我不是很懂，能麻烦您具体说说吗

作者
Author: 原子分子物理 时间: 2020-10-20 10:03

biogon 发表于 2020-10-20 09:54
NO基态简并了

那产生的chk文件能用吗

作者
Author: zjxitcc 时间: 2020-10-20 10:05
你打算如何做CASPT2计算呢？高斯不支持CASPT2，你要到OpenMolcas, Molpro之类的软件里去算CASPT2的。先拿高斯算单点练练手？

作者
Author: biogon 时间: 2020-10-20 10:10

原子分子物理发表于 2020-10-20 10:03
那产生的chk文件能用吗

先看轨道图形跑出来的是什么东西

作者
Author: 原子分子物理 时间: 2020-10-20 10:11

zjxitcc 发表于 2020-10-20 10:05
你打算如何做CASPT2计算呢？高斯不支持CASPT2，你要到OpenMolcas, Molpro之类的软件里去算CASPT2的。先拿高 ...

是的，我想先用高斯练练手，之后再用OpenMolcas算

作者
Author: zjxitcc 时间: 2020-10-20 10:19
条条大路通罗马，做法有很多种。我不用对称性（个人偏好），我用MOKIT（https://gitlab.com/jxzou/mokit）里的AutoMR程序自动做CASPT2计算的输入文件举个例子

%mem=2GB
%nprocshared=1
#p CASPT2/def2TZVP
mokit{}
0 2
N 0.0 0.0 0.0
O 0.0 0.0 1.136

复制代码

运行
automr NO_def2TZVP.gjf >& NO_def2TZVP.out &

整个过程就1 min时间。最终结果是基于CASSCF(3,3)的CASPT2，OpenMolcas的输入输出文件都有。可以在此基础上自己写OpenMolcas的关键词再算激发态等其他问题。当然，这个做法没法给出对称性信息，优点是可以迅速算、不需人为干预。

作者
Author: 原子分子物理 时间: 2020-10-20 10:36

zjxitcc 发表于 2020-10-20 10:19
条条大路通罗马，做法有很多种。我不用对称性（个人偏好），我用MOKIT（https://gitlab.com/jxzou/mokit） ...

好的，谢谢老师，我研究一下这个程序。

作者
Author: wzkchem5 时间: 2020-10-20 12:20

原子分子物理发表于 2020-10-20 10:01
是的。基态是X2pi,是简并的，“换用(state-averaged) CASSCF、CASPT2等可以描述简并基态的方法后，这个问 ...

就是比如你做state-averaged CASSCF，你是优化的两个简并基态的平均能量，所以两个态总是能保持简并，可以严格保持C∞v对称性。
但是DFT只能优化一个态的能量，所以会把两个简并态里的一个态的能量选择性拉低，破坏简并性

作者
Author: 原子分子物理 时间: 2020-10-20 14:57

wzkchem5 发表于 2020-10-20 12:20
就是比如你做state-averaged CASSCF，你是优化的两个简并基态的平均能量，所以两个态总是能保持简并，可 ...

理解了，感谢大神。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)