A microscopic model of the formation of an alloy on the interface has been constructed, which takes into account the exchange of atoms with the substrate atoms and the "floating up" of the latter into the upper layers in the process of epitaxial growth. The self-consistent calculations of atomic magnetic moments of spatially inhomogeneous structures obtained in this case are used for the interpretation of data of Mössbauer spectroscopy. The proposed scenario of mixing leads to the appearance of a preferred direction in the sample and the asymmetry of interfaces in the direction of epitaxial growth. In the multilayer M ₁/M ₂ (M _1,2 = Fe, Cr, V, Sn, or Ag) systems, this asymmetry makes it possible to understand the difference in the magnetic behavior of M ₁-on M ₂ and M ₂-on-M ₁ interfaces which has been observed experimentally. The correlation between the calculated distributions of magnetic moments and the measured distributions of hyperfine fields at iron atoms confirms the assumption about their proportionality for a broad class of metallic multilayer systems. However, a linear decrease of hyperfine fields at the ⁵⁷Fe nuclei with increasing number of impurity atoms among the nearest and next-nearest neighbors is not confirmed for Fe/Cr systems, although is correct in Fe/V superlattices. In the Fe/Cr multilayer systems, the experimentally measured value of magnetoresistance grows with increasing fraction of the "floated up" atoms of ⁵⁷Fe. Thus, it is the bulk scattering by impurity atoms that gives the basic contribution to the effect of giant magnetoresistance. The problem of the influence of mixing and adsorption of hydrogen in the vanadium layers on the state of the spin-density wave in V/Cr superlattices has been considered.