Higher order hybrid Monte Carlo at finite temperature