“…(21) are given, in theσ z,i basis, by: |1 = |+, + , |2 = (|+, − − |−, + )/ √ 2, |3 = (|+, − + |−, + )/ √ 2 and |4 = |−, − . The non-equilibrium thermal steady-state density operator for the pair of qubits turns out to be diagonal in this basis taking the form of a direct product as [15] ρ ss =ρ ss…”