Dynamical properties of dissipative XYZ Heisenberg lattices

We study dynamical properties of dissipative XYZ Heisenberg lattices where anisotropic spin-spin coupling competes with local incoherent spin flip processes. In particular, we explore a region of the parameter space where dissipative magnetic phase transitions for the steady state have been recently predicted by mean-field theories and exact numerical methods. We investigate the asymptotic decay rate towards the steady state both in 1D (up to the thermodynamical limit) and in finite-size 2D lattices, showing that critical dynamics does not occur in 1D, but it can emerge in 2D. We also analyze the behavior of individual homodyne quantum trajectories, which well reveal the nature of the transition

Variational neural network ansatz for steady states in open quantum systems

We present a general variational approach to determine the steady state of open quantum lattice systems via a neural network approach. The steady-state density matrix of the lattice system is constructed via a purified neural network ansatz in an extended Hilbert space with ancillary degrees of freedom. The variational minimization of cost functions associated to the master equation can be performed using a Markov chain Monte Carlo sampling. As a first application and proof-of-principle, we apply the method to the dissipative quantum transverse Ising model.Comment: 6 pages, 4 figures, 54 references, 5 pages of Supplemental Information

Photon transport in a dissipative chain of nonlinear cavities

We analyze a chain of coupled nonlinear optical cavities driven by a coherent source of light localized at one end and subject to uniform dissipation. We characterize photon transport by studying the populations and the photon correlations as a function of position. When complemented with input-output theory, these quantities provide direct information about photon transmission through the system. The position of single- and multi-photon resonances directly reflect the structure of the many-body energy levels. This shows how a study of transport along a coupled cavity array can provide rich information about the strongly correlated (many-body) states of light even in presence of dissipation. By means of a numerical algorithm based on the time-evolving block decimation scheme adapted to mixed states, we are able to simulate arrays up to sixty cavities.Comment: 12 pages, 14 figure

Linked cluster expansions for open quantum systems on a lattice

We propose a generalization of the linked-cluster expansions to study driven-dissipative quantum lattice models, directly accessing the thermodynamic limit of the system. Our method leads to the evaluation of the desired extensive property onto small connected clusters of a given size and topology. We first test this approach on the isotropic spin-1/2 Hamiltonian in two dimensions, where each spin is coupled to an independent environment that induces incoherent spin flips. Then we apply it to the study of an anisotropic model displaying a dissipative phase transition from a magnetically ordered to a disordered phase. By means of a Pad\'e analysis on the series expansions for the average magnetization, we provide a viable route to locate the phase transition and to extrapolate the critical exponent for the magnetic susceptibility.Comment: 10 pages, 5 figure

Energy transport between two integrable spin chains

We study the energy transport in a system of two half-infinite XXZ chains initially kept separated at different temperatures, and later connected and let free to evolve unitarily. By changing independently the parameters of the two halves, we highlight, through bosonisation and time-dependent matrix-product-state simulations, the different contributions of low-lying bosonic modes and of fermionic quasi-particles to the energy transport. In the simulations we also observe that the energy current reaches a finite value which only slowly decays to zero. The general pictures that emerges is the following. Since integrability is only locally broken in this model, a pre-equilibration behaviour may appear. In particular, when the sound velocities of the bosonic modes of the two halves match, the low-temperature energy current is almost stationary and described by a formula with a non-universal prefactor interpreted as a transmission coefficient. Thermalisation, characterized by the absence of any energy flow, occurs only on longer time-scales which are not accessible with our numerics.Comment: 15 pages, 14 figure

Stabilizing strongly correlated photon fluids with non-Markovian reservoirs

We introduce a novel 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 non-equilibrium 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 non-integer 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 non-equilibrium 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