Estimating optimal PAC-Bayes bounds with Hamiltonian Monte Carlo