A solvable model for strongly interacting nonequilibrium excitons

We study the driven-dissipative Bose-Hubbard model with an all-to-all hopping term in the system Hamiltonian, while subject to incoherent pumping and decay from the environment. This system is naturally probed in several recent experiments on excitons in WS₂/WSe₂ moiré systems, as well as quantum simulators. By positing a particular form of coupling to the environment, we derive the Lindblad jump operators and show that, in certain limits, the system admits a closed-form expression for the steady-state density matrix. Away from the exactly solvable regions, the steady state can be obtained numerically for 100s to 1,000s of sites. We study the nonequilibrium phase diagram and phase transitions, which qualitatively matches the equilibrium phase diagram, agreeing with the intuition that increasing the intensity of the light is equivalent to changing the bosonic chemical potential. However, the steady states are far from thermal states, and the nature of the phase transitions is changed.