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. 
H.-J. Stoeckmann