Color-octet nonrelativistic QCD matrix elements for heavy quarkonium decays in the refined Gribov-Zwanziger theory

Hee Sok Chung

DOI: https://doi.org/10.1103/physrevd.109.054001

IF: 5.407

2024-03-02

Physical Review D

Abstract:We determine color-octet nonrelativistic QCD matrix elements for quarkonium decays from moments of the two-point correlation function of the QCD field-strength tensor computed in the refined Gribov-Zwanziger theory. We find that a tree-level calculation in the refined Gribov-Zwanziger theory can give a suitable description of the QCD field-strength correlation function at both short and long distances, which leads to moments that are infrared finite and can be properly renormalized. By using the color-octet matrix elements we obtain, we quantitatively improve the nonrelativistic effective field theory description of quarkonium decay rates, especially for the χQJ and ηQ states, where Q=c or b . https://doi.org/10.1103/PhysRevD.109.054001 Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI. Funded by SCOAP 3 . Published by the American Physical Society

astronomy & astrophysics,physics, particles & fields

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to determine the octet non - relativistic QCD (NRQCD) matrix elements in heavy quarkonium decays by calculating the two - point correlation function of the quantum chromodynamics (QCD) field - strength tensor and its moments within the framework of the refined Gribov - Zwanziger (RGZ) theory. Specifically, the paper focuses on the following aspects:

1. **Calculation of the field - strength correlation function**: The author calculates the two - point correlation function of the QCD field - strength tensor in the RGZ theory and verifies its consistency with the lattice QCD results. These correlation functions perform well at both short and long distances and can provide infrared - finite and appropriately renormalizable moments.

2. **Determination of the octet matrix elements**: Through the moments of the calculated field - strength correlation function, the author determines the octet NRQCD matrix elements. These matrix elements are particularly important for describing the decay rates of heavy quarkonium, especially in the decay of P - wave heavy quarkonium, where the octet contribution appears at the leading order of the non - relativistic expansion.

3. **Improvement of the theoretical description**: Using these octet matrix elements, the author quantitatively improves the description of the decay rates of heavy quarkonium by the non - relativistic effective field theory, especially for the decay rates of the $\chi_{QJ}$ and $\eta_Q$ states.

### Formula display

1. **Definition of the field - strength correlation function**:
\[
D_{\mu\nu\rho\sigma}(x) = T_F \langle \Omega | g G^a_{\mu\nu}(x) \Phi^{ab}(x,0) g G^b_{\rho\sigma}(0) | \Omega \rangle
\]
where $\langle \Omega \rangle$ is the QCD vacuum, $G^a_{\mu\nu} = \partial_\mu A^a_\nu - \partial_\nu A^a_\mu + g f^{abc} A^b_\mu A^c_\nu$ is the QCD field - strength tensor, $A_\mu$ is the gauge field, $g$ is the strong coupling constant, $T_F=\frac{1}{2}$, and $\Phi^{ab}(x,0)$ is the straight Wilson line.

2. **General form of the field - strength correlation function**:
\[
D_{\mu\nu\rho\sigma}(x) = (\delta_{\mu\rho} \delta_{\nu\sigma} - \delta_{\mu\sigma} \delta_{\nu\rho}) \left[ D(x^2)+D_1(x^2)+(\delta_{\mu\rho} x_\nu x_\sigma - \delta_{\mu\sigma} x_\nu x_\rho - \delta_{\nu\rho} x_\mu x_\sigma + \delta_{\nu\sigma} x_\mu x_\rho) \frac{\partial D_1(x^2)}{\partial x^2} \right]
\]

3. **Invariant functions under dimensional regularization**:
\[
D_E(x^2)=\frac{1}{d - 1}D_{0i0i}(x)=D(x^2)+D_1(x^2)+x^2\frac{\partial D_1(x^2)}{\partial x^2}
\]
\[
D_B(x^2)=\frac{1}{(d - 1)(d - 2)}D_{ijij}(x)=D(x^2)+D_1(x^2)
\]

4. **Tree - level gluon propagator**:
\[
\langle \Omega | A^a_\mu(x) A^b_\nu(y) | \Omega \rangle_{\text{RGZ}}=\delta^{ab}\int\frac{d^d k}{(2\pi)^d}e^{