请选择 进入手机版 | 继续访问电脑版
第13届北京科音初级量子化学培训班将于10月5~8日于北京举办,请点击此链接查看详情。这是新人一次性正确、完整学习量子化学计算的最好、最快机会,能少走无数弯路,欢迎参加并相互转告!(已报满)

计算化学公社

 找回密码
 现在注册!
查看: 145|回复: 7

[Molpro] molpro用CASSCF算NO势能曲线求助

[复制链接]

37

帖子

0

威望

95

eV
积分
132

Level 2 能力者

发表于 2020-9-1 11:42:01 | 显示全部楼层 |阅读模式
本人Molpro初学者:用casscf计算NO,下面是我编辑的molpro输入文件!结果总是说全局处理器错误,而且也出不来图片,请大家帮忙看一下,谢谢!

molpro 算NO 势能曲线

molpro 算NO 势能曲线

18

帖子

0

威望

406

eV
积分
424

Level 3 能力者

发表于 2020-9-1 14:05:20 | 显示全部楼层
建议将out文件里报错信息贴出来

37

帖子

0

威望

95

eV
积分
132

Level 2 能力者

 楼主| 发表于 2020-9-1 15:56:55 | 显示全部楼层
绫波レイ 发表于 2020-9-1 14:05
建议将out文件里报错信息贴出来

你好,下面是输出文件,它还是说全局错误
Primary working directories    : /tmp
Secondary working directories  : /tmp
Wavefunction directory         : /public1/home/sc31181/wfu/
Main file repository           : /tmp/

ARCHNAME  : Linux/x86_64
FC        : /opt/intel/compilerpro-12.0.2.137/bin/intel64/ifort
FCVERSION : 12.0.2
BLASLIB   :
id        : tech

Nodes          nprocs
c3801.para.cstc   24
ga_uses_ma=false, calling ma_init with nominal heap.
GA-space will be limited to   8.0 MW (determined by -G option)

Using customized tuning parameters: mindgm=4; mindgv=4; mindgc=1; mindgr=1; noblas=0; mincuda=1000; minvec=7
default implementation of scratch files=sf  

***,NO merge
r=[1,2,3,4,5] ang
  
i=1
geometry={N
           O,1,r(i)}
  
basis=av5z
hf
  
do i=1,#r
{multi
closed,2,0,0,0
occ,   6,2,2,0
wf;state,2
canonical,ci}
  
{rs2,MIX=2,INIT
wf;state,1,1}
  
e1_caspt2(i)=energy
  
{rs2,MIX=2
wf;state,1,2}
  
e2_caspt2(i)=energy
e1_mscaspt2(i)=msenergy(1)
e2_mscaspt2(i)=msenergy(2)
enddo
  
{table,r,e1_caspt2,e2_caspt2,e1_mscaspt2,e2_mscaspt2
title,SS-SR-CASPT2 for NO
plot,file='no_sr_mscaspt2.plot'
}

Variables initialized (660), CPU time= 0.01 sec
Commands  initialized (462), CPU time= 0.03 sec, 486 directives.
Default parameters read. Elapsed time= 1.13 sec

Checking input...
Passed
1


                                         ***  PROGRAM SYSTEM MOLPRO  ***
                         Copyright, University College Cardiff Consultants Limited, 2008

                                    Version 2010.1 linked 15 Sep 2011 12:01:52


**********************************************************************************************************************************
LABEL *   NO merge                                                                     
Linux-3.10.0-1127.18.2.el7.x86_64/c3801.para.cstc(x86_64) 64 bit mpp version            DATE: 01-Sep-20          TIME: 15:52:01  
**********************************************************************************************************************************

SHA1:             c82fc9886d3ef081defdb91dc08d0bb97a8ad3e0
**********************************************************************************************************************************
SETTING R(1)           =         1.00000000  ANG                             
SETTING R(2)           =         2.00000000  ANG                             
SETTING R(3)           =         3.00000000  ANG                             
SETTING R(4)           =         4.00000000  ANG                             
SETTING R(5)           =         5.00000000  ANG                             
SETTING I              =         1.00000000                                 

Variable memory set to    8000000 words,  buffer space   230000 words

SETTING BASIS          =    AV5Z


Recomputing integrals since basis changed


Using spherical harmonics

Library entry N      S aug-cc-pV5Z          selected for orbital group  1
Library entry N      P aug-cc-pV5Z          selected for orbital group  1
Library entry N      D aug-cc-pV5Z          selected for orbital group  1
Library entry N      F aug-cc-pV5Z          selected for orbital group  1
Library entry N      G aug-cc-pV5Z          selected for orbital group  1
Library entry N      H aug-cc-pV5Z          selected for orbital group  1
Library entry O      S aug-cc-pV5Z          selected for orbital group  2
Library entry O      P aug-cc-pV5Z          selected for orbital group  2
Library entry O      D aug-cc-pV5Z          selected for orbital group  2
Library entry O      F aug-cc-pV5Z          selected for orbital group  2
Library entry O      G aug-cc-pV5Z          selected for orbital group  2
Library entry O      H aug-cc-pV5Z          selected for orbital group  2

1PROGRAM * SEWARD (Integral evaluation for generally contracted gaussian basis sets)     Author: Roland Lindh, 1990

Geometry written to block  1 of record 700

Orientation using atomic masses  
Molecule type: Linear
Symmetry elements: X,Y
Rotational constants:      67.6686037     67.6686037      0.0000000 GHz  (calculated with average atomic masses)

Point group  C2v



ATOMIC COORDINATES

NR  ATOM    CHARGE       X              Y              Z

   1  N       7.00    0.000000000    0.000000000   -1.007611261
   2  O       8.00    0.000000000    0.000000000    0.882114870

Bond lengths in Bohr (Angstrom)

1-2  1.889726131
     ( 1.000000000)

NUCLEAR CHARGE:                   15
NUMBER OF PRIMITIVE AOS:         398
NUMBER OF SYMMETRY AOS:          288
NUMBER OF CONTRACTIONS:          254   (  92A1  +  62B1  +  62B2  +  38A2  )
NUMBER OF CORE ORBITALS:           2   (   2A1  +   0B1  +   0B2  +   0A2  )
NUMBER OF VALENCE ORBITALS:        8   (   4A1  +   2B1  +   2B2  +   0A2  )


NUCLEAR REPULSION ENERGY   29.63392370

EXTRA SYMMETRY OF AOS IN SYMMETRY 1:   1 1 1 1 1 1 1 1 1 1   1 1 1 1 2 1 2 1 2 1   2 1 2 1 2 1 2 1 2 1   2 1 3 2 1 3 2 1 3 2
                                        2 3 1 2 3 1 1 1 1 1   1 1 1 1 1 1 1 1 1 1   2 1 2 1 2 1 2 1 2 1   2 1 2 1 2 1 2 1 3 2
                                        1 3 2 1 3 2 2 3 1 2   3 1
EXTRA SYMMETRY OF AOS IN SYMMETRY 2:   1 1 1 1 1 1 1 1 1 1   1 1 2 1 2 1 2 1 2 1   2 1 2 1 2 1 2 3 1 2   3 1 1 1 1 1 1 1 1 1
                                        1 1 1 2 1 2 1 2 1 2   1 2 1 2 1 2 1 2 3 1   2 3
EXTRA SYMMETRY OF AOS IN SYMMETRY 3:   1 1 1 1 1 1 1 1 1 1   1 1 2 1 2 1 2 1 2 1   2 1 2 1 2 1 2 3 1 2   3 1 1 1 1 1 1 1 1 1
                                        1 1 1 2 1 2 1 2 1 2   1 2 1 2 1 2 1 2 3 1   2 3
EXTRA SYMMETRY OF AOS IN SYMMETRY 4:   1 1 1 1 1 1 1 1 1 1   2 1 2 1 2 2 1 2 1 1   1 1 1 1 1 1 1 1 1 2   1 2 1 2 2 1 2 1

Eigenvalues of metric

         1 0.831E-05 0.630E-04 0.928E-04 0.226E-03 0.245E-03 0.546E-03 0.718E-03 0.834E-03
         2 0.657E-04 0.529E-03 0.113E-02 0.186E-02 0.419E-02 0.555E-02 0.106E-01 0.142E-01
         3 0.657E-04 0.529E-03 0.113E-02 0.186E-02 0.419E-02 0.555E-02 0.106E-01 0.142E-01
         4 0.834E-03 0.649E-02 0.248E-01 0.282E-01 0.648E-01 0.882E-01 0.943E-01 0.159E+00


Contracted 2-electron integrals neglected if value below      1.0D-12
###############################################################################
###############################################################################
###                                                                         ###
###                                                                         ###
###     Allocation failed for Work3                                         ###
###                                                                         ###
###                                                                         ###
###############################################################################
###############################################################################
                 37892                 66072                 66072
                 30371                     0                     0
                     1                     1                     1
                     1                     4                     4
                     1                     1
xquit called with rc= 128

GLOBAL ERROR fehler on processor   0                                         

37

帖子

0

威望

95

eV
积分
132

Level 2 能力者

 楼主| 发表于 2020-9-2 16:57:31 | 显示全部楼层
顶一下,大神可以帮我看一下到底是哪里出错了

964

帖子

0

威望

3381

eV
积分
4345

Level 6 (一方通行)

发表于 2020-9-2 17:00:39 | 显示全部楼层
wulalalana 发表于 2020-9-2 16:57
顶一下,大神可以帮我看一下到底是哪里出错了

提示 Allocation failed for Work3,大概率是内存不足。你设置内存再算
欢迎关注“量子化学”公众号,也欢迎投稿。

37

帖子

0

威望

95

eV
积分
132

Level 2 能力者

 楼主| 发表于 2020-9-3 10:53:15 | 显示全部楼层
zjxitcc 发表于 2020-9-2 17:00
提示 Allocation failed for Work3,大概率是内存不足。你设置内存再算

您好,我做了一些修改之后,出现了下面的问题,它这里是说我的 wf 设置的不对吗,我已经对wf 做了多次修改,仍然是下面的这种情况,麻烦您看一下,谢谢

Inconsistent input for SPIN (MS2):

NELEC= 15  MS2=0  CHARGE= 0
? Error
? Wrong spin
? The problem occurs in get_wf

GLOBAL ERROR fehler on processor   0      

964

帖子

0

威望

3381

eV
积分
4345

Level 6 (一方通行)

发表于 2020-9-3 11:03:05 | 显示全部楼层
wulalalana 发表于 2020-9-3 10:53
您好,我做了一些修改之后,出现了下面的问题,它这里是说我的 wf 设置的不对吗,我已经对wf 做了多次修 ...

你这MS2是何含义?自旋多重度?alpha-beta电子数?NO的自旋多重度是2.
欢迎关注“量子化学”公众号,也欢迎投稿。

10

帖子

0

威望

43

eV
积分
53

Level 2 能力者

发表于 4 天前 | 显示全部楼层
你的问题最终解决了么?
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2020-9-20 21:02 , Processed in 0.154878 second(s), 28 queries .

快速回复 返回顶部 返回列表