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.

Investigation of Low-Energy Carbon Cluster Depositions on Surfaces by a Molecular Dynamics Simulation

Theoretical studies of the structure of C2H2 adsorbed on Si(001)

Theoretical studies of the structure of C 2H 2 adsorbed on Si(0 0 1)

Structure Character in Small-Carbon-cluster Deposition on Diamond Surface

Low energy deposition of size-selected Si clusters onto graphite

Memory Effect in the Deposition of C-20 Fullerenes on A Diamond Surface

Deposition of Small Clusters on Surface: A Molecular Dynamics Simulation

Molecular Dynamics Simulation of Energetic Cu55 Clusters Deposition on a Fe (0 0 1) Surface

Impact-induced Chemisorption of C2H2 on Diamond(001) Surfaces: a Molecular Dynamics Simulation

The Effect of Incident Energy and Cluster Size on the Deposition Cu Clusters on the Fe (001) Surface

Effect of Charge Transfer on the Geometric Structure of a C70monolayer on the Surface of Ag(111)

Carbon superatom thin films

Adsorption of a Carbon Atom on the Ni 38 Magic Cluster and Three Low-Index Nickel Surfaces: A Comparative First-Principles Study

Ejection of clusters from solid surface by impact of size-selected cluster ion

Deposition and diffusion of size-selected (Ag400+) clusters on a stepped graphite surface

Chemisorption Ofc28fullerene Onc(4×4)reconstructed GaAs(001) Surface: A Density Functional Theory Study

Exploring the Adsorption Mechanism of Tetracene on Ag(110) by STM and Dispersion-Corrected DFT

Aggregation of Organic Molecules on Silver Surface with the Balance Between Molecule–substrate Interaction and Intermolecular Interaction

Tight-binding molecular-dynamics study of amorphous carbon deposits over silicon surfaces

The Structural and Vibrational Properties of C-60 Adsorbed on the Graphite(0001) Surface

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