24小时热门版块排行榜    

查看: 1989  |  回复: 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

木虫 (正式写手)

引用回帖:
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的回帖
查看全部 6 个回答

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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[职场人生] 40+中年人失业后,还能找到好工作吗? +11 xqy2009 2026-10-13 13/650 2026-10-15 21:32 by bear2007
[有机交流] 用NBS溴化时溶剂除了CCl4外还有什么 +3 18834002860 2026-10-12 6/300 2026-10-15 14:30 by 长沙摩尔信息科�
[基金申请] 申请科研项目纯靠实力能成功吗? 200+7 David1982 2026-10-13 11/550 2026-10-15 10:30 by zyqchem
[硕博家园] 极限拉扯,如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
信息提示
请填处理意见