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.

DFT-based Mesoscopic Interactions of N Adatoms on Cu(100)

DFT Total Energy Study on the Atomic Geometry and Adsorption of Cu(100) C(2 × 2)-N Surface

Elastic interactions of Fe adatoms on Cu(111) at the mesoscale

A First-Principles and Mesoscopic Model Analysis for Pair, Trio, and Quarto Interactions of Au Adatoms on Ag(100)

A Density-Functional Study On The Atomic Geometry And Adsorption Of The Cu(100) C(2 X 2)/N Surface

Self‐organization of square nitrogen islands on copper (100) surfaces: Island pair distributions based on an ab‐initio DFT and mesoscopic interaction model

Density Functional Theory Study on the Atomic Structure and Electronic States of Cu(100) (2×22)R45°-O Surface

Density functional theory study on the atomic structure and electronic states of Cu(l00) (√2 × 2 √2)R45°-O surface

First principles study of adsorbed Cu_n (n=1-4) microclusters on MgO(100): structural and electronic properties

${\it Ab \: initio}$ density-functional studies of 13-atom Cu and Ag clusters

Molecular assembly of tetracene on the (2×1)O reconstructed Cu(110) surface

Ab Initio Calculation of Diffusion Barriers for Cu Adatom Hopping on Cu(1 0 0) Surface and Evolution of Atomic Configurations

Theoretical Study on the Atomic Structure, Electronic States, and STM Image of C(2×2) Cu(001)/O Surface

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

Dynamical Behaviour of Ag and Cu Double-Layer Islands on Fcc (111), Surfaces

Density Functional Theory Calculations of Tetracene on Low Index Surfaces of Copper Crystal

Atomic Superlattice Formation Mechanism Revealed by Scanning Tunneling Microscopy and Kinetic Monte Carlo Simulations

Theoretical Study on the Aggregation of Copper Clusters on a Liquid Surface.

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

Diffusion rates of Cu adatoms on Cu(111) in the presence of an adisland nucleated at FCC or HCP sites

Structures and Magic Numbers of Adatom Clusters on Metal Fcc (001) Surfaces