Some six-coordinate iron (II) coordination compounds exhibit thermal-, optical-, electrical-, magnetic- and pressure-induced switching between the diamagnetic low-spin (LS, S=0) and the paramagnetic high-spin (HS; S=2) states [1]. This may lead to potential application of these complexes in molecular devices such as temperature and pressure sensors [2]. An Ising-like model has been proposed to explain the occurrence of the thermal hysteresis behaviour [3,4] of this switchable solids. In this contribution, the local mean field approximation is applied to solve the Hamiltonian modelling interactions pertaining to 2D nanoparticles embedded in a magnetically-inactive matrix.