  The dynamics of simple qubit systems in a chain configuration coupled at both
ends to separate bosonic baths at different temperatures is studied. An exact
analytical solution of the master equation in the Born-Markov approximation for
the reduced density matrix of the qubit system is constructed. The unique
non-equilibrium stationary state for the long time behavior of the reduced
density matrix in obtained. Dynamical and steady state properties of the
concurrence between the first and the last spin are studied.
