We present a simulation algorithm for dynamical fermions that combines
the multiboson technique with the hybrid Monte Carlo algorithm. We fi
nd that the algorithm gives a substantial gain over the standard metho
ds in practical simulations. We point out the ability of the algorithm
to treat fermion zero modes in a clean and controllable manner. (C) 1
997 Elsevier Science B.V.