Abstract:Cu has been widely used as a substrate material for graphene growth. To understand the atomistic mechanism of growth, an efficient and accurate method for describing Cu-C interactions is necessary, which is the prerequisite of any possible large-scale molecular simulation studies. The semi-empirical density-functional tight-binding (DFTB) method has a solid basis from the density functional theory (DFT) and is believed to be a good tool for achieving a balance between efficiency and accuracy. However, existing DFTB parameters cannot provide a reasonable description of the Cu surface structure. At the same time, DFTB parameters for Cu-C interactions are not available. Therefore, it is highly desirable to develop a set of DFTB parameters that can describe the Cu-C system, especially for surface reactions. In this study, a parametrization for Cu-C systems within the self-consistent-charge DFTB (SCC-DFTB) framework is performed. One-center parameters, including on-site energy, Hubbard, and spin parameters, are obtained from DFT calculations on free atoms. Two-center parameters can be calculated based on atomic wavefunctions. The remaining repulsive potential is obtained as the best compromise to describe different kinds of systems. Test calculations on Cu surfaces and Cu- or C atom-adsorbed Cu surfaces indicate that the obtained parameters can generate reasonable geometric structures and energetics. Based on this parameter set, carbon dimerization on the Cu(111) surface has been investigated via molecular dynamics simulations. Since they are the feeding species for graphene growth, it is important to understand how carbon dimers are formed on the Cu surface. It is difficult to observe carbon dimerization in brute-force MD simulations even at high temperatures, because of the surface structure distortion. To study the dimerization mechanism, metadynamics simulations are performed. Our simulations suggest that carbon atoms will rotate around the bridging Cu atom after a bridging metal structure is formed, which eventually leads to the dimer formation. The free energy barrier for dimerization at 1300 K is about 0.9 eV. The results presented here provide useful insights for understanding graphene growth.

Thermodynamically Driven Formation of Intercalated Cu Carpets from Supported Cu Pyramids on MoS2.

Unusual flat and extended morphology of intercalated Cu under MoS2

Edge Reconstruction-Dependent Growth Kinetics of MoS2

Adsorption, Intercalation, Diffusion, and Adhesion of Cu at the 2H-Mos_{2} (0001) Surface from First-Principles Calculations

DFT calculations of the structure and stability of copper clusters on MoS2

On-Surface Self-Assembly Kinetic Study of Cu-Hexaazatriphenylene 2D Conjugated Metal-Organic Frameworks on Coinage Metals and MoS2 Substrates

Ab initio characterization of hybrid MOF-MXenes surfaces: The case of Cu-pyridyl on Ti2CO2

Constant-Potential Machine Learning Molecular Dynamics Simulations Reveal Potential-Regulated Cu Cluster Formation on MoS$_{2}$

Constant-Potential Machine Learning Molecular Dynamics Simulations Reveal Potential-Regulated Cu Cluster Formation on MoS_2

The dynamics of copper intercalated molybdenum ditelluride

Quantum-confinement-induced Periodic Surface States in Two-Dimensional Metal-Organic Frameworks

Formation of Multilayer Cu Islands Embedded Beneath the Surface of Graphite: Characterization and Fundamental Insights

A Molecular Dynamics Study of Carbon Dimerization on Cu(111) Surface with Optimized DFTB Parameters

First-Principles Thermodynamics of Graphene Growth on Cu Surfaces

Spontaneous self-intercalation of copper atoms into transition metal dichalcogenides

Studies of two-dimensional MoS2 on enhancing the electrical performance of ultrathin copper films

Strong thermal transport along polycrystalline transition metal dichalcogenides revealed by multiscale modelling for MoS2

Theoretical realization of two-dimensional M3(C6X6)2(M= Co, Cr, Cu, Fe, Mn, Ni, Pd, Rh and X= O, S, Se) metal-organic frameworks

Concentric Advancing Front Corrugations and Multiple Ordered Growth of 2D Mo2C Crystals

Multi-terminal electrical transport measurements of molybdenum disulphide using van der Waals heterostructure device platform