posic/publications/c_defects_in_si.tex

   1 \documentclass[prb,twocolumn,superscriptaddress,a4paper,showkeys,showpacs]{revtex4}\r
   2 \usepackage{graphicx}\r
   3 \usepackage{subfigure}\r
   4 \usepackage{dcolumn}\r
   5 \usepackage{booktabs}\r
   6 \usepackage{units}\r
   7 \usepackage{amsmath}\r
   8 \usepackage{amsfonts}\r
   9 \usepackage{amssymb}\r
  10 \r
  11 \r
  12 \begin{document}\r
  13 \r
  14 \title{Description of Defects in Carbon implanted Silicon -- a comparison of classical potentials and first principles methods}\r
  15 \author{F. Zirkelbach} \author{B. Stritzker}\r
  16 \affiliation{Experimentalphysik IV, Universit\"at Augsburg, D-86153 Augsburg, Germany}\r
  17 \author{K. Nordlund}\r
  18 \affiliation{Accelerator Laboratory, University of Helsinki, 00014 Helsinki, Finland}\r
  19 \author{J. K. N. Lindner}\r
  20 \affiliation{Experimentelle Physik, Universit\"at Paderborn, 33095 Paderborn, Germany}\r
  21 \author{W. G. Schmidt} \author{E. Rauls}\r
  22 \affiliation{Theoretische Physik, Universit\"at Paderborn, 33095 Paderborn, Germany}\r
  23 \r
  24 \begin{abstract}\r
  25 We present a comparative theoretical investigation of carbon interstitials in silicon.\r
  26 Calculations using classical potentials are put aside first principles density functional theory calculations of the geometries, formation and activation energies of the carbon dumbbell interstitial, showing the importance of a quantum mechanical description of this system.\r
  27 \end{abstract}\r
  28 \r
  29 \keywords{point defects, migration, interstitials, first principles calculations, classical potentials }\r
  30 \pacs{ find out later... }\r
  31 \maketitle\r
  32 \r
  33 %  --------------------------------------------------------------------------------\r
  34 \section{Introduction}\r
  35 \r
  36 %Frank:  Idea: description of 3C-SiC-precipitation in C-implanted silicon.\\ \r
  37 %        cite and describe briefly experimental work - why is this material important/better than other SiC).  \\ \r
  38 %        Describe the precipitation process in brief.\\ \r
  39 %        Sum up literature where classical potentials have been used (a) successful, and (b) failed. Also add citations of Nordlunds papers. Not only on silicon or SiC!\\\r
  40 \r
  41 % there should be a short motivation for the material system!\r
  42 Silicon carbide (SiC) has a number of remarkable physical and chemical properties.\r
  43 The wide band gap semiconductor (2.3 eV - 3.3 eV) exhibiting a high breakdown field, saturated electron drift velocity and thermal conductivity in conjunction with its unique thermal and mechanical stability as well as radiation hardness is a suitable material for high-temperature, high-frequency and high-power devices\cite{wesch96}, which are moreover deployable in harsh environments\cite{capano97}.\r
  44 % there are different polytpes with different properties and 3c-sic in special\r
  45 SiC, which forms fourfold coordinated covalent bonds, tends to crystallize into many different modifications, which solely differ in the one-dimensional stacking sequence of identical, close-packed SiC bilayers\cite{fischer90}.\r
  46 Different polytypes exhibit different properties, in which the cubic phase (3C-SiC) shows increased values for the thermal conductivity and breakdown field compared to other polytypes\cite{wesch96}, which is thus most effective for high-performance electronic devices.\r
  47 \r
  48 % (thin films of) 3c-sic can be produced by ibs\r
  49 Next to 3C-SiC epitaxial layer growth by chemical vapor deposition (CVD) on 6H-SiC substrates\cite{powell90} or molecular beam epitaxy on silicon\cite{} as well as 6H-SiC\cite{}, ion beam synthesis (IBS) constitutes a promising method to produce 3C-SiC epitaxial layers of high quality in crystalline silicon\cite{lindner02} (c-Si).\r
  50 Highly energetic carbon ions are ...\r
  51 \r
  52 The relevant structures are with $\approx$ 20000 atoms/nanocrystal way too large to be completely be described with high accuracy \r
  53 quantum mechanical methods. Modelling the processes described above require the use of less accurate methods, like e.g. classical \r
  54 potentials (Erhard/Albe\cite{albe},Stillinger-Weber\cite{stillinger},...).  Whether such potentials are appropriate for the description of the \r
  55 physical problem has, however, to be verified first by applying both methods to relevant processes that can be treated by both methods.     \r
  56 In this work, we have implemented and compared the applicability of several (?) classical potentials to ab initio results for one \r
  57 of the most important processes of our original question.   \r
  58 \r
  59 In the following we will present a comparative investigation of density functional theory (DFT) studies and \r
  60 classical potential calculations of the structure, energetics and mobility of carbon defects in silicon. \r
  61 \r
  62 %  --------------------------------------------------------------------------------\r
  63 \section{Methodology}\r
  64 % ----- DFT ------\r
  65 The first-principles DFT calculations were performed with the plane-wave based \r
  66 Vienna Ab-initio Simulation Package (VASP)\cite{kresse96}. The Kohn-Sham equations were solved \r
  67 using the generalized-gradient XC-functional approximation proposed by Perdew and \r
  68 Wang (GGA-PW91)\cite{perdew92}. \r
  69 The electron-ion interaction was described by the projector-augmented wave (PAW) method\cite{bloechel94,kresse99}. \r
  70 In the PAW data scalar relativistic corrections are contained. Throughout this work an \r
  71 energy cut-off of \unit[300]{eV} was used to expand the wave functions into the plane-wave basis. \r
  72 For the sampling of the Brillouin zone, only the $\Gamma$-point was used. \r
  73 The defect structures and the migration paths were modelled in cubic supercells containing 216 Si-atoms. Spin polarization has been \r
  74 fully accounted for.    \r
  75 \r
  76 % ------ Albe potential ---------\r
  77 %% Frank:  Setup/short description of the potential ?\r
  78 For the calculations with the classical potentials... \r
  79 \r
  80 \r
  81 \section{Results}\r
  82 \r
  83 After ion implantation, carbon interstitials are the most common defects in the silicon sample. Their mobility is the \r
  84 crucial quantity to be investigated.   We thus started our comparative investigations by comparing the stability and the \r
  85 mobility of an isolated carbon interstitial in silicon bulk in the various possible geometries it can  take.  \r
  86 \r
  87 \subsection{Carbon interstitials in various geometries}\r
  88 \r
  89 Several geometries have been calculated to be stable for the carbon interstitial. Fig.\ref{fig:interstitials} shows all these \r
  90 structures. However, there are some discrepancies between the results from classical potential calculations and those obtained \r
  91 from first principles.    \r
  92 Table \ref{table:formation} summarizes the formation energies of the interstitial geometries for both methods used in this work \r
  93 and compares the results to literature values. (...check references for more data, ..) \r
  94 \r
  95 % Tables:  like in the talk, but add further literature data and give the references/citations (also to bibliography \r
  96 % at the end!) \r
  97 %\begin{figure}\r
  98 %\includegraphics[width=1.0\columnwidth]{models.eps}\r
  99 %\caption{\label{fig:interstitials} Molecular model of the possible carbon interstitials. }\r
 100 %\end{figure} \r
 101 \r
 102 While the Albe potential predicts ... as stable, DFT does not. ...(further comparisons, trend "too high/low" E-formation,...)...    \r
 103  Nevertheless, both methods predict the (110) dumb bell configuration to  be the most stable...  (?) \r
 104 \r
 105 \r
 106 \subsection{Mobility}\r
 107 A measure for the mobility of the interstitial carbon is the activation energy for the migration path from one stable \r
 108 position to another. The stable defect geometries have been discussed in the previous subsection.   We now investigate \r
 109 the migration from the most stable structure (...should be named somehow...) on one site of the silicon host lattice to \r
 110 a neighbored site.  \r
 111 On the lowest energy path (first principles), the carbon atom starts to move along (110)..(check that!)...   The center of the line connecting \r
 112 initial and final structure has been found to be a local minimum and not a saddle point as could be expected. The two \r
 113 saddle points shortly before and behind this local minimum are slightly displaced out of the (110) plane by ... {\AA}.  ..(check that!)..\r
 114 This path is not surprising -- a similar behavior was e.g. found earlier for the carbon split interstitial \cite{rauls03a} and the phosphorus \r
 115 interstitial \cite{rauls03b,gerstmann03} in SiC.  However, an interesting effect is the change of the spin state from zero at the (110) dumb bell \r
 116 configuration to one at the local minimum.  By this, the energy of the local minimum is lowered by 0.3 eV  (... check it!!..). \r
 117 %\begin{figure}\r
 118 %\includegraphics[width=1.0\columnwidth]{path-DFT.eps}\r
 119 %\caption{\label{fig:path-DFT} Energy of the carbon interstitial during migration from ... to ... calculated from first principles. The \r
 120 %         activation energy of 0.9 eV (?) agrees well with experimental findings (0.7-0.9 eV?). }\r
 121 %\end{figure} \r
 122 Fig.\ref{fig:path-DFT} shows the energy along this lowest energy migration path.  The activation energy of 0.9 eV (?) agrees well \r
 123 with experimental findings (0.7-0.9 eV?).\r
 124 \r
 125 Calculations with the Albe potential yield a different picture.  \r
 126 %\begin{figure}\r
 127 %\includegraphics[width=1.0\columnwidth]{path-Albe.eps}\r
 128 %\caption{\label{fig:path-Albe} Energy of the carbon interstitial during migration from ... to ... calculated using the classical potential \r
 129 %          method. Here, the activation energy is 2.2 eV (?). }\r
 130 %\end{figure} \r
 131 Fig.\ref{fig:path-Albe} shows the energy along the lowest energy migration path found by this method. The activation energy of 2.2 eV (?) \r
 132 is way too high to explain the experimental findings (0.7-0.9 eV?).  (...further discussion...) \r
 133 \r
 134 \section{Discussion}\r
 135 The first principles results are in good agreement to previous work on this subject \cite{joannopoulos,xyz}  (...add some references!...).   With an \r
 136 activation energy of 0.9 eV, the carbon interstitial can be expected to be mobile at temperatures in the range of... (?).    \r
 137 The description of the same processes obviously fails if we use the classical potential method.  \r
 138 Already the geometry of the most stable dumb bell configuration differs considerably from that of the first principles calculated \r
 139 structure.  (..... add description, the two main angles and bond lengths and an explanation...)\r
 140 Formation energies are throughout too high... (...reason?...)\r
 141 \r
 142 A reason for this failure of the classical description is most likely...  (cut-off, neglect of quantum mechanical effects,...)\r
 143 \r
 144 \r
 145 \section{Summary}\r
 146 In summary, we have shown that ab initio calculations are very close to the results expected from experimental data.  \r
 147 Furthermore, they agree well with other theoretical results. (...some results - later...)  The classical potentials, however, fail to describe the \r
 148 selected processes. This has been shown to have two reasons, i.e. the artificial cut-off of the next nearest neighbor \r
 149 interaction on the one hand, on the other hand the quantum mechanical effects which are crucial in the problem under study. \r
 150 \r
 151 % ----------------------------------------------------\r
 152 \section*{Acknowledgment}\r
 153 The calculations were done using grants of computer time from the \r
 154 Paderborn Center for Parallel Computing (PC$^2$) and the \r
 155 H\"ochstleistungs-Rechenzentrum Stuttgart. The Deutsche \r
 156 Forschungsgemeinschaft is acknowledged for financial support.\r
 157 \r
 158 % --------------------------------- references -------------------\r
 159 \begin{thebibliography}{99}\r
 160 \bibitem{kresse96} G. Kresse and J. Furthm\"uller, \r
 161   Comput. Mater. Sci. {\bf 6}, 15 (1996). \r
 162 \bibitem{perdew92} J. P. Perdew, J. A. Chevary, S. H. Vosko, K. A. Jackson, M. R. Pederson, D. J. Singh and C. Fiolhais, \r
 163   Phys. Rev. B {\bf 46}, 6671 (1992). \r
 164 \bibitem{ceperley80} D. M. Ceperley and B. J. Alder, \r
 165   Phys. Rev. Lett. {\bf 45}, 556 (1980). \r
 166 \bibitem{perdew81} J. P. Perdew and A. Zunger, \r
 167   Phys. Rev. B {\bf 23}, 5048 (1981). \r
 168 \bibitem{bloechel94} P. E. Bl\"ochl, \r
 169   Phys. Rev. B {\bf 50}, 17953 (1994).\r
 170 \bibitem{kresse99} G. Kresse and D. Joubert, \r
 171   Phys. Rev. B {\bf 59}, 1758 (1999).\r
 172 \bibitem{monk76} H. J. Monkhorst and J. D. Pack, \r
 173   Phys. Rev. B {\bf 13}, 5188 (1976). \r
 174 \bibitem{albe}  Albe potential\r
 175 \bibitem{stillinger} Stillinger-Weber potential \r
 176 \bibitem{joannopoulos} Joannopoulos\r
 177 \bibitem{xyz} who else?  \r
 178 \bibitem{rauls03a} E. Rauls, A. Gali, P. De´ak, and Th. Frauenheim, Phys. Rev. B, 68, 155208 (2003).\r
 179 \bibitem{rauls03b} E. Rauls, U. Gerstmann, H. Overhof, and Th. Frauenheim, Physica B, Vols. 340-342, p. 184-189 (2003). \r
 180 \bibitem{gerstmann03} U. Gerstmann, E. Rauls, Th. Frauenheim, and H. Overhof, Phys. Rev. B, 67, 205202, (2003).\r
 181 \r
 182 \end{thebibliography}\r
 183 \r
 184 \r
 185 \end{document}\r
 186 \r