24小时热门版块排行榜    

查看: 1276  |  回复: 0

qcl龍

铜虫 (初入文坛)

[求助] lammps fix rigid 命令模拟CO分子运动

我用下面的输入文件运行第二步就出现原子丢失,请教下是什么原因,fix rigid 针对气体分子的体系要怎么用,ERROR: Lost atoms: original 200 current 192 (../thermo.cpp:430)

units     metal
boundary   p p p
atom_style  full
              
read_data   CO.data

kspace_style ewald 0.00001
special_bonds lj/coul 0.0 0.0 1.0           
pair_style   lj/cut/coul/long 8
bond_style   harmonic

bond_coeff   1 1400  1.13      #C-O
pair_coeff    1 1   0.0011   3.83 8 #C-C
pair_coeff    2 2   0.0069   3.12 8 #O-O

neighbor   2.0 bin
neigh_modify every 1 delay 0 check yes

velocity      all create 300 56789 rot yes mom yes dist gaussian units box
fix      1 all rigid/nvt/small molecule temp 300 300 100


thermo    1
thermo_style  custom step temp pe etotal press vol
timestep   1
dump     1 all custom 1 dump.lammpstrj id type x y z
dump_modify  1 sort id
run      2000000
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qcl龍 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 ttH3f1oA6ace 2026-11-04 8/400 2026-11-05 13:44 by My20PsvIExYO
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-05 13:33 by My20PsvIExYO
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 kvewwGtKktqB 2026-11-04 6/300 2026-11-05 13:22 by My20PsvIExYO
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 9/450 2026-11-05 13:22 by My20PsvIExYO
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-05 13:13 by My20PsvIExYO
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 ncTVBM6AOJHF 2026-11-04 9/450 2026-11-05 13:11 by My20PsvIExYO
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ncTVBM6AOJHF 2026-11-04 9/450 2026-11-05 13:11 by My20PsvIExYO
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 13:11 by My20PsvIExYO
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 ncTVBM6AOJHF 2026-11-04 7/350 2026-11-05 13:00 by My20PsvIExYO
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ncTVBM6AOJHF 2026-11-04 8/400 2026-11-05 13:00 by My20PsvIExYO
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 6/300 2026-11-05 12:16 by My20PsvIExYO
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 ncTVBM6AOJHF 2026-11-04 7/350 2026-11-05 11:12 by PYZmcHr5kSaq
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 11:10 by PYZmcHr5kSaq
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 ncTVBM6AOJHF 2026-11-04 8/400 2026-11-05 11:01 by PYZmcHr5kSaq
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ncTVBM6AOJHF 2026-11-04 8/400 2026-11-05 10:59 by PYZmcHr5kSaq
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-05 10:48 by PYZmcHr5kSaq
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 10:48 by PYZmcHr5kSaq
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-05 10:39 by PYZmcHr5kSaq
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 hNpp4Y2CGvSC 2026-11-04 9/450 2026-11-05 09:20 by PYZmcHr5kSaq
[基金申请] 各位基础老师们,看完中央《从严治党》的决定后,对学术界更加有信心了吧 +5 iaeyuan 2026-11-01 5/250 2026-11-04 14:43 by jinchengzi
信息提示
请填处理意见