Thank you for that inside.
What do you mean by softmax making the model unidentified? I can reproduce the multimodal results only if I replace one component by tt.zeros((nd, 1)) as you suggested. With the original model I get
pm.traceplot(trace, var_names=['decomp'], compact=False)
