应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » LAMMPS运行结果出现nan
北京石油化工学院2026年研究生招生接收调剂公告
211/3返回列表123下一页
查看: 4557  |  回复: 20
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

WJ-MS

新虫 (小有名气)

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

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

Screenshot.png
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
加油！
1楼 2012-12-07 11:12:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

Bessel

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
jiaoyixiong: 金币+3, 鼓励交流 2012-12-09 23:11:22
WJ-MS: 金币+7, ★★★很有帮助 2012-12-17 09:23:36
引用回帖:
5楼: Originally posted by WJ-MS at 2012-12-09 14:28:05
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 bo ...

create_atoms        1 random 106 16 regbox units box
create_atoms        2 random 106 16 regbox units box

可能问题出在这了，两种粒子都是random分布的，有的可能重合在一起了。所以能量变得无穷大。

在第一次run之前minimize一下。
赞一下(1人) 回复此楼
9楼2012-12-09 22:19:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

★ ★ ★ ★
感谢参与，应助指数 +1
zh1987hs: 金币+1, 谢谢 2012-12-11 21:30:01
WJ-MS: 金币+3, ★★★很有帮助 2012-12-17 09:26:38
我记得这样是初始模型建立的好不引起的。还有把你的 in文件贴出来，让大家看看

你有木有先进行能量最小化？
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

WJ-MS
2楼2012-12-07 14:04:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Bessel

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
jiaoyixiong: 金币+2, 鼓励交流 2012-12-08 20:08:34
初始构型的问题，在run之前先minimize一下。

http://lammps.sandia.gov/doc/minimize.html
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

WJ-MS
3楼2012-12-08 04:58:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tracel

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
jiaoyixiong: 金币+2, 鼓励交流 2012-12-08 20:08:51
如前几楼说的，构型不好，导致能量太大，超出了lammps所能记录的数字大小。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

WJ-MS
4楼2012-12-08 19:47:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WJ-MS

新虫 (小有名气)
引用回帖:
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
赞一下 回复此楼
加油！
5楼2012-12-09 21:28:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WJ-MS

新虫 (小有名气)
送鲜花一朵
引用回帖:
3楼: Originally posted by Bessel at 2012-12-08 04:58:17
初始构型的问题，在run之前先minimize一下。

http://lammps.sandia.gov/doc/minimize.html

好的 我试一试哈
回复此楼
加油！
6楼2012-12-09 21:28:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WJ-MS

新虫 (小有名气)
送鲜花一朵
引用回帖:
4楼: Originally posted by tracel at 2012-12-08 19:47:03
如前几楼说的，构型不好，导致能量太大，超出了lammps所能记录的数字大小。

我回复了2楼in文件了，麻烦帮我看看呢，哪儿出了问题了……新手上路，还请多多直接哈
赞一下 回复此楼
加油！
7楼2012-12-09 21:30:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WJ-MS

新虫 (小有名气)
送鲜花一朵
引用回帖:
2楼: Originally posted by jiaoyixiong at 2012-12-07 14:04:29
我记得这样是初始模型建立的好不引起的。还有把你的 in文件贴出来，让大家看看

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

新手上路，还请多多指教呀！！！
赞一下 回复此楼
加油！
8楼2012-12-09 21:31:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

★ ★ ★
zh1987hs: 金币+3, 谢谢 2012-12-12 11:02:20
引用回帖:
9楼: Originally posted by Bessel at 2012-12-09 22:19:11
create_atoms        1 random 106 16 regbox units box
create_atoms        2 random 106 16 regbox units box

可能问题出在这了，两种粒子都是random分布的，有的可能重合在一起了。所以能量变得无穷大。
...

楼上所言甚是！

如果我没有理解错的话，楼主是想模拟氯化钠晶体的融化吧？

如果是这样的话，你的模型貌似有问题，
create_atoms        1 random 106 16 regbox units box
这个命令建立的原子位置是随机的。而氯化钠晶体中氯和钠的位置可不是这样的。你得根据晶体类型和晶格常数来生成初始位置吧？？？
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

WJ-MS
10楼2012-12-09 23:14:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 WJ-MS 的主题更新
211/3返回列表123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 318求调剂 +12 ykyhsa 2026-04-05 14/700 2026-04-06 17:46 by fuyu_
[考研] 285求调剂 +5 mapmath 2026-04-06 6/300 2026-04-06 17:18 by 蓝云思雨
[考研] 287分求调剂 有专利国奖一志愿哈工大085406 +3 白易辰 2026-04-06 4/200 2026-04-06 15:24 by going home
[考研] 一志愿武汉理工大学080200机械工程308分，求调剂 +4 终不似从前 2026-04-05 4/200 2026-04-06 11:46 by 考研学校招点人
[考研] 材料与化工371求调剂 +14 陪琳看海 2026-04-04 15/750 2026-04-06 06:59 by houyaoxu
[考研] 材料调剂 +14 壹贰贰亿 2026-04-04 14/700 2026-04-05 23:31 by 来看流星雨10
[考研] 301求调剂 +3 XYPLR 2026-04-05 4/200 2026-04-05 19:07 by XYPLR
[考研] 复试调剂 +8 春日来信- 2026-04-03 8/400 2026-04-05 18:58 by 蓝云思雨
[考研] 一志愿同济大学323分（080500）求调剂 +8 yikeniu 2026-04-01 8/400 2026-04-05 18:15 by cql1109
[考研] 301求调剂 +18 骆驼男人 2026-04-02 18/900 2026-04-04 20:33 by 蓝云思雨
[考研] 22408，264求调剂 +3 ywh729 2026-04-03 4/200 2026-04-04 11:04 by ywh729
[考研] 0856调剂 +8 曲听筠 2026-03-30 8/400 2026-04-04 08:46 by tianyyysss
[考研] 387求调剂 +4 爱吃片豆土 2026-04-03 5/250 2026-04-04 08:10 by 岸上的一条鱼
[考研] 材料科学与工程考研 +10 拯救皮特托先生 2026-04-02 10/500 2026-04-03 23:57 by userper
[考研] 281求调剂 +10 aaawhy 2026-04-03 10/500 2026-04-03 21:42 by lbsjt
[考研] 338求调剂 +4 zzz，，r 2026-04-03 4/200 2026-04-03 16:39 by lijunpoly
[考研] 一志愿复旦材料，英一专硕，总分357调剂 +4 1050389037 2026-04-02 5/250 2026-04-02 21:40 by dongzh2009
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +9 双马尾痞老板2 2026-04-01 9/450 2026-04-02 12:14 by oooqiao
[考研] 379求调剂 +3 ?苦瓜不苦 2026-04-01 3/150 2026-04-01 20:09 by 暮云清寒
[考研] 材料调剂 +10 Eujd1 2026-03-31 11/550 2026-04-01 11:23 by ivanqyq
信息提示
关闭
请填处理意见
关闭