Thanks @Datseris. I will prepare a minimal working example tonight and post it here.
The problem is not the map_pairwise! function. That function correctly maps a given function on pairs of points within the given distance, generated by CellList. The problem is that in some situations of my simulations the list does not get updated. So the problem should be in UpdateCellList!. I will send the minimal example tonight.