Dear users,
I’m having trouble with the airebo potential. I want to simulate diffusion around a dynamic CNT.
However the CNT is destroyed when inserting LJ potentials of diffusing atoms. Temperature has no effect, here is my simplified input file:
I would have expected only little deformation of the tube, is there anything wrong with the input or are the results ok?
dimension 3
boundary p p p
atom_style bond
units metal
#(10,10 tube with 90 atoms around each end)
read_data run.data
mass 1 12
mass 2 16
group CNT type 1
group CH4 type 2
neighbour 2 bin
pair_style hybrid lj/cut 15 airebo 2 0 0
pair_coeff * * airebo CH.airebo C NULL
pair_coeff 1 2 lj/cut 0.005549 3.605
pair_coeff 2 2 lj/cut 0.01276 3.81
fix 1 all nvt 1400 1400 0.02 drag 0.5
timestep 0.0002
thermo_style multi
thermo 200
dump 1 all atom 200 result.lammpstrj
run 15000
thanks a lot in advance
Thomas