我在sympy中使用以下代码：
from sympy.physics.secondquant import F, Fd, NO, Commutator
from sympy import symbols
a, b, c, d = symbols("a,b,c,d")
comm = NO(Commutator( Fd(a) * F(b), Fd(c) * F(d) ).doit())
print(comm)
# gives me 0
很显然，这不应为零。
好吧，也许我听不懂sympy。这是我要计算的：
[F†_a F_b，F†_c F_d]的索引a，b，c和d不一定相等。