  We use the antenna subtraction method to isolate the double real radiation
infrared singularities present in gluonic scattering amplitudes at
next-to-next-to-leading order. The antenna subtraction framework has been
successfully applied to the calculation of NNLO corrections to the 3-jet cross
section and related event shape distributions in electron-positron
annihilation. Here we consider processes with two coloured particles in the
initial state, and in particular two-jet production at hadron colliders such as
the Large Hadron Collider (LHC). We construct a subtraction term that describes
the single and double unresolved contributions from the six-gluon tree-level
process using antenna functions with initial state partons and show numerically
that the subtraction term correctly approximates the matrix elements in the
various single and double unresolved configurations.
