Loop-level gluon OPEs in celestial holography

We compute one-loop corrections to the OPE of gluons in the celestial conformal field theory corresponding to Yang-Mills coupled to arbitrary matter. We exploit universal hard/soft factorization to derive an IR finite OPE for the hard gluon operators. This OPE involves logarithms and operators that resemble logarithmic partners of primary operators. We derive an exact all-loop OPE in a limit of the Higgs-regulated planar $\mathcal{N}=4$ super Yang-Mills theory.