We study the radiative decay widths of vector quarkonia for the process of / ( ) → ( ) and Υ( ) → ( ) (for principal quantum numbers = 1, 2, 3) in the framework of Bethe-Salpeter equation under the covariant instantaneous ansatz using a 4 × 4 form of BSE. The parameters of the framework were determined by a fit to the mass spectrum of ground states of pseudoscalar and vector quarkonia, such as , , / , and Υ. These input parameters so fixed were found to give good agreements with data on mass spectra of ground and excited states of pseudoscalar and vector quarkonia, leptonic decay constants of pseudoscalar and vector quarkonia, two-photon decays, and two-gluon decays of pseudoscalar quarkonia in our recent paper. With these input parameters so fixed, the radiative decay widths of ground (1 ) and excited (2 , 3 ) states of heavy vector quarkonia ( /Ψ and Υ) are calculated and found to be in reasonable agreement with data.