24小时热门版块排行榜    

查看: 6658  |  回复: 101
本帖产生 3 个 模拟EPI ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

qphll

金虫 (正式写手)

[交流] 【求助】CuBTC MOF + CO2 MuSiC 问题 已有12人参与

在用MuSiC测试体系: CO2+CuBTC的时候, 发现了两个问题:
(1) 采用与文献相同模拟, 相同参数, 相同软件的情况下, 吸附量差别很大;
(2) 改变压力做吸附等温线, 不同压力下的吸附量竟然差别不大.

应该是哪里的设置有些问题. 让我细细道来, 大家评论一下.

拿来测试的文献是: Molecular Simulations and Experimental Studies of CO2, CO and N2 adsorption in Metal-Organic Frameworks, J. Phys. Chem. C, 114, 15735, 2010

更具体一点, 试图重复的数据是该文献Figure-2中, 298K下CO2在CuBTC中的吸附. 在GetData的帮助下, Figure-2中CO2在CuBTC的吸附量大致如下:

Pressure, KPa              Loading, mol/kg
41.28                          2.22
120.67                        5.40
321.05                       10.93
530.09                       12.94
720.96                       14.03
1149.34                      14.95

以41.28KPa这个压力点为例, 贴出所有我计算中用到的文件.
回复此楼
Life, Love, Laugh.
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)

★ ★ ★
zh1987hs(金币+3):谢谢 2010-12-05 17:22:54
update,

前1E7扔掉, 分析了后1E7的能量变化. 看起来, 在1.5E7以后, 体系差不多到达平衡态啦.

Life, Love, Laugh.
84楼2010-12-05 13:19:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)

★ ★ ★ ★
zh1987hs(金币+4):谢谢 2010-12-05 17:23:05
引用回帖:
Originally posted by gjh123 at 2010-12-03 23:05:45:
问楼主个问题:当考虑库仑相互作用时,ALPHA@0.2  这个参数是什么意思啊?谢谢楼主

回复晚啦, 望见谅.

你的这个问题, 我前面有讲过:

KAPPA 是和 EWALD SUM里面的Gaussian function有关的一个参数, 它定义了 Gaussian Function的带宽.
如果你看Frenkel&Smit 的第十二章, 这里程序中的KAPPA参数值等于Frenkel&Smit Chap.12中的alpha^0.5.

请参见 本帖 69楼 我的某次update, 原则上讲, 不怎么需要担心这个参数. 但是弄明白些总是好的.
Life, Love, Laugh.
85楼2010-12-05 13:23:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)

★ ★
zh1987hs(金币+2):谢谢 2011-01-06 22:50:06
引用回帖:
Originally posted by ghcacj at 2010-12-11 10:45:22:
问LZ一个问题,目前文献中计算吸附等温线,一般分两种,一种是直接用music之类的MC软件跑等温线,也就是说在每个压力点下,吸附剂起始状态都是空的,然后根据uVT计算吸附量。另外一种方法是除了第一个压力点用空吸 ...

正解!

其实两种方法都是正确的. 只要每种方法的精度都达到要求, 那么应该是给出相同的结果. 其实这也是为什么MC可以续算的原因.

另外, 我推荐第二种方法. 确实能省不少时间.
Life, Love, Laugh.
88楼2010-12-11 11:28:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)

多谢版主帮忙回复, 才过完年就跑出去了. 没有时间来做更新和答复, 大家多多包涵.
Life, Love, Laugh.
95楼2011-01-05 11:11:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qphll 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见