This is a 2D-native analog of the 3D shell-growth model: instead of a Fibonacci-sphere projection of the 3D cluster, atoms here are placed on an actual 2D triangular (hexagonal-packed) lattice — the same close-packing every {111} plane of an FCC metal crystal shows in cross-section. Ring k (k = 1, 2, 3…) around the central atom contains exactly
atoms in ring k = 6k
the classic "centered hexagonal number" sequence (6, 12, 18, 24… giving totals 7, 19, 37, 61…). Summing rings 1..K plus the center atom gives the total atom count:
N(K) = 1 + Σ(6k) for k=1..K
= 1 + 3K(K+1)
Only the outermost ring sits at the disk's edge, so the surface fraction is 6K / N(K). Because N grows as K² while the edge ring grows only as K, shrinking the disk makes edge atoms dominate — at K=1 (7 atoms) 86% of atoms are on the rim; by K=8 (217 atoms) that falls to about 22%, and a macroscopic 2D crystal (K → ∞) tends to 0%. This is the exact same volume-vs-surface scaling argument as the 3D model, just one dimension down (perimeter ~ K vs area ~ K² here, instead of surface ~ K² vs volume ~ K³ in 3D) — the underlying physics (why small clusters are surface-dominated) is identical, computed independently from the hex-ring geometry rather than read off the 3D scene.
Interior atoms on a 2D triangular lattice have 6 nearest neighbours; rim atoms are missing roughly 2 of those, leaving about 4. That reduced coordination is why edge/surface atoms are more reactive and more mobile — the same physical reason real 3D nanoparticles are better catalysts and melt at a lower temperature. The simulator estimates this with:
avg. coordination ≈ (6·N_interior + 4·N_surface) / N
Positions are exact axial hex-grid coordinates (pointy-top rings, generated by walking 6 lattice directions), not an approximation — the atom counts per ring and resulting surface fraction / coordination numbers follow directly from the combinatorial formulas above. The physical spacing d uses the same per-metal nearest-neighbour distances (≈ FCC lattice constant / √2) as the 3D model: disk extent ≈ (2K + 1)·d.