Scalable Hamiltonian Monte Carlo Via Surrogate Methods