“…In the traditional non-relativisitic bound state calculation, the two-photon widths for the P -wave quarkonium state depend on the derivative at the origin of the spatial wave function which has to be extracted from potential models [5]. Though the physics of quarkonium decay seems to be better understood within the conventional framework of QCD [6], unlike the two-photon width of S-wave η c and η b quarkonia which can be predicted from the corresponding J/ψ and Υ leptonic widths using HQSS, there is no similar prediction for the Pwave χ c and χ b states and all the existing theoretical values for the decay rates are based on potential model calculations [5,7,8,9,10,11,12,13,14,15,16,17]. To have a prediction for the two-photon width of P -wave quarkonia, one need to express the decay amplitude in terms of the matrix element of a heavy quark field local operator extracted from some known physical processes or computed in an essentially model-independent manner, such as QCD sum rules technique [18,19] or lattice simulations [20].…”