Package-X is hosted by Hepforge, IPPP Durham

Hiren H. Patel

Version 2.0 of Package-X is now available for download (link provided below).


Package-X is a Mathematica package to obtain analytic expressions for one-loop integrals in higher order calculations of relativistic quantum field theory in 4 spacetime dimensions.

               Technical features
  • Computes dimensionally regulated one-loop integrals with up to four distinct propagators of arbitrarily high rank.
  • Calculates traces of Dirac matrices in d dimensions for closed fermion loops, or carries out Dirac algebra for open fermion lines
  • Generates analytic results for any kinematic configuration (e.g. at zero external momentum, at physical threshold, etc.) for real masses and external invariants.
  • Provides analytic expressions for UV-divergent, IR-divergent and finite parts either separately or all together.
  • Constructs Taylor series expansions for one-loop integrals
  • Evaluates with either machine precision or arbitrary precision arithmetic, and is always consistent with the + prescription
               User interface features
  • Complete documentation available within the Wolfram Documentation Center which includes
    • over 300 basic usage examples with detailed information for every symbol defined by Package-X,
    • 3 project-scale tutorials,
    • and instructions on linking Package-X to FeynCalc and LoopTools.
  • Readable and intuitive input (see Featured code snippets below)
  • Output easily manipulatable and readily numerically evaluatable

Featured code snippets

Compute the transverse part of the QED vacuum polarization tensor:
Its series expansion near q2 = 0:
Compute the QED vertex function:
Compute the UV-divergent part of the QCD vacuum polarization function in general covariant Rξ gauges:
Compute the massless one loop scalar four-point integral in 6 dimensions:


Package-X was developed after several long years of tedious work. Please include the following citation if Package-X was useful for your publication.
    Hiren H. Patel, Comput. Phys. Commun. 197, 276 (2015),
    ePrint: arXiv:1503.01469


Package-X is closed source and distributed under the standard CPC license agreement.

Latest release

v2.0.1: (October 25, 2016) Download
  • Fixed a bug where FermionLineExpand and other tensor manipulation routines would drop the quadruply contracted Levi-Civita symbol (thanks, Matthew Kirk).
  • Fixed a bug where FermionLineProduct would fail to automatically expand over deeply nested sums over FermionLine objects.
  • Fixed a bug where on rare occasions Series would incorrectly expand DiscB around normal threshold.
v2.0.0: (October 14, 2016) Download

Backwards compatibility

The following changes from v1.0.4 to v2.0.0 may affect existing notebooks that were created with Package-X v1.0:
  • Passarino-Veltman functions pvA, pvB, pvC, pvC0, and pvC0IR6 are now named PVA, PVB, PVC, ScalarC0, and ScalarC0IR6, respectively, consistent with the naming convention of capitalizing the first letter of pre-defined symbols.
  • LoopIntegrate now normalizes its integration measure so that e-γ ϵ is factored out instead of rΓ. This is to prepare Package-X for computing two loop integrals in the future. This only changes the output of LoopRefine for integrals exhibiting 1/ϵ2 poles, arising from overlapping soft and collinear divergences, by an amount proportional to ~π2/12 to the finite part.
  • Order of arguments of PVC (was pvC) and ScalarC0 (was pvC0) has been changed to match that of FeynCalc, LoopTools, and essentially all other authors in the literature (Sorry to v1.0 users who struggled to match against literature!). The relations are PVC[r,n1,n2,s1,s12,s2,m0,m1,m2] = pvC[r,n1,n2,s1,s2,s12,m2,m1,m0] and ScalarC0[s1,s12,s2,m0,m1,m2] = pvC0[s1,s2,s12,m2,m1,m0]
  • ScalarC0 is now a function defined manifestly in four dimensions, and will return ComplexInfinity for cases that are IR divergent. The related function PVC[0,0,0, args] lives in d dimensions, and will yield 1/ϵ poles by LoopRefine for IR divergent cases.
  • Package contexts X`IndexAlg`, X`Spur`, and X`OneLoop` are deprecated; all package symbols now belong to a common context X`
  • Auxiliary Passarino-Veltman function pvb is now obsolete, and is covered by higher weight PVB functions.
  • Dimensional regularization ('t Hooft) scale μR is now denoted with µ (\[Micro]), and entered with the keyboard alias [esc][m][m][esc].
LoopIntegrate generates capitalized PVA...PVD, and the three-point functions PVC and ScalarC0 are generated with the updated ordering for their arguments. LoopRefine generates expressions with µ denoting the dimensional regularization scale, and with the updated normalization affecting integrals with overlapping soft and collinear IR divergences. Since legacy syntax of LoopIntegrate is still possible and the syntax of LoopRefine is unchanged, rerunning code with v2.0 in notebooks that used v1.0 should update expressions reflecting the changes. You may download v1.0.4 here, but is no longer supported.


Any computer running Mathematica
Minimum: Mathematica 8.0.4 or higher.
Recommended: Mathematica 10.1.0 or higher.

Installation Instructions

Uncompress downloaded file, and copy the Package-X folder (/X) with all its contents to the folder $UserBaseDirectory/Applications/. See the tutorial for details.


Get started with the tutorial (October 14, 2016).