| 查看: 4870 | 回复: 20 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
WJ-MS新虫 (小有名气)
|
[求助]
LAMMPS运行结果出现nan
|
|
改写了一段小的程序，运行计算的时候，屏幕上热力学输出都只显示nan，见附图，这是怎么回事阿？ Screenshot.png |
回复此楼» 猜你喜欢
面上合作单位盖章
已经有5人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有8人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
申请删除本帖
已经有7人回复
-
基金不中，共勉
已经有10人回复
-
国自然面上复盘~欢迎讨论
已经有11人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
怎么看青基中了没有啊
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于towhee 回答zhang_jaj的问题
已经有44人回复
- 我这样运行LAMMPS软件中自带的例子有错误吗？
已经有6人回复
- 【求助】LAMMPS-运行报错
已经有9人回复
- 【讨论】lammps运行eam/fs出错
已经有38人回复
tracel
木虫 (正式写手)
- 应助: 48 (小学生)
- 金币: 4116.7
- 散金: 1052
- 红花: 1
- 帖子: 314
- 在线: 2048.9小时
- 虫号: 831941
- 注册: 2009-08-21
- 专业: 理论和计算化学
WJ-MS4楼2012-12-08 19:47:03
2楼2012-12-07 14:04:29
Bessel
木虫 (正式写手)
- 模拟EPI: 2
- 应助: 47 (小学生)
- 金币: 1506.5
- 散金: 4755
- 红花: 8
- 帖子: 609
- 在线: 234.1小时
- 虫号: 370174
- 注册: 2007-05-13
- 专业: 凝聚态物性I:结构、力学和
【答案】应助回帖
★ ★
感谢参与，应助指数 +1
jiaoyixiong: 金币+2, 鼓励交流 2012-12-08 20:08:34
感谢参与，应助指数 +1
jiaoyixiong: 金币+2, 鼓励交流 2012-12-08 20:08:34
|
初始构型的问题，在run之前先minimize一下。 http://lammps.sandia.gov/doc/minimize.html |
» 本帖已获得的红花（最新10朵）
3楼2012-12-08 04:58:17
WJ-MS
新虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1559.2
- 散金: 100
- 帖子: 194
- 在线: 92.5小时
- 虫号: 1697144
- 注册: 2012-03-17
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
|
3d NaCl melt units metal dimension 3 boundary p p p atom_style charge region regbox block 0 18.8 0 18.8 0 18.8 units box create_box 2 regbox create_atoms 1 random 106 16 regbox units box create_atoms 2 random 106 16 regbox units box mass 1 23.0 mass 2 35.5 set type 1 charge 1.0 set type 2 charge -1.0 pair_style born/coul/long 10.0 8.0 pair_coeff 1 1 0.2640625 0.317 2.340 1.05 0.5 pair_coeff 1 2 0.21125 0.317 2.755 7 8.6875 pair_coeff 2 2 0.1584375 0.317 3.170 72.5 145.625 kspace_style ewald 1.0e-5 neighbor 0.1 bin neigh_modify every 5 delay 0 check no timestep 0.005 compute 1 all msd velocity all create 1073.0 87287 dist gaussian run 0 thermo 5000 thermo_style custom step etotal temp c_1[4] thermo_modify flush yes dump 1 all custom 100000 dump.* id type x y z dump 2 all cfg 100000 a*.cfg id type xs ys zs dump_modify 1 element Na1 Cl2 fix 1 all npt temp 1073 1073 4 aniso 0 0 4 run 500000 unfix 1 fix 1 all nvt temp 1073 1073 0.1 drag 0.2 run 500000 |
5楼2012-12-09 21:28:05
10