计算化学公社

标题: 求助：使用OpenMolcas的GUGACI模块，算出来显然不合理的结果？ [打印本页]

---

作者
Author: 海棠    时间: 2023-7-16 16:02
标题: 求助：使用OpenMolcas的GUGACI模块，算出来显然不合理的结果？


这是输入文件 h2.input

1. &GATEWAY
   
2.   coord=h2.xyz
   
3.   basis=6-31g
   
4.   group=nosym
   
5. &SEWARD
   
6. &SCF
   
7.   charge=-2
   
8. &MOTRA
   
9. &GUGADRT
   
10.   electrons=4
    
11.   inactive=1
    
12.   active=3
    
13.   ciall
    
14. &GUGACI
    
15.   nrroots=4

复制代码

坐标文件 h2.xyz

1. 2
   
2. title
   
3. H 0.0 0.0 0.0
   
4. H 0.0 0.0 1.0

复制代码

SCF 模块的结果

1. ::    Total SCF energy                                 -0.3906595868
   
2.       One-electron energy                              -3.2904119394
   
3.       Two-electron energy                               2.3705751417
   
4.       Nuclear repulsion energy                          0.5291772109
   
5.       Kinetic energy (interpolated)                     1.6082712508
   
6.       Virial theorem                                    0.2429065287
   
7.       Total spin, S(S+1)                                0.0000000000
   
8.       Total spin, S                                     0.0000000000
   
9.       Max non-diagonal density matrix element           0.0000003295
   
10.       Max non-diagonal Fock matrix element              0.0000001593

复制代码

GUGACI 模块的结果中出现了

1. nroot      ci energy          dav energy        coef
   
2.    1         -3.614647078       -3.614690004   0.996155
   
3.    2         -2.977814407       -2.977819501   0.998547
   
4.    3         -2.782980266       -2.783085263   0.993457
   
5.    4         -2.456700695       -2.456784615   0.993578

复制代码

为什么 CI 算得的能量比 SCF 高呢？
另外，运行程序到 GUGACI 模块时，有报错：

1. double free or corruption (out)

复制代码

另附输出信息在附件。

---

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

Powered by Discuz! X3.3