北京石油化工学院2026年研究生招生接收调剂公告
查看: 2272  |  回复: 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

木虫 (著名写手)

引用回帖:
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的回帖
查看全部 11 个回答

beefly

专家顾问 (职业作家)

地沟油冶炼专家

1,试试这个输入:
#P uhf/3-21g guess=alter iop(5/15=2)

test

0 1
h
h 1 0.7

1 2

2,需要用大基组,做基组外推
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
2楼2013-09-11 09:26:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gmy1990

荣誉版主 (著名写手)

优秀版主优秀版主

引用回帖:
2楼: Originally posted by beefly at 2013-09-11 09:26:22
1,试试这个输入:
#P uhf/3-21g guess=alter iop(5/15=2)

test

0 1
h
h 1 0.7

1 2

2,需要用大基组,做基组外推

也不行
按道理是hf的结果就是随着两原子距离增加而能量一直递减。还是无法得到限制性和非限制性的差别。
关于Ar 二聚体的情况是这样的,我换了He二聚体可以得到满意的结果,就是换成Ar原子后就出现这样的问题,基组已经用到aug-cc-pvtz.
这样一来,所以,我就减小了步长到0.005,方才可以解决了这个问题。是一个极小的energy minimum point

[ Last edited by gmy1990 on 2013-9-11 at 19:45 ]
3楼2013-09-11 19:32:17
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +6 白QF 2026-03-21 7/350 2026-03-28 00:42 by 白QF
[考研] 285求调剂 +4 AZMK 2026-03-27 7/350 2026-03-27 20:59 by AZMK
[考研] 学硕274求调剂 +5 Li李鱼 2026-03-26 5/250 2026-03-27 20:51 by 热情沙漠
[考研] 考研化学308分求调剂 +10 你好明天你好 2026-03-23 12/600 2026-03-27 14:43 by shangxh
[考研] 22408 359分调剂 +3 Qshers 2026-03-27 3/150 2026-03-27 12:22 by wxiongid
[考研] 304求调剂 +3 曼殊2266 2026-03-27 3/150 2026-03-27 10:17 by guoweigw
[考研] 材料学硕,求调剂 6+5 糖葫芦888ll 2026-03-22 10/500 2026-03-27 08:18 by hypershenger
[考研] 321求调剂 +5 材料cailiao 2026-03-21 5/250 2026-03-26 20:41 by fmesaito
[考研] 材料调剂 5+4 想要一壶桃花水 2026-03-25 10/500 2026-03-26 19:56 by 不吃魚的貓
[考研] 一志愿厦门大学化学学硕307求调剂 +8 y7czhao 2026-03-26 8/400 2026-03-26 19:51 by 不吃魚的貓
[考研] 材料科学与工程 317求调剂 +4 JKSOIID 2026-03-26 4/200 2026-03-26 15:58 by 不吃魚的貓
[考研] 263求调剂 +6 yqdszhdap- 2026-03-22 10/500 2026-03-26 13:11 by 公瑾逍遥
[考研] 一志愿天津大学339材料与化工求调剂 +3 江往卖鱼 2026-03-26 3/150 2026-03-26 09:42 by 王小欠i
[考研] 080500求调剂 +3 zzzzfan 2026-03-24 3/150 2026-03-24 16:38 by barlinike
[考研] 284求调剂 +3 yanzhixue111 2026-03-23 6/300 2026-03-23 22:58 by pswait
[考研] 361求调剂 +3 Glack 2026-03-22 3/150 2026-03-23 22:03 by fuyu_
[考研] 333求调剂 +3 ALULU4408 2026-03-23 3/150 2026-03-23 19:04 by macy2011
[考研] 工科0856求调剂 +5 沐析汀汀 2026-03-21 5/250 2026-03-23 17:56 by 海瑟薇-
[考研] 275求调剂 +6 shansx 2026-03-22 8/400 2026-03-22 15:27 by barlinike
[考研] 求调剂 +3 .m.. 2026-03-21 4/200 2026-03-21 16:25 by barlinike
信息提示
请填处理意见