24小时热门版块排行榜    

查看: 1683  |  回复: 5

小猪女魔头

铜虫 (初入文坛)

[求助] 第一次用MOLCAS算H2O2在激发态时的能量出现错误已有1人参与

各位虫友,我是新手,目前正在用新学的软件进行H2O2的固定结构优化和激发态的能量计算,但输出文件中出现错误,请各位高手帮助我一下,我金币不多,请各位谅解,谢谢!
以下是我的输入文件和输出文件

>>> DO WHILE <<<
  &SEWARD &END
    Title
    OO11OH15-CASPT2
    AMFI
    Basis set
    O.ANO-L...5s4p3d2f.
    O1      -0.502120000     -0.206942000      0.058689000   /Angstrom
    O2       0.551150000      0.045617000     -0.133239000   /Angstrom
    End of basis
    Basis set
    H.ANO-L...3s2p1d.
    H1      -1.434006000      0.960071000     -0.081556000   /Angstrom
    H2       1.041768000      0.330529000      0.677958000   /Angstrom
    End of basis
    End of Input
  &SCF
  &RASSCF &END
    spin
      1
    NACTEL
      14 0 0
    INACTIVE
      2
    RAS2
      10
    CIRoot
      2 2 1
    LUMORB
    Iterations
      200,200
    End of Input
>>COPY $Project.JobIph $Project.JobOld
  &CASPT2 &END
    Frozen
      2
    Multistate
      2 1 2
    IPEA
      0.0
    imag
      0.05
    MAXIterations
      100
    RLXRoot
      1
    End of Input
  &ALASKA &END
    End of Input
  &SLAPAF
    constraints
    B1 = Bond O1 O2
    B2 = Bond O1 H1
    Value
    B1 = 1.1  Angstrom
    B2 = 1.5  Angstrom
    End of Constraints
    Iterations
      100
    Thrs
      1.0D-05 3.0D-03
    End of Input
>>> EndDo <<<

  &RASSCF &END
    spin
      1
    NACTEL
      14 0 0
    INACTIVE
      2
    RAS2
      10
    CIRoot
      4 4 1
    LUMORB
    Iterations
      200,200
    End of Input
>>COPY $Project.JobIph $Project.JobOld
    !ln -sf $Project.JobMix001 JOBMIX
  &CASPT2 &END
    Frozen
      2
    Multistate
      4 1 2 3 4
    IPEA
      0.0
    imag
      0.05
    MAXIterations
      100
    End of Input
    !rm JOBMIX
    !rm JOBIPH
  &RASSCF &END
    spin
      3
    NACTEL
      14 0 0
    INACTIVE
      2
    RAS2
      10
    CIRoot
      4 4 1
    LUMORB
    Iterations
      200,200
    End of Input
>>COPY $Project.JobIph $Project.JobOld
    !ln -sf $Project.JobMix002 JOBMIX
  &CASPT2 &END
    Frozen
      2
    Multistate
      4 1 2 3 4
    IPEA
      0.0
    imag
      0.05
    MAXIterations
      100
    End of Input
    !rm JOBMIX
    !rm JOBIPH
    !ln -sf $Project.JobMix001 JOB001
    !ln -sf $Project.JobMix002 JOB002

  &RASSI &END
    Nr Of jobIph
      2 4 4
      1 2 3 4
      1 2 3 4
    spinorbit
    EJOB
    End of Input


输出文件

                          RASSCF iterations: Energy and convergence statistics
                          ----------------------------------------------------

      Iter CI   SX   CI       RASSCF       Energy    max ROT   max BLB   max BLB  Level Ln srch  Step   QN    Time(min)
          iter iter root      energy       change     param    element    value    shift minimum  type update     CPU
David: nBasVec less than lRoots
nBasvec, lRoots =                      3                     4
CIREstart was used. Check the number of roots in the previous calculation
--- Stop Module:  rasscf at Fri Apr  4 15:52:18 2014 /rc= _INTERNAL_ERROR_ ---
--- Module rasscf spent 43 seconds

         Some internal inconsistency of the code was detected

Non-zero return code - check program input/output
--- Stop Module:  auto at Fri Apr  4 15:52:20 2014 /rc= _INTERNAL_ERROR_ ---
--- Module auto spent 44 minutes and 37 seconds
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiangning198511

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
小猪女魔头: 金币+5, ★★★很有帮助 2014-04-05 20:10:06
先将自带的例子测试一下,成功后说明你的程序没问题。
其次,将任务分解,先将CASSCF 跑通,再做CASPT2
PS : MOLCAS的输入比较麻烦,建议使用MOLPRO。
2楼2014-04-05 15:20:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小猪女魔头

铜虫 (初入文坛)

引用回帖:
2楼: Originally posted by jiangning198511 at 2014-04-05 15:20:36
先将自带的例子测试一下,成功后说明你的程序没问题。
其次,将任务分解,先将CASSCF 跑通,再做CASPT2
PS : MOLCAS的输入比较麻烦,建议使用MOLPRO。

你好,谢谢你的回答!
我是研一的学生,这个输入文件是博三的师兄帮助写出来的,我只需要负责固定O-O、O-H键长改坐标就行,所以对于你说的并不特别理解,如果不耽误你的时间,麻烦你说的详细一点好吗?
感觉楼主是大牛,不知道楼主在使用MOLCAS方面有没有好的建议,师兄告诉我要看说明书,但我基础很薄弱,目前处于看不懂的状态,如果楼主对于MOLCAS的学习有好的建议,还请你多多赐教!
万分感谢!
3楼2014-04-05 20:15:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiangning198511

木虫 (正式写手)

引用回帖:
3楼: Originally posted by 小猪女魔头 at 2014-04-05 20:15:17
你好,谢谢你的回答!
我是研一的学生,这个输入文件是博三的师兄帮助写出来的,我只需要负责固定O-O、O-H键长改坐标就行,所以对于你说的并不特别理解,如果不耽误你的时间,麻烦你说的详细一点好吗?
感觉楼主 ...

如果是你师兄写的,最好的办法就是让他给你讲解一下,嘿嘿
PS:刚开始就接触这么高级的计算确实不容易,建议看看MOLCAS网站上辅导资料
4楼2014-04-06 10:00:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangshaohong

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by jiangning198511 at 2014-04-05 15:20:36
先将自带的例子测试一下,成功后说明你的程序没问题。
其次,将任务分解,先将CASSCF 跑通,再做CASPT2
PS : MOLCAS的输入比较麻烦,建议使用MOLPRO。

你好!我现在在用MOLCAS的CASSCF做结构优化,但是实验室里没有人用过这个软件,我想请问一下输入文件怎么写,谢谢楼主
5楼2017-01-04 10:48:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangshaohong

新虫 (初入文坛)

你好!我也是研一的,现在在用MOLCAS的CASSCF做结构优化,但是实验室里没有人用过这个软件,我想请问一下输入文件怎么写,谢谢楼主,万分感谢
6楼2017-01-04 10:50:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小猪女魔头 的主题更新
信息提示
请填处理意见