We explore theoretically the nonequilibrium photonic phases of an array of coupled cavities in presence of incoherent driving and dissipation. In particular, we consider a Hubbard model system where each site is a Kerr nonlinear resonator coupled to a two-level emitter, which is pumped incoherently. Within a Gutzwiller mean-field approach, we determine the steady-state phase diagram of such a system. We find that, at a critical value of the intercavity photon hopping rate, a second-order nonequilibrium phase transition associated with the spontaneous breaking of the U(1) symmetry occurs. The transition from an incompressible Mott-like photon fluid to a coherent delocalized phase is driven by commensurability effects and not by the competition between photon hopping and optical nonlinearity. The essence of the mean-field predictions is corroborated by finite-size simulations obtained with matrix product operators and corner-space renormalization methods

We report a theoretical study of a quantum optical model consisting of an array of strongly nonlinear cavities incoherently pumped by an ensemble of population-inverted two-level atoms. Projective methods are used to eliminate the atomic dynamics and write a generalized master equation for the photonic degrees of freedom only, where the frequency-dependence of gain introduces non-Markovian features. In the simplest single cavity configuration, this pumping scheme gives novel optical bistability effects and allows for the selective generation of Fock states with a well-defined photon number. For many cavities in a weakly non-Markovian limit, the non-equilibrium steady state recovers a Grand-Canonical statistical ensemble at a temperature determined by the effective atomic linewidth. For a two-cavity system in the strongly nonlinear regime, signatures of a Mott state with one photon per cavity are found.

We introduce a frequency-dependent incoherent pump scheme with a square-shaped spectrum as a way to study strongly correlated photons in arrays of coupled nonlinear resonators. This scheme can be implemented via a reservoir of population-inverted two-level emitters with a broad distribution of transition frequencies. Our proposal is predicted to stabilize a nonequilibrium steady state sharing important features with a zero-temperature equilibrium state with a tunable chemical potential. We confirm the efficiency of our proposal for the Bose-Hubbard model by computing numerically the steady state for finite system sizes: first, we predict the occurrence of a sequence of incompressible Mott-insulator-like states with arbitrary integer densities presenting strong robustness against tunneling and losses. Secondly, for stronger tunneling amplitudes or noninteger densities, the system enters a coherent regime analogous to the superfluid state. In addition to an overall agreement with the zero-temperature equilibrium state, exotic nonequilibrium processes leading to a finite entropy generation are pointed out in specific regions of parameter space. The equilibrium ground state is shown to be recovered by adding frequency-dependent losses. The promise of this improved scheme in view of quantum simulation of the zero-temperature many-body physics is highlighted

No abstract