The sampler by default adds an initial jittering around the initval. You can try to set init="adapt_diag" in pm.sample. But the bigger problem is that the model is still numerically unstable (even though the current parameterization is unbounded). Can you try to reduce the mean of v? That may be getting very large for the later trials if the sequence is large enough.