Mean square displacement of intruders in freely cooling multicomponent granular mixtures

The mean square displacement (MSD) of intruders (tracer particles) immersed in a multicomponent granular mixture made up of smooth inelastic hard spheres in a homogeneous cooling state is explicitly computed. The multicomponent granular mixture is constituted by $s$ species with different masses, diameters, and coefficients of restitution. In the hydrodynamic regime, the time decay of the granular temperature of the mixture gives rise to a time decay of the intruder's diffusion coefficient $D_0$. The corresponding MSD of the intruder is determined by integrating the corresponding diffusion equation. As expected from previous works on binary mixtures, we find a logarithmic time dependence of the MSD which involves the coefficient $D_0$. To analyze the dependence of the MSD on the parameter space of the system, the diffusion coefficient is explicitly determined by considering the so-called second Sonine approximation (two terms in the Sonine polynomial expansion of the intruder's distribution function). The theoretical results for $D_0$ are compared with those obtained by numerically solving the Boltzmann equation by means of the direct simulation Monte Carlo method. We show that the second Sonine approximation improves the predictions of the first Sonine approximation, especially when the intruders are much lighter than the particles of the granular mixture. In the long-time limit, our results for the MSD agree with those recently obtained by Bodrova [Phys. Rev. E \textbf{109}, 024903 (2024)] when $D_0$ is determined by considering the first Sonine approximation.