The formalism of the preceding paper is applied to work out the theory of the Gorsky effect, or anelastic relaxation due to the long-range diffusion of interstitials in a host lattice, for non-interacting (low-concentration) interstitials (e.g., H in Nb). It is shown how linear response theory (LRT) provides a number of advantages that simplify the solution of the problem and permit the handling of complications due to specimen geometry and stress inhomogeneity. The multiple-relaxation time creep function of Alefeldet al is first re-derived. Next, the dynamic responseand the short-time behaviour of the creep function are deduced exactly, and theω^−1/2 fall-off of the internal friction at high frequencies is exhibited. Finally, it is pointed out that the true asymptotic behaviour of the dynamic response must be found by going beyond the diffusion equation model. A two-state random walk analysis is used to predict a cross-over to a trueω⁻¹ asymptotic behaviour, and the physical reasons for this phenomenon are elucidated.