Disentangling the effects of many-body forces on depletion interactions

When considering effective interactions in a many-body or thermodynamic system, it is important to recognize that the inherent character of the effective forces between two bodies already incorporates many-body effects that appear unified in a nontrivial way. Working backward to disentangle these effects into the contributions of pairs, triplets, or larger groups is generally complicated and seldom accomplished. However, this could offer essential insights into the structural architecture and self-assembly of complex systems, such as soft materials. In this contribution, we tackle this problem by employing a simulation-based method, here termed contraction of the bare forces (CBF), which, on the one hand, permits a precise evaluation of the effective interactions between colloids [de los Santos-López et al., J. Chem. Phys. 157, 074903 (2022)0021-960610.1063/5.0099919] and, on the other hand, makes it possible to quantify the many-body contributions to these effective interactions. The CBF approach is applied to calculate the depletion forces in asymmetric binary mixtures of hard spheres; the approach is sufficiently sensitive to quantify the effects of higher-order terms on the depletion forces between large particles. To assess the accuracy of our approach, we first explore and confirm the long-established prediction founded on purely geometric considerations, which states that for the particle size ratio, q, smaller than q