Poor convergence using interpolated distribution for observed data

Questions
3

Using Metropolis sampling give the same bad results…

The page where I get the function does say this (in the middle of the page)

It is still possible to continue using NUTS sampling method because Interpolated class implements calculation of gradients that are necessary for Hamiltonian Monte Carlo samplers.

So the NUTS sampler should work ok… But it doesn’t ;(