The infrared divergences of QCD have always rendered computations in perturbation theory dicult, making integrating the corresponding regions of phase-space delicate. In the current eorts to bring the state-of-the-art computations in QCD to N3LO, it is then crucial to develop a subtraction scheme at that order. Such a project requires the knowledge of the soft and collinear behaviour of QCD amplitudes, which is made easier by the universal factorisations that happen in these limits. For squared amplitudes in the collinear limit, the universal pieces are known as splitting functions. In this thesis, we compute the quadruple collinear splitting functions, adding an essential piece to the eort of the community to build a subtraction scheme at N3LO in QCD. We review how to compute such functions, and detail their symmetries as well as their colour structure. In addition, we check the soft and collinear limits of the splitting functions as both a way to verify our results and to explore these specific regions of phase-space.