Conjugation of powers of subscripted constants in mathematica

[joshkamm]

Description

Not sure if this is the right place to post about the SNEG library underlying this library but I didn't see a specific repository for it.

There seems to be a quirk in the low level manipulation of Dirac notation.

Here is a minimal reproducing code example which is a simpler version of the more complicated manipulation I was originally performing which got stuck simplifying.

Code example

<< "sneg-1.250/sneg.m"
snegrealconstants[c, Subscript[c, 1]]
conj[c*VACUUM]
conj[c*c*VACUUM]
conj[Subscript[c, 1]*VACUUM]
conj[Subscript[c, 1]*Subscript[c, 1]*VACUUM]

Output:

Expected behavior: SNEG recognizes a power of a subscripted constant like the other constant above

Actual behavior: SNEG appears to not recognize the power of a real constant containing a subscript.

Versions

Mathematica 13.0 SNEG 1.250

[rokzitko]

I'll reply here, since there is indeed no separate repository for SNEG. The best solution is to simply define the required properties:

snegnonopQ[Subscript[c, ___]] := True;
Conjugate[x : Subscript[c, ___]] ^= x;

The first line specifies that c with a subscript is not an operator, and the second one declares it to be a real quantity.

[rokzitko]

Just for the record, there is now a github repo for sneg: https://github.com/rokzitko/sneg