Heavy quarkonium decay <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mi>V</mml:mi> <mml:mo stretchy="false">→</mml:mo> <mml:mi>g</mml:mi> <mml:mi>g</mml:mi> <mml:mi>g</mml:mi> </mml:math> with both relativistic and Q

In the heavy quarkonium decay process <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"> <a:mi>V</a:mi> <a:mo stretchy="false">→</a:mo> <a:mi>g</a:mi> <a:mi>g</a:mi> <a:mi>g</a:mi> </a:math> ( <d:math xmlns:d="http://www.w3.org/1998/Math/MathML" display="inline"> <d:mi>V</d:mi> <d:mo>=</d:mo> <d:mi>J</d:mi> <d:mo>/</d:mo> <d:mi>ψ</d:mi> <d:mo>,</d:mo> <d:mi mathvariant="normal">ϒ</d:mi> </d:math> ), making a definite prediction including relativistic corrections has so far remained a significant challenge. In this work, we study this decay process by taking into account the relativistic corrections in the Bethe-Salpeter formalism, where the relativistic bound-state wave function of quarkonium is obtained by solving the Bethe-Salpeter equation under the covariant instantaneous ansatz. Through analytical calculation, we find that some polarized decay widths vanish due to the helicity selection rule, which suppresses the corresponding helicity amplitudes. Owing to helicity flip symmetry and phase space symmetry, the nonvanishing polarized decay widths are not all independent; they are related through a set of symmetry relations. Then we obtain the unpolarized decay width formula <g:math xmlns:g="http://www.w3.org/1998/Math/MathML" display="inline"> <g:mi mathvariant="normal">Γ</g:mi> <g:mo stretchy="false">(</g:mo> <g:mi>V</g:mi> <g:mo stretchy="false">→</g:mo> <g:mi>g</g:mi> <g:mi>g</g:mi> <g:mi>g</g:mi> <g:mo stretchy="false">)</g:mo> <g:mo>=</g:mo> <g:mfrac> <g:mrow> <g:mn>80</g:mn> <g:mo stretchy="false">(</g:mo> <g:msup> <g:mi>π</g:mi> <g:mn>2</g:mn> </g:msup> <g:mo stretchy="false">−</g:mo> <g:mn>9</g:mn> <g:mo stretchy="false">)</g:mo> <g:msubsup> <g:mi>α</g:mi> <g:mi>s</g:mi> <g:mn>3</g:mn> </g:msubsup> <g:msubsup> <g:mi>N</g:mi> <g:mi>V</g:mi> <g:mn>2</g:mn> </g:msubsup> <g:msubsup> <g:mi>β</g:mi> <g:mi>V</g:mi> <g:mn>3</g:mn> </g:msubsup> </g:mrow> <g:mrow> <g:mn>81</g:mn> <g:msup> <g:mi>π</g:mi> <g:mrow> <g:mn>9</g:mn> <g:mo stretchy="false">/</g:mo> <g:mn>2</g:mn> </g:mrow> </g:msup> <g:mi>M</g:mi> </g:mrow> </g:mfrac> <g:mo stretchy="false">(</g:mo> <g:mn>1</g:mn> <g:mo>−</g:mo> <g:mi>κ</g:mi> <g:mfrac> <g:msubsup> <g:mi>β</g:mi> <g:mi>V</g:mi> <g:mn>2</g:mn> </g:msubsup> <g:msup> <g:mi>M</g:mi> <g:mn>2</g:mn> </g:msup> </g:mfrac> <g:mo stretchy="false">)</g:mo> </g:math> , where the factor <s:math xmlns:s="http://www.w3.org/1998/Math/MathML" display="inline"> <s:mi>κ</s:mi> <s:mfrac> <s:msubsup> <s:mi>β</s:mi> <s:mi>V</s:mi> <s:mn>2</s:mn> </s:msubsup> <s:msup> <s:mi>M</s:mi> <s:mn>2</s:mn> </s:msup> </s:mfrac> </s:math> arises from the relativistic corrections with <u:math xmlns:u="http://www.w3.org/1998/Math/MathML" display="inline"> <u:mi>κ</u:mi> <u:mo>≡</u:mo> <u:mfrac> <u:mrow> <u:mn>3</u:mn> <u:mo stretchy="false">(</u:mo> <u:mn>112</u:mn> <u:mo stretchy="false">+</u:mo> <u:mn>25</u:mn> <u:msup> <u:mi>π</u:mi> <u:mn>2</u:mn> </u:msup> <u:mo stretchy="false">)</u:mo> </u:mrow> <u:mrow> <u:mn>16</u:mn> <u:mo stretchy="false">(</u:mo> <u:msup> <u:mi>π</u:mi> <u:mn>2</u:mn> </u:msup> <u:mo stretchy="false">−</u:mo> <u:mn>9</u:mn> <u:mo stretchy="false">)</u:mo> </u:mrow> </u:mfrac> </u:math> . Furthermore, including both relativistic and QCD radiative corrections within the factorization assumption, our predictions of <cb:math xmlns:cb="http://www.w3.org/1998/Math/MathML" display="inline"> <cb:mi mathvariant="script">B</cb:mi> <cb:mo stretchy="false">(</cb:mo> <cb:mi>V</cb:mi> <cb:mo stretchy="false">→</cb:mo> <cb:mi>g</cb:mi> <cb:mi>g</cb:mi> <cb:mi>g</cb:mi> <cb:mo stretchy="false">)</cb:mo> </cb:math> and <ib:math xmlns:ib="http://www.w3.org/1998/Math/MathML" display="inline"> <ib:mi mathvariant="script">B</ib:mi> <ib:mo stretchy="false">(</ib:mo> <ib:mi>V</ib:mi> <ib:mo stretchy="false">→</ib:mo> <ib:msup> <ib:mi>e</ib:mi> <ib:mo>+</ib:mo> </ib:msup> <ib:msup> <ib:mi>e</ib:mi> <ib:mo>−</ib:mo> </ib:msup> <ib:mo stretchy="false">)</ib:mo> </ib:math> agree well with their experimental data. As a crossing check, with the experimental value of the ratio <ob:math xmlns:ob="http://www.w3.org/1998/Math/MathML" display="inline"> <ob:msub> <ob:mi>R</ob:mi> <ob:mi>V</ob:mi> </ob:msub> <ob:mo>=</ob:mo> <ob:mfrac> <ob:mrow> <ob:mi mathvariant="normal">Γ</ob:mi> <ob:mo stretchy="false">(</ob:mo> <ob:mi>V</ob:mi> <ob:mo stretchy="false">→</ob:mo> <ob:mi>g</ob:mi> <ob:mi>g</ob:mi> <ob:mi>g</ob:mi> <ob:mo stretchy="false">)</ob:mo> </ob:mrow> <ob:mrow> <ob:mi mathvariant="normal">Γ</ob:mi> <ob:mo stretchy="false">(</ob:mo> <ob:mi>V</ob:mi> <ob:mo stretchy="false">→</ob:mo> <ob:msup> <ob:mi>e</ob:mi> <ob:mo stretchy="false">+</ob:mo> </ob:msup> <ob:msup> <ob:mi>e</ob:mi> <ob:mo stretchy="false">−</ob:mo> </ob:msup> <ob:mo stretchy="false">)</ob:mo> </ob:mrow> </ob:mfrac> </ob:math> and our result for <ac:math xmlns:ac="http://www.w3.org/1998/Math/MathML" display="inline"> <ac:msub> <ac:mi>R</ac:mi> <ac:mi>V</ac:mi> </ac:msub> </ac:math> , we extract <cc:math xmlns:cc="http://www.w3.org/1998/Math/MathML" display="inline"> <cc:msub> <cc:mi>α</cc:mi> <cc:mi>s</cc:mi>