Problem with pymc softmax transformation producing unexpected results

Can I maybe ask a follow up question?
Is there any obvious reason why the parameters would not be recovered by the model correctly? There is a pattern that the parameter are systematically underestiamted by approximately the same magnitue.

Is this a magnitude to worry about? Shouldn’t they be recovered with the peak at the very value I used for simulation?