--- /dev/null
+\pdfoutput=0
+%\documentclass[doublespacing]{elsart}
+\documentclass[]{elsart}
+\usepackage[activate]{pdfcprot}
+\usepackage{verbatim}
+\usepackage{a4}
+\usepackage{a4wide}
+\usepackage[german,english]{babel}
+\usepackage[latin1]{inputenc}
+\usepackage[T1]{fontenc}
+\usepackage{amsmath}
+\usepackage{latexsym}
+\usepackage{ae}
+\usepackage{aecompl}
+\usepackage[dvips]{graphicx}
+\graphicspath{{../img/}}
+\usepackage{color}
+\usepackage{pstricks}
+\usepackage{pst-node}
+\usepackage{rotating}
+
+%\bibliographystyle{h-physrev3.bst}
+
+\selectlanguage{english}
+
+\begin{document}
+
+\begin{frontmatter}
+
+\title{Molecular dynamics simulation of defect formation and precipitation
+ in heavily carbon doped silicon}
+
+\author[augsburg]{F. Zirkelbach\corauthref{cor}}
+\author[augsburg]{J. K. N. Lindner}
+\author[helsinki]{K. Nordlund}
+\author[augsburg]{B. Stritzker}
+
+\address[augsburg]{Experimentalphysik IV, Institut f\"ur Physik,
+ Universit\"at Augsburg,\\
+ Universit\"atsstr. 1, D-86135 Augsburg, Germany}
+\address[helsinki]{Accelerator Laboratory, Department of Physical Sciences,
+ University of Helsinki,\\
+ Pietari Kalmink. 2, 00014 Helsinki, Finland}
+
+\corauth[cor]{Corresponding author.\\
+ Tel.: +49-821-5983008; fax: +49-821-5983425.\\
+ E-mail address: frank.zirkelbach@physik.uni-augsburg.de
+ (Frank Zirkelbach)}
+
+\begin{abstract}
+The precipitation process of silicon carbide in heavily carbon doped silicon is not yet fully understood.
+High resolution transmission electron microscopy observations suggest that in a first step carbon atoms form C-Si dumbbells on regular Si lattice sites which agglomerate into large clusters.
+In a second step, when the cluster size reaches a radius of a few nm, the high interfacial energy due to the SiC/Si lattice misfit of almost 20\% is overcome and the precipitation occurs.
+By simulation, details of the precipitation process can be obtained on the atomic level.
+A newly parametrized Tersoff-like bond order potential is used to model the system appropriately.
+First results gained by molecular dynamics simulations using this potential are presented.
+\end{abstract}
+
+\begin{keyword}
+Silicon, carbon, silicon carbide, nucleation, defect formation,
+molecular dynamics simulations
+\end{keyword}
+
+\end{frontmatter}
+
+\section{Introduction}
+Understanding the precipitation process of cubic silicon carbide (3C-SiC) in heavily carbon doped silicon will enable significant technological progress in thin film formation of the important wide band gap semiconductor material SiC \cite{1}.
+On the other hand it will likewise offer perspectives for processes which rely upon prevention of precipitation events, e.g. the fabrication of strained, pseudomorphic $\text{Si}_{1-y}\text{C}_y$ heterostructures \cite{2}.
+
+Epitaxial growth of 3C-SiC films is achieved either by ion beam synthesis (IBS) \cite{3}, chemical vapour deposition (CVD) or molecular beam epitaxy (MBE) techniques.
+While in CVD and MBE surface effects need to be taken into account, SiC formation during IBS takes place in the bulk of the Si crystal.
+In the present work the simulation assumes conditions similar to those of the ion implantation process.
+
+First of all a picture of the supposed precipitation event is presented.
+Then the applied simulation sequences are introduced.
+Finally first simulation results are presented.
+
+\section{Supposed Si to 3C-SiC conversion mechanism}
+Crystalline silicon (c-Si) has diamond structure and thus is composed of two fcc lattices which are displaced by one quarter of the volume diagonal.
+3C-SiC grows in zincblende structure, i.e. is also composed of two fcc lattices one out of which is occupied by Si and the other by C atoms.
+The length of four lattice constants of c-Si is approximately equal to the length of five 3C-SiC lattice constants ($4a_{\text{Si}}\approx 5a_{\text{3C-SiC}}$) resulting in a lattice misfit of almost 20\%.
+Due to this the silicon atomic density of 3C-SiC is slightly lower than the one of c-Si.
+
+There is a supposed conversion mechanism of heavily carbon doped c-Si into 3C-SiC \cite{4}.
+As concluded from high resolution transmission electron microscopy (HREM) carbon atoms introduced into c-Si form C-Si dumbbells on regular Si lattice sites.
+The dumbbells agglomerate into large clusters, called embryos.
+Finally, when the cluster size reaches a critical radius of 2 to 4 nm, the high interfacial energy due to the 3C-SiC/c-Si lattice misfit is overcome and precipitation occurs.
+In HREM this transformation manifestsitself by the disappearance of patches with dark background in favour of high resolution Moir\'e fringes due to 3C-SiC precipitates embedded in c-Si \cite{3,4}.
+Due to the slightly lower silicon density of 3C-SiC excessive silicon atoms exist which will most probably end up as self-interstitials in the c-Si matrix since there is more space than in 3C-SiC.
+
+Thus, in addition to the precipitation event itself, knowledge of C and Si interstitials in Si are of great interest in order to investigate the precipitation of 3C-SiC in heavily C doped c-Si.
+
+\section{Simulation sequences}
+A molecular dynamics simulation approach is used to examine the steps involved in the precipitation process.
+For integrating the equations of motion the velocity verlet algorithm \cite{5} with a timestep of 1 fs is adopted.
+The interaction of the silicon and carbon atoms is realized by a newly parametrized Tersoff-like bond order potential \cite{6}.
+Since temperature and pressure of the system is kept constant in experiment the isothermal-isobaric NPT ensemble is chosen for the simulation.
+Coupling to the heat bath is achieved by the Berendsen thermostat \cite{7} with a time constant of 100 fs.
+The pressure is scaled by the Berendsen barostat \cite{7} again using a timeconstant of 100 fs and a bulk modulus of 100 GPa for silicon.
+To exclude surface effects periodic boundary conditions are applied.
+In order to study the behaviour of carbon in c-Si, two different simulation sequences are used.
+
+To investigate the interstitial configurations of C and Si atoms in c-Si, a simulation volume of $9\times9\times9$ silicon unit cells is used.
+The temperature is set to 0 K.
+In separated simulation runs a carbon or a silicon atom is inserted in the central unit cell at the tetrahedral (0,0,0) (${\color{red}\triangleleft}$), a hexagonal (-1/8,-1/8,1/8) (${\color{green}\triangleright}$), a nearly dumbbell (-1/8,-1/8,-1/4) (${\color{magenta}\Box}$) or at a random position (in units of the silicon lattice constant) where the origin is located in the centre of the unit cell (Fig. 1).
+In order to avoid too high potential energies in the case of the dumbbell configuration the nearest silicon neighbour atom is manually shifted to (-3/8,-3/8,-1/4) ($\circ$).
+The energy introduced into the system by the insertion of atoms is scaled out within a relaxation phase of 2 ps.
+
+For the simulations aiming to reproduce a precipitation process the simulation volume is 31 silicon lattice constants in each direction.
+The system temperature is set to $450\, ^{\circ} \textrm{C}$ like in IBS \cite{3}.
+$6000$ carbon atoms (the number necessary to form a 3C-SiC precipitate with a radius of $\sim 3$ nm) are consecutively inserted in a way to keep the system temperature constant.
+Carbon is inserted statistically distributed over either the whole simulation volme (insertion volume $V_1$), an insertion volume $V_2$ corresponding to the minimum size of a 3C-SiC precipitate or the volume $V_3$ containing the number of Si atoms necessary for the formation of such a minimum precipitate.
+The two latter insertion volumes are considered since no diffusion of carbon atoms is observed at this temperature.
+Following the insertion procedure the system is cooled down to $20\, ^{\circ} \textrm{C}$.
+
+\section{Results}
+The insertion of individual Si atoms leads to stable tetrahedral and <110> dumbbell self-interstitial configurations as observed earlier in \cite{6}.
+The formation energies are 3.4 eV and 4.4 eV, respectively.
+However, the hexagonal dumbbell configuration is not stable opposed to results in \cite{6}.
+The Si interstitial atom moves towards an energetically more favorable position very close to the tetrahedral one but slightly displaced along the three coordinate axes.
+The formation energy of 4.0 eV for this type of interstitials equals the result obtained in \cite{6} for the hexagonal one.
+The same type of interstitial arises using random insertions.
+In addition, variations exist in which the displacement is only along two <100> axes ($E_f=3.8\, \text{eV}$) or along a single <100> axis ($E_f=3.6\, \text{eV}$) succesively approximating the tetrahedral configuration and formation energy.
+
+In carbon insertion simulations the tetrahedral and <110> dumbbell carbon interstitial configurations are observed to be stable.
+The formation energies are 2.7 eV and 1.8 eV, respectively.
+Again, the hexagonal configuration is found to be unstable.
+The carbon interstitial atom moves to the more favorable <100> dumbbell position which has a formation energy of 0.5 eV.
+This interstitial configuration is shown in Fig. 2.
+There is experimental evidence \cite{8} of the existence of this configuration.
+It is frequently generated and has the lowest formation energy of all the defects observed in all the simulation runs in which carbon is inserted at random positions.
+Fig. 3 schematically displays the <100> dumbbell configuration including the atomic displacements relative to the atoms initial positions together with resulting new Si-Si and C-Si pair distances.
+
+Fig. 4 shows resulting pair correlation functions of the simulation runs targeting the observation of precipitation events.
+The contributions of Si-C and C-C pairs are presented separately.
+The three curves represent results for the insertion of 6000 C atoms into the three different volumes $V_1$, $V_2$ and $V_3$, as explained above.
+There is no significant difference between carbon insertions into $V_2$ and $V_3$.
+
+The number of C-C bonds for $V_1$ is much smaller than for $V_2$ and $V_3$ since carbon atoms are spread over the total simulation volume, which means that there is only 0.2 carbon atoms per silicon unit cell on average.
+The first C-C peak appears at about 0.15 nm.
+This is comparable to the nearest neighbour distance of graphite or diamond.
+It is assumed that these carbon atoms form strong C-C bonds which is supported by a decrease of the total energy during carbon insertion for the $V_2$ and $V_3$ in contrast to the $V_1$ simulation.
+
+The C-C peak at $\sim 0.31$ nm perfectly matches the distance of two carbon atoms in the 3C-SiC lattice. In 3C-SiC the same distance is also expected for Si-Si bonds.
+The Si-Si portion of the pair correlation function for the case of 3000 C atoms inserted into $V_2$ is displayed in Fig. 5, together with a reference function for pure c-Si.
+After insertion of C atoms the Si-Si pair correlation function in fact shows non-zero g(r) values around distance values of 0.31 nm while the amount of Si pairs at the regular distances of 0.24 and 0.38 nm decreases.
+However, no clear peak is observed but the interval of enhanced g(r) values corresponds to the width of the C-C g(r) peak.
+Analyses of randomly chosen configurations in which distances around 0.3 nm appear, identify <100> C-Si dumbbells to be responsible for stretching the Si-Si next neighbour distance for low C concentrations, i.e. for the $V_1$ and early stages of $V_2$ and $V_3$ simulation runs.
+This excellently agrees with the calculation for a single <100> dumbbell ($r(13)$ in Fig. 3).
+For higher C concentrations the defect concentration is likewise increased and a considerable amount of damage is introduced into the insertion volume.
+Damage and superposition of defects generate new displacement arrangements which become hard to categorize and trace and obviously lead to a broader distribution of pair distances.
+The step-like increase of Si-Si pairs at 0.29 nm is probably due to the Si-Si cut-off radius of 0.296 nm in the used bond order potential \cite{6}.
+The cut-off function causes artificial forces pushing the Si atoms out of the cut-off region.
+By visualizing the C-C atom pairs with distances of 0.31 nm, concatenated, differently oriented <100> dumbbell interstitials are frequently observed.
+This leads to the assumption that C atoms introduced in that way are the first elements arranged at distances expected for 3C-SiC which at a later stage pull the Si atoms into the right configuration.
+
+The Si-C pair correlation (Fig. 4) shows the formation of Si-C atom pairs displaced approximately by 0.187 nm.
+This corresponds quite well to the expected next neighbour distance of 0.189 nm for Si and C atoms in 3C-SiC.
+For low concentrations again <100> dumbbell configurations can be identified to contribute to this pair distance (r(3C) in Fig. 3).
+Further more, the bump at 0.175 (r(1C) in Fig. 3) and the peaks at about 0.333, 0.374, 0.386, 0.431 and 0.464 nm are comprehensible by comparing with calculated values in Fig. 3 and the dashed vertical lines in Fig. 4.
+The dashed lines mark further calculated C-Si atom pair distances appearing in the <100> C-Si dumbbell interstitial configuration which are not displayed in Fig. 3.
+The distinct peak at 0.26 nm, which exactly matches the cut-off radius of Si-C interaction, is again a potential artefact.
+For high carbon concentrations again superposition of defects and high amounts of damage characterize the pair distances.
+Short range order indeed is observed but only hardly visible is the long range order.
+This suggests to use higher temperatures particularly with regard to $V_2$ and $V_3$ simulations, i.e. simulations with high carbon concentrations, so that the system is able to get over an energy barrier escaping a local minimum configuration.
+An increase of the simulation temperature should not be in conflict with experimental conditions since the target temperature of $450\, ^{\circ} \textrm{C}$ measured during the IBS of 3C-SiC \cite{3} at the substrat holder is most certainly lower than the temperature predominating near the target surface, where 3C-SiC formation takes place.
+
+\section{Summary}
+The conversion of heavily C doped c-Si into 3C-SiC precipitates is investigated.
+Molecular dynamics simulation sequences to analyze interstitial configurations and the precipitation of 3C-SiC are introduced.
+The <100> C-Si dumbbel reproduced by simulation is observed the energetically most favorable configuration and explained in detail.
+First results of the simulations suggest that in the precipitation process C atoms introduced as differently oriented C-Si dumbbells into the c-Si matrix are the first elements properly arranged for the 3C-SiC formation.
+Furthermore, the observation of high amounts of damage particularly for high carbon concentrations demands for elevated system temperatures to achieve the precipitation event.
+
+%\bibliography{../../bibdb/bibdb}
+\begin{thebibliography}{8}
+ \bibitem{1} J. H. Edgar, J. Mater. Res. 7 (1992) 235.
+ \bibitem{2} J. W. Strane, S. R. Lee, H. J. Stein, S. T. Picraux,
+ J. K. Watanabe, J. W. Mayer, J. Appl. Phys. 79 (1996) 637.
+ \bibitem{3} J. K. N. Lindner, Appl. Phys. A 77 (2003) 27.
+ \bibitem{4} P. Werner, R. Koegler, W. Skorupa, and D. Eichler:
+ in Proc. 11th Int. Conf. on Ion Implantation Technology,
+ ed. by E. Ishidida, S. Banerjee, S. Mehta, T. C. Smith,
+ M. Current, L. Larson, A. Tasch, T. Romig
+ (Austin, TX, USA, June 16-21, 1996) pp. 675
+ \bibitem{5} L. Verlet, Phys. Rev. 159 (1967) 98.
+ \bibitem{6} P. Erhart and K. Albe, Phys. Rev. B 71 (2005) 35211.
+ \bibitem{7} H. J. C. Berendsen, J. P. M. Postma, W. F. van Gunsteren,
+ A. DiNola, J. R. Haak, J. Chem. Phys. 81 (1984) 3684.
+ \bibitem{8} G. D. Watkins and K. L. Brower, Phys. Rev. Lett. 36 (1976) 1329.
+ \end{thebibliography}
+
+\newpage
+\section*{Figures}
+
+\begin{figure}[!h]
+ \begin{center}
+ \includegraphics[width=8cm]{unit_cell_s.eps}
+ \caption{Insertion positions for the tetrahedral (${\color{red}\triangleleft}$), hexagonal (${\color{green}\triangleright}$) and <110> dumbbell (${\color{magenta}\Box}$) interstitial configuration.}
+ \end{center}
+\end{figure}
+
+\begin{figure}[!h]
+ \begin{center}
+ \includegraphics[width=6cm]{c_in_si_int_001db_0.eps}
+ \caption{Position of a <100> dumbbell carbon interstitial in silicon.
+ Only bonds of the carbon interstitial atom are shown.}
+ \end{center}
+\end{figure}
+
+\begin{figure}[!h]
+ \begin{center}
+ \includegraphics[width=15cm]{100-c-si-db_s.eps}
+ \caption{Schematic of the <100> C-Si dumbbell configuration.
+ Displacements of the atoms relative to their initial position are given.
+ The displacement of the carbon atom is relative to the initial position of atom 1.
+ New resulting Si-Si and C-Si pair distances r(ij) for the atoms (ij) shown in the schematic are recorded.}
+ \end{center}
+\end{figure}
+
+\begin{figure}[!h]
+ \begin{center}
+ \includegraphics[]{pc_si-c_c-c.eps}
+ \caption{Pair correlation functions for Si-C and C-C bonds.
+ The three curves represent results for the three different insertion volumes $V_1$, $V_2$ and $V_3$, as explained in the text.
+ The dashed vertical lines mark further calculated C-Si atom pair distances appearing in the <100> C-Si dumbbell interstitial configuration, which are not displayed in Fig. 3.}
+ \end{center}
+\end{figure}
+
+\begin{figure}[!h]
+ \begin{center}
+ \includegraphics[]{pc_si-si.eps}
+ \caption{Si-Si pair correlation function for pure Si and Si with 3000 inserted C atoms.
+ The inset shows a magnified region between 0.28 and 0.36 nm.}
+ \end{center}
+\end{figure}
+
+\end{document}
+++ /dev/null
-\pdfoutput=0
-\documentclass[a4paper,11pt]{article}
-\usepackage[activate]{pdfcprot}
-\usepackage{verbatim}
-\usepackage{a4}
-\usepackage{a4wide}
-\usepackage[german,english]{babel}
-\usepackage[latin1]{inputenc}
-\usepackage[T1]{fontenc}
-\usepackage{amsmath}
-\usepackage{latexsym}
-\usepackage{ae}
-\usepackage{aecompl}
-\usepackage[dvips]{graphicx}
-\graphicspath{{../img/}}
-\usepackage{color}
-\usepackage{pstricks}
-\usepackage{pst-node}
-\usepackage{rotating}
-
-\bibliographystyle{h-physrev3}
-
-\setlength{\headheight}{0mm} \setlength{\headsep}{0mm}
-\setlength{\topskip}{-10mm} \setlength{\textwidth}{17cm}
-\setlength{\oddsidemargin}{-10mm}
-\setlength{\evensidemargin}{-10mm} \setlength{\topmargin}{-1cm}
-\setlength{\textheight}{26cm} \setlength{\headsep}{0cm}
-
-%\linespread{2.0}
-
-\selectlanguage{english}
-
-\begin{document}
-
-% header
-\begin{center}
- {\LARGE {\bf Molecular dynamics simulation
- of defect formation and precipitation
- in heavily carbon doped silicon
- }\\}
- \vspace{16pt}
- \textsc{\Large F. Zirkelbach$^1$, J. K. N. Lindner$^1$,
- K. Nordlund$^2$, B. Stritzker$^1$}\\
- \vspace{16pt}
- $^1$ Experimentalphysik IV, Institut f\"ur Physik, Universit\"at Augsburg,\\
- Universit\"atsstr. 1, D-86135 Augsburg, Germany\\
- \vspace{16pt}
- $^2$ Accelerator Laboratory, Department of Physical Sciences,
- University of Helsinki,\\
- Pietari Kalmink. 2, 00014 Helsinki, Finland\\
- \vspace{16pt}
- {\scriptsize Corresponding author: Frank Zirkelbach
- <frank.zirkelbach@physik.uni-augsburg.de>}
-\end{center}
-
-%\vspace{24pt}
-
-\section*{Abstract}
-The precipitation process of silicon carbide in heavily carbon doped silicon is not yet fully understood.
-High resolution transmission electron microscopy observations suggest that in a first step carbon atoms form C-Si dumbbells on regular Si lattice sites which agglomerate into large clusters.
-In a second step, when the cluster size reaches a radius of a few $nm$, the high interfacial energy due to the SiC/Si lattice misfit of almost 20\% is overcome and the precipitation occurs.
-By simulation, details of the precipitation process can be obtained on the atomic level.
-A newly parametrized Tersoff-like bond order potential is used to model the system appropriately.
-First results gained by molecular dynamics simulations using this potential are presented.
-\\\\
-{\bf Keywords:} Silicon, carbon, silicon carbide, nucleation, defect formation,
- molecular dynamics simulations
-
-\section*{Introduction}
-Understanding the precipitation process of cubic silicon carbide (3C-SiC) in heavily carbon doped silicon will enable significant technological progress in thin film formation of the important wide band gap semiconductor material SiC \cite{edgar92}.
-On the other hand it will likewise offer perspectives for processes which rely upon prevention of precipitation events, e.g. the fabrication of strained, pseudomorphic $\text{Si}_{1-y}\text{C}_y$ heterostructures \cite{}.
-
-Epitaxial growth of 3C-SiC films is achieved either by ion beam synthesis (IBS) \cite{lindner02} and chemical vapour deposition (CVD) or molecular beam epitaxy (MBE) techniques.
-While in CVD and MBE surface effects need to be taken into account, SiC formation during IBS takes place in the bulk of the Si crystal.
-In the present work the simulation tries to realize conditions which hold for the ion implantation process.
-
-First of all a picture of the supposed precipitation event is presented.
-Afterwards the applied simulation sequences are discussed.
-Finally first results gained by simulation are presented.
-
-\section*{Supposed conversion mechanism}
-Silicon has diamond structure and thus is composed of two fcc lattices which are displaced by one quarter of the volume diagonal.
-3C-SiC grows in zincblende structure, i.e. is also composed of two fcc lattices out of which one is occupied by Si the other by C atoms.
-The length of four lattice constants of Si is approximately equal to the length of five 3C-SiC lattice constants ($4a_{\text{Si}}\approx 5a_{\text{3C-SiC}}$) resulting in a lattice misfit of almost 20\%.
-Due to this the silicon atomic density of 3C-SiC is slightly lower than the one of pure Si.
-
-%\begin{figure}[!h]
-% \begin{center}
-% \begin{minipage}{5.5cm}
-% \includegraphics[width=5cm]{sic_prec_seq_01_s.eps}
-% \end{minipage}
-% \begin{minipage}{5.5cm}
-% \includegraphics[width=5cm]{sic_prec_seq_02_s.eps}
-% \end{minipage}
-% \begin{minipage}{5.5cm}
-% \includegraphics[width=5cm]{sic_prec_seq_03_s.eps}
-% \end{minipage}
-% \caption{Schematic of the supposed conversion mechanism of highly C (${\color{red}\Box}$) doped Si (${\color{black}\bullet}$) into SiC ($_{\color{black}\bullet}^{{\color{red}\Box}}$) and residual Si atoms ($\circ$). The figure shows the dumbbell formation (left), the agglomeration into clusters (middle) and the situation after precipitation (right).}
-% \end{center}
-%\end{figure}
-There is a supposed conversion mechanism of heavily carbon doped Si into SiC \cite{werner97}.
-As concluded by high resolution transmission electron microscopy introduced carbon atoms form C-Si dumbbells on regular Si lattice sites.
-The dumbbells agglomerate into large clusters, called embryos.
-Finally, when the cluster size reaches a critical radius of 2 to 4 nm, the high interfacial energy due to the 3C-SiC/Si lattice misfit is overcome and precipitation occurs.
-Due to the slightly lower silicon density of 3C-SiC excessive silicon atoms exist which will most probably end up as self-interstitials in the silicon matrix since there is more space than in 3C-SiC.
-
-Thus, in addition to the precipitation event itself, knowledge of C and Si interstitials in Si are of great interest in order to investigate the precipitation of heavily C doped Si into SiC.
-%Additionaly the influence of interstitials on atomic diffusion is investigated.
-
-\section*{Simulation sequences}
-A molecular dynamics simulation approach is used to examine the steps involved in the precipitation process.
-For integrating the equations of motion the velocity verlet algorithm \cite{verlet67} with a timestep of $1\, fs$ is adopted.
-The interaction of the silicon and carbon atoms is realized by a newly parametrized Tersoff-like bond order potential \cite{albe_sic_pot}.
-Since temperature and pressure of the system is kept constant in experiment the isothermal-isobaric NPT ensemble is chosen for the simulation.
-Coupling to the heat bath is achieved by the Berendsen thermostat \cite{berendsen84} with a time constant $\tau_T=100\, fs$.
-The pressure is scaled by the Berendsen barostat \cite{berendsen84} again using a timeconstant of $\tau_P=100\, fs$ and a bulk modulus of $100\, GPa$ for silicon.
-To exclude surface effects periodic boundary conditions are applied.
-
-\begin{figure}[!h]
- \begin{center}
- \includegraphics[width=8cm]{unit_cell_s.eps}
- \caption{Insertion positions for the tetrahedral (${\color{red}\triangleleft}$), hexagonal (${\color{green}\triangleright}$) and <110> dumbbell (${\color{magenta}\Box}$) interstitial configuration.}
- \end{center}
-\end{figure}
-To investigate the interstitial configurations of C and Si in Si, a simulation volume of 9 silicon unit cells in each direction is used.
-The temperature is set to $T=0\, K$.
-The insertion positions are illustrated in Fig. 2.
-In separated simulation runs a carbon and a silicon atom respectively is inserted at the tetrahedral $(0,0,0)$ (${\color{red}\triangleleft}$), hexagonal $(-1/8,-1/8,1/8)$ (${\color{green}\triangleright}$), nearby dumbbell $(-1/8,-1/8,-1/4)$ (${\color{magenta}\Box}$) and at random positions (in units of the silicon lattice constant) where the origin is located in the centre of the unit cell.
-In order to avoid too high potential energies in the case of the dumbbell configuration the nearest silicon neighbour atom is shifted to $(-3/8,-3/8,-1/4)$ ($\circ$).
-The energy introduced into the system is scaled out within a relaxation phase of $2\, ps$.
-
-%The same volume is used to investigate diffusion.
-%Different amounts of silicon atoms are inserted at random positions within a centered region of $11 \,\textrm{\AA}$ in each direction.
-%Insertion events are carried out step by step maintaining a constant system temperature of $450\, ^{\circ} \textrm{C}$.
-%Finally a single carbon atom is inserted at a random position within the unit cell located in the middle of the simulation volume.
-%The simulation is proceeded for another $30\, ps$.
-
-For the simulations aiming to reproduce a precipitation process the volume is 31 silicon lattice constants in each direction.
-The system temperature is set to $450\, ^{\circ} \textrm{C}$.
-$6000$ carbon atoms (the amount necessary to form a 3C-SiC precipitate with a radius of 3 nm) are consecutively inserted in a way to keep constant the system temperature.
-Precipitation is examined for three insertion volumes which differ in size.
-The whole simulation volume $V_1$, the volume corresponding to the size of a minimal SiC precipitate $V_2$ and the volume containing the amount of silicon necessary for the formation of such a minimal precipitate $V_3$ are examined.
-The two latter ones are accomplished since no long range diffusion of the carbon atoms is expected at this temperature.
-Following the insertion procedure the system is cooled down to $20\, ^{\circ} \textrm{C}$.
-
-\section*{Results}
-
-The tetrahedral and the <110> dumbbell self-interstitial configurations can be reproduced as observed in \cite{albe_sic_pot}.
-The formation energies are $3.4\, eV$ and $4.4\, eV$ respectively.
-However the hexagonal one is not stable opposed to what is presented in \cite{albe_sic_pot}.
-The atom moves towards an energetically more favorable position very close to the tetrahedral one but slightly displaced along the three coordinate axes.
-The formation energy of $4.0\, eV$ of this type of interstitial equals the result obtained in the reference for the hexagonal one.
-The same type of interstitial may arise using random insertions.
-In addition variations exist in which the displacement is only along two axes ($E_f=3.8\, eV$) or along a single axis ($E_f=3.6\, eV$) succesively approximating the tetrahedral configuration and formation energy.
-
-The tetrahedral and <110> dumbbell carbon interstitial configurations are stable.
-The formation energies are $2.7\, eV$ and $1.8\, eV$ respectively.
-Again the hexagonal one is found to be unstable.
-The interstitial atom moves to the more favorable <100> dumbbell position which has a formation energy of $0.5\, eV$.
-The interstitial configuration is shown in Fig. 2.
-There is experimental evidence \cite{watkins76} of the existence of this configuration.
-It is frequently generated and has the lowest formation energy of all the defects observed in all the simulation runs in which carbon is inserted at random positions.
-Fig. 3 schematically displays the <110> dumbbell configuration including the displacements relative to their initial positions and resulting new Si-Si and C-Si pair distances.
-
-\begin{figure}[!h]
- \begin{center}
- \includegraphics[width=8cm]{c_in_si_int_001db_0.eps}
- \caption{Position of a <100> dumbbell carbon interstitial in silicon.
- Only bonds of the carbon interstitial atom are shown.}
- \end{center}
-\end{figure}
-\begin{figure}[!h]
- \begin{center}
- \includegraphics[width=16cm]{100-c-si-db.eps}
- \caption{Schematic of the <100> C-Si dumbbell configuration.
- Displacements of the atoms relative to their initial position are given.
- The displacement of the carbon atom is relative to the initial position of atom 1.
- New resulting Si-Si and C-Si pair distances for the atoms shown in the schematic and the distances to Si' atoms outside of the displayed region are recorded.}
- \end{center}
-\end{figure}
-
-%\begin{figure}[!h]
-% \begin{center}
-% \includegraphics[width=12cm]{../plot/diff_dep.ps}
-% \caption{Diffusion coefficients of a single carbon atom for different amount of Si selft interstitials}
-% \end{center}
-%\end{figure}
-%The influence of Si self-interstitials on the diffusion of a single carbon atom is displayed in Fig. 3.
-%Diffusion coefficients for different amount of Si self-interstitials are shown.
-%A slight increase is first observed in the case of 30 interstitial atoms.
-%Further increasing the amount of interstitials leads to a tremendous decay of the diffusion coeeficient.
-%Generally there is no long range diffusion of the carbon atom for a temperature of $450\, ^{\circ} \textrm{C}$.
-%The maximal displacement of the carbon atom relativ to its insertion position is between 0.5 and 0.7 \AA.
-
-\begin{figure}[!h]
- \begin{center}
- \includegraphics[width=12cm]{pc_si-c_c-c.ps}
- \caption{Pair correlation functions for Si-C and C-C bonds.
- Carbon atoms are introduced into the whole simulation volume $V_1$, the region which corresponds to the size of a minimal SiC precipitate $V_2$ and the volume which contains the necessary amount of silicon for such a minimal precipitate $V_2$ respectively.}
- \end{center}
-\end{figure}
-\begin{figure}[!h]
- \begin{center}
- \includegraphics[width=12cm]{pc_si-si.ps}
- \caption{Si-Si pair correlation function for pure Si and Si with 3000 inserted C atoms.
- The inset shows a magnified region between 0.28 and 0.36 nm.}
- \end{center}
-\end{figure}
-Fig. 4 shows resulting pair correlation functions of the simulation runs targeting the observation of precipitation events.
-The contributions of Si-C and C-C pairs are presented separately each of them displaying the pair correlation for the three different volumes $V_1$, $V_2$ and $V_3$ (as explained above) exposed to carbon insertion.
-Results show no signigicant difference between $V_1$ and $V_2$.
-Si-Si pairs for the case of 3000 inserted C atoms inserted into $V_2$ and a reference function for pure Si are displayed in Fig. 5.
-
-The amount of C-C bonds for $V_1$ are much smaller than for $V_2$ and $V_3$ since carbon atoms are spread over the total simulation volume which means that there are only 0.2 carbon atoms per silicon unit cell on average.
-The first C-C peak appears at about 0.15 nm.
-This is comparable to the nearest neighbour distance for graphite or diamond.
-It is assumed that these carbon atoms form strong C-C bonds, which is supported by a decrease of the total energy during carbon insertion for the $V_2$ and $V_3$ in contrast to the $V_3$ simulation.
-
-The peak at 0.31 nm perfectly matches the distance of two carbon atoms in the SiC lattice which in SiC is also expected for the Si-Si bonds.
-After insertion of carbon atoms the Si-Si pair correlation function in fact shows non-zero values in the range of the C-C peak width while the amount of Si pairs at the regular distances at 0.24 and 0.38 nm decreases.
-However no clear peak is observed and random analyses of configurations in which distances around 0.3 nm appear, i.e. visualization of such atom pairs, identify <100> C-Si dumbbells to be responsible for stretching the Si-Si next neighbour distance for low concentrations of carbon, i.e. for the $V_1$ and early stages of $V_2$ and $V_3$ simulation runs.
-This excellently agrees with the calculation for a single <100> dumbbell ($r(13)$ in Fig. 4).
-For higher carbon concentrations the defect concentration is likewise increased and a considerable amount of damage is introduced into the inserted volume.
-Damage and superposition of defects generate new displacement arrangements which become hard to categorize and trace and obviously lead to a broader distribution of pair distances.
-The slightly higher amount and intense increase of Si-Si pairs at distances smaller 0.31 nm is probably due to the Si-Si cutoff radius of 0.296 nm.
-The cutoff function causes artificial forces pushing the Si atoms out of the cutoff region.
-By again visualizing the C-C atom pairs with distances of 0.31 nm concatenated, differently oriented <100> dumbbell interstitials are frequently observed.
-Since dumbbells of this type with different orientations are perpendicularly aligned the C atoms are displaced along the plane diagonal of the original lattice.
-One might now assume for the precipitation process, that C atoms are arrenged first and at a later point pull the Si atoms into the right configuration.
-In this way the hkl planes of the SiC in Si and the Si matrix would have equal orientations which is supported by experimental transmission electron microscopy data \cite{}.
-\\\\
-Ab hier weiter ...
-\\\\
-On the other hand the Si-C pair correlation function indicates formation of SiC bonds with an increased crystallinity for the simulation in which carbon is inserted into the whole simulation volume.
-There is more carbon forming Si-C bonds than C-C bonds.
-This gives suspect to the competition of Si-C and C-C bond formation in which the predominance of either of them depends on the method handling carbon insertion.
-
-\section*{Summary}
-The supposed conversion mechanism of heavily carbon doped silicon into silicon carbide is presented.
-Molecular dynamics simulation sequences to investigate interstitial configurations
-%, the influence of interstitials on the atomic diffusion
-and the precipitation of SiC are explained.
-The <100> C-Si dumbbel is reproduced and is the energetically most favorable configuration observed by simulation.
-%The influence of silicon self-interstitials on the diffusion of a single carbon atom is demonstrated.
-Two competing bond formations, either Si-C or C-C, seem to coexist, where the strength of either of them depends on the size of the region in which carbon is introduced.
-
-\bibliography{../../bibdb/bibdb}
-
-\end{document}
projects / lectures / latex.git / commitdiff