I later found out that this inconsistency in the single point energy does not affect the result of MD simulation.