-for $R_{ij} < r_{ij} < S_{ij}$ and otherwise zero.
-The derivations of the repulsive and attractive part are:
-\begin{eqnarray}
-\partial_{x_n^i} f_R(r_{ij}) & = & - \lambda_{ij} A_{ij} \exp (-\lambda_{ij} r_{ij})\\
-\partial_{x_n^i} f_A(r_{ij}) & = & \mu_{ij} B_{ij} \exp (-\mu_{ij} r_{ij}) \textrm{ .}
-\end{eqnarray}
-The angle $\theta_{ijk}$ can be expressed by the atom distances with the law of cosines:
-\begin{eqnarray}
-\theta_{ijk} & = & \arccos \Big( (r_{ij}^2 + r_{ik}^2 - r_{jk}^2)/(2 r_{ij} r_{ik}) \Big) \\
-\partial_{x^i_n} \theta_{ijk} & = &
-\frac{-1}{\sqrt{1 - ((r_{ik}^2+r_{ij}^2-r_{jk}^2)/2r_{ik}r_{ij})^2}} \times \nonumber\\
- & & \times \Big( \frac{4 r_{ik}r_{ij} (2 x^i_n - x^k_n - x^j_n) + 2(x^j_n - x^i_n)\frac{r_{ik}}{r_{ij}} + 2(x^k_n - x^i_n)\frac{r_{ij}}{r_{ik}} }{4 r^2_{ik} r^2_{ij}}\Big) \label{eq:d_theta}
-\end{eqnarray}
-Using the expressions \eqref{eq:d_cutoff} and \eqref{eq:d_theta} the derivation of $b_{ij}$ with respect to $x^i_n$ can be written as:
-\begin{eqnarray}
-\partial_{x^i_n} b_{ij} & = &
-- \frac{1}{2n_i} \chi_{ij} \Bigg( 1 + \beta_i^{n_i} \Bigg[ \sum_{k \ne i,j} \bigg( f_C(r_{ij}) \omega_{ik} \Big( 1 + \frac{c_i^2}{d_i^2} - \frac{c_i^2}{d_i^2 + (h_i - \cos \theta_{ijk})^2} \Big) \bigg)^{n_i} \Bigg] \Bigg)^{-\frac{1}{2n_i} - 1} \times \nonumber\\
-&& \times n_i \beta_i^{n_i} \Bigg[ \sum_{k \ne i,j} f_C(r_{ik}) \omega_{ik} \Big( 1 \frac{c_i^2}{d_i^2} - \frac{c_i^2}{d_i^2 + (h_i - \cos \theta_{ijk})^2} \Big) \Bigg]^{n_i -1} \times \nonumber\\
-&& \times \sum_{k \ne i,j} \Bigg( \omega_{ik} \Big( 1 + \frac{c_i^2}{d_i^2} - \frac{c_i^2}{d_i^2 + (h_i - \cos \theta_{ijk})^2} \Big) \partial_{x^i_n} f_C(r_{ik}) + \nonumber\\
-&& + f_C(r_{ik}) \omega_{ik} (-1) \frac{c_i^2}{(d_i^2 + (h_i - \cos \theta_{ijk})^2)^2} \times \nonumber\\
-&& \times 2 \Big( h_i - \cos \theta_{ijk} \Big) \sin \theta_{ijk} \partial_{x^i_n} \theta_{ijk} \Bigg)
-\end{eqnarray}
+The volume of the synthesized material can hardly be controlled in experiment.
+Instead the pressure can be adjusted.
+Holding constant the pressure in addition to the temperature of the system its states are represented by the isothermal-isobaric $NpT$ ensemble.
+The expression for the pressure of a system derived from the equipartition theorem is given by
+\begin{equation}
+pV=Nk_{\text{B}}T+\langle W\rangle\text{, }W=-\frac{1}{3}\sum_i\vec{r}_i\nabla_{\vec{r}_i}U
+\text{, }
+\label{eq:basics:ps}
+\end{equation}
+where $W$ is the virial and $U$ is the configurational energy.
+
+Berendsen~et~al.~\cite{berendsen84} proposed a method, which is easy to implement, to couple the system to an external bath with constant temperature $T_0$ or pressure $p_0$ with adjustable time constants $\tau_T$ and $\tau_p$ determining the strength of the coupling.
+Control of the respective variable is based on the relations given in equations \eqref{eq:basics:ts} and \eqref{eq:basics:ps}.
+The thermostat is achieved by scaling the velocities of all atoms in every time step $\delta t$ from $\vec{v}_i$ to $\lambda \vec{v}_i$, with
+\begin{equation}
+\lambda=\left[1+\frac{\delta t}{\tau_T}(\frac{T_0}{T}-1)\right]^\frac{1}{2}
+\text{ ,}
+\end{equation}
+where $T$ is the current temperature according to equation \eqref{eq:basics:ts}.
+The barostat adjusts the pressure by changing the virial through scaling of the particle positions $\vec{r}_i$ to $\mu \vec{r}_i$ and the volume $V$ to $\mu^3 V$, with
+\begin{equation}
+\mu=\left[1-\frac{\beta\delta t}{\tau_p}(p_0-p)\right]^\frac{1}{3}\text{ ,}
+\end{equation}
+where $\beta$ is the isothermal compressibility and $p$ corresponds to the current pressure, which is determined by equation \eqref{eq:basics:ps}.
+
+Using this method the system does not behave like a true $NpT$ ensemble.
+On average $T$ and $p$ correspond to the expected values.
+For large enough time constants, i.e. $\tau > 100 \delta t$, the method shows realistic fluctuations in $T$ and $p$.
+The advantage of the approach is that the coupling can be decreased to minimize the disturbance of the system and likewise be adjusted to suit the needs of a given application.
+It provides a stable algorithm that allows smooth changes of the system to new values of temperature or pressure, which is ideal for the investigated problem.
+
+\section{Denstiy functional theory}
+\label{section:dft}
+
+Dirac declared that chemistry has come to an end, its content being entirely contained in the powerul equation published by Schr\"odinger in 1926 \cite{schroeder26} marking the beginning of wave mechanics.
+Following the path of Schr\"odinger the problem in quantum-mechanical modeling of describing the many-body problem, i.e. a system of a large amount of interacting particles, is manifested in the high-dimensional Schr\"odinger equation for the wave function $\Psi({\vec{R}},{\vec{r}})$ that depends on the coordinates of all nuclei and electrons.
+The Schr\"odinger equation contains the kinetic energy of the ions and electrons as well as the electron-ion, ion-ion and electron-electron interaction.
+This cannot be solved exactly and finding approximate solutions requires several layers of simplification in order to reduce the number of free parameters.
+Approximations that consider a truncated Hilbert space of single-particle orbitals yield promising results, however, with increasing complexity and demand for high accuracy the amount of Slater determinats to be evaluated massively increases.
+
+In contrast, instead of using the description by the many-body wave function, the key point in density functional theory (DFT) is to recast the problem to a description utilizing the charge density $n(\vec{r})$, which constitutes a quantity in real space depending only on the three spatial coordinates.
+In the following sections the basic idea of DFT will be outlined.
+As will be shown, DFT can formally be regarded as an exactification of the Thomas Fermi theory \cite{thomas27,fermi27} and the self-consistent Hartree equations \cite{hartree28}.
+A nice review is given in the Nobel lecture of Kohn \cite{kohn99}, one of the inventors of DFT.
+
+\subsection{Born-Oppenheimer approximation}
+
+Born and Oppenheimer proposed a simplification enabling the effective decoupling of the electronic and ionic degrees of freedom \cite{born27}.
+Within the Born-Oppenheimer (BO) approximation the light electrons are assumed to move much faster and, thus, follow adiabatically to the motion of the heavy nuclei, if the latter are only slightly deflected from their equilibrium positions.
+Thus, on the timescale of electronic motion the ions appear at fixed positions.
+On the other way round, on the timescale of nuclear motion the electrons appear blurred in space adding an extra term to the ion-ion potential.
+The simplified Schr\"odinger equation no longer contains the kinetic energy of the ions.
+The momentary positions of the ions enter as fixed parameters and, therefore, the ion-ion interaction may be regarded as a constant added to the electronic energies.
+The Schr\"odinger equation describing the remaining electronic problem reads
+\begin{equation}
+\left[-\frac{\hbar^2}{2m}\sum_j\nabla^2_j-
+\sum_{j,l} \frac{Z_le^2}{|\vec{r}_j-\vec{R}_l|}+
+\frac{1}{2}\sum_{j\neq j'}\frac{e^2}{|\vec{r}_j-\vec{r}_{j'}|}
+\right] \Psi = E \Psi
+\text{ ,}
+\end{equation}
+where $Z_l$ are the atomic numbers of the nuclei and $\Psi$ is the many-electron wave function, which depends on the positions and spins of the electrons.
+Accordingly, there is only a parametrical dependence on the ionic coordinates $\vec{R}_l$.
+However, the remaining number of free parameters is still too high and need to be further decreased.