24小时热门版块排行榜    

查看: 2655  |  回复: 5

missingwudi

新虫 (初入文坛)

[求助] 用lammps做CH4的MSD均方位移过程中出现问题 已有1人参与

in文件是这样的:
units         real
boundary         p p p
atom_style         full
angle_style         harmonic
bond_style         harmonic

pair_style  lj/cut 10.0
pair_modify mix arithmetic  

read_data        data.final_0.3
replicate    1 1 1

group         A_C type 3
group   CH4 type 1 2

timestep  0.5

neighbor         2.0 bin
neigh_modify   delay 5 every 1 check yes

velocity CH4 create 300.0 4928459
fix   freeze A_C setforce 0.0 0.0 0.0

minimize 1.0e-4 1.0e-6 100 1000

fix    1  CH4 rigid single
fix    2  CH4 nvt temp 300.0 300.0 0.01

compute methane_temp   CH4 temp

thermo_style custom step temp pe etotal press vol c_methane_temp
thermo    1000

dump  1 all custom 1000 dump.lammpstrj id type x y z
dump_modify 1 sort 1 first no

log                 log.lammps
run                 1000000

undump  1
unfix   2

fix   3   CH4 nve
fix   4  CH4 temp/rescale 100 1.0 1.1 0.02 0.5

dump  2 all custom 10000 nve.lammpstrj id type xu yu zu
dump_modify 2 sort 1
dump 3 all xyz 10000 nve.xyz
dump_modify 3 sort 1

compute                msd_CH4 CH4 msd com yes
fix                5 CH4 ave/time 1 1 10000 c_msd_CH4 file msd_CH4 start 1000000 mode vector

run                 1000000
问题:导出来的文件出现了错误
dump.lammpstrj文件:
ITEM: TIMESTEP
1000
ITEM: NUMBER OF ATOMS
3310
ITEM: BOX BOUNDS pp pp pp
0 39.9
0 39.9
0 39.9
ITEM: ATOMS id type x y z
1 1 -nan -nan -nan
2 2 -nan -nan -nan
3 2 -nan -nan -nan
4 2 -nan -nan -nan
5 2 -nan -nan -nan
6 1 -nan -nan -nan
7 2 -nan -nan -nan
8 2 -nan -nan -nan
9 2 -nan -nan -nan
10 2 -nan -nan -nan
11 1 -nan -nan -nan
12 2 -nan -nan -nan
13 2 -nan -nan -nan
14 2 -nan -nan -nan
15 2 -nan -nan -nan
16 1 -nan -nan -nan
17 2 -nan -nan -nan
18 2 -nan -nan -nan
19 2 -nan -nan -nan
20 2 -nan -nan -nan
21 1 -nan -nan -nan
22 2 -nan -nan -nan
23 2 -nan -nan -nan
24 2 -nan -nan -nan
25 2 -nan -nan -nan
26 1 -nan -nan -nan
27 2 -nan -nan -nan
28 2 -nan -nan -nan
29 2 -nan -nan -nan
30 2 -nan -nan -nan
31 1 -nan -nan -nan
32 2 -nan -nan -nan
33 2 -nan -nan -nan
34 2 -nan -nan -nan
35 2 -nan -nan -nan
36 1 -nan -nan -nan
37 2 -nan -nan -nan
38 2 -nan -nan -nan
39 2 -nan -nan -nan
40 2 -nan -nan -nan
41 1 -nan -nan -nan
42 2 -nan -nan -nan
43 2 -nan -nan -nan
44 2 -nan -nan -nan
45 2 -nan -nan -nan
46 1 -nan -nan -nan
47 2 -nan -nan -nan
48 2 -nan -nan -nan
49 2 -nan -nan -nan
50 2 -nan -nan -nan
51 1 -nan -nan -nan
52 2 -nan -nan -nan
53 2 -nan -nan -nan
54 2 -nan -nan -nan
55 2 -nan -nan -nan
56 1 -nan -nan -nan
57 2 -nan -nan -nan
58 2 -nan -nan -nan
59 2 -nan -nan -nan
60 2 -nan -nan -nan
61 1 -nan -nan -nan
62 2 -nan -nan -nan
63 2 -nan -nan -nan
64 2 -nan -nan -nan
65 2 -nan -nan -nan
66 1 -nan -nan -nan
67 2 -nan -nan -nan
68 2 -nan -nan -nan
69 2 -nan -nan -nan
70 2 -nan -nan -nan
71 1 -nan -nan -nan
72 2 -nan -nan -nan
73 2 -nan -nan -nan
74 2 -nan -nan -nan
75 2 -nan -nan -nan
76 1 -nan -nan -nan
77 2 -nan -nan -nan
78 2 -nan -nan -nan
79 2 -nan -nan -nan
80 2 -nan -nan -nan
81 1 -nan -nan -nan
82 2 -nan -nan -nan
83 2 -nan -nan -nan
log文件:
WARNING: One or more atoms are time integrated more than once (../modify.cpp:245)
WARNING: Inconsistent image flags (../domain.cpp:645)
Memory usage per processor = 6.25594 Mbytes
Step Temp PotEng TotEng Press Volume methane_
      32   0.10282812    3612.6324     3612.927    16623.972    63521.199     98.81782
    1000         -nan         -nan         -nan         -nan    63521.199         -nan
    2000         -nan         -nan         -nan         -nan    63521.199         -nan
    3000         -nan         -nan         -nan         -nan    63521.199         -nan
    4000         -nan         -nan         -nan         -nan    63521.199         -nan
    5000         -nan         -nan         -nan         -nan    63521.199         -nan
    6000         -nan         -nan         -nan         -nan    63521.199         -nan
    7000         -nan         -nan         -nan         -nan    63521.199         -nan
    8000         -nan         -nan         -nan         -nan    63521.199         -nan
回复此楼

» 猜你喜欢

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

lsloneil

专家顾问 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
missingwudi: 金币+50, ★★★★★最佳答案, 好像是可以的,程序还在跑,等结果出来了再向您讨教 2015-08-25 15:00:13
minimize 1.0e-4 1.0e-6 100 1000
这个先去掉, minimize和fix nvt不能同时用。
2楼2015-08-25 11:24:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

missingwudi

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by lsloneil at 2015-08-25 11:24:33
minimize 1.0e-4 1.0e-6 100 1000
这个先去掉, minimize和fix nvt不能同时用。

开始计算时好的,算着算着又出现问题了,还是同样的问题
1100 2 -nan -nan -nan
1101 1 -nan -nan -nan
1102 2 -nan -nan -nan
1103 2 -nan -nan -nan
1104 2 -nan -nan -nan
1105 2 -nan -nan -nan
1106 1 -nan -nan -nan
1107 2 -nan -nan -nan
1108 2 -nan -nan -nan
1109 2 -nan -nan -nan
1110 2 -nan -nan -nan
1111 1 -nan -nan -nan
1112 2 -nan -nan -nan
1113 2 -nan -nan -nan
1114 2 -nan -nan -nan
1115 2 -nan -nan -nan
1116 1 -nan -nan -nan
1117 2 -nan -nan -nan
1118 2 -nan -nan -nan
1119 2 -nan -nan -nan
1120 2 -nan -nan -nan
1121 1 -nan -nan -nan
1122 2 -nan -nan -nan
1123 2 -nan -nan -nan
1124 2 -nan -nan -nan
1125 2 -nan -nan -nan
1126 1 -nan -nan -nan
1127 2 -nan -nan -nan
1128 2 -nan -nan -nan
1129 2 -nan -nan -nan
1130 2 -nan -nan -nan
1131 1 -nan -nan -nan
1132 2 -nan -nan -nan
1133 2 -nan -nan -nan
1134 2 -nan -nan -nan
1135 2 -nan -nan -nan
1136 1 -nan -nan -nan
1137 2 -nan -nan -nan
1138 2 -nan -nan -nan
1139 2 -nan -nan -nan
1140 2 -nan -nan -nan
1141 1 -nan -nan -nan
1142 2 -nan -nan -nan
1143 2 -nan -nan -nan
1144 2 -nan -nan -nan
1145 2 -nan -nan -nan
1146 1 -nan -nan -nan
1147 2 -nan -nan -nan
1148 2 -nan -nan -nan
1149 2 -nan -nan -nan
1150 2 -nan -nan -nan
1151 1 -nan -nan -nan
1152 2 -nan -nan -nan
1153 2 -nan -nan -nan
1154 2 -nan -nan -nan
1155 2 -nan -nan -nan
1156 1 -nan -nan -nan
1157 2 -nan -nan -nan
1158 2 -nan -nan -nan
1159 2 -nan -nan -nan
1160 2 -nan -nan -nan
1161 1 -nan -nan -nan
3楼2015-08-25 15:32:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)

引用回帖:
3楼: Originally posted by missingwudi at 2015-08-24 19:32:53
开始计算时好的,算着算着又出现问题了,还是同样的问题
1100 2 -nan -nan -nan
1101 1 -nan -nan -nan
1102 2 -nan -nan -nan
1103 2 -nan -nan -nan
1104 2 -nan -nan -nan
1105 2 -nan -nan -nan
...

我不是很清楚,可以尝试以下

1. 采用更小的timestep
2. 先把体系minimize之后再用nvt平衡
3. fix    1  CH4 rigid single
fix    2  CH4 nvt temp 300.0 300.0 0.01
这个用法按照lammps手册上的说法有点可疑,或许可以把fix rigid去掉后试试看,或者把这两个fix结合起来,直接用fix rigid/nvt.

4. 当你体系的温度变得无限大,可以把你lammpstrj 文件在VMD里显示出来,看是否有很奇怪的构型出现。
4楼2015-08-26 08:17:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

missingwudi

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by lsloneil at 2015-08-26 08:17:56
我不是很清楚,可以尝试以下

1. 采用更小的timestep
2. 先把体系minimize之后再用nvt平衡
3. fix    1  CH4 rigid single
fix    2  CH4 nvt temp 300.0 300.0 0.01
这个用法按照lammps手册上的说法有点可 ...

好的,我把fix rigid去掉之后可以计算了,我有个问题想问一下,利用lammps自身的msd语句compute  msd_CH4 CH4 msd com yes计算出来的结果我该怎么处理才能得到扩散系数
5楼2015-08-28 13:05:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)

【答案】应助回帖

引用回帖:
5楼: Originally posted by missingwudi at 2015-08-27 17:05:42
好的,我把fix rigid去掉之后可以计算了,我有个问题想问一下,利用lammps自身的msd语句compute  msd_CH4 CH4 msd com yes计算出来的结果我该怎么处理才能得到扩散系数...

http://muchong.com/bbs/viewthread.php?tid=6899039&fpage=1
6楼2015-08-29 07:32:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 missingwudi 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿武理材料工程348求调剂 +5  ̄^ ̄゜汗 2026-03-19 7/350 2026-03-22 19:44 by 公瑾逍遥
[考研] 289材料与化工(085600)B区求调剂 +3 这么名字咋样 2026-03-22 4/200 2026-03-22 17:56 by 云民大李老师
[考研] 一志愿中南化学(0703)总分337求调剂 +9 niko- 2026-03-19 10/500 2026-03-22 16:08 by ColorlessPI
[考研] 305分求调剂(食品工程) +4 Sxy112 2026-03-21 6/300 2026-03-22 15:26 by 无懈可击111
[考研] 求调剂 +7 Auroracx 2026-03-22 7/350 2026-03-22 12:38 by 素颜倾城1988
[基金申请] 山东省面上项目限额评审 +4 石瑞0426 2026-03-19 4/200 2026-03-22 08:50 by Wei_ren
[考研] 化学调剂 +5 yzysaa 2026-03-21 5/250 2026-03-21 22:12 by peike
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +11 吃吃吃才有意义 2026-03-19 11/550 2026-03-21 18:23 by 学员8dgXkO
[考研] 332求调剂 +3 凤凰院丁真 2026-03-20 3/150 2026-03-21 10:27 by luoyongfeng
[考研] 机械专硕299求调剂至材料 +3 kkcoco25 2026-03-16 4/200 2026-03-21 03:52 by JourneyLucky
[考研] 070300化学319求调剂 +7 锦鲤0909 2026-03-17 7/350 2026-03-21 03:46 by JourneyLucky
[考研] 085700资源与环境308求调剂 +12 墨墨漠 2026-03-18 13/650 2026-03-21 01:42 by JourneyLucky
[考研] 一志愿华南师大 070300(化学)304分求调剂 +3 0703武芊慧雪304 2026-03-18 3/150 2026-03-21 00:48 by JourneyLucky
[考研] 294求调剂材料与化工专硕 +15 陌の森林 2026-03-18 15/750 2026-03-20 23:28 by JourneyLucky
[考研] 工科材料085601 279求调剂 +7 困于星晨 2026-03-17 9/450 2026-03-20 17:38 by 无懈可击111
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
信息提示
请填处理意见