Exploring a $\Sigma_{c}\bar{D}$ state: with focus on $P_{c}(4312)^{+}$

TL;DR

This paper investigates the nature of the $P_{c}(4312)^{+}$ pentaquark candidate as a $ ext{Sigma}_c ar{D}$ molecular state using QCD sum rules, providing a mass estimate consistent with experimental observations.

Contribution

The study applies QCD sum rules to interpret $P_{c}(4312)^{+}$ as a $ ext{Sigma}_c ar{D}$ state, offering a theoretical mass range aligning with experimental data.

Findings

Mass range 4.07–4.97 GeV supports the $ ext{Sigma}_c ar{D}$ interpretation.

Results are consistent with the observed $P_{c}(4312)^{+}$ mass.

QCD sum rules effectively describe the pentaquark candidate.

Abstract

Stimulated by the new discovery of $P_{c}(4312)^{+}$ by LHCb Collaboration, we endeavor to perform the study of $P_{c}(4312)^{+}$ as a $\Sigma_{c}\bar{D}$ state in the framework of QCD sum rules. Taking into account the results from two sum rules, a conservative mass range $4.07\sim4.97~\mbox{GeV}$ is presented for the $\Sigma_{c}\bar{D}$ hadronic system, which agrees with the experimental data of $P_{c}(4312)^{+}$ and could support its interpretation as a $\Sigma_{c}\bar{D}$ state.

Equations

$$\Pi (q^2)=i\int d^{4}x\text{mbox}{e}^{iq.x}⟨0|T[j(x)\overline{j}(0)]|0⟩.$$

$$j={ϵ}_{abe}(q_a^{T}C{\gamma }_{\mu }{q}_b){\gamma }^{\mu }{\gamma }_5{c}_e{\bar{c}}_{f}i{\gamma }_5{q}_f.$$

$$\Pi (q^2)={\Pi }_1(q^2)+to0.0pt/\text{hss}q{\Pi }_2(q^2).$$

$$\Pi (q^2)={\lambda }_H^2\frac{M_H+to0.0pt/\text{hss}q}{M_H^2-q^2}+\frac{1}{\pi }{\int }_{s_0}^{\infty }ds\frac{\text{mbox}Im{\Pi }_1^{\text{mbox}phen}+to0.0pt/\text{hss}q\text{mbox}Im{\Pi }_2^{\text{mbox}phen}}{s-q^2}+...$$

$$\Pi (q^2)={\int }_{4m_c^2}^{\infty }ds\frac{{\rho }_1}{s-q^2}+to0.0pt/\text{hss}q{\int }_{4m_c^2}^{\infty }ds\frac{{\rho }_2}{s-q^2}+...$$

$${\lambda }_H^2{M}_H{e}^{-M_H^2/M^2}$$

$${\lambda }_H^2{e}^{-M_H^2/M^2}$$

$$M_H^2$$

$${\rho }_i={\rho }_i^{\text{mbox}pert}+{\rho }_i^{⟨\bar{q}q⟩}+{\rho }_i^{⟨g^2{G}^2⟩}+{\rho }_i^{⟨g\bar{q}\sigma \cdot Gq⟩}+{\rho }_i^{⟨\bar{q}q{⟩}^2}+{\rho }_i^{⟨g^3{G}^3⟩}+{\rho }_i^{⟨\bar{q}q⟩⟨g^2{G}^2⟩}+{\rho }_i^{⟨\bar{q}q⟩⟨g\bar{q}\sigma \cdot Gq⟩}$$

$${\rho }_1^{\text{mbox}pert}$$

$${\rho }_1^{⟨\bar{q}q⟩}$$

$${\rho }_1^{⟨g^2{G}^2⟩}$$

$${\rho }_1^{⟨g\bar{q}\sigma \cdot Gq⟩}$$

$${\rho }_1^{⟨\bar{q}q{⟩}^2}$$

$${\rho }_1^{⟨g^3{G}^3⟩}$$

$${\rho }_1^{⟨\bar{q}q⟩⟨g^2{G}^2⟩}$$

$${\rho }_1^{⟨\bar{q}q⟩⟨g\bar{q}\sigma \cdot Gq⟩}$$

$${\rho }_2^{\text{mbox}pert}$$

$${\rho }_2^{⟨\bar{q}q⟩}$$

$${\rho }_2^{⟨g^2{G}^2⟩}$$

$${\rho }_2^{⟨g\bar{q}\sigma \cdot Gq⟩}$$

$${\rho }_2^{⟨\bar{q}q{⟩}^2}$$

$${\rho }_2^{⟨g^3{G}^3⟩}$$

$${\rho }_2^{⟨\bar{q}q⟩⟨g^2{G}^2⟩}$$

$${\rho }_2^{⟨\bar{q}q⟩⟨g\bar{q}\sigma \cdot Gq⟩}$$

Full text

Keywords: QCD sum rules; $P_{c}(4312)^{+}$

I Introduction

Very recently, LHCb Collaboration reported the discovery of a narrow state $P_{c}(4312)^{+}$ with a statistical significance of $7.3\sigma$ in a data sample of $\Lambda_{b}^{0}\rightarrow J/\psi pK^{-}$ decays LHCb . Moreover, $P_{c}(4450)^{+}$ formerly announced by LHCb is confirmed and observed to consist of two narrow overlapping peaks, $P_{c}(4440)^{+}$ and $P_{c}(4457)^{+}$. Soon after the LHCb’s new observation, many works th-1 ; th-2 ; th-3 ; th-4 ; th-5 ; th-6 ; th-7 ; th-8 ; th-9 ; th-10 ; th-11 ; th-12 ; th-13 ; th-14 ; th-15 ; th-16 have been promptly triggered. Among these new experimental results, the most exciting point should attribute to the freshly discovered $P_{c}(4312)^{+}$. After all, there already have existed plenty of researches on $P_{c}(4450)^{+}$ others (one can also see a recent review e.g. Penta-overview ). Besides, $P_{c}(4312)^{+}$ is narrow and below the $\Sigma_{c}^{+}\bar{D}^{0}$ threshold within a plausible hadron-hadron binding energy, hence it provides the strongest experimental evidence to date for the existence of a $\Sigma_{c}\bar{D}$ bound state LHCb . Meanwhile, some different opinion has also appeared in Ref. th-15 , in which the authors could find evidence for the attractive effect of the $\Sigma_{c}^{+}\bar{D}^{0}$ channel, however not strong enough to form a bound state and they infer that the $P_{c}(4312)^{+}$ peak is more likely to be a virtual (unbound) state instead. Whether or no, to realize the nature of $P_{c}(4312)^{+}$, it certainly requires more theoretical scrutiny.

In this work, we focus all our attention on the newly discovered $P_{c}(4312)^{+}$ and would investigate the possibility of $P_{c}(4312)^{+}$ being a $\Sigma_{c}\bar{D}$ state, if the $\Sigma_{c}\bar{D}$ state does exist. While studying a baryon-meson state, one inevitably has to confront and treat nonperturbative QCD problem. As one reliable method for evaluating nonperturbative effects, the QCD sum rule svzsum is an analytic formalism firmly established on QCD theory and has been successfully applied to different hadronic systems overview1 ; overview2 ; overview3 ; reinders ; overview4 . As a matter of fact, there have appeared some related works on these $P_{c}$ hadrons basing on baryon-meson configuration QCD sum rules th-1 ; Peta-SR1 ; Peta-SR2 ; Peta-SR3 ; Peta-SR4 . In QCD sum rule analysis, it is of great importance to carefully inspect both the operator product expansion (OPE) convergence and the pole dominance in order to ensure the extracted result authentic. In practice, one could note that some condensate may play an important role in some multiquark cases Zs0 ; Zs1 ; Zs2 ; Zs , which causes that it is of difficulty to find conventional work windows. Specially for the four-quark condensate, a general factorization $\langle\bar{q}q\bar{q}q\rangle=\varrho\langle\bar{q}q\rangle^{2}$ has been hotly discussed A ; B , where $\varrho$ is a constant, which may be equal to $1$, to $2$, or be smaller than $1$. Moreover, the factorization parameter $\varrho$ could be about $3\sim 4$ C . Compromisingly, the parameter $\varrho$ is taken as $2$ in this work.

The rest paper is organized as follows. In Sec. II, $P_{c}(4312)^{+}$ is studied as a $\Sigma_{c}\bar{D}$ state through the QCD sum rule approach. Numerical analysis and discussions are given in Sec. III. The last part is a brief summary.

II QCD sum rule study of $P_{c}(4312)^{+}$ as a $\Sigma_{c}\bar{D}$ state

Mass sum rules for a $\Sigma_{c}\bar{D}$ state can be derived from the two-point correlator

[TABLE]

To represent the $\Sigma_{c}\bar{D}$ state, one can construct its interpolating current $j$ from baryon-meson type of fields adopting currents for the heavy baryon baryon-current and for the heavy meson reinders . Concretely, the current can be written as

[TABLE]

Here $q$ could be the light $u$ or $d$ quark, $c$ denotes the heavy charm quark, $T$ means matrix transposition, $C$ is the charge conjugation matrix, and the subscript $a$, $b$, $e$, and $f$ are color indices.

Lorentz covariance implies that the two-point correlator (1) has the general form

[TABLE]

In phenomenology, it can be expressed as

[TABLE]

where $M_{H}$ is the hadron’s mass, and $\lambda_{H}$ denotes the coupling of the current to the hadron $\langle 0|j|H\rangle=\lambda_{H}u(p,s)$. In the OPE side, one can write the correlator as

[TABLE]

where spectral densities are $\rho_{i}=\frac{1}{\pi}\mbox{Im}\Pi_{i}^{\mbox{OPE}}$, with $i=1,2$. After equating the two expressions, applying quark-hadron duality, and making a Borel transform, the sum rules are

[TABLE]

and

[TABLE]

where $M^{2}$ indicates the Borel parameter. Taking the derivative of Eq. (5) or (6) with respect to $1/M^{2}$ and dividing the equation itself, one can obtain mass sum rules

[TABLE]

and

[TABLE]

In the deriving of spectral densities, one can utilize the similar techniques as Refs. e.g. overview4 ; Zhang . The heavy-quark propagator in momentum-space reinders can be used to keep the heavy-quark mass finite, and the correlator’s light-quark part can be obtained in the coordinate space, which is then Fourier-transformed to the $D$ dimension momentum space. The resulting light-quark part is combined with the heavy-quark part before it is dimensionally regularized at $D=4$. As follows, we concretely present spectral densities $\rho_{i}$ deduced from $\Pi_{i}(q^{2})$ and put them forward to further numerical analysis, with

[TABLE]

up to dimension $8$. In detail,

[TABLE]

and

[TABLE]

in which the general $\langle\bar{q}q\bar{q}q\rangle=\varrho\langle\bar{q}q\rangle^{2}$ factorization has been used. The integration limits are \alpha_{min}=\Big{(}1-\sqrt{1-4m_{c}^{2}/s}\Big{)}/2, \alpha_{max}=\Big{(}1+\sqrt{1-4m_{c}^{2}/s}\Big{)}/2, and $\beta_{min}=\alpha m_{c}^{2}/(s\alpha-m_{c}^{2})$. Those condensates higher than dimension $8$ are not involved here, as one could expect that kind of high dimension contributions may not radically influence the OPE’s character Zhang1 ; Zhang2 .

III numerical analysis and discussions

In this part, we firstly perform the numerical analysis of sum rule (8) to extract the value of $M_{H}$, and take $m_{c}$ as the running charm quark mass $1.275_{-0.035}^{+0.025}~{}\mbox{GeV}$ PDG along with other input parameters as $\langle\bar{q}q\rangle=-(0.24\pm 0.01)^{3}~{}\mbox{GeV}^{3}$, $\langle g\bar{q}\sigma\cdot Gq\rangle=m_{0}^{2}~{}\langle\bar{q}q\rangle$, $m_{0}^{2}=0.8\pm 0.1~{}\mbox{GeV}^{2}$, $\langle g^{2}G^{2}\rangle=0.88\pm 0.25~{}\mbox{GeV}^{4}$, and $\langle g^{3}G^{3}\rangle=0.58\pm 0.18~{}\mbox{GeV}^{6}$ svzsum ; overview2 . Steering a middle course, the factorization parameter $\varrho$ is set to be $2$. According to a standard procedure, both the OPE convergence and the pole dominance should be considered to find appropriate work windows for the threshold $\sqrt{s_{0}}$ and the Borel parameter: the lower bound of $M^{2}$ is gained by analyzing the OPE convergence, and the upper one is obtained by viewing that the pole contribution should be larger than QCD continuum contribution. Besides, $\sqrt{s_{0}}$ characterizes the beginning of continuum states and should not be taken at will. It is correlated to the next excited state energy and empirically $400\sim 600~{}\mbox{MeV}$ above the eventually achieved value $M_{H}$.

In FIG. 1, the relative contributions of various OPE in sum rule (6) are compared as a function of $M^{2}$ for the $\Sigma_{c}\bar{D}$ state. Visually, there four main condensate contributions could play an important role on the OPE side, i.e. the two-quark condensate $\langle\bar{q}q\rangle$, the mixed condensate $\langle g\bar{q}\sigma\cdot Gq\rangle$, the four-quark condensate $\langle\bar{q}q\rangle^{2}$, and the $\langle\bar{q}q\rangle\langle g\bar{q}\sigma\cdot Gq\rangle$ condensate. The direct consequence is that it is not easy to find the standard Borel window, in which the low dimension condensate contribution should be bigger than the high dimension one. To say the least, these four main condensates could cancel each other out to some extent. In this way, the perturbative term still plays an important role on the OPE side and the OPE’s convergence could be under control at the relatively low value of $M^{2}$. Thus, the lower bound of $M^{2}$ is taken as $2.0~{}\mbox{GeV}^{2}$ for the sum rule (6).

Phenomenologically, a comparison between pole contribution and continuum contribution of sum rule (6) for $\sqrt{s_{0}}=4.8~{}\mbox{GeV}$ is shown in FIG. 2, which manifests that the relative pole contribution is about $50\%$ at $M^{2}=2.7~{}\mbox{GeV}^{2}$ and decreases with $M^{2}$. In a similar way, the upper bounds of Borel parameters are $M^{2}=2.6~{}\mbox{GeV}^{2}$ for $\sqrt{s_{0}}=4.7~{}\mbox{GeV}$ and $M^{2}=2.9~{}\mbox{GeV}^{2}$ for $\sqrt{s_{0}}=4.9~{}\mbox{GeV}$. Thus, Borel windows are taken as $2.0\sim 2.6~{}\mbox{GeV}^{2}$ for $\sqrt{s_{0}}=4.7~{}\mbox{GeV}$, $2.0\sim 2.7~{}\mbox{GeV}^{2}$ for $\sqrt{s_{0}}=4.8~{}\mbox{GeV}$, and $2.0\sim 2.9~{}\mbox{GeV}^{2}$ for $\sqrt{s_{0}}=4.9~{}\mbox{GeV}$. The mass $M_{H}$ of $\Sigma_{c}\bar{D}$ is shown in FIG. 3 as a function of $M^{2}$ from sum rule (8). In the chosen work windows, $M_{H}$ is calculated to be $4.35\pm 0.07~{}\mbox{GeV}$. Furthermore, in view of the uncertainty due to the variation of quark masses and condensates, we have $4.35\pm 0.07^{+0.55}_{-0.21}~{}\mbox{GeV}$ (the first error is resulted from the variation of $\sqrt{s_{0}}$ and $M^{2}$, and the second error reflects the uncertainty rooting in the variation of QCD parameters) or briefly $4.35^{+0.62}_{-0.28}~{}\mbox{GeV}$ for $\Sigma_{c}\bar{D}$.

Furthermore, one could put forward the numerical analysis of sum rule (7) analogously. In FIG. 4, the relative contributions of various OPE in sum rule (5) are shown as a function of $M^{2}$ for $\sqrt{s_{0}}=4.8~{}\mbox{GeV}$. Similarly, four main condensates (i.e. $\langle\bar{q}q\rangle$, $\langle g\bar{q}\sigma\cdot Gq\rangle$, $\langle\bar{q}q\rangle^{2}$, and $\langle\bar{q}q\rangle\langle g\bar{q}\sigma\cdot Gq\rangle$) could cancel each other out to some extent. For the sum rule (5), the lower bound of $M^{2}$ is taken as $2.2~{}\mbox{GeV}^{2}$ at which the OPE’s convergence could still be controllable. In FIG. 5, a comparison between pole and continuum contribution of sum rule (5) is shown for $\sqrt{s_{0}}=4.8~{}\mbox{GeV}$, which indicates that the relative pole contribution is about $50\%$ at $M^{2}=2.9~{}\mbox{GeV}^{2}$ and decreases with $M^{2}$. Thereby, the ranges of $M^{2}$ are fixed as $2.2\sim 2.9~{}\mbox{GeV}^{2}$ for $\sqrt{s_{0}}=4.7~{}\mbox{GeV}$, $2.2\sim 3.1~{}\mbox{GeV}^{2}$ for $\sqrt{s_{0}}=4.8~{}\mbox{GeV}$, and $2.2\sim 3.2~{}\mbox{GeV}^{2}$ for $\sqrt{s_{0}}=4.9~{}\mbox{GeV}$. The mass of $\Sigma_{c}\bar{D}$ state is shown in FIG. 6 as a function of $M^{2}$ from sum rule (7). In the chosen work windows, $M_{H}$ is calculated to be $4.38\pm 0.09~{}\mbox{GeV}$. In view of the uncertainty due to the variation of quark masses and condensates, we have $4.38\pm 0.09^{+0.13}_{-0.07}~{}\mbox{GeV}$ (the first error is resulted from the variation of $\sqrt{s_{0}}$ and $M^{2}$, and the second error reflects the uncertainty rooting in the variation of QCD parameters) or briefly $4.38^{+0.22}_{-0.16}~{}\mbox{GeV}$ for $\Sigma_{c}\bar{D}$.

In the end, combining the eventual results from both (7) and (8), one could arrive at a conservative mass range $4.07\sim 4.97~{}\mbox{GeV}$ for the $\Sigma_{c}\bar{D}$ state, which is consistent with the data of $P_{c}(4312)^{+}$ and could support its explanation as a $\Sigma_{c}\bar{D}$ state.

IV Summary

Motivated by LHCb’s new discovery of $P_{c}(4312)^{+}$, we study that whether $P_{c}(4312)^{+}$ could be a $\Sigma_{c}\bar{D}$ state in QCD sum rules. In order to insure the quality of sum rule analysis, contributions of condensates up to dimension $8$ have been computed to test the OPE convergence. We find that some condensates, i.e. the two-quark condensate, the mixed condensate, the four-quark condensate, and the $\langle\bar{q}q\rangle\langle g\bar{q}\sigma\cdot Gq\rangle$ condensate are of importance to the OPE side. Not bad, those main condensates could cancel each other out to some extent, which brings that the OPE convergence is still controllable. By combining those results from two sum rules, we finally obtain that a conservative mass range for $\Sigma_{c}\bar{D}$ is $4.07\sim 4.97~{}\mbox{GeV}$, which is in agreement with the experimental value of $P_{c}(4312)^{+}$. This result supports that $P_{c}(4312)^{+}$ could be explained as a $\Sigma_{c}\bar{D}$ state.

In the future, one can expect that further experimental observations may shed more light on the nature of $P_{c}(4312)^{+}$ and the inner structure of $P_{c}(4312)^{+}$ could be further revealed by continual efforts in both experiment and theory.

Acknowledgements.

This work was supported by the National Natural Science Foundation of China under Contract Nos. 11475258, 11105223, and 11675263, and by the project for excellent youth talents in NUDT.