We calculate the gg → γ γ amplitude by including the tt bound-state effects near their mass threshold. In terms of the non-relativistic expansion of the amplitude, the LO contribution is an energy-independent term in the oneloop amplitude. We include the NLO contribution described by the non-relativistic Green function and part of the NNLO contribution. Despite a missing NLO piece which can be accomplished with the two-loop-level amplitude via massive quarks, the shape of the diphoton mass spectrum is predicted with a good accuracy. Thanks to the simple and clean nature of the observable, its experimental measurement can be a direct method to determine the short-distance mass of the top quark at hadron colliders.At the LHC, a diphoton mass spectrum dσ/dm γ γ has attracted broad attention for observations of the properties of the Higgs boson in the standard model (SM) [1][2][3][4][5] and searches for new phenomena beyond the SM [6][7][8][9][10][11]. At hadron colliders, pairs of high-p T photon are produced by qq annihilation and gluon-fusion mechanisms [12][13][14], and processes which involve fragmentation-photon contributions [15]. The gg → γ γ process, which is the main focus of this letter, is described by loop diagrams with quarks in the SM. The analytic expression of the one-loop amplitude has been known for a long time for both the massless-and the massive-quark loops [16][17][18][19][20][21][22]. The two-loop amplitude has been calculated only for the massless-quark loops [23][24][25].The threshold structure of the massive-quark-loop amplitude deserves particular interests [21,26,27] where the massive quark is regarded as the top quark or a hypothetical particle beyond the SM. Beyond the one-loop level, the amplitude receives large QCD corrections due to the Coulomb-gluon exchanges between the nearly on-shell and low-velocity heavy quarks in s-channel. Thus, the description of the amplitude requires an elaborate treatment based on the non-relativistic QCD formalism. For gg → γ γ process, a e-mail: skawabata@seoultech.ac.kr b e-mail: hyokoya@kias.re.kr such a study cannot be found in the literature. The aim of this letter is to compile present knowledge of the non-relativistic QCD theory for the description of the bound-state effects in the massive-quark-loop amplitude, and to present a dedicated and quantitative study on the diphoton mass spectrum near the tt threshold. Our framework follows the preceding studies on h → γ γ [28,29], and some of our numerical results overlap with that in Ref. [26].We discuss further to utilize the predicted mass spectrum for a precise determination of the top-quark mass, which is one of the fundamental parameters in the SM. Although the top-quark mass has been measured with an error of sub-GeV level [30], its interpretation in terms of well-defined mass parameters is not settled yet in perturbative QCD. It is known that the well-defined mass parameters can be determined by using the threshold scan method at future e + e − colliders [31][32][33]. We show that the ...