|
|
[交流]
【求助】求吸附能,吸附热的问题
我用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) 相差十万八千里,
所以想请教一下大家。 |
» 猜你喜欢
售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
» 本主题相关价值贴推荐,对您同样有帮助:
» 抢金币啦!回帖就可以得到:
查看全部散金贴
|