Unable to find maximum with NUTS. Too many local maxima maybe?

This actually solved my problem! I think what helped the most was switching from an MvNormal with a fixed value for sigma to and a Normal with a prior for sigma.