“…( 12) we have introduced g, the coupling strength between two-level systems and the boson field, and the parameter 0 ≤ α ≤ 1. When α = 0 or α = 1, the Dicke Hamiltonian reduces to the Tavis-Cummings model, which is integrable and has an additional conserved quantity, the total number of excitations in the system, Mk = Ĵ + Ĵz + b † b, see; otherwise, for 0 < α < 1, the model is in the chaotic regime [22,[34][35][36]38]. Thus, we can analyse the behaviour of this system in the integrable and non-integrable limits simply by considering cases with α ∈ {0, 1} and α ∈ {0, 1}, respectively.…”