应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Molpro/Molcas » multi mrci计算问题，态平均？
191/2返回列表12下一页
查看: 3533  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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的态各自计算后取其能量平均值？
请高手指点！
回复此楼

» 本帖已获得的红花（最新10朵）

18766315007 18766315007

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-07-04 23:20:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

kent1022

金虫 (著名写手)
引用回帖:
13楼: Originally posted by 18766315007 at 2017-12-11 18:43:23
好的，谢谢！还有一个问题。就是OCC卡后面的轨道是怎么确定的？像同核双原子分子D2h点群，ag,b3u,b2u,b1g,b1u,b2g,b3g,au对应分子的什么轨道。也就是说，给定一个双原子分子的电子态，怎么确定occ卡后面的数值？...

让程序计算hf，输出里面就有默认的轨道
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

18766315007
14楼2017-12-11 21:59:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

kent1022

金虫 (著名写手)
2楼2013-07-05 21:42:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kent1022

金虫 (著名写手)
3楼2013-07-11 22:33:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kent1022

金虫 (著名写手)
4楼2013-07-12 11:08:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

fatpig8832

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
kent1022: 金币+20 2013-07-16 10:36:13
引用回帖:
6楼: Originally posted by kent1022 at 2013-07-13 22:19:43
rhf和multi还都是按照一篇文章里的设置的，看来有些发表了的东西也不可靠！本人目前还在摸索阶段，所以先拿本人文章的例子练习。
我是要算1Delta态，如果把A2加到2，输入改为
{multi
occ,16,6,6,2
wf,50,1,0   ...

首先请把active space搞清楚，为什么你要选16,6,6,2而不是15,6,6,2或14,7,7,2或其他的呢？也许16,6,6,2是合理的，只要你搞清楚每个轨道所代表的是哪个具体轨道，能够说服自己就行...

然后，“态平均”和“分开算再平均”是两回事...举例说明：假设你用“态平均”算了两个态，能量分别是-100.100000和-100.200000，这两个能量就是每个态的“态平均”能量，不需要再加起来除以2...
赞一下 回复此楼
8楼2013-07-14 13:45:22
已阅   回复此楼   关注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的回帖

18766315007

新虫 (著名写手)
送红花一朵
大神，我初学，遇到了同样的问题，你弄明白了吗？能不能讲一下？
赞一下 回复此楼
10楼2017-12-11 16:45:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kent1022 的主题更新
191/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +55 医学老男孩 2026-08-13 123/6150 2026-08-19 23:44 by stan_et
[基金申请] 今天基金会出结果吗？20260819 +13 kkkl_v 2026-08-19 13/650 2026-08-19 22:37 by 宝贝虫子
[基金申请] 时间戳变了，能看出什么问题？ +12 基诺咪客 2026-08-17 16/800 2026-08-19 22:04 by kk 的小木虫
[基金申请] 今天放榜没戏了吧 +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
[基金申请] 什么时候开奖？ +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 重要消息，中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[基金申请] 明天放榜？ +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
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护，明天很可能放榜！ +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
信息提示
关闭
请填处理意见
关闭