| 查看: 985 | 回复: 0 | ||
[求助]
刚做分子模拟的研究，求大神帮忙下！
|
|
下面是我做的计算乙炔在AL2O3孔道中的扩散系数的源文件，但是用LAMMPS运行的时候老是出现LOST ATOMS，能帮忙看看是哪里出问题了么 # description c2h2 diffusion in AL2O3 # used as template for the single atom molecules diffusion in AL2O3. variable temperature equal 298.0 variable dumpstep equal 10000 variable thermostep equal 10000 variable comstep equal 10000 variable restartstep equal 10000000 variable runrelaxstep equal 1000000 variable runstep equal 20000000 units real atom_style atomic read_data AL2O3.c2h2.data.lammps05 pair_style lj/cut 17.0 #pair_coeff * * none # parameter from Dreiding for materials: generally the type equals to the atom number. 4-Zinc # 2-He from Peng X, 3-Li from UFF. pair_coeff 1 1 0.5050 2.2495 pair_coeff 2 2 0.2104 1.7210 #parameter from Peng Xuan for species: 9-Methane, 10-Hydrogen2, 11-Nitrogen2, 12-Oxygen2 pair_coeff 3 3 0.1588 3.5742 pair_modify shift yes mix arithmetic label loop variable a loop 2 group fluid type 3 group AL2O3 subtract all fluid # fix freeze AL2O3 setforce 0.0 0.0 0.0 velocity fluid create ${temperature} 3654321 dist gaussian neighbor 2.0 bin neigh_modify exclude group AL2O3 AL2O3 timestep 0.1 # attention: for the new lammps exe, the fix nvt changed as fix 1 fluid nvt temp 298.0 298.0 100.0 #fix 1 fluid nvt temp ${temperature} ${temperature} 100.0 #fix 1 fluid nvt ${temperature} ${temperature} 100.0 run ${runrelaxstep} reset_timestep 0 # in the new lammps exe file, this two fix have been discarded. fix 2 fluid msd 10000 AL2O3.c2h2.diff.out #fix 3 fluid com 10000 AL2O3.c2h2.com compute fluidtem fluid temp #compute fluidpress all pressure fluidtem #compute peratom fluid pe/atom #compute pe fluid reduce sum c_peratom #dump 1 fluid custom 100 dump.AL2O3.c2h2.fluid.wrap.$t id type x y z ix iy iz dump 2 fluid custom 1000 dump.AL2O3.c2h2.fluid.unwrap.$a id type xu yu zu #dump 3 all custom 1000 dump.AL2O3.c2h2.all id type x y z ix iy iz #dump 4 fluid custom 100 dump.AL2O3.c2h2.fluid.v.$t id type vx vy vz # dump_modify id scale no image yes thermo ${thermostep} thermo_style custom step temp c_fluidtem press pe etotal #thermo_style custom step temp etotal press pe c_pe # thermo_modify temp fluidtem press fluidpress restart ${restartstep} AL2O3.c2h2.restart.$a run ${runstep} clear # use the restart file as the initial configuration of the md every loop # because in the velocity command, the system has been set a new velocity # so it may be independent with the former configuration run, i think. read_restart AL2O3.c2h2.restart.$a.${runstep} next a jump in.lammps loop 这个程序是别人拿来做另外一种研究的，跟我的差不多，它的就能运行，怎么我的不行呢 |
回复此楼» 猜你喜欢
某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！
已经有5人回复
-
26/27申博自荐
已经有9人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有7人回复
-
河北省自然科学基金
已经有7人回复
-
揭秘青基评审内幕：几个A才能顺利中标
已经有4人回复
-
青B发送上会通知了吗
已经有7人回复
-
博士申请
已经有3人回复
-
投稿求助，期刊
已经有8人回复
-
献血感触
已经有13人回复
» 本主题相关商家推荐: (我也要在这里推广)
找到一些相关的精华帖子，希望有用哦~
关于优先吸附的问题，求大神指点！
已经有7人回复
- 2015年考研求助，求大神帮忙
已经有15人回复
- 各位大神求帮助，高分子的一些小问题
已经有7人回复
- 投了篇文章，竟然有四个外审，求大神帮忙分析下情况
已经有5人回复
- 求大神帮忙解答！
已经有8人回复
- 请问有没有哪位环境类的大神做分子模拟之类的
已经有5人回复
- 求高数大神 u帮忙做下 题
已经有12人回复
- 求两个值中最大的值但编写的代码不知道哪里错了，求大神帮忙。
已经有10人回复
- 求各路大神们帮忙鉴定一下下面几种植物,谢谢了
已经有5人回复
- 求积分，不会啊求大神帮帮忙！
已经有6人回复
- 求师哥师姐大神帮忙翻译一下 不胜感激！！
已经有6人回复
- amber MD过程中出现的问题 求各位大神解答
已经有7人回复
- 求做聚酰亚胺的大神，帮忙翻译一下几个专业名词
已经有9人回复
- 求各位大神帮忙下一篇文章
已经有3人回复
- 求大神帮忙研究一下关于铝合金的抗弯强度的问题！！！！！！！
已经有13人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
10