【答案】应助回帖

» 猜你喜欢

中国地质大学（武汉）—国家级青年人才杨明教授组-招收博士-新能源材料化学及催化材料已经有32人回复
求光催化产过氧化氢的PPT ！！！毕业答辩需要已经有1人回复
物理化学论文润色/翻译怎么收费? 已经有127人回复
[电子教材]天津大学物理化学（第七版）上、下册已经有1人回复
湖南大学招收锂离子电池、离子液体、高分子与凝胶材料、电催化方向博士生多名已经有5人回复
半透性容器进行高湿或者低湿研究已经有0人回复
关于五水硝酸铋还原得到单质铋的问题已经有1人回复
阴离子交换膜电解二氧化碳还原已经有0人回复
山东省青年基金什么时候拨款已经有30人回复

» 本主题相关价值贴推荐，对您同样有帮助:

二聚体的分离---分子量1000 已经有6人回复
血红蛋白四聚体二聚体高效液相和飞行时间质谱已经有2人回复
【sobereva个人文集】乱弹DFT-D（从sobereva的博客中转载过来的）已经有27人回复

1楼2013-09-11 06:54:26

beefly

专家顾问 (职业作家)

地沟油冶炼专家

专家经验: +457
QC强帖: 12
应助: 86 (初中生)
金币: 12754.9
散金: 27307
红花: 300
帖子: 4663
在线: 925.3小时
虫号: 408372
注册: 2007-06-21
性别: GG
专业: 理论和计算化学
管辖: 计算模拟

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

gmy1990

荣誉版主 (著名写手)

QC强帖: 1
应助: 163 (高中生)
贵宾: 0.782
金币: 13422.7
散金: 526
红花: 42
沙发: 4
帖子: 1627
在线: 1019.7小时
虫号: 745920
注册: 2009-04-11
专业: 理论和计算化学
管辖: 量子化学

引用回帖:

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

virtualzx

木虫 (著名写手)

QC强帖: 2
应助: 263 (大学生)
金币: 7161.3
红花: 54
帖子: 1605
在线: 317.6小时
虫号: 2069080
注册: 2012-10-18
性别: GG
专业: 理论和计算化学

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +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

gmy1990

荣誉版主 (著名写手)

QC强帖: 1
应助: 163 (高中生)
贵宾: 0.782
金币: 13422.7
散金: 526
红花: 42
沙发: 4
帖子: 1627
在线: 1019.7小时
虫号: 745920
注册: 2009-04-11
专业: 理论和计算化学
管辖: 量子化学

引用回帖:

4楼: Originally posted by virtualzx at 2013-09-11 21:22:14
楼上的输入应该没问题的；你再仔细查看下是否少抄了一行。不行的话你还可以再试试guess=mix。

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

如果你算Ar2 ...

对于H2的解离，是否应该要用非限制性方法？但是由于都是基于闭壳层参考态，所以，得到的结果就是限制和非限制得到的结果都一样？

Ar的二聚体问题已经解决了，由于作用力很小，在设置合适的步长扫描后就没问题了。这个不考虑限制和非限制的问题。

5楼2013-09-11 21:38:08

virtualzx

木虫 (著名写手)

QC强帖: 2
应助: 263 (大学生)
金币: 7161.3
红花: 54
帖子: 1605
在线: 317.6小时
虫号: 2069080
注册: 2012-10-18
性别: GG
专业: 理论和计算化学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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朵）

gmy1990

6楼2013-09-12 03:28:43

gmy1990

荣誉版主 (著名写手)

QC强帖: 1
应助: 163 (高中生)
贵宾: 0.782
金币: 13422.7
散金: 526
红花: 42
沙发: 4
帖子: 1627
在线: 1019.7小时
虫号: 745920
注册: 2009-04-11
专业: 理论和计算化学
管辖: 量子化学

送红花一朵

引用回帖:

6楼: Originally posted by virtualzx at 2013-09-12 03:28:43
H2用各种方法计算出来的势能面你可以在szabo&ostlund的书上看到，uhf和rhf是在227页，ump2和ump3是在375页。

H2在解离时相关能应该是0，uhf,ump2,casscf和ccsd都应该给出精确解-1a.u.；你可以以此来查看你u ...

第一个图就是做的H2的解离，casscf和ccsd都给出了很好的结果，但是uhf和ump2得出来的结果是不对的，而且限制性和非限制性得出来的都一个样。

7楼2013-09-12 04:23:37

virtualzx

木虫 (著名写手)

QC强帖: 2
应助: 263 (大学生)
金币: 7161.3
红花: 54
帖子: 1605
在线: 317.6小时
虫号: 2069080
注册: 2012-10-18
性别: GG
专业: 理论和计算化学

引用回帖:

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