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.
We propose a simple scheme to construct composition-dependent interatomic potentials for multicomponent systems that when superposed onto the potentials for the pure elements can reproduce not only the heat of mixing of the solid solution in the enti
re concentration range but also the energetics of a wider range of configurations including intermetallic phases. We show that an expansion in cluster interactions provides a way to systematically increase the accuracy of the model, and that it is straightforward to generalise this procedure to multicomponent systems. Concentration-dependent interatomic potentials can be built upon almost any type of potential for the pure elements including embedded atom method (EAM), modified EAM, bond-order, and Stillinger-Weber type potentials. In general, composition-dependent N-body terms in the total energy lead to explicit (N+1)-body forces, which potentially renders them computationally expensive. We present an algorithm that overcomes this problem and that can speed up the calculation of the forces for composition-dependent pair potentials in such a way as to make them computationally comparable in efficiency and scaling behaviour to standard EAM potentials. We also discuss the implementation in Monte-Carlo simulations. Finally, we exemplarily review the composition-dependent EAM model for the Fe-Cr system [PRL 95, 075702 (2005)].
