Manipulating energy and spin currents in non-equilibrium systems of interacting qubits

We consider a generic interacting chain of qubits, which are coupled at the edges to baths of fixed polarizations. We can determine the non-equilibrium steady states, described by the fixed point of the Lindblad master equation. Under rather general assumptions about local pumping and interactions, symmetries of the reduced density matrix are revealed. The symmetries drastically restrict the form of the steady density matrices in such a way that an exponentially large subset of one-point and many-point correlation functions are found to vanish. As an example we show how in a Heisenberg spin chain a suitable choice of the baths can completely switch off either the spin or the energy current, or both of them, despite the presence of large boundary gradients.