Thanks alot for your reply!
Removing the as_op decorator gives me errors and passing compile_kwargs=dict(backend="NUMBA") to pm.sample, has no appreciable effects on the time of convergence.
For the >=12 parameter versions of the model, the algorithm does not converge even after running for 2 days. Is there something else I can do?