查看: 3288  |  回复: 6

hukd

银虫 (小有名气)


[交流] 【求助】求吸附能,吸附热的问题

我用sorption做丁烷吸附在MFI的模拟。想求吸附热。我看了文献和论坛上交流去做但结果不对。下面是我做的过程:

1.   import 丁烷,用 discover 优化,并记下优化好的能量 -5.52kcal/mol

2.  import MFI,  做 P1,  用discover 优化,并记下优化好的能量-5697.72 (单位它没写清楚但我推测一般discover能量是用kcal/mol)

3.  用sorption 做fixed locating, loading 分子数为1, 温度 300K, 力场 PCFF,其他基本按默认值,平衡500000步,产出也用了 500000步。

4.  吸附完后,我把轨迹文件中的lowest energy的3D 文件打开,并用discover计算了能量(我没有再对其优化),得到能量为-1949.51kcal/mol

5. 根据吸附能=吸附质能量+吸附体能量-吸附复合体能量
=-5.52 kcal/mol+(-5697.72 kcal/mol )- (-1949.51 kcal/mol) =3753.73kcal/mol , 和文献中实验以及其他模拟的值(50 kj/mol) 相差十万八千里,

所以想请教一下大家。
回复此楼

» 猜你喜欢

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

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hukd

银虫 (小有名气)


没有人知道吗? 请高手现身
2楼2011-03-28 16:20:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxrhuanan

金虫 (著名写手)


★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
franch(金币+1): 谢谢回帖交流 2011-06-06 21:02:17
你这样子算差的肯定很远,看out文件
3楼2011-06-06 20:39:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kmw.8668

金虫 (正式写手)


★
小木虫(金币+0.5):给个红包,谢谢回帖
是不是计算所用的软件不一样?
4楼2011-06-07 10:29:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

charityqi

木虫 (小有名气)


★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
zzy870720z(金币+1): 鼓励交流 2011-06-22 10:04:19
你是要算化学吸附吗?印象里adsorption得到的物理吸附的,所以分子距离离得都很远
5楼2011-06-21 22:10:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qi-cheng

铜虫 (小有名气)


★
小木虫(金币+0.5):给个红包,谢谢回帖
引用回帖:
1楼: Originally posted by hukd at 2011-03-28 11:06:23:
我用sorption做丁烷吸附在MFI的模拟。想求吸附热。我看了文献和论坛上交流去做但结果不对。下面是我做的过程:

1.   import 丁烷,用 discover 优化,并记下优化好的能量 -5.52kcal/mol

2.  import MFI,  做 ...

我想请教一下,怎么把实验值和用MS得到的计算值进行对比?
比如说用DMOL一般是算一个分子在一个面上的吸附,而实验一般是一定量的气体在一定量的吸附剂上的吸附,这样两个得到的结果有可比性吗?我不知道怎么把实验和计算结合起来,请教。。。
6楼2011-09-08 13:58:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

026li

金虫 (小有名气)


★
小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
15433楼: Originally posted by hukd at 2011-03-28 11:06:23
我用sorption做丁烷吸附在MFI的模拟。想求吸附热。我看了文献和论坛上交流去做但结果不对。下面是我做的过程:

1.   import 丁烷,用 discover 优化,并记下优化好的能量 -5.52kcal/mol

2.  import MFI,  做  ...

你想得到吸附能,用量化的方法计算更精确,主要是化学吸附,也可能得到物理吸附,比如形成氢键的弱相互作用。你用MS作的是物理吸附,这和实验文献会差很多
7楼2012-11-30 15:04:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hukd 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 12/600 2026-10-13 20:40 by Merryya
[硕博家园] 极限拉扯,如何抉择 +9 Isak0320 2026-10-12 13/650 2026-10-13 11:57 by Isak0320
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见