24小时热门版块排行榜    

查看: 1183  |  回复: 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龍 的主题更新
信息提示
请填处理意见