24小时热门版块排行榜    

查看: 1933  |  回复: 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的回帖

小猪女魔头

铜虫 (初入文坛)

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

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 可靠消息,8月21号(明天)放榜 +11 医学老男孩 2026-08-20 13/650 2026-08-20 10:50 by anjeeshine
[基金申请] 时间戳变了,能看出什么问题? +13 基诺咪客 2026-08-17 17/850 2026-08-20 10:43 by wangyue2015
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 21/1050 2026-08-20 10:37 by Vivilian
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +15 Ldrop2023 2026-08-13 19/950 2026-08-20 10:30 by 陈晨晨陈啊
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +55 医学老男孩 2026-08-13 125/6250 2026-08-20 09:02 by 医学老男孩
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见