I have a system composed by 2 infinite layers that form an interface.
To every step I need to compute the normal vector for each atom.
Considering the distance vectors between the atom and his first neighbors and mediating the vector products between these vectors we obtain the normal vector.
First of all: is there a routine in lammps that compute the normal vector in this way?
Otherwise is there a routine in lammps that finds the first neighbors? And if yes how to use that in another routine?
thank you very much!