Numerically exact treatment of many body self-organization in a cavity

We investigate the full quantum evolution of ultracold interacting bosonic atoms on a chain and coupled to an optical cavity. Extending the time-dependent matrix product state techniques and the many-body adiabatic elimination technique to capture the global coupling to the cavity mode and the open nature of the cavity, we examine the long time behavior of the system beyond the mean-field elimination of the cavity field. We investigate the many body steady states and the self-organization transition for a wide range of parameters. We show that in the self-organized phase the steady state consists in a mixture of the mean-field predicted density wave states and excited states with additional defects. In particular, for large dissipation strengths a steady state with a fully mixed atomic sector is obtained crucially different from the predicted mean-field state.