One possibility is to simply skip the force contributions containing the derivatives of the cut-off function, which was successfully applied to reproduce the brittle propagation of fracture in SiC at zero temperature~\cite{mattoni2007}.
Another one is to use variable cut-off values scaled by the system volume, which properly describes thermomechanical properties of 3C-SiC~\cite{tang95} but might be rather ineffective for the challenge inherent to this study.
-To conclude the obstacle needed to get passed is twofold.
+To conclude, the obstacle needed to get passed is twofold.
The sharp cut-off of the employed bond order model potential introduces overestimated high forces between next neighbored atoms enhancing the problem of slow phase space propagation immanent to MD simulations.
This obstacle could be referred to as {\em potential enhanced slow phase space propagation}.
Due to this, pushing the time scale to the limits of computational resources or applying one of the above mentioned accelerated dynamics methods exclusively will not be sufficient enough.
\section{Long time scale simulations at maximum temperature}
-As discussed in section~\ref{section:md:limit} and~\ref{section: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.
+As discussed in section~\ref{section:md:limit} and~\ref{section: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.
Furthermore, these results indicate that increased temperatures are necessary to drive the system out of equilibrium enabling conditions needed for the formation of a metastable cubic polytype of SiC.
A maximum temperature to avoid melting is determined in section~\ref{section:md:tval} to be 120 \% of the Si melting point but due to defects lowering the transition point a maximum temperature of 95 \% of the Si melting temperature is considered useful.
Next to the employment of longer time scales and a maximum temperature a few more changes are applied.
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 first MD simulations.
-To speed up the simulation the initial simulation volume is reduced to 21 Si unit cells in each direction and 5500 inserted C atoms in either the whole volume or in a sphere with a radius of 3 nm corresponding to the size of a precipitate consisting of 5500 C atoms.
+To speed up the simulation, the initial simulation volume is reduced to 21 Si unit cells in each direction and 5500 inserted C atoms in either the whole volume or in a sphere with a radius of 3 nm corresponding to the size of a precipitate consisting of 5500 C atoms.
The \unit[100]{ps} sequence after C insertion intended for structural evolution is exchanged by a \unit[10]{ns} sequence, which is hoped to result in the occurrence of infrequent processes and a subsequent phase transition.
The return to lower temperatures is considered separately.
projects / lectures / latex.git / blobdiff