Hello, I stumbled upon your project while searching for existing codes to deal with products of fermionic and bosonic operators. I am not quite familiar with Python programming, but I think I found an inconsistency in your code-base (not quite sure if this is fixed in secondquant). From a physical point of view powers of a fermionic operator other than op^1 do not make sense. From Pauli's exclusion principle one would expect something like a^n=S.Zero for n>=2 if a is any fermionic operator.
Another point inside sympy.physics.quantum.operatorordering: Since one has the fermionic anticommutator-relation: Switching two 'independent' fermionic operators should switch the sign of the total expression, i.e. a*b = -b*a for any two independent fermionic operators a and b. In addition to the aforementioned points I have implemented a few tweaks to your library to deal with indexed FermionOp's. I am just learning how to use git so I don't know how to properly set up pull requests. Nevertheless I would like to share my results. Kind regars Peter -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To post to this group, send email to [email protected]. Visit this group at https://groups.google.com/group/sympy. To view this discussion on the web visit https://groups.google.com/d/msgid/sympy/fc2c5524-5b32-4f27-bcde-c76b214da572%40googlegroups.com. For more options, visit https://groups.google.com/d/optout.
