Dear lammps users,
i am trying to simulate the radiation damage in Nickle.initially i taken a domain of size 25a025a025a0.and chosen a atom with PKA energy of 500eV moving initially in 111 direction.also i am using fix dt/reset command to account the energy decrease due to collision.
the problem is, temp of complete domain is reaching around 5000K.and all the atoms are moving randomly.
and also i tried to put the rescaling of velocity at boundary to remove the heat, but it is not working.
what to do?
please help me.
if u have any sample input file of radiation damage simulation, i request u to please send a copy of that.