The work performed by the biasing potential for the formation of a surface defect was calculated for each surface/solvent combination in order to extract quantitative thermodynamic information for the process. The findings for each case are reported in Fig. . Overall, the absolute value of the work associated with a removal of a surface molecule is dominated by the crystal surface morphology. In particular, the formation of a surface vacancy is least probable on the {100} apolar crystal facets in all solvent cases, while the probability is found to be highest on its polar counter layer. These trends correlate with the observations made in the discussion of the transition times in the previous section, suggesting the presence of intrinsic face-dependent factors controlling the absolute process kinetics and thermodynamics. At the same time, it can be noted that while the values for the work for the case of the {100} apolar crystal face vary only marginally, those calculated for the rest of the crystal facets exhibit a much more significant variation when changing the solvent. This observation indicates that the polar {100}, {002} and {011} crystal faces are much more prone to specificsurface solvent interactions, subsequently affecting the ease of molecule association/dissociation at the crystal surface.