应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » MS » 分子模拟中MSD的问题讨论
261/1返回列表
查看: 7139  |  回复: 25
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zhl1987091

木虫 (小有名气)

[交流] 分子模拟中MSD的问题讨论

首先我先交代一下背景哈：我现在用Materials Studio 做分子模拟，研究的主要是高压下的膜气体分离的问题；
然后说一下细节：在做气体在膜中的扩散模拟的时候，需要用到（均方位移曲线）MSD，大多数的文献中对于得到MSD的模拟都是用的Forcite中的NVT动力学模拟再分析得到；
问题来了：
1、但是因为我的模拟需要考虑高压的问题，NVT并没有压力的设定啊，是不是可以认为压力为0呢？
2、请问我可以用NPT来模拟这个过程吗？如果用NPT模拟的话，膜晶格的密度肯定会变化啊（现实中膜材料的密度肯定是不变的）
这个问题实在是让我纠结啊！
3、如果都不行的话，这个压力的因素怎样和扩散联系起来呢？
4、计算自由体积时，分子探针半径大小该怎么选择呢？


上次发问题好多回帖的是来蹭金币的，这次改为手动发放金币了，希望有明白的同仁，能给指点一二，只要问题能解决，金币我不在乎，可以无限追加！

[ Last edited by zhl1987091 on 2014-1-6 at 09:50 ]
回复此楼

» 收录本帖的淘帖专辑推荐

考博

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-01-02 09:07:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

笑笑愁

金虫 (小有名气)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
月只蓝: 金币+1, 鼓励交流! 2014-01-02 10:43:30
zhl1987091: 金币+10, 本来是想给30个金币的，无奈最多只能给10个 2014-01-02 15:34:07
做扩散肯定是要NVT系综的，楼主的三个问题我并在一起回答。
文献里的经常有扩散系数与压力的关系曲线，实际上是先做一定压力下的GCMC，然后把得到的构型跑MD，原来的一定压力就是你扩散MSD对应的压力哦；我们组做了很多这方面的。希望对你有帮助。
赞一下(1人) 回复此楼
3楼2014-01-02 09:38:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

笑笑愁

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
4楼: Originally posted by zhl1987091 at 2014-01-02 15:33:01
在前面做吸附溶解模拟的过程中有用到这个GCMC，比如说我做扩散要用到的压力是50bar的话，是不是就是要先做一下在50bar的压力下，晶格的溶解吸附能力，然后再利用得到的这个含有气体分子的模型，进而模拟扩散过程， ...

对的，就是这样滴
回帖评分哪里可以多次给金币的吧
赞一下(1人) 回复此楼
5楼2014-01-02 15:47:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhl1987091

木虫 (小有名气)
引用回帖:
3楼: Originally posted by 笑笑愁 at 2014-01-02 09:38:16
做扩散肯定是要NVT系综的，楼主的三个问题我并在一起回答。
文献里的经常有扩散系数与压力的关系曲线，实际上是先做一定压力下的GCMC，然后把得到的构型跑MD，原来的一定压力就是你扩散MSD对应的压力哦；我们组做了 ...

在前面做吸附溶解模拟的过程中有用到这个GCMC，比如说我做扩散要用到的压力是50bar的话，是不是就是要先做一下在50bar的压力下，晶格的溶解吸附能力，然后再利用得到的这个含有气体分子的模型，进而模拟扩散过程，就算是考虑到压力的影响了呢？
赞一下 回复此楼
4楼2014-01-02 15:33:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhl1987091

木虫 (小有名气)
引用回帖:
5楼: Originally posted by 笑笑愁 at 2014-01-02 15:47:27
对的，就是这样滴
回帖评分哪里可以多次给金币的吧...

我先试一下，但是呢，高压下的吸附量是比较大的啊，这样是不是会很费机时，有没有可以取代的方法？（金币好像不能超过10个啊，没关系，如果成功了，我可以买成红花送给你啊）我可能还有后续的问题，还要多多帮助啊！
赞一下 回复此楼
6楼2014-01-02 16:06:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

笑笑愁

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
6楼: Originally posted by zhl1987091 at 2014-01-02 16:06:00
我先试一下，但是呢，高压下的吸附量是比较大的啊，这样是不是会很费机时，有没有可以取代的方法？（金币好像不能超过10个啊，没关系，如果成功了，我可以买成红花送给你啊）我可能还有后续的问题，还要多多帮助啊 ...

MD需要动力学积分，很耗机时。除非减小盒子大小，没有好方法
赞一下 回复此楼
7楼2014-01-02 16:08:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhl1987091

木虫 (小有名气)
引用回帖:
7楼: Originally posted by 笑笑愁 at 2014-01-02 16:08:25
MD需要动力学积分，很耗机时。除非减小盒子大小，没有好方法...

还有个问题就是对于MSD来说，一个晶格中吸附溶解的气体分子数目是否会影响到曲线值的大小？
赞一下 回复此楼
8楼2014-01-02 16:23:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

笑笑愁

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
8楼: Originally posted by zhl1987091 at 2014-01-02 16:23:23
还有个问题就是对于MSD来说，一个晶格中吸附溶解的气体分子数目是否会影响到曲线值的大小？...

肯定影响啊，这个就是压力与扩散系数的关系，可以做一条曲线，文献里经常有的
赞一下 回复此楼
9楼2014-01-02 16:27:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhl1987091

木虫 (小有名气)
引用回帖:
9楼: Originally posted by 笑笑愁 at 2014-01-02 16:27:22
肯定影响啊，这个就是压力与扩散系数的关系，可以做一条曲线，文献里经常有的...

一开始你们做的那个GCMC是怎么把气体小分子放进去的呢？我试了一下，得到的晶格中只有吸附位，并没有相应的气体分子被填充进去啊？你能否说一下怎么把相应压力下的气体分子放到晶格中去的呢？
赞一下 回复此楼
10楼2014-01-02 16:32:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

笑笑愁

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
10楼: Originally posted by zhl1987091 at 2014-01-02 16:32:19
一开始你们做的那个GCMC是怎么把气体小分子放进去的呢？我试了一下，得到的晶格中只有吸附位，并没有相应的气体分子被填充进去啊？你能否说一下怎么把相应压力下的气体分子放到晶格中去的呢？...

我们是用软件MUSIC做的，分子文件做好文件会自动添加，MS怎么做我不知道
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

zhl1987091
11楼2014-01-02 16:37:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhl1987091

木虫 (小有名气)
送红花一朵
引用回帖:
11楼: Originally posted by 笑笑愁 at 2014-01-02 16:37:33
我们是用软件MUSIC做的，分子文件做好文件会自动添加，MS怎么做我不知道...

好吧！还是谢谢你，我再试试别的办法。
赞一下 回复此楼
12楼2014-01-02 16:40:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mzhyan

至尊木虫 (文坛精英)
blessing
回复此楼
14楼2014-01-02 17:26:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mlanqiang

木虫之王 (文学泰斗)
blessing
回复此楼
15楼2014-01-02 17:37:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lusostty

金虫 (小有名气)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
zhl1987091: 金币+2 2014-01-03 15:26:13
引用回帖:
10楼: Originally posted by zhl1987091 at 2014-01-02 16:32:19
一开始你们做的那个GCMC是怎么把气体小分子放进去的呢？我试了一下，得到的晶格中只有吸附位，并没有相应的气体分子被填充进去啊？你能否说一下怎么把相应压力下的气体分子放到晶格中去的呢？...

你做定压吸附 就能得到含有小分子的结构了
赞一下 回复此楼
16楼2014-01-03 10:15:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhl1987091

木虫 (小有名气)
引用回帖:
16楼: Originally posted by lusostty at 2014-01-03 10:15:58
你做定压吸附 就能得到含有小分子的结构了...

我用fox loading做的 ，然后就得到了
赞一下 回复此楼
17楼2014-01-03 15:25:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhl1987091

木虫 (小有名气)
引用回帖:
16楼: Originally posted by lusostty at 2014-01-03 10:15:58
你做定压吸附 就能得到含有小分子的结构了...

你说的应该是用fix pressure ，这个我没找到，我先做定压吸附，然后得到要吸附溶解的气体个数，然后fix loading 进去相应的个数，你说我的这个过程和你说的等价吗？
赞一下 回复此楼
18楼2014-01-03 15:31:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lusostty

金虫 (小有名气)
★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
zhl1987091: 金币+5 2014-01-04 10:29:17
引用回帖:
18楼: Originally posted by zhl1987091 at 2014-01-03 15:31:28
你说的应该是用fix pressure ，这个我没找到，我先做定压吸附，然后得到要吸附溶解的气体个数，然后fix loading 进去相应的个数，你说我的这个过程和你说的等价吗？...

loading主要考察吸附位点 而你主要目的是得到吸附稳定构型 没必要多此一举吧
赞一下 回复此楼
19楼2014-01-03 21:20:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhl1987091

木虫 (小有名气)
引用回帖:
19楼: Originally posted by lusostty at 2014-01-03 21:20:39
loading主要考察吸附位点 而你主要目的是得到吸附稳定构型 没必要多此一举吧...

哦，知道了，还有一个问题，就是在及算自由体积的时候，用到分子探针，这个探针的半径大小是根据什么来定的呢？一般多大？
赞一下 回复此楼
20楼2014-01-04 10:26:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lusostty

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
20楼: Originally posted by zhl1987091 at 2014-01-04 10:26:11
哦，知道了，还有一个问题，就是在及算自由体积的时候，用到分子探针，这个探针的半径大小是根据什么来定的呢？一般多大？...

每个探针分子大小都不一样 这个在文献里有的 我记得N2是1.8吧 记不清了 你自己查阅下文献吧
赞一下 回复此楼
21楼2014-01-04 11:30:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhl1987091

木虫 (小有名气)
引用回帖:
21楼: Originally posted by lusostty at 2014-01-04 11:30:01
每个探针分子大小都不一样 这个在文献里有的 我记得N2是1.8吧 记不清了 你自己查阅下文献吧...

比如说我要模拟的是N2的扩散，那么在算自由体积的时候就要用氮气做探针分子吗？如果是这样的话，算两种气体一起扩散时，用谁做探针分子呢？
赞一下 回复此楼
22楼2014-01-04 11:37:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lusostty

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
22楼: Originally posted by zhl1987091 at 2014-01-04 11:37:41
比如说我要模拟的是N2的扩散，那么在算自由体积的时候就要用氮气做探针分子吗？如果是这样的话，算两种气体一起扩散时，用谁做探针分子呢？...

这个问题还真不知道 不同的探针分子测量的体积是不一样的 你只能查查文献了
赞一下 回复此楼
24楼2014-01-05 12:30:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Awcse

新虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
3楼: Originally posted by 笑笑愁 at 2014-01-02 09:38:16
做扩散肯定是要NVT系综的，楼主的三个问题我并在一起回答。
文献里的经常有扩散系数与压力的关系曲线，实际上是先做一定压力下的GCMC，然后把得到的构型跑MD，原来的一定压力就是你扩散MSD对应的压力哦；我们组做了 ...

你好，想算MSD，请问我用VASP中MD跑完之后如何处理数据啊？如果方便能留下联系方式么，想向你请教一下
赞一下 回复此楼
25楼2015-11-29 16:54:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

elfyu

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
25楼: Originally posted by Awcse at 2015-11-29 16:54:29
你好，想算MSD，请问我用VASP中MD跑完之后如何处理数据啊？如果方便能留下联系方式么，想向你请教一下...

您好，我也是用VASP，不知你的问题解决了吗？能否告诉我？不胜感激
赞一下 回复此楼
26楼2019-10-26 19:59:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
sjmr12212楼
2014-01-02 09:26   回复  
ohcesc13楼
2014-01-02 17:06   回复  
12
lusostty23楼
2014-01-05 12:28   回复  
引用回帖:
22楼: Originally posted by zhl1987091 at 2014-01-04 11:37:41 比如说我要模拟的是N2的扩散，那么在算自由体积的时候就要用氮气做探针分子吗？如果是这样的话，算两种气体一起扩散时，用谁做探针分子呢？...

相关版块跳转 我要订阅楼主 zhl1987091 的主题更新
261/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么到现在没收到通知？ +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 29号明天会评吗 +3 笨笨唐 2026-08-28 3/150 2026-08-30 09:29 by cww8181
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金未中，这种答复是模板吗？ +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
关闭
请填处理意见
关闭