Using selected-area low-energy electron diffraction analysis, we showed strict orientational alignment of monolayer hexagonal boron nitride (h-BN) crystallites with Cu(100) surface lattices of Cu foil substrates during atmospheric pressure chemical vapor deposition. In sharp contrast, the graphene-Cu(100) system is well-known to assume a wide range of rotations despite graphene's crystallographic similarity to h-BN. Our density functional theory calculations uncovered the origin of this surprising difference: The crystallite orientation is determined during nucleation by interactions between the cluster's edges and the substrate. Unlike the weaker B-and N-Cu interactions, strong C-Cu interactions rearrange surface Cu atoms, resulting in the aligned geometry not being a distinct minimum in total energy. The discovery made in this specific case runs counter to the conventional wisdom that strong epilayer-substrate interactions enhance orientational alignment in epitaxy and sheds light on the factors that determine orientational relation in van der Waals epitaxy of 2D materials.two-dimensional materials | van der Waals epitaxy | hexagonal boron nitride | graphene | orientational relation R esearch on van der Waals (vdW) heterostructures formed by stacking up various 2D crystals is an emerging field (1) because the numerous possible interfaces may lead to new physics not necessarily associated with the constituent 2D materials and, therefore, to novel applications. Most pioneering works on vdW heterostructures are based on mechanical placement of exfoliated 2D crystal flakes (2-4), where uncertainty in the orientational relation between layers is inevitable. For some purposes, e.g., to enhance graphene's transport by placing it on hexagonal boron nitride (3) (h-BN, or simply BN hereafter), the orientational relation is not important. In other cases, however, a definitive orientational relation is necessary for certain physical properties to arise. Examples include the tunneling behavior of a graphene-insulator-graphene junction (5) and the band structure of graphene placed on top of BN (2).The only practical way to achieve certainty in orientational relation, or azimuthal order, in vdW heterostructures is epitaxial growth. Actually, vdW epitaxy has been studied for decades as a method to overcome lattice mismatch (6, 7) enabled by the relatively weak and flexible vdW interactions between the epilayer and the substrate. Although vdW epitaxy [or quasiepitaxy (8)] is typically incommensurate, there can be a well-defined orientational relation (8), which has been observed in many examples of vdW epitaxy (6, 7, 9-11). In some cases, however, a definitive epilayer-substrate orientational relation is lacking (12, 13). At the moment when vdW epitaxy is shaping into a new thrust area in 2D crystal research, it is imperative to uncover the mechanisms that determine the orientational relation.Epilayer-substrate interactions lead to registry in commensurate epitaxy, which naturally guarantees orientational alignment between the...