Velocity/Momentum dependent potential

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

So, long story short: I need to simulate interaction between particles
depending on their relative momentum. Can I do that on LAMMPS?

The DPD pairstyle has force components that depend on the relative

momentum of particles.
Carlos

Awesome, I totally missed that! I’ll take a closer look at it, thanks Carlos!