This is a generalisation of Electrodynamics. Here, the gauge fields (gluons) are matrices and do not commute, i.e., they will be self-interacting, unlike photons.
,
where are traceless Hermitian matrices corresponding to generators of .
- corresponds to QCD, describing strong interactions.
and .
Now, with .
.
Now, the gauge symmetry can be generalised to where .
Note: For , e.g., , we recover the gauge symmetry from electromagnetism since is just a function.
The action is invariant under the above symmetry. Further,
.
Under a gauge transformation,
,
, using the cyclicity of trace and the fact that is unitary.
We can fix the gauge using the Fadeev-Popov procedure:
, where is the gauge-fixing functional and we will neglect the ghosts since we are interested only in the tree-amplitudes.
We will choose the Gervais-Neveu gauge, which is convenient for amplitudes,
.
The Lagrangian then becomes .