Shape Comparison Between Non-equilibrium Droplets and Equilibrium Crystal Shapes for the 2D Ising Model

We use the Boost Graph Library to identify droplets in the single-droplet region of metastable decay for the 2D Ising model (using a single-spin flip algorithm), and likewise for equilibrium crystals (using Kawasaki dynamics). As a quantitative way of comparing shapes, we measure the number of broken bonds for the droplets/crystals and also (using the \texttt{igraph} library) the graph-theoretic diameter. Both measurements show that the droplets and crystals have the same shapes on average. This is consistent with earlier research, which used the nucleation rate of droplets to compare their surface free energy with the free energy of equilibrium crystals. The advantage of the present method is that it can be extended to nucleation and growth on irregular lattices.