| 查看: 984 | 回复: 0 | |||
[求助]
LAMMPS中dpd模块的参数如何设置，求大家给予帮助
|
|
小弟刚刚学习DPD，使用LAMMPS中的DPD模块功能进行大分子自组装模拟，在重复文献中的模拟时，不能将文献中的参数完全对应的放进去，尤其是timestep，差了一个量级，文献中是0.02,而我只能到0.002，稍微设置大一点儿就出现bond atoms missing 或者angle atom missing这样的错误。。。不知道有没有人能帮忙改一改参数设置, 或者有什么指点的地方，感谢。文献是这篇 Martin Bates and Martin Walker,Soft Matter, 2009, 5, 346–353, Dissipative particle dynamics simulation of T- and X-shaped polyphilic molecules exhibiting honeycomb columnar phases. boundary p p p units lj atom_style full bond_style harmonic angle_style harmonic variable T equal 0.7 pair_style dpd ${T} 1.0 34387 read_data resultTshape3.data mass * 1 bond_coeff 1 4.0 0. angle_coeff 1 4.0 180 angle_coeff 2 4.0 140 angle_coeff 3 1.0 180 angle_coeff 4 4.0 120 special_bonds lj 1.0 1.0 1.0 #default is 0 0 0. variable gamma equal 3.67^2/2./${T} pair_coeff * * 0 0 pair_coeff 6 6 5 ${gamma} #6:polar; 7:core1; 9:core2; 8:side chain. Sigma=sqrt(2*T*gamma)=3.67(given in the paper) pair_coeff 6 7 15 ${gamma} #pair_coeff specifies A and gamma, do not support mixing pair_coeff 6 8 40 ${gamma} #coefficients for all I,J pairs must be specified explicitly. pair_coeff 6 9 10 ${gamma} #LAMMPS sets ij and ji to be the same values. pair_coeff 7 7 15 ${gamma} #pair_coeff command can override a previous setting. pair_coeff 7 8 40 ${gamma} pair_coeff 7 9 15 ${gamma} pair_coeff 8 8 20 ${gamma} pair_coeff 8 9 40 ${gamma} pair_coeff 9 9 10 ${gamma} communicate multi vel yes #velocites are stored by ghost atoms for dpd style. velocity all create ${T} 49863 dist uniform neighbor 2.0 bin neigh_modify delay 0 every 1 timestep 0.001 # timestep 0.02, 1M steps was found to be sufficient for N=5000 system thermo 100000 thermo_style one restart 100000 restart dump d1 all xyz 1000000 dump.fluid.*.xyz dump_modify d1 flush yes fix 1 all nvt temp 0.7 1.7 1.0 run 5000000 unfix 1 fix 1 all nvt temp 1.7 1.7 1.0 run 50000000 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : Computer_Simulation-4.pdf
2014-06-09 15:00:16, 1.77 M
» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有159人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关商家推荐: (我也要在这里推广)
5