多体Pauli算子乘积与对易式化简程序优化求助

我在Physics Stack Exchange上看到一个未迁移到计算科学板块的问题，需求是化简涉及多体Pauli矩阵乘积与对易式的表达式，示例如下：

我尝试使用SymPy进行处理，代码如下：

from sympy import *
from sympy.physics.paulialgebra import Pauli, evaluate_pauli_product
from sympy import I
from sympy.physics.quantum import TensorProduct as tp
from sympy.physics.quantum import Commutator as cm
from sympy.physics.quantum import tensor_product_simp as tpsimp
Id = 1
X = Pauli(1)
Y = Pauli(2)
Z = Pauli(3)
term1 = tp(X,Y)
term2 = tp(Id,cm(Y,X).doit())
result = tpsimp(cm(term1,term2).doit())
print(evaluate_pauli_product(result))

但输出结果2*I*(1xsigma3*sigma1xsigma2 - sigma1xsigma2*1xsigma3)远未达到简化效果，也不知道该如何指定张量积对应的索引，求更好的解决方法。

内容的提问来源于stack exchange，提问作者Andrew Hardy