Wigner-Dyson statistics for a class of integrable models

We construct an ensemble of second-quantized Hamiltonians with two bosonic degrees of freedom, whose members display with probability one Gaussian orthogonal ensemble (GOE) or Gaussian unitary ensemble (GUE) statistics. Nevertheless, these Hamiltonians have a second integral of motion, namely, the boson number, and thus are integrable. To construct this ensemble we use some "reverse engineering" starting from the fact that n bosons in a two-level system with random interactions have an integrable classical limit by the old Heisenberg association of boson operators to actions and angles. By choosing an n-body random interaction and degenerate levels we end up with GOE or GUE Hamiltonians. Ergodicity of these ensembles completes the example.