24小时热门版块排行榜    

查看: 1417  |  回复: 4

sevenci5460

金虫 (小有名气)

[求助] MRCI 如何指定参考态? 已有1人参与

计算势能面时候,输入如下:
......
{multi;occ,9,2;closed,2;pspace,2d0;config,csf;
wf,15,1,1;state,6;
maxiter,35,600;
natorb,2140.2;}
ecas1=energy(1)
ecas2=energy(2)
ecas3=energy(3)
ecas4=energy(4)
ecas5=energy(5)
ecas6=energy(6)

{mrci;occ,9,2;closed,2;wf,15,1,1;state,4;maxiter,400,1200;}
......
——————————————————————————————————————————————————

部分结果文件如下:

********************************************************************************************************************
SETTING ECAS1          =      -109.62135417  HARTREE
SETTING ECAS2          =      -109.37561668  HARTREE
SETTING ECAS3          =      -109.33520972  HARTREE
SETTING ECAS4          =      -109.31983119  HARTREE
SETTING ECAS5          =      -109.28870659  HARTREE
SETTING ECAS6          =      -109.25822546  HARTREE

1PROGRAM * CI (Multireference internally contracted CI)     Authors: H.-J. Werner, P.J. Knowles, 1987
......
Reference wavefunction optimized for reference space (refopt=1)

State     Reference Energy
   1        -109.62135417
   2        -109.37561668
   3        -109.31983119
————————————————————————————————
可以看出来,计算MRCI时选取的reference energy 并非CASSCF算出来的前三个态的能量,而是第1,2,4个,导致MRCI结果有误,
请问 如何指定相应的参考态能量呢?
用refstate?,不知具体怎么写,求教!
回复此楼
做一个和自己赛跑的人
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★
gkf高: 金币+2 2014-02-27 21:43:33
你不是写了state,4;吗?怎么会只有3个态?
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
2楼2014-02-27 20:59:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sevenci5460

金虫 (小有名气)

引用回帖:
2楼: Originally posted by beefly at 2014-02-27 20:59:20
你不是写了state,4;吗?怎么会只有3个态?

不好意思,应该是 state,3,那是我后来改成state,4的时候,就能算对了,可是这样计算时间会长一点
做一个和自己赛跑的人
3楼2014-02-28 14:24:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
sevenci5460: 金币+10, ★★★★★最佳答案 2014-02-28 20:46:17
gmy1990: 金币+2 2014-03-04 05:55:47
{ci;....
option,nstati=6;
wf,....;state,3;}
这样就不会出现丢失某个态的问题了

如果不关心跃迁偶极矩,第二行可以改成option,ITRDM=-1,nstati=6;,计算还能再快一点
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
4楼2014-02-28 15:43:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sevenci5460

金虫 (小有名气)

引用回帖:
4楼: Originally posted by beefly at 2014-02-28 15:43:28
{ci;....
option,nstati=6;
wf,....;state,3;}
这样就不会出现丢失某个态的问题了

如果不关心跃迁偶极矩,第二行可以改成option,ITRDM=-1,nstati=6;,计算还能再快一点

谢谢! 恩 这样就对了!
做一个和自己赛跑的人
5楼2014-02-28 20:46:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sevenci5460 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 u1xkze8ybe 2026-07-30 7/350 2026-07-31 12:10 by syhhffkrym
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 u1xkze8ybe 2026-07-30 3/150 2026-07-31 12:08 by syhhffkrym
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 u1xkze8ybe 2026-07-30 6/300 2026-07-31 12:05 by syhhffkrym
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +6 锦衣卫寒战 2026-07-28 6/300 2026-07-31 11:27 by syl200707
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 小木虫看见有人已经知道结果了 +15 1234567wang 2026-07-25 19/950 2026-07-31 08:30 by ZJTJZ
[基金申请] 2026年国自然面上资助率 +8 布布和一二 2026-07-30 8/400 2026-07-31 01:51 by alongwaytogo
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 微信指数没变化,科研之友没阅读 +14 wangze12014 2026-07-28 17/850 2026-07-30 20:10 by kissu88
[考博] 申请2027年材料/化学类博士 10+3 考研老狗? 2026-07-25 5/250 2026-07-30 19:33 by 白色木鱼
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[基金申请] 今年的WR进展到哪一步了? +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 时间戳他又来了 +14 晓晓爱翠翠 2026-07-26 17/850 2026-07-29 16:30 by 呵呵二哥
[有机交流] 这个自发加氧反应的机理是什么? 50+4 YUAN2273 2026-07-25 7/350 2026-07-29 12:27 by xiaobenpu
[基金申请] 没消息就是被刷了呗 +13 lambert2014 2026-07-24 14/700 2026-07-29 10:38 by msh145800
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
[分析] 水蒸气蒸馏法乳化现象 +3 月光娇娇 2026-07-26 3/150 2026-07-26 23:25 by caiyun
信息提示
请填处理意见