Hi
I cannot offer a full answer but I noticed that the sampling stops after 10,000 before the convergence call back kicks in. I increased the n iter to 100,000 and the convergence halts sampling after around 50k or 60k steps.
Now the alpha is closer to 3.4 with approx 25% of the mass above 4.0
