**Next message:**Jongping Hsu: "Re: Question about Subscript[E, 0](0,0,Overscript[p, _]^2,0,Overscript[p, _]^2,0,Overscript[p, _]^2,Overscript[p, _]^2,Overscript[p, _]^2,Overscript[p, _]^2,0,0,0,0,0)"**Previous message:**Vladyslav Shtabovenko: "Re: Question about TID in FC9.2.0"**Maybe in reply to:**Xiu-Lei Ren: "Question about TID in FC9.2.0"**Next in thread:**Jongping Hsu: "Re: Question about Subscript[E, 0](0,0,Overscript[p, _]^2,0,Overscript[p, _]^2,0,Overscript[p, _]^2,Overscript[p, _]^2,Overscript[p, _]^2,Overscript[p, _]^2,0,0,0,0,0)"**Reply:**Jongping Hsu: "Re: Question about Subscript[E, 0](0,0,Overscript[p, _]^2,0,Overscript[p, _]^2,0,Overscript[p, _]^2,Overscript[p, _]^2,Overscript[p, _]^2,Overscript[p, _]^2,0,0,0,0,0)"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ]**Mail actions:**[ respond to this message ] [ mail a new topic ]

*>> Unfortunately, at the moment one cannot evaluate it numerically directly
*

*>> from FeynHelpers (while developing the add-on my main focus were
*

*>> symbolic evaluations). However, you can easily do something like
*

*>>
*

*>> exp = PaXDiLog[2.7, 1] /. PaXDiLog -> X`DiLog
*

*>> Export["exp.m", exp]
*

*>>
*

*>> Quit[]
*

*>> << X`
*

*>> Import["exp.m"]
*

*>>
*

*>> to obtain the numerical value from Package-X. I'll contact the developer
*

*>> of Package-X to see if we can find a better solution...
*

Sorry, ignore this part. Somehow I completely forgot that the 1-loop library of Package-X is loaded not immediately but during the first call of PaXEvaluate.

So if you have already done some calculations with PaXEvaluate on a running kernel, then

exp = PaXDiLog[2.7, 1] /. PaXDiLog -> X`DiLog

is sufficient. On a fresh kernel just call PaXEvaluate once (with any input) and then it will work as well:

$LoadAddOns = {"FeynHelpers"};

<< FeynCalc`

PaXEvaluate[1]

exp = PaXDiLog[2.7, 1] /. PaXDiLog -> X`DiLog

Cheers,

Vladyslav

Am 25.12.2016 um 20:33 schrieb Vladyslav Shtabovenko:

*> Dear Xiu-Lei,
*

*>
*

*> actually, when doing the expansion by the means of Package-X, the LO
*

*> coefficient of the 1/mN expansion turns out to be zero (provided that I
*

*> got your example right):
*

*>
*

*> SPD[p4, p4] = mN^2;
*

*> XC0 = C0[SPD[p4], SPD[q], SPD[p4 + q], mN^2, mpi^2, mpi^2] //
*

*> ExpandScalarProduct;
*

*> XC0Re = PaXEvaluate[XC0, PaXC0Expand -> True,
*

*> PaXSeries -> {{mN, Infinity, 0}}, PaXAnalytic -> True] // Normal
*

*>
*

*> Doing the expansion with Series afterwards
*

*>
*

*> SPD[p4, p4] = mN^2;
*

*> XC0 = C0[SPD[p4], SPD[q], SPD[p4 + q], mN^2, mpi^2, mpi^2] //
*

*> ExpandScalarProduct;
*

*> XC0Re = PaXEvaluate[XC0, PaXC0Expand -> True] // Normal;
*

*> Series[XC0Re, {mN, Infinity, 0}] // Normal // Simplify
*

*>
*

*> produces several suspicious terms, like Sqrt[-SPD[q,q]]. As you probably
*

*> know, Mathematica is not always careful when choosing the branch cuts of
*

*> logs and square roots and does not really provide options to control
*

*> that consistently, so I would rather trust the output of Package-X
*

*> (which takes care of those things in a special way internally) than the
*

*> output of Series.
*

*>
*

*> By the way, the author of Package-X has released several fixes
*

*> in the meantime (ver 2.0.3 being the most current). One can update
*

*> Package-X manually, by downloading the tarball from
*

*> packagex.hepforge.org or via FeynHelpers' installer
*

*>
*

*> Import["https://raw.githubusercontent.com/FeynCalc/feynhelpers/master/\
*

*> install.m"]
*

*> InstallPackageX[]
*

*>
*

*> As to the second part of your question:
*

*>
*

*> PaXDiLog is just a placeholder that for the DiLog of Package-X. Its
*

*> relation to PolyLog is described in 1503.01469, Sec VI.
*

*>
*

*> Unfortunately, at the moment one cannot evaluate it numerically directly
*

*> from FeynHelpers (while developing the add-on my main focus were
*

*> symbolic evaluations). However, you can easily do something like
*

*>
*

*> exp = PaXDiLog[2.7, 1] /. PaXDiLog -> X`DiLog
*

*> Export["exp.m", exp]
*

*>
*

*> Quit[]
*

*> << X`
*

*> Import["exp.m"]
*

*>
*

*> to obtain the numerical value from Package-X. I'll contact the developer
*

*> of Package-X to see if we can find a better solution...
*

*>
*

*> I agree that PaXDiLog[Complex[-1,-6],-0.2] does not look correct at all.
*

*> Could you provide a minimal working code example that generates this
*

*> weird expression?
*

*>
*

*> I also wish you happy holidays.
*

*>
*

*> Cheers,
*

*> Vladyslav
*

*>
*

*> Am 25.12.2016 um 10:51 schrieb Xiu-Lei Ren:
*

*>> Dear Vladyslav,
*

*>>
*

*>> Thank you very much for your quick reply. It helps a lot.
*

*>>
*

*>> However, when i try to obtain the analytic expressions of triangle
*

*>> diagram mentioned in the previous email, I also encountered two
*

*>> questions about PaXDiLog.
*

*>>
*

*>> In order to avoid unexpected results when performing Dimension -> 4, I
*

*>> use the recommended FeynHelper--Package-X.
*

*>>
*

*>> When I do this, the treatment of pave coefficient C0 is necessary.
*

*>> In my case, (I am handling the two-nucleon scattering with two-pion
*

*>> exchange.
*

*>> mN, mpi deonte as nucleon and pion masses, p4 is the momentum of
*

*>> outgoing nucleon, q is the transfer momentum between two nucleons.)
*

*>>
*

*>> XC0 = C0[p4^2, q^2, (p4+q)^2, mN^2, mpi^2, mpi^2]
*

*>>
*

*>> should be replaced by using
*

*>>
*

*>> XC0Re = PaXEvaluate[XC0, PaXC0Expand -> True]//Normal
*

*>>
*

*>> Apparently, the output is lengthy with conditions.
*

*>>
*

*>> Then, perform the 1/mN expansion,
*

*>>
*

*>> Series[XC0Re, {mN, infty, 0}]//Normal
*

*>>
*

*>> The result always contains Li2 functions (PaXDiLog).
*

*>>
*

*>> 1) How one can transfer PaXDiLog to PolyLog?
*

*>>
*

*>> Furthermore, when I do the numerical evaluation for checking,
*

*>> I also find another problem about PaXDiLog.
*

*>>
*

*>> 2) e.g. PaXDiLog[Complex[-1,-6],-0.2], it cannot give a numerical value.
*

*>>
*

*>> Could you kindly let me know how to handle these problem?
*

*>>
*

*>> Merry Christmas and happy new year.
*

*>>
*

*>> Cheers,
*

*>> Xiu-Lei
*

*>>
*

*>
*

**Next message:**Jongping Hsu: "Re: Question about Subscript[E, 0](0,0,Overscript[p, _]^2,0,Overscript[p, _]^2,0,Overscript[p, _]^2,Overscript[p, _]^2,Overscript[p, _]^2,Overscript[p, _]^2,0,0,0,0,0)"**Previous message:**Vladyslav Shtabovenko: "Re: Question about TID in FC9.2.0"**Maybe in reply to:**Xiu-Lei Ren: "Question about TID in FC9.2.0"**Next in thread:**Jongping Hsu: "Re: Question about Subscript[E, 0](0,0,Overscript[p, _]^2,0,Overscript[p, _]^2,0,Overscript[p, _]^2,Overscript[p, _]^2,Overscript[p, _]^2,Overscript[p, _]^2,0,0,0,0,0)"**Reply:**Jongping Hsu: "Re: Question about Subscript[E, 0](0,0,Overscript[p, _]^2,0,Overscript[p, _]^2,0,Overscript[p, _]^2,Overscript[p, _]^2,Overscript[p, _]^2,Overscript[p, _]^2,0,0,0,0,0)"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ]**Mail actions:**[ respond to this message ] [ mail a new topic ]

*
This archive was generated by hypermail 2b29
: 02/18/18-12:20:01 AM Z CET
*