This post Permanent magnet attraction from Maxwell's equations / Lorentz force presents a similar question to mine. I believe it was closed before it was fully answered.
The problem:
In classical electromagnetism (so, putting aside spin).
We model a magnet bar as a current loop and assume a constant external magnetic field perpendicular to the loop surface. We then apply the Lorenz force equation and find that the field applies torque on the current loop but no linear force. So, no pulling or pushing, only rotation.
It seems that the model above should describe the interaction between two bars of magnet when placed south to north. Still, we find that two magnetic bars attract each other.
I suggest the following answer:
Instead of a constant external magnetic field, lets consider the field of a second current loop (this second current loop would model another magnet bar). Now, if we look at the field lines "at and around" the north pole of the current loop we see that they have an orthogonal component facing "north" but also a radial component facing from the axis of the loop outwards.
If we place a current loop A above another current loop B, close to B's north pole, then at every point along the loop A there is a component of the magnetic field facing outwards and orthogonal to the current through A. Appling Lorentz force equation will give us a force perpendicular to loop A surface downwards (towards the loop B).
Does this explain the attraction between two magnets bars for classical electromagnetism?