In this paper, we attempt to discuss the topics mentioned in the title by scrutinizing and improving the 0 −+ pseudoscalar di-gluonia/glueballs sum rules within the standard SVZ-expansion at N2LO without instantons.First, we reconsider the estimate of the slope of the topological charge χ (0)(Q 2 = 2 GeV 2 ) = 23.9(3.4) MeV from low degree moments which imply 1 0 dx g P 1 (x) = 0.144(5) [data : 0.145(13)] for the first moment of the polarized proton structure function (proton spin) and G (0) A (10 GeV 2 ) ≡ ∆u+∆d+∆s = 0.337(50) [data = 0.330(39)] for the singlet form factor of the axial current. Second, we work with high degree moments and parametrize the spectral function beyond the minimal duality ansatz: "One resonance ⊕ QCD continuum" to get the 0 −+ pseudoscalar di-gluonia spectra. Then, we obtain three groups of gluonia : The familiar light η 1 [singlet gluon component of the η (958)] with [M η1 , f η1 ] = [870, 881(68)] MeV which is important for understanding the U (1) A anomaly; The two new medium gluonia with M P1a = 1328(112) MeV and [M P 1b = 1474(117) MeV or their mean [M P1 , f P1 ] = [1398(81), 576(144)] MeV which support the gluonium nature of the excellent experimental candidate η(1405) and may bring a small gluon piece to the η(1295); their corresponding 1st radial excitations : M P 1a = 1481(216) MeV and M P 1b = 1556(138) MeV with their mean: [M P 1 , f P 1 ] = [1524(116), 230(260)] MeV which may be identified (up to some eventual mixings with qq states) with the observed η(1475, 1700) states; The heavy gluonium with the mean mass: [M P2 , f P2 ] = [2778(133), 500(43)] MeV which is comparable with the lattice results. One can remark the (natural) one to one correspondence between the pseudoscalar gluonia and their chiral scalar analogue from Ref. [1P 2 which is mainly due to the importance of the QCD PT contributions in the sum rule analysis that are almost equal in these two channels. The results for the spectrum are summarized in Table 2.