The dynamics of interstitial dopants govern the properties of a wide variety of doped crystalline materials. To describe the hopping dynamics of such interstitial impurities, classical approaches often assume that dopant particles do not interact and travel through a static potential energy landscape. Here we show, using computer simulations, how these assumptions and the resulting predictions from classical Eyring-type theories break down in entropically stabilized body-centered cubic (BCC) crystals due to the thermal excitations of the crystalline matrix. Deviations are particularly severe close to melting where the lattice becomes weak and dopant dynamics exhibit strongly localized and heterogeneous dynamics. We attribute these anomalies to the failure of both assumptions underlying the classical description: (i) The instantaneous potential field experienced by dopants becomes largely disordered due to thermal fluctuations and (ii) elastic interactions cause strong dopant-dopant interactions even at low doping fractions. These results illustrate how describing nonclassical dopant dynamics requires taking the effective disordered potential energy landscape of strongly excited crystals and dopantdopant interactions into account.anomalous dynamics | doping | crystals D oping pure crystalline solids with small amounts of interstitial impurities is a widely used method to enhance material properties such as heat and electric conductivity (1-4) or to tailor mechanical properties (5). Prototypical examples include the introduction of carbon atoms in iron crystals to make steel or the doping of plastic crystals with Li ions to create solid-state batteries (4). To ensure longevity of doped materials, it is essential that the spatial homogeneity and transport dynamics of the dopants within the crystal are well controlled and understood. Although theories and models are abundant (5-11), it remains unclear how large thermal excitations of the matrix lattice affect the dynamics of dopants. This becomes of particular interest during the processing of doped crystals, where they are heated close to or beyond their melting point. For example in body-centered cubic (BCC) iron doped with carbon, significant deviations from the exponential increase of diffusivity with temperature, expected from Arrhenius' law, are observed close to the melting temperature where lattice excitations are strong (12). Whereas doping is typically performed to tailor material properties at the macroscopic scale, these enhanced properties emerge from the dynamics and interactions between dopants at the scale of individual atoms (13). In classical theories for dopant dynamics, impurity particles are described as hopping through a potential energy landscape that is set by a perfect lattice symmetry, with transition rates governed by the energy barriers between adjacent interstitial sites and their occupancy (6,7,14). In reality, thermal fluctuations of atoms away from their equilibrium lattice positions will randomize the instantaneous potential energy l...