Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

FermionSpinSum fails with spinor of combined momenta #167

Open
ChiMaoShuPhy opened this issue Apr 13, 2022 · 1 comment
Open

FermionSpinSum fails with spinor of combined momenta #167

ChiMaoShuPhy opened this issue Apr 13, 2022 · 1 comment

Comments

@ChiMaoShuPhy
Copy link

ChiMaoShuPhy commented Apr 13, 2022

  • Your Mathematica version

```12.3.1 for Linux x86 (64-bit) (June 24, 2021)``

  • Your FeynCalc version

9.3.1

  • Did you try to reinstall FeynCalc (stable version) using the automatic installer to make sure that you have the latest bugfixes?

Yes

  • Does your Mathematica initialization file contain statements that might influence the behavior of FeynCalc? Sometimes external packages may modify `init.m` in unusual ways, causing troubles for other codes.

no

  • The ```FermionSpinSum``` fails if the spinor contains combination of momenta such as p+k, p/2+k

SpinorUBar[p,m].SpinorU[p,m]//FermionSpinSum works
SpinorUBar[p+k,m].SpinorU[p+k,m]//FermionSpinSum reports error FermionSpinSum: Error! FermionSpinSum encountered a fatal problem and must abort the computation. The problem reads: The input contains Spinor objects with incorrect syntax.

  • Same error was reported for SpinorV and SpinorVBar
@vsht
Copy link
Member

vsht commented Apr 13, 2022

I'm aware of this behavior, but I don't really see a way to fix this.

Functions that deal with spinor chains expect each spinor to depend on
a single momentum like p or q. This is mainly related to the way how FeynCalc
figures out whether the given Spinor is a u, v, ubar or vbar. The sign of the momentum
is crucial here, cf.

https://feyncalc.github.io/FeynCalcBookDev/Spinor.html

Now if you make a spinor depend on a linear combination of momenta, all
this machinery breaks down and one runs into very weird bugs, cf. #76
The worst thing is that you'd get incorrect results without even noticing this.
This is why I explicitly added a check to forbid such input.

For nonrelativistic expansions of spinor chains I advocate the approach where
one first rewrites everything in terms of Pauli matrices and stuff (using FeynOnium routines)
and then substitutes the decomposition of the heavy particles momenta with FCReplaceMomenta.
This is how I handle this in my calculations.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants