Analysis of the hidden-bottom tetraquark mass spectrum with the QCD sum rules

TL;DR

This paper uses QCD sum rules to systematically analyze the mass spectrum of ground state hidden-bottom tetraquarks, aiming to compare predictions with future experimental data to understand the nature of Z_b states.

Contribution

It extends previous work by providing a systematic QCD sum rule analysis of hidden-bottom tetraquark masses, highlighting effective diquark configurations.

Findings

Scalar, axialvector, and tensor diquark components are effective in forming the lowest tetraquark states.

Predicted tetraquark masses can be tested against future experimental data.

The study offers insights into the internal structure of Z_b states.

Abstract

In this article, we extend our previous work to study the mass spectrum of the ground state hidden-bottom tetraquark states with the QCD sum rules in a systematic way. The predicted hidden-bottom tetraquark masses can be confronted to the experimental data in the future to diagnose the nature of the $Z_b$ states. In calculations, we observe that the scalar diquark states, the axialvector diquark states and the axialvector components of the tensor diquark state are all good diquarks in building the lowest tetraquark states.

Tables

Table 1. Table 1: The quark structures and corresponding current operators for the hidden-bottom tetraquark states.

$Z_b$	$J^{PC}$	Currents
${\left[ub\right]}_S{\left[\overline{db}\right]}_S$	$0^{++}$	$J_{SS}\left(x\right)$
${\left[ub\right]}_A{\left[\overline{db}\right]}_A$	$0^{++}$	$J_{AA}\left(x\right)$
${\left[ub\right]}_{\tilde{A}}{\left[\overline{db}\right]}_{\tilde{A}}$	$0^{++}$	$J_{\tilde{A}\tilde{A}}\left(x\right)$
${\left[ub\right]}_{\tilde{V}}{\left[\overline{db}\right]}_{\tilde{V}}$	$0^{++}$	$J_{\tilde{V}\tilde{V}}\left(x\right)$
${\left[ub\right]}_S{\left[\overline{db}\right]}_A-{\left[ub\right]}_A{\left[\overline{db}\right]}_S$	$1^{+-}$	$J_{-,\mu }^{SA}\left(x\right)$
${\left[ub\right]}_A{\left[\overline{db}\right]}_A$	$1^{+-}$	$J_{-,\mu \nu }^{AA}\left(x\right)$
${\left[ub\right]}_{\tilde{A}}{\left[\overline{db}\right]}_A-{\left[ub\right]}_A{\left[\overline{db}\right]}_{\tilde{A}}$	$1^{+-}$	$J_{-,\mu }^{\tilde{A}A}\left(x\right)$
${\left[ub\right]}_S{\left[\overline{db}\right]}_{\tilde{A}}-{\left[ub\right]}_{\tilde{A}}{\left[\overline{db}\right]}_S$	$1^{+-}$	$J_{-,\mu \nu }^{S\tilde{A}}\left(x\right)$
${\left[ub\right]}_S{\left[\overline{db}\right]}_A+{\left[ub\right]}_A{\left[\overline{db}\right]}_S$	$1^{++}$	$J_{+,\mu }^{SA}\left(x\right)$
${\left[ub\right]}_{\tilde{V}}{\left[\overline{db}\right]}_V-{\left[ub\right]}_V{\left[\overline{db}\right]}_{\tilde{V}}$	$1^{++}$	$J_{+,\mu }^{\tilde{V}V}\left(x\right)$
${\left[ub\right]}_{\tilde{A}}{\left[\overline{db}\right]}_A+{\left[ub\right]}_A{\left[\overline{db}\right]}_{\tilde{A}}$	$1^{++}$	$J_{+,\mu }^{\tilde{A}A}\left(x\right)$
${\left[ub\right]}_A{\left[\overline{db}\right]}_A$	$2^{++}$	$J_{+,\mu \nu }^{AA}\left(x\right)$

Table 2. Table 2: The Borel parameters, continuum threshold parameters, energy scales of the QCD spectral densities, pole contributions, and the contributions of the vacuum condensates of dimension 10 10 10 for the ground state hidden-bottom tetraquark states.

$Z_b$	$J^{PC}$	$T^2\left({\mathrm{GeV}}^2\right)$	$\sqrt{s_0}\left(\mathrm{GeV}\right)$	$\mu \left(\mathrm{GeV}\right)$	pole	$\left|D\left(10\right)\right|$
${\left[ub\right]}_S{\left[\overline{db}\right]}_S$	$0^{++}$	$7.0-8.0$	$11.16\pm 0.10$	$2.40$	$\left(44-66\right)\%$	$\le 3\%$
${\left[ub\right]}_A{\left[\overline{db}\right]}_A$	$0^{++}$	$6.4-7.4$	$11.14\pm 0.10$	$2.30$	$\left(44-68\right)\%$	$\le 11\%$
${\left[ub\right]}_{\tilde{A}}{\left[\overline{db}\right]}_{\tilde{A}}$	$0^{++}$	$7.2-8.2$	$11.17\pm 0.10$	$2.40$	$\left(45-66\right)\%$	$\le 4\%$
${\left[ub\right]}_{\tilde{V}}{\left[\overline{db}\right]}_{\tilde{V}}$	$0^{++}$	$11.4-12.8$	$12.22\pm 0.10$	$5.40$	$\left(44-61\right)\%$	$\ll 1\%$
${\left[ub\right]}_S{\left[\overline{db}\right]}_A-{\left[ub\right]}_A{\left[\overline{db}\right]}_S$	$1^{+-}$	$7.0-8.0$	$11.16\pm 0.10$	$2.40$	$\left(44-66\right)\%$
${\left[ub\right]}_A{\left[\overline{db}\right]}_A$	$1^{+-}$	$7.1-8.1$	$11.17\pm 0.10$	$2.40$	$\left(44-65\right)\%$	$\le 4\%$
${\left[ub\right]}_{\tilde{A}}{\left[\overline{db}\right]}_A-{\left[ub\right]}_A{\left[\overline{db}\right]}_{\tilde{A}}$	$1^{+-}$	$6.9-7.9$	$11.17\pm 0.10$	$2.40$	$\left(44-66\right)\%$	$\le 7\%$
${\left[ub\right]}_S{\left[\overline{db}\right]}_{\tilde{A}}-{\left[ub\right]}_{\tilde{A}}{\left[\overline{db}\right]}_S$	$1^{+-}$	$7.1-8.1$	$11.17\pm 0.10$	$2.40$	$\left(44-66\right)\%$	$\le 4\%$
${\left[ub\right]}_S{\left[\overline{db}\right]}_A+{\left[ub\right]}_A{\left[\overline{db}\right]}_S$	$1^{++}$	$7.1-8.1$	$11.18\pm 0.10$	$2.45$	$\left(44-65\right)\%$	$\le 3\%$
${\left[ub\right]}_{\tilde{V}}{\left[\overline{db}\right]}_V-{\left[ub\right]}_V{\left[\overline{db}\right]}_{\tilde{V}}$	$1^{++}$	$6.8-7.8$	$11.19\pm 0.10$	$2.50$	$\left(44-66\right)\%$	$\le 4\%$
${\left[ub\right]}_{\tilde{A}}{\left[\overline{db}\right]}_A+{\left[ub\right]}_A{\left[\overline{db}\right]}_{\tilde{A}}$	$1^{++}$	$9.7-11.1$	$11.99\pm 0.10$	$4.90$	$\left(44-63\right)\%$	$\ll 1\%$
${\left[ub\right]}_A{\left[\overline{db}\right]}_A$	$2^{++}$	$7.2-8.2$	$11.19\pm 0.10$	$2.50$	$\left(44-65\right)\%$

Table 3. Table 3: The masses and pole residues of the ground state hidden-bottom tetraquark states.

$Z_b$	$J^{PC}$	$M_Z\left(\mathrm{GeV}\right)$	${\lambda }_Z\left({\mathrm{GeV}}^5\right)$
${\left[ub\right]}_S{\left[\overline{db}\right]}_S$	$0^{++}$	$10.61\pm 0.09$	$\left(1.10\pm 0.17\right)\times {10}^{-1}$
${\left[ub\right]}_A{\left[\overline{db}\right]}_A$	$0^{++}$	$10.60\pm 0.09$	$\left(1.61\pm 0.25\right)\times {10}^{-1}$
${\left[ub\right]}_{\tilde{A}}{\left[\overline{db}\right]}_{\tilde{A}}$	$0^{++}$	$10.61\pm 0.09$	$\left(1.81\pm 0.27\right)\times {10}^{-1}$
${\left[ub\right]}_{\tilde{V}}{\left[\overline{db}\right]}_{\tilde{V}}$	$0^{++}$	$11.66\pm 0.12$	$3.03\pm 0.31$
${\left[ub\right]}_S{\left[\overline{db}\right]}_A-{\left[ub\right]}_A{\left[\overline{db}\right]}_S$	$1^{+-}$	$10.61\pm 0.09$	$\left(1.08\pm 0.16\right)\times {10}^{-1}$
${\left[ub\right]}_A{\left[\overline{db}\right]}_A$	$1^{+-}$	$10.62\pm 0.09$	$\left(1.07\pm 0.16\right)\times {10}^{-1}$
${\left[ub\right]}_{\tilde{A}}{\left[\overline{db}\right]}_A-{\left[ub\right]}_A{\left[\overline{db}\right]}_{\tilde{A}}$	$1^{+-}$	$10.62\pm 0.09$	$\left(2.12\pm 0.31\right)\times {10}^{-1}$
${\left[ub\right]}_S{\left[\overline{db}\right]}_{\tilde{A}}-{\left[ub\right]}_{\tilde{A}}{\left[\overline{db}\right]}_S$	$1^{+-}$	$10.62\pm 0.09$	$\left(1.08\pm 0.16\right)\times {10}^{-1}$
${\left[ub\right]}_S{\left[\overline{db}\right]}_A+{\left[ub\right]}_A{\left[\overline{db}\right]}_S$	$1^{++}$	$10.63\pm 0.09$	$\left(1.17\pm 0.17\right)\times {10}^{-1}$
${\left[ub\right]}_{\tilde{V}}{\left[\overline{db}\right]}_V-{\left[ub\right]}_V{\left[\overline{db}\right]}_{\tilde{V}}$	$1^{++}$	$10.63\pm 0.09$	$\left(1.22\pm 0.20\right)\times {10}^{-1}$
${\left[ub\right]}_{\tilde{A}}{\left[\overline{db}\right]}_A+{\left[ub\right]}_A{\left[\overline{db}\right]}_{\tilde{A}}$	$1^{++}$	$11.45\pm 0.14$	$\left(8.52\pm 1.02\right)\times {10}^{-1}$
${\left[ub\right]}_A{\left[\overline{db}\right]}_A$	$2^{++}$	$10.65\pm 0.09$	$\left(1.72\pm 0.24\right)\times {10}^{-1}$

Equations

$$\mu$$

$$\hat{P}{\epsilon }^{abc}{q}^{Tb}(x)C{\sigma }_{\mu \nu }{\gamma }_5{Q}^c(x){\hat{P}}^{-1}$$

$$\hat{P}{\epsilon }^{abc}{q}^{Tb}(x)C{\sigma }_{\mu \nu }{Q}^c(x){\hat{P}}^{-1}$$

$$\hat{P}{\epsilon }^{abc}{q}^{Tb}(x)C{\sigma }_{\mu \nu }^t{\gamma }_5{Q}^c(x){\hat{P}}^{-1}$$

$$\hat{P}{\epsilon }^{abc}{q}^{Tb}(x)C{\sigma }_{\mu \nu }^v{Q}^c(x){\hat{P}}^{-1}$$

$$\hat{P}{\epsilon }^{abc}{q}^{Tb}(x)C{\sigma }_{\mu \nu }^t{Q}^c(x){\hat{P}}^{-1}$$

$$\hat{P}{\epsilon }^{abc}{q}^{Tb}(x)C{\sigma }_{\mu \nu }^v{\gamma }_5{Q}^c(x){\hat{P}}^{-1}$$

$$\Pi (p)$$

$${\Pi }_{\mu \nu }(p)$$

$${\Pi }_{\mu \nu \alpha \beta }(p)$$

$$J_{SS}(x)$$

$$J_{AA}(x)$$

$$J_{\tilde{A}\tilde{A}}(x)$$

$$J_{\tilde{V}\tilde{V}}(x)$$

$$J_{-,\mu }^{SA}(x)$$

$$J_{-,\mu \nu }^{AA}(x)$$

$$J_{-,\mu }^{\tilde{A}A}(x)$$

$$J_{-,\mu \nu }^{S\tilde{A}}(x)$$

$$J_{+,\mu }^{SA}(x)$$

$$J_{+,\mu }^{\tilde{V}V}(x)$$

$$J_{+,\mu }^{\tilde{A}A}(x)$$

$$J_{+,\mu \nu }^{AA}(x)$$

$$\Pi (p)$$

$${\Pi }_{\mu \nu }(p)$$

$${\Pi }_{\mu \nu \alpha \beta }^{AA,-}(p)$$

$${\Pi }_{\mu \nu \alpha \beta }^{S\tilde{A},-}(p)$$

$${\Pi }_{\mu \nu \alpha \beta }^{AA,+}(p)$$

$$⟨0|J(0)|Z_b^{+}(p)⟩$$

$$⟨0|J_{\mu }(0)|Z_b^{+}(p)⟩$$

$$⟨0|J_{-,\mu \nu }^{S\tilde{A}}(0)|Z_b^{-}(p)⟩$$

$$⟨0|J_{-,\mu \nu }^{S\tilde{A}}(0)|Z_b^{+}(p)⟩$$

$$⟨0|J_{-,\mu \nu }^{AA}(0)|Z_b^{+}(p)⟩$$

$$⟨0|J_{-,\mu \nu }^{AA}(0)|Z_b^{-}(p)⟩$$

$$⟨0|J_{+,\mu \nu }^{AA}(0)|Z_b^{+}(p)⟩$$

$${\lambda }_{Z^{+}}^2\exp \left(-\frac{M_{Z^{+}}^2}{T^2}\right)={\int }_{4m_b^2}^{s_0}ds\rho (s)\exp \left(-\frac{s}{T^2}\right).$$

$$M_{Z^{+}}^2$$

$$⟨\bar{q}q⟩(\mu )$$

$$⟨\bar{q}{g}_s\sigma Gq⟩(\mu )$$

$$m_b(\mu )$$

$${\alpha }_s(\mu )$$

$$\mathrm{PC}$$

$$D(n)$$

$$Z_b^{+}(1^{+-})$$

$$Z_b^{+}(0^{++})$$

$$Z_b^{+}(1^{++})$$

Full text

PACS number: 12.39.Mk, 12.38.Lg

Key words: Tetraquark state, QCD sum rules

1 Introduction

In 2011, the Belle collaboration observed the $Z^{\pm}_{b}(10610)$ and $Z_{b}^{\pm}(10650)$ in the $\pi^{\pm}\Upsilon({\rm 1,2,3S})$ and $\pi^{\pm}h_{b}({\rm 1,2P})$ mass spectrum for the first time, the favored quantum numbers (Isospin, $G$-parity, Spin, Parity) are $I^{G}(J^{P})=1^{+}(1^{+})$ [1]. Later, the Belle collaboration updated the values of the masses and widths $M_{Z_{b}(10610)}=(10607.2\pm 2.0)\,\rm{MeV}$, $M_{Z_{b}(10650)}=(10652.2\pm 1.5)\,\rm{MeV}$, $\Gamma_{Z_{b}(10610)}=(18.4\pm 2.4)\,\rm{MeV}$ and $\Gamma_{Z_{b}(10650)}=(11.5\pm 2.2)\,\rm{MeV}$ [2], which are adopted in the Review of Particle Physics by the Particle Data Group now [3]. The possible assignments of the $Z^{\pm}_{b}(10610)$ and $Z_{b}^{\pm}(10650)$ are the tetraquark states [4, 5, 6], molecular states [7, 8, 9, 10, 11], threshold cusps [12], re-scattering effects [13], etc. For more literatures on the $Z_{b}$ states, one can consult the old review [14] and the recent review [15].

In 2013, the BESIII collaboration (also the Belle collaboration) observed the $Z_{c}^{\pm}(3900)$ in the $\pi^{\pm}J/\psi$ mass spectrum [16, 17]. Later, the BESIII collaboration observed the $Z^{\pm}_{c}(4025)$ near the $(D^{*}\bar{D}^{*})^{\pm}$ threshold [18], and the $Z_{c}^{\pm}(4020)$ in the $\pi^{\pm}h_{c}$ mass spectrum [19]. Now the $Z_{c}^{\pm}(4025)$ and $Z_{c}^{\pm}(4020)$ are taken to be the same particle, and are denoted as the $X(4020)$ in the Review of Particle Physics [3].

The $Z_{c}^{\pm}(3900)$ and $Z_{c}^{\pm}(4020)$ ($Z^{\pm}_{b}(10610)$ and $Z_{b}^{\pm}(10650)$) are charged charmonium-like states (bottomonium-like states), their quark constituents must be $c\bar{c}u\bar{d}$ or $c\bar{c}d\bar{u}$ ($b\bar{b}u\bar{d}$ or $b\bar{b}d\bar{u}$), irrespective of the diquark-antidiquark type or meson-meson type substructures. The $Z_{b}(10610)$, $Z_{b}(10650)$, $Z_{c}(3900)$ and $Z_{c}(4020)$ are observed in the analogous decays to the final states $\pi^{\pm}h_{b}({\rm 1,2P})$, $\pi^{\pm}\Upsilon({\rm 1,2,3S})$, $\pi^{\pm}J/\psi$, $\pi^{\pm}h_{c}$, and should have analogous structures. In Refs.[6, 20, 21], we assign the $Z_{b}(10610)$, $Z_{b}(10650)$, $Z_{c}(3900)$ and $Z_{c}(4020)$ to be the diquark-antidiquark type axialvector tetraquark states, and study their masses with the QCD sum rules in details. Furthermore, we explore the energy scale dependence of the hidden-charm and hidden-bottom tetraquark states for the first time [20], and suggest a formula,

[TABLE]

with the effective heavy mass ${\mathbb{M}}_{Q}$ to determine the optimal energy scales $\mu$ of the QCD spectral densities [6, 21]. The experimental values of the masses can be well reproduced. In Ref.[22], we study the two-body strong decays $Z_{c}^{\pm}(3900)\to J/\psi\pi^{\pm}$, $\eta_{c}\rho^{\pm}$, $(D\bar{D}^{*})^{\pm}$ with the QCD sum rules in details. We take into account both the connected and disconnected Feynman diagrams, and pay special attentions to matching the hadron side with the QCD side of the correlation functions to obtain solid duality, the predicted width $\Gamma_{Z_{c}(3900)}=54.2\pm 29.8\,\rm{MeV}$ is consistent with the experimental data [16, 17], and supports assigning the $Z_{c}^{\pm}(3900)$ to be the diquark-antidiquark type axialvector tetraquark state [20]. In Ref.[6], we use the method proposed in Ref.[23] to study the two-body strong decays $Z_{b}^{\pm}(10610)\to\Upsilon\pi^{\pm}$, $\Upsilon\rho^{\pm}$ with the QCD sum rules by taking into account only the connected Feynman diagrams. Although the predictions are good, the subtractions of the higher resonances and continuum states are introduced by hand, the contaminations cannot be subtracted completely. The widths of the $Z_{b}^{\pm}(10610)$ and $Z_{b}^{\pm}(10650)$ should be studied in a consistent way to make the assignments more robust. A updated analysis of the masses and widths with the QCD sum rules is needed.

If the $Z_{b}^{\pm}(10610)$ and $Z_{b}^{\pm}(10650)$ are the diquark-antidiquark type axialvector hidden-bottom tetraquark states, there should exist a holonomic spectrum for the scalar, axialvector and tensor hidden-bottom tetraquark states without introducing an additional P-wave. Now we extend our previous work to study the mass spectrum of the hidden-bottom tetraquark states in a systematic way. Those hidden-bottom tetraquark states may be observed at the LHCb, Belle II, CEPC (Circular Electron Positron Collider), FCC (Future Circular Collider), ILC (International Linear Collider) in the future, and shed light on the nature of the exotic $X$, $Y$, $Z$ particles.

We usually take the diquarks (in color antitriplet) and antidiquarks (in color triplet) as the basic building blocks to construct the tetraquark states. The diquarks (or diquark operators) $\varepsilon^{abc}q^{T}_{b}C\Gamma q^{\prime}_{c}$ have five structures in Dirac spinor space, where $C\Gamma=C\gamma_{5}$, $C$, $C\gamma_{\mu}\gamma_{5}$, $C\gamma_{\mu}$ and $C\sigma_{\mu\nu}$ (or $C\sigma_{\mu\nu}\gamma_{5}$) for the scalar ($S$), pseudoscalar ($P$), vector ($V$), axialvector ($A$) and tensor ($T$) diquarks, respectively, the $a$, $b$, $c$ are color indexes. The scalar, pseudoscalar, vector and axialvector diquark states have been studied with the QCD sum rules in details, the good diquark correlations in building the lowest tetraquark states are the scalar and axialvector diquark states [24, 25], the axialvector diquark states are not bad diquark states.

Under parity transform $\widehat{P}$, the tensor diquark operators have the properties,

[TABLE]

where $x^{\mu}=(t,\vec{x})$ and $\tilde{x}^{\mu}=(t,-\vec{x})$. The tensor diquark states have both $J^{P}=1^{+}$ and $1^{-}$ components, we introduce the four vector $t^{\mu}=(1,\vec{0})$ and project out the $1^{+}$ and $1^{-}$ components explicitly,

[TABLE]

where \sigma^{t}_{\mu\nu}=\frac{i}{2}\Big{[}\gamma^{t}_{\mu},\gamma^{t}_{\nu}\Big{]}, \sigma^{v}_{\mu\nu}=\frac{i}{2}\Big{[}\gamma^{v}_{\mu},\gamma^{t}_{\nu}\Big{]}, $\gamma^{v}_{\mu}=\gamma\cdot tt_{\mu}$, $\gamma^{t}_{\mu}=\gamma_{\mu}-\gamma\cdot tt_{\mu}$ [26]. Thereafter, we will denote the axialvector diquark operators $\varepsilon^{abc}q^{Tb}(x)C\sigma^{v}_{\mu\nu}Q^{c}(x)$, $\varepsilon^{abc}q^{Tb}(x)C\sigma^{t}_{\mu\nu}\gamma_{5}Q^{c}(x)$ as $\widetilde{A}$, and the vector diquark operators $\varepsilon^{abc}q^{Tb}(x)C\sigma^{t}_{\mu\nu}Q^{c}(x)$, $\varepsilon^{abc}q^{Tb}(x)C\sigma^{v}_{\mu\nu}\gamma_{5}Q^{c}(x)$ as $\widetilde{V}$.

In this article, we take the scalar ($S$), axialvector ($A$, $\widetilde{A}$), vector ($V$, $\widetilde{V}$) diquark operators and antidiquark operators as the basic building blocks to construct the tetraquark operators to study the mass spectrum of the hidden-bottom tetraquark states with the QCD sum rules in a systematic way.

The article is arranged as follows: we derive the QCD sum rules for the masses and pole residues of the hidden-bottom tetraquark states in section 2; in section 3, we present the numerical results and discussions; section 4 is reserved for our conclusion.

2 QCD sum rules for the hidden-bottom tetraquark states

We write down the two-point correlation functions $\Pi(p)$, $\Pi_{\mu\nu}(p)$ and $\Pi_{\mu\nu\alpha\beta}(p)$ in the QCD sum rules,

[TABLE]

where the currents $J(x)=J_{SS}(x)$, $J_{AA}(x)$, $J_{\widetilde{A}\widetilde{A}}(x)$, $J_{\widetilde{V}\widetilde{V}}(x)$, $J_{\mu}(x)=J^{SA}_{-,\mu}(x)$, $J_{-,\mu}^{\widetilde{A}A}(x)$, $J^{SA}_{+,\mu}(x)$, $J_{+,\mu}^{\widetilde{V}V}(x)$, $J_{+,\mu}^{\widetilde{A}A}(x)$, $J_{\mu\nu}(x)=J^{AA}_{-,\mu\nu}(x)$, $J^{S\widetilde{A}}_{-,\mu\nu}(x)$, $J^{AA}_{+,\mu\nu}(x)$,

[TABLE]

the $i$, $j$, $k$, $m$, $n$ are color indexes, the $C$ is the charge conjugation matrix, the subscripts $\pm$ denote the positive charge conjugation and negative charge conjugation, respectively. The diquark operators $\varepsilon^{ijk}q^{Tj}(x)C\sigma_{\mu\nu}\gamma_{5}b^{j}(x)$ and $\varepsilon^{ijk}q^{Tj}(x)C\sigma_{\mu\nu}b^{k}(x)$ have both positive parity and negative parity components, we project out the $J^{P}=1^{+}$ and $1^{-}$ components unambiguously with suitable diquark operators to obtain the current operators $J_{-,\mu}^{\widetilde{A}A}(x)$, $J_{+,\mu}^{\widetilde{V}V}(x)$ and $J_{+,\mu}^{\widetilde{A}A}(x)$ with the $J^{P}=1^{+}$. The current operators $J^{AA}_{-,\mu\nu}(x)$ and $J^{S\widetilde{A}}_{-,\mu\nu}(x)$ couple potentially to both the positive parity and negative parity tetraquark states, we separate those contributions explicitly to obtain reliable QCD sum rules. In Table 1, we present the quark structures and corresponding interpolating currents for the hidden-bottom tetraquark states. The four vector $t^{\mu}=(1,\vec{0})$ breaks down Lorentz covariance, the currents $J_{\widetilde{A}\widetilde{A}}(x)$, $J_{\widetilde{V}\widetilde{V}}(x)$, $J_{-,\mu}^{\widetilde{A}A}(x)$, $J^{S\widetilde{A}}_{-,\mu\nu}(x)$, $J_{+,\mu}^{\widetilde{V}V}(x)$ and $J_{+,\mu}^{\widetilde{A}A}(x)$ are not Lorentz covariant, it is the shortcoming of the present method, the calculations can be understood as carried out at a particular (or given) coordinate system, which cannot impair the predictive ability.

At the hadron side, we insert a complete set of intermediate hadronic states with the same quantum numbers as the current operators $J(x)$, $J_{\mu}(x)$ and $J_{\mu\nu}(x)$ into the correlation functions $\Pi(p)$, $\Pi_{\mu\nu}(p)$ and $\Pi_{\mu\nu\alpha\beta}(p)$ to obtain the hadronic representation [27, 28], and isolate the ground state hidden-bottom tetraquark contributions,

[TABLE]

where $\widetilde{g}_{\mu\nu}=g_{\mu\nu}-\frac{p_{\mu}p_{\nu}}{p^{2}}$, the superscripts (subscripts) $\pm$ in the hidden-bottom tetraquark states $Z_{b}^{\pm}$ (components $\Pi_{\pm}(p^{2})$, $\widetilde{\Pi}_{\pm}(p^{2})$) denote the positive-parity and negative parity, respectively. The pole residues $\lambda_{Z^{\pm}}$ are defined by

[TABLE]

the $\varepsilon_{\mu/\alpha}$ and $\varepsilon_{\mu\nu}$ are the polarization vectors of the hidden-bottom tetraquark states. In this article, we choose the components $\Pi_{+}(p^{2})$ and $p^{2}\widetilde{\Pi}_{+}(p^{2})$ to study the scalar, axialvector and tensor hidden-bottom tetraquark states with the QCD sum rules.

At the QCD side, we carry out the operator product expansion for the correlation functions $\Pi(p)$, $\Pi_{\mu\nu}(p)$ and $\Pi_{\mu\nu\alpha\beta}(p)$ up to the vacuum condensates of dimension $10$ in a consistent way, and obtain the QCD spectral densities $\rho(s)$ through dispersion relation. We match the hadron side with the QCD side of the correlation functions below the continuum threshold $s_{0}$ and perform Borel transform with respect to $P^{2}=-p^{2}$ to obtain the QCD sum rules:

[TABLE]

The explicit expressions of the QCD spectral densities $\rho(s)$ are available upon request by contacting me via E-mail. For the technical details, one can consult Refs.[20, 21].

We derive Eq.(8) with respect to $\tau=\frac{1}{T^{2}}$, and obtain the QCD sum rules for the masses of the scalar, axialvector and tensor hidden-bottom tetraquark states $Z_{b}$ through a ratio,

[TABLE]

3 Numerical results and discussions

We take the standard values of the vacuum condensates $\langle\bar{q}q\rangle=-(0.24\pm 0.01\,\rm{GeV})^{3}$, $\langle\bar{q}g_{s}\sigma Gq\rangle=m_{0}^{2}\langle\bar{q}q\rangle$, $m_{0}^{2}=(0.8\pm 0.1)\,\rm{GeV}^{2}$, $\langle\frac{\alpha_{s}GG}{\pi}\rangle=(0.33\,\rm{GeV})^{4}$ at the energy scale $\mu=1\,\rm{GeV}$ [27, 28, 29], and take the $\overline{MS}$ mass $m_{b}(m_{b})=(4.18\pm 0.03)\,\rm{GeV}$ from the Particle Data Group [3], and set $m_{u}=m_{d}=0$. Furthermore, we take into account the energy-scale dependence of the input parameters at the QCD side,

[TABLE]

where $t=\log\frac{\mu^{2}}{\Lambda^{2}}$, $b_{0}=\frac{33-2n_{f}}{12\pi}$, $b_{1}=\frac{153-19n_{f}}{24\pi^{2}}$, $b_{2}=\frac{2857-\frac{5033}{9}n_{f}+\frac{325}{27}n_{f}^{2}}{128\pi^{3}}$, $\Lambda=210\,\rm{MeV}$, $292\,\rm{MeV}$ and $332\,\rm{MeV}$ for the flavors $n_{f}=5$, $4$ and $3$, respectively [3, 30], and evolve all the input parameters to the optimal energy scales $\mu$ with the flavor $n_{f}=5$ to extract the tetraquark masses.

In all the QCD sum rules for the hidden-charm tetraquark states [21, 31], hidden-charm pentaquark states [32] and hidden-bottom tetraquark states [6], we search for the optimal Borel parameters $T^{2}$ and continuum threshold parameters $s_{0}$ to satisfy the four criteria:

$\bf 1.$ Pole dominance at the hadron side;

$\bf 2.$ Convergence of the operator product expansion at the QCD side;

$\bf 3.$ Appearance of the Borel platforms;

$\bf 4.$ Satisfying the energy scale formula,

via try and error.

The pole contributions (PC) or ground state tetraquark contributions are defined by

[TABLE]

this definition is adopted in all the QCD sum rules. The $Z_{b}(10610)$, $Z_{b}(10650)$, $Z_{c}(3900)$ and $Z_{c}(4020)$ have analogous properties, the $Z_{c}(4430)$ can be tentatively assigned to be the first radial excited state of the $Z_{c}(3900)$ [31, 33, 34]. The relevant mass gaps $M_{Z_{c}(4430)}-M_{Z_{c}(3900)}=591\,\rm{MeV}$, $M_{\psi^{\prime}}-M_{J/\psi}=589\,\rm{GeV}$, $M_{\Upsilon^{\prime}}-M_{\Upsilon}=563\,\rm{MeV}$ from the Particle Data Group [3]. The energy gaps at the charmonium section have the relation $M_{Z^{\prime}_{c}}-M_{Z_{c}}=M_{\psi^{\prime}}-M_{J/\psi}$, we expect such a relation survives in the bottom section $M_{Z^{\prime}_{b}}-M_{Z_{b}}=M_{\Upsilon^{\prime}}-M_{\Upsilon}=0.55\,\rm{GeV}$. In this article, we choose the continuum threshold parameters $\sqrt{s_{0}}=Z_{b}+0.55\pm 0.10\,\rm{GeV}$ as a constraint.

To judge the convergence of the operator product expansion, we calculate the contributions of the vacuum condensates $D(n)$ with the formula,

[TABLE]

rather than with the formula,

[TABLE]

The definition in Eq.(13) works only when all the contributions at the hadron side are included, such as the ground state, first radial excited state, second radial excited state, $\cdots$, continuum states, where the index $n$ denotes the dimension of the vacuum condensates.

The energy scale formula $\mu=\sqrt{M^{2}_{X/Y/Z}-(2{\mathbb{M}}_{Q})^{2}}$ for the QCD spectral densities can enhance the pole contributions remarkably and improve the convergence of the operator product expansion considerably. In Ref.[6], we take the energy scale formula $\mu=\sqrt{M^{2}_{X/Y/Z}-(2{\mathbb{M}}_{b})^{2}}$ with the effective $b$-quark mass ${\mathbb{M}}_{b}=5.13\,\rm{GeV}$ to determine the ideal energy scales of the QCD spectral densities. After its publication, we re-checked the numerical calculations and found that there existed a small error involving the mixed condensates. We corrected the small error and observed that the Borel windows were modified slightly and the predicted masses were improved slightly, but the conclusions survived, the updated effective $b$-quark mass was ${\mathbb{M}}_{b}=5.17\,\rm{GeV}$ [35]. In this article, we take the updated value ${\mathbb{M}}_{b}=5.17\,\rm{GeV}$. Moreover, we recalculate the high dimensional vacuum condensates using the formula $t^{a}_{ij}t^{a}_{mn}=-\frac{1}{6}\delta_{ij}\delta_{mn}+\frac{1}{2}\delta_{jm}\delta_{in}$, and obtain slightly different expressions compared to the old calculations, where $t^{a}=\frac{\lambda^{a}}{2}$, the $\lambda^{a}$ is the Gell-Mann matrix.

We obtain the Borel parameters, continuum threshold parameters, energy scales of the QCD spectral densities, pole contributions, and the contributions of the vacuum condensates of dimension $10$, which are shown explicitly in Table 2. From the Table, we can see that the pole contributions are about $(44-66)\%$, the pole dominance condition is well satisfied.

In calculations, we observe that the contributions of the vacuum condensates of dimension $10$ $|D(10)|\leq 4\%$ for the most QCD sum rules, the operator product expansion is well convergent. In the QCD sum rules for the tetraquark state $[ub]_{A}[\overline{db}]_{A}$ with $J^{PC}=0^{++}$, the $|D(10)|=(4\sim 11)\%$, which is somewhat large. If we choose Borel window $T^{2}=(6.7-7.7)\,\rm{GeV}^{2}$, we can obtain a lightly smaller pole contribution, ${\rm PC}=(40-63)\%$, and slightly smaller $D(10)$, $|D(10)|=(3\sim 8)\%$, which is acceptable. As the predicted mass changes slowly with variation of the Borel parameter, the Borel windows $T^{2}=(6.7-7.7)\,\rm{GeV}^{2}$ and $(6.4-7.4)\,\rm{GeV}^{2}$ lead to the same tetraquark mass. The two basic criteria of the QCD sum rules are satisfied.

We take into account all the uncertainties of the input parameters and obtain the masses and pole residues of the scalar, axialvector, tensor hidden-bottom tetraquark states, which are shown explicitly in Table 3 and in Figs.1–4. From Tables 2–3, we can see that the energy scale formula $\mu=\sqrt{M^{2}_{X/Y/Z}-(2{\mathbb{M}}_{b})^{2}}$ is well satisfied. In Figs.1–4, we plot the masses and pole residues of the scalar, axialvector, tensor hidden-bottom tetraquark states with variations of the Borel parameters at much larger ranges than the Borel widows, the regions between the two perpendicular lines are the Borel windows, where the Borel platforms appear. Now the four criteria of the QCD sum rules are all satisfied.

In Figs.1–2, we also present the experimental values of the masses of the $Z_{b}(10610)$ and $Z_{b}(10650)$ [2, 3]. From the figures, we can see that the masses of all the $[ub]_{S}[\overline{db}]_{A}-[ub]_{A}[\overline{db}]_{S}$, $[ub]_{A}[\overline{db}]_{A}$, $[ub]_{\widetilde{A}}[\overline{db}]_{A}-[ub]_{A}[\overline{db}]_{\widetilde{A}}$, $[ub]_{S}[\overline{db}]_{\widetilde{A}}-[ub]_{\widetilde{A}}[\overline{db}]_{S}$ hidden-bottom tetraquark states with the $J^{PC}=1^{+-}$ are in excellent agreements with the experimental values $M_{Z_{b}(10610)}=(10607.2\pm 2.0)\,\rm{MeV}$ and $M_{Z_{b}(10650)}=(10652.2\pm 1.5)\,\rm{MeV}$ within uncertainties [2, 3]. We cannot assign the $Z_{b}(10610)$ and $Z_{b}(10650)$ unambiguously with the mass alone, we should study the partial decay widths exclusively to obtain a more robust assignment.

From Table 3, we can see that the scalar tetraquark states $[ub]_{S}[\overline{db}]_{S}$, $[ub]_{A}[\overline{db}]_{A}$, $[ub]_{\tilde{A}}[\overline{db}]_{\tilde{A}}$, axialvector tetraquark states $[ub]_{S}[\overline{db}]_{A}-[ub]_{A}[\overline{db}]_{S}$, $[ub]_{A}[\overline{db}]_{A}$, $[ub]_{\widetilde{A}}[\overline{db}]_{A}-[ub]_{A}[\overline{db}]_{\widetilde{A}}$, $[ub]_{S}[\overline{db}]_{\widetilde{A}}-[ub]_{\widetilde{A}}[\overline{db}]_{S}$, $[ub]_{S}[\overline{db}]_{A}+[ub]_{A}[\overline{db}]_{S}$, $[ub]_{\widetilde{V}}[\overline{db}]_{V}-[ub]_{V}[\overline{db}]_{\widetilde{V}}$, and tensor tetraquark state $[ub]_{A}[\overline{db}]_{A}$ have almost degenerated masses, i.e. about $10.6\,\rm{GeV}$. The calculations based on the QCD sum rules indicate that the scalar ($S$) and axialvector ($A$) bottom diquark states have degenerated masses [25], furthermore, the spin-spin interaction is proportional to $\frac{\vec{S}_{i}\cdot\vec{S}_{j}}{m_{i}m_{j}}$, the mass of the $b$-quark is large, so it is reasonable that the lowest scalar, axialvector, tensor hidden-bottom tetraquark have almost degenerated masses. The scalar diquark states $S$ and axialvector diquark states $A$, $\widetilde{A}$ are all good diquark states in building the lowest tetraquark states.

The hidden-bottom tetraquark masses obtained in the present work can be confronted to the experimental data at the LHCb, Belle II, CEPC, FCC, ILC in the future, and shed light on the nature of the exotic $X$, $Y$, $Z$ particles. We can take the pole residues as input parameters to study the two-body strong decays of those hidden-bottom tetraquark states

[TABLE]

with the three-point QCD sum rules, and obtain the partial decay widths, and compare them to the experimental data in the future to diagnose the nature of the $Z_{b}$ states, as different quark structures lead to different partial decay widths. In Ref.[22], we tentatively assign the $Z_{c}^{\pm}(3900)$ to be the diquark-antidiquark type axialvector tetraquark state, study the two-body strong decays $Z_{c}^{+}(3900)\to J/\psi\pi^{+}$, $\eta_{c}\rho^{+}$, $D^{+}\bar{D}^{*0}$, $\bar{D}^{0}D^{*+}$ with the QCD sum rules based on solid quark-hadron quality, and obtain the total width of the $Z_{c}^{\pm}(3900)$, which supports assigning the $Z_{c}^{\pm}(3900)$ to be the diquark-antidiquark type axialvector tetraquark state. In Ref.[36], we extend the method to study the two-body strong decays of the $Y(4660)$ as a diquark-antidiquark type vector tetraquark state, and illustrate how to study the relevant hadronic coupling constants based on solid quark-hadron quality. The new method can be applied to study the two-body strong decays of the $Z_{b}$ tetraquark states directly.

The ratios of the partial widths of the decays $Z_{c}(3900/4020)\to\eta_{c}\rho\,,\,J/\psi\pi$ at $\sqrt{s}=4.23\,\rm{GeV}$ measured by the BESIII collaboration are

[TABLE]

at the $90\%$ C.L. [37]. More precise experimental data are still needed to examine the theoretical calculations. As far as the hidden-bottom axialvector tetraquark candidates $Z_{b}(10610)$ and $Z_{b}(10650)$ are concerned, the partial widths have not been measured yet, the experimental data are scarce. More theoretical and experimental works are still needed to assign the $Z_{b}(10610)$ and $Z_{b}(10650)$ unambiguously.

4 Conclusion

In this article, we take the scalar, axialvector and vector bottom (anti)diquark operators as the basic building blocks, and construct the scalar, axialvector and tensor hidden-bottom tetraquark currents to study the mass spectrum of the ground state hidden-bottom tetraquark states with the QCD sum rules in a systematic way by carrying out the operator product expansion up to vacuum condensates of dimension $10$ consistently. In calculations, we use the energy scale formula $\mu=\sqrt{M^{2}_{X/Y/Z}-(2{\mathbb{M}}_{b})^{2}}$ to determine the ideal energy scales of the QCD spectral densities and choose the continuum threshold parameters $\sqrt{s_{0}}=Z_{b}+0.55\pm 0.10\,\rm{GeV}$ as a constraint to extract the masses and pole residues from the QCD sum rules. The predicted masses $10.61\pm 0.09\,\rm{GeV}$ and $10.62\pm 0.09\,\rm{GeV}$ for the $1^{+-}$ tetraquark states supports assigning the $Z_{b}(10610)$ and $Z_{b}(10650)$ to be the axialvector hidden-bottom tetraquark states, more theoretical and experimental works are still needed to assign the $Z_{b}(10610)$ and $Z_{b}(10650)$ unambiguously according to the partial decay widths. The predicted tetraquark masses can be confronted to the experimental data in the future at the LHCb, Belle II, CEPC, FCC, ILC. The pole residues can be taken as input parameters to study the two-body strong decays of those hidden-bottom tetraquark states with the three-point QCD sum rules. Furthermore, we observe that the scalar diquark states $S$ and axialvector diquark states $A$, $\widetilde{A}$ are all good diquark states in building the lowest tetraquark states.

Acknowledgements

This work is supported by National Natural Science Foundation, Grant Number 11775079.