Geometric control of the moire twist angle in heterobilayer flakes

We demonstrate a finite twist-angle stabilization mechanism in lattice-mismatched 2D heterobilayers, which results from the geometric alignment between the flake edges and its moire pattern. Using atomistic simulations of graphene on hexagonal boron nitride flakes with diameters of up to $\sim 2500$Å, we identify robust metastable angles at $\sim {0.61}^{\circ }$ for armchair and $\sim {1.89}^{\circ }$ for zigzag-edged flakes, tunable via in-plane heterostrain. This locking mechanism, which relies on energy barriers that are an order of magnitude larger than those of nearby metastable twist angles, provides a geometric route to precision twist-angle control of two-dimensional heterostructures and to understand the self-orientation of macroscopic flakes.