There may be a few ways to speed up the model itself, but the biggest speed-up will come from calling Turing.setrdcache(true). With just this change, on my machine, Turing sampled 4 chains in 13 minutes. See Automatic Differentiation for a description of when this is safe to use.