The characterization of semiconductors applied in photovoltaic devices by means of ab-initio methods has become an emerging field over recent years. This is not only because of the increasing availability of more powerful computers, but mostly due to the development of refined numerical and theoretical methods.Explaining and predicting solid-state properties requires an understanding of the behavior of electrons in solids. First-principles or ab-initio calculations are those that start directly at the Hamiltonian of the Schr€ odinger equation containing the kinetic energy operators and the potential energy due to electrostatic interactions of electrons and nuclei. The value of ab-initio methods lies in the fact that material properties can be calculated from scratch without experimental input or empirical parameters. Obtaining materials data from ab-initio calculations may serve various purposes. Firstly, they allow for validating experimental results and therefore contribute to a better understanding of material properties. Secondly, calculated data can be used for identifying new properties or mechanisms. Thirdly, by calculating specific properties for a large database of various structures, new materials with optimized properties may be designed. In semiconductor research, ab-initio methods are now frequently used for calculating structural, thermomechanical and electronic properties. In practice, however, the feasibility of an ab-initio approach always depends on the problem at hand and may be limited for conceptual reasons. Most importantly, the computational costs which scale with the system size, that is, the number of electrons necessary to calculate a certain property of a material, restrict the applicability of ab-initio methods. Moreover, it is often very crucial to balance accuracy and computational cost by choosing the appropriate method for the problem at hand. The most accurate approach of treating a complex system consisting of many interacting particles, such as a semiconducting crystal, is ab-initio methods solving the many-body Schr€ odinger equation directly. This includes quantum Monte Carlo [1], configuration interaction (CI), and coupled cluster (CC) methods [2], which are computationally expensive and therefore more often applied for studying small