Although the Large Hadron Collider (LHC) has not observed supersymmetric (SUSY) partners of the Standard Model particles, their existence is not ruled out yet. One recently explored scenario in which there are light SUSY partners that have evaded current bounds from the LHC is that of a light long-lived stop quark. In this paper we consider light stop pair production at the LHC when the stop mass is between 200 and 400 GeV. If the stops are long-lived they can form a bound state, stoponium, which then undergoes two-body decays to Standard Model particles. By considering the near-threshold production of such a pair through the gluon-gluon fusion process and taking into account the strong Coulombic interactions responsible for the formation of this bound state, we obtain factorization theorems for the stop pair inclusive and differential production cross sections. We also perform a resummation of large threshold logarithms up to next-to-next-to-leading logarithmic accuracy using well-established renormalization group equations in an effective field theory methodology. These results are used to calculate the invariant mass distributions of two photons or two Z bosons coming from the decay of the stoponium at the LHC. For our choices of SUSY model parameters, the stoponium is not detectable above Standard Model backgrounds in γγor ZZ at 8 TeV, but will be visible with 400 fb^(-1) of accumulated data if its mass is below 500 GeV when the LHC runs at 14 TeV.

Figures 3, 4, 5 and Tables 1, 2 are revised after fixing a coding error. Conclusions are slightly changed acoordingly. A version to appear in PRD