Thanks a lot for the detailed perspective!

The potential fields are a curious phenomenon. I feel like - and I haven't verified this in any way yet - the smoothness is dictated by the fact that all values are within a not-so-broad range. Then as you "peel" the hexagon from the outer layer using those 6-rings, building the potential field, each next layer inwards shouldn't changed too much.

I had also explored a case of n->inf where the sums turn into integrals, and consecutiveness turns into the uniform distribution. Then it's not so hard to find a solution. It's not in the article, because it turned out to be a dead end for solving the discrete case, and I didn't want to make the article too heavy - but if I remember correctly, there I also observed the smoothness.