| 查看: 955 | 回复: 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
» 猜你喜欢
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有141人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有1人回复
-
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有5人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
-
如何从铁电相到顺电相。
已经有1人回复
» 本主题相关商家推荐: (我也要在这里推广)
