| 查看: 2548 | 回复: 8 | |||
luqing6879木虫 (著名写手)
|
[求助]
请问mrci 和mrci+q的区别
|
|
我用molpro计算势能面，当用mrci的能量值进行拟合时，得到的误差符合所需精度，但是用mrciq进行拟合时，拟合误差就很大。 请问mrci 和mrci+q 的区别是什么？mrci+q“错误”地考虑了什么因素造成了这么大的误差？ 谢谢。^_^ |
回复此楼» 收录本帖的淘帖专辑推荐
 我的收藏 |
» 猜你喜欢
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有11人回复
-
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有10人回复
-
物理化学论文润色/翻译怎么收费?
已经有267人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
已知面上挂了，小学校实验方向无平台，感觉太难了
已经有26人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
求助各位大佬，球球了，在这科萨头
已经有0人回复
-
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340
已经有2人回复
-
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
-
分子动力学模拟合作需求
已经有0人回复
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 12859.4
- 散金: 27307
- 红花: 300
- 帖子: 4667
- 在线: 925.7小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
2楼2012-09-30 01:20:18
luqing6879
木虫 (著名写手)
- 应助: 8 (幼儿园)
- 金币: 3565.5
- 散金: 249
- 红花: 3
- 帖子: 1907
- 在线: 807.8小时
- 虫号: 614833
- 注册: 2008-09-28
- 专业: 理论和计算化学
|
谢谢回复。 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就没事。 |
3楼2012-09-30 02:27:17
fmklz
木虫 (正式写手)
搞电器工业化的化学打工仔
- 应助: 2 (幼儿园)
- 金币: 2874.3
- 散金: 50
- 红花: 4
- 帖子: 352
- 在线: 383.4小时
- 虫号: 770620
- 注册: 2009-05-14
- 性别: GG
- 专业: 能源化工
4楼2012-10-10 22:07:45
luqing6879
木虫 (著名写手)
- 应助: 8 (幼儿园)
- 金币: 3565.5
- 散金: 249
- 红花: 3
- 帖子: 1907
- 在线: 807.8小时
- 虫号: 614833
- 注册: 2008-09-28
- 专业: 理论和计算化学
5楼2012-10-11 07:40:42
fmklz
木虫 (正式写手)
搞电器工业化的化学打工仔
- 应助: 2 (幼儿园)
- 金币: 2874.3
- 散金: 50
- 红花: 4
- 帖子: 352
- 在线: 383.4小时
- 虫号: 770620
- 注册: 2009-05-14
- 性别: GG
- 专业: 能源化工
6楼2012-10-15 17:17:38
virtualzx
木虫 (著名写手)
- QC强帖: 2
- 应助: 263 (大学生)
- 金币: 7161.3
- 红花: 54
- 帖子: 1605
- 在线: 317.6小时
- 虫号: 2069080
- 注册: 2012-10-18
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
luqing6879: 金币+10, ★★★很有帮助, 谢谢。^_^ 2012-10-18 07:01:47
luqing6879: 金币+10, ★★★很有帮助, 谢谢。^_^ 2012-10-18 07:01:47
|
+Q就是davidson correction，是一个非常简单快捷的微扰理论估计四阶激发对能量贡献的方法。不需要额外的计算，只用参考态(scf或mcscf)能量，ci能量，和ci波函数与参考态波函数的交叠三个数字，就可以直接外推出来。 对mrci，一般有三种算法 (1-c0^2)*(E[ci]-E[ref]) (1-c0^2)*(E[ci]-E[ref])/(c0^2) (1-c0^2)*(E[ci]-E[ref])/(2*c0^2-1) 这不是一个很干净的矫正，但是经常有效，而且又不花时间，所以比较常用。当然，因为不干净，不准的情况也有发生。 如果矫正结果不连续，有几个可能性。首先，这个矫正是基于微扰假定的，所以如果你的mcscf和mrci的解已经出现定性差别，这个矫正就完全不成立了。也就是如果c0比较小，结果就不可信。一般c0应该在0.8以上 再有可能是c0不连续造成的。有可能是你的mcscf不连续，出现了某个点轨道忽然翻转的情况；同时因为要把mcscf的波函数和mrci的波函数一一对应，如果对应找错了就很容易出现问题。你可以在输出文件中找到每个mrci解和所有参考态的交叠，看一下是不是选错了。 |
7楼2012-10-18 03:39:16
8楼2012-10-19 09:22:36
luqing6879
木虫 (著名写手)
- 应助: 8 (幼儿园)
- 金币: 3565.5
- 散金: 249
- 红花: 3
- 帖子: 1907
- 在线: 807.8小时
- 虫号: 614833
- 注册: 2008-09-28
- 专业: 理论和计算化学
9楼2012-10-19 23:39:15