Brief Tutorial On Feyncalc: Vladyslav Shtabovenko

Brief Tutorial on FeynCalc
based on
arXiv:1601.01167, arXiv:1611.06793
Vladyslav Shtabovenko
Technische Universität München, Germany
Tutorial at Methods of EFT

& Lattice Field Theory School
3 RD OF J ULY, 2017
Physik-Department T30f
V. Shtabovenko (TUM) @ EFT School, 3.07.2017 FeynCalc 1 / 14

Outline
1 Motivation And Brief History of FeynCalc
2 Using FeynCalc in your research
3 FeynHelpers: FeynCalc on Steroids
4 Summary and Outlook

Motivation And Brief History of FeynCalc
A generic perturbative QFT calculation There are many ways to automatize such
may involve many steps calculations:
I Feynman diagrams
Self-written codes
I Feynman rules from L
I Diagram generation Write private codes in
I Amplitudes FORM [Vermaseren, 2000], M ATHEMATICA,
I ... M APLE, R EDUCE, . . .
I Dirac algebra
I Simplification of γ-matrix chains Separate packages
I Dirac traces Automatize each step separately using
I Fierz identities standalone packages (e. g. F EYN A RTS,
I ... L OOP TOOLS [Hahn & Perez-Victoria, 1999],
I Loop integrals F EYN RULES [Christensen & Duhr, 2008],
QG RAF [Nogueira, 1993],
I Tensor reduction T RACER [Jamin & Lautenbacher, 1993],
I Partial fractioning F ORM T RACER [Cyrol et al., 2016],
I Mapping of topologies F ORCER [Ruijl et al., 2017],
I IBP-Reduction PY S EC D EC [Borowka et al., 2017], . . . ).
I Numerical evaluation
I ...
I Phase space integrals All-in-one packages
I ... Employ all-in-one packages that handle
most of these steps in one framework.

Two big categories of all-in-one packages

I Fully-automatic (F ORM C ALC [Hahn & Perez-Victoria, 1999], G O S AM [Cullen et al., 2014],
GRACE [Belanger et al., 2006], D IANA [Tentyukov & Fleischer, 2000], FDC [Wang, 2004],
C ALC HEP [Belyaev et al., 2012], C OMP HEP [Boos et al., 2004] . . . )
I Semi-automatic (F EYN C ALC [Mertig et al., 1991, Shtabovenko et al., 2016],
HEPM ATH [Wiebusch, 2014], PACKAGE -X [Patel, 2015], . . . )
Fully-automatic tools Semi-automatic tools

Blackbox: Require only minimal user Toolbox: Combine different tools with
input and provide a small set of options. many options to get the computation done
The code takes care of the rest. in your way.
Easy to use You must know what you are doing

Foolproof Easy to make mistakes
Constantly good performance The performance depends on your skills
Saves your time Writing codes may take quite some time
Limited number of templated calculations Very broad range of applications
Difficult to extend/modify for your needs Extendable with user-defined objects
Not easy to obtain intermediate results Intermediate results at each step

F EYN C ALC
I M ATHEMATICA package for symbolic
semi-automatic QFT calculations
I Open source (GPLv3) and publicly available
I Can process Feynman diagrams and
standalone QFT expressions
I Widely used in Effective Field Theory
calculations
A bit of history ...

1991 • F EYN C ALC 1.0 by R. Mertig [Mertig et al., 1991]
1991-1996 • New features and improvements
1997 • Release of TARCER [Mertig & Scharf, 1998]
1997-2000 • Many contributions from F. Orellana (e. g. PHI [Orellana, 2002] for ChPT)
2012 • F EYN C ALC F ORM L INK [Feng & Mertig, 2012]
2016 • F EYN C ALC 9.0 [Shtabovenko et al., 2016]
2017 • F EYN H ELPERS [Shtabovenko, 2017]

Current status of the project

I (Again) under active development since 2014
I Regular releases since 2016: F EYN C ALC 9.0.0,
F EYN C ALC 9.0.1, F EYN C ALC 9.1.0,
F EYN C ALC 9.2.0
I Easy installation
I Online documentation: https:
//feyncalc.github.io/reference
I Many sample tree and 1-loop calculations (QED, EW, QCD) included
I Public source code repository: https://www.github.com/feyncalc
I hotfix-stable branch: stable version + bugfixes
I master branch: development version (new features, less stable)
I Extensive unit testing framework (based on MUnit, over 4000 tests)
F EYN C ALC developer team

I Rolf Mertig (GluonVision GmbH): original author of the package, first release 1991
I Frederik Orellana (Technical University of Denmark): joined 1997
I VS (Technical University of Munich, soon Zhejiang University (China)): joined 2014

Using FeynCalc in your research General picture
What can F EYN C ALC do?

I Standard feature set: Lorentz algebra, Dirac algebra, Color algebra
I Tensor reduction of 1-loop integrals (Passarino–Veltman)
I Basic support for manipulating multi-loop integrals
I Built-in interface to F EYN A RTS [Hahn & Perez-Victoria, 1999] (also works with custom
F EYN A RTS model from F EYN RULES [Christensen & Duhr, 2008])
I Optional interfaces to other HEP tools
I Extensive typesetting for better readability
I Extendable with self-written M ATHEMATICA codes
When is F EYN C ALC useful? What are the limitations of F EYN C ALC?
I Small or medium-sized I Cannot be used without
calculations, too specific for fully M ATHEMATICA
automatic packages I Inherits M ATHEMATICA’s
I F EYN C ALC as a “calculator” for performance problems with large
QFT expressions number of terms
I Cross-check results from other I Not really suited for very large and
people complex calculations
I Extensive manipulations on the I Much slower than FORM
level of the amplitudes I Only algebraic manipulations.

Using FeynCalc in your research Brief overview of selected functions
Lorentz algebra (most used functions)

I ChangeDimension: p̄ · q̄ → p̂ · q̂
I Contract: pµ qµ → p · q
I ExpandScalarProduct: a · (b + c) → a · b + a · c
I FourDivergence: ∂µ
I FourLaplacian: ∂µ ∂ µ
I MomentumCombine: pµ + qµ → (p + q)µ
I Uncontract: p · q → pµ qµ
New functions
I FCRenameDummyIndices (v. 9.0): pµ qµ → pα qα
I FCCanonicalizeDummyIndices (v. 9.1): pµ qµ + pν qν → 2pα qα
I p − /̄
FCGetDimensions (v. 9.2): / p → {D, D − 4}

Dirac algebra (most used functions)

I Chisholm: γ̄ µ γ̄ ν γ̄ λ → ḡµν γ̄ λ + ḡνλ γ̄ µ − ḡµλ γ̄ ν + iε̄µνλσ γ¯σ γ 5
I pu(p) → mu(p)
DiracEquation: /
I DiracGammaCombine: /p+/q → γ µ (p + q)µ
I DiracGammaExpand: γ µ (p + q)µ → /
p+/
q
ν µ µν µ ν
I DiracOrder: γ γ → 2g − γ γ
I DiracReduce: γ̄ µ γ̄ ν → ḡµν − iσ µν
I DiracSigmaExplicit: σ µν → 2i [γ µ , γ ν ]
I DiracTrace: Tr(γ µ γ ν ) → 4gµν
I DiracSimplify: Applies all available simplifications
New functions
γ µ ( p +m)γ e2
I FCDiracIsolate (v. 9.1): e2 (q2 −m/2 )(q2 −pµ2 ) → (q2 −m2 )(q2 −p2 )
[γ µ (/
p + m)γµ ]
I Most changes in the Dirac algebra related functions happened under the hood
(performance improvements, cleaner code . . . )

Loop calculations (most used functions)

1 1
I FeynAmpDenominatorSimplify: (q+p)2 −m2
→ q2 −m2
qµ qν m2 gµν
I TID: q2 −m2
→ D q2 −m2
m4
I PaVeReduce: A0000 (m2 ) → A (m2 )
D(D+2) 0
New functions
I ToPaVe (v. 9.0): 1
q2 −m2
→ iπ 2 A0 (m2 )
q/1 q/2 1 m2 1
I FCMultiLoopTID (v. 9.0): q21 q22 [(q1 −q2 )2 −m2 ]
→ q21 [(q1 −q2 )2 −m2 ]
− 2 q2 q2 [(q1 −q2 )2 −m2 ]
1 2
1 1 1
I ApartFF (v. 9.0): q2 (q−p)2 (q+p)2
→ p2 q2 (q−p)2
− p2 q2 (q−2p)2
I FCLoopBasisIncompleteQ (v. 9.0): Do the propagators of the given loop
integral form a basis?
I FCLoopBasisOverDeterminedQ (v. 9.0): Is the propagator basis
overdetermined?
I FCLoopBasisFindCompletion (v. 9.0): Which propagators need to be added to
have the complete basis?

FeynHelpers: FeynCalc on Steroids An example from physics of heavy quarkonia
We have to go deeper ...

I Often, the final F EYN C ALC result is not the most useful one
I No automatic evaluation of Passarino–Veltman functions
I No automatic IBP-reduction of (multi-)loop integrals a
I A lot of good HEP tools are publicly available, no need to reinvent the wheel
I However, combining them with F EYN C ALC is not always easy
I Well-tested interfaces supported by the developers are better than private hacks
a
TARCER [Mertig & Scharf, 1998] is limited to 2-loop propagator-type integrals
The first step into the right direction: F EYN H ELPERS [Shtabovenko, 2017]
I Seamless integration of PACKAGE -X [Smirnov, 2015] and FIRE [Patel, 2015] into
F EYN C ALC
I Access the PACKAGE -X library of analytic results for 1-loop integrals and the
IBP-reduction engine of FIRE directly from a F EYN C ALC session
I All conversions (syntax, normalization) are automatic and require no user
intervention
I Open-source interface code: https://github.com/FeynCalc/feynhelpers
I F EYN H ELPERS comes with many neat examples from QED, QCD, EW, EFTs etc.

Consider
lµ lν
Z
I µν = , with P2 = 4m2Q
l l2 [(l + P/2)2 − mQ ][(l − P/2)2 − mQ ]2
F EYN C ALC can handle this integral, but the result is not particularly enlightening
In fact, one can show that

(D − 4)(D − 2) Z 1
I µν = (D − 6)P µ ν
P + 4m 2 µν
Q g .
64(D − 5)(D − 3)m6Q 2
l l − mQ
2
Is it possible to get this using F EYN C ALC in an easy way?

With F EYN C ALC and F EYN H ELPERS such calculations become almost trivial...
The performance is, however, still an issue and will be improved in the subsequent
versions (by switching to the C++ backend of FIRE)
Summary and Outlook
I Modern HEP theory and phenomenology heavily rely on programs for automatic
calculations.
I The market for QFT software is growing rapidly, new tools appear almost monthly.
I Different projects require different tools and F EYN C ALC is one of them.
I The development of F EYN C ALC focuses on versatility, flexibility, convenience and
performance.
I F EYN C ALC is particularly useful for “nonstandard” calculations that are not covered
by other (more automatic) tools.
I In the last 3 years we tried to achieve several goals:
I Refactor the source code, fix old (and new) bugs, make things faster
I Facilitate manipulations of tree and loop amplitudes by introducing new functions and
improving the existing ones
I Interface F EYN C ALC with other tools useful for HEP calculations
I Yet, much still remains to be done:
I There is still a lot of room for improvements regarding stability, the number of available
function and their performance.
I We need more interfaces to other HEP packages, e. g. to L ITE R ED [Lee, 2012], S@M
[Maitre & Mastrolia, 2008], F ORM T RACER [Cyrol et al., 2016], QG RAF [Nogueira, 1993].

Backup
Tensor reduction a la Passarino–Veltman

I Fairly old technique for dealing with tensor 1-loop integrals [Passarino & Veltman, 1979]
I Still widely used in many loop calculations.
I Main idea: convert all the tensor integrals into scalar ones (Passarino–Veltman
coefficient functions)
I Evaluation of any 1-loop integral can be reduced to the evaluation of the resulting
coefficient functions
I A lot of tools for numerical evaluation: FF [van Oldenborgh, 1991],
L OOP TOOLS [Hahn & Perez-Victoria, 1999], QCDL OOP [Carrazza et al., 2016],
O NE LO OP [van Hameren, 2011], G OLEM 95C [Cullen et al., 2011],
PJF RY [Fleischer & Riemann, 2011], C OLLIER [Denner et al., 2017], . . .
I Where to get analytic results for singular kinematics or zero Gram determinants?
Often needed for renormalization, EFTs, . . .
I Most of the results can be found somewhere in the literature.

Backup
PACKAGE -X
I Recent [Patel, 2015] M ATHEMATICA package for semi-automatic 1-loop calculations
(closed-source freeware)
I Unique feature: Library of analytic expressions for Passarino–Veltman functions
with up to 4 legs and almost arbitrary kinematics.
I Can also extract UV- and IR-parts and expand coefficient functions in their
arguments.
I Someone indeed has collected all those results from the literature!
Interface to PACKAGE -X
I Main function: PaXEvaluate
I Works: on scalar 1-loop integrals (unit numerators) and Passarino–Veltman
coefficient functions A, B, C and D
I Takes two arguments (plus options): input expression, loop momentum.
I Use PaXEvaluateUV(PaXEvaluateIR) to get the UV(IR)-divergent part of the
result
I PaXEvaluateUVIRSplit returns the full result with the explicit distinction
between UV and IR .
I All four functions share the same set of options

Backup
dD q 1
R
Let us compute (2π)D q2 −m2
In[1]:= →1/(2Pi)^D]
int=PaXEvaluate[FAD[{q,m}],q,PaXImplicitPrefactor→
µ2
im2 −log

im2 +γ−1−log(4π)
m2
Out[1]= −
16π 2 ε 16 π 2
Make the result look more compact (∆ ≡ 1/ − γE + log(4π)) using FCHideEpsilon
In[2]:= int//FCHideEpsilon
2 µ2
i∆m2 im log m2 +1
Out[2]= +
16π 2 16π 2
Evaluation of Passarino–Veltman functions:
In[3]:= PaXEvaluate[B0[SPD[p,p],0,m^2]]
2 m2 log m2

m2

1 µ m2 −p2
Out[3]= +log 2
− 2
+log 2 2
−γ+2
ε πm p m −p

Backup
We can also expand coefficient functions in their parameters (masses or external

momenta). To expand B0 (p2 , 0, m2 ) around p2 = m2 up to first order with PaXEvaluate
we first need to assign an arbitrary symbolic value to the scalar product p2 , e.g. pp
In[4]:= SPD[p,p]=pp;
Then use the option PaXSeries to specify the expansion parameters and activate the
option PaXAnalytic
In[5]:= →{{pp,m^2,1}},PaXAnalytic→
PaXEvaluate[B0[SPD[p,p],0,m^2],PaXSeries→ →True]
2 µ2
3 m2 −pp (3 m −pp) −log m2 +γ−2+log(π)
Out[5]= −
2εm2 2m2
Get only in the UV-part of this series: PaXEvaluate with PaXEvaluateUV
PaXEvaluateUV[B0[SPD[p,p],0,m^2],PaXSeries→ →True]
1
Out[6]=
εUV
The IR-part is equally easy
PaXEvaluateIR[B0[SPD[p,p],0,m^2],PaXSeries→ →True]
m2 −pp
Out[7]=
2m2 εIR
Full result with the explicit distinction between UV and IR singularities
PaXEvaluateUVIRSplit[B0[SPD[p,p],0,m^2],PaXSeries→ →True]
2 µ2
m2 −pp (3 m −pp) −log m2 +γ−2+log(π) 1
Out[8]= − +
2m2 εIR 2 m2 εUV

Backup
Integration-by-parts identities
I Reduction of scalar loop integrals using integration-by-parts (IBP) identities
[Chetyrkin & Tkachov, 1981] is a standard technique in modern loop calculations.
I Many publicly available IBP-packages on the market: FIRE [Smirnov & Smirnov, 2013],
L ITE RED [Lee, 2012], R EDUZE [Studerus, 2009], AIR [Anastasiou & Lazopoulos, 2004], . . .
I Expected input: loop integrals with propagators that form a basis.
I What about integrals with an incomplete or overdetermined basis?
I FCLoopBasisIncompleteQ detects integrals that require a basis completion
I FCLoopBasisFindCompletion gives a list of propagators (with zero exponents)
required to complete the basis
I FCLoopBasisOverdeterminedQ checks if the propagators are linearly dependent.
Such integrals can be decomposed further using ApartFF.

Backup
Interface to FIRE
I Main function: FIREBurn
I Reduces scalar multi-loop integrals to simpler ones using IBP-techniques.
I Takes three arguments (plus options): input expression, list of loop momenta and
the list of external momenta.
I Automatically adds propagators to integrals with incomplete bases of propagators
I Automatically detects integrals with linearly dependent propagators
Current limitations
I No recognition of integral families
I Each loop integral is evaluated separately
I Hence, rather inefficient ...

Backup
dD l
R
IBP-reduce the 1-loop integral [l2 ]2 [(l−p)2 −m2 ]2
In[9]:= FIREBurn[FAD[{l,0,2},{l−p,m,2}],{l},{p}]
(D−2)(2 D m2 −9 m2 −pp) (D−3)(D m2 +D pp−4m2 −6 pp)

Out[9]=
3
− 3
2 2 2 2
2m (m −pp) ((l−p) −m ) (m2 −pp) l2 .((l−p)2 −m2 )
No dependence on external momenta → supply an empty list for the third argument. For
D D
q2 d D q3
example, for [q2 −m2 ]2 [(q d+qq1)d2 −m
R
2 ][(q −q )2 ][q2 ]2
1 1 3 2 3 2
In[10]:= FIREBurn[FAD[{q1,m,2},{q1+q3,m},{q2−q3},{q2,0,2}],{q1,q2,q3},{}]
(D−3)(3 D−10)(3 D−8)
Out[10]= −
16(2 D−7)m4 (q12 −m2 ).q22 .(q2−q3)2 .((q1+q3)2 −m2 )

Backup
Check out the included examples for more nice calculations

I ABJ axial anomaly (following Chapter 19.2 in Peskin and Schroeder)
QEDABJAxialAnomaly.m
I Electron’s g − 2 in QED at 1-loop
QEDElectronGMinusTwoOneLoop.m
I 1-loop renormalization of QED in MS, MS and OS schemes (full gauge dependence)
QEDRenormalizationOneLoop.m
I 1-loop gluon self-energy (full gauge dependence) in QCD
QCDQuarkSelfEnergyOneLoop.m
I 1-loop quark self-energy (full gauge dependence) in QCD
QCDGluonSelfEnergyOneLoop.m
I LO matching between QED and Euler-Heisenberg’s EFT (development version)
QEDToEulerHeisenbergLagrangianMatching.m
I UV part of the 3-gluon vertex at 1-loop in QCD (development version)
QCDThreeGluonVertexOneLoop.m
I H → gg decay via a top-quark loop
EWHiggsToTwoGluonsOneLoop.m

Backup
Anastasiou, C. & Lazopoulos, A. (2004).

Automatic Integral Reduction for Higher Order Perturbative Calculations.
JHEP, 0407, 046.
Belanger, G., Boudjema, F., Fujimoto, J., Ishikawa, T., Kaneko, T., Kato, K., & Shimizu, Y. (2006).
GRACE at ONE-LOOP: Automatic calculation of 1-loop diagrams in the electroweak theory with gauge
parameter independence checks.
Phys. Rept., 430, 117–209.
Belyaev, A., Christensen, N. D., & Pukhov, A. (2012).

CalcHEP 3.4 for collider physics within and beyond the Standard Model.
Boos, E., Bunichev, V., Dubinin, M., Dudko, L., Ilyin, V., Kryukov, A., Edneral, V., Savrin, V., Semenov, A., &
Sherstnev, A. (2004).
CompHEP 4.4: Automatic computations from Lagrangians to events.
Nucl. Instrum. Meth., A534, 250–259.
Borowka, S., Heinrich, G., Jahn, S., Jones, S. P., Kerner, M., Schlenk, J., & Zirke, T. (2017).
pySecDec: a toolbox for the numerical evaluation of multi-scale integrals.
Carrazza, S., Ellis, R. K., & Zanderighi, G. (2016).

QCDLoop: a comprehensive framework for one-loop scalar integrals.
Comput. Phys. Commun., 209, 134–143.
Chetyrkin, K. & Tkachov, F. (1981).

Integration by parts: The algorithm to calculate β-functions in 4 loops.
Nucl. Phys. B, 192(1), 159–204.
Backup
Christensen, N. D. & Duhr, C. (2008).

FeynRules - Feynman rules made easy.
Cullen, G., Guillet, J. P., Heinrich, G., Kleinschmidt, T., Pilon, E., Reiter, T., & Rodgers, M. (2011).
Golem95C: A library for one-loop integrals with complex masses.
Cullen, G., van Deurzen, H., Greiner, N., Heinrich, G., Luisoni, G., Mastrolia, P., Mirabella, E., Ossola, G.,
Peraro, T., Schlenk, J., von Soden-Fraunhofen, J. F., & Tramontano, F. (2014).
GoSam-2.0: a tool for automated one-loop calculations within the Standard Model and beyond.
Eur. Phys. J. C, 74, 8, 3001.
Cyrol, A. K., Mitter, M., & Strodthoff, N. (2016).

FormTracer - A Mathematica Tracing Package Using FORM.
Denner, A., Dittmaier, S., & Hofer, L. (2017).

Collier: a fortran-based Complex One-Loop LIbrary in Extended Regularizations.
Feng, F. & Mertig, R. (2012).

FormLink/FeynCalcFormLink : Embedding FORM in Mathematica and FeynCalc.
Fleischer, J. & Riemann, T. (2011).

A complete algebraic reduction of one-loop tensor Feynman integrals.
Phys. Rev. D, 83, 073004.
Backup
Hahn, T. & Perez-Victoria, M. (1999).

Automatized One-Loop Calculations in 4 and D dimensions.
Jamin, M. & Lautenbacher, M. E. (1993).

TRACER version 1.1.
Comput. Phys. Commun., 74(2), 265–288.
Lee, R. N. (2012).
Presenting LiteRed: a tool for the Loop InTEgrals REDuction.
Maitre, D. & Mastrolia, P. (2008).

S@M, a Mathematica Implementation of the Spinor-Helicity Formalism.
Mertig, R., Böhm, M., & Denner, A. (1991).

Feyn Calc - Computer-algebraic calculation of Feynman amplitudes.
Mertig, R. & Scharf, R. (1998).

TARCER - A Mathematica program for the reduction of two-loop propagator integrals.
Nogueira, P. (1993).
Automatic Feynman graph generation.
J. Comput. Phys., 105, 279–289.
Backup
Orellana, F. (2002).
Mesonic Final State Interactions.
PhD thesis, University of Zurich.
Passarino, G. & Veltman, M. (1979).

One Loop Corrections for e+ e− Annihilation Into µ+ µ− in the Weinberg Model.
Nucl. Phys., B160, 151.
Patel, H. H. (2015).
Package-X: A Mathematica package for the analytic calculation of one-loop integrals.
Ruijl, B., Ueda, T., & Vermaseren, J. A. M. (2017).

Forcer, a FORM program for the parametric reduction of four-loop massless propagator diagrams.
Shtabovenko, V. (2017).
FeynHelpers: Connecting FeynCalc to FIRE and Package-X.
Shtabovenko, V., Mertig, R., & Orellana, F. (2016).

New Developments in FeynCalc 9.0.
Smirnov, A. V. (2015).
FIRE5: a C++ implementation of Feynman Integral REduction.
Backup
Smirnov, A. V. & Smirnov, V. A. (2013).

FIRE4, LiteRed and accompanying tools to solve integration by parts relations.
Studerus, C. (2009).
Reduze - Feynman Integral Reduction in C++.
Tentyukov, M. & Fleischer, J. (2000).

A Feynman Diagram Analyser DIANA.
van Hameren, A. (2011).

OneLOop: for the evaluation of one-loop scalar functions.
van Oldenborgh, G. (1991).

FF - a package to evaluate one-loop Feynman diagrams.
Vermaseren, J. A. M. (2000).
New features of FORM.
Wang, J.-X. (2004).

Progress in FDC project.
Nucl. Instrum. Meth. A, 534(1-2), 241–245.
Backup
Wiebusch, M. (2014).
HEPMath: A Mathematica Package for Semi-Automatic Computations in High Energy Physics.

Brief Tutorial On Feyncalc: Vladyslav Shtabovenko

Uploaded by

Copyright:

Available Formats

Brief Tutorial On Feyncalc: Vladyslav Shtabovenko

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Brief Tutorial On Feyncalc: Vladyslav Shtabovenko

Uploaded by

Copyright:

Available Formats

Brief Tutorial on FeynCalc

Technische Universität München, Germany

Tutorial at Methods of EFT

V. Shtabovenko (TUM) @ EFT School, 3.07.2017 FeynCalc 1 / 14

1 Motivation And Brief History of FeynCalc

2 Using FeynCalc in your research

3 FeynHelpers: FeynCalc on Steroids

4 Summary and Outlook

V. Shtabovenko (TUM) @ EFT School, 3.07.2017 FeynCalc 2 / 14

V. Shtabovenko (TUM) @ EFT School, 3.07.2017 FeynCalc 3 / 14

Two big categories of all-in-one packages

Fully-automatic tools Semi-automatic tools

Easy to use You must know what you are doing

V. Shtabovenko (TUM) @ EFT School, 3.07.2017 FeynCalc 4 / 14

A bit of history ...

V. Shtabovenko (TUM) @ EFT School, 3.07.2017 FeynCalc 5 / 14

Current status of the project

F EYN C ALC developer team

V. Shtabovenko (TUM) @ EFT School, 3.07.2017 FeynCalc 6 / 14

What can F EYN C ALC do?

V. Shtabovenko (TUM) @ EFT School, 3.07.2017 FeynCalc 7 / 14

Lorentz algebra (most used functions)

V. Shtabovenko (TUM) @ EFT School, 3.07.2017 FeynCalc 8 / 14

Dirac algebra (most used functions)

V. Shtabovenko (TUM) @ EFT School, 3.07.2017 FeynCalc 9 / 14

Loop calculations (most used functions)

V. Shtabovenko (TUM) @ EFT School, 3.07.2017 FeynCalc 10 / 14

We have to go deeper ...

V. Shtabovenko (TUM) @ EFT School, 3.07.2017 FeynCalc 11 / 14

In fact, one can show that

Is it possible to get this using F EYN C ALC in an easy way?

V. Shtabovenko (TUM) @ EFT School, 3.07.2017 FeynCalc 14 / 14

Tensor reduction a la Passarino–Veltman

V. Shtabovenko (TUM) @ EFT School, 3.07.2017 FeynCalc 15 / 14

V. Shtabovenko (TUM) @ EFT School, 3.07.2017 FeynCalc 16 / 14

V. Shtabovenko (TUM) @ EFT School, 3.07.2017 FeynCalc 17 / 14

We can also expand coefficient functions in their parameters (masses or external

V. Shtabovenko (TUM) @ EFT School, 3.07.2017 FeynCalc 18 / 14

V. Shtabovenko (TUM) @ EFT School, 3.07.2017 FeynCalc 19 / 14

V. Shtabovenko (TUM) @ EFT School, 3.07.2017 FeynCalc 20 / 14

(D−2)(2 D m2 −9 m2 −pp) (D−3)(D m2 +D pp−4m2 −6 pp)

V. Shtabovenko (TUM) @ EFT School, 3.07.2017 FeynCalc 21 / 14

Check out the included examples for more nice calculations

V. Shtabovenko (TUM) @ EFT School, 3.07.2017 FeynCalc 22 / 14

Anastasiou, C. & Lazopoulos, A. (2004).

Belyaev, A., Christensen, N. D., & Pukhov, A. (2012).

Carrazza, S., Ellis, R. K., & Zanderighi, G. (2016).

Chetyrkin, K. & Tkachov, F. (1981).

Christensen, N. D. & Duhr, C. (2008).

Cyrol, A. K., Mitter, M., & Strodthoff, N. (2016).

Denner, A., Dittmaier, S., & Hofer, L. (2017).

Feng, F. & Mertig, R. (2012).

Fleischer, J. & Riemann, T. (2011).

Hahn, T. & Perez-Victoria, M. (1999).

Jamin, M. & Lautenbacher, M. E. (1993).

Maitre, D. & Mastrolia, P. (2008).

Mertig, R., Böhm, M., & Denner, A. (1991).

Mertig, R. & Scharf, R. (1998).

Passarino, G. & Veltman, M. (1979).