Graphene layer of hybrid graphene/hexagonal boron nitride model upon heating

Hybrid graphene/h-BN model is studied via molecular dynamics simulation to observe the evolution of graphene layer upon heating. Model containing 20,064 atoms is heated up from 50 to 8000 K via Tersoff and Lennard-Jones potentials. Various thermodynamic quantities, structural characteristics, and the occurrence of liquid-like atoms are studied. The Lindemann criterion for 2D case is calculated and used to observe the appearance of liquid-like atoms. The atomic mechanism of structural evolution upon heating is analyzed on the basis of the occurrence/growth of liquid-like atoms, the formation of clusters, the coordination number, and the ring statistics. The liquid-like atoms tend to form clusters and the largest cluster increases slightly in order to form a single largest cluster of liquid-like atoms. The other models such as free-standing graphene, zigzag GNR, and armchair GNR are also presented to have an entire picture about the evolution of graphene upon heating in different models. Note that the largest clusters of free-standing graphene as well as zigzag GNR, and armchair GNR tend to decrease to form a ring-like 2D liquid carbon.