where i is the number of the simulation performed and V i (s) is the total bias accumulated for the duration of the simulation. The work reported in subsequent sections refers to the average value of W i computed over the set of simulations performed for a single system. The estimate of W is robust with respect to the simulation setup, and it has been used to evaluate the solvent-dependence on the thermodynamics of the defect formation for all crystal facets studied.