The plot looks the same for 4 chains. I ran 30 to see, if there is a difference. The rhat is quite large (>3), but I assumed that was due to the different chains converging onto the different modes.
I also ran the same model with pm.Metropolis() and 10 chains and the result is quite similar still.
Any suggestions what could be happening?