应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » Gaussian 计算出的能量有问题
121/2返回列表12下一页
查看: 1292  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

dybf0428

银虫 (初入文坛)

[求助] Gaussian 计算出的能量有问题

用 Gaussian计算CO2的分子轨道，键长为1.163埃，采用3-21G基组，计算出最高的几个分子轨道的能量为-0.5211，-0.71881，-0.72059。

但查文献，这几个轨道的电离势应该是-0.507，-0.636，-0.665.计算出的数值都偏高，即使换其他的基计算，算出来的能量仍然对不上？

尝试计算N2分子，能量也对不上，请问是什么原因？
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-11-14 21:24:10
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

scpchaos

新虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
gmy1990: 金币+2 2012-11-14 22:04:39
dybf0428: 金币+10, ★有帮助 2012-11-16 12:35:01
1.3-21G太小了。。。换大一点的
2.理论值和实验值总是会有点差别的啊（虽然在这里应该不属于这种情况）
赞一下 回复此楼
2楼2012-11-14 22:00:43
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

dybf0428

银虫 (初入文坛)
引用回帖:
2楼: Originally posted by scpchaos at 2012-11-14 22:00:43
1.3-21G太小了。。。换大一点的
2.理论值和实验值总是会有点差别的啊（虽然在这里应该不属于这种情况）

基组换成 6-311++G(2df,2p)后，算的的能量为-0.54258，-0.71086，-0.74301，仍然差别很大？

这个基组行么，还是得换更大的，应该用哪个比较合适？
赞一下 回复此楼
3楼2012-11-14 23:07:27
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

gujing136

禁虫 (文坛精英)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
dybf0428: 金币+10, ★有帮助 2012-11-16 12:35:14
一个是方法，一个是基组，都是要考虑的——你用的什么方法和什么样的基组？

准确的来说，想要把计算值算的和实验值一模一样是不大可能的，因为量化的理论都还只是处在模拟阶段——只有某种方法适合某种体系，没有通行的方法的。
赞一下 回复此楼
叫老猫的兔子
4楼2012-11-15 21:36:58
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

dybf0428

银虫 (初入文坛)
引用回帖:
4楼: Originally posted by gujing136 at 2012-11-15 21:36:58
一个是方法，一个是基组，都是要考虑的——你用的什么方法和什么样的基组？

准确的来说，想要把计算值算的和实验值一模一样是不大可能的，因为量化的理论都还只是处在模拟阶段——只有某种方法适合某种体系，没有 ...

方法用的是HF，基组试过3-21G和6-311++G(2df,2p)，能量都差别很大
赞一下 回复此楼
5楼2012-11-16 12:34:36
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

gauss98

禁虫 (正式写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
dybf0428: 金币+10, ★★★很有帮助 2012-11-17 12:44:42
本帖内容被屏蔽
6楼2012-11-16 13:06:58
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

gujing136

禁虫 (文坛精英)

【答案】应助回帖

引用回帖:
5楼: Originally posted by dybf0428 at 2012-11-16 12:34:36
方法用的是HF，基组试过3-21G和6-311++G(2df,2p)，能量都差别很大...

HF？？？？天，这个方法，只能是计算大体系，用来定性的——你这个起码都要用MP2……
赞一下 回复此楼
叫老猫的兔子
7楼2012-11-16 13:56:20
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

dybf0428

银虫 (初入文坛)
引用回帖:
6楼: Originally posted by gauss98 at 2012-11-16 13:06:58
轨道能量和电离能只是近似关系

直接算  B3LYP/6-311++G(2df,2p)

E(CO2) -E (CO2+) = 0.5082(弛豫)   0.5089（不弛豫）
与 0.507的电离势 只相差 0.0012 Hartree， 约合 0.7 kcal/mol
相当准确了

轨道能量不等于电离能么？
请问您用的是不是DFT方法，在计算E(CO2) -E (CO2+) = 0.5082(弛豫)   0.5089（不弛豫）时，E(CO2)和E(CO2)是多少，在哪里看？
赞一下 回复此楼
8楼2012-11-17 12:46:50
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

dybf0428

银虫 (初入文坛)
引用回帖:
7楼: Originally posted by gujing136 at 2012-11-16 13:56:20
HF？？？？天，这个方法，只能是计算大体系，用来定性的——你这个起码都要用MP2……...

用MP4算了，轨道能量仍然是-0.54，差别较大
赞一下 回复此楼
9楼2012-11-17 12:47:37
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

cxneng123

银虫 (小有名气)

【答案】应助回帖

★ ★
cenwanglai: 金币+2, 谢谢回复~ 2012-11-18 16:34:57
我个人觉得机组可能有一定的影响，但是泛函的影响更大的，你可以改变泛函在试试，不过你的已经很准了。像我计算吸收的时候有时都差100nm的。
赞一下 回复此楼
努力，就有希望，相信未来！
10楼2012-11-17 23:52:57
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dybf0428 的主题更新
121/2返回列表12下一页
信息提示
关闭
请填处理意见
关闭