Non-Gaussian superradiant transition via three-body ultrastrong coupling

We introduce a class of quantum optical Hamiltonians characterized by three-body couplings and propose a circuit-QED scheme based on state-of-the-art technology that implements the considered model. Unlike two-body light-matter interactions, this three-body coupling Hamiltonian is exclusively composed of terms which do not conserve the particle number. We explore the three-body ultrastrong-coupling regime, showing the emergence of a superradiant phase transition which is of first order, is characterized by the breaking of Z2 x Z2 symmetry, and has a strongly non-Gaussian nature. Indeed, in contrast to what is observed in any two-body-coupling model, in proximity of the transition the ground state exhibits a divergent coskewness, i.e., quantum correlations that cannot be captured within semiclassical and Gaussian approximations. Furthermore, we demonstrate the robustness of our findings by including dissipative processes in the model, showing that the steady state of the system inherits from the ground states the most prominent features of the transition.