Computational approach to (ZnS)_{i} nanoclusters in ionic liquids

Unique and attractive properties have been predicted for II-VI-type semiconductor nanoclusters within the field of nanotechnology. However, the low reaction kinetics within the usual solvents gives only thermodynamic control during their production process, making the obtention of different metastable polymorphs extremely difficult. The use of ionic liquids as solvents has been proposed to overcome this problem. Identifying how these nanoclusters are solvated within ionic liquids is fundamental if this strategy is to be pursued. While computational chemistry tools are best suited for this task, the complexity and size of the system requires a careful design of the simulation protocol, which is put forward in this work. Taking as reference the (ZnS)_{12} nanocluster and the [EMIM][EtSO_{4}] ionic liquid, we characterize the interactions between the nanoparticle and first solvation shell by density functional theory calculations, considering most of the solvent implicitly. The DFT results are consistent through different theory levels showing a strong interaction between the Zn atoms of the nanocluster and the [EtSO_{4}^{-}] anion of the ionic liquid. A more realistic representation of the system is obtained by classical MD calculations, for which various classical force fields were considered and several atomic interactions parameterized. This new set of parameters correctly describes the interaction of different (ZnS) nanoclusters, supporting its transferability. The resulting MD simulation shows the formation of a structured ionic liquid solvation shell around the nanocluster with no exchange of ions for at least 5 ns, in agreement with the strong interactions observed in the density functional theory calculations.