Digluon contribution to J/psi production

In this paper we study the contribution of the double parton distributions of gluons to the 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 enhancement of gluonic densities in the small Bjorken $x_{B}$ limit. This contribution is not suppressed at large $J/\psi$ momenta $p_{T}$ and thus presents one of the complementary mechanisms of charmonia production in this kinematics.