| 查看: 3339 | 回复: 5 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
关于npt系综设定？
|
|||
|
最近学习用LAMMPS做lipid bilayer的分子动力学模拟，对几个问题感到困惑，有没有做过这方面模拟的童鞋，请教一下： 1.模型构建，我看文献上的模型构建都写的不怎么详细，就自己写了段小程序，大体是这么进行的，参考文献中的每个lipid分子的平均横向面积，按照正方形晶格排一层头基珠子，然后在头基珠子上再接tail。然后复制成双层膜，形成的初始构型，这种做法做出的初始构型类似于有序的晶体结构。请教有经验的同学，这种做法是否合理？ 2.我看文献中跑双层膜的过程，其中文献中是这样说的。 The initial configuration is equilibrated to a zero-lateral-tension state using the NPT ensemble. This typically took around 104τ .The equilibrium box length for the uncharged system was found to be L ≈ 50σ, while for the most charged 1 : 3 system,the equilibrium box size slightly shrank to L ≈ 48.5σ.Upon equilibration, the system was simulated for an additional 105τ in the NVT ensemble using standard Nosé-Hoover thermostat. 这里我用npt系综平衡至a zero-lateral-tension state，这个命令遇到了麻烦 我写的命令是 fix 1 all npt temp 0.7 0.7 10.0 aniso 0.0 0.0 100.0， 如果用压强是0.0的话，系统体积在不断的膨胀，也不能形成文件中的双层膜，而是整的盒子都充满了lipid，不像文献中的那样，只在盒子中间形成了很好的双层膜，用压强是一的话，盒子总是变小，也形不成膜，想请教下，怎么设定npt系综平衡至a zero-lateral-tension state ps：参考文献 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : Sknepnek_JChemPhys_137_104905_2012.pdf
2014-08-11 17:55:45, 1.05 M
» 猜你喜欢
收到复试调剂但是去不了
已经有7人回复
-
260求调剂
已经有5人回复
-
急需调剂
已经有10人回复
-
申博/考博
已经有4人回复
-
化工学硕294分，求导师收留
已经有37人回复
-
求调剂
已经有11人回复
-
一志愿华中农业071010，320求调剂
已经有19人回复
-
304求调剂
已经有7人回复
-
求博导｜生物质基多孔碳/超级电容方向，已有相关成果，寻能源材料/碳材料方向老师
已经有3人回复
-
二苯甲酮酸类衍生物
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 碳纳米管MD模拟中积分时间步长选取时存在的问题
已经有15人回复
- 关于NPT控压器的疑问
已经有8人回复
- NPT系综下判断是否平衡的标志
已经有12人回复
- 扩散系数的计算及体系平衡的判断
已经有4人回复
- 关于NVT系统下的压力设定
已经有10人回复
- 关于lammps fix npt/deform做拉伸时温度的变化
已经有6人回复
- NPT系综动力学
已经有19人回复
- gromacs 模拟NPT系宗出问题
已经有8人回复
- gromacs NPT 系综如何模拟？
已经有8人回复
- 第一原理动力学的系综选择，和各种系综的输入文件Example
已经有13人回复
- Gibbs系综和NPT系综的区别
已经有4人回复
- nvt系综，系统总能量震荡厉害
已经有9人回复
- 【原创】EMASS是CPMD最重要的参数
已经有20人回复
- 【求助】NPT分子动力学程序
已经有5人回复
- 【求助】跑discover出现的问题
已经有5人回复
- 【求助】md跑几百步之后out文件出现NaN（无报错）
已经有11人回复
- 【讨论】系综问题
已经有32人回复
- 【讨论】关于md中系统稳定的问题
已经有10人回复
- 【讨论】高温高压的体系需要注意的环节
已经有18人回复
4楼2014-08-12 11:49:31
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
本帖内容被屏蔽 |
2楼2014-08-12 00:30:17
阿Q~~
至尊木虫 (文坛精英)
- 应助: 77 (初中生)
- 金币: 33343.6
- 红花: 107
- 帖子: 36186
- 在线: 1028.1小时
- 虫号: 3237212
- 注册: 2014-05-27
- 性别: MM
- 专业: 宇宙学
3楼2014-08-12 06:35:57
5楼2014-08-13 17:27:09
20