We investigate the basic quantum mechanical processes behind non-proportional response of scintillators to incident radiation responsible for reduced resolution. For this purpose, we conduct a comparative first principles study of quasiparticle spect
ra on the basis of the $G_0W_0$ approximation as well as absorption spectra and excitonic properties by solving the Bethe-Salpeter equation for two important systems, NaI and SrI$_2$. The former is a standard scintillator material with well-documented non-proportionality while the latter has recently been found to exhibit a very proportional response. We predict band gaps for NaI and SrI$_2$ of 5.5 and 5.2 eV, respectively, in good agreement with experiment. Furthermore, we obtain binding energies for the groundstate excitons of 216 meV for NaI and 195$pm$25 meV for SrI$_2$. We analyze the degree of exciton anisotropy and spatial extent by means of a coarse-grained electron-hole pair-correlation function. Thereby, it is shown that the excitons in NaI differ strongly from those in SrI$_2$ in terms of structure and symmetry, even if their binding energies are similar. Furthermore, we show that quite unexpectedly the spatial extents of the highly anisotropic low-energy excitons in SrI$_2$ in fact exceed those in NaI by a factor of two to three in terms of the full width at half maxima of the electron-hole pair-correlation function.
Excitonic effects in optical spectra and electron-hole pair excitations are described by solutions of the Bethe-Salpeter equation (BSE) that accounts for the Coulomb interaction of excited electron-hole pairs. Although for the computation of excitoni
c optical spectra in an extended frequency range efficient methods are available, the determination and analysis of individual exciton states still requires the diagonalization of the electron-hole Hamiltonian $hat{H}$. We present a numerically efficient approach for the calculation of exciton states with quadratically scaling complexity, which significantly diminishes the computational costs compared to the commonly used cubically scaling direct-diagonalization schemes. The accuracy and performance of this approach is demonstrated by solving the BSE numerically for the Wannier-Mott two-band model in {bf k} space and the semiconductors MgO and InN. For the convergence with respect to the $vk$-point sampling a general trend is identified, which can be used to extrapolate converged results for the binding energies of the lowest bound states.
