as the composition of the sample changes from Cu-rich to Cu-poor, so that in the Cupoor samples only one peak is detected. Although only the Cu-poor material is important for actual devices Cu-rich samples give indispensable information for defect identification.First-principles calculations based on the density functional theory (DFT) can be used to obtain important, complementary information about point defects such as formation energies and charge transition levels. [6] A plethora of studies concerning defects in CuInSe 2 and CuGaSe 2 has been published over the last two decades. [7][8][9][10][11][12][13] The most recent DFT investigations have employed hybrid functionals, which can overcome the energy band gap problem plaguing the older studies and can thus also provide information about the defect level positions within the band gap. [9][10][11][12][13] The results of different calculations agree well with respect to general trends in formation energies of the most important defects, such as the copper vacancy (V Cu ), indium antisite on copper place (In Cu ), and copper interstitial (Cu int ). However, the results differ in some important cases. For instance, there are clearly different values for the ionization levels within the band gap for the copper antisite on the indium place (Cu In ) and indium vacancy (V In ). Based on the formation energy calculations, V Cu and Cu In are abundant acceptors, and are most probably responsible for some of the above-mentioned PL peaks, but it is unclear whether any native defect can be responsible for the third acceptor level.One important goal of the present work was to gain a perspective on the present unsatisfactory situation in modeling point defects in CIGSe and to approach the ultimate accuracy by which DFT is able to predict the properties of bulk crystalline materials. [14] First we carried out a detailed benchmarking of the first-principles computational scheme used. We checked effects due to the supercell size and shape, as well as those of the finite-size supercell correction scheme. We used also two very different implementations of the first-principles DFT method, which differ in describing valence-core electron interaction and electron wave functions (see below). After finding the computational parameters yielding accurate results, we calculated formation energies and charge transition levels for different acceptor candidates in CuInSe 2 . By carefully considering the relevant chemical potential limits, we were able to draw conclusions about the abundances of different defects. In addition to simple native defects, we have also considered a set of complexes formed by them. Our paper is organized as follows.