Exact Nonequilibrium Steady State of an Open Hubbard Chain

We discuss current carrying non-equilibrium steady state of an open fermionic Hubbard chain that is strongly driven by markovian incoherent processes localized at the chain ends. An explicit form of exact many-body density operator for any value of the coupling parameter is presented. The structure of a matrix product form of the solution is encoded in terms of a novel diagrammatic technique which should allow for generalization to other integrable non-equillibrium models.