Molecular dynamics simulation of graphene sinking during chemical vapor deposition growth on semi-molten Cu substrate

Copper foil is the most promising catalyst for the synthesis of large-area, high-quality monolayer graphene. Experimentally, it has been found that the Cu substrate is semi-molten at graphene growth temperatures. In this study, based on a self-developed C-Cu empirical potential and density functional theory (DFT) methods, we performed systematic molecular dynamics simulations to explore the stability of graphene nanostructures, i.e., carbon nanoclusters and graphene nanoribbons, on semi-molten Cu substrates. Many atomic details observed in the classical MD simulations agree well with those seen in DFT-MD simulations, confirming the high accuracy of the C-Cu potential. Depending on the size of the graphene island, two different sunken-modes are observed: (i) graphene island sinks into the first layer of the metal substrate and (ii) many metal atoms surround the graphene island. Further study reveals that the sinking graphene leads to the unidirectional alignment and seamless stitching of the graphene islands, which explains the growth of large single-crystal graphene on Cu foil. This study deepens our physical insights into the CVD growth of graphene on semi-molten Cu substrate with multiple experimental mysteries well explained and provides theoretic references for the controlled synthesis of large-area single-crystalline monolayer graphene.