| 查看: 1019 | 回复: 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 这个程序是别人拿来做另外一种研究的，跟我的差不多，它的就能运行，怎么我的不行呢 |
回复此楼» 猜你喜欢
梦想
已经有7人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
面上合作单位盖章
已经有8人回复
-
面上函评意见出来了，像什么等级？
已经有15人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
» 本主题相关商家推荐: (我也要在这里推广)
找到一些相关的精华帖子，希望有用哦~
关于优先吸附的问题，求大神指点！
已经有7人回复
- 2015年考研求助，求大神帮忙
已经有15人回复
- 各位大神求帮助，高分子的一些小问题
已经有7人回复
- 投了篇文章，竟然有四个外审，求大神帮忙分析下情况
已经有5人回复
- 求大神帮忙解答！
已经有8人回复
- 请问有没有哪位环境类的大神做分子模拟之类的
已经有5人回复
- 求高数大神 u帮忙做下 题
已经有12人回复
- 求两个值中最大的值但编写的代码不知道哪里错了，求大神帮忙。
已经有10人回复
- 求各路大神们帮忙鉴定一下下面几种植物,谢谢了
已经有5人回复
- 求积分，不会啊求大神帮帮忙！
已经有6人回复
- 求师哥师姐大神帮忙翻译一下 不胜感激！！
已经有6人回复
- amber MD过程中出现的问题 求各位大神解答
已经有7人回复
- 求做聚酰亚胺的大神，帮忙翻译一下几个专业名词
已经有9人回复
- 求各位大神帮忙下一篇文章
已经有3人回复
- 求大神帮忙研究一下关于铝合金的抗弯强度的问题！！！！！！！
已经有13人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
5