text: The kinetic Monte Carlo model for two-dimensional dewetting nanoﬂuids [33] was ﬁrst proposed 
  
  in Ref. [35] and extended to include next-nearest neighbour interactions in [37]. The two key 
  
  assumptions used are: (i) the relevant processes can be mapped on to a two-dimensional lattice 
  
  gas model, thereby neglecting continuous changes in the thickness of the evaporating ﬁlm, and (ii) 
  
  all relevant dynamics results from diffusing nanoparticles and evaporating/condensing solvent. 
  
  The model builds on an Ising-type model for the liquid-gas phase transition. The surface is divided 
  
  up into a regular array of lattice sites whose size is dictated by the nanoparticles. One then con- 
  
  siders each lattice site to be occupied either by a nanoparticle, liquid or vapour. This effectively 
  
  maps the system onto a two-dimensional two-component lattice gas having two ﬁelds n and l. The 
  
  resulting three possible states of a cell are: liquid (l = 1, n = 0), nanoparticle (l = 0, n = 1), 
  
  and vapour (l = 0, n = 0, i.e., cell empty). The energy of an overall conﬁguration is given by the 