Anatomy of the newly observed hidden-charm pentaquark states: $P_c(4312)$, $P_c(4440)$ and $P_c(4457)$

TL;DR

This paper investigates the internal structure and composition of the newly observed hidden-charm pentaquark states $P_c(4312)$, $P_c(4440)$, and $P_c(4457)$ using effective-range expansion and resonance compositeness relations, providing insights into their molecular components.

Contribution

It introduces a detailed analysis of the pentaquark states' scattering properties and compositeness, revealing their likely molecular nature and channel contributions.

Findings

Calculated scattering lengths and effective ranges for $Σ_c\bar{D}$ and $Σ_c\bar{D}^*$ channels.

Estimated probabilities of $J/\psi p$ and $Σ_c\bar{D}$ components inside $P_c(4312)$.

Determined partial decay widths and compositeness coefficients for $P_c(4440)$ and $P_c(4457)$.

Abstract

We study the newly reported hidden-charm pentaquark candidates $P_c(4312)$, $P_c(4440)$ and $P_c(4457)$ from the LHCb Collaboration, in the framework of the effective-range expansion and resonance compositeness relations. The scattering lengths and effective ranges from the $S$-wave $\Sigma_c\bar{D}$ and $\Sigma_c\bar{D}^*$ scattering are calculated by using the experimental results of the masses and widths of the $P_c(4312)$, $P_c(4440)$ and $P_c(4457)$. Then we calculate the couplings between the $J/\psi p,\,\Sigma_c\bar{D}$ channels and the pentaquark candidate $P_c(4312)$, with which we further estimate the probabilities of finding the $J/\psi p$ and $\Sigma_c\bar{D}$ components inside $P_c(4312)$. The partial decay widths and compositeness coefficients are calculated for the $P_c(4440)$ and $P_c(4457)$ states by including the $J/\psi p$ and $\Sigma_c\bar{D}^*$ channels. Similar studies are also carried out for the three $P_c$ states by including the $\Lambda_c\bar{D}^{*}$ and $\Sigma_c\bar{D}^{(*)}$ channels.

Tables

Table 1. Table 1: The values of the scattering lengths and effective ranges of the S 𝑆 S -wave amplitudes for different channels. The uncertainties for a 𝑎 a and r 𝑟 r are determined by adding in quadrature the resulting ones from the systematic and statistical errors of the masses and widths of the P c subscript 𝑃 𝑐 P_{c} states. The errors of the different thresholds are negligible in comparison with the uncertainties of the masses and widths of the P c subscript 𝑃 𝑐 P_{c} states.

Resonance	Mass	Width	Threshold	$a$	$r$
	(MeV)	(MeV)	(MeV)	(fm)	(fm)
$P_c\left(4312\right)$	$4311.9\pm 6.8$	$9.8\pm 5.2$	${\Sigma }_c^{+}{\overline{D}}^0$ (4317.7)	$-2.9\pm 0.8$	$-1.7\pm 0.7$
			${\Sigma }_c^{++}{D}^{-}$ (4323.6)	$-2.4\pm 0.6$	$-1.2\pm 0.3$
$P_c\left(4440\right)$	$4440.3\pm 4.9$	$20.6\pm 11.2$	${\Sigma }_c^{+}{\overline{D}}^{\ast 0}$ (4459.8)	$-1.7\pm 0.2$	$-0.9\pm 0.1$
			${\Sigma }_c^{++}{D}^{\ast -}$ (4464.2)	$-1.6\pm 0.2$	$-0.8\pm 0.1$
$P_c\left(4457\right)$	$4457.3\pm 4.1$	$6.4\pm 6.0$	${\Sigma }_c^{+}{\overline{D}}^{\ast 0}$ (4459.8)	$-3.8\pm 1.6$	$-2.3\pm 1.3$
			${\Sigma }_c^{++}{D}^{\ast -}$ (4464.2)	$-3.0\pm 0.7$	$-1.6\pm 0.4$

Table 2. Table 2: Results obtained with X = X 1 + X 2 = 1 𝑋 subscript 𝑋 1 subscript 𝑋 2 1 X=X_{1}+X_{2}=1 . The J / ψ ​ p 𝐽 𝜓 𝑝 J/\psi p and Σ c ​ D ¯ ( ∗ ) subscript Σ 𝑐 superscript ¯ 𝐷 \Sigma_{c}\bar{D}^{(*)} channels, which are labeled as 1 1 1 and 2 2 2 respectively, are included.

Resonance	$\left|g_1\right|$	$\left|g_2\right|$	${\Gamma }_1$	${\Gamma }_2$	$X_1$	$X_2$
	(GeV)	(GeV)	(MeV)	(MeV)
$P_c\left(4312\right)$
$m_{{\Sigma }_c^{+}}+m_{{\overline{D}}^0}$	${2.1}_{-2.1}^{+0.8}$	${10.9}_{-2.9}^{+2.1}$	${6.5}_{-6.5}^{+4.9}$	${3.3}_{-3.3}^{+10.5}$	${0.006}_{-0.006}^{+0.005}$	${0.994}_{-0.005}^{+0.006}$
$m_{{\Sigma }_c^{++}}+m_{D^{-}}$	${2.5}_{-0.9}^{+0.6}$	${12.6}_{-2.6}^{+1.6}$	${8.5}_{-4.6}^{+4.7}$	${1.3}_{-1.3}^{+6.1}$	${0.008}_{-0.005}^{+0.005}$	${0.992}_{-0.005}^{+0.005}$
$P_c\left(4440\right)$
$m_{{\Sigma }_c^{+}}+m_{{\overline{D}}^{\ast 0}}$	${3.2}_{-0.9}^{+0.6}$	${14.9}_{-1.4}^{+1.2}$	${16.3}_{-7.4}^{+6.7}$	${4.3}_{-4.3}^{+9.2}$	${0.010}_{-0.004}^{+0.005}$	${0.990}_{-0.005}^{+0.004}$
$m_{{\Sigma }_c^{++}}+m_{D^{\ast -}}$	${3.3}_{-0.9}^{+0.6}$	${15.6}_{-1.1}^{+1.0}$	${17.7}_{-8.2}^{+6.9}$	${2.9}_{-2.9}^{+8.3}$	${0.011}_{-0.005}^{+0.005}$	${0.989}_{-0.005}^{+0.005}$
$P_c\left(4457\right)$
$m_{{\Sigma }_c^{+}}+m_{{\overline{D}}^{\ast 0}}$	${1.5}_{-1.0}^{+0.7}$	${9.5}_{-5.1}^{+2.2}$	${3.5}_{-3.5}^{+4.2}$	${2.9}_{-2.9}^{+9.5}$	${0.002}_{-0.002}^{+0.003}$	${0.998}_{-0.003}^{+0.002}$
$m_{{\Sigma }_c^{++}}+m_{D^{\ast -}}$	${1.8}_{-0.9}^{+0.6}$	${11.2}_{-2.5}^{+1.6}$	${5.4}_{-4.0}^{+4.2}$	${1.0}_{-1.0}^{+6.1}$	${0.003}_{-0.002}^{+0.003}$	${0.997}_{-0.003}^{+0.002}$

Table 3. Table 3: Results obtained for X = 0.8 𝑋 0.8 X=0.8 and X = 0.5 𝑋 0.5 X=0.5 by including the J / ψ ​ p 𝐽 𝜓 𝑝 J/\psi p (labeled as 1) and Σ c ​ D ¯ ( ∗ ) subscript Σ 𝑐 superscript ¯ 𝐷 \Sigma_{c}\bar{D}^{(*)} (labeled as 2) channels. The values in the table are calculated by using the masses Σ c + superscript subscript Σ 𝑐 \Sigma_{c}^{+} and D ¯ ( ∗ ) ​ 0 superscript ¯ 𝐷 0 \bar{D}^{(*)0} .

Resonance	$\left|g_1\right|$	$\left|g_2\right|$	${\Gamma }_1$	${\Gamma }_2$	$X_1$	$X_2$
	(GeV)	(GeV)	(MeV)	(MeV)
$P_c\left(4312\right)$
$X=0.8$	${2.3}_{-1.8}^{+0.7}$	${9.8}_{-2.5}^{+1.8}$	${7.1}_{-6.8}^{+5.0}$	${2.7}_{-2.7}^{+7.3}$	${0.007}_{-0.007}^{+0.005}$	${0.793}_{-0.005}^{+0.007}$
$X=0.5$	${2.4}_{-1.2}^{+0.7}$	${7.7}_{-2.0}^{+1.5}$	${8.1}_{-6.2}^{+5.1}$	${1.7}_{-1.7}^{+5.1}$	${0.008}_{-0.006}^{+0.005}$	${0.492}_{-0.005}^{+0.006}$
$P_c\left(4440\right)$
$X=0.8$	${3.2}_{-0.9}^{+0.7}$	${13.3}_{-1.3}^{+1.0}$	${17.2}_{-8.2}^{+7.6}$	${3.4}_{-3.4}^{+7.4}$	${0.011}_{-0.005}^{+0.005}$	${0.789}_{-0.005}^{+0.005}$
$X=0.5$	${3.4}_{-1.0}^{+0.7}$	${10.5}_{-1.0}^{+0.7}$	${18.5}_{-9.3}^{+9.0}$	${2.1}_{-2.1}^{+4.5}$	${0.012}_{-0.006}^{+0.006}$	${0.488}_{-0.006}^{+0.006}$
$P_c\left(4457\right)$
$X=0.8$	${1.6}_{-1.5}^{+0.7}$	${8.5}_{-4.5}^{+2.0}$	${4.1}_{-4.1}^{+4.6}$	${2.3}_{-2.3}^{+7.9}$	${0.002}_{-0.002}^{+0.003}$	${0.798}_{-0.003}^{+0.003}$
$X=0.5$	${1.7}_{-1.6}^{+0.8}$	${6.7}_{-3.3}^{+1.5}$	${5.0}_{-5.0}^{+5.1}$	${1.4}_{-1.4}^{+5.0}$	${0.003}_{-0.003}^{+0.003}$	${0.497}_{-0.003}^{+0.003}$

Table 4. Table 4: Results obtained when including the Λ c + ​ D ¯ ∗ 0 superscript subscript Λ 𝑐 superscript ¯ 𝐷 absent 0 \Lambda_{c}^{+}\bar{D}^{*0} (labeled as 1) and Σ c + ​ D ¯ ∗ 0 superscript subscript Σ 𝑐 superscript ¯ 𝐷 absent 0 \Sigma_{c}^{+}\bar{D}^{*0} (labeled as 2) channels for X = 1.0 𝑋 1.0 X=1.0 , 0.8 0.8 0.8 and 0.5 0.5 0.5 .

Resonance	$\left|g_1\right|$	$\left|g_2\right|$	${\Gamma }_1$	${\Gamma }_2$	$X_1$	$X_2$
	(GeV)	(GeV)	(MeV)	(MeV)
$P_c\left(4312\right)$
$X=1.0$	${4.0}_{-3.8}^{+2.0}$	${10.5}_{-2.5}^{+1.3}$	${6.8}_{-6.8}^{+5.4}$	${3.0}_{-3.0}^{+10.6}$	${0.09}_{-0.09}^{+0.16}$	${0.91}_{-0.16}^{+0.09}$
$X=0.8$	${4.2}_{-3.4}^{+2.0}$	${9.2}_{-2.0}^{+1.2}$	${7.5}_{-7.2}^{+5.5}$	${2.3}_{-2.3}^{+8.1}$	${0.10}_{-0.10}^{+0.16}$	${0.70}_{-0.16}^{+0.10}$
$X=0.5$	${4.5}_{-2.5}^{+2.0}$	${6.8}_{-1.2}^{+0.9}$	${8.5}_{-6.5}^{+5.7}$	${1.3}_{-1.3}^{+4.3}$	${0.11}_{-0.09}^{+0.17}$	${0.39}_{-0.17}^{+0.09}$
$P_c\left(4440\right)$
$X=1.0$	${3.8}_{-1.0}^{+0.7}$	${14.8}_{-1.3}^{+1.0}$	${16.4}_{-7.5}^{+6.8}$	${4.2}_{-4.2}^{+9.1}$	${0.03}_{-0.02}^{+0.01}$	${0.97}_{-0.01}^{+0.02}$
$X=0.8$	${3.9}_{-1.1}^{+0.8}$	${13.1}_{-1.1}^{+0.9}$	${17.3}_{-8.3}^{+7.7}$	${3.3}_{-3.3}^{+7.2}$	${0.03}_{-0.02}^{+0.01}$	${0.77}_{-0.01}^{+0.02}$
$X=0.5$	${4.0}_{-1.2}^{+1.0}$	${10.2}_{-0.8}^{+0.6}$	${18.6}_{-9.4}^{+9.2}$	${2.0}_{-2.0}^{+4.3}$	${0.03}_{-0.01}^{+0.02}$	${0.47}_{-0.02}^{+0.01}$
$P_c\left(4457\right)$
$X=1.0$	${1.7}_{-1.6}^{+0.9}$	${9.4}_{-5.0}^{+2.3}$	${3.5}_{-3.5}^{+3.7}$	${2.9}_{-2.9}^{+9.5}$	${0.005}_{-0.005}^{+0.007}$	${0.995}_{-0.007}^{+0.005}$
$X=0.8$	${1.9}_{-1.9}^{+0.8}$	${8.4}_{-4.4}^{+2.0}$	${4.1}_{-4.1}^{+4.6}$	${2.3}_{-2.3}^{+7.9}$	${0.006}_{-0.006}^{+0.008}$	${0.794}_{-0.008}^{+0.006}$
$X=0.5$	${2.0}_{-2.0}^{+0.9}$	${6.6}_{-3.2}^{+1.6}$	${5.0}_{-5.0}^{+5.1}$	${1.4}_{-1.4}^{+4.9}$	${0.008}_{-0.008}^{+0.008}$	${0.492}_{-0.008}^{+0.008}$

Equations

$$V(k)=-\frac{1}{a}+\frac{1}{2}r{k}^2+O(k^2),$$

$$k=\sqrt{2\mu (E-m_{\mathrm{th}})},$$

$$T(E)=\frac{1}{-\frac{1}{a}+\frac{1}{2}r{k}^2-ik},$$

$$\mathrm{Im}T(E{)}^{-1}=-k,(E>m_{\mathrm{th}}).$$

$$a\to -\frac{m_{\mathrm{th}}-M_{\mathrm{CDD}}}{g_{\mathrm{CDD}}},r\to -\frac{g_{\mathrm{CDD}}}{\mu (m_{\mathrm{th}}-M_{\mathrm{CDD}}{)}^2},$$

$$E_R=M_R-i{\Gamma }_R/2,$$

$$T_{II}(E)=\frac{1}{-\frac{1}{a}+\frac{1}{2}r{k}^2+ik}.$$

$$-\frac{1}{a}+\frac{1}{2}r{k}_R^2+i{k}_R=0,$$

$$k_R=\sqrt{\mu (E_R-m_{\mathrm{th}})}.$$

$$a$$

$$r$$

$$X_j=|g_j{|}^2\left|\frac{\partial G_j(s_R)}{\partial s}\right|,$$

$$G(s)=i\int \frac{{\mathrm{d}}^{4}q}{(2\pi {)}^4}\frac{1}{(q^2-m_1^2+iϵ)[(P-q{)}^2-m_2^2+iϵ]},s\equiv P^2.$$

$$G(s)$$

$$\sigma (s)=\sqrt{[s-(m_1+m_2{)}^2][s-(m_1-m_2{)}^2]}.$$

$${\Gamma }_1=|g_1{|}^2\frac{q(M_R^2)}{8\pi M_R^2},$$

$$q(M_R^2)=\frac{\sqrt{[M_R^2-(m_1+m_2{)}^2][M_R^2-(m_1-m_2{)}^2]}}{2M_R}.$$

$${\Gamma }_2=|g_2{|}^2{\int }_{m_{\mathrm{th}}}^{M_R+2{\Gamma }_R}dw\frac{q(w^2)}{16{\pi }^2{w}^2}\frac{{\Gamma }_R}{(M_R-w{)}^2+{\Gamma }_R^2/4}.$$

$$|g_1{|}^2\frac{q_1(M_R^2)}{8\pi M_R^2}+|g_2{|}^2{\int }_{m_{\mathrm{th}}}^{M_R+2{\Gamma }_R}dw\frac{q_2(w^2)}{16{\pi }^2{w}^2}\frac{{\Gamma }_R}{(M_R-w{)}^2+{\Gamma }_R^2/4}$$

$$|g_1{|}^2\left|\frac{\partial G_1^{II}(s_R)}{\partial s}\right|+|g_2{|}^2\left|\frac{\partial G_2(s_R)}{\partial s}\right|$$

Full text

1 Introduction

The first discovery of the hidden-charm pentaquark states $P_{c}(4380)$ and $P_{c}(4450)$ [1] has triggered a plethora of in-depth theoretical studies [2]. Very recently, the LHCb Collaboration has reported updated results on the pentaquark states based on the combinations of the Run 1 + Run 2 data [3]. The first notable finding from the updated measurements is that a new hidden-charm pentaquark state $P_{c}(4312)$ is observed with the mass $4311.9\pm 0.7^{+6.8}_{-0.6}$ MeV and the width $9.8\pm 2.7^{+3.7}_{-4.5}$ MeV. The second notable and intriguing observation is that the previous single state $P_{c}(4450)$ is superseded by two nearby states $P_{c}(4440)$ and $P_{c}(4457)$, with their masses $4440.3\pm 1.3^{+4.1}_{-4.7}$ MeV and $4457.3\pm 0.6^{+4.1}_{-1.7}$ MeV, respectively, and their widths $20.6\pm 4.9^{+8.7}_{-10.1}$ MeV and $6.4\pm 2.0^{+5.7}_{-1.9}$ MeV, respectively. The previous peak around the $P_{c}(4380)$ state now becomes less clear and its existence needs to be confirmed further by the experimental analysis. The new measurements have already attracted attention from many groups [4, 5, 6, 7, 8, 9, 10, 11, 12, 13].

All of the three new states are observed in the $J/\psi p$ invariant mass distributions from the $\Lambda_{b}\to J/\psi\,p\,K^{-}$ decay. One of the common features of the newly measured pentaquark states is that they all have small widths. Another important common feature is that they lie quite close to the thresholds of two underlying hadrons. In the following discussions we take a conservative way to estimate the experimental values of the masses and widths for the $P_{c}(4312),P_{c}(4440)$ and $P_{c}(4457)$ [3]. To be more specific, we take the larger values in magnitude of the upper or lower limits for the systematic uncertainties, and add them quadratically to the statistical ones as the total uncertainties. The resulting masses and widths are summarized in the second and third columns of Table 1. The differences between the mass of the $P_{c}(4312)$ and the $\Sigma_{c}^{+}\bar{D}^{0}$ and $\Sigma_{c}^{++}D^{-}$ thresholds are $5.8\pm 6.8$ MeV and $11.7\pm 6.8$ MeV, respectively. The mass of $P_{c}(4440)$ lie $19.5\pm 4.9$ MeV and $23.9\pm 4.9$ MeV below the $\Sigma_{c}^{+}\bar{D}^{*0}$ and $\Sigma_{c}^{++}D^{*-}$ thresholds, respectively. For the $P_{c}(4457)$, the differences between its mass and the $\Sigma_{c}^{+}\bar{D}^{*0}$ and $\Sigma_{c}^{++}D^{*-}$ thresholds are $2.5\pm 4.1$ MeV and $6.9\pm 4.1$ MeV, respectively. Taking into account the uncertainties of the experimental measurements of the $P_{c}(4312)$, we notice that its mass can be either below or above the $\Sigma_{c}^{+}\bar{D}^{0}$ threshold, but it is always below the $\Sigma_{c}^{++}D^{-}$ threshold. For the mass of the $P_{c}(4457)$ a similar situation occurs, so that, within the present experimental uncertainties [3], its mass can be also either below or above the $\Sigma_{c}^{+}\bar{D}^{*0}$ threshold, but it is always below the $\Sigma_{c}^{++}D^{*-}$ threshold. As a result one would expect that the isospin breaking effects could be visible [6]. In order to quantify the possible isospin breaking effects, we shall distinguish the elastic scattering with different thresholds involving $\Sigma_{c}^{+}$ or $\Sigma_{c}^{++}$ in a later study.

In this work, our key aim is to quantify the possibilities of the $P_{c}(4312)$ as the $S$-wave $\Sigma_{c}\bar{D}$, and the $P_{c}(4440)$ and $P_{c}(4457)$ as the $S$-wave $\Sigma_{c}\bar{D}^{*}$ molecular states. The effective range expansion (ERE) approach offers a reliable tool to analyze the dynamics around the threshold energy region. The combinations of the analyticity, unitarity and ERE have been demonstrated to be successful in analyzing the heavy-flavor exotic hadrons near thresholds [14, 15, 16, 17]. Another powerful tool that can help to reveal the inner structures of the hadrons is the Weinberg’s compositeness relation [18], which is extended to the resonance case in Refs. [19, 20, 21]. Other forms of generalization for other compositeness relation to address the resonances can be also found in Refs. [22, 23, 24, 25, 26]. In the current work we combine analyticity, unitarity, the ERE and the resonance compositeness relation to study the three newly measured pentaquark states.

2 Effective-range-expansion study of the pentaquark states

The ERE approach relies on the power series expansion of the $K$-matrix $V(k)$ at around threshold

[TABLE]

where $a$ is the scattering length, $r$ denotes the effective range and $k$ stands for the magnitude of three-momentum in the center of mass (CM) frame. For a two-particle system with masses $m_{1}$ and $m_{2}$, in the non-relativistic limit the three-momentum $k$ is related to the CM energy $E$ through

[TABLE]

with the threshold $m_{\rm th}=m_{1}+m_{2}$ and the reduced mass $\mu=\frac{m_{1}m_{2}}{m_{1}+m_{2}}$.

For the $\Sigma_{c}\bar{D}$ scattering near the $P_{c}(4312)$ and the $\Sigma_{c}\bar{D}^{*}$ scattering near the $P_{c}(4440)$ and $P_{c}(4457)$ energy regions, the magnitudes of the three-momenta of the two-particle systems can range from 0 to 250 MeV, after taking into account the experimental uncertainties of the masses of the $P_{c}$ states [3]. For the scattering of two heavy-flavor hadrons, it is plausible that the pion exchanges can be treated perturbatively [27, 28, 29, 30, 31]. For the heavier vector-resonance exchanges, their contributions can be effectively included via contact interactions, since their masses are clearly larger than the scale of the relevant three-momenta. Therefore we take the point of view from the pionless effective field theory, which only needs to include the local contact interactions [32].

Under these circumstances only the unitarity/right-hand cut enters and there is no crossed-channel dynamics. The elastic $S$-wave scattering amplitude around threshold that results from Eq. (1) (without the crossed-channel cuts) can be written as

[TABLE]

which satisfies the unitarity condition

[TABLE]

The formula $T(E)$ in Eq. (3) generally works well in the energy region near threshold even when resonances appear, except in the special situation that an underlying Castillejo-Dalitz-Dyson (CDD) pole sits on top of the threshold. In the latter case, one has to explicitly include the CDD pole in Eq. (3) and we refer to Ref. [17] for further details. It is difficult to know whether there is a CDD pole near threshold a priori. Nevertheless in Refs. [17, 15] it is proved that when a CDD pole approaches to the threshold one has the following behaviors for the scattering length and effective range

[TABLE]

with $M_{\rm CDD}$ the bare CDD pole mass and $g_{\rm CDD}$ the residue. According to Eq. (5), one can infer that there exists a CDD pole near the threshold only for the situations with $|a|\ll 1$ fm and $|r|\gg 1$ fm. In this situation, one should use the formalism developed in Ref. [17] to proceed, instead of Eq. (3).

In the present work we first blindly use the ERE formalism in Eq. (3). If the resulting $a$ and $r$ have natural values of the long-range hadronic scale at $1/m_{\pi}\sim\,1$ fm, one could then safely conclude that the formalism in Eq. (3) is applicable in our study (with no indication of a near-threshold CDD pole). We demonstrate below that the resulting values of $a$ and $r$ from the $S$-wave $\Sigma_{c}\bar{D}$ scattering around $P_{c}(4312)$ and the $S$-wave $\Sigma_{c}\bar{D}^{*}$ scattering around the $P_{c}(4440)$ and $P_{c}(4457)$, indeed have typical long-range hadronic scale around $1$ fm. Another issue that needs to be clarified is that we implicitly assume a definite isospin number for the $P_{c}$ states (although we do not need to specify it), otherwise we had to use a coupled-channel scattering formalism in the ERE study. Regarding the quantum numbers of $J^{P}$, the negative $P$ parity can be uniquely fixed in the $S$-wave $\Sigma_{c}\bar{D}^{(*)}$ scattering. Similarly, the total angular momentum is $J=1/2$ for $\Sigma_{c}\bar{D}$ $S$-wave scattering, while for the analogous $\Sigma_{c}\bar{D}^{*}$ case there are two possibilities, $J=1/2$ or $3/2$, which can not be pinned down from our study.

For a resonance pole, its position $E_{R}$ is denoted as

[TABLE]

where $M_{R}$ is the resonance mass and $\Gamma_{R}$ denotes its width. The resonance poles lie on the second Riemann sheet (RS) of the scattering amplitude $T_{II}(E)$, which is given by

[TABLE]

We mention that the convention ${\rm Im}k>0$ should be taken in Eqs. (3) and (7). Given the mass and width of the resonance, we can determine the scattering length $a$ and effective range $r$ by requiring that $T_{II}(E_{R})^{-1}=0$, i.e.

[TABLE]

where $k_{R}$ is the corresponding three-momentum at the pole position

[TABLE]

By solving Eq. (8), it is straightforward to determine the values of $a$ and $r$ once the masses and widths of the resonances are given, with the result [15]

[TABLE]

where $k_{r}={\rm Re}\,k_{R}$ and $k_{i}={\rm Im}\,k_{R}$. As mentioned above in the Introduction, we distinguish the different charged states $\Sigma_{c}^{+}$ and $\Sigma_{c}^{++}$, in order to quantify the isospin breaking effects. The thresholds of the different charged states are explicitly given in the fourth column of Table 1. The results for the scattering lengths $a$ and effective ranges $r$ with uncertainties are collected in the fifth and sixth columns of Table 1, respectively.

According to the values obtained for $a$ and $r$ in Table 1, although we see some discrepancies in the central values for the channels with different charged states, they are compatible after taking into account the uncertainties. It implies that the isospin breaking effects in the three $P_{c}$ states seem mild and further experimental reduction of the uncertainties could help to identify the roles of the isospin breaking.

All of the resulting scattering lengths $a$ and effective ranges $r$ in Table 1 have natural values of the order of $1$ fm, indicating that indeed there is no need for introducing CDD poles near the thresholds. Let us notice that this outcome is consistent with the application of Eq. (3) in our study. Furthermore, the natural values of the $a$ and $r$ allow us to qualitatively conclude that the $P_{c}(4312)$ can be described as an $S$-wave $\Sigma_{c}\bar{D}$ molecular state, and the $P_{c}(4440)$ and $P_{c}(4457)$ are $S$-wave $\Sigma_{c}\bar{D}^{*}$ composite states. Nevertheless, in the ERE approach we can not use the prescription in Ref. [19] to give a quantitative estimate of the probabilities of the $\Sigma_{c}\bar{D}$ component in the $P_{c}(4312)$ and of the $\Sigma_{c}\bar{D}^{*}$ component in the $P_{c}(4440)$ and $P_{c}(4457)$ resonances. In Ref. [19] it has been demonstrated that one can only give the probabilistic interpretation of the compositeness coefficients when the resonance pole $s_{R}=E_{R}^{2}$ lies in an unphysical RS that is directly connected to the physical one in the region $s_{k}<s<s_{k+1}$, such that $s_{k}<{\rm Re}s_{R}<s_{k+1}$, with $s_{k}$ and $s_{k+1}$ the two nearby thresholds. In the single-channel scattering case, it requires that the resonance pole mass should lie above the threshold in the second RS. However in most of cases the pole positions of the $P_{c}$ states in Table 1 are below the thresholds. This fact refrains us from discussing the probabilities of finding the two-particle components in the $P_{c}$ states in the ERE approach.

In the present formalism we are assuming that the whole width of a resonance is due to the corresponding $\Sigma_{c}\bar{D}^{(*)}$ channel, and the resulting $a$ and $r$ are real. On general grounds, because of the presence of the inelastic channels below threshold, like the $J/\psi p$ one to which these resonances decay, the ERE parameters $a$ and $r$ are complex. One possible way to proceed is to include explicitly the inelastic channels below the $\Sigma_{c}\bar{D}^{(*)}$ channel, such as the aforementioned $J/\psi p$. However, in the coupled-channel scattering case, there would be needed extra scattering input which is beyond the scope of the present study. In order to give quantitative information of the inner structures of the $P_{c}$ states, and take into account at least one inelastic channel, we proceed the study by relating the compositeness coefficients with the partial decay widths in next section.

3 Compositeness relations and the partial widths

As mentioned previously, we can not access the quantitative information of the constituents inside the $P_{c}(4312)$ in the elastic scattering $\Sigma_{c}\bar{D}$ and the $P_{c}(4440)$ and $P_{c}(4457)$ from the elastic $\Sigma_{c}\bar{D}^{*}$ scattering. A straightforward extension is to include the additional $J/\psi p$ channel, in which invariant-mass distribution the different $P_{c}$ resonances are actually detected  [3]. For the two-channel $J/\psi p$ and $\Sigma_{c}\bar{D}^{(*)}$ systems, it is natural to assume that the $P_{c}$ resonances lie in the second RS, which now allows us to exploit the formalism in Ref. [19] to calculate the probabilities of the two-particle components in the $P_{c}$. Analogous study has been carried out for the obsolete $P_{c}(4450)$ state by including the $J/\psi p$ and $\chi_{c1}p$ channels in Ref. [21].

The essential prescription of Ref. [19] to calculate the partial compositeness coefficient $X_{j}$ of a resonance $R$ contributed by the $j$th channel is given by

[TABLE]

where $g_{j}$ denotes the coupling between the two-particle state and the resonance $R$, and the one-loop two-point function $G(s)$ is given by

[TABLE]

This expression can be explicitly integrated out by using a once-subtracted dispersion relation or dimensional regularization (replacing the divergence by a constant), which then reads  [33]

[TABLE]

where

[TABLE]

The evaluation of $G_{j}(s)$ for the $j$th channel in Eq. (11) requires to use the proper masses $m_{1}$ and $m_{2}$ in that channel. In this equation $\partial G_{j}(s_{R})/\partial s$ denotes the partial derivative evaluated at the resonance pole position $s_{R}=E_{R}^{2}=(M_{R}-i\Gamma_{R}/2)^{2}$. Notice that $\partial G_{j}(s_{R})/\partial s$ is independent on the subtraction constant $a(\mu_{g})$ and the regularization scale $\mu_{g}$ in Eq. (13).

In order to fix the two couplings $g_{i=1,2}$, we impose that the decay widths $\Gamma_{R}$ of the $P_{c}$ states are saturated by the two channels $J/\psi p$ and $\Sigma_{c}\bar{D}^{(*)}$. The partial decay width $\Gamma_{1}$ to $J/\psi p$ takes the standard form [34]

[TABLE]

where the relativistic three-momentum $q(M_{R}^{2})$ is

[TABLE]

Since in many cases the masses of the $P_{c}$ resonances are below the thresholds of $\Sigma_{c}\bar{D}^{(*)}$, we introduce a Lorentzian mass distribution to calculate the partial width $\Gamma_{2}$ to the $\Sigma_{c}\bar{D}^{(*)}$ channel as

[TABLE]

To restrict the discussion to the resonance energy region, we set the upper integration limit at $M_{R}+2\Gamma_{R}$ in Eq. (17), as in Ref. [21]. After taking into account Eqs. (15) and (17), the saturation condition of the $P_{c}$ decay widths by the $J/\psi p$ and $\Sigma_{c}\bar{D}^{(*)}$ channels gives

[TABLE]

with $q_{1}$ and $q_{2}$ the three-momenta of the $J/\psi p$ and $\Sigma_{c}\bar{D}^{(*)}$ channels, respectively.

For the resonance poles in the second RS in the coupled-channel $J/\psi p$ and $\Sigma_{c}\bar{D}^{(*)}$ scattering, one can identify the compositeness coefficient $X_{j}$ in Eq. (11) as the probability to find the two-particle state from the $j$th channel in the considered resonance. We mention that within the uncertainties of the masses of the $P_{c}(4312)$ and $P_{c}(4457)$, a tiny portion of their poles lies in the third RS (in which the three-momenta of the two channels flip sign) so that they are continuously connected with the physical RS above the $\Sigma_{c}\bar{D}^{(*)}$ threshold. Nevertheless, due to their closeness to the thresholds, their effects can be covered by the large uncertainties of the $P_{c}$ masses. Therefore we shall only focus on the poles on the second RS in the following.

As a clarification remark, let us notice that in Eq. (11) the coupling is taken constant in the range of masses of the resonance along its Lorentzian mass distribution because of the finite width of the resonance, cf. Eq. (17). In this way, there is a smooth transition in the calculation of $X_{2}$ as the value of the nominal resonance pole mass $M_{R}$ varies from above to below the threshold. This allows us some flexibility in order to bypass the strict requirement that the resonance mass should lie above the thresholds of the channels for which $X_{j}$ is calculated. However, in the elastic ERE approach discussed in Sec. 2, the whole width is accounted for only by the channel explicitly taken into account (the second one in the present coupled-channel study), and the situation is more stringent in this respect [15].

The total compositeness $X$ is the sum of $X_{1}$ and $X_{2}$, with $X_{1}$ the partial compositeness coefficient of the $J/\psi p$ and $X_{2}$ the coefficient of $\Sigma_{c}\bar{D}^{(*)}$. By using Eq. (11), we can obtain

[TABLE]

where $G_{1}^{II}(s)$ stands for the $G(s)$ function on the second RS and it is related to the expression in Eq. (13) through $G^{II}(s)=G(s)+i\sigma(s)/(8\pi s)$.

For a given value of the total compositeness $X$ contributed by the $J/\psi p$ and $\Sigma_{c}\bar{D}^{(*)}$ channels, we can determine the couplings $|g_{1}|$ and $|g_{2}|$ by combining Eqs. (18) and (19). In this way, we can further calculate the partial compositeness coefficients $X_{1,2}$ using Eq. (11) and the partial decay widths $\Gamma_{1,2}$ via Eqs. (15) and (17). In principle the partial widths $\Gamma_{2}$ consist of combinations of the $\Sigma_{c}^{+}\bar{D}^{(*)0}$ and $\Sigma_{c}^{++}D^{(*)-}$ channels, depending on the isospin of the pentaquark states $P_{c}$. Nevertheless, we point out that the method employed is not sensitive to whether we assume a definite isospin for the $P_{c}$ states or not, as long as the same masses of the $\Sigma_{c}\bar{D}^{(*)}$ are taken in Eqs. (17) and (18). The reason is because the couplings squared of the different charged states simply add together in these equations. In order to check the isospin breaking effects, we separately solve Eqs. (18) and (19) by using either the masses of $\Sigma_{c}^{+}\bar{D}^{0}(\bar{D}^{*0})$ or $\Sigma_{c}^{++}D^{-}(D^{*-})$.

Concerning the value of $X$ in Eq. (19), we distinguish three different scenarios. In the first scenario, we assume that the compositeness of the $P_{c}$ states is completely saturated by the $J/\psi p$ and $\Sigma_{c}\bar{D}^{(*)}$channels, that is, we first assume that $X=1$. For each $P_{c}$ state, we separately perform the calculations by using either the masses of $\Sigma_{c}^{+}\bar{D}^{0}(\bar{D}^{*0})$ or $\Sigma_{c}^{++}D^{-}(D^{*-})$. The resulting values of the couplings $|g_{1}|$ and $|g_{2}|$, the partial widths $\Gamma_{1}$ and $\Gamma_{2}$, and the partial compositeness coefficients $X_{1}$ and $X_{2}$ are summarized in Table 2. The first lesson we learn from Table 2 is that the $P_{c}$ couplings $|g_{1}|$ to the $J/\psi p$ channel are much smaller than the couplings $|g_{2}|$ to the $\Sigma_{c}\bar{D}^{(*)}$ channel. The situation for the partial decay widths becomes less clear, since many of them have large uncertainties. In all the cases, the overwhelmingly dominant components of the $P_{c}$ states are found to be the $\Sigma_{c}\bar{D}^{(*)}$, in agreement with our qualitative understanding in Sec. 2 from the values of $a$ and $r$ given in Table 1.

In the next two scenarios, we set the compositeness $X=0.8$ and $0.5$ in Eq. (19). In order not to overload the table, we only show the values obtained by using the masses of $\Sigma_{c}^{+}$ and $\bar{D}^{(*)0}$ in Table 3. The results by using the masses of $\Sigma_{c}^{++}$ and $D^{(*)-}$ are quantitatively similar. All the values in Table 3 show quite similar trends as those in Table 2, with $X_{2}\gg X_{1}$.

The previous discussions rely on the assumption that the decay widths of the $P_{c}$ states are saturated by the $J/\psi p$ and $\Sigma_{c}\bar{D}^{(*)}$ channels. Other decay patterns are also predicted, such as those in Refs. [35, 36], which suggest that the partial decay widths of the pentaquark states to the $\Lambda_{c}\bar{D}^{(*)}$ channels could be more important than to the $J/\psi p$. In order to check the robustness of our conclusion, we include the $\Lambda_{c}\bar{D}^{(*)}$ and $\Sigma_{c}\bar{D}^{(*)}$ channels to perform a similar study. To be specific, we give the results in Table 4 by using the masses of $\Lambda_{c}^{+}\bar{D}^{*0}$ and $\Sigma_{c}^{+}\bar{D}^{*0}$. It is verified that to use the masses of other charged states leads to quantitatively similar results. Since to replace the $\Lambda_{c}\bar{D}^{*}$ channel by the $\Lambda_{c}\bar{D}$ does not lead to qualitatively new trends, we do not explicitly show the corresponding results. Comparing the numbers in Tables 2, 3 and those in Table 4, not only the partial decay widths of the two different sets of dynamical channels are quite similar, but also the compositeness coefficients in the different cases are compatible within uncertainties.

Summarizing, we have studied the newly discovered hidden-charm exotic states $P_{c}(4312)$, $P_{c}(4440)$ and $P_{c}(4457)$ by the LHCb Collaboration [3]. We have first applied elastic effective-range expansion in the $\Sigma_{c}\bar{D}^{(*)}$ channel with the scattering length and the effective range fixed by reproducing the mass and width of every resonance separately. In all the cases one obtains values for these parameters of ${\cal O}(1)$ fm, which supports their interpretation as composite resonances of $\Sigma_{c}\bar{D}^{(*)}$. We have also employed another coupled-channel approach involving the two channels $J/\psi p$ and $\Sigma_{c}\bar{D}^{(*)}$ for each resonance, so that we require the saturation of the total width of the resonance. By assuming some values for the total compositeness coefficients for these two channels, ranging from 0.5 to 1, we conclude that the weight of the $\Sigma_{c}\bar{D}^{(*)}$ channel is much larger than the one for $J/\psi p$, in agreement with the ERE approach. We have also performed similar studies by including alternatively the $\Lambda_{c}\bar{D}^{*}$ and $\Sigma_{c}\bar{D}^{(*)}$ as dynamical channels. The conclusions are basically the same as those obtained in the $J/\psi p$ and $\Sigma_{c}\bar{D}^{(*)}$ channels. Needless to say that more thorough studies are needed, e.g., to disentangle the dynamics giving rise to the two nearby $P_{c}(4440)$ and $P_{c}(4457)$ resonances around the $\Sigma_{c}^{+}\bar{D}^{*0}$ and $\Sigma_{c}^{++}D^{*-}$ thresholds and its possible connection with isospin breaking (more likely in the case of composite resonances).

Acknowledgements

This work is funded in part by the Natural Science Foundation of China under Grants No. 11575052, the Natural Science Foundation of Hebei Province under Contract No. A2015205205, and the MINECO (Spain) and EU grant FPA2016-77313-P.