24小时热门版块排行榜    

查看: 3053  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

icemantoby

木虫 (初入文坛)

[交流] 【讨论】GCMC系综中化学势如何确定 已有6人参与

向大家请教一下:
我用GCMC系综模拟水物体表面的吸附情况,也就是UVT系综,那么体系的化学势是通过什么方法确定的。
我目前了解了几种方法:
第一种方法: 随意假定一个化学势,进行GCMC的模拟,就可以得到体系的压强,同样地得到了水分子的个数,假定一系列的化学势做一些列的模拟,得到压强和水分子个数,做吸附等温线。
存在的疑问:由于是液态水,假定化学势后计算得到的压强是非常大的,有几千甚至上万个bar,所以不知道这种方法是不是可靠,而且也无法得到特定压强下的水分子个数。

第二种方法:听说有这样的公式,可以通过知道压强,而计算得到水或整个体系的化学势,如果这样的话就能得到特定压强下的水分子个数。
存在的疑问:是不是存在这样一个公式?是什么样子的?

第三种方法:在文献中看到有Widom方法来计算化学势。
存在的疑问:我对Widom本身不太了解,另外我用的是MCCCS Towhee来做MC的,这个软件里面基本没有介绍Widom方法的使用,而且好像也没有相关参数的设置。

一般你们都用哪种方法来确定体系的化学势?上述三种方法是不是都可行,还有没有其他的方法可以确定体系的化学势,最好是针对水的。
Thank you in advance!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

icemantoby

木虫 (初入文坛)


zh1987hs(金币+1):谢谢参与讨论 2010-09-23 19:00:03
引用回帖:
Originally posted by ghcacj at 2010-09-22 19:16:23:

不是这样的,我记得是先将固定的分子放入一个空盒里,跑MC,会计算出一个单位为K的化学势,这个就是第二步计算吸附时的化学势输入量(此时已经将吸附剂加入盒子)。

你是指固定的分子数码? 也就是先跑NVT系综的MC,此时能得到体系的化学势,然后再将这个化学势作为输入量跑UVT系综的MC。可是第一步的分子数的值怎么确定?因为这个值必定跟体系的压强有关。
5楼2010-09-23 09:28:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 icemantoby 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见