first peak shift to larger distances in comparison to the pmf minima. With the increase of the alkyl chain of the IL cation (two bottom rows of Figure ), the radial distribution function curves start to display several peaks that match the positions of the corresponding pmf minima. For the hydrophilic NPs dispersion in C8, no aggregation was observed up to 1000 ns of simulation, a time at which every NP was bonded to at least one other in the corresponding case in C4, and even at 2500 only 2 dimers were formed from a system with 64 NPs. This results in only very small variations of the NP-NP interaction energy when compared to the dispersion of the same NPs in C4 (Figure ). Although this simulation isn't long enough to reach an equilibrium state, at which certainly every particle would be aggregated since the Δ ag G is very negative, the g(r) patterns change very slowly after 1000 ns, indicating that no huge change would be expected until simulation times orders of magnitude larger than performed here. Notice that simulations with coarse grained force fields as the ones performed here tend to result in faster dynamics than in real world due to the smooth potential surface in the coarse grained models, which is also noticed by the underestimated zero shear viscosity of the C4 obtained with the same force field. Hence, in the real world, the kinetic stabilization of similar NPs dispersions is likely to be even greater.