Theoretical description of the decays $\Lambda_b \to \Lambda^{(\ast)}(\frac12^\pm,\frac32^\pm) + J/\psi$
Thomas Gutsche, Mikhail A. Ivanov, J\"urgen G. K\"orner, Valery E., Lyubovitskij, Vladimir V. Lyubushkin, Pietro Santorelli

TL;DR
This paper provides a theoretical analysis of the decay amplitudes for $ ext{Lambda}_b$ baryon decays into various $ ext{Lambda}^*$ states plus $J/ ext{psi}$, aiding in the identification of hidden charm pentaquark states.
Contribution
It offers a covariant quark model calculation of decay amplitudes for multiple $ ext{Lambda}^*$ states, including excited states with different spins and parities.
Findings
Helicity amplitudes for certain excited states are suppressed compared to ground states.
The analysis helps distinguish between different $ ext{Lambda}^*$ contributions in decay chains.
Supports identification of pentaquark states in $ ext{Lambda}_b$ decays.
Abstract
We calculate the invariant and helicity amplitudes for the transitions where the are -type ground and excited states with quantum numbers . The calculations are performed in the framework of a covariant confined quark model previously developed by us. We find that the values of the helicity amplitudes for the and the are suppressed compared with those for the ground state and the excited state . This analysis is important for the identification of the hidden charm pentaquark states and which were discovered in the decay chain because the cascade decay chain…
| state | current | |
|---|---|---|
| scalar diquark | ||
| pseudoscalar diquark | ||
| vector diquark | ||
| axial-vector diquark |
| 1116 | 1405 | 1890 | 1520 | |
| 0 | 0 | |||
| 0.53 | 0.51 | 0.54 | ||
| 0.47 | 0.47 | 0.44 | ||
| 0 | 0 |
| 1116 | 1405 | 1890 | 1520 | |
| 8.00 | 7.07 | 0.45 | 0.19 |
| , | , | |||
|---|---|---|---|---|
| Quantity | Our results | LHCb Aaij:2013oxa | ATLAS Aad:2014iba | CMS CMS:2016iaf |
| , | , | , | , | , | ||||
|---|---|---|---|---|---|---|---|---|
| Quantity | Our results | LHCb Jurik:2016bdm | Our results | LHCb Jurik:2016bdm | Our results | LHCb Jurik:2016bdm | Our results | LHCb Jurik:2016bdm |
| , | , | , | , | |||
|---|---|---|---|---|---|---|
| Quantity | Our results | LHCb Jurik:2016bdm | Our results | LHCb Jurik:2016bdm | Our results | LHCb Jurik:2016bdm |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
MITP/17-036 (Mainz)
Theoretical description of the decays
Thomas Gutsche
Institut für Theoretische Physik, Universität Tübingen,
Kepler Center for Astro and Particle Physics,
Auf der Morgenstelle 14, D-72076, Tübingen, Germany
Mikhail A. Ivanov
Bogoliubov Laboratory of Theoretical Physics,
Joint Institute for Nuclear Research, 141980 Dubna, Russia
Jürgen G. Körner
PRISMA Cluster of Excellence, Institut für Physik, Johannes Gutenberg-Universität,
D-55099 Mainz, Germany
Valery E. Lyubovitskij
Institut für Theoretische Physik, Universität Tübingen,
Kepler Center for Astro and Particle Physics,
Auf der Morgenstelle 14, D-72076, Tübingen, Germany
Departamento de Física y Centro Científico Tecnológico de Valparaíso (CCTVal), Universidad Técnica Federico Santa María, Casilla 110-V, Valparaíso, Chile
Department of Physics, Tomsk State University, 634050 Tomsk, Russia
Laboratory of Particle Physics, Mathematical Physics Department, Tomsk Polytechnic University, 634050 Tomsk, Russia
Vladimir V. Lyubushkin
Dzhelepov Laboratory of Nuclear Problems,
Joint Institute for Nuclear Research, 141980 Dubna, Russia
Pietro Santorelli
Dipartimento di Fisica, Università di Napoli Federico II, Complesso Universitario di Monte S. Angelo, Via Cintia, Edificio 6, 80126 Napoli, Italy
Istituto Nazionale di Fisica Nucleare, Sezione di Napoli, 80126 Napoli, Italy
Abstract
We calculate the invariant and helicity amplitudes for the transitions , where the are -type ground and excited states with quantum numbers . The calculations are performed in the framework of a covariant confined quark model previously developed by us. We find that the values of the helicity amplitudes for the and the are suppressed compared with those for the ground state and the excited state . This analysis is important for the identification of the hidden charm pentaquark states and which were discovered in the decay chain because the cascade decay chain involves the same final state.
relativistic quark model, light and bottom baryons, charmonium, decay rates and asymmetries
pacs:
12.39.Ki,13.30.Eg,14.20.Jn,14.20.Mr
I Introduction
Recently the LHCb Collaboration has performed an angular analysis of the decay , where the ’s are produced in collisions at TeV at the LHC (CERN) Aaij:2013oxa . They reported on the measurement of the relative magnitude of the helicity amplitudes in the decay by a fit to several asymmetry parameters in the cascade decay distribution and .
In an earlier paper Gutsche:2013oea we have performed a detailed analysis of the decay process within a covariant quark model. We have worked out two variants of the threefold joint angular decay distributions in the cascade decay for polarized and unpolarized decays. We have further listed results on helicity amplitudes which determine the rate and the asymmetry parameters in the decay processes and .
In this paper we calculate the corresponding invariant and helicity amplitudes in the transitions where the are -type ground and excited states with quantum numbers . The calculations are performed in the framework of our covariant confined quark model developed previously by us. We find that the values of the helicity amplitudes for the transitions are suppressed compared with those for the transitions to the ground state also calculated in Gutsche:2013oea and the excited state . This analysis is important for the identification of the hidden charm pentaquark states and since the cascade decay involves the same final states as the decay . The subject of the hidden charm pentaquark states has been intensively discussed in the literature (see e.g. Roca:2015tea -Bayar:2016ftu ).
Our paper is structured as follows. In Sec. II, we give explicit expressions for the hadronic matrix elements in terms of dimensionless invariant form factors . The corresponding vector and axial helicity amplitudes are linearly related to the invariant form factors. The linear relations are explicitly calculated and listed. The helicity amplitudes are the basic building blocks in the calculation of the rate and in the construction of the full angular decay distributions for the cascade decays.
In Sec. III, we construct local interpolating three-quark currents corresponding to the states with parity . We then use nonlocal variants of the local interpolating currents to evaluate all invariant amplitudes in the framework of the covariant confined quark model. In Sec. IV, we give numerical results for the normalized helicity amplitudes and branching ratios. Finally, in Sec. V, we summarize our findings.
II The decays
: matrix element and helicity amplitudes
The matrix element of the exclusive decay is defined by (in the present application the vector meson label stands for the )
[TABLE]
where and are the mass and the leptonic decay constant of the vector meson , and and are Cabibbo-Kabayashi-Maskawa (CKM) matrix elements. The coefficient stands for the combination of Wilson coefficients
[TABLE]
The color factor will be set to zero such that we only keep the leading term in the expansion. We take the numerical values of the Wilson coefficients from Altmannshofer:2008dz :
[TABLE]
The hadronic matrix element is expressed in terms of six and eight, respectively, dimensionless invariant form factors viz.
Transition :
[TABLE]
Transition :
[TABLE]
Transition :
[TABLE]
Transition :
[TABLE]
where and all matrices are defined as in the text book by Bjorken-Drell. We use the same notation for the form factors in all transitions even though their numerical values differ. For completeness we have kept the form factors in the transitions and in the transitions although they do not contribute to the decay since .
Next we express the vector and axial helicity amplitudes in terms of the invariant form factors , where and are the helicity components of the vector meson and the daughter baryon , respectively. Note again that the time-component helicity amplitudes do not contribute since the is a spin 1 meson. We need to calculate the expressions
[TABLE]
where we split the helicity amplitudes into their vector and axial parts. We shall work in the rest frame of the parent baryon with the daughter baryon moving in the positive -direction: , and . In this case . Following Ref. Faessler:2009xn one has
Transition : and .
[TABLE]
Transition : and .
[TABLE]
Transition : and .
[TABLE]
Transition : and .
[TABLE]
We use the abbreviations , , .
For the decay width one finds
[TABLE]
The sum over helicities includes all helicities satisfying the angular momentum constraint . Compared to Eq. (11) of Gutsche:2013oea we have dropped a factor containing the lepton mass in the rate expression.
Using the helicity amplitudes one can write down angular decay distributions in the cascade decays where and are the final baryon (, , etc.) and meson (, , etc.) states. Note that the decays are strong and therefore parity conserving while the decay is a weak decay and therefore parity violating. The angular decay distribution involving the strong decays can be obtained from that involving the weak decay by setting the relevant asymmetry parameter to zero. When the is polarized the angular decay distributions are characterized by three polar angles and two azimuthal angles. The full five-fold angular decay distribution can be found in Lednicky:1985zx ; Bialas:1992ny ; Kadeer:2005aq . Corresponding three-fold polar angle distributions for polarized decay and a threefold joint decay distribution for unpolarized decay can be obtained from the full five-fold decay distributions written down in Bialas:1992ny ; Kadeer:2005aq by appropriate angular integrations or by setting the polarization of the to zero.
As mentioned in the introduction there are two variants of how the angular decay distributions of such cascade decay processes can be presented. The unprocessed form the angular decay distribution is written down directly from the traces of the production and the rotated decay spin density matrices. In the present case describes the relative orientation of the decay (or ), the relative orientation of the leptonic decay and the polar orientation of the polarization of the . In the normal form one subtracts off unity corresponding to the normalized total rate, i.e. one writes
[TABLE]
where . The normal form of the threefold angular decay distribution can be written in terms of three linear combinations of normalized squared helicity amplitudes which are Aaij:2013oxa ; Gutsche:2013oea
[TABLE]
where . The last bracketed contribution in (8) only comes in for the case.
III The form factors in the
covariant quark model
We employ generic three-quark currents to describe the states:
[TABLE]
Here or , the color index is denoted by and is the charge conjugation matrix. The Dirac matrices and are chosen in such a way to provide the correct P-parity for the -baryons. A set of currents for the flavor-antisymmetric diquark states is shown in Table 1 by analogy with the classification given in Ref. Nielsen:2009uh .
One can then construct local three-quark currents with the appropriate quantum numbers of the the states. They are given by
[TABLE]
Note that we do not employ derivative couplings in our interpolating currents. It would be interesting to find out whether the use of derivative couplings would change our results.
The covariant quark model employs nonlocal renditions of the local three-quark currents in Eq. (11). The nonlocal Lagrangian describing the couplings of the baryons () with their constituent quarks is then given by
[TABLE]
where and is the mass of the quark placed at the space-time point .
First of all, one has to calculate the mass functions (or self-energy functions) arising from the interactions of the baryons with the constituent quarks as written down in Eq. (12). Then one can determine the coupling constants by using the so-called compositeness condition.
The Fourier-transforms of the mass functions are given by
[TABLE]
The calculation of the Fourier-transforms of the vertex functions can be done in a straightforward way by using Jacobi coordinates. One arrives at the following expressions
[TABLE]
where the sign stands for the final baryon states with and and the sign stands for the final baryon states with and . We have omitted some unnecessary indices and self-explanatory notation.
In the numerical calculations we choose a simple Gaussian form for the vertex functions (for both mesons and baryons):
[TABLE]
where is a size parameter describing the distribution of the quarks inside a given hadron and . We emphasize that the Minkowskian momentum variable turns into the Euclidean form needed for the appropriate falloff behavior of the correlation function (15) in the Euclidean region.
The compositeness condition implies that the renormalization constant of the hadron wave function is set equal to zero. This condition has been suggested by Weinberg Weinberg:1962hj and Salam Salam:1962ap (for a review, see Hayashi:1967hk ) and extensively used in our approach (for details, see Efimov:1993ei ). In the case the compositeness condition may be written in the form
[TABLE]
where is the derivative of the mass function taken on the mass shell . In the case one has to account for the Rarita-Schwinger conditions and . This can be done by splitting off a scalar function in the form . The compositeness condition for the case reads
[TABLE]
In practice, it is more convenient to use a form equivalent to Eqs. (16) and (17) by writing
[TABLE]
The loop integrations in Eq. (14) are performed by using the Fock-Schwinger representations of the quark propagators. The tensorial loop integrations and the manipulations with Dirac matrices are performed with the help of FORM Vermaseren:2000nd . The final relations needed for the determination of the coupling constants may be symbolically cast in the form
[TABLE]
The infrared cutoff parameter provides for the absence of all constituent quark threshold singularities. The threefold integrals are calculated by a FORTRAN code using the NAG library.
The matrix elements of the transitions finally read
[TABLE]
where or , and . The reduced quark masses are defined by
[TABLE]
Below we show the different Dirac structures and in (20) for the different final state baryons:
[TABLE]
The expressions for the scalar form factors are represented by the fourfold integrals
[TABLE]
The model parameters are the constituent quark masses and the infrared cutoff parameter responsible for quark confinement. They are taken from a new fit done and used in our papers on the semileptonic decays Ivanov:2017mrj ; Ivanov:2016qtw ; Ivanov:2015woa ; Ivanov:2015tru , rare decays Issadykov:2015iba ; Dubnicka:2015iwg ; Dubnicka:2016nyy , the semileptonic decays , Gutsche:2015mxa ; Gutsche:2015rrt and for the calculation of nucleon tensor form factors Gutsche:2016xff . The best fit values for the constituent quark masses and the infrared cutoff parameter are
[TABLE]
The dimensional-size parameters of the ground-state and baryons have been determined by a fit to the semileptonic decays and . The resulting values are , GeV. The values of the size parameters of the final states are set equal to the size parameter of the ground state . For the size parameter of the we take GeV as determined from our most recent fit described above.
IV Numerical results
In Tables 2 and 3 we list our predictions for the normalized helicity amplitudes and branching fractions.
The helicity amplitudes of the produced states are clearly dominated by the helicity configuration as in the quark level transition . For the spin states in the transition this implies that the two states are almost purely left-handed.
It is also apparent that the branching ratios involving the excited states are suppressed relative to those of the ground state and the excited state with .
Next, we discuss the experimental results on the pentaquark evidence in light of our theoretical findings. Large data samples of the bottom baryon state were collected by the LHC experiments from collisions during Run I. An angular analysis of the decay was first done by the LHCb Collaboration. Aaij:2013oxa and was then repeated by the ATLAS Aad:2014iba and CMS CMS:2016iaf Collaborations. From Table 4, one can assess the present-day accuracy of the measured helicity amplitudes for the transition .
The dominant production mechanism of the ’s at the LHC proceeds via the strong interactions. Therefore the longitudinal polarization of the produced vanishes because of parity conservation. Contrary to this a transverse polarization component is not forbidden by parity. The component depends strongly on the Feynman variable , where is the longitudinal momentum of the and is the collision center-of-mass energy. For the collisions of identical unpolarized initial state particles one has by virtue of the invariance under the rotation of the coordinate system through an angle of about the normal to the reaction plane Abramov:2005sb . This implies that . Taking into account the very small value for ’s produced at the LHC at TeV, is estimated to be in Hiller:2007ur . The value measured in Aaij:2013oxa ; Aad:2014iba ; CMS:2016iaf is consistent with the expected value of zero. We therefore treat to be unpolarized in our further analysis.
Large samples of the decay decay have been collected by the LHCb experiment Aaij:2015tga . This decay was expected to be dominated by resonances decaying into final states. The measured fit fractions of the and states are and , respectively. It was also found that the data cannot be satisfactorily described without the inclusion of two Breit-Wigner resonances decaying strongly to . These new pentaquark states, called and , have large fit fractions of and of the total sample, respectively.
The presence of various conventional resonances and exotic pentaquark states , which can interfere with each other, makes the analysis of experimental data particularly difficult and challenging. There are 32 complex parameters ( amplitudes) which describe the transition . The 32 parameters were determined in a six-dimensional fit Aaij:2015tga to the angular decay distribution of the cascade decay process. Their numerical values for the default fit variant can be found in the Appendix G of Jurik:2016bdm .
The helicity amplitudes used in the present approach are linearly related to the -amplitudes used in Aaij:2015tga (see e.g. ms1970 ; chung2014 ). The coefficients of this linear relation can be obtained with the help of angular momentum Clebsch-Gordan coefficients (see Eq. (2) in Aaij:2015tga ). There are two important remarks that one has to make here. First, the amplitudes from the fit have already been redefined in order to account for the helicity couplings from the strong decays , i.e. the experiment reports values for the products . It is obvious that the additional factor still allows one to compare the moduli of normalized helicity amplitudes with theoretical expectations. The second remark concerns the unusual sign of the helicity in the decay chain in Jurik:2016bdm leading to a redefinition .
In Tables 5 and 6, we compare our predictions for the moduli squared of normalized helicity amplitudes with values recalculated from a fit to experimental data Jurik:2016bdm . Unfortunately, the absence of the correlation matrix does not allow us to estimate the error bands for those values; nevertheless, they are expected to be rather large. We found that the predicted values of do not vary significantly with the invariant mass and the leading contribution should come from the and helicity amplitudes for both the and cases. No such pattern could be found in the experimental data (may be except for ).
At this stage we cannot decide definitely whether our theoretical approach contradicts the experimental analysis or not. A refit of the LHCb data using our theoretical restrictions for the helicity amplitudes is needed to obtain an unambiguous conclusion.
Despite this fact, we propose a simple check. For this purpose we have generated a sample of 750k events containing the decay using the PYTHIA 8.1 Sjostrand:2007gs Monte Carlo generator. The decay products are distributed isotropically in the rest frame. We have used an event-by-event reweighting to take into account the presence of resonances as well as possible interference effects between them. The weight is given by the matrix element squared and is calculated for each simulated decay using the 4-momenta of the outgoing particles. We have tried to reproduce the default LHCb fit which is known to describe the experimental data well by using the corresponding matrix elements (Eq. (8) from Aaij:2015tga ) with the appropriate constants given in Jurik:2016bdm .
The invariant mass distributions for the and invariant masses are shown in Fig. 1. We have concentrated on events with where this subsample contains about 39k of the decays. The contribution from the state should be dominant here while the influence of pentaquark states can be safely neglected.
Next we have tried to consider the helicity angle distributions for , and (our definitions for helicity angles are identical with those of the LHCb analysis). should be trivial for the decay of a unpolarized particle. The shape of the helicity angle distribution is the same as for the single resonance: , with a coefficient .
The most interesting distribution is shown on the left plot in Fig. 2. The black dotted curve corresponds to the full LHCb model. The curve does not look like the expected even function of the cosine of the helicity angle. The main reason is the strong interference between different intermediate resonances. We have tried to describe this distribution by taking into account only the three lowest states (i.e. , and ). We have further neglected all helicity amplitudes except for and ; the overall fraction of each resonance is also fixed to its LHCb values.
In our approach all helicity amplitudes are real. Complex phases can result from the decay through final state interactions. To start with we assign an identical complex phase to each helicity amplitude. The result is shown on Fig. 2 by a solid blue line. Obviously, it does not agree with the reference plot. One can achieve a reasonable agreement by varying the 5 complex phases of the different helicity amplitudes and keeping the moduli of the helicity amplitudes unchanged (see the solid red line in Fig. 2). We take this as evidence that experimental data can in fact be described within our theoretical approach which includes quite strong constraints for the moduli of the helicity amplitudes.
V Summary
We have calculated the invariant and helicity amplitudes in the transitions where the states are -resonances with quantum numbers . The calculations were performed in the framework of a covariant confined quark model previously developed by us. We have found that the values of the helicity amplitudes for the transitions into the and states are suppressed compared with those for the transitions into the ground state and the excited state .
We have compared our numerical results for the helicity amplitudes and decay assymmetry parameters for the set of resonances with those recalculated from the LHCb fit. We have shown that the helicity angle distributions for the low-mass resonances () can be reproduced using our predicted values for the normalized helicity amplitudes and .
This analysis is important for the identification of the hidden charm pentaquark states and since the cascade decay involves the same final states as the decay .
Acknowledgements.
This work was supported by the German Bundesministerium für Bildung und Forschung (BMBF) under Project 05P2015 - ALICE at High Rate (BMBF-FSP 202): “Jet- and fragmentation processes at ALICE and the parton structure of nuclei and structure of heavy hadrons”, by CONICYT (Chile) PIA/Basal FB0821, by the Tomsk State University Competitiveness Improvement Program and the Russian Federation program “Nauka” (Contract No. 0.1764.GZB.2017). M.A.I. acknowledges the support from the PRISMA cluster of excellence (Mainz Uni.). M.A.I. and J.G.K. thank the Heisenberg-Landau Grant for the partial support. P.S. acknowledges support by the INFN (QFT-HEP project).
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1(1) R. Aaij et al. (LH Cb Collaboration), Phys. Lett. B 724 , 27 (2013) [ar Xiv:1302.5578 [hep-ex]].
- 2(2) T. Gutsche, M. A. Ivanov, J. G. Körner, V. E. Lyubovitskij and P. Santorelli, Phys. Rev. D 88 , no. 11, 114018 (2013) [ar Xiv:1309.7879 [hep-ph]].
- 3(3) L. Roca, M. Mai, E. Oset and U. G. Meißner, Eur. Phys. J. C 75 , no. 5, 218 (2015) [ar Xiv:1503.02936 [hep-ph]].
- 4(4) A. Feijoo, V. K. Magas, A. Ramos and E. Oset, Phys. Rev. D 92 , no. 7, 076015 (2015) Erratum: [Phys. Rev. D 95 , no. 3, 039905 (2017)] [ar Xiv:1507.04640 [hep-ph]].
- 5(5) H. X. Chen, L. S. Geng, W. H. Liang, E. Oset, E. Wang and J. J. Xie, Phys. Rev. C 93 , no. 6, 065203 (2016) [ar Xiv:1510.01803 [hep-ph]].
- 6(6) E. Wang, H. X. Chen, L. S. Geng, D. M. Li and E. Oset, Phys. Rev. D 93 , no. 9, 094001 (2016) [ar Xiv:1512.01959 [hep-ph]].
- 7(7) A. Feijoo, V. K. Magas, A. Ramos and E. Oset, Eur. Phys. J. C 76 , no. 8, 446 (2016) [ar Xiv:1512.08152 [hep-ph]].
- 8(8) J. X. Lu, E. Wang, J. J. Xie, L. S. Geng and E. Oset, Phys. Rev. D 93 , 094009 (2016) [ar Xiv:1601.00075 [hep-ph]].
