“…During past years, considerable effort has been dedicated to solve this intricate and challenging problem via different numerical and methodological routes [3,4,5,6,7]. In another contribution [8] we have demonstrated that a suitable combination of two complementary numerical approaches is able to provide a highly satisfactory answer to this yet open issue. To be more specific, we have combined the following two methods: (i) on one hand, an optimization technique, which employs ideas of evolutionary algorithms [9] that is able to predict efficiently and with high reliability ordered equilibrium structures at vanishing temperature; (ii) on the other hand, suitably developed Monte Carlo simulations [10,11] which allow to evaluate accurately, via thermodynamic integration, the thermodynamic properties of a particular ordered structure formed by patchy particles at finite temperature.…”