巨正则蒙特卡罗模拟

zhaoxin · 发表于 Post on 2016-12-17 21:32:05

关于吸附问题，很多研究工作者用巨正则蒙特卡罗模拟方法，但是我查了一些资料后，只是了解到时保持化学势，温度，体积不变，具体的很多操作什么都不明白，有没有大神可以推荐下一些教程或者网址可以学习这个模拟方法啊？

winterzen · 发表于 Post on 2016-12-17 21:54:08

可以看一下分子模拟的专著并结合软件学习
介绍MC的一本很有名的专著是Frenkel & Smit《Understanding Molecular Simulation: From Algorithms to Applications 》
GCMC常用软件：Towhee、RASPA、MUSIC、MS Sorption等

amani73986 · 发表于 Post on 2016-12-18 12:09:09

补充一下MUSIC是免费软件用的比较多（美国西北大学snurr教授课题组开发，相应的算例也比较多）、MS Sorption（你懂的，有图形化界面操作简单，就是有这个软件版权的比较少，如果资金充足可以试试超算中心，有的中心有MS版权）

ruanyang · 发表于 Post on 2016-12-18 17:07:37

可试试最新的RASPAhttp://muchong.com/html/201509/9443327.html

zhaoxin · 发表于 Post on 2016-12-20 08:47:25

谢谢大家的指导

akun · 发表于 Post on 2024-5-25 15:28:36

winterzen 发表于 2016-12-17 21:54
可以看一下分子模拟的专著并结合软件学习
介绍MC的一本很有名的专著是Frenkel & Smit《Understanding Mole ...

您好，请问您有这本书吗？有的话能麻烦您发给我下吗？祝您一切顺利!
我的邮箱：3198349661@qq.com

xiaowei2022 · 发表于 Post on 2025-5-28 19:26:20

Lammps也可以，看Lammps手册

wmg166 · 发表于 Post on 2026-4-23 14:26:35

大神您在研究巨正则公式吗？绘制过这种图吗有代码吗求分析好像就是利用巨正则公式在MATLAB里面绘制的？

yuzc · 发表于 Post on 2026-4-28 10:36:48

wmg166 发表于 2026-4-23 14:26
大神您在研究巨正则公式吗？绘制过这种图吗有代码吗求分析好像就是利用巨正则公式在MATLAB里面绘制的 ...

没有所谓的“巨正则公式”，你理解有问题。这个很明显就是对参数空间进行扫描然后画了一个3D的热图，以吸附数为染色标准。

wmg166 · 发表于 Post on 2026-4-28 15:23:06

yuzc 发表于 2026-4-28 10:36
没有所谓的“巨正则公式”，你理解有问题。这个很明显就是对参数空间进行扫描然后画了一个3D的热图，以 ...

这里的表1第一列数据代入公式8，公式8好像就是巨正则公式然后绘制的哪个三维图，我想用Python绘制绘制出来的图片和这个近似 P轴方向相反不晓得怎样调换？

yuzc · 发表于 Post on 2026-4-29 09:19:39

wmg166 发表于 2026-4-28 15:23
这里的表1第一列数据代入公式8，公式8好像就是巨正则公式然后绘制的哪个三维图，我想用Python绘制绘制 ...

这叫配分函数。每个系综种类都有对应的配分函数。它写的不是原始形式，你最好看一下出处。

这里的唯一门槛就是，你作为化学系的人，要稍微学一点平衡态统计力学，学到系综法结束。

最后说到绘图，因为凡GCMC，最后要通过H2的 EOS，将温度、压力+H2（一般要提供临界温度，临界压力，偏心因子这三个数据，通过P-R EOS将温度，压力折合成逸度）性质转换成目标化学势，所以你才有了这个热图。

算法上，你可以自己写Metropolis，对着这个配分函数做采样。相当于你有了相空间密度ρ。

总的来讲，研究细节你去看参考文献53，54就好了，看一下使用的是自主软件/脚本还是towhee，raspa等软件。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[蒙特卡罗] 巨正则蒙特卡罗模拟

评分 Rate

评分 Rate

评分 Rate

浏览过的版块