“…12 Molecular dynamics simulations, Monte Carlo simulations and quantum chemical calculations had been employed for some years to elucidate the organization and interactions of ILs close to NPs, 17,24,39,40,41,42,43 nanotubes, 44,45 and solid substrates. 29,44,46,47,48,49 However, to the best of our knowledge, no computer simulation work provided a direct calculation of how the domain segregation and the interaction between IL and NP modulates the aggregation thermodynamics and the kinetic stability of NPs dispersions since most of them consist of simulations with only one NP or only one solid substrate interacting with the IL. Some exceptions are the works of Frost et al 41 In the present work, the potential of mean force was computed between pairs of NPs in two ILs that differ by the length of the cation alkyl chain: 1-butyl-3-methyl-imidazolium tetrafluoroborate ([BMIM][BF4] or C4) and 1-octyl-3-methyl-imidazolium tetrafluoroborate ([OMIM][BF4] or C8) and concentrated dispersions with several NPs were simulated with microseconds time scale simulations.…”