Roper-like resonances with various flavor contents and their two-pion emission decays

A. J. Arifi1, H. Nagahiro1,2, A. Hosaka1,3, and K. Tanida3 1Research Center for Nuclear Physics (RCNP), Osaka University, Ibaraki, Osaka 567-0047, Japan 2Department of Physics, Nara Women’s University, Nara 630-8506, Japan 3Advanced Science Research Center, Japan Atomic Energy Agency, Tokai, Ibaraki 319-1195, Japan (Dated: April 17, 2020)

∗ We study the three-body decay of the newly observed bottom baryon Λb (6072) by LHCb; ∗ Λb (6072) → Λbππ. Its mass about 500 MeV above the ground state Λb and a broad width im- ply that the state could be an analogue of the Roper resonance of the nucleon N(1440). In terms ∗ of sequential processes going through Σb and Σb , we find that the observed invariant mass distri- bution is reproduced assuming its spin and parity J P = 1/2+. We discuss that the ratio of the two sequential processes and angular correlation of two pions are useful for the determination of spin and parity. We suggest further studies for the Roper resonance analogue in various flavor contents, raising an interesting and important question in baryon spectroscopy. Keywords: Roper-like resonance, flavor-independent nature, two-pion emission decay

In this letter, motivated by the recent observation of the spin and parity of these charmed and bottom candi-

Λb∗(6072) by CMS and LHCb collaborations [1,2], with dates are not yet determined. Observing this situation, mass and width the M = 6072.3 MeV and Γ = 72 MeV we show possible candidate states in Fig.1 with their ex- measured by LHCb [2], we discuss it as a state analogous citation energies. They are similar not only in masses but to the Roper resonance N(1440) with spin and parity also in decays through two pion emission. These rather J P = 1/2+ in heavy flavor sectors. It is shown that universal features in various flavor sectors may suggest the three-body decay of two pion emission is particularly important dynamics of low energy QCD. useful to determine its unknown spin and parity. The Roper resonance is the first excited state of the In the present letter, focusing on Λb∗(6072), we study → nucleon of J P = 1/2+ [3]. It has been a mysterious state in detail its three-body decay Λb∗(6072) Λbππ. In our because its properties such as the mass and level order- previous preprint [11], we have shown that the Dalitz ing with negative parity nucleons have not been easily plot analysis is particularly useful for the determination explained by the conventional quark model. Recently, of spin and parity of Λc∗(2765) when the decay is domi- detailed studies by a dynamical model for meson-baryon nated by sequential processes going through Σc and Σc∗. scatterings [4,5] and by the measurement of the transi- The relevant diagrams are shown in Fig.2 (for charmed tion form factor in a wide range of momentum transfer baryons, replace the label b for bottom by c). Follow- Q2 [6] have convinced that it is a radial excitation of the ing this, we investigate the sequential decays of Λb∗(6072) nucleon [7]. with an intermediate state Σb or Σb∗. In principle, there is another one of the direct process where the two pions are emitted from a single vertex. As we will see later, however, the sequential processes can explain the so far N(1440) (1600) (1660) c(2765) c(2970) b(6072) ] 0.5

V observed data very well. e G [

M

N c c b + + + 0 ⇡ ⇡ ⇡ ⇡ ⇡ ⇡ arXiv:2004.07423v1 [hep-ph] 16 Apr 2020 (udu) (uds) (uds) (udc) (usc) (udb) + + Quark content 0 0 0 0 0 0 ⇤b⇤ ⇤b ⇤b⇤ ⇤b ⇤b⇤ ⇤b ( ) ( )+ ⌃b⇤ ⌃b⇤ FIG. 1. Excitation energies of the first excited J P = 1/2+

baryons and candidates with various flavor contents. Red bars ∗0 0 + − FIG. 2. Sequential decay processes of Λb → Λb π π going are for those with undetermined spin and parity. (∗)− (∗)+ through Σb and Σb . On the other hand, similar states of J P = 1/2+ with almost the same excitation energies have been known for some time for hyperons systematically [8,9]. The recent To compute these diagrams, let us first introduce

observation of Λb∗(6072) could add another candidate in amplitudes for various vertices. Because the bot- the list of the analogous states. In fact, there is also a tom baryons are sufficiently heavy, we employ a non- candidate in the charm sector, Λc∗(2765), which was re- relativistic method. Assuming that the spin and parity + cently identified as an isoscalar particle [10]. However, of Λb∗ equal 1/2 , the amplitudes at the first vertex of 2

Fig.2 are given by s-wave decay. Therefore, we find R  1 which does not seem consistent with the experimental data. This ∗ ∗ −iTΛ Σbπ = g1 χ† (σ · p)χΛ , b → Σb b demonstrates the advantage in the use of the ratio R for ∗ ∗ ∗ −iTΛ Σ π = g2 χΣ† ∗ (S† · p)χΛ , (1) spin and parity. b → b b b Now we discuss Dalitz plots and other related quan- and for the second vertex, tities. In what follows, we assume again J P = 1/2+ for Λ∗(6072). For decays of unpolarized particles, Dalitz −iTΣb Λbπ = g3 χ† (σ · p) χΣb , b → Λb plots are expressed as functions of two kinematical vari- ∗ ∗ −iTΣ Λbπ = g4 χ† (S · p) χΣ . (2) b → Λb b ables formed by the three decaying particles. In fact, they In Eqs. (1) and (2), g ’s are the coupling constants and also depend on the mass of the initial particle, which i ∼ χ’s are the two-component spinors for the baryons as distributes over a finite range of width 72 MeV for indicated by lower indices. Λb∗(6072). In experiments with low statistics, plots are By using these vertex amplitudes, the amplitude of often made by integrating the signals from such a finite mass range. Therefore, to compare directly theory pre- the first diagram of Fig.2 with the intermediate Σ b− is expressed as dictions with the actual data, this feature is properly treated. TΣ− Λ0π− TΛ∗0 Σ−π+   b → b b → b In the upper three panels of Fig.3 shown are diagonal −iT Σb− = −i i , (3) m23 − m − + Γ − views of the Dalitz plots for decay probabilities at three Σb 2 Σb different fixed masses of Λb∗(6072) around the peak value, where m23 is the invariant mass of the subsystem m(Λbππ) = 6.052, 6.072 and 6.092 GeV, as functions 0 2 2 2 2 π−(particle 2) and Λb (particle 3). The amplitudes with of invariant mass squares m and m , P (m; m , m ). + 12 23 12 23 the Σb intermediate state and those with Σb∗∓ are com- Note that the invariant mass m23 is for Λbπ− that couples ( ) puted similarly. The total amplitude is then the coherent to Σ ∗ −, which will be relevant in the discussions below. sum of the four of them. b The unknown coupling constants gi’s are constrained by the heavy quark spin symmetry of QCD [12]; the ra- tios of g1 to g2 and of g3 to g4 are computed when Σb and Σb∗ are regarded as heavy quark spin partners. In the heavy quark limit, the spin J of a baryon is com- posed of a spin j of the brown muck of light degrees of freedom [13] and spin 1/2 of the heavy quark, such that J = j ± 1/2.

For a spin 1/2 Λb∗, the brown muck spin can take two values, j = 0, 1. Assume further that the brown muck ( ) of Σb∗ has spin j = 1 as in the quark model, their total spin J can be 1/2 or 3/2. Having this information, we can compute the ratio as g √ 1 2 = 2 for j = 0 and √ for j = 1. (4) g1 2 Combined with two-body phase space factors, when the same partial L-wave is possible for the two final states

πΣb and πΣb∗, we find the ratio Γ(Λ (6072) → Σ π) g2 p2L+1 R = b∗ b∗ = 2 2 . (5) → 2 2L+1 Γ(Λb∗(6072) Σbπ) g1 p1 For J P = 1/2+, the two decays are in p-wave, L = 1. Using p1 = 215 MeV and p2 = 192 MeV, we obtain R = 1.43 for j = 0 and 0.36 for j = 1 (6) As we will see shortly, this agrees well with the observed P data from LHCb. On the other hand, if J = 1/2− the partial waves of πΣb and of πΣ∗ are s and d-waves, FIG. 3. The Dalitz and invariant mass plots at three differ- b ∗ respectively. For low momentum pion as in the present ent initial masses of Λb (6072) (upper three panels) and the case, the d-wave decay is suppressed as compared to the convoluted one (bottom panel). 3

In each layer, the solid red line shown in the back panel is the projected strength, namely the invariant mass plot. In all the three plots, we find four resonance bands corre- + + sponding to Σb−,Σb∗−,Σb∗ and Σb from the left to right. The first two are for the first diagram of Fig.2 showing the genuine resonance band, while the latter two corre- spond to the second one that are the so called kinematical reflection. Now what we need is the sum of the plots like these three ones. More precisely, the convolution is performed over the finite mass range weighted by the Breit-Wigner function,

Z 2 2 ∗ 2 2 1 P (m; m12, m23) dm FIG. 5. The convoluted square Dalitz plot for Λb (6072) and ∗− P (m12, m23) = 2 2 , (7) its corresponding angular correlation along Σ with a mass − ∗ ∗ b N (m MΛb ) + ΓΛ /4 b cut M ∗− ± Γ ∗− . Σb Σb where the factor N is for normalization and m for the Λbππ mass. The result is shown in the bottom panel of Fig.3 as well as the invariant mass plot on the back states. If only Λbπ− is included, the left peak is slightly panel. As shown there, in the convoluted plot the peaks higher than that of the right one as shown in the bot- of kinematical reflections in the upper three panels are tom panel of Fig.3 and by the green line in Fig.4. This smeared out, while the real resonance peaks remain. analysis so far supports the spin and parity of Λb∗(6072) + Because the experimental data is shown for the sum is 1/2 , and their decay is dominated by the two pion + emission of the sequential processes. This is one of the of signals of Λbπ and Λbπ−, we have shown in Fig.4 the corresponding one of our calculation (red solid line) main conclusions of the present study. If there will be compared with the experimental data. Note that exper- further high statistics data, it is interesting to compare imental resolution of around 1 MeV is not considered in the Dalitz plots at each different energies. the calculation. For theory side, we have also shown the Let us further look at the angular correlation (depen- + dence) along the resonance band. The angle is the one strengths of Λbπ− and Λbπ separately by blue and green solid lines, respectively. We find that our theory calcula- between the two pions θ12 in the intermediate resonance tion agrees remarkably well with experimental data not rest frame. To see the correlation better, it is convenient only in overall shape of the peaks and background. The to make a square plot as a function of m23 and cos θ12 relative strengths of the peaks is determined by the ratio as shown in the lower plot of Fig.5. As seen there the R, while the background shape is reproduced by the con- angular correlation along the Σb band (near side) is flat voluted kinematical reflection. It is interesting to see that due to its spin 1/2. On the other hand, the angular cor- even such a detailed structure, the relative height of the relation along the Σb∗ band (far side) has a concave shape due to its spin 3/2. two resonance peaks of Σb and Σb∗ are well reproduced. This is achieved by taking the sum of the two charged In general the angular correlation along the Σb∗ band decaying from a particle of spin J contains the two terms weighted by the helicity amplitudes Ah, where h is the conserved helicity in the decay, 0 270 b 2 2 0 + W (θ12) ∝ A (Λ∗ → Σ∗π) × (1 + 3 cos θ12) b 1/2 b b Total 2 2 + A3/2(Λb∗ → Σb∗π) × 3 sin θ12. (8) 180 LHCb When the initial Λb∗ has spin 1/2, the helicity amplitude A3/2 vanishes. Thus, the angular correlation will have 90 2 a 1 + 3 cos θ12 dependence (concave shape). In reality, interference between the resonance of a finite width and background (kinematical reflection in the present case) 0 contaminates that angular correlation. Fitting the pro- 5.75 5.80 5.85 5.90 5.95 6.00 jected angular correlation as plotted in the upper panel m( 0 ±) [GeV] b of Fig.3, our prediction for the angular correlation is 2 0 − 0 + W (θ12) ∝ 1 + 3.3 cos θ12. (9) FIG. 4. Invariant mass plots of Λb π (blue), Λb π (green) and their sum (red). The data from LHCb is the sum of these So far we have focused our discussions on Λb∗(6072). two [2] and compared with the red line. As anticipated in the beginning of this paper, there are 4 other candidate baryons of similar nature, in particular Ξ and by Grants-in Aid for Scientific Research on Innova- baryons. One is Ξc∗(2970) with excitation energy around tive Areas (No. 18H05407). 500 MeV and decays into Ξcππ [14]. If it is regarded as a Roper-like resonance, its decay into ΛcK is forbidden due to the brown muck selection rule, namely the transition + + [1] A. M. Sirunyan et al. [CMS Collaboration], Phys. Lett. of 0 (initial brown muck) → 0 (final brown muck) +0− B 803, 135345 (2020) (kaon) does not occur. This selection rule applies to the [2] R. Aaij et al. [LHCb Collaboration], arXiv:2002.05112 other candidates and is a unique feature of the Roper like [hep-ex] resonance. [3] L. D. Roper, Phys. Rev. Lett. 12, 340 (1964) In recent analysis, LHCb reported several resonances [4] N. Suzuki, B. Julia-Diaz, H. Kamano, T. Lee, A. Mat- in Λ K [15]. One of them is Ξ (2965) with similar excita- suyama and T. Sato, Phys. Rev. Lett. 104, 042302 (2010) c c∗ [5] H. Kamano, S. Nakamura, T. Lee and T. Sato, Phys. tion energy, but it has significantly smaller decay width Rev. C 81, 065207 (2010) than Ξc∗(2970). Therefore, Ξc∗(2965) could be a particle of [6] I. Aznauryan et al. [CLAS], Phys. Rev. C 80, 055203 different nature that can be identified as a negative parity (2009) excitation of Ξc0 as pointed out in Ref. [16–18]. Likewise, [7] V. D. Burkert and C. D. Roberts, Rev. Mod. Phys. 91, 011003 (2019) in the bottom sector, the Ξb∗(6227) resonance observed in [8] M. Tanabashi et al. [Particle Data Group], Phys. Rev. D ΛbK and Ξbπ decays [19] would be this analogue of neg- 98, 030001 (2018) and 2019 update ative parity excitation of Ξb0 [20–24]. We expect a Roper P + [9] M. Takayama, H. Toki and A. Hosaka, Prog. Theor. Phys. like Ξb∗ of J = 1/2 decaying into Ξbππ at around 6.2 101, 1271 (1999) GeV. [10] A. Abdesselam et al. [Belle Collaboration], We may expect further candidates in the strangeness arXiv:1908.06235 [hep-ex]. sector. For example, the resonance Ξ∗(1820) listed in [11] A. J. Arifi, H. Nagahiro, A. Hosaka and K. Tanida, arXiv:2003.08202 [hep-ph] PDG [8] may be identified with 3/2− resonance as ob- served in ΛK invariant mass [25]. However, a resonance [12] N. Isgur and M. B. Wise, Phys. Rev. Lett. 66, 1130 (1991) with a similar excitation energy is also observed in Ξππ [13] M. Neubert, Phys. Rept. 245, 259 (1994) invariant masses with a larger decay width [26]. This is [14] J. Yelton et al. [Belle], Phys. Rev. D 94, 052011 (2016) a candidate of the Roper like resonance. [15] R. Aaij et al. [LHCb], arXiv:2003.13649 [hep-ex] Here we summarize several common features of the [16] H. Yang, H. Chen and Q. Mao, arXiv:2004.00531 [hep- Roper-like resonances in our observation: ph] [17] Q. Lu, arXiv:2004.02374 [hep-ph] • The excitation energy is around 500 MeV which is [18] K. Wang, L. Xiao and X. Zhong, arXiv:2004.03221 [hep- ph] rather flavor-independent. [19] R. Aaij et al. [LHCb], Phys. Rev. Lett. 121, 072002 (2018) • The decay width is rather large with a significant [20] B. Chen, K. Wei, X. Liu and A. Zhang, Phys. Rev. D 98, coupling to two-pion emission decay. 031502 (2018) [21] E. Cui, H. Yang, H. Chen and A. Hosaka, Phys. Rev. D • Their two-pion emission decays are dominated by 99, 094021 (2019) the sequential processes going through 1/2+ reso- [22] K. Wang, Q. L and X. Zhong, Phys. Rev. D 99, 014011 (2019) e.g. + e.g. nance, Σb and 3/2 resonance, Σb∗, or [23] H. Yang, H. Chen, E. Cui, A. Hosaka and Q. Mao, Eur. their analogues, indicating that the direct contri- Phys. J. C 80, 80 (2020) butions including from f0(500) are insignificant as [24] D. Jia, W. Liu and A. Hosaka, Phys. Rev. D 101, 034016 observed also for other baryons [27–29]. (2020) [25] S. Biagi, M. Bourquin, R. Brown, H. Burckhart, P. Ex- • The ratio R and angular correlation can be used to termann, M. Gailloud, C. P. Gee, W. Gibson, P. Jacot- determine their spin and parity 1/2+. Guillarmod, J. Perrier, K. Ragan, P. Rosselet, P. Schi- rato, H. Siebert, V. Smith, K. Streit, J. Thresher, A. Wood and C. Yanagisawa, Z. Phys. C 34, 175 (1987) Observations in the heavy flavor baryons and revisit to [26] S. Apsell, N. Barash-Schmidt, L. Kirsch, P. Schmidt, the light flavor baryons could provide new insights into C. Chang, R. Hemingway, B. Khoury, A. Stottlemyer, unique flavor independent dynamics of low energy QCD. H. Whiteside, G. Yodh, S. Glickman, M. Goldberg, S. Ja- We thank support from the Reimei Research Promo- cobs, K. Jaeger, C. Mccarthy, B. Meadows, G. Moneti, tion project (Japan Atomic Energy Agency) in comple- J. Bartley, R. Dowd, J. Schneps and G. Wolsky, Phys. tion of this work. A. J. A. is supported by a scholar- Rev. Lett. 24, 777 (1970) [27] S. Prakhov et al. [Crystal Ball], Phys. Rev. C 69, 045202 ship from the Ministry of Education, Culture, Science (2004) and Technology of Japan. H. N. is supported in part [28] S. Prakhovet al. [Crystal Ball], Phys. Rev. C 69, 042202 by Grants-in Aid for Scientific Research, Grants No. (2004) 17K05443(C) and A. H. by Grants No. 17K05441(C) [29] R. Mizuk et al. [Belle Collaboration], Phys. Rev. Lett. 98, 262001 (2007)