24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1311  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 286求调剂 +10 PolarBear11 2026-03-26 10/500 2026-03-28 07:52 by baoball
[考研] 压国家一区线,求导师收留,有恩必谢! +3 迷人的哈哈 2026-03-28 3/150 2026-03-28 07:34 by ACS Nano——
[考研] 材料与化工085600,总分304,本科有两篇sci参与,求调剂 +10 幸运的酱酱 2026-03-22 12/600 2026-03-27 16:08 by muchong357
[考研] 317求调剂 +5 十闲wx 2026-03-24 5/250 2026-03-27 13:48 by 杨杨杨紫
[考研] 一志愿北京化工大学材料与化工 264分各科过A区国家线 +10 哈哈157349 2026-03-21 10/500 2026-03-27 13:06 by zzll406
[考研] 298调剂 +3 jiyingjie123 2026-03-27 3/150 2026-03-27 11:57 by wxiongid
[考研] 求调剂 +3 刘柯@ 2026-03-24 4/200 2026-03-27 11:28 by shangxh
[考研] 一志愿吉大071010,316分求调剂 +3 xgbiknn 2026-03-27 3/150 2026-03-27 10:36 by guoweigw
[考研] 求调剂 +6 林之夕 2026-03-24 6/300 2026-03-27 08:38 by hypershenger
[考研] 289求调剂 +17 硕星赴 2026-03-23 17/850 2026-03-26 16:18 by 不吃魚的貓
[考研] 263求调剂 +6 yqdszhdap- 2026-03-22 10/500 2026-03-26 13:11 by 公瑾逍遥
[考研] 309求调剂 +4 gajsj 2026-03-25 5/250 2026-03-26 00:27 by Dyhoer
[考研] 302求调剂 +4 锦衣卫藤椒 2026-03-25 4/200 2026-03-25 16:29 by 功夫疯狂
[考研] 材料考研调剂生 +3 黄粱一梦千年 2026-03-24 3/150 2026-03-24 17:00 by barlinike
[考研] 307求调剂 +3 余意卿 2026-03-21 6/300 2026-03-24 15:03 by 余意卿
[考研] 一志愿国科过程所081700,274求调剂 +3 三水研0水立方 2026-03-23 3/150 2026-03-23 23:11 by MajorWen
[考研] 接收2026硕士调剂(学硕+专硕) +4 allen-yin 2026-03-23 6/300 2026-03-23 15:04 by 汪!?!
[考研] 一志愿东华大学化学070300,求调剂 +7 2117205181 2026-03-21 8/400 2026-03-22 22:55 by chixmc
[考研] 306求调剂 +5 来好运来来来 2026-03-22 5/250 2026-03-22 16:17 by BruceLiu320
[考研] 一志愿深大,0703化学,总分302,求调剂 +4 七月-七七 2026-03-21 4/200 2026-03-21 18:20 by 学员8dgXkO
信息提示
请填处理意见