| 查看: 3147 | 回复: 18 | ||
[求助]
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» 猜你喜欢
国内青年破解费曼难题，从分子层面阐释液体蒸发（气液相变）微观过程
已经有5人回复
-
研究发现一个新的临界系数，独立于临界压缩因子
已经有7人回复
-
物理化学论文润色/翻译怎么收费?
已经有187人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
已知面上挂了，小学校实验方向无平台，感觉太难了
已经有26人回复
-
求助各位大佬，球球了，在这科萨头
已经有0人回复
-
七嗪类物质合成求助
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- molpro说明书的例子运行不了？？
已经有4人回复
- molpro试做旋轨耦合,怎么在输出中找omega及其对应的能量
已经有5人回复
- 【求助】Molpro如何执行旋轨耦合计算 如何执行态平均计算?
已经有10人回复
14楼2017-12-11 21:59:53
 |
2楼2013-07-05 21:42:06
|
|
3楼2013-07-11 22:33:37
|
|
4楼2013-07-12 11:08:34
fatpig8832
铁杆木虫 (著名写手)
- 应助: 16 (小学生)
- 金币: 7420.2
- 散金: 31
- 红花: 28
- 帖子: 1438
- 在线: 610.9小时
- 虫号: 770403
- 注册: 2009-05-14
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
kent1022: 金币+30 2013-07-13 22:30:34
gmy1990: 金币+3 2013-07-14 07:21:27
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
|
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
7楼2013-07-13 22:25:43
fatpig8832
铁杆木虫 (著名写手)
- 应助: 16 (小学生)
- 金币: 7420.2
- 散金: 31
- 红花: 28
- 帖子: 1438
- 在线: 610.9小时
- 虫号: 770403
- 注册: 2009-05-14
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
kent1022: 金币+20 2013-07-16 10:36:13
kent1022: 金币+20 2013-07-16 10:36:13
|
首先请把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
fatpig8832
铁杆木虫 (著名写手)
- 应助: 16 (小学生)
- 金币: 7420.2
- 散金: 31
- 红花: 28
- 帖子: 1438
- 在线: 610.9小时
- 虫号: 770403
- 注册: 2009-05-14
9楼2013-07-14 13:50:27
18766315007
新虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 6318.5
- 散金: 100
- 红花: 6
- 帖子: 1261
- 在线: 263.5小时
- 虫号: 2979660
- 注册: 2014-02-20
- 专业: 原子和分子物理
10楼2017-12-11 16:45:20
