I am trying to use the SNEG Library for evaluating some commutators. For instance I am trying to evaluate $$ \left[n_1^2,a^\dagger_1\right]$$ with the code in Mathematica
<< sneg`; snegbosonoperators[a]; nSquare = (nc[number[a], number[a]]); aCR1 = a[CR, 1]; Comm1 = komutator[nSquare, aCR1] // SnegFullSimplify
The result using SNEG is $$ a_1^\dagger + 2 a_1^\dagger \cdot a_1^\dagger \cdot a_1 $$.
However I would like to simplify it to $$ a_1^\dagger + 2 a_1^\dagger \cdot n_1 $$, automatically.
Is there any way to tell SNEG to simplify the expression in this way?