This work is an extensive study of the spectral statistics of three representative classically integrable systems, namely rectangle, torus and circle billiards. We analyse the statistics and focus on the related level spacing distribution P(S) and the delta statistics . The agreement with the Poissonian model is typically found to be perfect up to the outer (unfolded) energy scale , beyond which the saturation is observed, in agreement with Berry's dynamical theory of the spectral rigidity, where as .
The untypical systems are those, where for example P(S) is not a smooth distribution but a sum of the delta functions, due to the `granularity' of the energy scale, for example in the rectangle with the rational squared sides ratio. However, even there we find reasonable trend towards Poissonian statistics for large ranges L but . We describe theoretically and numerically the broadening of the delta spikes when the rectangular billiard is slightly distorted away from a rational to an irrational shape and find excellent agreement. Also, in irrational rectangle billiards we show and explain the existence of large fluctuations, by one order of magnitude bigger than the statistical ones, whose origin is in the closeness to some rational billiard shape. These fluctuations and their amplitude are independent of the energy if the bin size shrinks inversely with energy.
Finally we tested the mode distribution (i.e. distribution of the reduced mode fluctuation number W) and found that it was not Poissonian, in agreement with Steiner's conjecture, and in fact follows the prediction by Bleher et al, that its tail behaves as . The general reason for non-universal behaviour of P(W) is that at largest energy scales we are always in the saturation regime .