| 查看: 930 | 回复: 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
» 猜你喜欢
德国亥姆霍兹Hereon中心招收两位医用镁合金腐蚀与LPSO相变方向2026公派博士生
已经有0人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，推荐试试
已经有11人回复
-
物理学I论文润色/翻译怎么收费?
已经有259人回复
-
2026-CJ开始申报了
已经有1人回复
-
西安电子科学大学杭州研究院刘丽香教授招收智能多模态传感器和微型储能器件方向博士
已经有12人回复
-
重庆交通大学光子学微结构与器件课题组招收2026年硕士研究生信息
已经有0人回复
-
一志愿郑大材料学硕298分，求调剂
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
10