24小时热门版块排行榜

>论坛更新日志 (4579)
>文献求助 (470)
>虫友互识 (376)
>导师招生 (311)
>考博 (185)
>招聘信息布告栏 (164)
>硕博家园 (157)
>休闲灌水 (115)
>论文道贺祈福 (114)
>博后之家 (100)
>考研 (88)
>论文投稿 (78)
>基金申请 (67)
>教师之家 (67)
>公派出国 (62)
>找工作 (55)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

WJ-MS

新虫 (小有名气)

应助: 2 (幼儿园)
金币: 1559.2
散金: 100
帖子: 194
在线: 92.5小时
虫号: 1697144
注册: 2012-03-17
性别: MM
专业: 凝聚态物性I:结构、力学和

[求助] LAMMPS运行结果出现nan

改写了一段小的程序，运行计算的时候，屏幕上热力学输出都只显示nan，见附图，这是怎么回事阿？

Screenshot.png

回复此楼

» 猜你喜欢

真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
论文投稿，期刊推荐已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
孩子确诊有中度注意力缺陷已经有14人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

加油！

1楼2012-12-07 11:12:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

WJ-MS

新虫 (小有名气)

应助: 2 (幼儿园)
金币: 1559.2
散金: 100
帖子: 194
在线: 92.5小时
虫号: 1697144
注册: 2012-03-17
性别: MM
专业: 凝聚态物性I:结构、力学和

引用回帖:

2楼: Originally posted by jiaoyixiong at 2012-12-07 14:04:29
我记得这样是初始模型建立的好不引起的。还有把你的 in文件贴出来，让大家看看

你有木有先进行能量最小化？

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