24小时热门版块排行榜    

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

gmy1990

荣誉版主 (著名写手)

优秀版主优秀版主

[求助] H2键解离和Ar二聚体两个原子间的弱作用计算

采用Gaussian 03
H2的键解离过程应该用非限制性方法计算,但是为什么得到的和限制性方法的结果是一样的呢?MP2计算出来的曲线怎么又会如此奇葩?
对于Ar dimer的计算的话,势能曲线应该是在4左右的位置有一个最低点(Figure 3),而不是现在的曲线形状(Figure 2)。
@sobereva
@zhou2009
@yjcmwgk
@beefly
H2键解离和Ar二聚体两个原子间的弱作用计算
Figure1.jpg


H2键解离和Ar二聚体两个原子间的弱作用计算-1
Figure2.jpg


H2键解离和Ar二聚体两个原子间的弱作用计算-2
Figure3.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
gmy1990: 金币+5, ★★★★★最佳答案, 非常感谢! 2013-09-11 21:38:25
楼上的输入应该没问题的;你再仔细查看下是否少抄了一行。不行的话你还可以再试试guess=mix。

在解离部分你的rmp2和ump2实际上都是基于rhf参考态,参考态本身有定性问题,大概微扰理论发散了吧。

如果你算Ar2我觉得限制和非限制就是应该没有差别的。你用hf和mp2的话本身就不能准确描述色散相互作用,即使找到极小也没有什么参考价值。尽量要用ccsd或者casscf吧
4楼2013-09-11 21:22:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
gmy1990: 金币+5, ★★★★★最佳答案 2013-09-12 15:06:20
gmy1990: 金币+5, ★★★★★最佳答案, 非常感谢! 2013-09-12 15:06:51
gmy1990: 金币+15, ★★★★★最佳答案, 非常感谢! 2013-09-25 19:05:19
引用回帖:
5楼: Originally posted by gmy1990 at 2013-09-11 08:38:08
对于H2的解离,是否应该要用非限制性方法?但是由于都是基于闭壳层参考态,所以,得到的结果就是限制和非限制得到的结果都一样?

Ar的二聚体问题已经解决了,由于作用力很小,在设置合适的步长扫描后就没问题了 ...

H2用各种方法计算出来的势能面你可以在szabo&ostlund的书上看到,uhf和rhf是在227页,ump2和ump3是在375页。

H2在解离时相关能应该是0,uhf,ump2,casscf和ccsd都应该给出精确解-1a.u.;你可以以此来查看你uhf的解是否正确,如果离-1很远说明可能有问题。

对于Ar2,即使找到一个极小点,这个极小点并没有意义;不论mp2和hf都不描述色散相互作用,只能描述排斥力。你得到的极小很可能只是基组重叠误差造成的,不反映实际物理相互作用。要得到有物理意义的势能面还是需要换方法

» 本帖已获得的红花(最新10朵)

6楼2013-09-12 03:28:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)

引用回帖:
7楼: Originally posted by gmy1990 at 2013-09-11 15:23:37
第一个图就是做的H2的解离,casscf和ccsd都给出了很好的结果,但是uhf和ump2得出来的结果是不对的,而且限制性和非限制性得出来的都一个样。...

我把这两幅图扫下来了,只有纸板书,手机照相不太清楚
uhf,rhf和精确解
http://muchong.com/bbs/wpic.php?id=12650
uhr,ump2,ump3和精确解
http://muchong.com/bbs/wpic.php?id=12651

UHF大概从bohr2.3左右开始就应该开始低于RHF,然后逐渐收敛到和ccsd一样的能量。
8楼2013-09-12 10:03:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gmy1990 的主题更新
信息提示
请填处理意见