  We present a set of related Hybrid Monte Carlo methods to simulate an
arbitrary number of dynamical overlap fermions. Each fermion is represented by
a chiral pseudo-fermion field. The new algorithm reduces critical slowing down
in the chiral limit and for sectors of nontrivial topology.
