Quantization of a generic chaotic 3D billiard with smooth boundary. I. Energy level statistics

Numerical calculation and analysis of extremely high-lying energy spectra, containing thousands of levels with sequential quantum number up to 62000 per symmetry class, of a generic chaotic 3D quantum billiard is reported. The shape of the billiard is given by a simple and smooth deformation of a unit sphere, which gives rise to nearly fully chaotic classical dynamics. We present an analysis of (i) the quantum length spectrum whose smooth part agrees with the 3D Weyl formula and whose oscillatory part is peaked around the periods of classical periodic orbits, (ii) the nearest neighbor level spacing distribution and (iii) the number variance. Although the chaotic classical dynamics quickly and uniformly explores almost the entire energy shell, while the measure of the regular part of phase space is insignificantly small, we find small but significant deviations from GOE statistics which are explained in terms of localization of eigenfunctions onto lower dimensional classically invariant manifolds.