| 查看: 1685 | 回复: 4 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
Lammps模拟苯酚和水在PDMS膜中的扩散，msd没有增加
|
||||
|
我做水和苯酚在聚二甲基硅氧烷膜中的扩散对比，以下是in文件中设置的部分，为什么水和苯酚的msd根本没有增加，只是稍微有波动？（本人是初学者，刚接触几个月） 体系中有3个聚二甲基硅氧烷分子，苯酚和水各五个 pair_style hybrid lj/charmm/coul/charmm 8.0 10.0 lj/cut/tip4p/long 6 7 5 8 0.15 12.0 bond_style harmonic angle_style charmm kspace_style pppm/tip4p 0.001 neighbor 2.0 bin neigh_modify delay 0 every 1 check yes timestep 1.0 min_style fire minimize 1.0e-16 1.0e-16 50000 100000 fix limit all nve fix speed all langevin 293 293 10 482799 run 500000 unfix limit fix relax all nvt temp 293 293 100 run 1000000 |
回复此楼» 猜你喜欢
生物学308分求调剂（一志愿华东师大）
已经有3人回复
-
求调剂：一志愿：南京大学 专业：0705 总分320 ，本科985，四六级已过
已经有3人回复
-
环境工程 085701，267求调剂
已经有14人回复
-
308求调剂
已经有12人回复
-
求调剂：085600材料与化工，考材科基，总分319
已经有21人回复
-
311求调剂
已经有8人回复
-
294分080500材料科学与工程求调剂
已经有13人回复
-
288求调剂 一志愿哈工大 材料与化工
已经有19人回复
-
288资源与环境专硕求调剂，不限专业，有学上就行
已经有23人回复
-
一志愿085600中科院宁波所276分求调剂
已经有17人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助LAMMPS关于MSD计算的问题
已经有17人回复
|
我是想模拟水和苯酚的布朗运动所以用了fix langevin，我对模拟的过程了解的不太多，msd计算用的是 compute m1 phenol msd/molecule compute m2 h2o msd/molecule variable msd1 equal c_m1[1][4] variable msd2 equal c_m1[2][4] variable msd3 equal c_m1[3][4] variable msd4 equal c_m1[4][4] variable msd5 equal c_m1[5][4] variable msd6 equal c_m2[1][4] variable msd7 equal c_m2[2][4] variable msd8 equal c_m2[3][4] variable msd9 equal c_m2[4][4] variable msd10 equal c_m2[5][4] variable step equal step fix msd all print 10000 "${step} ${msd1} ${msd2} ${msd3} ${msd4} ${msd5} ${msd6} ${msd7} ${msd8} ${msd9} ${msd10}" file all.msd screen no |
3楼2015-04-27 14:37:43
lsloneil
无虫
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
【答案】应助回帖
★ ★ ★ ★
感谢参与，应助指数 +1
wzjing103: 金币+2, ★有帮助 2015-04-27 16:13:19
wzjing103(abinitio代发): 金币+2, 鼓励交流 2015-04-28 07:42:17
感谢参与，应助指数 +1
wzjing103: 金币+2, ★有帮助 2015-04-27 16:13:19
wzjing103(abinitio代发): 金币+2, 鼓励交流 2015-04-28 07:42:17
|
很难说，你的设置看上去有点奇怪。正常体系的MD模拟，一般是先minimize，然后NVT下平衡，然后NVE下采样。 而你的步骤，在minimize之后，却用了fix langevin + fix nve 你这样做有什么特殊的理由吗？是别人文献里也这么做的？ fix langevin一般是对隐式溶剂(implicit solvent）下的模拟使用的，也就是说在模拟中没有水分子，所有水分子被当成均一的介质，用随机力来模拟溶质分子与水分子的碰撞。但我看你的pair_style里有 lj/cut/tip4p，看起来是有水分子存在的，所以我觉得用fix langevin不是非常合理。 另外, fix langevin不能和fix nvt混用，所以你在unfix nve后也需要unfix langevin。 你的MSD计算是怎么设置的？ |
2楼2015-04-27 13:26:14
silence-jin
版主
|
本帖内容被屏蔽 |
4楼2015-12-16 21:36:20
没落的华丽
管理员
- 应助: 3 (幼儿园)
- 金币: 586.6
- 红花: 6
- 帖子: 96
- 在线: 259.1小时
- 虫号: 2131493
- 注册: 2012-11-17
- 性别: GG
- 专业: 生化分析及生物传感
5楼2015-12-18 19:02:19
10