ColorMath

ColorMath is a Mathematica package for symbolically performing color summed calculations in SU(N_c). It is based on advanced pattern matching and uses a syntax which is very similar to how QCD color structure is written on paper.

For a quick start, it is recommended to use the Mathematica notebook ColorMathTutorial1.0.nb. However, the ColorMath paper is recommended for more details and better understanding. ColorMath has been tested with Mathematica 7, 8 and 9, but also seems to work with version 10 and 11. ColorMath is published in EPJC as well as on the Wolfram Library Archive.

Downloads

The ColorMath paper, package and tutorial can be downloaded here:

Color bases

For many applications, it is desirable to use color bases for treating the color structure. Below is a selection of multiplet bases as Mathematica .m packages, which can be read in to Mathematica and used with ColorMath. The files are also fairly human readable. More information about the bases can be found by reading the top sections in the respective .m-files.

Gluon bases

Below are the orthogonal normalized multiplet color bases needed for processes with four, five and six gluons, for any N_c. They come in two versions, one charge conjugation invariant version, and one general version. The invariant bases use that in any pure gluon amplitude, the color structures related by swapping quarks and antiquarks always appear together. These bases are thus sufficient for gluon amplitudes. The larger, non-invariant, versions can be used for building bases with quarks (a gluon can be traded for a quark-antiquark pair) and for more straight forward calculations using Wigner 6j coefficients. Smaller N_c=3 versions of the bases are obtained by crossing out the vectors which do not appear for N_c=3. The bases are valid for any order in perturbation theory.

CMCNMBq0g4.m: Charge conjugation invariant four gluon basis (published in 1503.00530).
CMCNMBq0g5.m: Charge conjugation invariant five gluon basis (published in 1503.00530).
CMCNMBq0g6.m: Charge conjugation invariant six gluon basis (published in 1503.00530).
CMNMBq0g4.m: Four gluon basis.
CMNMBq0g5.m: Five gluon basis (published in 1507.03814).
CMNMBq0g6.m: Six gluon basis (published in 1507.03814).

Bases with quarks

Will be added shortly... rsn

Wigner 6j coefficients

ColorMath is based on rewriting color structures into small representations, and then contracting indices to evaluate invariants. An alternative way of performing calculations in color space is to use Wigner 6j and 3j coefficients, however, only some Wigner 6j coefficients have been calculated, and ColorMath cannot (currently) perform contractions using Wigner coefficients. A set of Wigner 6j coefficients, sufficient up to six gluon amplitudes up to NLO or tree level amplitudes with up to four incoming plus four outgoing color lines is supplied below. The 3j coefficients, which correspond to vertex normalizations, are normalized to one.

Wigner6jCoefficientsWithQuarks.m: The set of Wigner 6j coefficients published in 1809.05002.
UsingWignerCoefficients.nb: A Mathematica notebook illustrating how to access the 6j coefficients.