I have a map set as

$ Phi (X) = a ^ 2 Tr[X] | 0 rangle langle 0 | + b ^ 2 Tr[X] | 1 rangle langle 1 | + ab Tr[sigma_z X] | 0 rangle langle 1 | + ab Tr[sigma_z X] | 1 row langle 0 |.

$

Right here $ | 0 rangle = begin {pmatrix} 1 \ 0 end {pmatrix} $ and $ | 1 rangle = begin {pmatrix} 0 \ 1 end {pmatrix} $. I want to perform the following operation

$ sum limits_ {j = 0} ^ {1} sum limits_ {k = 0} ^ {1} Psi ( j rangle langle k |) otimes | j rangle langle k |. $

How can this operation be done on Mathematica?