]> hackdaworld.org Git - lectures/latex.git/blob - posic/publications/emrs2012.tex
removed openstack, i didn't do much!
[lectures/latex.git] / posic / publications / emrs2012.tex
1 \documentclass[pss]{wiley2sp}
2 \usepackage{graphicx}
3 \usepackage{units}
4 \usepackage{amsmath}
5 \usepackage{amsfonts}
6 \usepackage{amssymb}
7
8 % additional stuff
9 \usepackage{miller}
10
11 % roman numbers
12 \newcommand{\RM}[1]{\MakeUppercase{\romannumeral #1{}}}
13
14 \begin{document}
15
16 \title{First-principles and empirical potential simulation study of intrinsic
17        and carbon-related defects in silicon}
18
19 \titlerunning{First-principles and empirical potential simulation study
20               of intrinsic and carbon-related defects in silicon}
21
22 \author{%
23  F. Zirkelbach\textsuperscript{\Ast,\textsf{\bfseries 1}},
24  B. Stritzker\textsuperscript{\textsf{\bfseries 1}},
25  K. Nordlund\textsuperscript{\textsf{\bfseries 2}},
26  W. G. Schmidt\textsuperscript{\textsf{\bfseries 3}},
27  E. Rauls\textsuperscript{\textsf{\bfseries 3}},
28  J. K. N. Lindner\textsuperscript{\textsf{\bfseries 3}}
29 }
30
31 \authorrunning{F. Zirkelbach et al.}
32
33 \mail{e-mail
34   \textsf{frank.zirkelbach@physik.uni-augsburg.de}
35 }
36
37 \institute{%
38   \textsuperscript{1}\, Experimentalphysik IV, Universit\"at Augsburg,
39                         86135 Augsburg, Germany\\
40   \textsuperscript{2}\, Department of Physics, University of Helsinki,
41                         00014 Helsinki, Finland\\
42   \textsuperscript{3}\, Department Physik, Universit\"at Paderborn,
43                         33095 Paderborn, Germany}
44
45 \received{XXXX, revised XXXX, accepted XXXX}
46
47 \published{XXXX}
48
49 \keywords{Silicon, carbon, silicon carbide, defect formation, defect migration,
50           density functional theory, empirical potential, molecular dynamics.}
51 %\pacs{61.72.J-,61.72.Yx,61.72.uj,66.30.J-,79.20.Rf,31.15.A-}
52
53 \abstract{%
54 Results of atomistic simulations aimed at understanding precipitation of the highly attractive wide band gap semiconductor material silicon carbide in silicon are presented.
55 The study involves a systematic investigation of intrinsic and carbon-related defects as well as defect combinations and defect migration by both, quantum-mechanical first-principles as well as empirical potential methods.
56 Comparing formation and activation energies, ground-state structures of defects and defect combinations as well as energetically favorable agglomeration of defects are predicted.
57 Moreover, accurate {\em ab initio} calculations unveil limitations of the analytical method based on a Tersoff-like bond order potential.
58 A work-around is proposed in order to subsequently apply the highly efficient technique on large structures not accessible by first-principles methods.
59 The outcome of both types of simulation provides a basic microscopic understanding of defect formation and structural evolution particularly at non-equilibrium conditions strongly deviated from the ground state as commonly found in SiC growth processes.
60 A possible precipitation mechanism, which conforms well to experimental findings and clarifies contradictory views present in the literature is outlined.
61 }
62
63 \maketitle
64
65 \section{Introduction}
66
67 Silicon carbide (SiC) is a promising material for high-temperature, high-power and high-frequency electronic and optoelectronic devices, which can operate under extreme conditions
68 % shorten
69 % \cite{edgar92,morkoc94,wesch96,capano97,park98}.
70 \cite{edgar92,capano97,park98}.
71 Ion beam synthesis (IBS) consisting of high-dose carbon implantation into crystalline silicon (c-Si) and subsequent or in situ annealing is a promising technique to fabricate nano-sized precipitates and thin films of the favorable cubic SiC (3C-SiC) polytype topotactically aligned to and embedded in the silicon host
72 % shorten
73 % \cite{borders71,lindner99,lindner01,lindner02}.
74 \cite{borders71,lindner01}.
75 However, the process of formation of SiC precipitates in Si during C implantation is not yet fully understood and controversial ideas exist in the literature.
76 Based on experimental high resolution transmission electron microscopy (HREM) studies 
77 % shorten
78 % \cite{werner96,werner97,eichhorn99,lindner99_2,koegler03}
79 \cite{werner96,lindner99_2,koegler03}
80 it is assumed that incorporated C atoms form C-Si dimers (dumbbells) on regular Si lattice sites.
81 The highly mobile C interstitials agglomerate into large clusters followed by the formation of incoherent 3C-SiC nanocrystallites once a critical size of the cluster is reached.
82 In contrast, a couple of other studies \cite{strane94,nejim95,serre95} suggest initial coherent SiC formation by agglomeration of substitutional instead of interstitial C followed by the loss of coherency once the increasing strain energy surpasses the interfacial energy of the incoherent 3C-SiC precipitate and the c-Si substrate.
83
84 To solve this controversy and in order to understand the effective underlying processes on a microscopic level atomistic simulations are performed.
85 % ????
86 A lot of theoretical work has been done on intrinsic point defects in Si
87 % shorten
88 % \cite{bar-yam84,bar-yam84_2,car84,batra87,bloechl93,tang97,leung99,colombo02,goedecker02,al-mushadani03,hobler05,sahli05,posselt08,ma10} 
89 \cite{bar-yam84,car84,bloechl93,tang97,leung99,al-mushadani03,hobler05,sahli05,posselt08,ma10} 
90 and C defects and defect reactions in Si 
91 % shorten
92 %\cite{tersoff90,dal_pino93,capaz94,burnard93,leary97,capaz98,zhu98,mattoni2002,park02}.
93 \cite{tersoff90,dal_pino93,capaz94,burnard93,leary97,capaz98,mattoni2002}.
94 However, none of the mentioned studies consistently investigates entirely the relevant defect structures and reactions concentrated on the specific problem of 3C-SiC formation in C implanted Si.
95 % ????
96
97 In the present study, an accurate first-principles treatment is utilized to systematically investigate relevant intrinsic as well as carbon related defect structures and defect mobilities in silicon, which allow to draw conclusions on the mechanism of SiC precipitation in Si.
98 These findings are compared to empirical potential results, which, by taking into account the drawbacks of the less accurate though computationally efficient method enabling molecular dynamics (MD) simulations of large structures, support and complete previous findings on SiC precipitation based on the quantum-mechanical treatment.
99
100 \section{Methodology}
101
102 The plane-wave based Vienna {\em ab initio} simulation package (VASP) \cite{kresse96} is used for the first-principles calculations based on density functional theory (DFT).
103 Exchange and correlation is taken into account by the generalized-gradient approximation \cite{perdew86,perdew92}.
104 Norm-con\-ser\-ving ultra-soft pseudopotentials \cite{hamann79} as implemented in VASP \cite{vanderbilt90} are used to describe the electron-ion interaction.
105 A kinetic energy cut-off of \unit[300]{eV} is employed.
106 Defect structures and migration paths were modelled in cubic supercells with a side length of \unit[1.6]{nm} containing $216$ Si atoms.
107 These structures are large enough to restrict sampling of the Brillouin zone to the $\Gamma$-point and formation energies and structures are reasonably converged.
108 The ions and cell shape are allowed to change in order to realize a constant pressure simulation realized by the conjugate gradient algorithm.
109 Spin polarization is fully accounted for.
110
111 Migration and recombination pathways are investigated utilizing the constraint conjugate gradient relaxation technique (CRT) \cite{kaukonen98}.
112 The defect formation energy $E-N_{\text{Si}}\mu_{\text{Si}}-N_{\text{C}}\mu_{\text{C}}$ is defined by choosing SiC as a particle reservoir for the C impurity, i.e. the chemical potentials are determined by the cohesive energies of a perfect Si and SiC supercell after ionic relaxation.
113 The binding energy of a defect pair is given by the difference of the formation energy of the complex and the sum of the two separated defect configurations.
114 Accordingly, energetically favorable configurations result in binding energies below zero while unfavorable configurations show positive values for the binding energy.
115 The interaction strength, i.e. the absolute value of the binding energy, approaches zero for increasingly non-interacting isolated defects.
116
117 Within the empirical approach, defect structures are modeled in a supercell of nine Si lattice constants in each direction consisting of 5832 Si atoms.
118 Reproducing SiC precipitation is attempted by successive insertion of 6000 C atoms to form a minimal 3C-SiC precipitate with a radius of about \unit[3.1]{nm} within the Si host consisting of 31 unit cells (238328 atoms) in each direction.
119 At constant temperature 10 atoms are inserted on statistically distributed positions at a time.
120 A Tersoff-like bond order potential by Erhart and Albe (EA) \cite{albe_sic_pot} has been utilized, which accounts for nearest neighbor interactions realized by a cut-off function dropping the interaction to zero in between the first and second nearest neighbor distance.
121 The Berendsen barostat and thermostat \cite{berendsen84} with a time constant of \unit[100]{fs} enables the isothermal-isobaric ensemble.
122 The velocity Verlet algorithm \cite{verlet67} and a fixed time step of \unit[1]{fs} is used to integrate the equations motion.
123 Structural relaxation of defect structures is treated by the same algorithms at zero temperature.
124
125 \section{Defect configurations in silicon}
126
127 Table~\ref{tab:defects} summarizes the formation energies of relevant defect structures for the EA and DFT calculations, which are shown in Fig.~\ref{fig:intrinsic_def} and \ref{fig:carbon_def}.
128 \begin{table*}
129 \centering
130 \begin{tabular}{l c c c c c c c c c}
131 \hline
132 $E_{\text{f}}$ [eV] & Si$_{\text{i}}$ \hkl<1 1 0> DB & Si$_{\text{i}}$ H & Si$_{\text{i}}$ T & Si$_{\text{i}}$ \hkl<1 0 0> DB & V & C$_{\text{s}}$ & C$_{\text{i}}$ \hkl<1 0 0> DB & C$_{\text{i}}$ \hkl<1 1 0> DB & C$_{\text{i}}$ BC \\
133 \hline
134 VASP & 3.39 & 3.42 & 3.77 & 4.41 & 3.63 & 1.95 & 3.72 & 4.16 & 4.66 \\
135 Erhart/Albe & 4.39 & 4.48$^*$ & 3.40 & 5.42 & 3.13 & 0.75 & 3.88 & 5.18 & 5.59$^*$ \\
136 \hline
137 \end{tabular}
138 \caption{Formation energies of C and Si point defects in c-Si given in eV. T denotes the tetrahedral, H the hexagonal and BC the bond-centered interstitial configuration. V corresponds to the vacancy. Subscript i and s indicates the interstitial and substitutional configuration. Dumbbell configurations are abbreviated by DB. Formation energies of unstable configurations are marked by an asterisk.}
139 \label{tab:defects}
140 \end{table*}
141 \begin{figure}
142 \subfloat[Intrinsic Si point defects.]{%
143 \begin{minipage}{0.9\columnwidth}
144 \centering
145 \begin{minipage}[t]{0.43\columnwidth}
146 \centering
147 \underline{Si$_{\text{i}}$ \hkl<1 1 0> DB}\\
148 \includegraphics[width=0.8\columnwidth]{si110_bonds.eps}
149 \end{minipage}
150 \begin{minipage}[t]{0.43\columnwidth}
151 \centering
152 \underline{Si$_{\text{i}}$ hexagonal}\\
153 \includegraphics[width=0.8\columnwidth]{sihex_bonds.eps}
154 \end{minipage}\\
155 \begin{minipage}[t]{0.43\columnwidth}
156 \centering
157 \underline{Si$_{\text{i}}$ tetrahedral}\\
158 \includegraphics[width=0.8\columnwidth]{sitet_bonds.eps}
159 \end{minipage}
160 \begin{minipage}[t]{0.43\columnwidth}
161 \centering
162 \underline{Si$_{\text{i}}$ \hkl<1 0 0> DB}\\
163 \includegraphics[width=0.8\columnwidth]{si100_bonds.eps}
164 \end{minipage}
165 %\caption{Configurations of intrinsic Si point defects. Dumbbell configurations are abbreviated by DB.}
166 \end{minipage}
167 \label{fig:intrinsic_def}
168 }\\
169 \subfloat[C point defects in Si.]{%
170 \begin{minipage}{0.9\columnwidth}
171 \centering
172 \begin{minipage}[t]{0.43\columnwidth}
173 \centering
174 \underline{C$_{\text{s}}$}\\
175 \includegraphics[width=0.8\columnwidth]{csub_bonds.eps}
176 \end{minipage}
177 \begin{minipage}[t]{0.43\columnwidth}
178 \centering
179 \underline{C$_{\text{i}}$ \hkl<1 0 0> DB}\\
180 \includegraphics[width=0.8\columnwidth]{c100_bonds.eps}
181 \end{minipage}\\
182 \begin{minipage}[t]{0.43\columnwidth}
183 \centering
184 \underline{C$_{\text{i}}$ \hkl<1 1 0> DB}\\
185 \includegraphics[width=0.8\columnwidth]{c110_bonds.eps}
186 \end{minipage}
187 \begin{minipage}[t]{0.43\columnwidth}
188 \centering
189 \underline{C$_{\text{i}}$ bond-centered}\\
190 \includegraphics[width=0.8\columnwidth]{cbc_bonds.eps}
191 \end{minipage}
192 %\caption{Configurations of carbon point defects in silicon. Silicon and carbon atoms are illustrated by yellow and gray spheres respectively. Dumbbell configurations are abbreviated by DB.}
193 \end{minipage}
194 \label{fig:carbon_def}
195 }
196 \caption{Defect configurations in Si. Si and C atoms are illustrated by yellow and gray spheres respectively. Dumbbell configurations are abbreviated by DB.}
197 \end{figure}
198
199 Regarding intrinsic defects in Si, classical potential and {\em ab initio} methods predict energies of formation that are within the same order of magnitude.
200 The EA potential does not reproduce the correct ground state, i.e. the interstitial Si (Si$_{\text{i}}$) \hkl<1 1 0> dumbbell (DB), which is consensus for Si$_{\text{i}}$ \cite{leung99,al-mushadani03}.
201 Instead, the tetrahedral configuration is favored, a limitation assumed to arise due to the sharp cut-off as has already been discussed by Tersoff \cite{tersoff90}.
202
203 In the case of C impurities, although discrepancies exist, classical potential and first-principles methods depict the correct order of the formation energies.
204 Next to the substitutional C (C$_{\text{s}}$) configuration, which is not an interstitial configuration since the C atom occupies an already vacant Si lattice site, the interstitial C (C$_{\text{i}}$) \hkl<1 0 0> DB constitutes the energetically most favorable configuration, in which the C and Si dumbbell atoms share a regular Si lattice site.
205 This finding is in agreement with several theoretical \cite{dal_pino93,capaz94,burnard93,leary97} and experimental \cite{watkins76,song90} investigations, which all predict this configuration to be the ground state.
206 It is worth to note that the bond-centered (BC) configuration constitutes a real local minimum in spin polarized calculations in contrast to results \cite{capaz94} without spin predicting a saddle point configuration as well as to the empirical description, which shows a relaxation into the C$_{\text{i}}$ \hkl<1 0 0> DB ground-state configuration.
207
208 \section{Mobility of the carbon defect}
209
210 The migration barriers of the ground-state C defect are investigated by both, first-principles as well as the empirical method.
211 The migration pathways are shown in Fig.~\ref{fig:mig}.
212
213 \begin{figure}
214 \subfloat[Transition path obtained by first-principles methods.]{%
215 \includegraphics[width=\columnwidth]{path2_vasp_s.ps}
216 }\\
217 %\caption{Migration barrier and structures of the C$_{\text{i}}$ \hkl[0 0 -1] DB (left) to the \hkl[0 -1 0] DB (right) transition as obtained by first principles methods.}
218 \subfloat[Transition involving the {\hkl[1 1 0]} DB (center) configuration within the EA description.]{%
219 \includegraphics[width=\columnwidth]{110mig.ps}
220 }
221 \caption{Migration barriers and structures of the C$_{\text{i}}$ \hkl[0 0 -1] DB (left) to the hkl[0 -1 0] DB (right) transition.}
222 %\caption{Migration barrier and structures of the C$_{\text{i}}$ \hkl[0 0 -1] DB (left) to the hkl[0 -1 0] DB (right) transition involving the \hkl[1 1 0] DB (center) configuration within EA description. Migration simulations were performed utilizing time constants of \unit[1]{fs} (solid line) and \unit[100]{fs} (dashed line) for the Berendsen thermostat.}
223 \label{fig:mig}
224 \end{figure}
225
226 In qualitative agreement with the results of Capaz~et~al.\  \cite{capaz94}, the lowest migration barrier of the ground-state C$_{\text{i}}$ defect within the quantum-mechanical treatment is found for the path, in which a C$_{\text{i}}$ \hkl[0 0 -1] DB migrates to a C$_{\text{i}}$ \hkl[0 -1 0] DB located at the neighbored Si lattice site in \hkl[1 1 -1] direction.
227 Calculations in this work reinforce this path by an additional improvement of the quantitative conformance of the barrier height of \unit[0.90]{eV} to experimental values (\unit[0.70-0.87]{eV}) \cite{song90,lindner06,tipping87}.
228
229 In contrast, the empirical approach does not reproduce the same path.
230 Related to the above mentioned instability of the BC configuration, a pathway involving the C$_{\text{i}}$ \hkl<1 1 0> DB as an intermediate configuration must be considered most plausible \cite{zirkelbach11}.
231 Considering a two step diffusion process and assuming equal preexponential factors, a total effective migration barrier 3.5 times higher than the one obtained by first-principles methods is obtained.
232 A more detailed description can be found in previous studies \cite{zirkelbach11,zirkelbach10}.
233
234 \section{Defect combinations}
235
236 The implantation of highly energetic C atoms results in a multiplicity of possible point defects and respective combinations.
237 Thus, defect combinations of an initial C$_{\text{i}}$ DB and further types of defects created at certain neighbor positions have been investigated \cite{zirkelbach11} exclusively by DFT calculations.
238 Some of the most important results are presented in the following.
239
240 \begin{figure}
241 \includegraphics[width=\columnwidth]{db_along_110_cc_n.ps}
242 \caption{Minimum binding energy of dumbbell combinations separated along \hkl[1 1 0] with respect to the C-C distance. The blue line is a guide for the eye and the green curve corresponds to the most suitable fit function consisting of all but the first data point.}
243 \label{fig:dc_110}
244 \end{figure}
245
246 The agglomeration of interstitial C is found to be energetically favorable.
247 As can be seen in Fig.~\ref{fig:dc_110}, which shows the binding energies of DB combinations separated along the \hkl[1 1 0] chain, a capture radius clearly exceeding \unit[1]{nm} is observed.
248 The interaction is proportional to the reciprocal cube of the C-C distance for extended separations of the defects.
249 However, the interpolated graph shows the disappearance of attractive forces corresponding to the slope of the graph in between the two lowest separation distances, which clearly indicates a preferable C agglomeration but the absence of C clustering.
250
251 In IBS, configurations may arise, in which the impinging C atom creates a vacant site (V) near a C$_{\text{i}}$ DB, but does not occupy it.
252 All these structures were found to be energetically preferred compared to isolated largely separated defects \cite{zirkelbach11} showing an entirely attractive interaction between defects of these types.
253 The ground-state configuration is obtained for a V located right next to the C atom of the DB.
254 The C atom moves towards the vacant site forming a stable C$_{\text{s}}$ configuration resulting in the release of a huge amount of energy.
255 The second most favorable configuration is accomplished for a V located right next to the Si atom of the DB structure.
256 This is due to the reduction of compressive strain of the Si DB atom and its two upper Si neighbors present in the isolated C$_{\text{i}}$ DB configuration.
257 This configuration is followed by the structure, in which the V is created at one of the neighbored lattice sites below one of the Si atoms that are bound to the C atom of the initial DB.
258 Relaxed structures of the latter two defect combinations are shown in the bottom left of Fig.~\ref{fig:314-539} and \ref{fig:059-539} respectively together with their energetics during transition into the ground state.
259 \begin{figure}
260 \subfloat[V created right next to the Si atom of the initial DB. Activation energy: {\unit[0.1]{eV}}.]{%
261 \includegraphics[width=\columnwidth]{314-539.ps}
262 \label{fig:314-539}
263 }\\
264 %\caption{Migration barrier and structures of the transition of the initial C$_{\text{i}}$ \hkl[0 0 -1] DB and a V created right next to the Si atom of the initial DB (left) into a C$_{\text{s}}$ configuration (right). An activation energy of \unit[0.1]{eV} is observed.}
265 \subfloat[V created next to one of the Si atoms that is bound to the C atom of the initial DB. Activation energy: {\unit[0.6]{eV}}.]{%
266 \includegraphics[width=\columnwidth]{059-539.ps}
267 \label{fig:059-539}
268 }
269 %\caption{Migration barrier and structures of the transition of the initial C$_{\text{i}}$ \hkl[0 0 -1] DB and a V created next to one of the Si atoms that is bound to the C atom of the initial DB (left) into a C$_{\text{s}}$ configuration (right). An activation energy of \unit[0.6]{eV} is observed.}
270 \caption{Migration barrier and structures of transitions of an initial C$_{\text{i}}$ \hkl[0 0 -1] DB and a V (left) into a C$_{\text{s}}$ configuration (right).}
271 \end{figure}
272 These transitions exhibit activation energies as low as \unit[0.1]{eV} and \unit[0.6]{eV}.
273 In the first case the Si and C atom of the DB move towards the vacant and initial DB lattice site respectively.
274 In the second case Si number 1, which is substituted by the C$_{\text{i}}$ atom, migrates towards the vacant site.
275 In both cases, the formation of additional bonds is responsible for the vast gain in energy rendering almost impossible the reverse processes.
276 Considering the small activation energies, a high probability for the formation of stable C$_{\text{s}}$ must be concluded.
277
278 In addition, it is instructive to investigate combinations of C$_{\text{s}}$ and Si$_{\text{i}}$, which can be created in IBS by highly energetic C atoms that kick out a Si atom from its lattice site, resulting in a Si self-interstitial accompanied by a vacant site, which might get occupied by another C atom that lost almost all of its kinetic energy.
279 The most favorable configuration, which is C$_{\text{s}}$ located right next to the ground-state Si$_{\text{i}}$ defect, i.e.\  the Si$_{\text{i}}$ \hkl<1 1 0> DB, and the transition of this structure into the ground-state configuration, i.e. the C$_{\text{i}}$ \hkl<1 0 0> DB is displayed in Fig.~\ref{fig:162-097}
280 \begin{figure}
281 \includegraphics[width=\columnwidth]{162-097.ps}
282 \caption{Transition of a \hkl[1 1 0] Si$_{\text{i}}$ DB next to C$_{\text{s}}$ (right) into the C$_{\text{i}}$ \hkl[0 0 -1] DB configuration (left).}
283 \label{fig:162-097}
284 \end{figure}
285 Due to the low barrier of \unit[0.12]{eV}, the C$_{\text{i}}$ \hkl<1 0 0> DB configuration is very likely to occur.
286 However, the barrier of only \unit[0.77]{eV} for the reverse process indicates the possibility to form a C$_{\text{s}}$ and Si$_{\text{i}}$ DB out of the ground state activated without much effort either thermally or by introduced energy of the implantation process.
287 \begin{figure}
288 \includegraphics[width=\columnwidth]{c_sub_si110.ps}
289 %\includegraphics[width=\columnwidth]{c_sub_si110_data.ps}
290 %\caption{Binding energies of combinations of a C$_{\text{s}}$ and a Si$_{\text{i}}$ DB with respect to the separation distance.}
291 \caption{Binding energies of combinations of a C$_{\text{s}}$ and a Si$_{\text{i}}$ DB with respect to the separation distance. The interaction strength of the defect pairs are well approximated by a Lennard-Jones 6-12 potential, which is used for curve fitting.}
292 \label{fig:dc_si-s}
293 \end{figure}
294 Furthermore, the interaction strength quickly drops to zero with increasing separation distance as can be seen in Fig.~\ref{fig:dc_si-s}.
295 The interaction of the defects is well approximated by a Lennard-Jones (LJ) 6-12 potential.
296 Unable to model possible positive values of the binding energy, i.e. unfavorable configurations, the LJ fit should rather be thought as a guide for the eye describing the decrease of the interaction strength, i.e. the absolute value of the binding energy, with increasing separation distance.
297 The LJ fit estimates almost zero interaction already at \unit[0.5-0.6]{nm}, indicating a low interaction capture radius of the defect pair.
298 In IBS separations exceeding this capture radius are easily produced.
299 For these reasons, it must be concluded that configurations of C$_{\text{s}}$ and Si$_{\text{i}}$ instead of the thermodynamically stable C$_{\text{i}}$ \hkl<1 0 0> DB play a decisive role in IBS, a process far from equilibrium.
300 Indeed, in a previous study, an {\em ab initio} molecular dynamics run at \unit[900]{$^{\circ}$C} results in a separation of the C$_{\text{s}}$ and Si$_{\text{i}}$ DB located right next to each other \cite{zirkelbach11}.
301
302 To summarize, these obtained  results suggest an increased participation of C$_{\text{s}}$ already in the initial stages of precipitation under IBS conditions.
303
304 \section{Large scale empirical potential MD results}
305
306 Results of the MD simulations at \unit[450]{$^{\circ}$C}, an operative and efficient temperature in IBS \cite{lindner01}, indicate the formation of C$_{\text{i}}$ \hkl<1 0 0> DBs if C is inserted into the total simulation volume.
307 However, no agglomeration is observed within the simulated time, which was increased up to several nanoseconds.
308 This is attributed to the drastically overestimated migration barrier of the C defect, which hampers C agglomeration.
309 To overcome this obstacle, the simulation temperature is successively increased up to \unit[2050]{$^{\circ}$C}.
310 Fig.~\ref{fig:tot} shows the resulting radial distribution functions of Si-C bonds for various elevated temperatures.
311 \begin{figure}
312 \includegraphics[width=\columnwidth]{tot_pc_thesis.ps}
313 \caption{Radial distribution function for Si-C pairs for C insertion at various elevated temperatures. Si-C distances of a single C$_{\text{s}}$ defect configuration are plotted.}
314 \label{fig:tot}
315 \end{figure}
316 Although not intended, a transformation from a structure dominated by C$_{\text{i}}$ into a structure consisting of C$_{\text{s}}$ with increasing temperature can clearly be observed if compared with the radial distribution of C$_{\text{s}}$ in c-Si.
317
318 Thus, the C$_{\text{s}}$ defect and resulting stretched coherent structures of SiC must be considered to play an important role in the IBS at elevated temperatures.
319 This, in fact, satisfies experimental findings of annealing experiments \cite{strane94,nejim95,serre95} as well as previous DFT results, which suggest C$_{\text{s}}$ to be involved at higher temperatures and in conditions that deviate the system out of the thermodynamic ground state.
320
321 \section{Summary and discussion}
322
323 Although investigations of defect combinations show the agglomeration of C$_{\text{i}}$ DBs to be energetically most favorable, configurations that may arise during IBS were presented, their dynamics indicating C$_{\text{s}}$ to play an important role particularly at high temperatures.
324 This is supported by the classical MD results, which show an increased participation of C$_{\text{s}}$ at increased temperatures that allow the system to deviate from the ground state.
325
326 Based on these findings, it is concluded that in IBS at elevated temperatures, SiC conversion takes place by an initial agglomeration of C$_{\text{s}}$ into coherent, tensily strained structures of SiC followed by precipitation into incoherent SiC once a critical size is reached and the increasing strain energy of the coherent structure surpasses the interfacial energy of the incoherent precipitate.
327 Rearrangement of stable C$_{\text{s}}$ is enabled by excess Si$_{\text{i}}$, which not only acts as a vehicle for C but also as a supply of Si atoms needed elsewhere to form the SiC structure and to reduce possible strain at the interface of coherent SiC precipitates and the Si host.
328
329 %It is worth to point out that the experimentally observed alignment of the \hkl(h k l) planes of precipitate and substrate is satisfied by this mechanism.
330 %In contrast, the topotactic orientation of the SiC precipitate originating from an agglomerate consisting exclusively of C-Si dimers would necessarily involve a much more profound change in structure.
331
332 \begin{acknowledgement}
333 We gratefully acknowledge financial support by the Bayerische Forschungsstiftung (Grant No. DPA-61/05) and the Deutsche Forschungsgemeinschaft (Grant No. DFG SCHM 1361/11).
334 \end{acknowledgement}
335
336 %\bibliography{../../bibdb/bibdb}{}
337 %\bibliographystyle{pss.bst}
338
339
340 \providecommand{\WileyBibTextsc}{}
341 \let\textsc\WileyBibTextsc
342 \providecommand{\othercit}{}
343 \providecommand{\jr}[1]{#1}
344 \providecommand{\etal}{~et~al.}
345
346
347 \begin{thebibliography}{[10]}
348
349 \bibitem{edgar92}% article
350  \textsc{J.\,H. Edgar}\iffalse Prospects for device implementation of wide band
351   gap semiconductors\fi,
352  \jr{J. Mater. Res.} \textbf{7}(January), 235 (1992).
353
354
355 \bibitem{capano97}% article
356  \textsc{M.\,A. Capano} and  \textsc{R.\,J. Trew}\iffalse Silicon carbide
357   electronic materials and devices\fi,
358  \jr{MRS Bull.} \textbf{22}(3), 19--22 (1997).
359
360
361 \othercit
362 \bibitem{park98}% book
363  \textsc{Y.\,S. Park},
364 Si{C} Materials and Devices (Academic Press, San Diego, 1998).
365
366
367 \bibitem{borders71}% article
368  \textsc{J.\,A. Borders},  \textsc{S.\,T. Picraux},  and
369   \textsc{W.~Beezhold}\iffalse {FORMATION} {OF} si{C} {IN} {SILICON} {BY} {ION}
370   {IMPLANTATION}\fi,
371  \jr{Appl. Phys. Lett.} \textbf{18}(11), 509--511 (1971).
372
373
374 \bibitem{lindner01}% article
375  \textsc{J.\,K.\,N. Lindner}\iffalse Ion beam synthesis of buried si{C} layers
376   in silicon: Basic physical processes\fi,
377  \jr{Nucl. Instrum. Methods Phys. Res. B} \textbf{178}(1-4), 44--54 (2001).
378
379
380 \othercit
381 \bibitem{werner96}% inproceedings
382  \textsc{P.~Werner},  \textsc{R.~K{\"{o}}gler},  \textsc{W.~Skorupa},  and
383   \textsc{D.~Eichler},
384 {TEM} investigation of {C}-si defects in carbon implanted silicon,
385  in: Proceedings of the 11th International Conference on Ion Implantation
386   Technology.,  (June 1996),  pp.\,675--678.
387
388
389 \bibitem{lindner99_2}% article
390  \textsc{J.\,K.\,N. Lindner} and  \textsc{B.~Stritzker}\iffalse Mechanisms in
391   the ion beam synthesis of si{C} layers in silicon\fi,
392  \jr{Nucl. Instrum. Methods Phys. Res. B} \textbf{148}(1-4), 528--533 (1999).
393
394
395 \bibitem{koegler03}% article
396  \textsc{R.~K{\"{o}}gler},  \textsc{F.~Eichhorn},  \textsc{J.\,R. Kaschny},
397   \textsc{A.~M{\"{u}}cklich},  \textsc{H.~Reuther},  \textsc{W.~Skorupa},
398   \textsc{C.~Serre},  and  \textsc{A.~Perez-Rodriguez}\iffalse Synthesis of
399   nano-sized si{C} precipitates in si by simultaneous dual-beam implantation of
400   {C}+ and si+ ions\fi,
401  \jr{Appl. Phys. A} \textbf{76}(March), 827--835 (2003).
402
403
404 \bibitem{strane94}% article
405  \textsc{J.\,W. Strane},  \textsc{H.\,J. Stein},  \textsc{S.\,R. Lee},
406   \textsc{S.\,T. Picraux},  \textsc{J.\,K. Watanabe},  and  \textsc{J.\,W.
407   Mayer}\iffalse Precipitation and relaxation in strained si[sub 1 - y]{C}[sub
408   y]/si heterostructures\fi,
409  \jr{J. Appl. Phys.} \textbf{76}(6), 3656--3668 (1994).
410
411
412 \bibitem{nejim95}% article
413  \textsc{A.~Nejim},  \textsc{P.\,L.\,F. Hemment},  and
414   \textsc{J.~Stoemenos}\iffalse Si{C} buried layer formation by ion beam
415   synthesis at 950 [degree]{C}\fi,
416  \jr{Appl. Phys. Lett.} \textbf{66}(20), 2646--2648 (1995).
417
418
419 \bibitem{serre95}% article
420  \textsc{C.~Serre},  \textsc{A.~P\'{e}rez-Rodr\'{\i}guez},
421   \textsc{A.~Romano-Rodr\'{\i}guez},  \textsc{J.\,R. Morante},
422   \textsc{R.~K{\"{o}}gler},  and  \textsc{W.~Skorupa}\iffalse Spectroscopic
423   characterization of phases formed by high-dose carbon ion implantation in
424   silicon\fi,
425  \jr{J. Appl. Phys.} \textbf{77}(7), 2978--2984 (1995).
426
427
428 \bibitem{bar-yam84}% article
429  \textsc{Y.~Bar-Yam} and  \textsc{J.\,D. Joannopoulos}\iffalse Barrier to
430   migration of the silicon self-interstitial\fi,
431  \jr{Phys. Rev. Lett.} \textbf{52}(13), 1129--1132 (1984).
432
433
434 \bibitem{car84}% article
435  \textsc{R.~Car},  \textsc{P.\,J. Kelly},  \textsc{A.~Oshiyama},  and
436   \textsc{S.\,T. Pantelides}\iffalse Microscopic theory of atomic diffusion
437   mechanisms in silicon\fi,
438  \jr{Phys. Rev. Lett.} \textbf{52}(20), 1814--1817 (1984).
439
440
441 \bibitem{bloechl93}% article
442  \textsc{P.\,E. Bl{\"o}chl},  \textsc{E.~Smargiassi},  \textsc{R.~Car},
443   \textsc{D.\,B. Laks},  \textsc{W.~Andreoni},  and  \textsc{S.\,T.
444   Pantelides}\iffalse First-principles calculations of self-diffusion constants
445   in silicon\fi,
446  \jr{Phys. Rev. Lett.} \textbf{70}(16), 2435--2438 (1993).
447
448
449 \bibitem{tang97}% article
450  \textsc{M.~Tang},  \textsc{L.~Colombo},  \textsc{J.~Zhu},  and  \textsc{T.\,D.
451   de\,la Rubia}\iffalse Intrinsic point defects in crystalline silicon:
452   Tight-binding molecular dynamics studies of self-diffusion,
453   interstitial-vacancy recombination, and formation volumes\fi,
454  \jr{Phys. Rev. B} \textbf{55}(21), 14279--14289 (1997).
455
456
457 \bibitem{leung99}% article
458  \textsc{W.\,K. Leung},  \textsc{R.\,J. Needs},  \textsc{G.~Rajagopal},
459   \textsc{S.~Itoh},  and  \textsc{S.~Ihara}\iffalse Calculations of silicon
460   self-interstitial defects\fi,
461  \jr{Phys. Rev. Lett.} \textbf{83}(12), 2351--2354 (1999).
462
463
464 \bibitem{al-mushadani03}% article
465  \textsc{O.\,K. Al-Mushadani} and  \textsc{R.\,J. Needs}\iffalse Free-energy
466   calculations of intrinsic point defects in silicon\fi,
467  \jr{Phys. Rev. B} \textbf{68}(23), 235205 (2003).
468
469
470 \bibitem{hobler05}% article
471  \textsc{G.~Hobler} and  \textsc{G.~Kresse}\iffalse Ab initio calculations of
472   the interaction between native point defects in silicon\fi,
473  \jr{Mater. Sci. Eng., B} \textbf{124-125}, 368--371 (2005),
474 EMRS 2005, Symposium D - Materials Science and Device Issues for Future
475   Technologies.
476
477
478 \bibitem{sahli05}% article
479  \textsc{B.~Sahli} and  \textsc{W.~Fichtner}\iffalse Ab initio molecular
480   dynamics simulation of self-interstitial diffusion in silicon\fi,
481  \jr{Phys. Rev. B} \textbf{72}(24), 245210 (2005).
482
483
484 \bibitem{posselt08}% article
485  \textsc{M.~Posselt},  \textsc{F.~Gao},  and  \textsc{H.~Bracht}\iffalse
486   Correlation between self-diffusion in si and the migration mechanisms of
487   vacancies and self-interstitials: An atomistic study\fi,
488  \jr{Phys. Rev. B} \textbf{78}(3), 035208 (2008).
489
490
491 \bibitem{ma10}% article
492  \textsc{S.~Ma} and  \textsc{S.~Wang}\iffalse Ab initio study of self-diffusion
493   in silicon over a wide temperature range: Point defect states and migration
494   mechanisms\fi,
495  \jr{Phys. Rev. B} \textbf{81}(19), 193203 (2010).
496
497
498 \bibitem{tersoff90}% article
499  \textsc{J.~Tersoff}\iffalse Carbon defects and defect reactions in silicon\fi,
500  \jr{Phys. Rev. Lett.} \textbf{64}(15), 1757--1760 (1990).
501
502
503 \bibitem{dal_pino93}% article
504  \textsc{A.~{Dal Pino}},  \textsc{A.\,M. Rappe},  and  \textsc{J.\,D.
505   Joannopoulos}\iffalse Ab initio investigation of carbon-related defects in
506   silicon\fi,
507  \jr{Phys. Rev. B} \textbf{47}(19), 12554--12557 (1993).
508
509
510 \bibitem{capaz94}% article
511  \textsc{R.\,B. Capaz},  \textsc{A.~{Dal Pino}},  and  \textsc{J.\,D.
512   Joannopoulos}\iffalse Identification of the migration path of interstitial
513   carbon in silicon\fi,
514  \jr{Phys. Rev. B} \textbf{50}(11), 7439--7442 (1994).
515
516
517 \bibitem{burnard93}% article
518  \textsc{M.\,J. Burnard} and  \textsc{G.\,G. DeLeo}\iffalse Interstitial carbon
519   and the carbon-carbon pair in silicon: Semiempirical electronic-structure
520   calculations\fi,
521  \jr{Phys. Rev. B} \textbf{47}(16), 10217--10225 (1993).
522
523
524 \bibitem{leary97}% article
525  \textsc{P.~Leary},  \textsc{R.~Jones},  \textsc{S.~{\"O}berg},  and
526   \textsc{V.\,J.\,B. Torres}\iffalse Dynamic properties of interstitial carbon
527   and carbon-carbon pair defects in silicon\fi,
528  \jr{Phys. Rev. B} \textbf{55}(4), 2188--2194 (1997).
529
530
531 \bibitem{capaz98}% article
532  \textsc{R.\,B. Capaz},  \textsc{A.~{Dal Pino}},  and  \textsc{J.\,D.
533   Joannopoulos}\iffalse Theory of carbon-carbon pairs in silicon\fi,
534  \jr{Phys. Rev. B} \textbf{58}(15), 9845--9850 (1998).
535
536
537 \bibitem{mattoni2002}% article
538  \textsc{A.~Mattoni},  \textsc{F.~Bernardini},  and
539   \textsc{L.~Colombo}\iffalse Self-interstitial trapping by carbon complexes in
540   crystalline silicon\fi,
541  \jr{Phys. Rev. B} \textbf{66}(19), 195214 (2002).
542
543
544 \bibitem{kresse96}% article
545  \textsc{G.~Kresse} and  \textsc{J.~Furthm{\"{u}}ller}\iffalse Efficiency of
546   ab-initio total energy calculations for metals and semiconductors using a
547   plane-wave basis set\fi,
548  \jr{Comput. Mater. Sci.} \textbf{6}(1), 15--50 (1996).
549
550
551 \bibitem{perdew86}% article
552  \textsc{J.\,P. Perdew} and  \textsc{Y.~Wang}\iffalse Accurate and simple
553   density functional for the electronic exchange energy: Generalized gradient
554   approximation\fi,
555  \jr{Phys. Rev. B} \textbf{33}(12), 8800--8802 (1986).
556
557
558 \bibitem{perdew92}% article
559  \textsc{J.\,P. Perdew},  \textsc{J.\,A. Chevary},  \textsc{S.\,H. Vosko},
560   \textsc{K.\,A. Jackson},  \textsc{M.\,R. Pederson},  \textsc{D.\,J. Singh},
561   and  \textsc{C.~Fiolhais}\iffalse Atoms, molecules, solids, and surfaces:
562   Applications of the generalized gradient approximation for exchange and
563   correlation\fi,
564  \jr{Phys. Rev. B} \textbf{46}(11), 6671--6687 (1992).
565
566
567 \bibitem{hamann79}% article
568  \textsc{D.\,R. Hamann},  \textsc{M.~Schl{\"u}ter},  and
569   \textsc{C.~Chiang}\iffalse Norm-conserving pseudopotentials\fi,
570  \jr{Phys. Rev. Lett.} \textbf{43}(20), 1494--1497 (1979).
571
572
573 \bibitem{vanderbilt90}% article
574  \textsc{D.~Vanderbilt}\iffalse Soft self-consistent pseudopotentials in a
575   generalized eigenvalue formalism\fi,
576  \jr{Phys. Rev. B} \textbf{41}(11), 7892--7895 (1990).
577
578
579 \bibitem{kaukonen98}% article
580  \textsc{M.~Kaukonen},  \textsc{P.\,K. Sitch},  \textsc{G.~Jungnickel},
581   \textsc{R.\,M. Nieminen},  \textsc{S.~P{\"o}ykk{\"o}},  \textsc{D.~Porezag},
582   and  \textsc{T.~Frauenheim}\iffalse Effect of {N} and {B} doping on the
583   growth of {CVD} diamond $(100):{H}(2\ifmmode\times\else\texttimes\fi{}1)$
584   surfaces\fi,
585  \jr{Phys. Rev. B} \textbf{57}(16), 9965--9970 (1998).
586
587
588 \bibitem{albe_sic_pot}% article
589  \textsc{P.~Erhart} and  \textsc{K.~Albe}\iffalse Analytical potential for
590   atomistic simulations of silicon, carbon, and silicon carbide\fi,
591  \jr{Phys. Rev. B} \textbf{71}(3), 035211 (2005).
592
593
594 \bibitem{berendsen84}% article
595  \textsc{H.\,J.\,C. Berendsen},  \textsc{J.\,P.\,M. Postma},  \textsc{W.\,F.
596   van Gunsteren},  \textsc{A.~DiNola},  and  \textsc{J.\,R. Haak}\iffalse
597   Molecular dynamics with coupling to an external bath\fi,
598  \jr{J. Chem. Phys.} \textbf{81}(8), 3684--3690 (1984).
599
600
601 \bibitem{verlet67}% article
602  \textsc{L.~Verlet}\iffalse Computer {"}experiments{"} on classical fluids.
603   {I}. thermodynamical properties of lennard-jones molecules\fi,
604  \jr{Phys. Rev.} \textbf{159}(1), 98 (1967).
605
606
607 \bibitem{watkins76}% article
608  \textsc{G.\,D. Watkins} and  \textsc{K.\,L. Brower}\iffalse {EPR} observation
609   of the isolated interstitial carbon atom in silicon\fi,
610  \jr{Phys. Rev. Lett.} \textbf{36}(22), 1329--1332 (1976).
611
612
613 \bibitem{song90}% article
614  \textsc{L.\,W. Song} and  \textsc{G.\,D. Watkins}\iffalse {EPR} identification
615   of the single-acceptor state of interstitial carbon in silicon\fi,
616  \jr{Phys. Rev. B} \textbf{42}(9), 5759--5764 (1990).
617
618
619 \bibitem{lindner06}% article
620  \textsc{J.\,K.\,N. Lindner},  \textsc{M.~H{\"a}berlen},
621   \textsc{G.~Thorwarth},  and  \textsc{B.~Stritzker}\iffalse On the balance
622   between ion beam induced nanoparticle formation and displacive precipitate
623   resolution in the {C}-si system\fi,
624  \jr{Mater. Sci. Eng., C} \textbf{26}(5-7), 857--861 (2006),
625 Current Trends in Nanoscience - from Materials to Applications.
626
627
628 \bibitem{tipping87}% article
629  \textsc{A.\,K. Tipping} and  \textsc{R.\,C. Newman}\iffalse The diffusion
630   coefficient of interstitial carbon in silicon\fi,
631  \jr{Semicond. Sci. Technol.} \textbf{2}(5), 315--317 (1987).
632
633
634 \bibitem{zirkelbach11}% article
635  \textsc{F.~Zirkelbach},  \textsc{B.~Stritzker},  \textsc{K.~Nordlund},
636   \textsc{J.\,K.\,N. Lindner},  \textsc{W.\,G. Schmidt},  and
637   \textsc{E.~Rauls}\iffalse Combined \textit{ab initio} and classical potential
638   simulation study on silicon carbide precipitation in silicon\fi,
639  \jr{Phys. Rev. B} \textbf{84}(August), 064126 (2011).
640
641
642 \bibitem{zirkelbach10}% article
643  \textsc{F.~Zirkelbach},  \textsc{B.~Stritzker},  \textsc{K.~Nordlund},
644   \textsc{J.\,K.\,N. Lindner},  \textsc{W.\,G. Schmidt},  and
645   \textsc{E.~Rauls}\iffalse Defects in carbon implanted silicon calculated by
646   classical potentials and first-principles methods\fi,
647  \jr{Phys. Rev. B} \textbf{82}(9), 094110 (2010).
648
649
650 \end{thebibliography}
651
652
653 \end{document}