Axialvector tetraquark candidates for the $Z_c(3900)$, $Z_c(4020)$, $Z_c(4430)$, $Z_c(4600)$

TL;DR

This study systematically investigates axialvector tetraquark states with specific quantum numbers using QCD sum rules, identifying candidates for several observed exotic mesons.

Contribution

The paper constructs axialvector and tensor currents and applies QCD sum rules to identify tetraquark candidates for known $Z_c$ states, providing new theoretical insights.

Findings

Identifies one axialvector tetraquark candidate for $Z_c(3900)$ and $Z_c(4430)$.

Finds two candidates for $Z_c(4020)$.

Finds three candidates for $Z_c(4600)$.

Abstract

In this paper, we construct the axialvector and tensor current operators to investigate the ground state tetraquark states and the first radially excited tetraquark states with the quantum numbers $J^{PC}=1^{+-}$ via the QCD sum rules systematically, and observe that there are one axialvector tetraquark candidate for the $Z_c(3900)$ and $Z_c(4430)$, two axialvector tetraquark candidates for the $Z_c(4020)$, three axialvector tetraquark candidates for the $Z_c(4600)$.

Tables

Table 1. Table 1: The Borel parameters, continuum threshold parameters, energy scales of the QCD spectral densities and pole contributions for the QCDSR I.

$Z_c$	$T^2\left({\mathrm{GeV}}^2\right)$	$s_0$	$\mu \left(\mathrm{GeV}\right)$	pole
${\left[uc\right]}_S{\left[\overline{d}\overline{c}\right]}_A-{\left[uc\right]}_A{\left[\overline{d}\overline{c}\right]}_S$	$2.7-3.1$	${\left(4.4\pm 0.1\mathrm{GeV}\right)}^2$	$1.4$	$\left(40-63\right)\%$
${\left[uc\right]}_{\tilde{A}}{\left[\overline{d}\overline{c}\right]}_A-{\left[uc\right]}_A{\left[\overline{d}\overline{c}\right]}_{\tilde{A}}$	$3.2-3.6$	$21.0\pm 1.0{\mathrm{GeV}}^2$	$1.7$	$\left(40-60\right)\%$
${\left[uc\right]}_{\tilde{V}}{\left[\overline{d}\overline{c}\right]}_V+{\left[uc\right]}_V{\left[\overline{d}\overline{c}\right]}_{\tilde{V}}$	$3.7-4.1$	${\left(5.25\pm 0.10\mathrm{GeV}\right)}^2$	$2.9$	$\left(41-60\right)\%$
${\left[uc\right]}_A{\left[\overline{d}\overline{c}\right]}_A$	$3.2-3.6$	$21.0\pm 1.0{\mathrm{GeV}}^2$	$1.7$	$\left(41-61\right)\%$

Table 2. Table 2: The Borel parameters, continuum threshold parameters, energy scales of the QCD spectral densities and pole contributions for the QCDSR II.

$Z_c+Z_c^{\prime }$	$T^2\left({\mathrm{GeV}}^2\right)$	$s_0$	$\mu \left(\mathrm{GeV}\right)$	pole ($Z_c$)
${\left[uc\right]}_S{\left[\overline{d}\overline{c}\right]}_A-{\left[uc\right]}_A{\left[\overline{d}\overline{c}\right]}_S$	$2.7-3.1$	${\left(4.85\pm 0.10\mathrm{GeV}\right)}^2$	$2.6$	$\left(72-88\right)\%$ ($\left(35-52\right)\%$)
${\left[uc\right]}_{\tilde{A}}{\left[\overline{d}\overline{c}\right]}_A-{\left[uc\right]}_A{\left[\overline{d}\overline{c}\right]}_{\tilde{A}}$	$3.2-3.6$	${\left(4.95\pm 0.10\mathrm{GeV}\right)}^2$	$2.8$	$\left(64-80\right)\%$ ($\left(30-44\right)\%$)
${\left[uc\right]}_A{\left[\overline{d}\overline{c}\right]}_A$	$3.2-3.6$	${\left(4.95\pm 0.10\mathrm{GeV}\right)}^2$	$2.8$	$\left(64-81\right)\%$ ($\left(29-43\right)\%$)

Table 3. Table 3: The masses and pole residues of the ground state tetraquark states Z c subscript 𝑍 𝑐 Z_{c} from the QCDSR I, where the superscripts ± plus-or-minus \pm represent the positive parity and negative parity constituents of the tensor diquark states, respectively.

$Z_c$	$|S_{uc},S_{\overline{d}\overline{c}};S⟩$	$M_Z\left(\mathrm{GeV}\right)$	${\lambda }_Z\left({\mathrm{GeV}}^5\right)$
${\left[uc\right]}_S{\left[\overline{d}\overline{c}\right]}_A-{\left[uc\right]}_A{\left[\overline{d}\overline{c}\right]}_S$	$|0,1;1⟩-|1,0;1⟩$	$3.90\pm 0.08$	$\left(2.09\pm 0.33\right)\times {10}^{-2}$
${\left[uc\right]}_{\tilde{A}}{\left[\overline{d}\overline{c}\right]}_A-{\left[uc\right]}_A{\left[\overline{d}\overline{c}\right]}_{\tilde{A}}$	$|1^{+},1;1⟩-|1,1^{+};1⟩$	$4.01\pm 0.09$	$\left(5.96\pm 0.94\right)\times {10}^{-2}$
${\left[uc\right]}_{\tilde{V}}{\left[\overline{d}\overline{c}\right]}_V+{\left[uc\right]}_V{\left[\overline{d}\overline{c}\right]}_{\tilde{V}}$	$|1^{-},1;1⟩+|1,1^{-};1⟩$	$4.66\pm 0.10$	$\left(1.18\pm 0.22\right)\times {10}^{-1}$
${\left[uc\right]}_A{\left[\overline{d}\overline{c}\right]}_A$	$|1,1;1⟩$	$4.00\pm 0.09$	$\left(2.91\pm 0.46\right)\times {10}^{-2}$

Table 4. Table 4: The masses and pole residues of the ground state tetraquark states Z c subscript 𝑍 𝑐 Z_{c} and the first radially excited tetraquark states Z c ′ superscript subscript 𝑍 𝑐 ′ Z_{c}^{\prime} from the QCDSR II.

$Z_c+Z_c^{\prime }$	$M_Z\left(\mathrm{GeV}\right)$	${\lambda }_Z\left({\mathrm{GeV}}^5\right)$	$M_{Z^{\prime }}\left(\mathrm{GeV}\right)$	${\lambda }_{Z^{\prime }}\left({\mathrm{GeV}}^5\right)$
${\left[uc\right]}_S{\left[\overline{d}\overline{c}\right]}_A-{\left[uc\right]}_A{\left[\overline{d}\overline{c}\right]}_S$	$3.81$	$1.77\times {10}^{-2}$	$4.47\pm 0.09$	$\left(6.02\pm 0.80\right)\times {10}^{-2}$
${\left[uc\right]}_{\tilde{A}}{\left[\overline{d}\overline{c}\right]}_A-{\left[uc\right]}_A{\left[\overline{d}\overline{c}\right]}_{\tilde{A}}$	$3.78$	$3.94\times {10}^{-2}$	$4.60\pm 0.09$	$\left(1.35\pm 0.18\right)\times {10}^{-1}$
${\left[uc\right]}_A{\left[\overline{d}\overline{c}\right]}_A$	$3.73$	$1.76\times {10}^{-2}$	$4.58\pm 0.09$	$\left(6.55\pm 0.85\right)\times {10}^{-2}$

Equations

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

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

$$J_{\mu }^1(x)$$

$$J_{\mu }^2(x)$$

$$J_{\mu }^3(x)$$

$$J_{\mu \nu }(x)$$

$$\hat{C}{J}_{\mu }(x){\hat{C}}^{-1}$$

$$\hat{C}{J}_{\mu \nu }(x){\hat{C}}^{-1}$$

$$\hat{P}{J}_{\mu }(x){\hat{P}}^{-1}$$

$$\hat{P}{J}_{\mu \nu }(x){\hat{P}}^{-1}$$

$$\hat{P}{\epsilon }^{abc}{q}_b^T(x)C{\sigma }_{jk}{\gamma }_5{Q}_c(x){\hat{P}}^{-1}$$

$$\hat{P}{\epsilon }^{abc}{q}_b^T(x)C{\sigma }_{0j}{Q}_c(x){\hat{P}}^{-1}$$

$$\hat{P}{\epsilon }^{abc}{q}_b^T(x)C{\sigma }_{0j}{\gamma }_5{Q}_c(x){\hat{P}}^{-1}$$

$$\hat{P}{\epsilon }^{abc}{q}_b^T(x)C{\sigma }_{jk}{Q}_c(x){\hat{P}}^{-1}$$

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

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

$$⟨0|J_{\mu }(0)|Z_c(p)⟩$$

$$⟨0|J_{\mu \nu }(0)|Z_c(p)⟩$$

$$⟨0|J_{\mu \nu }(0)|Y(p)⟩$$

$$\sum _{\lambda }{\epsilon }_{\mu }^{\ast }(\lambda ,p){\epsilon }_{\nu }(\lambda ,p)$$

$$J_{\mu }^1(x)$$

$$J_{\mu }^1(x,ϵ)$$

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

$${\Sigma }_{D{D}^{\ast }}(p)$$

$${\Sigma }_{J/\psi \pi }(p)$$

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

$${\lambda }_Z^2\exp \left(-\frac{M_Z^2}{T^2}\right)={\int }_{4m_c^2}^{s_0}ds\rho (s)\exp \left(-\frac{s}{T^2}\right),$$

$$\rho (s)$$

$${\rho }_3(s)$$

$${\rho }_4(s)$$

$${\rho }_5(s)$$

$${\rho }_6(s)$$

$${\rho }_7(s)$$

$${\rho }_8(s)$$

$${\rho }_{10}(s)$$

$$M_Z^2$$

$${\lambda }_Z^2\exp \left(-\frac{M_Z^2}{T^2}\right)+{\lambda }_{Z^{\prime }}^2\exp \left(-\frac{M_{Z^{\prime }}^2}{T^2}\right)$$

$${\lambda }_1^2\exp \left(-\tau M_1^2\right)+{\lambda }_2^2\exp \left(-\tau M_2^2\right)$$

Full text

PACS number: 12.39.Mk, 12.38.Lg

Key words: Tetraquark state, QCD sum rules

1 Introduction

In 2019, the LHCb collaboration performed an angular analysis of the weak decays $B^{0}\to J/\psi K^{+}\pi^{-}$ using proton-proton collision data, examined the $m(J/\psi\pi^{-})$ versus the $m(K^{+}\pi^{-})$ plane, and observed two possible resonant structures in the vicinity of the energies $m(J/\psi\pi^{-})=4200\,\rm{MeV}$ and $4600\,\rm{MeV}$, respectively [1]. There have been two tentative assignments of the structure $Z_{c}(4600)$ in the vicinity of $m(J/\psi\pi^{-})=4600\,\rm{MeV}$, the $[dc]_{P}[\bar{u}\bar{c}]_{A}-[dc]_{A}[\bar{u}\bar{c}]_{P}$ type vector tetraquark state with $J^{PC}=1^{--}$ [2] and the first radially excited $[dc]_{T}[\bar{u}\bar{c}]_{A}-[dc]_{A}[\bar{u}\bar{c}]_{T}$ type tetraquark state with $J^{PC}=1^{+-}$ [3]. In this paper, we use the subscripts $P$, $S$, $V$, $A$ and $T$ to represent the pseudoscalar, scalar, vector, axialvector and tensor color-antitriplet diquark states, respectively.

In 2013, the BESIII collaboration observed the charged charmonium-like resonance $Z^{\pm}_{c}(3900)$ in the $\pi^{\pm}J/\psi$ invariant mass spectrum in the process $e^{+}e^{-}\to J/\psi\pi^{+}\pi^{-}$ with $M_{Z_{c}}=(3899.0\pm 3.6\pm 4.9)\,\rm{MeV}$ and $\Gamma_{Z_{c}}=(46\pm 10\pm 20)\,\rm{MeV}$, respectively [4]. The Belle collaboration also observed the $Z^{\pm}_{c}(3900)$ in the same process [5], furthermore, the CLEO collaboration confirmed the existence of the $Z_{c}^{\pm}(3900)$ [6]. Almost at the same time, the BESIII collaboration observed the charmonium-like resonance $Z^{\pm}_{c}(4025)$ in the vicinity of the threshold $(D^{*}\bar{D}^{*})^{\pm}$ in the electron-positron scattering process $e^{+}e^{-}\to(D^{*}\bar{D}^{*})^{\pm}\pi^{\mp}$ [7]. Moreover, the BESIII collaboration observed the charmonium-like resonance $Z_{c}^{\pm}(4020)$ in the $\pi^{\pm}h_{c}$ invariant mass spectrum in the electron-positron collisions $e^{+}e^{-}\to\pi^{+}\pi^{-}h_{c}$ [8]. Now the $Z^{\pm}_{c}(4020)$ and $Z^{\pm}_{c}(4025)$ are listed in The Review of Particle Physics as the same particle [9]. In 2014, the LHCb collaboration performed a four-dimensional fit of the scattering amplitude for the decay $B^{0}\to\psi^{\prime}\pi^{-}K^{+}$ in proton-proton collisions, and obtained the first independent confirmation of the charmonium-like resonance $Z^{-}_{c}(4430)$ and determined its quantum numbers to be $J^{P}=1^{+}$ [10]. In 2017, the BESIII collaboration established the charmonium-like resonance $Z_{c}(3900)$’s quantum numbers to be $J^{P}=1^{+}$ [11].

There have been several possible explanations for the exotic states $Z_{c}(3900)$ and $Z_{c}(4020)$, such as the molecular states (from the heavy quark symmetries [12, 13], the QCD sum rules [14, 15], the light-front quark model [16], the one-pion exchange model [17], and the phenomenological Lagrangian approach [18]), the tetraquark states (from the diquark model with the effective Hamiltonian [19], the QCD sum rules [20, 21, 22, 23], and the potential model [24]), the triangle singularities (in the rescattering amplitudes) [25], the threshold effects [26], etc.

We can tentatively assign the hidden-charm resonances $Z_{c}(3900)$ and $Z_{c}(4430)$ as the ground state tetraquark state and the first radially excited tetraquark state respectively considering the similar decays, $Z_{c}(3900)^{\pm}\to J/\psi\pi^{\pm}$, $Z_{c}(4430)^{\pm}\to\psi^{\prime}\pi^{\pm}$, and the almost equal energy gaps $M_{Z_{c}(4430)}-M_{Z_{c}(3900)}=591\,\rm{MeV}$, $M_{\psi^{\prime}}-M_{J/\psi}=589\,\rm{MeV}$ [27, 28]. In Ref.[29], we adopt the method invented in Ref.[30] for the conventional quarkonium to study the $Z^{\pm}_{c}(3900)$ as the ground state axialvector tetraquark state and the $Z^{\pm}_{c}(4430)$ as the first radially excited axialvector tetraquark state respectively, and employ the energy scale formula $\mu=\sqrt{M^{2}_{X/Y/Z}-(2{\mathbb{M}}_{c})^{2}}$ to select the best energy scales of the spectral densities at the QCD side of the QCD sum rules with the effective (or constituent) charm quark mass ${\mathbb{M}}_{c}$ [31]. In Ref.[32], this subject is studied with the QCD sum rules adopting another parameter system. In Refs.[22, 23], we observe that we can assign the $Z_{c}(4020/4025)$ to be the ground state axialvector $[uc]_{A}[\bar{d}\bar{c}]_{A}$ tetraquark state with the quantum numbers $J^{PC}=1^{+-}$ according to the QCD sum rules calculations. If the $Z_{c}(4600)$ is the first radial excitation of the hidden-charm tetraquark candidate $Z_{c}(4020/4025)$, its preferred decay mode is $Z_{c}(4600)\to\psi^{\prime}\pi$ rather than $Z_{c}(4600)\to J/\psi\pi$.

In this paper, we intend to perform a detailed and updated analysis of the ground states and the first radially excited states of the charged hidden-charm tetraquark states with the QCD sum rules, and explore the possible assignments of the $Z_{c}(4600)$ state in the scenario of the axialvector hidden-charm tetraquark states with the quantum numbers $J^{PC}=1^{+-}$.

The paper is arranged as follows: in Sect.2 we obtain the analytical expressions of the QCD sum rules for the hidden-charm axialvector tetraquark states $Z_{c}$; in Sect.3 we provide the numerical results for the masses and pole residues of the $Z_{c}$ states and detailed discussions; in Sect.4 we reach the conclusion.

2 The QCD sum rules for the axialvector tetraquark states

In the following, let us write down the two-point Green functions (or correlation functions) $\Pi_{\mu\nu}(p)$ and $\Pi_{\mu\nu\alpha\beta}(p)$ as the first step,

[TABLE]

where the four-quark current operators $J_{\mu}(x)=J^{1}_{\mu}(x)$, $J^{2}_{\mu}(x)$, $J^{3}_{\mu}(x)$,

[TABLE]

the superscripts $i$, $j$, $k$, $m$ and $n$ are color indexes with values obeying the antisymmetric tensor $\varepsilon$, the charge conjugation matrix $C=i\gamma^{2}\gamma^{0}$. If we preform charge conjugation (and parity) transform $\widehat{C}$ (and $\widehat{P}$), the axialvector current operators $J_{\mu}(x)$ and tensor current operator $J_{\mu\nu}(x)$ have the following properties,

[TABLE]

where the coordinates $x^{\mu}=(t,\vec{x})$ and $\tilde{x}^{\mu}=(t,-\vec{x})$.

The diquark operators $\varepsilon^{ijk}q^{T}_{j}C\Gamma Q_{k}$ in the attractive color-antitriplet $\bar{3}_{c}$ channel have five spinor structures, 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}$ correspond to the scalar, pseudoscalar, vector, axialvector and tensor diquark operators, respectively. The favorable quark-quark correlations are the scalar diquark and axialvector diquark in the color-antitriplet $\bar{3}_{c}$ channel from the QCD sum rules [33]. If we introduce a relative P-wave between the light quark and heavy quark, we can obtain the pseudoscalar diquark operator $\varepsilon^{ijk}q^{T}_{j}C\gamma_{5}\underline{\gamma_{5}}Q_{k}$ and vector diquark operator $\varepsilon^{ijk}q^{T}_{j}C\gamma_{\mu}\underline{\gamma_{5}}Q_{k}$ without introducing the additional P-wave explicitly, as multiplying a $\gamma_{5}$ can change the parity, the P-wave effect is embodied in the underlined $\gamma_{5}$. The pseudoscalar diquark states and vector diquark states (or the P-wave diquark states) have larger masses compared with the scalar diquark states and axialvector diquark states, we would like to choose the scalar diquark and axialvector diquark to construct the four-quark current operators to interpolate the lower tetraquark states.

The tensor heavy diquark operators $\varepsilon^{abc}q^{T}_{b}(x)C\sigma_{\mu\nu}\gamma_{5}Q_{c}(x)$ and $\varepsilon^{abc}q^{T}_{b}(x)C\sigma_{\mu\nu}Q_{c}(x)$ have both axialvector constituents and vector constituents,

[TABLE]

where the space indexes $j$, $k=1$, $2$, $3$. The tensor diquark operators also play an important role in constructing the tetraquark current operators [34]. We multiply the tensor diquark (antidiquark) operators with the axialvector or vector antidiquark (diquark) operators so as to project out the axialvector constituents and vector constituents to construct the four-quark current operators $J^{2}_{\mu}(x)$ and $J^{3}_{\mu}(x)$. Thereafter, we will use the $\tilde{V}$ and $\tilde{A}$ to represent the vector constituent and axialvector constituent of the tensor diquark operators, respectively.

The four-quark current operators $J^{1}_{\mu}(x)$, $J^{2}_{\mu}(x)$ and $J^{3}_{\mu}(x)$ couple potentially to the $[uc]_{S}[\bar{d}\bar{c}]_{A}-[uc]_{A}[\bar{d}\bar{c}]_{S}$ type, $[uc]_{\tilde{A}}[\bar{d}\bar{c}]_{A}-[uc]_{A}[\bar{d}\bar{c}]_{\tilde{A}}$ type and $[uc]_{\tilde{V}}[\bar{d}\bar{c}]_{V}+[uc]_{V}[\bar{d}\bar{c}]_{\tilde{V}}$ type axialvector hidden-charm tetraquark states with the spin-parity-charge-conjugation $J^{PC}=1^{+-}$, respectively. While the current operator $J_{\mu\nu}(x)$ couples potentially to both the $[uc]_{A}[\bar{d}\bar{c}]_{A}$ type axialvector tetraquark state with $J^{PC}=1^{+-}$ and vector tetraquark state with $J^{PC}=1^{--}$. Thereafter, we will not distinguish the negative or positive electric-charge of the $Z_{c}$ tetraquark states, as they have degenerate masses.

We insert a complete set of hadron states which have nonvanishing couplings with the four-quark current operators $J_{\mu}(x)$ and $J_{\mu\nu}(x)$ into the Green functions $\Pi_{\mu\nu}(p)$ and $\Pi_{\mu\nu\alpha\beta}(p)$ to get the hadronic representation [35, 36]. Then we separate the ground state axialvector and vector tetraquark state contributions from other contributions, such as the higher excited tetraquark states and continuum states, to obtain the results,

[TABLE]

where the $Z$ represents the axialvector tetraquark states, the $Y$ represents the vector tetraquark states, the $\lambda_{Z}$, $\widetilde{\lambda}_{Z}$ and $\widetilde{\lambda}_{Y}$ are the pole residues or current-tetraquark coupling constants,

[TABLE]

the antisymmetric tensor $\varepsilon_{0123}=-1$, the $\varepsilon_{\mu}(\lambda,p)$ are the polarization vectors of the axialvector and vector tetraquark states satisfy the summation formula,

[TABLE]

The diquark-antidiquark type four-quark current operators $J_{\mu}(x)$ and $J_{\mu\nu}(x)$ couple potentially to the diquark-antidiquark type hidden-charm tetraquark states. We can perform Fierz rearrangements to those currents in both the spinor space and color space to obtain a series of color-singlet-color-singlet (or meson-meson) type current operators, for example,

[TABLE]

while the constituents, such as $\bar{c}(x)i\gamma_{5}c(x)\,\bar{d}(x)\gamma^{\mu}u(x)$, $\bar{c}(x)\gamma^{\mu}c(x)\,\bar{d}(x)i\gamma_{5}u(x)$, etc, couple potentially to the meson-meson type scattering states or tetraquark molecular states.

However, we should be careful in performing the Fierz rearrangements, the rearrangements in the spinor space and color space are quite non-trivial, the scenarios of the diquark-antidiquark type tetraquark states and meson-meson type molecular states are quite different.

According to the arguments of Selem and Wilczek, a diquark-antidiquark type tetraquark can be plausibly described by two diquarks trapped in a double potential well, the two potential wells are separated apart by a barrier [37]. At long distances, the diquark and antiquark serve as point color charges respectively, and attract each other strongly just like in the quark and antiquark bound states. However, when the two diquarks approach each other, the attractions between quark and antiquark in different diquarks decrease the bound energy of the diquarks and tend to destroy the diquarks. Those effects (beyond the naive one-gluon exchange force) increase when the distance between the diquark and antidiquark decreases, and a repulsive interaction between the diquark and antidiquark emerges, if large enough, it will lead to a barrier between the diquark and antidiquark [38]. The two potential wells which are separated apart by a barrier can give successful descriptions of the diquark-antidiquark type tetraquark states [38].

While in the dynamical picture of the tetraquark states, the large spatial separation between the diquark and antidiquark leads to small wave-function overlap between the quark-antiquark pair [39], the rearrangements in the spinor space and color space are highly suppressed.

In practical calculations, it is difficult to account for the non-local effects between the diquark and antidiquark pair in the four-quark currents $J_{\mu}(x)$ and $J_{\mu\nu}(x)$ directly, for example, the current $J^{1}_{\mu}(x)$ can be modified to

[TABLE]

to account for the non-locality by adding a finite $\epsilon$, however, it is very difficult to deal with the finite $\epsilon$ both at the hadron side and at the QCD side in a consistent way, we expend the current $J^{1}_{\mu}(x,\epsilon)$ in terms of Taylor series of $\epsilon$,

[TABLE]

then expend the correlation function $\Pi_{\mu\nu}(p)$ also in terms of Taylor series of $\epsilon$,

[TABLE]

where the components $\Pi_{\mu\nu}\left(\mathcal{O}(\epsilon^{i})\right)$ with $i=0$, $1$, $2$, $\cdots$ stand for the contributions of the order $\mathcal{O}(\epsilon^{i})$. In this article, we study the leading order contributions $J^{1}_{\mu}(x)=J^{1}_{\mu}(x,0)$ and $\Pi_{\mu\nu}(p)=\Pi_{\mu\nu}\left(\mathcal{O}(\epsilon^{0})\right)$, the effects beyond the leading order frustrate the Fierz rearrangements of the diquark-antidiquark type currents into a series of color-singlet-color-singlet (meson-meson) type currents freely.

We can use the Feynman diagram drawn in Fig.1 to describe the lowest order contributions in the correlation functions for the diquark-antidiquark type four-quark currents, and use the Feynman diagrams drawn in Fig.2 to describe the corresponding lowest order contributions in the correlation functions for the color-singlet-color-singlet type four-quark currents. The Feynman diagram drawn in Fig.1 cannot be factorized into the two Feynman diagrams drawn in Fig.2 freely due to the barrier (or spatial separation) between the diquark and antidiquark [38, 39]. When a quark (antiquark) in the diquark (antidiquark) penetrates the barrier, the Feynman diagram drawn in Fig.1 is factorizable in color space. In this case, the non-factorizable diagrams start at the order $\mathcal{O}(\alpha_{s}^{2})$ [40].

In Ref.[40], Lucha, Melikhov, and Sazdjiand argue that the diquark-antidiquark type four-quark currents can be changed into the color-singlet-color-singlet (meson-meson) type currents through Fierz transformation, the Feynman diagrams which make contributions to the quark-gluon operators of the order $\mathcal{O}(\alpha_{s}^{0})$ and $\mathcal{O}(\alpha_{s})$ in accomplishing the operator product expansion are factorizable and are canceled out by the contributions of the two-meson scattering states at the phenomenological side, furthermore, the factorizable parts (in color space) of the Feynman diagrams of the order $\mathcal{O}(\alpha_{s}^{2})$ are also canceled out by the contributions come from the two-meson scattering states (or more precisely, the free two-meson states), the relevant non-factorizable contributions start at the order $\mathcal{O}(\alpha_{s}^{2})$. We do not agree with their viewpoint, as there exists a repulsive barrier [37, 38] or a large spatial separation [39], which are embodied in the non-local effects, to prevent performing the Fierz transformation freely, although at the present time we cannot take into account the non-local effects in the QCD sum rules, and have to take the leading order approximations $J^{1}_{\mu}(x)=J^{1}_{\mu}(x,0)$ and $\Pi_{\mu\nu}(p)=\Pi_{\mu\nu}\left(\mathcal{O}(\epsilon^{0})\right)$. Our viewpoint is that the relevant contributions begin at the order $\mathcal{O}(\alpha_{s}^{0})$, it is not necessary to perform or it is very difficult to perform the Fierz transformation to separate the factorizable and non-factorizable contributions in the color space, we should take into account both the factorizable and nonfactorizable Feynman diagrams for the diquark-antidiquark type currents.

When the quark or antiquark penetrates the barrier, we can perform the Fierz rearrangements, and study the effects of the scattering states. Now let us begin to explore the contributions of the meson-meson type scattering states (in other words, the two-meson loops) to the Green function $\Pi_{\mu\nu}(p)$ for the four-quark current $J^{1}_{\mu}(x)$ as a representative example,

[TABLE]

where

[TABLE]

$\widetilde{g}_{\mu\nu}(p)=-g_{\mu\nu}+\frac{p_{\mu}p_{\nu}}{p^{2}}$, the $G_{ZDD^{*}}$ and $G_{ZJ/\psi\pi}$ are the hadronic coupling constants. We resort to the bare quantities $\widehat{\lambda}_{Z}$ and $\widehat{M}_{Z}$ so as to absorb the divergent terms which appear in the integrals in calculating the self-energies $\Sigma_{DD^{*}}(p)$, $\Sigma_{J/\psi\pi}(p)$, etc. The self-energies after renormalization result in a finite energy-dependent width to modify the dispersion relation,

[TABLE]

the experimental value of the total decay width $\Gamma_{Z_{c}(3900)}(M_{Z}^{2})=(46\pm 10\pm 20)\,\rm{MeV}$ [4] (or $(28.2\pm 2.6)\,\rm{MeV}$ [9]), the zero width approximation in the spectral densities at the phenomenological side are approved reasonable [41]. In this paper, we neglect the contributions of the meson-meson type scattering states or the two-meson loops, the predictions are still robust.

We calculate all the Feynman diagrams in performing the operator product expansion to obtain the QCD spectral representation of the Green functions (or correlation functions) $\Pi_{\mu\nu}(p)$ and $\Pi_{\mu\nu\alpha\beta}(p)$. In analytical calculations, we take into account the vacuum condensates (by selecting the quark-gluon operators of the orders $\mathcal{O}(\alpha_{s}^{k})$ with $k\leq 1$) up to dimension $10$ consistently and factorize the higher dimensional condensates into the lower dimensional condensates by assuming vacuum saturation. After getting the analytical expressions of the Green functions at the quark-gluon level, we obtain the spectral representation via dispersion relation. Now we match the hadronic representation with the QCD representation of the Green functions (or correlation functions) $\Pi_{Z}(p^{2})$ below the continuum threshold parameters $s_{0}$ and carry out the Borel transformation in regard to $P^{2}=-p^{2}$ to get the QCD sum rules:

[TABLE]

[TABLE]

$\lambda_{Z}=\widetilde{\lambda}_{Z}M_{Z}$, the $T^{2}$ is the Borel parameter, the subscripts $i$ in the components of the QCD spectral densities $\rho_{i}(s)$ represent the dimensions of the vacuum condensates,

[TABLE]

We neglect the cumbersome analytical expressions of the spectral densities at the quark-gluon level in saving printed pages. We can refer to Ref.[20] for the technical details in calculating the Feynman diagrams. On the other hand, we can refer to Refs.[20, 23] for the explicit expressions of the spectral densities at the quark-gluon level for the axialvector current $J^{1}_{\mu}(x)$ and tensor current $J_{\mu\nu}(x)$. In this work, we recalculate those QCD spectral densities, and use the formula $t^{a}_{ij}t^{a}_{mn}=-\frac{1}{6}\delta_{ij}\delta_{mn}+\frac{1}{2}\delta_{jm}\delta_{in}$ with $t^{a}=\frac{\lambda^{a}}{2}$ to deal with the higher dimensional vacuum condensates, where the $\lambda^{a}$ is the Gell-Mann matrix. This routine leads to slight but neglectful differences compared to the old calculations. For the currents $J^{2}_{\mu}(x)$ and $J^{3}_{\mu}(x)$, we neglect the tiny contributions of the $4\pi\alpha_{s}\langle\bar{q}q\rangle^{2}$, which originate from the operators like $\langle\bar{q}_{j}\gamma_{\mu}q_{i}g_{s}D_{\nu}G^{a}_{\alpha\beta}t^{a}_{mn}\rangle$.

We derive Eq.(16) in regard to $\tau=\frac{1}{T^{2}}$, then reach the QCD sum rules for the tetraquark masses by eliminating the pole residues $\lambda_{Z}$ through a fraction,

[TABLE]

Thereafter, we will refer the QCD sum rules in Eq.(16) and Eq.(19) as QCDSR I.

If we take into account the contributions of the first radially excited tetraquark states $Z_{c}^{\prime}$ in the hadronic representation, we can obtain the QCD sum rules,

[TABLE]

where the $s_{0}^{\prime}$ is continuum threshold parameter, then we introduce the notations $\tau=\frac{1}{T^{2}}$, $D^{n}=\left(-\frac{d}{d\tau}\right)^{n}$, and resort to the subscripts $1$ and $2$ to represent the ground state tetraquark state $Z_{c}$ and the first radially excited tetraquark state $Z_{c}^{\prime}$ respectively for simplicity. We rewrite the QCD sum rules as

[TABLE]

here we introduce the subscript $QCD$ to represent the QCD representation. We derive the QCD sum rules in Eq.(21) in regard to $\tau$ to get

[TABLE]

From Eqs.(21)-(22), we obtain the QCD sum rules,

[TABLE]

where the indexes $i\neq j$. Let us derive the QCD sum rules in Eq.(23) in regard to $\tau$ to get

[TABLE]

The squared masses $M_{i}^{2}$ satisfy the equation,

[TABLE]

where

[TABLE]

the indexes $i=1,2$ and $j,k=0,1,2,3$. Finally we solve above equation analytically to obtain two solutions [30],

[TABLE]

[TABLE]

From now on, we will denote the QCD sum rules in Eq.(20) and Eqs.(27)-(28) as QCDSR II. In calculations, we observe that if we specify the energy scales of the spectral densities in the QCD representation, only one solution satisfies the energy scale formula $\mu=\sqrt{M^{2}_{X/Y/Z}-(2{\mathbb{M}}_{c})^{2}}$ in the QCDSR II, we have to abandon the other solution. In this paper, we retain the mass $M_{2}$ ($M_{Z^{\prime}}$) and discard the mass $M_{1}$ ($M_{Z}$).

3 Numerical results and discussions

We take the standard values or conventional 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 typical energy scale $\mu=1\,\rm{GeV}$ [35, 36, 42], and take the modified minimal subtraction mass of the charm quark $m_{c}(m_{c})=(1.275\pm 0.025)\,\rm{GeV}$ from the Particle Data Group [9]. We should evolve the quark condensate, mixed quark condensate and modified minimal subtraction mass to a special energy scale to warrant the parameters in the QCD spectral densities having the same energy scale. Now let us account for the energy-scale dependence of the input parameters at the quark-gluon level,

[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}}$, with the values $\Lambda=210\,\rm{MeV}$, $292\,\rm{MeV}$ and $332\,\rm{MeV}$ for the quark flavors $n_{f}=5$, $4$ and $3$, respectively [9, 43]. As we explore the hidden-charm tetraquark states, we choose the flavor $n_{f}=4$ and search for the best energy scales $\mu$.

The Okubo-Zweig-Iizuka supper-allowed decays

[TABLE]

are expected to take place easily. The energy gaps maybe have the relations $M_{Z^{\prime}}-M_{Z}=m_{\psi^{\prime}}-m_{J/\psi}$ and $M_{Z^{\prime\prime}}-M_{Z^{\prime}}=m_{\psi^{\prime\prime}}-m_{\psi^{\prime}}$. The charmonium masses are $m_{J/\psi}=3.0969\,\rm{GeV}$, $m_{\psi^{\prime}}=3.686097\,\rm{GeV}$ and $m_{\psi^{\prime\prime}}=4.039\,\rm{GeV}$ from the Particle Data Group [9], $m_{\psi^{\prime}}-m_{J/\psi}=0.59\,\rm{GeV}$, $m_{\psi^{\prime\prime}}-m_{J/\psi}=0.94\,\rm{GeV}$, we can choose the continuum threshold parameters to be $\sqrt{s_{0}}=M_{Z}+0.59\,\rm{GeV}$ and $\sqrt{s_{0}^{\prime}}=M_{Z}+0.95\,\rm{GeV}$ tentatively and vary the continuum threshold parameters and Borel parameters to satisfy the following four criteria:

$\bf 1.$ The ground state tetraquark state or single-pole term makes dominant contribution at the hadron side;

$\bf 2.$ The operator product expansion is convergent below the continuum thresholds, and the higher dimensional vacuum condensates make minor contribution;

$\bf 3.$ The Borel platforms appear in both the lineshapes of the tetraquark masses and pole residues with variations of the Borel parameters;

$\bf 4.$ The masses of the tetraquark states satisfy the energy-scale formula.

After trial and error, we reach the feasible continuum threshold parameters and Borel windows, we also acquire the best energy scales of the spectral densities at the quark-gluon level and the contributions of the ground state tetraquark states for the QCDSR I, see Table 1. In general, for the continuum threshold parameters $s_{0}$, we can take any values satisfy the relation $M_{gr}<\sqrt{s_{0}}\leq M_{gr}+\Delta$, where the subscript $gr$ denotes the ground states, as there exists an energy gap $\Delta$ between the ground state and the first radial excited state. For the conventional S-wave quark-antiquark mesons, the energy gaps $\Delta$ vary from $m_{m_{K^{*}(1410)}}-m_{K^{*}(892)}=522\,\rm{MeV}$ to $m_{\pi(1300)}-m_{\pi}=1160\,\rm{MeV}$, i.e. $\Delta=522\sim 1160\,\rm{MeV}$ [9]. In the QCD sum rules for the conventional quark-antiquark mesons, we usually choose the values $\sqrt{s_{0}}=M_{gr}+(0.4\sim 0.7)\,\rm{GeV}$ [42]. In Table 1, the continuum threshold parameters $s_{0}$ satisfy the relation $\sqrt{s_{0}}=M_{Z_{c}}+(0.4\sim 0.6)\,\rm{GeV}$ or $M_{Z_{c}}+(0.5\sim 0.7)\,\rm{GeV}$, it is reasonable, as the values $\exp\left(-s^{0}_{max}/T^{2}_{max}\right)=(1\sim 2)\%$, where the subscript $max$ denotes the maximum values, the contributions of the $Z_{c}^{\prime}$ are greatly suppressed if there are any. In Table 1, we write the continuum threshold parameters as $s_{0}=21.0\pm 1.0\,\rm{GeV}^{2}$ rather than as $s_{0}=(4.58\pm 0.11\,\rm{GeV})^{2}$ for the $[uc]_{\tilde{A}}[\bar{d}\bar{c}]_{A}-[uc]_{A}[\bar{d}\bar{c}]_{\tilde{A}}$ and $[uc]_{A}[\bar{d}\bar{c}]_{A}$ tetraquark states to retain the same form as in our previous work [23]. In Ref.[23], we study the axialvector $[uc]_{A}[\bar{d}\bar{c}]_{A}$ tetraquark state and choose the continuum threshold parameters as $s_{0}=21.0\pm 1.0\,\rm{GeV}^{2}$.

Again we obtain the corresponding parameters for the QCDSR II using trial and error, see Table 2. In this paper, we employ the energy scale formula $\mu=\sqrt{M^{2}_{X/Y/Z}-(2{\mathbb{M}}_{c})^{2}}$ with the effective charm quark mass (or constituent charm quark mass) ${\mathbb{M}}_{c}$ to restrain the tetraquark masses and the energy scales of the spectral densities [31]. The energy scale formula can enhance the contributions of the ground state tetraquark states remarkably at the hadron representation and improve the convergent behaviors of the operator product expansion remarkably at the QCD representation by augmenting the contributions of the lower dimensional vacuum condensates, and is feasible for the hidden-charm tetraquark states and hidden-charm pentaquark states [44].

From Table 1 and Table 2, we can see that the contributions of the single-pole terms (the ground state tetraquark states) are about $(40-60)\%$ for the QCDSR I, the contributions of the two-pole terms (the ground state tetraquark states plus the first radially excited tetraquark states) are about $(70-80)\%$ for the QCDSR II, which satisfy the pole dominance criterion very well. In the QCDSR II, the contributions of the ground state tetraquark states are about $(30-45)\%$, which are much less than the corresponding ground state tetraquark contributions in the QCDSR I, for the ground state tetraquark masses and pole residues, we prefer the predictions from the QCDSR I. In numerical calculations, we find that the contributions of the vacuum condensates of dimension 10 (the largest dimension) are of percent level at the Borel widows for both the QCDSR I and QCDSR II, the minor contributions warrant good convergent behaviors of the operator product expansion.

Now let us take into account all the uncertainties of input parameters, and reach the numerical values of the masses and pole residues of the ground state tetraquark states $Z_{c}$ and the first radially excited tetraquark states $Z_{c}^{\prime}$, which are shown in Table 3 and Table 4. From those Tables, we can see that the ground state tetraquark masses from the QCDSR I and the radially excited tetraquark masses from the QCDSR II satisfy the energy scale formula $\mu=\sqrt{M^{2}_{X/Y/Z}-(2{\mathbb{M}}_{c})^{2}}$, where the updated value of the effective charm quark mass (or constituent charm quark mass) ${\mathbb{M}}_{c}=1.82\,\rm{GeV}$ is adopted [23]. In Table 4, we also present the central values of the ground state tetraquark masses and pole residues extracted from the QCDSR II at the ideal energy scales which are shown in Table 1. We examine Table 4 and observe that the ground state tetraquark masses cannot satisfy the energy scale formula, so we will discard those values. This is the shortcoming of the QCDSR II.

In Fig.3, we plot the ground state tetraquark masses from the QCDSR I and the first radially excited tetraquark masses from the QCDSR II in regard to variations of the Borel parameters in much larger regions than the Borel windows, which are shown in Table 1 and Table 2. From the Fig.3, we find that there indeed appear very flat platforms in the Borel windows for the $[uc]_{S}[\bar{d}\bar{c}]_{A}-[uc]_{A}[\bar{d}\bar{c}]_{S}$ type, $[uc]_{\tilde{A}}[\bar{d}\bar{c}]_{A}-[uc]_{A}[\bar{d}\bar{c}]_{\tilde{A}}$ type and $[uc]_{A}[\bar{d}\bar{c}]_{A}$ type axialvector tetraquark states. For the $[uc]_{\tilde{V}}[\bar{d}\bar{c}]_{V}+[uc]_{V}[\bar{d}\bar{c}]_{\tilde{V}}$ type tetraquark state, we only plot the ground state tetraquark mass, as the ground state tetraquark mass is large enough. From the Fig.3, we also find that the platform in the Borel window is not flat enough, at the region $T^{2}<3.6\,\rm{GeV}^{2}$, the mass increases quickly and monotonously along with the increase of the value of Borel parameter, the platform appears approximately only at the region $T^{2}>3.6\,\rm{GeV}^{2}$.

The predicted mass $M_{Z}=3.90\pm 0.08\,\rm{GeV}$ for the ground state tetraquark state $[uc]_{S}[\bar{d}\bar{c}]_{A}-[uc]_{A}[\bar{d}\bar{c}]_{S}$ exhibits very good agreement with the experimental value $M_{Z(3900)}=(3899.0\pm 3.6\pm 4.9)\,\rm{MeV}$ from the BESIII collaboration [4], which is in favor of assigning the $Z_{c}(3900)$ to be the ground state tetraquark state $[uc]_{S}[\bar{d}\bar{c}]_{A}-[uc]_{A}[\bar{d}\bar{c}]_{S}$ with the quantum numbers $J^{PC}=1^{+-}$ [20]. In Ref.[45], we study the non-leptonic decays $Z_{c}^{+}(3900)\to J/\psi\pi^{+}$, $\eta_{c}\rho^{+}$, $D^{+}\bar{D}^{*0}$, $\bar{D}^{0}D^{*+}$ with the three-point QCD sum rules. In analytical calculations, we take into account both the factorizable and nonfactorizable Feynman diagrams, match the hadronic representation with the QCD representation according to solid quark-hadron duality, and get the total decay width $\Gamma_{Z_{c}}=54.2\pm 29.8\,\rm{MeV}$, which agrees with the experimental value $(46\pm 10\pm 20)\,\rm{MeV}$ very good considering the uncertainties [4].

The predicted mass $M_{Z}=4.47\pm 0.09\,\rm{GeV}$ for the first radially excited tetraquark state $[uc]_{S}[\bar{d}\bar{c}]_{A}-[uc]_{A}[\bar{d}\bar{c}]_{S}$ exhibits very good agreement with the experimental value $M_{Z(4430)}=(4475\pm 7\,{{}_{-25}^{+15}})\,\rm{MeV}$ from the LHCb collaboration [10], which is in favor of assigning the $Z_{c}(4430)$ to be the first radially excited tetraquark state $[uc]_{S}[\bar{d}\bar{c}]_{A}-[uc]_{A}[\bar{d}\bar{c}]_{S}$ with the quantum numbers $J^{PC}=1^{+-}$. We can investigate its non-leptonic decays with the three-point QCD sum rules to make more reasonable assignment.

The predicted mass $M_{Z}=4.01\pm 0.09\,\rm{GeV}$ for the ground state tetraquark state $[uc]_{\tilde{A}}[\bar{d}\bar{c}]_{A}-[uc]_{A}[\bar{d}\bar{c}]_{\tilde{A}}$ and $M_{Z}=4.00\pm 0.09\,\rm{GeV}$ for the ground state tetraquark state $[uc]_{A}[\bar{d}\bar{c}]_{A}$ both exhibit very good agreement with the experimental values $M_{Z(4020/4025)}=(4026.3\pm 2.6\pm 3.7)\,\rm{MeV}$ [7] and $(4022.9\pm 0.8\pm 2.7)\,\rm{MeV}$ [8] from the BESIII collaboration. There are two axialvector tetraquark state candidates with the quantum numbers $J^{PC}=1^{+-}$ for the $Z_{c}(4020)$. Again the two-body strong decays should be studied to make the assignment more reasonably.

The predicted mass $M_{Z}=4.60\pm 0.09\,\rm{GeV}$ for the first radially excited tetraquark state $[uc]_{\tilde{A}}[\bar{d}\bar{c}]_{A}-[uc]_{A}[\bar{d}\bar{c}]_{\tilde{A}}$ and $M_{Z}=4.58\pm 0.09\,\rm{GeV}$ for the first radially excited tetraquark state $[uc]_{A}[\bar{d}\bar{c}]_{A}$ both exhibit very good agreement with the experimental value $M_{Z(4600)}=4600\,\rm{MeV}$ from the LHCb collaboration [1]. On the other hand, the predicted mass $M_{Z}=4.66\pm 0.10\,\rm{GeV}$ for the ground state tetraquark state $[uc]_{\tilde{V}}[\bar{d}\bar{c}]_{V}+[uc]_{V}[\bar{d}\bar{c}]_{\tilde{V}}$ is also compatible with the experimental data $M_{Z(4600)}=4600\,\rm{MeV}$ from the LHCb collaboration [1]. Furthermore, the decay $Z_{c}(4600)\to J/\psi\pi$ can take place more easily for the ground state tetraquark state, which is in very good agreement with the observation of the $Z_{c}(4600)$ in the $J/\psi\pi$ invariant mass spectrum [1]. In summary, there are three axialvector tetraquark state candidates with $J^{PC}=1^{+-}$ for the $Z_{c}(4600)$, we still need more theoretical and experimental works to identify the $Z_{c}(4600)$ unambiguously.

In Ref.[2], we identify the $Z_{c}(4600)$ as the $[dc]_{P}[\bar{u}\bar{c}]_{A}-[dc]_{A}[\bar{u}\bar{c}]_{P}$ type vector tetraquark state tentatively according to the predicted mass $M_{Z}=(4.59\pm 0.08)\,\rm{GeV}$ from the QCD sum rules [46], and explore its non-leptonic decays $Z_{c}(4600)\to J/\psi\pi$, $\eta_{c}\rho$, $J/\psi a_{0}$, $\chi_{c0}\rho$, $D^{*}\bar{D}^{*}$, $D\bar{D}$, $D^{*}\bar{D}$ and $D\bar{D}^{*}$ with the QCD sum rules by matching the hadronic representation with the QCD representation with solid quark-hadron duality. The large partial decay width $\Gamma(Z_{c}^{-}(4600)\to J/\psi\pi^{-})=41.4^{+20.5}_{-14.9}\,\rm{MeV}$ exhibits very good agreement with the observation of the $Z_{c}(4600)$ in the $J/\psi\pi^{-}$ invariant mass spectrum.

In Table 3, we also present the diquark spin $S_{uc}$, antidiquark spin $S_{\bar{d}\bar{c}}$ and total spin $S$ of the hidden-charm tetraquark states. We examine the Table and find that the $[uc]_{\tilde{A}}[\bar{d}\bar{c}]_{A}-[uc]_{A}[\bar{d}\bar{c}]_{\tilde{A}}$ and $[uc]_{A}[\bar{d}\bar{c}]_{A}$ tetraquark states (which have the tetraquark structures $|1^{+},1;1\rangle-|1,1^{+};1\rangle$ and $|1,1;1\rangle$, respectively) have slightly larger masses than the $[uc]_{S}[\bar{d}\bar{c}]_{A}-[uc]_{A}[\bar{d}\bar{c}]_{S}$ tetraquark state (which has the structure $|0,1;1\rangle-|1,0;1\rangle$). It is reasonable, as the most favored diquark configurations or quark-quark correlations from the attractive interaction in the color-antitriplet (color-triplet) channel induced by one-gluon exchange are the scalar diquark (antidiquark) states. In previous works, the scalar, pseudoscalar, vector and axialvector diquarks states have been investigated with the QCD sum rules, the scalar and axialvector heavy diquark states in the color-antitriplet have almost degenerate masses or have almost the same typical quark-quark correlation lengths, the mass gaps between the scalar and axialvector heavy diquark states are very small or tiny [33]. Furthermore, it agrees with the predictions of the simple constituent diquark-antidiquark model [27].

The vector (or P-wave) diquark states $[uc]_{V}$ and $[uc]_{\tilde{V}}$ are expected to have larger masses than the axialvector (or S-wave) diquark states $[uc]_{A}$ and $[uc]_{\tilde{A}}$, as there exists a relative P-wave between the light quark and heavy quark. In the case of the traditional $c\bar{u}$ charmed mesons, the energy exciting a P-wave costs about $458\,\rm{MeV}$ from the Particle Data Group [9],

[TABLE]

If the energy exciting a P-wave in the $qc$ diquark systems also costs about $458\,\rm{MeV}$, the $[uc]_{\tilde{V}}[\bar{d}\bar{c}]_{V}+[uc]_{V}[\bar{d}\bar{c}]_{\tilde{V}}$ tetraquark state has the largest ground state mass, which is even larger than the masses of the first radially excited states of the hidden-charm tetraquark states $[uc]_{\tilde{A}}[\bar{d}\bar{c}]_{A}-[uc]_{A}[\bar{d}\bar{c}]_{\tilde{A}}$ and $[uc]_{A}[\bar{d}\bar{c}]_{A}$ with the quantum numbers $J^{PC}=1^{+-}$, as exciting two P-waves costs about $0.9\,\rm{GeV}$, which is larger than the energy gap $0.6\,\rm{GeV}$ between the ground state hidden-charm tetraquark state and the first radial excitation of the hidden-charm tetraquark states.

4 Conclusion

In this paper, we investigate the ground states and the first radially excited states of the $[uc]_{S}[\bar{d}\bar{c}]_{A}-[uc]_{A}[\bar{d}\bar{c}]_{S}$ type, $[uc]_{\tilde{A}}[\bar{d}\bar{c}]_{A}-[uc]_{A}[\bar{d}\bar{c}]_{\tilde{A}}$ type and $[uc]_{A}[\bar{d}\bar{c}]_{A}$ type tetraquark states and the ground state $[uc]_{\tilde{V}}[\bar{d}\bar{c}]_{V}+[uc]_{V}[\bar{d}\bar{c}]_{\tilde{V}}$ type tetraquark state with the quantum numbers $J^{PC}=1^{+-}$ via the QCD sum rules in a systematic way. The predicted tetraquark masses are in favor of assigning the $Z_{c}(3900)$ and $Z_{c}(4430)$ as the ground state and the first radially excited state of the $[uc]_{S}[\bar{d}\bar{c}]_{A}-[uc]_{A}[\bar{d}\bar{c}]_{S}$ type axialvector tetraquark states respectively; assigning the $Z_{c}(4020)$ as the ground state $[uc]_{\tilde{A}}[\bar{d}\bar{c}]_{A}-[uc]_{A}[\bar{d}\bar{c}]_{\tilde{A}}$ type axialvector tetraquark state or $[uc]_{A}[\bar{d}\bar{c}]_{A}$ type axialvector tetraquark state; assigning the $Z_{c}(4600)$ as the first radially excited $[uc]_{\tilde{A}}[\bar{d}\bar{c}]_{A}-[uc]_{A}[\bar{d}\bar{c}]_{\tilde{A}}$ type axialvector tetraquark state or $[uc]_{A}[\bar{d}\bar{c}]_{A}$ type axialvector tetraquark state, or the ground state $[uc]_{\tilde{V}}[\bar{d}\bar{c}]_{V}+[uc]_{V}[\bar{d}\bar{c}]_{\tilde{V}}$ type axialvector tetraquark state. We still need more experimental and theoretical works to identify the $Z_{c}(4600)$ unambiguously.

Acknowledgements

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