In my input script I have used this command:
thermo_style custom step temp vol pe evdwl ebond eangle edihed eimp etotal
by this command I am going to study each energy separately. When I run the simulation, pe and evdwl are equal, ebond, eangle, edihed, eimp are equal to zero which is strange. For pair_style and pair_coeff I have these lines in the input script:
pair_style airebo 2.5 1.0 1.0
pair_coeff * * …/…/potentials/CH.airebo C C
What do you think about this problem? Why pe and evdwl are equal and other energies are zero? Is this because of the pair_style I am using?