On the one hand this sheds light on characteristic values like the radial distribution function or the total amount of free energy for such a configuration that is aimed to be reproduced by simulation.
On the other hand, assuming a correct alignment of the precipitate with the c-Si matrix, properties of such precipitates and the surrounding as well as the interface can be investiagted.
Furthermore these investigations might establish the prediction of conditions necessary for the simulation of the precipitation process.
On the one hand this sheds light on characteristic values like the radial distribution function or the total amount of free energy for such a configuration that is aimed to be reproduced by simulation.
On the other hand, assuming a correct alignment of the precipitate with the c-Si matrix, properties of such precipitates and the surrounding as well as the interface can be investiagted.
Furthermore these investigations might establish the prediction of conditions necessary for the simulation of the precipitation process.
To play safe the precipitate configuration at 100 \% of the Si melting temperature is chosen and cooled down to $20\,^{\circ}\mathrm{C}$ with a cooling rate of $1\,^{\circ}\mathrm{C}/\text{ps}$.
{\color{blue}TODO: Wait for results and then compare structure (PC) and interface energy, maybe a energetically more favorable configuration arises.}
{\color{red}TODO: Mention the fact, that the precipitate is stable for eleveated temperatures, even for temperatures where the Si matrix is melting.}
To play safe the precipitate configuration at 100 \% of the Si melting temperature is chosen and cooled down to $20\,^{\circ}\mathrm{C}$ with a cooling rate of $1\,^{\circ}\mathrm{C}/\text{ps}$.
{\color{blue}TODO: Wait for results and then compare structure (PC) and interface energy, maybe a energetically more favorable configuration arises.}
{\color{red}TODO: Mention the fact, that the precipitate is stable for eleveated temperatures, even for temperatures where the Si matrix is melting.}
As discussed in section \ref{subsection:md:limit} and \ref{subsection:md:inct} a further increase of the system temperature might help to overcome limitations of the short range potential and accelerate the dynamics involved in structural evolution.
A maximum temperature to avoid melting was determined in section \ref{subsection:md:tval}, which is 120 \% of the Si melting point.
In the following simulations the system volume, the amount of C atoms inserted and the shape of the insertion volume are modified from the values used in the first MD simulations to now match the conditions given in the simulations of the self-constructed precipitate configuration for reasons of comparability.
To quantify, the initial simulation volume now consists of 21 Si unit cells in each direction and 5500 C atoms are inserted in either the whole volume or in a sphere with a radius of 3 nm.
Since the investigated temperatures exceed the Si melting point the initial Si bulk material is heated up slowly by $1\,^{\circ}\mathrm{C}/\text{ps}$ starting from $1650\,^{\circ}\mathrm{C}$ before the C insertion sequence is started.
As discussed in section \ref{subsection:md:limit} and \ref{subsection:md:inct} a further increase of the system temperature might help to overcome limitations of the short range potential and accelerate the dynamics involved in structural evolution.
A maximum temperature to avoid melting was determined in section \ref{subsection:md:tval}, which is 120 \% of the Si melting point.
In the following simulations the system volume, the amount of C atoms inserted and the shape of the insertion volume are modified from the values used in the first MD simulations to now match the conditions given in the simulations of the self-constructed precipitate configuration for reasons of comparability.
To quantify, the initial simulation volume now consists of 21 Si unit cells in each direction and 5500 C atoms are inserted in either the whole volume or in a sphere with a radius of 3 nm.
Since the investigated temperatures exceed the Si melting point the initial Si bulk material is heated up slowly by $1\,^{\circ}\mathrm{C}/\text{ps}$ starting from $1650\,^{\circ}\mathrm{C}$ before the C insertion sequence is started.
+{\color{red}TODO: Melting occurs, show and explain it and that it's due to the defects created.}
+
+{\color{red}TODO: Due to melting, after insertion, simulation is continued NVE, so melting hopefully will not occur, before it will be cooled down later on.}
+
+{\color{red}TODO: In additions simulations at 95 \% of the Si melting temperature are started again for longer times.}