I've been using my usual bag of tricks with str(mod), unclass(mod), coef(mod), etc. but, although I can get the fixed effects coefficients. I'd really like to get what's on the Pr(>|z|) column and that p value of 9.567e-10 for the variance component....