I implemented hybrid mcmc (Hamiltonian) in AD Model Builder's random effects module. So far as I know it has not been used much. It should be possible to access it for testing from Ben Bolker's glmmadmb package for R without too much trouble.