24小时热门版块排行榜    

查看: 3149  |  回复: 18
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

kent1022

金虫 (著名写手)

[求助] multi mrci计算问题,态平均?

计算一双原子分子体系的1Delta态或3Pi态,在C2v中分别是A1+A2和B1+B2.
若计算Delta态是否需要分别计算各个不同的对称性然后取平均?
先rhf,再multi,然后icMRCI,输入按以下
{rHF
occ,12,6,6,1      
wf,50, ,0
}      
{multi
occ,16,6,6,1
wf,50,1,0      !A1
lquant,2
wf,50,4,0      !A2  
lquant,2
}
{ci
occ,16,6,6,1
wf,50,1,0      !A1
wf,50,4,0      !A2    ! 1Delta   
}
计算分子核间距在r=1.5埃时没有问题,但是当r增加到1.55时,出现下列错误
SYYMETRY= 4  MS2= 0
? Error
? Open-shell singlet not implemented
? The problem occurs in scfpro

ERROR EXIT
CURRENT STACK:      MAIN

为什么计算r=1.5埃时没有出错?是输入哪里有问题?rhf中需要加入state,2么?

上述rhf,multi,mrci计算的是不是就是态平均的结果,是不是两个对称性分别为A1和A2的态各自计算后取其能量平均值?
请高手指点!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)

【答案】应助回帖

引用回帖:
7楼: Originally posted by kent1022 at 2013-07-13 22:25:43
那接下来的mrci计算用{ci;occ,16,6,6,2;wf,50, ,2;state,2}就可以了吧

还是也需要两个对称性分别计算
{ci;occ,16,6,6,2;wf,50,1,0}
{ci;occ,16,6,6,2;wf,50,4,0}

我不明白 wf,50,,2 (应该是wf,50,,0吧)中间为什么要省略...
你可以查一下manual,如果省略的话是不是默认为1...

这样的话,wf,50, ,2;state,2就相当于在算 wf,50,1,2 的ground state 和1st exited state而已...

另外,occ,16,6,6,1如果是文献上设置的,请搞清楚为啥这样设置...而不是看了文献和我的解释就乱选一个...
9楼2013-07-14 13:50:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答

fatpig8832

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
kent1022: 金币+30 2013-07-13 22:30:34
gmy1990: 金币+3 2013-07-14 07:21:27
楼主是单单要算1Delta态吧?

有两个问题...
首先,为什么multi选这个space?
16,6,6,1比起rhf的12,6,6,1来,只是在A1上增加了4个轨道,而A2的轨道只是1...那你算symmetry 4 时不出错才怪...1.5埃估计是碰巧算过去的...
建议:把每个轨道具体表示的是什么轨道搞清楚,在从中挑出需要的轨道组成你multi用的space。例如,如果你的计算需要用到3d轨道,那么必须把A1上的3d0,3d2+,B1上的3d1+,B2上的3d1-,以及A2上的3d2-都包含在你的active space才行,这样你的A2至少也要是2.

其次,你把两个态一起写,这就是在算态平均...得到的两个能量是已经态平均(默认50:50,比例可以用参数调)过的了,不需要你手动再平均一次。
rhf不能算态平均。

需要注意的是:active space必须要选好,要不然本来简并的两个态有可能算到不同能量;本来不简并的则没什么可能算到简并。
5楼2013-07-12 16:42:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kent1022

金虫 (著名写手)

引用回帖:
5楼: Originally posted by fatpig8832 at 2013-07-12 16:42:37
楼主是单单要算1Delta态吧?

有两个问题...
首先,为什么multi选这个space?
16,6,6,1比起rhf的12,6,6,1来,只是在A1上增加了4个轨道,而A2的轨道只是1...那你算symmetry 4 时不出错才怪...1.5埃估计是碰巧算过 ...

rhf和multi还都是按照一篇文章里的设置的,看来有些发表了的东西也不可靠!本人目前还在摸索阶段,所以先拿本人文章的例子练习。
我是要算1Delta态,如果把A2加到2,输入改为
{multi
occ,16,6,6,2
wf,50,1,0      !A1
lquant,2
wf,50,4,0      !A2  
lquant,2
}
这样计算的A1和A2的态平均结果是不是就是我要计算的1Delta态?
这个结果与我分别计算
{multi;occ,16,6,6,2;wf,50,1,0;lquant,2}
{multi;occ,16,6,6,2;wf,50,4,0;lquant,2}
得到的单点能能量的平均值应该是一样的吧?这就是所谓的态平均吧
6楼2013-07-13 22:19:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kent1022

金虫 (著名写手)

那接下来的mrci计算用{ci;occ,16,6,6,2;wf,50, ,2;state,2}就可以了吧

还是也需要两个对称性分别计算
{ci;occ,16,6,6,2;wf,50,1,0}
{ci;occ,16,6,6,2;wf,50,4,0}
7楼2013-07-13 22:25:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +5 @taotao 2026-03-19 5/250 2026-03-24 21:07 by greychen00
[考研] 一志愿北化315 求调剂 +3 akrrain 2026-03-24 3/150 2026-03-24 19:35 by 了了了了。。
[考研] 070300化学求调剂 +9 苑豆豆 2026-03-20 9/450 2026-03-24 17:15 by licg0208
[考研] 材料学硕,求调剂 6+3 糖葫芦888ll 2026-03-22 7/350 2026-03-24 17:11 by hello七七
[考研] 080500求调剂 +3 zzzzfan 2026-03-24 3/150 2026-03-24 16:38 by barlinike
[考研] 081700 调剂 267分 +9 迷人的哈哈 2026-03-23 9/450 2026-03-24 11:58 by 544594351
[考研] 【双一流院校新能源、环境材料,材料加工与模拟招收大量调剂】 +4 Higraduate 2026-03-22 7/350 2026-03-24 11:23 by 种大树
[考研] 284求调剂 +10 Zhao anqi 2026-03-22 10/500 2026-03-24 00:08 by Equinoxhua
[考研] 接收2026硕士调剂(学硕+专硕) +4 allen-yin 2026-03-23 6/300 2026-03-23 15:04 by 汪!?!
[考研] 306求调剂 +9 chuanzhu川烛 2026-03-18 9/450 2026-03-23 13:17 by luoyongfeng
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 280分求调剂 一志愿085802 +4 PUMPT 2026-03-22 7/350 2026-03-22 22:13 by 星空星月
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +11 枫桥ZL 2026-03-18 13/650 2026-03-22 20:26 by edmund7
[考研] 289求调剂 +7 怀瑾握瑜l 2026-03-20 7/350 2026-03-22 15:57 by ColorlessPI
[考研] 318求调剂 +4 plum李子 2026-03-21 7/350 2026-03-22 14:17 by ColorlessPI
[考研] 求助 +5 梦里的无言 2026-03-21 6/300 2026-03-21 17:51 by 学员8dgXkO
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 266求调剂 +3 哇呼哼呼哼 2026-03-20 3/150 2026-03-21 16:46 by barlinike
[考研] 296求调剂 +6 www_q 2026-03-18 10/500 2026-03-20 23:56 by JourneyLucky
[考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
信息提示
请填处理意见