I am using the quantum mathematica package and am trying to do computations that require me to use the outer (tensor) products of Pauli matrices. I would define this as, for example, Subscript[\[Sigma], 2]\[CircleTimes]Subscript[\[Sigma], 3] . However if I ask mathematica to compute this, it only finds the matrix product of the two operators, not the outer (tensor) product. Does anyone know of a way to do this calculation with Mathematica still recognising the operators as Pauli matrices, i.e. rather than writing them in matrix form or in Dirac notation?