**Next message:**Nikita Belyaev: "Re: Imaginary parts and Schouten identity"**Previous message:**Nikita Belyaev: "Imaginary parts and Schouten identity"**In reply to:**Nikita Belyaev: "Imaginary parts and Schouten identity"**Next in thread:**Nikita Belyaev: "Re: Imaginary parts and Schouten identity"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ]**Mail actions:**[ respond to this message ] [ mail a new topic ]

Hi Nikita,

actually I have some doubts that your expression should have no

imaginary part. At least, if I look at the pieces of the imaginary part

that are proportional to u^3:

SelectNotFree[FCE[TrA2B2], Complex];

u3Piece = SelectNotFree[%, u^3] // EpsEvaluate

-512 I u^3 SP[p1, p2] LC[][k1, k2, p, s] +

512 I u^3 SP[p, p2] LC[][k1, k2, p1, s] -

512 I u^3 SP[p, p1] LC[][k1, k2, p2, s] +

512 I u^3 SP[k1, p2] LC[][k2, p, p1, s] -

512 I u^3 SP[k1, p] LC[][k2, p1, p2, s]

Schouten[%]

then it is clear that they do not vanish by the Schouten identity.

However since these are the only terms that contain u^3, I don't think

that they can be cancelled by other terms in the imaginary part.

Cheers,

Vladyslav

Am 25.07.2015 um 18:39 schrieb Nikita Belyaev:

*> Hello, Vladyslav
*

*>
*

*> We've faced with a problem that we couldn't solve. It is related to imaginary parts of the traces. We have some of them to calculate and a lot of them were already calculated correctly but for one specific term we cannot get the result without imaginary part.
*

*> We've spent a lot of time using a FeynCalcFormLink procedure and also trying to apply the Scouten identity as you told us some time ago but we did not succeed because as I understand there is no common algorithm to check such kind of equivalence in the general case.
*

*>
*

*> Could you probably help us to check that the trace included contains no imaginary part?
*

*>
*

*> Here is the code of our calculation:
*

*> Clear["Global`*"];
*

*> <<HighEnergyPhysics`FeynCalc`
*

*> Needs["FeynCalcFormLink`"]
*

*> $LeviCivitaSign = -1;
*

*>
*

*> ScalarProduct[p,p]=m^2;
*

*> ScalarProduct[p1,p1]=m^2;
*

*> ScalarProduct[p2,p2]=m^2;
*

*> ScalarProduct[k1,k1]=0;
*

*> ScalarProduct[k2,k2]=0;
*

*> ScalarProduct[q,q] = u^2;
*

*> ScalarProduct[q,s] =0;
*

*>
*

*> Line13:= (GS[p2]-m).GA[\[Beta]1].(GS[p1]+m).GA[\[Beta]].(GS[p]-m).GA[\[Alpha]1].GS[k2].GA[\[Alpha]].(1-GA[5]);
*

*> Line14:= GS[k1].GA[\[Alpha]1].(GS[q]-GS[p1]-GS[p2]-u).GA[\[Beta]1].(GS[q]-u).(1+GA[5].GS[s]).GA[\[Beta]].(GS[q]-GS[p1]-GS[p]-u).GA[\[Alpha]].(1-GA[5]);
*

*> Line15:= (GS[p]-m).GA[\[Beta]1].(GS[p1]+m).GA[\[Beta]].(GS[p2]-m).GA[\[Alpha]1].GS[k2].GA[\[Alpha]].(1-GA[5]);
*

*> Line16:= GS[k1].GA[\[Alpha]1].(GS[q]-GS[p1]-GS[p]-u).GA[\[Beta]1].(GS[q]-u).(1+GA[5].GS[s]).GA[\[Beta]].(GS[q]-GS[p1]-GS[p2]-u).GA[\[Alpha]].(1-GA[5]);
*

*>
*

*> Tr13= DiracTrace[Line13];
*

*> Tr14= DiracTrace[Line14];
*

*> Tr15= DiracTrace[Line15];
*

*> Tr16= DiracTrace[Line16];
*

*>
*

*> TrA2B2=FeynCalcFormLink[Tr13.Tr14+Tr15.Tr16];
*

*>
*

*> ComplexExpand[FullSimplify[TrA2B2/1024]] //Schouten
*

*>
*

*> Thanks for your help!
*

*>
*

*> With the Best Regards,
*

*> Nikita Belyaev
*

*>
*

**Next message:**Nikita Belyaev: "Re: Imaginary parts and Schouten identity"**Previous message:**Nikita Belyaev: "Imaginary parts and Schouten identity"**In reply to:**Nikita Belyaev: "Imaginary parts and Schouten identity"**Next in thread:**Nikita Belyaev: "Re: Imaginary parts and Schouten identity"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ]**Mail actions:**[ respond to this message ] [ mail a new topic ]

*
This archive was generated by hypermail 2b29
: 03/23/17-07:00:03 AM Z CET
*