Exploring a ΣcD¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma _{c}\bar{D}$$\end{document} state: with focus on Pc(4312)+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{c}(4312)^{+}$$\end{document}

Stimulated by the new discovery of Pc(4312)+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{c}(4312)^{+}$$\end{document} by LHCb Collaboration, we endeavor to perform the study of Pc(4312)+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{c}(4312)^{+}$$\end{document} as a ΣcD¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma _{c}\bar{D}$$\end{document} state in the framework of QCD sum rules. Taking into account the results from two sum rules, a conservative mass range 4.07∼4.97GeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4.07{\sim }4.97~\text{ GeV }$$\end{document} is presented for the ΣcD¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma _{c}\bar{D}$$\end{document} hadronic system, which agrees with the experimental data of Pc(4312)+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{c}(4312)^{+}$$\end{document} and could support its interpretation as a ΣcD¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma _{c}\bar{D}$$\end{document} state.

In this work, we focus all our attention on the newly discovered P c (4312) + and would investigate the possibility of a e-mail: jrzhang@nudt.edu.cn P c (4312) + being a cD state, if the cD 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 [36][37][38] is an analytic formalism firmly established on QCD theory and has been successfully applied to different hadronic systems [39][40][41][42][43]. As a matter of fact, there have appeared some related works on these P c hadrons basing on baryon-meson configuration QCD sum rules [2,[44][45][46][47]. 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 [48][49][50][51][52], which causes that it is of difficulty to find conventional work windows. Specially for the four-quark condensate, a general factorization qqqq = qq 2 has been hotly discussed [53][54][55][56][57][58][59][60][61][62], where is a constant, which may be equal to 1, to 2, or be smaller than 1. Moreover, the factorization parameter could be about 3∼4 [63][64][65]. Compromisingly, the parameter is taken as 2 in this work.
The rest paper is organized as follows. In Sect. 2, P c (4312) + is studied as a cD state through the QCD sum rule approach. Numerical analysis and discussions are given in Sect. 3. The last part is a brief summary.

QCD sum rule study of P c (4312) + as a cD state
Mass sum rules for a cD state can be derived from the two-point correlator To represent the cD state, one can construct its interpolating current j from baryon-meson type of fields adopting currents for the heavy baryon [66,67] and for the heavy meson [42].
Concretely, the current can be written as j = abe (q T a Cγ μ q b )γ μ γ 5 c ec f iγ 5 q f . 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 In phenomenology, it can be expressed as where M H is the hadron's mass, and λ H denotes the coupling of the current to the hadron 0| j|H = λ H u( p, s). In the OPE side, one can write the correlator as where spectral densities are ρ i = 1 π Im OPE i , with i = 1, 2. After equating the two expressions, applying quark-hadron duality, and making a Borel transform, the sum rules are and where M 2 indicates the Borel parameter. Taking the derivative of Eqs. (5) or (6) with respect to 1/M 2 and dividing the equation itself, one can obtain mass sum rules and In the deriving of spectral densities, one can utilize the similar techniques as Refs. e.g. [43,[68][69][70]. The heavy-quark propagator in momentum-space [42] 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 heavyquark part before it is dimensionally regularized at D = 4. As follows, we concretely present spectral densities ρ i deduced from i (q 2 ) and put them forward to further numerical analysis, with in which the general qqqq = qq 2 factorization has been used. The integration limits are α min = 1 − . 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 [71][72][73][74].

Numerical analysis and discussions
In this part, we firstly perform the numerical analysis of sum rule (8) [36][37][38]40]. Steering a middle course, the factorization parameter 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 √ 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, √ 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∼600 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 cD state. Visually, there four main condensate contributions could play an important role on the OPE side, i.e. the two-quark condensate qq , the mixed condensate gqσ · Gq , the four-quark condensate qq 2 , and the qq gqσ · Gq 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 GeV 2 for the sum rule (6).
Phenomenologically, a comparison between pole contribution and continuum contribution of sum rule (6) for √ s 0 = 4.8 GeV is shown in Fig. 2, which manifests that the relative pole contribution is about 50% at M 2 = 2.7 GeV 2 −0.28 GeV for cD . 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 √ s 0 = 4.8 GeV. Similarly, four main condensates (i.e. qq , gqσ ·Gq , qq 2 , and qq gqσ ·Gq ) could cancel each other out to some extent. For the sum rule (5), the lower bound of M 2 is taken as 2.2 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 √ s 0 = 4.8 GeV, which indicates that the relative pole contribution is about 50% at M 2 = 2.9 GeV 2 and decreases with M 2 . Thereby, the ranges of M 2 are fixed as 2.2∼2.9 GeV 2 for √ s 0 = 4.7 GeV, 2.2∼3.1 GeV 2 for √ s 0 = 4.8 GeV, and 2.2∼3.2 GeV 2 for √ s 0 = 4.9 GeV. The mass of cD 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 ± 0.09 GeV. In view of the uncertainty due to the variation of quark masses and condensates, we  (7) and (8), one could arrive at a conservative mass range 4.07∼4.97 GeV for the cD state, which is consistent with the data of P c (4312) + and could support its explanation as a cD state.

Summary
Motivated by LHCb's new discovery of P c (4312) + , we study that whether P c (4312) + could be a cD 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 qq gqσ · Gq 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 con-trollable. By combining those results from two sum rules, we finally obtain that a conservative mass range for cD is 4.07∼4.97 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 cD 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.
Data Availability Statement This manuscript has no associated data or the data will not be deposited. [Authors' comment: All data generated during this study are included in this published article.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .