Within the framework of nonrelativistic QCD factorization formalism, we compute the helicity amplitude as well as the decay width of (, ) electromagnetic decay into two photons up to next-to-next-to-leading order in expansion. For the first time, we verify the validity of nonrelativistic QCD factorization for the -wave quarkonium decay at next-to-next-to-leading order. We find that the and corrections to the helicity amplitude are negative and moderate, nevertheless both corrections combine to suppress the leading-order prediction for the decay width significantly. By approximating the total decay width of as the sum of those for the hadronic decay and the electric transition, we obtain the branching ratios and . To explore the potential measurement on , we further evaluate the production cross section of at LHCb at the lowest order in expansion. With the kinematic constraint on the longitudinal rapidity and transverse momentum for , we find the cross section can reach 2–50 nb for , and 1–22 pb for . Considering the integrated luminosity at , 13 TeV, we estimate that there are several hundreds events of . Since the background is relatively clean, it is promising to reconstruct through its electromagnetic decay. On the other hand, due to the small branching ratio and production cross section, it is quite challenging to detect at LHCb.