24小时热门版块排行榜    

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

敏sunflower

金虫 (小有名气)

[求助] 计算扩散系数 已有2人参与

初学模拟,急需各位大侠的帮助!
我最近在模拟A、B两分子在聚合物中的扩散系数,但是遇到了几个问题,想向各位请教一下
模拟过程如下:1.建Cell——2.能量最小化——3.NPT——4.NVT平衡结构——5.NVT用于数据分析
(1)做完第5步NVT动力学,分析得到MSD-t曲线,从而得到扩散系数。但是发现对第4步结束时的最后一帧结构平行做两次NVT动力学,分析得到的扩散系数有时相差很大,这是什么原因那?
(2)有人说是扩散系数对初始结构很敏感,但重新建模之后进行上述1-5步模拟,发现得到的扩散系数也相差很大
面对上述情况,怎样取扩散系数的值那?
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

敏sunflower

金虫 (小有名气)

引用回帖:
11楼: Originally posted by wangfei天空 at 2014-11-23 16:39:40
您好   请问您计算扩散系数用的是forcite还是discover呢   我已经用sorption做完了两种分子AB在孔隙中的吸附  下面要进行AB各自扩散系数的计算   现在正在看各个帖子学习怎样计算扩散  有几个问题希望得到您的帮助
...

你好!由于我做的是聚合物吸附带苯环的较大分子,聚合物密度在1g/cm3左右,此时,用sorption模块吸附不上去,就没有做吸附过程。直接做的扩散,扩散模型中A、B分子的数量,是根据聚合物分子对A、B的饱和吸附量确定的。
计算扩散系数的话,forcite和discover都可以;
至于你提的第二个问题,我没有看懂;
第三个问题,我做NPT是因为建立初始模型的密度比实验值低很多,所以要通过NPT来增大模型密度。当密度达平衡后,还有通过NVT进行300PS的弛豫,此时能量已达到平衡,这样就进一步用NVT,相当于收集计算扩散系数的数据;
第四个问题,我的服务器24核,使用23核进行NVT的话,总原子数7000左右,需要接近80小时;
13楼2014-11-24 09:23:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 31 个回答

smutao

禁虫 (著名写手)

感谢参与,应助指数 +1
本帖内容被屏蔽

2楼2014-10-01 07:14:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

smutao

禁虫 (著名写手)

★ ★ ★ ★ ★
敏sunflower: 金币+5 2014-10-01 09:18:18
本帖内容被屏蔽

3楼2014-10-01 07:22:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

敏sunflower

金虫 (小有名气)

引用回帖:
3楼: Originally posted by smutao at 2014-10-01 07:22:17
请见此处
http://muchong.com/bbs/viewthread.php?tid=7330890&fpage=1

希望对你有用

之前看过的,立场参数的话,MS软件无法修改,至于动力学时间,我是用了300psNVT来平衡体系,再进行600psNVT收集数据,从输出的能量图来看,已经达到平衡,这样应该就可以排除动力学时间太短的问题了吧?
4楼2014-10-01 09:21:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见