24小时热门版块排行榜    

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

luqing6879

木虫 (著名写手)

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

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

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

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

luqing6879

木虫 (著名写手)

引用回帖:
8楼: Originally posted by lyq0120 at 2012-10-19 09:22:36
楼主好。我是刚刚接触势能面,现在在算一些数据,主要是算分子键长与能量的关系,也就是单点能。用origin画出了单点能的简单图像。但是怎么分析,以后再怎么做就不清楚了。比如说拟合,怎么拟合啊?请教楼主。。。。

我可能理解你的意思有误,不过只算一个点是不会得到势能面的。另外你是算近平衡点得势能面还是全局势能面呢?

你需要找到与粒子对称性合适的一组格点,然后算每个格点的单点能,最后用最小二乘或其他什么方法拟合。
cishan.goodweb.cn/慈善点击,点滴帮助。
9楼2012-10-19 23:39:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

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

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的回帖
信息提示
请填处理意见