I want to simulate electrons Pauli repulsion through a potential that depends not trivially on their relative momentum (Dorso et al, 1984 PRB) through a Gaussian interaction: V(p_{ij}) ~ exp(-(p_i - p_j)^2).
Is there around LAMMPS any velocity dependent potential? I realize it could be quite hard (different, at least) to make it evolve through any symplectic integrator, verlet-like (I must confess I cannot fully understand how to develop the actual algorithm, but I think I can do it with time and patience).
So, long story short: I need to simulate interaction between particles depending on their relative momentum. Can I do that on LAMMPS?
Thanks,
Pablo