A polynomial hybrid Monte Carlo algorithm

We present a simulation algorithm for dynamical fermions that combines the multiboson technique with the hybrid Monte Carlo algorithm. We find that the algorithm gives a substantial gain over the standard methods in practical simulations. We point out the ability of the algorithm to treat fermion zero modes in a clean and controllable manner.