A Supersymmetry approach to billiards with randomly distributed scatterers

The density of states for a chaotic billiard with randomly distributed point-like scatterers is calculated, doubly averaged over the positions of the impurities and the shape of the billiard. Truncating the billiard Hamiltonian to a N x N matrix, an explicit analytic expression is obtained for the case of broken time-reversal symmetry, depending on rank N of the matrix, number L of scatterers, and strength of the scattering potential. In the strong coupling limit a discontinuous change is observed in the density of states as soon as L exceeds N.