Pseudoscalar quarkonium exclusive decays to vector meson pair

The pseudoscalar quarkonia exclusive decays to light mesons still poses a challenge to the theoretical understanding of quarkonium properties in decay. In this work, we evaluate the processes of pseudoscalar heavy quarkonium decays into vector meson pairs, especially the helicity suppressed processes of ηb→J/ψJ/ψ and ηc→VV. In the framework of NRQCD, the branching fraction of Br[ηb→J/ψJ/ψ] are evaluated at the next-to-leading order of perturbative QCD; and within the light-cone distribution formalism, we calculate also the higher twist effects in these processes. Numerical results show that the higher twist terms contribute more than what from the NLO QCD corrections in the process of ηb→J/ψJ/ψ. It is found that the experimental results on ηc→VV are hard to be understood by merely the quark model and perturbative QCD calculation.