-\documentclass[prb,twocolumn,superscriptaddress,a4paper,showkeys,showpacs]{revtex4}\r
-\usepackage{graphicx}\r
-\usepackage{subfigure}\r
-\usepackage{dcolumn}\r
-\usepackage{booktabs}\r
-\usepackage{units}\r
-\usepackage{amsmath}\r
-\usepackage{amsfonts}\r
-\usepackage{amssymb}\r
-\r
-\r
-\begin{document}\r
-\r
-\title{Description of Defects in Carbon implanted Silicon -- a comparison of classical potentials and first principles methods}\r
-\author{F. Zirkelbach} \author{B. Stritzker}\r
-\affiliation{Experimentalphysik IV, Universit\"at Augsburg, D-86153 Augsburg, Germany}\r
-\author{K. Nordlund}\r
-\affiliation{Accelerator Laboratory, University of Helsinki, 00014 Helsinki, Finland}\r
-\author{J. K. N. Lindner}\r
-\affiliation{Experimentelle Physik, Universit\"at Paderborn, 33095 Paderborn, Germany}\r
-\author{W. G. Schmidt} \author{E. Rauls}\r
-\affiliation{Theoretische Physik, Universit\"at Paderborn, 33095 Paderborn, Germany}\r
-\r
-\begin{abstract}\r
-We present a comparative theoretical investigation of carbon interstitials in silicon.\r
-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
-\end{abstract}\r
-\r
-\keywords{point defects, migration, interstitials, first principles calculations, classical potentials }\r
-\pacs{ find out later... }\r
-\maketitle\r
-\r
-% --------------------------------------------------------------------------------\r
-\section{Introduction}\r
-\r
-%Frank: Idea: description of 3C-SiC-precipitation in C-implanted silicon.\\ \r
-% cite and describe briefly experimental work - why is this material important/better than other SiC). \\ \r
-% Describe the precipitation process in brief.\\ \r
-% 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
-\r
-% there should be a short motivation for the material system!\r
-Silicon carbide (SiC) has a number of remarkable physical and chemical properties.\r
-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,morkoc94}, which are moreover deployable in harsh and radiation-hard environments\cite{capano97}.\r
-% there are different polytpes with different properties and 3c-sic in special\r
-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
-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
-\r
-% (thin films of) 3c-sic can be produced by ibs\r
-Next to the fabrication of 3C-SiC layers by chemical vapor deposition (CVD) and molecular beam epitaxy (MBE) on hexagonal SiC\cite{powell90,fissel95,fissel95_apl} and Si\cite{nishino83,nishino87,kitabatake93,fissel95_apl} substrates, high-dose carbon implantation into crystalline silicon (c-Si) with subsequent or in situ annealing was found to result in SiC microcrystallites in Si\cite{borders71}.\r
-% maybe split CVD and MBE from IBS and explain remaining problems:\r
-% - on 6H-SiC: twin boundaries\r
-% - on Si: structural defects due to thermal conductivity and lattice mismatch\r
-Utilized and enhanced, ion beam synthesis (IBS) has become a promising method to form thin SiC layers of high quality exclusively of the 3C polytype embedded in and epitactically aligned to the Si host featuring a sharp interface\cite{lindner99,lindner01,lindner02}.\r
-% precipitation model\r
-However, only little is known on the SiC conversion in C implanted Si.\r
-High resolution transmission electron microscopy (HREM) studies\cite{werner96,werner97,lindner99_2} suggest the formation of C-Si dimers (dumbbells) on regular Si lattice sites, which agglomerate into large clusters indicated by dark contrasts and otherwise undisturbed Si lattice fringes in HREM.\r
-Once a critical radius of 2 nm to 4 nm is reached a topotactic transformation into a 3C-SiC precipitate occurs.\r
-The transformation is manifested by the disappearance of dark contrasts in favor of Moir\'e patterns due to the lattice mismatch of 20 \% of the 3C-SiC precipitate and c-Si.\r
-The insignificantly lower Si density of SiC ($\approx 4$ \%) compared to c-Si results in the emission of only a few excess Si atoms.\r
-% motivation to understand the precipitation and link to atomistic simulations\r
-A detailed understanding of the underlying processes will enable significant technological progress in 3C-SiC thin film formation and likewise offer perspectives for processes which rely upon prevention of precipitation events, e.g. the fabrication of strained pseudomorphic Si$_{1-y}$C$_y$ heterostructures\cite{strane96,laveant2002}.\r
-\r
-Atomistic simulations offer a powerful tool to study materials on a microscopic level providing detailed insight not accessible by experiment.\r
-Relevant structures consisting of $\approx 10^4$ atoms for the nanocrystal and even more atoms for a reasonably sized Si host matrix are too large to be completely described by high accuracy quantum mechanical methods.\r
-Directly modelling the dynamics of the processes mentioned above almost inevitably requires the atomic interaction to be described by less accurate though computationally more efficient classical potentials.\r
-The most common empirical potentials for covalent systems are the Stillinger-Weber\cite{stillinger85}, Brenner\cite{brenner90}, Tersoff\cite{tersoff_si3} and environment-dependent interatomic (EDIP)\cite{bazant96,bazant97,justo98} potential.\r
-From the mentioned potentials, until recently\cite{lucas10}, a parametrization to describe the C-Si multicomponent system did only exist for the Tersoff potential\cite{tersoff_m}.\r
-Description successful ... and failed ...\r
-Whether such potentials are appropriate for the description of the physical problem has, however, to be verified first by applying classical and quantum-mechanical methods to relevant processes that can be treated by both methods.\r
-In this work, the applicability of a Tersoff-like bond order potential\cite{albe_sic_pot} to some basic processes involved in the initially mentioned SiC precipitation mechanism is investigated by comparing results gained by classical and ab inito calculations.\r
-\r
-In the following a comparative investigation of density functional theory (DFT) studies and classical potential calculations of the structure, energetics and mobility of carbon defects in silicon is presented.\r
-\r
-% --------------------------------------------------------------------------------\r
-\section{Methodology}\r
-% ----- DFT ------\r
-The first-principles DFT calculations were performed with the plane-wave based \r
-Vienna Ab-initio Simulation Package (VASP)\cite{kresse96}. The Kohn-Sham equations were solved \r
-using the generalized-gradient XC-functional approximation proposed by Perdew and \r
-Wang (GGA-PW91)\cite{perdew92}. \r
-The electron-ion interaction was described by the projector-augmented wave (PAW) method\cite{bloechel94,kresse99}. \r
-In the PAW data scalar relativistic corrections are contained. Throughout this work an \r
-energy cut-off of \unit[300]{eV} was used to expand the wave functions into the plane-wave basis. \r
-For the sampling of the Brillouin zone, only the $\Gamma$-point was used. \r
-The defect structures and the migration paths were modelled in cubic supercells containing 216 Si-atoms. Spin polarization has been \r
-fully accounted for. \r
-\r
-% ------ Albe potential ---------\r
-%% Frank: Setup/short description of the potential ?\r
-For the calculations with the classical potentials... \r
-\r
-\r
-\section{Results}\r
-\r
-After ion implantation, carbon interstitials are the most common defects in the silicon sample. Their mobility is the \r
-crucial quantity to be investigated. We thus started our comparative investigations by comparing the stability and the \r
-mobility of an isolated carbon interstitial in silicon bulk in the various possible geometries it can take. \r
-\r
-\subsection{Carbon interstitials in various geometries}\r
-\r
-Several geometries have been calculated to be stable for the carbon interstitial. Fig.\ref{fig:interstitials} shows all these \r
-structures. However, there are some discrepancies between the results from classical potential calculations and those obtained \r
-from first principles. \r
-Table \ref{table:formation} summarizes the formation energies of the interstitial geometries for both methods used in this work \r
-and compares the results to literature values. (...check references for more data, ..) \r
-\r
-% Tables: like in the talk, but add further literature data and give the references/citations (also to bibliography \r
-% at the end!) \r
-%\begin{figure}\r
-%\includegraphics[width=1.0\columnwidth]{models.eps}\r
-%\caption{\label{fig:interstitials} Molecular model of the possible carbon interstitials. }\r
-%\end{figure} \r
-\r
-While the Albe potential predicts ... as stable, DFT does not. ...(further comparisons, trend "too high/low" E-formation,...)... \r
- Nevertheless, both methods predict the (110) dumb bell configuration to be the most stable... (?) \r
-\r
-\r
-\subsection{Mobility}\r
-A measure for the mobility of the interstitial carbon is the activation energy for the migration path from one stable \r
-position to another. The stable defect geometries have been discussed in the previous subsection. We now investigate \r
-the migration from the most stable structure (...should be named somehow...) on one site of the silicon host lattice to \r
-a neighbored site. \r
-On the lowest energy path (first principles), the carbon atom starts to move along (110)..(check that!)... The center of the line connecting \r
-initial and final structure has been found to be a local minimum and not a saddle point as could be expected. The two \r
-saddle points shortly before and behind this local minimum are slightly displaced out of the (110) plane by ... {\AA}. ..(check that!)..\r
-This path is not surprising -- a similar behavior was e.g. found earlier for the carbon split interstitial \cite{rauls03a} and the phosphorus \r
-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
-configuration to one at the local minimum. By this, the energy of the local minimum is lowered by 0.3 eV (... check it!!..). \r
-%\begin{figure}\r
-%\includegraphics[width=1.0\columnwidth]{path-DFT.eps}\r
-%\caption{\label{fig:path-DFT} Energy of the carbon interstitial during migration from ... to ... calculated from first principles. The \r
-% activation energy of 0.9 eV (?) agrees well with experimental findings (0.7-0.9 eV?). }\r
-%\end{figure} \r
-Fig.\ref{fig:path-DFT} shows the energy along this lowest energy migration path. The activation energy of 0.9 eV (?) agrees well \r
-with experimental findings (0.7-0.9 eV?).\r
-\r
-Calculations with the Albe potential yield a different picture. \r
-%\begin{figure}\r
-%\includegraphics[width=1.0\columnwidth]{path-Albe.eps}\r
-%\caption{\label{fig:path-Albe} Energy of the carbon interstitial during migration from ... to ... calculated using the classical potential \r
-% method. Here, the activation energy is 2.2 eV (?). }\r
-%\end{figure} \r
-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
-is way too high to explain the experimental findings (0.7-0.9 eV?). (...further discussion...) \r
-\r
-\section{Discussion}\r
-The first principles results are in good agreement to previous work on this subject \cite{joannopoulos,xyz} (...add some references!...). With an \r
-activation energy of 0.9 eV, the carbon interstitial can be expected to be mobile at temperatures in the range of... (?). \r
-The description of the same processes obviously fails if we use the classical potential method. \r
-Already the geometry of the most stable dumb bell configuration differs considerably from that of the first principles calculated \r
-structure. (..... add description, the two main angles and bond lengths and an explanation...)\r
-Formation energies are throughout too high... (...reason?...)\r
-\r
-A reason for this failure of the classical description is most likely... (cut-off, neglect of quantum mechanical effects,...)\r
-\r
-\r
-\section{Summary}\r
-In summary, we have shown that ab initio calculations are very close to the results expected from experimental data. \r
-Furthermore, they agree well with other theoretical results. (...some results - later...) The classical potentials, however, fail to describe the \r
-selected processes. This has been shown to have two reasons, i.e. the artificial cut-off of the next nearest neighbor \r
-interaction on the one hand, on the other hand the quantum mechanical effects which are crucial in the problem under study. \r
-\r
-% ----------------------------------------------------\r
-\section*{Acknowledgment}\r
-The calculations were done using grants of computer time from the \r
-Paderborn Center for Parallel Computing (PC$^2$) and the \r
-H\"ochstleistungs-Rechenzentrum Stuttgart. The Deutsche \r
-Forschungsgemeinschaft is acknowledged for financial support.\r
-\r
-% --------------------------------- references -------------------\r
-\begin{thebibliography}{99}\r
-\bibitem{kresse96} G. Kresse and J. Furthm\"uller, \r
- Comput. Mater. Sci. {\bf 6}, 15 (1996). \r
-\bibitem{perdew92} J. P. Perdew, J. A. Chevary, S. H. Vosko, K. A. Jackson, M. R. Pederson, D. J. Singh and C. Fiolhais, \r
- Phys. Rev. B {\bf 46}, 6671 (1992). \r
-\bibitem{ceperley80} D. M. Ceperley and B. J. Alder, \r
- Phys. Rev. Lett. {\bf 45}, 556 (1980). \r
-\bibitem{perdew81} J. P. Perdew and A. Zunger, \r
- Phys. Rev. B {\bf 23}, 5048 (1981). \r
-\bibitem{bloechel94} P. E. Bl\"ochl, \r
- Phys. Rev. B {\bf 50}, 17953 (1994).\r
-\bibitem{kresse99} G. Kresse and D. Joubert, \r
- Phys. Rev. B {\bf 59}, 1758 (1999).\r
-\bibitem{monk76} H. J. Monkhorst and J. D. Pack, \r
- Phys. Rev. B {\bf 13}, 5188 (1976). \r
-\bibitem{albe} Albe potential\r
-\bibitem{stillinger} Stillinger-Weber potential \r
-\bibitem{joannopoulos} Joannopoulos\r
-\bibitem{xyz} who else? \r
-\bibitem{rauls03a} E. Rauls, A. Gali, P. De´ak, and Th. Frauenheim, Phys. Rev. B, 68, 155208 (2003).\r
-\bibitem{rauls03b} E. Rauls, U. Gerstmann, H. Overhof, and Th. Frauenheim, Physica B, Vols. 340-342, p. 184-189 (2003). \r
-\bibitem{gerstmann03} U. Gerstmann, E. Rauls, Th. Frauenheim, and H. Overhof, Phys. Rev. B, 67, 205202, (2003).\r
-\r
-\end{thebibliography}\r
-\r
-\r
-\end{document}\r
-\r
+\documentclass[prb,twocolumn,superscriptaddress,a4paper,showkeys,showpacs]{revtex4}
+\usepackage{graphicx}
+\usepackage{subfigure}
+\usepackage{dcolumn}
+\usepackage{booktabs}
+\usepackage{units}
+\usepackage{amsmath}
+\usepackage{amsfonts}
+\usepackage{amssymb}
+
+% additional stuff
+\usepackage{miller}
+
+\begin{document}
+
+% ref mod: no capital letters (by editors)
+\title{Defects in carbon implanted silicon calculated by classical potentials and first principles methods}
+\author{F. Zirkelbach}
+\author{B. Stritzker}
+\affiliation{Experimentalphysik IV, Universit\"at Augsburg, 86135 Augsburg, Germany}
+\author{K. Nordlund}
+\affiliation{Department of Physics, University of Helsinki, 00014 Helsinki, Finland}
+\author{J. K. N. Lindner}
+\author{W. G. Schmidt}
+\author{E. Rauls}
+\affiliation{Department Physik, Universit\"at Paderborn, 33095 Paderborn, Germany}
+
+\begin{abstract}
+A comparative theoretical investigation of carbon interstitials in silicon is presented.
+% ref mod: language - put aside
+%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.
+Calculations using classical potentials are compared to 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.
+In contrast to previous studies, the present first principles calculations of the interstitial carbon migration path yield an activation energy that
+excellently matches the experiment.
+%ref mod: language - imperative
+%The bond-centered interstitial configuration shows a net magnetization of two electrons, illustrating the imperative of spin polarized calculations.
+The bond-centered interstitial configuration shows a net magnetization of two electrons, illustrating the need for spin polarized calculations.
+
+\end{abstract}
+
+\keywords{point defects, migration, interstitials, first principles calculations, classical potentials}
+\pacs{61.72.uf,66.30.-h,31.15.A-,34.20.Cf}
+\maketitle
+
+% --------------------------------------------------------------------------------
+\section{Introduction}
+
+Silicon carbide (SiC) has a number of remarkable physical and chemical properties.
+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,morkoc94}, which are moreover deployable in harsh and radiation-hard environments\cite{capano97}.
+SiC, which forms fourfold coordinated mostly 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}.
+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.
+
+Thin films of 3C-SiC can be fabricated by chemical vapor deposition (CVD) and molecular beam epitaxy (MBE) on hexagonal SiC\cite{powell90,fissel95,fissel95_apl} and Si\cite{nishino83,nishino87,kitabatake93,fissel95_apl} substrates.
+Next to these methods, high-dose carbon implantation into crystalline silicon (c-Si) with subsequent or in situ annealing was found to result in SiC microcrystallites in Si\cite{borders71}.
+% ref mod: language - omit first part (pretty much liked it!)
+%Utilized and enhanced, ion beam synthesis (IBS) has become a promising method to form thin SiC layers of high quality exclusively of the 3C polytype embedded in and epitactically aligned to the Si host featuring a sharp interface\cite{lindner99,lindner01,lindner02}.
+Ion beam synthesis (IBS) has become a promising method to form thin SiC layers of high quality exclusively of the 3C polytype embedded in and epitactically aligned to the Si host featuring a sharp interface\cite{lindner99,lindner01,lindner02}.
+However, only little is known about the SiC conversion in C implanted Si.
+High resolution transmission electron microscopy (HREM) studies\cite{werner96,werner97,lindner99_2} suggest the formation of C-Si dimers (dumbbells) on regular Si lattice sites, which agglomerate into large clusters indicated by dark contrasts and otherwise undisturbed Si lattice fringes in HREM.
+Once a critical radius of 2 nm to 4 nm is reached, a topotactic transformation into a 3C-SiC precipitate occurs.
+The transformation is manifested by the disappearance of dark contrasts in favor of Moir\'e patterns due to the lattice mismatch of \unit[20]{\%} of the 3C-SiC precipitate and c-Si.
+The insignificantly lower Si density of SiC ($\approx \unit[4]{\%}$) compared to c-Si results in the emission of only a few excess Si atoms.
+A detailed understanding of the underlying processes will enable significant technological progress in 3C-SiC thin film formation and likewise offer perspectives for processes which rely upon prevention of precipitation events, e.g. the fabrication of strained pseudomorphic Si$_{1-y}$C$_y$ heterostructures\cite{strane96,laveant2002}.
+
+Atomistic simulations offer a powerful tool to study materials on a microscopic level providing detailed insight not accessible by experiment.
+Relevant structures consisting of $\approx 10^4$ atoms for the nanocrystal and even more atoms for a reasonably sized Si host matrix are too large to be completely described by high accuracy quantum mechanical methods.
+Directly modelling the dynamics of the processes mentioned above almost inevitably requires the atomic interaction to be described by less accurate though computationally more efficient classical potentials.
+% ref mod: language EDIP
+%The most common empirical potentials for covalent systems are the Stillinger-Weber\cite{stillinger85} (SW), Brenner\cite{brenner90}, Tersoff\cite{tersoff_si3} and environment-dependent interatomic (EDIP)\cite{bazant96,bazant97,justo98} potential.
+The most common empirical potentials for covalent systems are the Stillinger-Weber\cite{stillinger85} (SW), Brenner\cite{brenner90}, Tersoff\cite{tersoff_si3} and environment-dependent interatomic potential (EDIP)\cite{bazant96,bazant97,justo98}.
+Until recently\cite{lucas10}, a parametrization to describe the C-Si multicomponent system within the mentioned interaction models did only exist for the Tersoff\cite{tersoff_m} and related potentials, e.g. the one by Gao and Weber\cite{gao02}.
+Whether such potentials are appropriate for the description of the physical problem has, however, to be verified first by applying classical and quantum-mechanical methods to relevant processes that can be treated by both methods.
+For instance, a comparison of empirical potential molecular dynamics (MD) and density functional theory (DFT) calculations showed that SW is best suited for simulations of dislocation nucleation processes\cite{godet03} and threshold displacement energy calculations\cite{holmstroem08} in Si important in ion implantation, while the Tersoff potential yielded a qualitative agreement for the interaction of Si self-interstitials with substitutional C\cite{mattoni2002}.
+Antisite pairs and defects in SiC have been investigated, both classically\cite{posselt06,gao04} employing the Gao/Weber potential\cite{gao02} and quantum-mechanically\cite{rauls03a,gao07,gao2001,bockstedte03}, which, both, agree very well with experimental results\cite{gali03,chen98,weber01}.
+An extensive comparison\cite{balamane92} concludes that each potential has its strengths and limitations and none of them is clearly superior to others.
+Despite their shortcomings these potentials are assumed to be reliable for large-scale simulations\cite{balamane92,huang95,godet03} on specific problems under investigation providing insight into phenomena that are otherwise not accessible by experimental or first principles methods.
+Remaining shortcomings have frequently been resolved by modifying the interaction\cite{tang95,gao02a,mattoni2007} or extending it\cite{devanathan98_2} with data gained from ab initio calculations\cite{nordlund97}.
+
+In this work, the applicability of a Tersoff-like bond order potential\cite{albe_sic_pot} to basic processes involved in the initially mentioned SiC precipitation mechanism has been investigated by comparing results gained by classical and ab initio calculations.
+In the following a comparative investigation of density functional theory studies and classical potential calculations of the structure, energetics and mobility of carbon defects in silicon is presented.
+
+% --------------------------------------------------------------------------------
+\section{Methodology}
+% ----- DFT ------
+The first principles DFT calculations have been performed with the plane-wave based Vienna Ab-initio Simulation Package (VASP)\cite{kresse96}.
+The Kohn-Sham equations were solved using the generalized-gradient XC-functional approximation proposed by Perdew and Wang (GGA-PW91)\cite{perdew86,perdew92}.
+The electron-ion interaction is described by norm-conserving ultra-soft pseudopotentials\cite{hamann79} as implemented in VASP\cite{vanderbilt90}.
+Throughout this work, an energy cut-off of \unit[300]{eV} was used to expand the wave functions into the plane-wave basis.
+Sampling of the Brillouin zone was restricted to the $\Gamma$-point. Some test calculations were performed with a $(2\times 2\times 2)$ k-point set.
+The defect structures and the migration paths have been modelled in cubic supercells containing 216 Si atoms.
+The ions and cell shape were allowed to change in order to realize a constant pressure simulation. Spin polarization has been fully accounted for.
+
+% ------ Albe potential ---------
+For the classical potential calculations, a supercell of 9 Si lattice constants in each direction consisting of 5832 Si atoms has been used.
+A Tersoff-like bond order potential by Erhart and Albe (EA)\cite{albe_sic_pot} has been utilized, which accounts for nearest neighbor interactions only realized by a cut-off function dropping the interaction to zero in between the first and second next neighbor distance.
+% ref mod: extension for short distances
+The potential was used as is, i.e. without any repulsive potential extension at short interatomic distances.
+% ref mod: time constants
+%Constant pressure simulations are realized by the Berendsen barostat\cite{berendsen84}.
+Constant pressure simulations are realized by the Berendsen barostat\cite{berendsen84} using a time constant of \unit[100]{fs} and a bulk modulus of \unit[100]{GPa} for Si.
+% ref mod: time constants + language (Verlet)
+%Structural relaxation in the MD run is achieved by the velocity Verlet algorithm\cite{verlet67} and the Berendsen thermostat\cite{berendsen84} with a time constant of \unit[1]{fs} resulting in direct velocity scaling and the temperature set to zero Kelvin.
+Structural relaxation in the MD run is achieved by the Velocity Verlet algorithm\cite{verlet67} and the Berendsen thermostat\cite{berendsen84} with a time constant of \unit[100]{fs} and the temperature set to zero Kelvin.
+Additionally, a time constant of \unit[1]{fs} resulting in direct velocity scaling was used for relaxation within the mobility calculations.
+% ref mod: time step
+A fixed time step of \unit[1]{fs} for integrating the equations of motion was used.
+
+\section{Results}
+
+According to the assumed SiC precipitation model described in the introductory part, carbon interstitial defects form and agglomerate into large clusters.
+Thus, it is of crucial importance to investigate the various possible structures of carbon defects and the mobility of the lowest energy, hence most probable, defect configuration in crystalline silicon.
+
+\subsection{Carbon interstitials in various geometries}
+
+Table~\ref{tab:defects} summarizes the formation energies of defect structures for the EA and DFT calculations performed in this work as well as further results from literature.
+The formation energy $E-N_{\text{Si}}\mu_{\text{Si}}-N_{\text{C}}\mu_{\text{C}}$ is defined in the same way as in the articles used for comparison\cite{tersoff90,dal_pino93} chosing SiC as a reservoir for the carbon impurity in order to determine $\mu_{\text{C}}$.
+Relaxed geometries are displayed in Fig.~\ref{fig:defects}.
+\begin{table*}
+\begin{ruledtabular}
+\begin{tabular}{l c c c c c c}
+ & T & H & \hkl<1 0 0> dumbbell & \hkl<1 1 0> dumbbell & S & BC \\
+\hline
+ Erhart/Albe & 6.09 & 9.05$^*$ & 3.88 & 5.18 & 0.75 & 5.59$^*$ \\
+ VASP & unstable & unstable & 3.72 & 4.16 & 1.95 & 4.66 \\
+ Tersoff\cite{tersoff90} & 3.8 & 6.7 & 4.6 & 5.9 & 1.6 & 5.3 \\
+ ab initio\cite{dal_pino93,capaz94} & - & - & x & - & 1.89 \cite{dal_pino93} & x+2.1 \cite{capaz94} \\
+\end{tabular}
+\end{ruledtabular}
+\caption{Formation energies of carbon point defects in crystalline silicon determined by classical potential and ab initio methods. The formation energies are given in eV. T denotes the tetrahedral, H the hexagonal and BC the bond-centered interstitial configuration. S corresponds to substitutional C. Formation energies for unstable configurations obtained by classical potential MD are marked by an asterisk and determined by using the low kinetic energy configuration shortly before the relaxation into the more favorable configuration starts.}
+\label{tab:defects}
+\end{table*}
+\begin{figure}
+\begin{minipage}[t]{0.32\columnwidth}
+\underline{Tetrahedral}\\
+\includegraphics[width=\columnwidth]{tet.eps}
+\end{minipage}
+\begin{minipage}[t]{0.32\columnwidth}
+\underline{Hexagonal}\\
+\includegraphics[width=\columnwidth]{hex.eps}
+\end{minipage}
+\begin{minipage}[t]{0.32\columnwidth}
+\underline{\hkl<1 0 0> dumbbell}\\
+\includegraphics[width=\columnwidth]{100.eps}
+\end{minipage}\\
+\begin{minipage}[t]{0.32\columnwidth}
+\underline{\hkl<1 1 0> dumbbell}\\
+\includegraphics[width=\columnwidth]{110.eps}
+\end{minipage}
+\begin{minipage}[t]{0.32\columnwidth}
+\underline{Substitutional}\\[0.05cm]
+\includegraphics[width=\columnwidth]{sub.eps}
+\end{minipage}
+\begin{minipage}[t]{0.32\columnwidth}
+\underline{Bond-centered}\\
+\includegraphics[width=\columnwidth]{bc.eps}
+\end{minipage}
+\caption{Configurations of carbon point defects in silicon. The silicon/carbon atoms and the bonds (only for the interstitial atom) are illustrated by yellow/grey spheres and blue lines. Bonds are drawn for atoms located within a certain distance and do not necessarily correspond to chemical bonds.}
+\label{fig:defects}
+\end{figure}
+
+Substitutional carbon (C$_{\text{sub}}$) occupying an already vacant Si lattice site, which is in fact not an interstitial defect, is found to be the lowest configuration with regard to energy for all potential models.
+DFT calculations performed in this work are in good agreement with results obtained by classical potential simulations by Tersoff\cite{tersoff90} and ab initio calculations done by Dal Pino et~al\cite{dal_pino93}.
+% ref mod: typo - underestim_a_tes
+%However, the EA potential dramatically underestimtes the C$_{\text{sub}}$ formation energy, which is a definite drawback of the potential.
+However, the EA potential dramatically underestimates the C$_{\text{sub}}$ formation energy, which is a definite drawback of the potential.
+
+Except for the Tersoff potential the \hkl<1 0 0> dumbbell (C$_{\text{i}}$) is the energetically most favorable interstital configuration, in which the C and Si dumbbell atoms share a Si lattice site.
+This finding is in agreement with several theoretical\cite{burnard93,leary97,dal_pino93,capaz94} and experimental\cite{watkins76,song90} investigations.
+Tersoff as well, considers C$_{\text{i}}$ to be the ground state configuration and believes an artifact due to the abrupt C-Si cut-off used in the potential to be responsible for the small value of the tetrahedral formation energy\cite{tersoff90}.
+It should be noted that EA and DFT predict almost equal formation energies.
+However, there is a qualitative difference: while the C-Si distance of the dumbbell atoms is almost equal for both methods, the position along \hkl[0 0 1] of the dumbbell inside the tetrahedron spanned by the four next neighbored Si atoms differs significantly.
+The dumbbell based on the EA potential is almost centered around the regular Si lattice site as can be seen in Fig.~\ref{fig:defects} whereas for DFT calculations it is translated upwards with the C atom forming an almost collinear bond to the two Si atoms of the top face of the tetrahedron and the bond angle of the Si dumbbell atom to the two bottom face Si atoms approaching \unit[120]{$^\circ$}.
+This indicates predominant sp and sp$^2$ hybridization for the C and Si dumbbell atom respectively.
+Obviously the classical potential is not able to reproduce the clearly quantum-mechanically dominated character of bonding.
+
+Both, EA and DFT reveal the hexagonal configuration unstable relaxing into the C$_{\text{i}}$ ground state structure.
+Tersoff finds this configuration stable, though it is the most unfavorable.
+Thus, the highest formation energy observed by the EA potential is the tetrahedral configuration, which turns out to be unstable in DFT calculations.
+The high formation energy of this defect involving a low probability to find such a defect in classical potential MD acts in concert with finding it unstable by the more accurate quantum-mechnical description.
+
+The \hkl<1 1 0> dumbbell constitutes the second most favorable configuration, reproduced by both methods.
+It is followed by the bond-centered (BC) configuration.
+However, even though EA yields the same difference in energy with respect to the \hkl<1 1 0> defect as DFT does, the BC configuration is found to be a saddle point within the EA description relaxing into the \hkl<1 1 0> configuration.
+Tersoff indeed predicts a metastable BC configuration.
+However, it is not in the correct order and lower in energy than the \hkl<1 1 0> dumbbell.
+Please note, that Capaz et~al.\cite{capaz94} in turn found this configuration to be a saddle point, which is about \unit[2.1]{eV} higher in energy than the C$_{\text{i}}$ configuration.
+This is assumed to be due to the neglection of the electron spin in these calculations.
+Another DFT calculation without fully accounting for the electron spin results in the smearing of a single electron over two non-degenerate states for the BC configuration.
+This problem is resolved by spin polarized calculations resulting in a net spin of one accompanied by a reduction of the total energy by \unit[0.3]{eV} and the transformation into a metastable local minimum configuration.
+All other configurations are not affected.
+
+To conclude, we observed discrepancies between the results from classical potential calculations and those obtained from first principles.
+Within the classical potentials EA outperforms Tersoff and is, therefore, used for further comparative studies.
+Both methods (EA and DFT) predict the \hkl<1 0 0> dumbbell interstitial configuration to be most stable.
+%ref mod: language - energetical order
+%Also the remaining defects and their energetical order are described fairly well.
+Also the remaining defects and their relative energies are described fairly well.
+It is thus concluded that -- so far -- modelling of the SiC precipitation by the EA potential might lead to trustable results.
+
+\subsection{Mobility}
+
+A measure for the mobility of the interstitial carbon is the activation energy for the migration path from one stable position to another.
+The stable defect geometries have been discussed in the previous subsection.
+In the following the migration of the most stable configuration, i.e. C$_{\text{i}}$, from one site of the Si host lattice to a neighboring site has been investigated by both, EA and DFT calculations utilizing the constraint conjugate gradient relaxation technique (CRT)\cite{kaukonen98}.
+Three migration pathways are investigated.
+The starting configuration for all pathways was the \hkl[0 0 -1] dumbbell interstitial configuration.
+In path~1 and 2 the final configuration is a \hkl[0 0 1] and \hkl[0 -1 0] dumbbell interstitial respectively, located at the next neighbored Si lattice site displaced by $\frac{a_{\text{Si}}}{4}\hkl[1 1 -1]$, where $a_{\text{Si}}$ is the Si lattice constant.
+In path~1 the C atom resides in the \hkl(1 1 0) plane crossing the BC configuration whereas in path~2 the C atom moves out of the \hkl(1 1 0) plane.
+Path 3 ends in a \hkl[0 -1 0] configuration at the initial lattice site and, for this reason, corresponds to a reorientation of the dumbbell, a process not contributing to long range diffusion.
+
+\begin{figure}
+\begin{center}
+\includegraphics[width=\columnwidth]{path2_vasp_s.ps}
+\end{center}
+\caption{Migration barrier and structures of the \hkl[0 0 -1] dumbbell (left) to the \hkl[0 -1 0] dumbbell (right) transition as obtained by first principles methods. The activation energy of \unit[0.9]{eV} agrees well with experimental findings of \unit[0.70]{eV}\cite{lindner06}, \unit[0.73]{eV}\cite{song90} and \unit[0.87]{eV}\cite{tipping87}.}
+\label{fig:vasp_mig}
+\end{figure}
+The lowest energy path (path~2) as detected by the first principles approach is illustrated in Fig.~\ref{fig:vasp_mig}, in which the \hkl[0 0 -1] dumbbell migrates towards the next neighbored Si atom escaping the $(1 1 0)$ plane forming a \hkl[0 -1 0] dumbbell.
+The activation energy of \unit[0.9]{eV} excellently agrees with experimental findings ranging from \unit[0.70]{eV} to \unit[0.87]{eV}\cite{lindner06,song90,tipping87}.
+
+\begin{figure}
+\begin{center}
+\includegraphics[width=\columnwidth]{path1_albe_s.ps}
+\end{center}
+\caption{Migration barrier and structures of the bond-centered (left) to \hkl[0 0 -1] dumbbell (right) transition utilizing the classical potential method. Two different pathways are obtained for different time constants of the Berendsen thermostat. The lowest activation energy is \unit[2.2]{eV}.}
+\label{fig:albe_mig}
+\end{figure}
+Calculations based on the EA potential yield a different picture.
+Fig.~\ref{fig:albe_mig} shows the evolution of structure and energy along the lowest energy migration path (path~1) based on the EA potential.
+Due to symmetry it is sufficient to merely consider the migration from the BC to the C$_{\text{i}}$ configuration.
+Two different pathways are obtained for different time constants of the Berendsen thermostat.
+With a time constant of \unit[1]{fs} the C atom resides in the \hkl(1 1 0) plane resulting in a migration barrier of \unit[2.4]{eV}.
+However, weaker coupling to the heat bath realized by an increase of the time constant to \unit[100]{fs} enables the C atom to move out of the \hkl(1 1 0) plane already at the beginning, which is accompanied by a reduction in energy, approaching the final configuration on a curved path.
+The energy barrier of this path is \unit[0.2]{eV} lower in energy than the direct migration within the \hkl(1 1 0) plane.
+It should be noted that the BC configuration is actually not a local minimum configuration in EA based calculations since a relaxation into the \hkl<1 1 0> dumbbell configuration occurs.
+However, investigating further migration pathways involving the \hkl<1 1 0> interstitial did not yield lower migration barriers.
+Thus, the activation energy should at least amount to \unit[2.2]{eV}.
+
+\section{Discussion}
+
+The first principles results are in good agreement to previous work on this subject\cite{burnard93,leary97,dal_pino93,capaz94}.
+The C-Si \hkl<1 0 0> dumbbell interstitial is found to be the ground state configuration of a C defect in Si.
+The lowest migration path already proposed by Capaz et~al.\cite{capaz94} is reinforced by an additional improvement of the quantitative conformance of the barrier height calculated in this work (\unit[0.9]{eV}) with experimentally observed values (\unit[0.70]{eV} -- \unit[0.87]{eV})\cite{lindner06,song90,tipping87}.
+However, it turns out that the bond-centered configuration is not a saddle point configuration as proposed by Capaz et~al.\cite{capaz94} but constitutes a real local minimum if the electron spin is properly accounted for.
+A net magnetization of two electrons, which is already clear by simple molecular orbital theory considerations on the bonding of the sp hybridized C atom, is settled.
+By investigating the charge density isosurface it turns out that the two resulting spin up electrons are localized in a torus around the C atom.
+With an activation energy of \unit[0.9]{eV} the C$_{\text{i}}$ carbon interstitial can be expected to be highly mobile at prevailing temperatures in the process under investigation, i.e. IBS.
+
+We found that the description of the same processes fails if classical potential methods are used.
+Already the geometry of the most stable dumbbell configuration differs considerably from that obtained by first principles calculations.
+The classical approach is unable to reproduce the correct character of bonding due to the deficiency of quantum-mechanical effects in the potential.
+%ref mod: language - energy / order
+%Nevertheless, both methods predict the same type of interstitial as the ground state configuration, and also the order in energy of the remaining defects is reproduced fairly well.
+Nevertheless, both methods predict the same type of interstitial as the ground state configuration.
+Furthermore, the relative energies of the other defects are reproduced fairly well.
+From this, a description of defect structures by classical potentials looks promising.
+% ref mod: language - changed
+%However, focussing on the description of diffusion processes the situation is changing completely.
+However, focussing on the description of diffusion processes the situation has changed completely.
+Qualitative and quantitative differences exist.
+First of all, a different pathway is suggested as the lowest energy path, which again might be attributed to the absence of quantum-mechanical effects in the classical interaction model.
+Secondly, the activation energy is overestimated by a factor of 2.4 compared to the more accurate quantum-mechanical methods and experimental findings.
+This is attributed to the sharp cut-off of the short range potential.
+As already pointed out in a previous study\cite{mattoni2007} the short cut-off is responsible for overestimated and unphysical high forces of next neighbor atoms.
+The overestimated migration barrier, however, affects the diffusion behavior of the C interstitials.
+By this artifact the mobility of the C atoms is tremendously decreased resulting in an inaccurate description or even absence of the dumbbell agglomeration as proposed by the precipitation model.
+
+\section{Summary}
+
+To conclude, we have shown that ab initio calculations on interstitial carbon in silicon are very close to the results expected from experimental data.
+The calculations presented in this work agree well with other theoretical results.
+So far, the best quantitative agreement with experimental findings has been achieved concerning the interstitial carbon mobility.
+For the first time, we have shown that the bond-centered configuration indeed constitutes a real local minimum configuration resulting in a net magnetization if spin polarized calculations are performed.
+Classical potentials, however, fail to describe the selected processes.
+This has been shown to have two reasons, i.e. the overestimated barrier of migration due to the artificial interaction cut-off on the one hand, and on the other hand the lack of quantum-mechanical effects which are crucial in the problem under study.
+% ref mod: language - being investigated
+%In order to get more insight on the SiC precipitation mechanism, further ab initio calculations are currently investigated.
+In order to get more insight on the SiC precipitation mechanism, further ab initio calculations are currently being performed.
+
+% ----------------------------------------------------
+\section*{Acknowledgment}
+We gratefully acknowledge financial support by the Bayerische Forschungsstiftung (DPA-61/05) and the Deutsche Forschungsgemeinschaft (DFG SCHM 1361/11).
+
+% --------------------------------- references -------------------
+
+%\bibliography{../../bibdb/bibdb}{}
+%\bibliographystyle{h-physrev3}
+
+\begin{thebibliography}{10}
+
+\bibitem{wesch96}
+W.~Wesch,
+\newblock Nucl. Instr. and Meth. B {\bf 116}, 305 (1996).
+
+\bibitem{morkoc94}
+H.~Morko\c{c}, S.~Strite, G.~B.~Gao, M.~E.~Lin, B.~Sverdlov, and M.~Burns,
+\newblock J. Appl. Phys. {\bf 76}, 1363 (1994).
+
+\bibitem{capano97}
+M.~A. Capano and R.~J. Trew,
+\newblock MRS Bull. {\bf 22}, 19 (1997).
+
+\bibitem{fischer90}
+G.~R. Fisher and P.~Barnes,
+\newblock Philosophical Magazine Part B {\bf 61}, 217 (1990).
+
+\bibitem{powell90}
+J.~A. Powell, D.~J.~Larkin, L.~G.~Matus, W.~J.~Choyke, J.~L.~Bradshaw, L.~Henderson, M.~Yoganathan, J.~Yang, and P.~Pirouz,
+\newblock Appl. Phys. Lett. {\bf 56}, 1353 (1990).
+
+\bibitem{fissel95}
+A.~Fissel, U.~Kaiser, E.~Ducke, B.~Schr{\"{o}}ter, and W.~Richter,
+\newblock J. Cryst. Growth {\bf 154}, 72 (1995).
+
+\bibitem{fissel95_apl}
+A.~Fissel, B.~Schr{\"{o}}ter, and W.~Richter,
+\newblock Appl. Phys. Lett. {\bf 66}, 3182 (1995).
+
+\bibitem{nishino83}
+S.~Nishino, J.~A. Powell, and H.~A. Will,
+\newblock Appl. Phys. Lett. {\bf 42}, 460 (1983).
+
+\bibitem{nishino87}
+S.~Nishino, H.~Suhara, H.~Ono, and H.~Matsunami,
+\newblock J. Appl. Phys. {\bf 61}, 4889 (1987).
+
+\bibitem{kitabatake93}
+M.~Kitabatake, M.~Deguchi, and T.~Hirao,
+\newblock J. Appl. Phys. {\bf 74}, 4438 (1993).
+
+\bibitem{borders71}
+J.~A. Borders, S.~T. Picraux, and W.~Beezhold,
+\newblock Appl. Phys. Lett. {\bf 18}, 509 (1971).
+
+\bibitem{lindner99}
+J.~K.~N. Lindner and B.~Stritzker,
+\newblock Nucl. Instr. and Meth. B {\bf 147}, 249 (1999).
+
+\bibitem{lindner01}
+J.~K.~N. Lindner,
+\newblock Nucl. Instr. and Meth. B {\bf 178}, 44 (2001).
+
+\bibitem{lindner02}
+J.~K.~N. Lindner,
+\newblock Appl. Phys. A {\bf 77}, 27 (2003).
+
+\bibitem{werner96}
+P.~Werner, R.~K{\"{o}}gler, W.~Skorupa, and D.~Eichler,
+\newblock {TEM} investigation of {C}-si defects in carbon implanted silicon,
+\newblock in {\em Ion Implantation Technology. Proceedings of the 11th
+ International Conference on}, pp. 675--678, 1996.
+
+\bibitem{werner97}
+P.~Werner, S.~Eichler, G.~Mariani, R.~K{\"{o}}gler, and W.~Skorupa,
+\newblock Appl. Phys. Lett. {\bf 70}, 252 (1997).
+
+\bibitem{lindner99_2}
+J.~K.~N. Lindner and B.~Stritzker,
+\newblock Nucl. Instr. and Meth. B {\bf 148}, 528 (1999).
+
+\bibitem{strane96}
+J.~W. Strane, S.~R.~Lee, H.~J.~Stein, S.~T.~Picraux, J.~K.~Watanabe, and J.~W.~Mayer,
+\newblock J. Appl. Phys. {\bf 79}, 637 (1996).
+
+\bibitem{laveant2002}
+P.~Lav\'eant, G.~Gerth, P.~Werner, and U.~G{\"{o}}sele,
+\newblock Mater. Sci. Eng. B {\bf 89}, 241 (2002).
+
+\bibitem{stillinger85}
+F.~H. Stillinger and T.~A. Weber,
+\newblock Phys. Rev. B {\bf 31}, 5262 (1985).
+
+\bibitem{brenner90}
+D.~W. Brenner,
+\newblock Phys. Rev. B {\bf 42}, 9458 (1990).
+
+\bibitem{tersoff_si3}
+J.~Tersoff,
+\newblock Phys. Rev. B {\bf 38}, 9902 (1988).
+
+\bibitem{bazant96}
+M.~Z. Bazant and E.~Kaxiras,
+\newblock Phys. Rev. Lett. {\bf 77}, 4370 (1996).
+
+\bibitem{bazant97}
+M.~Z. Bazant, E.~Kaxiras, and J.~F. Justo,
+\newblock Phys. Rev. B {\bf 56}, 8542 (1997).
+
+\bibitem{justo98}
+J.~F. Justo, M.~Z. Bazant, E.~Kaxiras, V.~V. Bulatov, and S.~Yip,
+\newblock Phys. Rev. B {\bf 58}, 2539 (1998).
+
+\bibitem{lucas10}
+G.~Lucas, M.~Bertolus, and L.~Pizzagalli,
+\newblock J. Phys.: Condens. Matter {\bf 22}, 035802 (2010).
+
+\bibitem{tersoff_m}
+J.~Tersoff,
+\newblock Phys. Rev. B {\bf 39}, 5566 (1989).
+
+\bibitem{gao02}
+F.~Gao and W.~J. Weber,
+\newblock Nucl. Instr. and Meth. B {\bf 191}, 504 (2002).
+
+\bibitem{godet03}
+J.~Godet, L.~Pizzagalli, S.~Brochard, and P.~Beauchamp,
+\newblock J. Phys.: Condens. Matter {\bf 15}, 6943 (2003).
+
+\bibitem{holmstroem08}
+E.~Holmstr{\"o}m, A.~Kuronen, and K.~Nordlund,
+\newblock Phys. Rev. B {\bf 78}, 045202 (2008).
+
+\bibitem{mattoni2002}
+A.~Mattoni, F.~Bernardini, and L.~Colombo,
+\newblock Phys. Rev. B {\bf 66}, 195214 (2002).
+
+\bibitem{posselt06}
+M.~Posselt, F.~Gao, and W.~J. Weber,
+\newblock Phys. Rev. B {\bf 73}, 125206 (2006).
+
+\bibitem{gao04}
+F.~Gao, W.~J. Weber, M.~Posselt, and V.~Belko,
+\newblock Phys. Rev. B {\bf 69}, 245205 (2004).
+
+\bibitem{rauls03a}
+E.~Rauls, T.~Frauenheim, A.~Gali, and P.~De\'ak,
+\newblock Phys. Rev. B {\bf 68}, 155208 (2003).
+
+\bibitem{gao07}
+F.~Gao, J.~Du, E.~J. Bylaska, M.~Posselt, and W.~J. Weber,
+\newblock Appl. Phys. Lett. {\bf 90}, 221915 (2007).
+
+\bibitem{gao2001}
+F.~Gao, E.~J. Bylaska, W.~J. Weber, and L.~R. Corrales,
+\newblock Phys. Rev. B {\bf 64}, 245208 (2001).
+
+\bibitem{bockstedte03}
+M.~Bockstedte, A.~Mattausch, and O.~Pankratov,
+\newblock Phys. Rev. B {\bf 68}, 205201 (2003).
+
+\bibitem{gali03}
+A.~Gali, P.~De\'ak, E.~Rauls, N.~T.~Son, I.~G.~Ivanov, F.~H.~C.~Carlsson, E.~Janz\'en, and W.~J.~Choyke,
+\newblock Phys. Rev. B {\bf 67}, 155203 (2003).
+
+\bibitem{chen98}
+J.~Chen, P.~Jung, and H.~Klein,
+\newblock J. Nucl. Mater. {\bf 258-263}, 1803 (1998).
+
+\bibitem{weber01}
+W.~J. Weber, W.~Jiang, and S.~Thevuthasan,
+\newblock Nucl. Instr. and Meth. B {\bf 175-177}, 26 (2001).
+
+\bibitem{balamane92}
+H.~Balamane, T.~Halicioglu, and W.~A. Tiller,
+\newblock Phys. Rev. B {\bf 46}, 2250 (1992).
+
+\bibitem{huang95}
+H.~Huang, N.~M. Ghoniem, J.~K. Wong, and M.~Baskes,
+\newblock Modelling and Simulation in Mater. Sci. Eng. {\bf
+ 3}, 615 (1995).
+
+\bibitem{tang95}
+M.~Tang and S.~Yip,
+\newblock Phys. Rev. B {\bf 52}, 15150 (1995).
+
+\bibitem{gao02a}
+F.~Gao and W.~J. Weber,
+\newblock Phys. Rev. B {\bf 66}, 024106 (2002).
+
+\bibitem{mattoni2007}
+A.~{Mattoni}, M.~{Ippolito}, and L.~{Colombo},
+\newblock Phys. Rev. B {\bf 76}, 224103 (2007).
+
+\bibitem{devanathan98_2}
+R.~Devanathan, T.~D. de~la Rubia, and W.~J. Weber,
+\newblock J. Nucl. Mater. {\bf 253}, 47 (1998).
+
+\bibitem{nordlund97}
+K.~Nordlund, N.~Runeberg, and D.~Sundholm,
+\newblock Nucl. Instr. and Meth. B {\bf 132}, 45 (1997).
+
+\bibitem{albe_sic_pot}
+P.~Erhart and K.~Albe,
+\newblock Phys. Rev. B {\bf 71}, 035211 (2005).
+
+\bibitem{kresse96}
+G.~Kresse and J.~Furthm{\"{u}}ller,
+\newblock Comp. Mater. Sci. {\bf 6}, 15 (1996).
+
+\bibitem{perdew86}
+J.~P. Perdew and W.~Yue,
+\newblock Phys. Rev. B {\bf 33}, 8800 (1986).
+
+\bibitem{perdew92}
+J.~P. Perdew, J.~A.~Chevary, S.~H.~Vosko, K.~A.~Jackson, M.~R.~Pederson, D.~J.~Singh, and C.~Fiolhais,
+\newblock Phys. Rev. B {\bf 46}, 6671 (1992).
+
+\bibitem{hamann79}
+D.~R. Hamann, M.~Schl{\"u}ter, and C.~Chiang,
+\newblock Phys. Rev. Lett. {\bf 43}, 1494 (1979).
+
+\bibitem{vanderbilt90}
+D.~Vanderbilt,
+\newblock Phys. Rev. B {\bf 41}, 7892 (1990).
+
+\bibitem{berendsen84}
+H.~J.~C. Berendsen, J.~P.~M. Postma, W.~F. van Gunsteren, A.~DiNola, and J.~R.
+ Haak,
+\newblock J. Chem. Phys. {\bf 81}, 3684 (1984).
+
+\bibitem{verlet67}
+L.~Verlet,
+\newblock Phys. Rev. {\bf 159}, 98 (1967).
+
+\bibitem{tersoff90}
+J.~Tersoff,
+\newblock Phys. Rev. Lett. {\bf 64}, 1757 (1990).
+
+\bibitem{dal_pino93}
+A.~{Dal Pino}, A.~M. Rappe, and J.~D. Joannopoulos,
+\newblock Phys. Rev. B {\bf 47}, 12554 (1993).
+
+\bibitem{capaz94}
+R.~B. Capaz, A.~{Dal Pino}, and J.~D. Joannopoulos,
+\newblock Phys. Rev. B {\bf 50}, 7439 (1994).
+
+\bibitem{burnard93}
+M.~J. Burnard and G.~G. DeLeo,
+\newblock Phys. Rev. B {\bf 47}, 10217 (1993).
+
+\bibitem{leary97}
+P.~Leary, R.~Jones, S.~{\"O}berg, and V.~J.~B. Torres,
+\newblock Phys. Rev. B {\bf 55}, 2188 (1997).
+
+\bibitem{watkins76}
+G.~D. Watkins and K.~L. Brower,
+\newblock Phys. Rev. Lett. {\bf 36}, 1329 (1976).
+
+\bibitem{song90}
+L.~W. Song and G.~D. Watkins,
+\newblock Phys. Rev. B {\bf 42}, 5759 (1990).
+
+\bibitem{kaukonen98}
+M.~Kaukonen, P.~K.~Sitch, G.~Jungnickel, R.~M.~Nieminen, S.~P{\"o}ykk{\"o}, D.~Porezag, and Th.~Frauenheim,
+\newblock Phys. Rev. B {\bf 57}, 9965 (1998).
+
+\bibitem{lindner06}
+J.~K.~N.~Lindner, M.~H{\"a}berlen, G.~Thorwarth, and B.~Stritzker,
+\newblock Mater. Sci. Eng. C {\bf 26}, 857 (2006).
+
+\bibitem{tipping87}
+A.~K. Tipping and R.~C. Newman,
+\newblock Semicond. Sci. Technol. {\bf 2}, 315 (1987).
+
+\end{thebibliography}
+
+\end{document}
projects / lectures / latex.git / blobdiff