Within the framework of nonrelativistic quark–antiquark Cornell potential model formalism, we study the annihilation of heavy quarkonia. We determine their annihilation widths resulting into digluon, dilepton, 3 γ , 3g and γ g g and compare our findings with the available theoretical results and experimental data. We also provide the charge radii and absolute square of radial Schrödinger wave function at zero quark–antiquark separation for heavy quarkonia and B c mesons. [ABSTRACT FROM AUTHOR]