Zolotarev proved a duality result that relates stable densities with different indices. In this paper, we show how Zolotarev duality leads to some interesting results on fractional diffusion. Fractional diffusion equations employ fractional derivatives in place of the usual integer order derivatives. They govern scaling limits of random walk models, with power law jumps leading to fractional derivatives in space, and power law waiting times between the jumps leading to fractional derivatives in time. The limit process is a stable Levy motion that models the jumps, subordinated to an inverse stable process that models the waiting times. Using duality, we relate the density of a spectrally negative stable process with index $1<alpha<2$ to the density of the hitting time of a stable subordinator with index $1/alpha$, and thereby unify some recent results in the literature. These results also provide a concrete interpretation of Zolotarev duality in terms of the fractional diffusion model.