\section{Silicon self-interstitials}
For investigating the \si{} structures a Si atom is inserted or removed according to Fig.~\ref{fig:basics:ins_pos} of section \ref{section:basics:defects}.
-The formation energies of \si{} configurations are listed in Table \ref{tab:defects:si_self} for both methods used in this work as well as results obtained by other {\em ab initio} studies \cite{al-mushadani03,leung99}.
+The formation energies of \si{} configurations are listed in Table \ref{tab:defects:si_self} for both methods used in this work as well as results obtained by other {\em ab initio} studies~\cite{al-mushadani03,leung99}.
\bibpunct{}{}{,}{n}{}{}
\begin{table}[tp]
\begin{center}
\textsc{vasp} & 3.39 & 3.42 & 3.77 & 4.41 & 3.63 \\
\textsc{posic} & 4.39 & 4.48$^*$ & 3.40 & 5.42 & 3.13 \\
\multicolumn{6}{c}{Other {\em ab initio} studies} \\
-Ref. \cite{al-mushadani03} & 3.40 & 3.45 & - & - & 3.53 \\
-Ref. \cite{leung99} & 3.31 & 3.31 & 3.43 & - & - \\
+Ref.~\cite{al-mushadani03} & 3.40 & 3.45 & - & - & 3.53 \\
+Ref.~\cite{leung99} & 3.31 & 3.31 & 3.43 & - & - \\
\hline
\hline
\end{tabular}
This is nicely reproduced by the DFT calculations performed in this work.
It has turned out to be very difficult to capture the results of quantum-mechanical calculations in analytical potential models.
-Among the established analytical potentials only the environment-dependent interatomic potential (EDIP) \cite{bazant97,justo98} and Stillinger-Weber \cite{stillinger85} potential reproduce the correct order in energy of the defects.
+Among the established analytical potentials only the environment-dependent interatomic potential (EDIP)~\cite{bazant97,justo98} and Stillinger-Weber~\cite{stillinger85} potential reproduce the correct order in energy of the defects.
However, these potentials show shortcomings concerning the description of other physical properties and are unable to describe the C-C and C-Si interaction.
In fact the EA potential calculations favor the tetrahedral defect configuration.
This limitation is assumed to arise due to the cut-off.
In the tetrahedral configuration the second neighbors are only slightly more distant than the first neighbors, which creates the particular problem.
-Indeed, an increase of the cut-off results in increased values of the formation energies \cite{albe_sic_pot}, which is most significant for the tetrahedral configuration.
-The same issue has already been discussed by Tersoff \cite{tersoff90} with regard to the description of the tetrahedral C defect using his potential.
+Indeed, an increase of the cut-off results in increased values of the formation energies~\cite{albe_sic_pot}, which is most significant for the tetrahedral configuration.
+The same issue has already been discussed by Tersoff~\cite{tersoff90} with regard to the description of the tetrahedral C defect using his potential.
While not completely rendering impossible further, more challenging empirical potential studies on large systems, the artifact has to be taken into account in the investigations of defect combinations later on in this chapter.
-The hexagonal configuration is not stable opposed to results of the authors of the potential \cite{albe_sic_pot}.
+The hexagonal configuration is not stable opposed to results of the authors of the potential~\cite{albe_sic_pot}.
In the first two picoseconds, while kinetic energy is decoupled from the system, the \si{} seems to condense at the hexagonal site.
The formation energy of \unit[4.48]{eV} is determined by this low kinetic energy configuration shortly before the relaxation process starts.
The \si{} atom then begins to slowly move towards an energetically more favorable position very close to the tetrahedral one but slightly displaced along the three coordinate axes.
-The formation energy of \unit[3.96]{eV} for this type of interstitial is equal to the result for the hexagonal one in the original work \cite{albe_sic_pot}.
+The formation energy of \unit[3.96]{eV} for this type of interstitial is equal to the result for the hexagonal one in the original work~\cite{albe_sic_pot}.
Obviously, the authors did not carefully check the relaxed results assuming a hexagonal configuration.
In Fig.~\ref{fig:defects:kin_si_hex} the relaxation process is shown on the basis of the kinetic energy plot.
\begin{figure}[tp]
\caption{Kinetic energy plot of the relaxation process of the hexagonal silicon self-interstitial defect simulation using the EA potential.}
\label{fig:defects:kin_si_hex}
\end{figure}
-To exclude failures in the implementation of the potential or the MD code itself the hexagonal defect structure was double-checked with the \textsc{parcas} MD code \cite{parcas_md}.
+To exclude failures in the implementation of the potential or the MD code itself the hexagonal defect structure was double-checked with the \textsc{parcas} MD code~\cite{parcas_md}.
The respective relaxation energetics are likewise plotted and look similar to the energetics obtained by \textsc{posic}.
In fact, the same type of interstitial arises using random insertions.
In addition, variations exist, in which the displacement is only along two \hkl<1 0 0> axes ($E_\text{f}=3.8\,\text{eV}$) or along a single \hkl<1 0 0> axes ($E_\text{f}=3.6\,\text{eV}$) successively approximating the tetrahedral configuration and formation energy.
The relaxed configurations are visualized in Fig.~\ref{fig:defects:c_conf}.
Again, the displayed structures are the results obtained by the classical potential calculations.
The type of reservoir of the C impurity to determine the formation energy of the defect is chosen to be SiC.
-This is consistent with the methods used in the articles \cite{tersoff90,dal_pino93}, which the results are compared to in the following.
+This is consistent with the methods used in the articles~\cite{tersoff90,dal_pino93}, which the results are compared to in the following.
Hence, the chemical potential of Si and C is determined by the cohesive energy of Si and SiC as discussed in section \ref{section:basics:defects}.
\begin{table}[tp]
\begin{center}
\textsc{posic} & 6.09 & 9.05$^*$ & 3.88 & 5.18 & 0.75 & 5.59$^*$ \\
\textsc{vasp} & Unstable & Unstable & 3.72 & 4.16 & 1.95 & 4.66 \\
Other studies & & & & & & \\
- Tersoff \cite{tersoff90} & 3.8 & 6.7 & 4.6 & 5.9 & 1.6 & 5.3 \\
- {\em Ab initio} \cite{dal_pino93,capaz94} & - & - & x & - & 1.89 & x+2.1 \\
+ Tersoff~\cite{tersoff90} & 3.8 & 6.7 & 4.6 & 5.9 & 1.6 & 5.3 \\
+ {\em Ab initio}~\cite{dal_pino93,capaz94} & - & - & x & - & 1.89 & x+2.1 \\
\hline
\hline
\end{tabular}
\end{figure}
\cs{} occupying an already vacant Si lattice site, which is in fact not an interstitial defect, is found to be the lowest configuration in energy for all potential models.
-An experimental value of the formation energy of \cs{} was determined by a fit to solubility data yielding a concentration of $3.5 \times 10^{24} \exp{(-2.3\,\text{eV}/k_{\text{B}}T)} \text{ cm}^{-3}$ \cite{bean71}.
+An experimental value of the formation energy of \cs{} was determined by a fit to solubility data yielding a concentration of $3.5 \times 10^{24} \exp{(-2.3\,\text{eV}/k_{\text{B}}T)} \text{ cm}^{-3}$~\cite{bean71}.
However, there is no particular reason for treating the prefactor as a free parameter in the fit to the experimental data.
It is simply given by the atomic density of pure silicon, which is $5\times 10^{22}\text{ cm}^{-3}$.
-Tersoff \cite{tersoff90} and Dal Pino et al. \cite{dal_pino93} pointed out that by combining this prefactor with the calculated values for the energy of formation ranging from \unit[1.6--1.89]{eV} an excellent agreement with the experimental solubility data within the entire temperature range of the experiment is obtained.
+Tersoff~\cite{tersoff90} and Dal Pino et al.~\cite{dal_pino93} pointed out that by combining this prefactor with the calculated values for the energy of formation ranging from \unit[1.6--1.89]{eV} an excellent agreement with the experimental solubility data within the entire temperature range of the experiment is obtained.
This reinterpretation of the solubility data, first proposed by Tersoff and later on reinforced by Dal~Pino~et~al. is in good agreement with the results of the quantum-mechanical calculations performed in this work.
Unfortunately the EA potential undervalues the formation energy roughly by a factor of two, which is a definite drawback of the potential.
Except for Tersoff's results for the tetrahedral configuration, the \ci{} \hkl<1 0 0> DB is the energetically most favorable interstitial configuration.
-As mentioned above, the low energy of formation for the tetrahedral interstitial in the case of the Tersoff potential is believed to be an artifact of the abrupt cut-off set to \unit[2.5]{\AA} (see Ref. 11 and 13 in \cite{tersoff90}) and the real formation energy is, thus, supposed to be located between \unit[3--10]{eV}.
+As mentioned above, the low energy of formation for the tetrahedral interstitial in the case of the Tersoff potential is believed to be an artifact of the abrupt cut-off set to \unit[2.5]{\AA} (see Ref. 11 and 13 in~\cite{tersoff90}) and the real formation energy is, thus, supposed to be located between \unit[3--10]{eV}.
Keeping these considerations in mind, the \ci{} \hkl<1 0 0> DB is the most favorable interstitial configuration for all interaction models.
This finding is in agreement with several theoretical~\cite{burnard93,leary97,dal_pino93,capaz94,jones04} and experimental~\cite{watkins76,song90} investigations, which all predict this configuration to be the ground state.
However, no energy of formation for this type of defect based on first-principles calculations has yet been explicitly stated in literature.
The highest energy is observed for the hexagonal interstitial configuration using classical potentials.
Quantum-mechanical calculations reveal this configuration to be unstable, which is also reproduced by the EA potential.
In both cases a relaxation towards the \ci{} \hkl<1 0 0> DB configuration is observed.
-Opposed to results of the first-principles calculations, Tersoff finds this configuration to be stable \cite{tersoff90}.
+Opposed to results of the first-principles calculations, Tersoff finds this configuration to be stable~\cite{tersoff90}.
In fact, the stability of the hexagonal interstitial could not be reproduced in simulations performed in this work using the unmodified Tersoff potential parameters.
Unfortunately, apart from the modified parameters, no more conditions specifying the relaxation process are given in Tersoff's study on C point defects in Si.
\label{subsection:100db}
As the \ci{} \hkl<1 0 0> DB constitutes the ground-state configuration of a C atom incorporated into otherwise perfect c-Si it is the most probable and, hence, one of the most important interstitial configurations of C in Si.
-The structure was initially suspected by IR local vibrational mode absorption \cite{bean70} and finally verified by electron paramagnetic resonance (EPR) \cite{watkins76} studies on irradiated Si substrates at low temperatures.
+The structure was initially suspected by IR local vibrational mode absorption~\cite{bean70} and finally verified by electron paramagnetic resonance (EPR)~\cite{watkins76} studies on irradiated Si substrates at low temperatures.
Fig.~\ref{fig:defects:100db_cmp} schematically shows the \ci{} \hkl<1 0 0> DB structure and Table \ref{tab:defects:100db_cmp} lists the details of the atomic displacements, distances and bond angles obtained by classical potential and quantum-mechanical calculations.
For comparison, the obtained structures for both methods are visualized in Fig.~\ref{fig:defects:100db_vis_cmp}.
\label{img:defects:bc_conf}
\end{figure}
In the BC interstitial configuration the interstitial atom is located in between two next neighbored Si atoms forming linear bonds.
-In a previous study this configuration was found to constitute an intermediate saddle point configuration determining the migration barrier of one possible migration path of a \ci{} \hkl<1 0 0> DB configuration into an equivalent one \cite{capaz94}.
+In a previous study this configuration was found to constitute an intermediate saddle point configuration determining the migration barrier of one possible migration path of a \ci{} \hkl<1 0 0> DB configuration into an equivalent one~\cite{capaz94}.
This is in agreement with results of the EA potential simulations, which reveal this configuration to be unstable relaxing into the \ci{} \hkl<1 1 0> configuration.
However, this fact could not be reproduced by spin polarized \textsc{vasp} calculations performed in this work.
Present results suggest this configuration to correspond to a real local minimum.
\label{fig:defects:00-1_0-10_mig}
\end{figure}
Fig.~\ref{fig:defects:00-1_0-10_mig} shows the migration barrier and structures of the \ci{} \hkl<0 0 -1> to \hkl<0 -1 0> DB transition.
-The resulting migration barrier of approximately \unit[0.9]{eV} is very close to the experimentally obtained values of \unit[0.70]{eV} \cite{lindner06}, \unit[0.73]{eV} \cite{song90} and \unit[0.87]{eV} \cite{tipping87}.
+The resulting migration barrier of approximately \unit[0.9]{eV} is very close to the experimentally obtained values of \unit[0.70]{eV}~\cite{lindner06}, \unit[0.73]{eV}~\cite{song90} and \unit[0.87]{eV}~\cite{tipping87}.
\begin{figure}[tp]
\begin{center}
\end{figure}
The third migration path, in which the DB is changing its orientation, is shown in Fig.~\ref{fig:defects:00-1_0-10_ip_mig}.
An energy barrier of roughly \unit[1.2]{eV} is observed.
-Experimentally measured activation energies for reorientation range from \unit[0.77]{eV} to \unit[0.88]{eV} \cite{watkins76,song90}.
+Experimentally measured activation energies for reorientation range from \unit[0.77]{eV} to \unit[0.88]{eV}~\cite{watkins76,song90}.
Thus, this pathway is more likely to be composed of two consecutive steps of the second path.
Since the activation energy of the first and last migration path is much greater than the experimental value, the second path is identified to be responsible as a migration path for the most likely C interstitial in Si explaining both, annealing and reorientation experiments.
The activation energy of roughly \unit[0.9]{eV} nicely compares to experimental values reinforcing the correct identification of the C-Si DB diffusion mechanism.
Slightly increased values compared to experiment might be due to the tightened constraints applied in the modified CRT approach.
-Nevertheless, the theoretical description performed in this work is improved compared to a former study \cite{capaz94}, which underestimates the experimental value by \unit[35]{\%}.
+Nevertheless, the theoretical description performed in this work is improved compared to a former study~\cite{capaz94}, which underestimates the experimental value by \unit[35]{\%}.
In addition, it is finally shown that the BC configuration, for which spin polarized calculations are necessary, constitutes a real local minimum instead of a saddle point configuration due to the presence of restoring forces for displacements in migration direction.
\begin{figure}[tp]
On the other hand, the activation energy obtained by classical potential simulations is tremendously overestimated by a factor of 2.4 to 3.5.
The overestimated barrier is due to the short range character of the potential, which drops the interaction to zero within the first and next neighbor distance.
Since the total binding energy is accommodated within a short distance, which according to the universal energy relation would usually correspond to a much larger distance, unphysical high forces between two neighbored atoms arise.
-This is explained in more detail in a previous study \cite{mattoni2007}.
+This is explained in more detail in a previous study~\cite{mattoni2007}.
Thus, atomic diffusion is wrongly described in the classical potential approach.
The probability of already rare diffusion events is further decreased for this reason.
However, agglomeration of C and diffusion of Si self-interstitials are an important part of the proposed SiC precipitation mechanism.
\caption[Relaxed structures of defect combinations obtained by creating {\hkl[1 0 0]} and {\hkl[0 -1 0]} DBs at position 1.]{Relaxed structures of defect combinations obtained by creating \hkl[1 0 0] (a) and \hkl[0 -1 0] (b) DBs at position 1.}
\label{fig:defects:comb_db_01}
\end{figure}
-Mattoni~et~al. \cite{mattoni2002} predict the ground-state configuration of \ci{} \hkl<1 0 0>-type defect pairs for a \hkl[1 0 0] or equivalently a \hkl[0 1 0] defect created at position 1 with both defects basically maintaining the as-isolated DB structure, resulting in a binding energy of \unit[-2.1]{eV}.
+Mattoni~et~al.~\cite{mattoni2002} predict the ground-state configuration of \ci{} \hkl<1 0 0>-type defect pairs for a \hkl[1 0 0] or equivalently a \hkl[0 1 0] defect created at position 1 with both defects basically maintaining the as-isolated DB structure, resulting in a binding energy of \unit[-2.1]{eV}.
In the present study, a further relaxation of this defect structure is observed.
The C atom of the second and the Si atom of the initial DB move towards each other forming a bond, which results in a somewhat lower binding energy of \unit[-2.25]{eV}.
The corresponding defect structure is displayed in Fig.~\ref{fig:defects:225}.
\end{figure}
Configuration A consists of a C$_{\text{i}}$ \hkl[0 0 -1] DB with threefold coordinated Si and C DB atoms slightly disturbed by the C$_{\text{s}}$ at position 3, facing the Si DB atom as a neighbor.
By a single bond switch, i.e.\ the breaking of a Si-Si in favor of a Si-C bond, configuration B is obtained, which shows a twofold coordinated Si atom located in between two substitutional C atoms residing on regular Si lattice sites.
-This configuration has been identified and described by spectroscopic experimental techniques \cite{song90_2} as well as theoretical studies \cite{leary97,capaz98}.
+This configuration has been identified and described by spectroscopic experimental techniques~\cite{song90_2} as well as theoretical studies~\cite{leary97,capaz98}.
Configuration B is found to constitute the energetically slightly more favorable configuration.
However, the gain in energy due to the significantly lower energy of a Si-C compared to a Si-Si bond turns out to be smaller than expected due to a large compensation by introduced strain as a result of the Si interstitial structure.
-Present results show a difference in energy of states A and B, which exactly matches the experimental value of \unit[0.02]{eV} \cite{song90_2}, reinforcing qualitatively correct results of previous theoretical studies on these structures.
+Present results show a difference in energy of states A and B, which exactly matches the experimental value of \unit[0.02]{eV}~\cite{song90_2}, reinforcing qualitatively correct results of previous theoretical studies on these structures.
% mattoni: A favored by 0.4 eV - NO, it is indeed B (reinforce Song and Capaz)!
%
% AB transition
-The migration barrier is identified to be \unit[0.44]{eV}, almost three times higher than the experimental value of \unit[0.16]{eV} \cite{song90_2} estimated for the neutral charge state transition in p- and n-type Si.
+The migration barrier is identified to be \unit[0.44]{eV}, almost three times higher than the experimental value of \unit[0.16]{eV}~\cite{song90_2} estimated for the neutral charge state transition in p- and n-type Si.
Keeping in mind the formidable agreement of the energy difference with experiment, the overestimated activation energy is quite unexpected.
Obviously, either the CRT algorithm fails to seize the actual saddle point structure or the influence of dopants has exceptional effect in the experimentally covered diffusion process being responsible for the low migration barrier.
% not satisfactory!
\caption[Migration barrier of the \si{} {\hkl[1 1 0]} DB into the hexagonal and tetrahedral configuration as well as the hexagonal \si{} to tetrahedral \si{} transition.]{Migration barrier of the \si{} \hkl[1 1 0] DB into the hexagonal (H) and tetrahedral (T) configuration as well as the hexagonal \si{} to tetrahedral \si{} transition.}
\label{fig:defects:si_mig2}
\end{figure}
-The obtained activation energies are of the same order of magnitude than values derived from other {\em ab initio} studies \cite{bloechl93,sahli05}.
+The obtained activation energies are of the same order of magnitude than values derived from other {\em ab initio} studies~\cite{bloechl93,sahli05}.
The low barriers indeed enable configurations of further separated \cs{} and \si{} atoms by the highly mobile \si{} atom departing from the \cs{} defect as observed in the previously discussed MD simulation.
% kept for nostalgical reason!
\ifnum1=0
-Obtained results for separated point defects in Si are in good agreement to previous theoretical work on this subject, both for intrinsic defects \cite{leung99,al-mushadani03} as well as for C point defects \cite{dal_pino93,capaz94}.
-The ground-state configurations of these defects, i.e.\ the Si$_{\text{i}}$ \hkl<1 1 0> and C$_{\text{i}}$ \hkl<1 0 0> DB, are reproduced and compare well to previous findings of theoretical investigations on Si$_{\text{i}}$ \cite{leung99,al-mushadani03} as well as theoretical \cite{dal_pino93,capaz94,burnard93,leary97,jones04} and experimental \cite{watkins76,song90} studies on C$_{\text{i}}$.
-A quantitatively improved activation energy of \unit[0.9]{eV} for a qualitatively equal migration path based on studies by Capaz et.~al.~\cite{capaz94} to experimental values \cite{song90,lindner06,tipping87} ranging from \unit[0.70--0.87]{eV} reinforce their derived mechanism of diffusion for C$_{\text{i}}$ in Si
+Obtained results for separated point defects in Si are in good agreement to previous theoretical work on this subject, both for intrinsic defects~\cite{leung99,al-mushadani03} as well as for C point defects~\cite{dal_pino93,capaz94}.
+The ground-state configurations of these defects, i.e.\ the Si$_{\text{i}}$ \hkl<1 1 0> and C$_{\text{i}}$ \hkl<1 0 0> DB, are reproduced and compare well to previous findings of theoretical investigations on Si$_{\text{i}}$~\cite{leung99,al-mushadani03} as well as theoretical~\cite{dal_pino93,capaz94,burnard93,leary97,jones04} and experimental~\cite{watkins76,song90} studies on C$_{\text{i}}$.
+A quantitatively improved activation energy of \unit[0.9]{eV} for a qualitatively equal migration path based on studies by Capaz et.~al.~\cite{capaz94} to experimental values~\cite{song90,lindner06,tipping87} ranging from \unit[0.70--0.87]{eV} reinforce their derived mechanism of diffusion for C$_{\text{i}}$ in Si
However, it turns out that the BC 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.
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 to 3.5 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.
+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 DB agglomeration as proposed by one of the precipitation models.
Quantum-mechanical investigations of two \ci{} of the \hkl<1 0 0>-type and varying separations and orientations state an attractive interaction between these interstitials.
-Obtained results for the most part compare well with results gained in previous studies \cite{leary97,capaz98,mattoni2002,liu02} and show an astonishingly good agreement with experiment \cite{song90}.
+Obtained results for the most part compare well with results gained in previous studies~\cite{leary97,capaz98,mattoni2002,liu02} and show an astonishingly good agreement with experiment~\cite{song90}.
%
Depending on orientation, energetically favorable configurations are found, in which these two interstitials are located close together instead of the occurrence of largely separated and isolated defects.
This is due to strain compensation enabled by the combination of such defects in certain orientations.
To conclude, the available results suggest precipitation by successive agglomeration of C$_{\text{s}}$.
However, the agglomeration and rearrangement of C$_{\text{s}}$ is only possible by mobile C$_{\text{i}}$, which has to be present at the same time.
Accordingly, the process is governed by both, C$_{\text{s}}$ accompanied by Si$_{\text{i}}$ as well as C$_{\text{i}}$.
-It is worth to mention that there is no contradiction to results of the HREM studies \cite{werner96,werner97,eichhorn99,lindner99_2,koegler03}.
+It is worth to mention that there is no contradiction to results of the HREM studies~\cite{werner96,werner97,eichhorn99,lindner99_2,koegler03}.
Regions showing dark contrasts in an otherwise undisturbed Si lattice are attributed to C atoms in the interstitial lattice.
However, there is no particular reason for the C species to reside in the interstitial lattice.
Contrasts are also assumed for Si$_{\text{i}}$.
projects / lectures / latex.git / blobdiff