  With a generally covariant equation of Dirac fields outside a black hole, we
develop a scattering theory for massive Dirac fields. The existence of modified
wave operators at infinity is shown by implementing a time-dependent
logarithmic phase shift from the free dynamics to offset a long-range mass
term. The phase shift we obtain is a matrix operator due to the existence of
both positive and negative energy wave components.
