24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1308  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 315分求调剂 +7 26考研上岸版26 2026-03-26 7/350 2026-03-28 04:05 by fmesaito
[考研] 求调剂 +6 白QF 2026-03-21 7/350 2026-03-28 00:42 by 白QF
[考研] 一志愿南师大0703化学 275求调剂 +4 Ripcord上岸 2026-03-27 4/200 2026-03-27 17:00 by zhyzzh
[考研] 274求调剂 +17 顾九笙要谦虚 2026-03-24 23/1150 2026-03-27 15:16 by caszguilin
[考研] 一志愿华东理工大学081700,初试分数271 +6 kotoko_ik 2026-03-23 7/350 2026-03-27 12:29 by 惠州彭于晏
[考研] 316求调剂 +5 Pigcasso 2026-03-24 5/250 2026-03-27 12:10 by zhshch
[考研] 342求调剂 +3 加油a李zs 2026-03-26 3/150 2026-03-27 00:29 by wxiongid
[考研] 一志愿郑州大学,080500学硕,总分317分求调剂 +4 举个栗子oi 2026-03-24 5/250 2026-03-26 23:15 by 不吃魚的貓
[考研] 321求调剂 +5 材料cailiao 2026-03-21 5/250 2026-03-26 20:41 by fmesaito
[考研] 309求调剂 +4 gajsj 2026-03-25 5/250 2026-03-26 00:27 by Dyhoer
[考研] 290分调剂求助 +3 吉祥止止陈 2026-03-25 3/150 2026-03-25 19:58 by barlinike
[考研] 26考研-291分-厦门大学(085601)-柔性电子学院材料工程专业求调剂 +3 min3 2026-03-24 4/200 2026-03-25 18:22 by xcjcqu
[考研] 考研一志愿苏州大学初始315(英一)求调剂 +3 sbdksD 2026-03-24 4/200 2026-03-25 18:16 by xcjcqu
[考研] 生物技术与工程 +3 1294608413 2026-03-25 4/200 2026-03-25 18:02 by 1294608413
[考研] 【2026考研调剂】制药工程 284分 求相关专业调剂名额 +4 袁奂奂 2026-03-25 8/400 2026-03-25 14:32 by lbsjt
[考研] 318求调剂 +5 plum李子 2026-03-21 8/400 2026-03-25 09:26 by aa331100
[考研] 340求调剂 +5 话梅糖111 2026-03-24 5/250 2026-03-25 06:53 by ilovexiaobin
[考研] 调剂 +4 13853210211 2026-03-24 4/200 2026-03-24 19:44 by ms629
[考研] 材料专硕找调剂 +5 哈哈哈吼吼吼哈 2026-03-23 5/250 2026-03-24 19:07 by 了了了了。。
[考研] 085404电子信息284分求调剂 +4 13659058978 2026-03-24 4/200 2026-03-24 12:15 by syl20081243
信息提示
请填处理意见