text: (cid:88)  (cid:88) 
  E = − 
  lilj − µ  li 
  
  i 
  
  where (cid:80) 
  
  <ij> denotes a sum over nearest neighbour pairs and εll, εnn and εnl are the liquid-liquid, 
  particle-particle and liquid-particle interaction energies, respectively. Fixing the three interaction 
  
  strength parameters εll, εnn, εnl and the effective chemical potential µ determines the equilibrium 
  
  state of the system. We choose εll as unit of energy – i.e. we set εll = 1. 
  
  The hamiltonian determines the equilibrium state and the energy landscape of the system. How- 
  
  ever, as the system ‘dries in’ during the course of the solvent evaporation, the ﬁnal nanoparticle 
  
  conﬁgurations do not necessarily represent equilibrium structures. This implies that the system 
  
  dynamics is of paramount importance. It is determined by the possible Monte Carlo moves, their 
  
  relative frequencies, and the probabilities for their acceptance. Two types of moves are allowed: (i) 
  
  evaporation/condensation of liquid and (ii) diffusion of nanoparticles within the liquid. A mobility 
  
  M corresponds to the ratio of cycles of particle and solvent moves and reﬂects the physical ratio of 