I wonder anything new for this topic? I am wondering, in np.random.exponential(lambda_param, random_size), lambda_param is actually the scale parameter, which is inverse of the actual lambda in exponential distribution. While in what pm.Exponential() takes is actually lambda itself. Not sure if I understand correctly about your original question, but is this possibly the actual reason?