Reactive force field potential for carbon deposition on silicon surfaces

In this paper a new interatomic potential based on the Kieffer force field and designed to perform molecular dynamics (MD) simulations of carbon deposition on silicon surfaces is implemented. This potential is a third-order reactive force field that includes a dynamic charge transfer and allows for the formation and breaking of bonds. The parameters for Si-C and C-C interactions are optimized using a genetic algorithm. The quality of the potential is tested on its ability to model silicon carbide and diamond physical properties as well as the formation energies of point defects. Furthermore, MD simulations of carbon deposition on reconstructed (100) silicon surfaces are carried out and compared to similar simulations using a Tersoff-like bond order potential. Simulations with both potentials produce similar results showing the ability to extend the use of the Kieffer potential to deposition studies. The investigation reveals the presence of a channelling effect when depositing the carbon at 45° incidence angle. This effect is due to channels running in directions symmetrically equivalent to the (110) direction. The channelling is observed to a lesser extent for carbon atoms with 30° and 60° incidence angles relative to the surface normal. On a pristine silicon surface, sticking coefficients were found to vary between 100 and 73%, depending on deposition conditions.

[1]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[2]  M. Shapiro,et al.  The influence of the ion-atom potential on molecular dynamics simulations of sputtering , 2004 .

[3]  Theoretical studies of SiC, AlN and their (110) surfaces , 2009 .

[4]  T. Wirtz,et al.  The Storing Matter technique: application to polymer samples using Ag collectors , 2011 .

[5]  M. Posselt,et al.  Improvement of the repulsive part of the classical interatomic potential for SiC , 2003 .

[6]  T. Wirtz,et al.  Storing Matter: A new quantitative and sensitive analytical technique , 2008 .

[7]  K. Nordlund,et al.  Dynamics of cluster induced sputtering in gold , 2007 .

[8]  Methfessel,et al.  Calculated elastic constants and deformation potentials of cubic SiC. , 1991, Physical review. B, Condensed matter.

[9]  Gerhard Hobler,et al.  On the useful range of application of molecular dynamics simulations in the recoil interaction approximation , 2001 .

[10]  Eric J. Bylaska,et al.  Ab Initio and Empirical Potential Studies of Defect Properties in 3C-SiC , 2001 .

[11]  Pantelides,et al.  Mechanism of self-diffusion in diamond. , 1988, Physical review letters.

[12]  John Kieffer,et al.  Molecular dynamics study of cristobalite silica using a charge transfer three-body potential: Phase transformation and structural disorder , 2003 .

[13]  R. Nieminen,et al.  Density-functional calculations of defect formation energies using supercell methods: Defects in diamond , 2005 .

[14]  I. Lyashenko,et al.  Statistical theory of the boundary friction of atomically flat solid surfaces in the presence of a lubricant layer , 2012 .

[15]  Liping Huang,et al.  Thermomechanical anomalies and polyamorphism in B 2 O 3 glass: A molecular dynamics simulation study , 2006 .

[16]  D. Sánchez-Portal,et al.  The SIESTA method for ab initio order-N materials simulation , 2001, cond-mat/0111138.

[17]  Kai Nordlund,et al.  Modelling of compound semiconductors: analytical bond-order potential for gallium, nitrogen and gallium nitride , 2003 .

[18]  W. J. Weber,et al.  Empirical potential approach for defect properties in 3C-SiC , 2002 .

[19]  Martins,et al.  Efficient pseudopotentials for plane-wave calculations. , 1991, Physical review. B, Condensed matter.

[20]  S. Lucas,et al.  Simulation at high temperature of atomic deposition, islands coalescence, Ostwald and inverse Ostwald ripening with a general simple kinetic Monte Carlo code , 2010 .

[21]  F. Priolo,et al.  CHANNELING IMPLANTS IN SILICON CRYSTALS , 1994 .

[22]  J. Kieffer,et al.  Low-energy oxygen bombardment of silicon by MD simulations making use of a reactive force field , 2011 .

[23]  K. Nordlund.  Computational materials science of ion irradiation , 2002 .

[24]  C. Mansilla,et al.  Application of the Storing Matter technique to the analysis of semiconductor materials , 2010 .

[25]  P. Erhart,et al.  Analytical potential for atomistic simulations of silicon, carbon, and silicon carbide , 2005 .

[26]  C. Mansilla,et al.  Storing Matter: A new quantitative and sensitive analytical technique based on sputtering and collection of sample material , 2009 .

[27]  A. K. Ramdas,et al.  Brillouin scattering in diamond , 1975 .

[28]  Tersoff.  Carbon defects and defect reactions in silicon. , 1990, Physical review letters.

[29]  Rafael Ramírez,et al.  Structural and thermodynamic properties of diamond: A path-integral Monte Carlo study , 2000 .

[30]  Ji‐an Xu,et al.  Total energy calculations of the lattice properties of cubic and hexagonal diamond , 1998 .

[31]  Liping Huang,et al.  Transformation pathways of silica under high pressure , 2006, Nature materials.

[32]  L. Angibaud,et al.  Parameter optimization in molecular dynamics simulations using a genetic algorithm , 2011 .

[33]  D. Sánchez-Portal,et al.  Numerical atomic orbitals for linear-scaling calculations , 2001, cond-mat/0104170.

[34]  M. Krisch,et al.  Phonon density of states probed by inelastic x-ray scattering , 2005 .

[35]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[36]  Kai Nordlund,et al.  Molecular dynamics simulation of ion ranges in the 1–100 keV energy range , 1995 .

[37]  Mark T. Robinson,et al.  Computer simulation of atomic-displacement cascades in solids in the binary-collision approximation , 1974 .