| 查看: 1644 | 回复: 9 | |||
[求助]
用MOLPRO中的基组外推发现结果有跃变？ 已有2人参与
|
|
计算了一体系的势能曲线，根据avqz，5z，6z的结果用molpro内置的外推方程外推，参考能量用E(n)=ECBS+A*exp(-C*n)，相关能用(n+ 0.5)^{-X} with X= 4.0。 发现qz，5z，6z各自的曲线都是光滑的，cbs结果除了一个小区域内的外推结果相对整个曲线的变化趋势有明显偏离外，其它地方均光滑而连续。如图，  Graph1.jpg  Graph2.jpg 有没有人也遇到过这种问题，应该怎样解决？ |
回复此楼» 猜你喜欢
淮北师范大学 化学与化工科学学院 电催化课题组还有1名硕士研究生招生名额
已经有0人回复
-
化学相关硕士研究生招收调剂-湖北师范大学 （快来捡漏，最新调剂信息，过国家线即可）
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有55人回复
-
河南理工大学化学化工学院招收学硕调剂生
已经有0人回复
-
黑龙江科技大学化工 环境工程调剂
已经有22人回复
-
莆田学院环境与生物工程学院资源与环境专硕（环境工程方向） 开放调剂招生
已经有13人回复
-
浙江海洋大学石油化工与环境学院招收工科调剂生（学硕、专硕均可，线上面试）
已经有0人回复
-
浙江海洋大学石油化工学院接收工科调剂（学硕、专硕均可，近期线上面试）
已经有8人回复
-
浙江海洋大学石油与天然气工程专业接收工科调剂（13日早上截止，线上面试）
已经有0人回复
-
招材料相关专业调剂生
已经有10人回复
-
考研调剂～湖南工程学院材料与化工学院
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用MOLPRO构造原子自然轨道（ANO）基函数
已经有3人回复
- 高精度量子化学计算软件Molpro 入门及理论精讲课程
已经有3人回复
- Molpro 讲习班 2012年第一期
已经有3人回复
luqing6879
木虫 (著名写手)
- 应助: 8 (幼儿园)
- 金币: 3566
- 散金: 249
- 红花: 3
- 帖子: 1908
- 在线: 808.5小时
- 虫号: 614833
- 注册: 2008-09-28
- 专业: 理论和计算化学
2楼2014-05-10 00:51:43
3楼2014-05-10 22:06:55
孙1986
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1046.9
- 散金: 1087
- 红花: 19
- 帖子: 455
- 在线: 243.1小时
- 虫号: 1405621
- 注册: 2011-09-17
- 专业: 原子和分子物理
4楼2014-05-11 12:48:32
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
kent1022: 金币+20 2014-05-12 16:14:37
感谢参与，应助指数 +1
kent1022: 金币+20 2014-05-12 16:14:37
|
correlation不光滑，本质上还是因为reference不光滑。 不过，通常即使reference不光滑，加上correlation之后整体也是光滑的。所以，参考能和相关能分别做外推，可能不是一个good idea. 比如使用RHF或者UHF计算闭壳层体系的UCCSD能量，键稍微长一点，RHF和UHF是不一样的，但最后总的UCCSD能量都是一样的。 一般的，reference到5Z也就收敛了，QZ也基本收敛，所以可以对总的能量，用外推correlation能量的方法外推两个最大基组的能量。 Science 342:1499，他们新的FH2势能面，就是外推AV5Z和AV6Z的能量，因为HF能量在AV5Z已经收敛，所以可以跟相关能放在一起做外推。 [ 发自手机版 http://muchong.com/3g ] |
5楼2014-05-11 13:55:08
luqing6879
木虫 (著名写手)
- 应助: 8 (幼儿园)
- 金币: 3566
- 散金: 249
- 红花: 3
- 帖子: 1908
- 在线: 808.5小时
- 虫号: 614833
- 注册: 2008-09-28
- 专业: 理论和计算化学
|
加入相关能后是可以使能量曲线看上去更平滑。但是如virtuallzx之前在楼主mrci里的那个帖子所说，看上去光滑的片段在放大后可能还是有间断。所以如4 5楼所判断的，问题还是出现在hf。 首先看看那些点的t1是多少吧。如果大的话可能说明在该点的附近有临近的激发态，做个2state casscf看看能量差多少。 至于基组外推，分开来推和合起来推我都试过，对于小分子来说，从绝对外推值来看，二者相差不大。 印象中分开推的人多一点。就我个人看，从本质上说，外推公式都是经验公式，公式的形式必须要拟合能量曲线的趋势，所以只要拟合的好，怎么推都可以。但问题是hf能和相关能的趋势往往是不尽相同的。所以用不同的公式分别拟合二者，才更符合科学精神。合起来推的结果，由于相关比hf能小几个数量级，所以相关能的信息完全被hf淹没了。这也是为什么用两种方法推差别不大的原因。 |
6楼2014-05-12 00:33:27
7楼2014-05-12 14:20:41
luqing6879
木虫 (著名写手)
- 应助: 8 (幼儿园)
- 金币: 3566
- 散金: 249
- 红花: 3
- 帖子: 1908
- 在线: 808.5小时
- 虫号: 614833
- 注册: 2008-09-28
- 专业: 理论和计算化学
8楼2014-05-12 14:59:39
9楼2014-05-12 20:24:12
10楼2014-05-12 22:52:23
20