| 查看: 7417 | 回复: 5 | |||
[求助]
求氮化硼热导率in文件出错 已有2人参与
|
|
因为是新人金币不多，希望好心人帮我看看lammps的in文件哪里出错了。我想求氮化硼的热导率，但是运行结果是这样的： ERROR on proc 0: Non-numeric atom coords - simulation unstable (../domain.cpp:510) 以下是我的in文件： units real dimension 3 boundary p p p atom_style full #read topology read_data bnsystem.data #Force-field parameters #B=1,N=1 pair_style lj/cut 2.5 bond_style harmonic angle_style harmonic dihedral_style charmm special_bonds lj/coul 0.0 0.0 0.0 pair_coeff * * 1.0 1.0 2.5 bond_coeff * 80.0 1.2 angle_coeff * 300.0 107.0 dihedral_coeff * 3.100 2 180 0.0 #NEMD calculate conductivity region box block -5 100 -1 100 -1 31.15 units box region begin block INF INF INF INF -1 7 units box region hot block INF INF INF INF 7 14 units box region mid block INF INF INF INF 14 21 units box region cold block INF INF INF INF 21 27 units box region end block INF INF INF INF 27 34 units box region boundary union 2 begin end units box group begin region begin group hot region hot group mid region mid group cold region cold group end region end group boundary region boundary group new subtract all boundary velocity new create 300 7986875 mom yes rot yes dist gaussian units box velocity boundary set 0.0 0.0 0.0 units box compute Thot hot temp compute Tcold cold temp compute Tmid mid temp #1st equilibration run fix 1 boundary setforce 0.0 0.0 0.0 fix 2 new nvt temp 300 300 0.025 thermo 20000 run 2000 velocity new scale 300 unfix 2 #2nd equilibration run fix 2 new nve fix hot all heat 1 50 region hot fix cold all heat 1 -50 region cold thermo_style custom step temp c_Thot c_Tcold c_Tmid thermo 20000 run 20000 #thermal conductivity calculation compute ke all ke/atom variable temp atom c_ke*1.6e+4/(1.5*1.38) fix 3 new ave/spatial 1 10000 10000 x (0,0,0) 2 v_temp file temp.profile thermo_style custom step temp c_Thot c_Tcold c_Tmid neighbor 1.0 bin neigh_modify every 1 delay 5 check yes timestep 0.001 dump 1 all atom 10 bn.lammpstrj run 500000 真的非常感谢！ |
回复此楼» 收录本帖的淘帖专辑推荐
 分子动力学-初级问题资源 |
» 猜你喜欢
明天E口面上会评
已经有5人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有3人回复
-
大龄残疾硕士的一点执念
已经有27人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
不要再数国自然申请书的 filecode 的分隔符个数了
已经有15人回复
-
祈祷青基必中
已经有16人回复
2楼2016-10-18 20:37:29
1195222016
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 383.7
- 散金: 15
- 红花: 1
- 帖子: 68
- 在线: 45.1小时
- 虫号: 3539198
- 注册: 2014-11-16
- 性别: GG
- 专业: 无机材料化学
3楼2016-10-27 09:20:38
4楼2018-06-10 11:31:48
wanziwen
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 407.9
- 散金: 10
- 帖子: 233
- 在线: 154.2小时
- 虫号: 3836618
- 注册: 2015-04-27
- 性别: GG
- 专业: 机械摩擦学与表面技术
5楼2018-08-03 16:38:57
6楼2018-11-22 09:42:34
50