Spin-1 models in the ultrastrong coupling regime of circuit QED

We propose a superconducting circuit platform for simulating spin-1 models. To this purpose we consider a chain of N ultrastrongly coupled qubit-resonator systems interacting through a grounded SQUID. The anharmonic spectrum of the qubit-resonator system and the selection rules imposed by the global parity symmetry allow us to activate well controlled two-body quantum gates via AC-pulses applied to the SQUID. We show that our proposal has the same simulation time for any number of spin-1 interacting particles. This scheme may be implemented within the state-of-the-art circuit QED in the ultrastrong coupling regime.