**Next message:**Vladyslav Shtabovenko: "Re: About ApartFF"**Previous message:**Jongping Hsu: "Fwd: 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)"**In reply to:**Jongping Hsu: "Fwd: 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)"**Next in thread:**Xiu-Lei Ren: "Re:Re:Re: Question about TID in FC9.2.0"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ]**Mail actions:**[ respond to this message ] [ mail a new topic ]

Hi Jongping,

*>> (Since your computer does not have these problem, I would like to ask
*

*>> you again to obtain the relation between E0 and B0 for me to check my
*

*>> result. Thank you.)
*

Sure:

PaVe[0, {0, 0, SPD[p, p], 0, SPD[p, p], 0, SPD[p, p], SPD[p, p],

SPD[p, p], SPD[p, p]}, {0, 0, 0, 0,

0}] == -(((-8 + D) (-5 + D) (-3 + D) B0[SPD[p, p], 0, 0])/(

2 SPD[p, p]^3))

I think that I also understand where your problems come from. In all the

notebooks you sent me I see the

FeynCalc is already loaded! To reload it, please restart the kernel.

$Aborted

message at the very beginning. This means that you tried to evaluate

those notebooks from a Mathematica session that already had FeynCalc loaded.

This is a problem, because add-ons like FeynHelpers can be activated

only when FeynCalc is loaded for the first time, but not at a later stage.

So if you have a notebook open, where FeynCalc was loaded alone,

like

<<FeynCalc`

and then you try to evaluate a notebook that begins with

$LoadAddOns={"FeynHelpers"};

<<FeynCalc`

then you will receive the "FeynCalc is already loaded!"-message and

FeynHelpers will not be loaded, so that you cannot FIREBurn and PaXEvaluate.

Hence, you should either open those notebooks in a new Mathematica

session, or restart the kernel via Evaluation -> Quit kernel -> Local

beforehand. Please make sure instead of

"FeynCalc is already loaded!" you see something like

Loading FeynCalc from /home/vs/.Mathematica/Applications/FeynCalc/

$PrePrint is set to FeynCalcForm. Use FI and FC to change the display

format.

FeynCalc 9.2.0. For help, use the documentation center, check out the

wiki or write to the mailing list.

See also the supplied examples. If you use FeynCalc in your research,

please cite

• V. Shtabovenko, R. Mertig and F. Orellana, Comput. Phys. Commun.,

207C, 432-444, 2016, arXiv:1601.01167

• R. Mertig, M. Böhm, and A. Denner, Comput. Phys. Commun., 64,

345-359, 1991.

FeynHelpers 1.0.0 loaded.

Have a look at the supplied examples. If you use FeynHelpers in your

research, please cite

• V. Shtabovenko, "FeynHelpers: Connecting FeynCalc to FIRE and

Package-X", TUM-EFT 75/15, arXiv:1611.06793

Furthermore, remember to cite the authors of the tools that you are

calling from FeynHelpers, which are

• FIRE by A. Smirnov, if you are using the function FIREBurn.

• Package-X by H. Patel, if you are using the function PaXEvaluate.

like in all the notebooks that I sent you. This way everything should

work properly.

Cheers,

Vladyslav

Am 27.12.2016 um 15:21 schrieb Jongping Hsu:

*> Sorry, I did not send the attachments last night. Here are the three
*

*> files (and the e-mail below) for your reference. Thanks. JP
*

*>
*

*> HSU Jongping,
*

*> Chancellor Professor
*

*> Department of Physics
*

*> Univ. of Massachusetts Dartmouth,
*

*> North Dartmouth, MA 02747. FAX (508)999-9115
*

*> http://www.umassd.edu/engineering/phy/people/facultyandstaff/jong-pinghsu/
*

*> recent monograph: Space-Time Symmetry and Quantum Yang–Mills Gravity
*

*> (https://sites.google.com/site/yangmillsgravity123/)
*

*>
*

*> ------------------------------------------------------------------------
*

*> *From: *"Jongping Hsu" <jhsu@umassd.edu>
*

*> *To: *"Vladyslav Shtabovenko" <dev@shtabovenko.df-kunde.de>
*

*> *Cc: *feyncalc@feyncalc.org
*

*> *Sent: *Monday, December 26, 2016 5:05:49 PM
*

*> *Subject: *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)
*

*>
*

*> Hi, Vladyslav,
*

*> Indeed, it is hard to believe those puzzling problems I told you in the
*

*> previous e-mails. I attached 3 nb files to let you see my puzzling
*

*> results. I have used both Mathematica 11 and 8 to run them and have the
*

*> same problems (shown in FIREBurn and "FCLoopSplit::fail").
*

*> 1. Your original file ( good for F0)
*

*> 2. The results in your original file changed after I ran it in my
*

*> computers (with Mathematica 11 and 8).
*

*> 3. I tried to follow your codes and failed to obtain the relation for E0.
*

*> Why?
*

*>
*

*> (Since your computer does not have these problem, I would like to ask
*

*> you again to obtain the relation between E0 and B0 for me to check my
*

*> result. Thank you.) JP
*

*>
*

*> HSU Jongping,
*

*> Chancellor Professor
*

*> Department of Physics
*

*> Univ. of Massachusetts Dartmouth,
*

*> North Dartmouth, MA 02747. FAX (508)999-9115
*

*> http://www.umassd.edu/engineering/phy/people/facultyandstaff/jong-pinghsu/
*

*> recent monograph: Space-Time Symmetry and Quantum Yang–Mills Gravity
*

*> (https://sites.google.com/site/yangmillsgravity123/)
*

*>
*

*> ------------------------------------------------------------------------
*

*> *From: *"Vladyslav Shtabovenko" <dev@shtabovenko.df-kunde.de>
*

*> *To: *feyncalc@feyncalc.org
*

*> *Sent: *Monday, December 26, 2016 3:10:58 PM
*

*> *Subject: *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)
*

*>
*

*> Hi Jongping,
*

*>
*

*> Am 26.12.2016 um 03:35 schrieb Jongping Hsu:
*

*>> Hi,Vladyslav,
*

*>> Happy Holidays!
*

*>
*

*> happy holidays to you too.
*

*>
*

*>> I need some more relations to check the consistency of the OneLoop result:
*

*>> What is the relation between B0(p^2,0,0) and
*

*>> F=E0(0,0,p^2,0,p^2,0,p^2,p^2,p^2,p^2,0,0,0,0,0) ?
*

*>
*

*> Just use the code that I sent you previously:
*

*>
*

*> https://feyncalc.org/forum/1156.html
*

*>
*

*> with
*

*>
*

*> FAD[{q, 0, 3}, {q + p, 0, 2}]/(I Pi^2)
*

*>
*

*> All your scalar integrals follow the same pattern, so
*

*> if you play a bit with the integers m and n in
*

*>
*

*> FAD[{q, 0, m}, {q + p, 0, n}]
*

*>
*

*> you can obtain all those relations for scalar integrals that you are
*

*> looking for.
*

*>
*

*>>
*

*>> When I used TID, e.g., "ampssPT = (TID[#, q, ToPaVe -> True,
*

*>> UsePaVeBasis -> True] & /@
*

*>> JPamp) ", I got almost 20 new and unfamiliar functions, E3,
*

*>> E4,...E444,..D002, D333, etc.
*

*>> Is there a source or reference that one can get the relation between
*

*>> them and, say, B(p^2,0,0)?
*

*>
*

*> The simplest way would be not to use the "ToPaVe" and "UsePaVeBasis"
*

*> options, so that you end up with scalar integrals that can be IBP
*

*> reduced. Otherwise I'm not aware of an existing symbolic codes to
*

*> evaluate higher order Passarino-Veltman functions, although this
*

*> probably should be using possible formulas from the work of Denner:
*

*>
*

*> https://arxiv.org/abs/0709.1075
*

*> https://arxiv.org/abs/hep-ph/0212259
*

*> https://arxiv.org/abs/hep-ph/0509141
*

*>
*

*> All those A,B,C,D,E,... PaVe functions are defined according to Denner.
*

*> The scalar integrals (those with 0 subscript) are quite simple, c.f. for
*

*> example the LoopTools manual
*

*> <http://www.feynarts.de/looptools/LT28Guide.pdf>
*

*>
*

*> and the appendix of the
*

*>
*

*> https://arxiv.org/pdf/1604.06792.pdf
*

*>
*

*> for the general formula.
*

*>
*

*> The coefficient functions (those with subscripts different from 0) are
*

*> more complicated. Depending on the kinematics it might be possible to
*

*> write them down only as Feynman parameter integrals, although others
*

*> might be reducible into scalar integrals.
*

*>
*

*> Cheers,
*

*> Vladyslav
*

*>
*

*>
*

*>> Thanks. JP
*

*>>
*

*>> HSU Jongping,
*

*>> Chancellor Professor
*

*>> Department of Physics
*

*>> Univ. of Massachusetts Dartmouth,
*

*>> North Dartmouth, MA 02747. FAX (508)999-9115
*

*>> http://www.umassd.edu/engineering/phy/people/facultyandstaff/jong-pinghsu/
*

*>> recent monograph: Space-Time Symmetry and Quantum Yang–Mills Gravity
*

*>> (https://sites.google.com/site/yangmillsgravity123/)
*

*>>
*

*>> ------------------------------------------------------------------------
*

*>> *From: *"Vladyslav Shtabovenko" <noreply@feyncalc.org>
*

*>> *To: *feyncalc@feyncalc.org
*

*>> *Sent: *Sunday, December 25, 2016 2:57:28 PM
*

*>> *Subject: *Re: Question about TID in FC9.2.0
*

*>>
*

*>>
*

*>>>> 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:**Vladyslav Shtabovenko: "Re: About ApartFF"**Previous message:**Jongping Hsu: "Fwd: 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)"**In reply to:**Jongping Hsu: "Fwd: 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)"**Next in thread:**Xiu-Lei Ren: "Re:Re:Re: Question about TID in FC9.2.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
: 06/24/17-07:00:06 AM Z CEST
*