Dislocations in stacking and commensurate-incommensurate phase transition in bilayer graphene and hexagonal boron nitride

Dislocations corresponding to a change of stacking in two-dimensional hexagonal bilayers, graphene and boron nitride, and associated with boundaries between commensurate domains are investigated using the two-chain Frenkel-Kontorova model on top of ab initio calculations. Structural transformations of bilayers in which the bottom layer is stretched and the upper one is left to relax freely are considered for gradually increased elongation of the bottom layer. Formation energies of dislocations, dislocation width, and orientation of the boundary between commensurate domains are analyzed depending on the magnitude and direction of elongation. The second-order phase transition from the commensurate phase to the incommensurate one with multiple dislocations is predicted to take place at some critical elongation. The order parameter for this transition corresponds to the density of dislocations, which grows continuously upon increasing the elongation of the bottom layer above the critical value. In graphene and metastable boron nitride with the layers aligned in the same direction, where elementary dislocations are partial, this transition, however, is preceded by formation of the first dislocation at the elongation smaller than the critical one. The phase diagrams including this intermediate state are plotted in coordinates of the magnitude and direction of elongation of the bottom layer.The authors also acknowledge the Russian Foundation of Basic Research (Grant No. 14-02-00739-a) and computational time on the Supercomputing Center of Lomonosov Moscow State University and the Multipurpose Computing Complex NRC “Kurchatov Institute.” I.L. acknowledges the financial support from Grupos Consolidados UPV/EHU del Gobierno Vasco (Grant No. IT578-13) and H2020-NMP-2014 project “MOSTOPHOS” (Project No. 646259).Peer reviewe