24小时热门版块排行榜    

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

luqing6879

木虫 (著名写手)

[求助] 请问mrci 和mrci+q的区别

我用molpro计算势能面,当用mrci的能量值进行拟合时,得到的误差符合所需精度,但是用mrciq进行拟合时,拟合误差就很大。

请问mrci 和mrci+q 的区别是什么?mrci+q“错误”地考虑了什么因素造成了这么大的误差?

谢谢。^_^
回复此楼
cishan.goodweb.cn/慈善点击,点滴帮助。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luqing6879

木虫 (著名写手)

引用回帖:
2楼: Originally posted by beefly at 2012-09-30 01:20:18
+Q部分大致对应于CISD与CCSD两种方法之间的差别。+Q有很多计算公式,有好有坏,但是一般差不太多。Molpro在2010版本以前,多根MRCI+Q计算存在BUG,可能导致错误的+Q能量校正。
http://www.molpro.net/info/2012.1/ ...

谢谢回复。
1
我看到过评论说“+Q部分大致对应于CISD与CCSD两种方法之间的差别”。 但是不是很理解为什么要进行这种校正。一个是单参考方法,一个是多参考方法。如果体系有multi-reference性质,为什么要关心ccsd呢?一个令我满意的解释是来自维基百科的:
It allows one to estimate the value of the full configuration interaction energy from a limited configuration interaction expansion result, although more precisely it estimates the energy of configuration interaction up to quadruple excitations (CISDTQ) from the energy of configuration interaction up to double excitations (CISD).
请指教。

另外搭车问一下,molpro是不是不支持逻辑运算符啊?例如:
if (k.eq10  .or.  k.eq.20) then
...
endif
总是报错。但是把上面的if分开来写成两个if就没事。
cishan.goodweb.cn/慈善点击,点滴帮助。
3楼2012-09-30 02:27:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

fmklz

木虫 (正式写手)

搞电器工业化的化学打工仔

【答案】应助回帖

if (k.eq10  .or.  k.eq.20) then
...
endif
支持逻辑运算符,k.eq10不知你的点是不是忘了k.eq.10
宁在一思进,莫在一思停
4楼2012-10-10 22:07:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luqing6879

木虫 (著名写手)

引用回帖:
4楼: Originally posted by fmklz at 2012-10-10 22:07:45
if (k.eq10  .or.  k.eq.20) then
...
endif
支持逻辑运算符,k.eq10不知你的点是不是忘了k.eq.10

Thanks for the reply.

Sorry for the typo. I have the dot in my original input, but it doesn't work.

Does the logical expression work for you?
cishan.goodweb.cn/慈善点击,点滴帮助。
5楼2012-10-11 07:40:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fmklz

木虫 (正式写手)

搞电器工业化的化学打工仔

【答案】应助回帖

★ ★ ★ ★ ★
luqing6879: 金币+5, ★★★很有帮助, thank you. 2012-10-17 00:50:35
I have used the command:
if (i.GE.5.AND.i.LT.9) THEN
...
endif

And this expresssion works for me.
Good luck
宁在一思进,莫在一思停
6楼2012-10-15 17:17:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:43 by OEbVnUOu01ol
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 03:31 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 8/400 2026-08-22 18:29 by 淀粉搬运工
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见