| 查看: 3469 | 回复: 13 | ||
[求助]
解离能怎么做 已有2人参与
|
||
|
我在做一个HgX2(x是卤素)的解离能 不太清楚具体怎么用高斯来计算 已经做了分子的opt+freq 接下来是不是分成片段Hg+X2来opt+freq然后再算前后的能差呢 求各位前辈大神教教我 都困扰了好几天 网上也没看到别人具体怎么做的 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
分子动力学模拟合作需求
已经有4人回复
-
中国科学院成都山地灾害与环境研究所2026年特别研究助理（博士后）招聘启事
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有225人回复
-
湖北文理学院材料与化工专硕还有调剂名额，欢迎大家推荐考生！
已经有1人回复
-
湖南大学刘巧玲课题组2026年博士研究生招生信息
已经有28人回复
-
南方科技大学表面催化反应动力学研究课题组招收2026届博士生
已经有1人回复
-
北京邮电大学物理学院老师招收博士生
已经有3人回复
-
Comsol 6.0 桶内肥料产热产湿过程模拟求助
已经有1人回复
-
微纳hub，微纳加工的宝库
已经有0人回复
-
微纳加工里的 lift-off 工艺
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 588 (博士)
- 贵宾: 0.01
- 金币: 26696.9
- 红花: 349
- 帖子: 2029
- 在线: 1461.6小时
- 虫号: 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
- 金币: 26696.9
- 红花: 349
- 帖子: 2029
- 在线: 1461.6小时
- 虫号: 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
