24小时热门版块排行榜    

查看: 3094  |  回复: 13

初衷。

新虫 (初入文坛)

[求助] 解离能怎么做 已有2人参与

我在做一个HgX2(x是卤素)的解离能    不太清楚具体怎么用高斯来计算   已经做了分子的opt+freq  接下来是不是分成片段Hg+X2来opt+freq然后再算前后的能差呢    求各位前辈大神教教我    都困扰了好几天   网上也没看到别人具体怎么做的

发自小木虫Android客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

paramecium86

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
scfslyzkf: 金币+5, 应助指数+1, 很详细,学习了 2017-04-19 15:01:23
conventional解离能可以直接opt算出complex的能量 和每一个monomer的能量 然后用楼上楼上给的公式BDE=E(complex)-E(monomer1)-E(monomer2)。

也可以考虑BSSE并且用counterpoise方法进行校正。
关键词就是counterpoise
input写法大概就是


# B3LYP/SDD Counterpoise=2 NoSymm Opt

HgX + X, optimization with counterpoise correction using ECP basis

0 1 0 2 0 2                              ! charge and multiplicity for  complex, monomer1, and monomer 2, respectively
Hg(Fragment=1)    0.00000000   0.00000000  2.000
X(Fragment=1)   0.00000000   0.00000000  1.000
X(Fragment=2)    0.00000000   0.00000000   3.000

正常算完的话找到输出文件里面类似这样的话
  Counterpoise corrected energy =    xxxxx
                  BSSE energy =       xxxxx
             sum of fragments =    xxxxxx
          complexation energy =     xxxxx kcal/mole (raw)
          complexation energy =     xxxxx kcal/mole (corrected)

基本上就能得到解离能了
4楼2017-04-19 03:13:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

含能材料

专家顾问 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
BDE=E(X-Hg-X)-E(X-Hg.)-E(.X),简单的理解就是断裂掉一根键,分子片段与整体的能量差
2楼2017-04-18 16:04:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

初衷。

新虫 (初入文坛)

谢谢   其实是能想通理论啦   就是不知道高斯要怎么做    那个关键词是用来算BDE的    要怎么做

发自小木虫Android客户端
3楼2017-04-18 20:07:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

初衷。

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by paramecium86 at 2017-04-19 03:13:43
conventional解离能可以直接opt算出complex的能量 和每一个monomer的能量 然后用楼上楼上给的公式BDE=E(complex)-E(monomer1)-E(monomer2)。

也可以考虑BSSE并且用counterpoise方法进行校正。
关键词就是counte ...

太感谢了   这个真的特详细    我还在sob的帖子里看到另一种感觉是不同的方法   先算整个分子HgX2的能量   然后减去Hg原子个X2分子的单点能    完了后两者想加减分子能量再加上热化学修正的量    这样算是不是也可以  那两种方法哪种的更可靠呢

发自小木虫Android客户端
5楼2017-04-19 09:53:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

零点校正 是应该加进去的 既然分解前后振动差距很大。做法就是 分别opt freq complex和两个monomer 得到零点校正以后的能量(zero-point corrected energy) 然后再用公式BDE=E(complex_corr)-E(mono1_corr)-E(mono2_corr)就可以了.这个和BSSE校正是不同类型的校正 BSSE校正主要是修正计算中人为加入的误差

发自小木虫IOS客户端
6楼2017-04-19 10:28:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

初衷。

新虫 (初入文坛)

引用回帖:
6楼: Originally posted by paramecium86 at 2017-04-19 10:28:55
零点校正 是应该加进去的 既然分解前后振动差距很大。做法就是 分别opt freq complex和两个monomer 得到零点校正以后的能量(zero-point corrected energy) 然后再用公式BDE=E(complex_corr)-E(mono1_corr)-E(mono2_ ...

明白了   十分感谢   我先试试算一下   看看结果怎么样再来请教

发自小木虫Android客户端
7楼2017-04-19 10:41:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scfslyzkf

木虫 (著名写手)

引用回帖:
2楼: Originally posted by 含能材料 at 2017-04-18 20:04:05
BDE=E(X-Hg-X)-E(X-Hg.)-E(.X),简单的理解就是断裂掉一根键,分子片段与整体的能量差

请问BDE的全称是?
8楼2017-04-19 15:00:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scfslyzkf

木虫 (著名写手)

引用回帖:
4楼: Originally posted by paramecium86 at 2017-04-19 07:13:43
conventional解离能可以直接opt算出complex的能量 和每一个monomer的能量 然后用楼上楼上给的公式BDE=E(complex)-E(monomer1)-E(monomer2)。

也可以考虑BSSE并且用counterpoise方法进行校正。
关键词就是counte ...

请问BSSE指的是?
9楼2017-04-19 15:01:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scfslyzkf

木虫 (著名写手)

引用回帖:
4楼: Originally posted by paramecium86 at 2017-04-19 07:13:43
conventional解离能可以直接opt算出complex的能量 和每一个monomer的能量 然后用楼上楼上给的公式BDE=E(complex)-E(monomer1)-E(monomer2)。

也可以考虑BSSE并且用counterpoise方法进行校正。
关键词就是counte ...

请问BSSE指的是?
10楼2017-04-19 15:01:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 初衷。 的主题更新
信息提示
请填处理意见