Next-to-leading-order QCD corrections to double prompt <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mrow><mml:mrow><mml:mi>J</mml:mi></mml:mrow><mml:mo>/</mml:mo><mml:mi>ψ</mml:mi></mml:mrow></mml:math> hadroproducti

Within the framework of nonrelativistic-QCD factorization, we perform a comprehensive investigation of the color singlet (CS) contribution to prompt <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"><a:mi>J</a:mi><a:mo>/</a:mo><a:mi>ψ</a:mi></a:math> pair production at the CERN Large Hadron Collider at next-to-leading order (NLO) in <c:math xmlns:c="http://www.w3.org/1998/Math/MathML" display="inline"><c:msub><c:mi>α</c:mi><c:mi>s</c:mi></c:msub></c:math>. Specifically, we compare our NLO predictions with measurements from the LHCb, CMS, and ATLAS collaborations. We find that the CS contribution itself can well describe the LHCb data in most of the experimental bins, except in the regions where the fixed-order calculation is spoiled by the emission of soft or hard-collinear gluons. In the CMS and ATLAS cases, however, the CS predictions greatly undershoot both the measured total and differential cross sections, despite sizable <e:math xmlns:e="http://www.w3.org/1998/Math/MathML" display="inline"><e:mi>K</e:mi></e:math> factors, of 2–3, for the total cross sections.