Hi Ju,
try the attached input/data sets and let me know how it goes.
The problem was the electron radial velocity variable should be ervel, instead of evel.
My fault, I just realized we need to update this in the velocities section of the read_data page in the online manual (electron atom-ID vx vy vz ervel)
http://lammps.sandia.gov/doc/read_data.html
Thanks for pointing this out.
Best,
Andres
data.li_plasma (4.16 KB)
in.li_plasma (745 Bytes)