pm.Bernoulli.dist(logit_p = 1.0, shape=2) works for me as well, @bhomass are you on master or the latest release?
pm.Bernoulli.dist(logit_p = 1.0, shape=2)