problem to simulate silicon nanotube

Dear all

I wanna simulate a silicon nanotube .my script is below and its good for carbon nanotubes but for silicon nanotube , the structure of nanotubes take apart after running.

would you please help me where is the problem ?

boundary s s s

atom_style atomic
units metal

read_data data_Sint.txt

mass 1 28.085

pair_style tersoff
pair_coeff * * SiO.tersoff Si

velocity all create 300 1234 dist gaussian mom yes rot yes

fix 1 all nve

timestep 0.001

thermo 100

thermo_style custom step etotal ke

dump 1 all xyz 10 dump_Sint22.xyz

run 5000

Best Regards
Hassan Shirzadi

Dear all

I wanna simulate a silicon nanotube .my script is below and its good for
carbon nanotubes but for silicon nanotube , the structure of nanotubes take
apart after running.

would you please help me where is the problem ?

could be lots of things: unsuitable geometry with high initial
potential energy, bad choice of parameters.

in any case, these are not LAMMPS problems, but problems in doing MD
and the science of your specific system. thus you should rather talk
to your adviser or the person teaching you how to do MD (if they are
not one and the same).

axel.