Abstract

In this paper we study the contribution of the double parton distributions of gluons to charmonium production. Despite being suppressed in the heavy quark mass limit, numerically, this contribution gives a sizeable correction to the leading order k(T) factorization result in LHC kinematics due to the enhancement of gluonic densities in the small Bjorken x(B) limit, and thus presents a complementary mechanism of charmonia production in this kinematics.