| 查看: 3184 | 回复: 13 | ||
[求助]
解离能怎么做 已有2人参与
|
||
|
我在做一个HgX2(x是卤素)的解离能 不太清楚具体怎么用高斯来计算 已经做了分子的opt+freq 接下来是不是分成片段Hg+X2来opt+freq然后再算前后的能差呢 求各位前辈大神教教我 都困扰了好几天 网上也没看到别人具体怎么做的 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
国内青年破解费曼难题，从分子层面阐释液体蒸发（气液相变）微观过程
已经有5人回复
-
研究发现一个新的临界系数，独立于临界压缩因子
已经有7人回复
-
物理化学论文润色/翻译怎么收费?
已经有255人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
已知面上挂了，小学校实验方向无平台，感觉太难了
已经有26人回复
-
求助各位大佬，球球了，在这科萨头
已经有0人回复
-
七嗪类物质合成求助
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 588 (博士)
- 贵宾: 0.01
- 金币: 26423.6
- 红花: 349
- 帖子: 2026
- 在线: 1442.8小时
- 虫号: 2466443
- 注册: 2013-05-16
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
scfslyzkf: 金币+5, 应助指数+1, 很详细，学习了 2017-04-19 15:01:23
感谢参与，应助指数 +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
2楼2017-04-18 16:04:05
3楼2017-04-18 20:07:49
|
太感谢了 这个真的特详细 我还在sob的帖子里看到另一种感觉是不同的方法 先算整个分子HgX2的能量 然后减去Hg原子个X2分子的单点能 完了后两者想加减分子能量再加上热化学修正的量 这样算是不是也可以 那两种方法哪种的更可靠呢 发自小木虫Android客户端 |
5楼2017-04-19 09:53:52
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 588 (博士)
- 贵宾: 0.01
- 金币: 26423.6
- 红花: 349
- 帖子: 2026
- 在线: 1442.8小时
- 虫号: 2466443
- 注册: 2013-05-16
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
|
零点校正 是应该加进去的 既然分解前后振动差距很大。做法就是 分别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
7楼2017-04-19 10:41:33
scfslyzkf
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 3827.2
- 散金: 957
- 红花: 12
- 帖子: 1239
- 在线: 296.2小时
- 虫号: 2853150
- 注册: 2013-12-06
- 性别: GG
- 专业: 理论和计算化学
8楼2017-04-19 15:00:32
scfslyzkf
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 3827.2
- 散金: 957
- 红花: 12
- 帖子: 1239
- 在线: 296.2小时
- 虫号: 2853150
- 注册: 2013-12-06
- 性别: GG
- 专业: 理论和计算化学
9楼2017-04-19 15:01:39
scfslyzkf
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 3827.2
- 散金: 957
- 红花: 12
- 帖子: 1239
- 在线: 296.2小时
- 虫号: 2853150
- 注册: 2013-12-06
- 性别: GG
- 专业: 理论和计算化学
10楼2017-04-19 15:01:42
20