Dear all,

I have a question but it may not be about the LAMMPS exactly.

We know that d/dt( 0.5 m Vi^2 ) is equal to Fi.Vi for each particle where

Fi is the force vector on atom i and Vi is its velocity vector.

When I compute the d/dt( 0.5 m Vi^2) numerically using the finite

difference method with a very small time step, I see it is not equal to

Fi.Vi . Is it for using finite difference method?

I am really confused. Any help would be appreciated.

Regards,

Ali Rajabpour