Depletion forces in dense mixtures of spheres and rods

We evaluate depletion forces in molecular dynamics simulation of a binary mixture of spheres (depleted particles) and rods (depletant particles) for a wide range of densities for both species. This evaluation was carried out using a recently proposed least squares fitting algorithm. We found that the restriction of the rods' rotational degrees of freedom, when the distance between two spheres is less than the rods length, creates a shallow, and apparently linear, attractive force ramp. For intersphere distances smaller than the rods' diameter, a much stronger attractive force is found, and a large repulsive barrier appears between these aforementioned regimes, roughly at the distance of the rods' thickness. The evaluated forces are validated via a comparison of the pairwise correlation functions obtained from molecular dynamics simulation of a mono-disperse sphere fluid, using the evaluated effective forces, against the original (full system) pairwise correlation functions. Agreement is excellent. We also record the angular pairwise correlation function, using the P2(x) Legendre polynomial, and find that for high densities of both species, a local nematic ordering starts to appear. This nematic order may be a factor in the small differences found between original and effective pairwise correlation functions at high densities of rods.