Improved bilinears in lattice QCD with nondegenerate quarks

We describe the extension of the improvement program for bilinear operators composed of Wilson fermions to nondegenerate dynamical quarks. We consider two, three and four flavors, and both flavor nonsinglet and singlet operators. We find that there are many more improvement coefficients than with degenerate quarks, but that, for three or four flavors, nearly all can be determined by enforcing vector and axial Ward identities. The situation is worse for two flavors, where many more coefficients remain undetermined