Я использую следующий код в 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.