1 Introduction
Charm baryon spectroscopy provides an excellent ground for studying the dynamics of light quarks in the environment of a heavy quark. In the past decade, many new excited charmed baryon states have been discovered by BaBar, Belle, CLEO and LHCb. B decays and the continuum are both very rich sources of charmed baryons. Many efforts have been made to identify the quantum numbers of these new states and understand their properties.
Consider the strong decays and , where . The mass differences between and and between and in the charmed baryon sector were found to be large enough that the strong decays of and are kinematically allowed. Consequently, the charmed baryon system offers a unique and excellent laboratory for testing the ideas and predictions of heavy quark symmetry of the heavy quark and chiral symmetry of the light quarks. These tests will have interesting implications for the low-energy dynamics of heavy baryons interacting with Goldstone bosons.
Theoretical interest in hadronic weak decays of charmed baryons peaked around the early 1990s and then faded away. To date, we still do not have a good phenomenological model, not mentioning the quantum chromodynamics (QCD)-inspired approach as in meson decays, to describe the complicated physics of baryon decays. We need cooperative efforts from both experimentalists and theorists to make progress in this arena.
This review is essentially an update of Ref. [
1], which described charmed baryon physics around 2007. The outline of the content is the same as that of Ref. [
1] except that we add discussions of the spectroscopy and lifetimes of doubly charmed baryons.
Several excellent review articles on charmed baryons can be found in Refs. [
2–
7].
2 Spectroscopy
2.1 Singly charmed baryons
The singly charmed baryon is composed of a charmed quark and two light quarks. Each light quark is a triplet of flavor SU(3). There are two different SU(3) multiplets of charmed baryons: a symmetric sextet
6 and an antisymmetric antitriplet
. The
and
form a
representation and they all decay weakly. The
,
,
and
form a
6 representation; among them, only
decays weakly. We follow the Particle Data Group’s (PDG’s) convention [
8] of using a prime to distinguish the
in the
6 from that in the
.
In the quark model, the orbital angular momentum of the light diquark can be decomposed into , where is the orbital angular momentum between the two light quarks, and is the orbital angular momentum between the diquark and the charmed quark. Although the separate spin angular momentum and orbital angular momentum of the light degrees of freedom are not well defined, they are included for guidance from the quark model. In the heavy quark limit, the spin of the charmed quark , and the total angular momentum of the two light quarks , are separately conserved. The total angular momentum is given by . It is convenient to use , , and to enumerate the spectrum of states. Moreover, one can define two independent relative momenta, and , from the two light quark momenta , and the heavy quark momentum . Denoting the quantum numbers and as the eigenvalues of and , respectively, the ρ-orbital momentum describes relative orbital excitations of the two light quarks, and the λ-orbital momentum describes orbital excitations of the center of mass of the two light quarks relative to the heavy quark (see Fig. 1). The p-wave heavy baryon can be in either the λ-state or the ρ-state. The orbital λ-state (ρ-state) is clearly symmetric (antisymmetric) under the interchange of and . In the following, we use the notation () to denote the states that are symmetric (antisymmetric) in the orbital wave functions under the exchange of two light quarks. The lowest-lying orbitally excited baryon states are the p-wave charmed baryons, the quantum numbers of which are listed in Table 1.
The next orbitally excited states are the positive-parity excitations with . There exist multiplets (e.g., and ) with the symmetric orbital wave function, corresponding to and (see Table 2). We use a hat to distinguish them. Because the orbital states are antisymmetric under the interchange of two light quarks, we use a tilde to denote them. Moreover, the notation is reserved for tilde states in the , as the quantum number is needed to distinguish different states.
The observed mass spectra and decay widths of charmed baryons are summarized in Table 3 (see also Fig. 2). Note that except for the parity of the lightest and the heavier one , none of the other JP quantum numbers given in Table 3 has been measured. One has to rely on the quark model to determine the spin-parity assignments.
In the following, we discuss some of the excited charmed baryon states:
2.1.1 states
The lowest-lying
p-wave
states are
,
,
and
. The doublet
is formed by
and
[
14]. The allowed strong decays are
and
. This explains why the width of
is narrower than that of
. Because of isospin conservation in strong decays,
is not allowed to decay into
.
is a broad state first seen in
by CLEO [
15]. However, it is still not known whether it is a
or a
and whether the width might be due to overlapping states. The Skyrme model [
16] and the quark model [
17,
18] suggest a
state with a mass of 2742 and 2775 MeV, respectively. Therefore,
could be the first positive-parity excitation of
. In the diquark model, it has also been proposed to be either the first radial (
) excitation of the
with
containing the light scalar diquark or the first orbital excitation (
) of the
with
containing the light axial vector diquark [
19].
The state
, first observed by CLEO [
15] in
, was also seen by BaBar in the
spectrum [
20]. Belle studied the experimental constraint on the
JP quantum numbers of
[
21] and found that
is favored by the angular analysis of
together with the ratio
, which was measured to be
In the quark model, the candidates for the parity-even spin- state are , , , , and (see Table 2). The first four candidates, with , decay to in an F wave and to in F and P waves. Neglecting the P-wave contribution for the moment,
where the factor of 4/5 follows from heavy quark symmetry. At first glance, this appears to be in good agreement with experiment. However, the channel is available via a P-wave and is enhanced by a factor of relative to the F-wave one. However, heavy quark symmetry cannot be applied to calculate the contribution of the channel to the ratio R, as the reduced matrix elements for the P-wave and F-wave modes differ. In this case, one has to rely on a phenomenological model to compute the ratio R. As for , it decays to , and all in F waves. Further,
Although this deviates from the experimental measurement (2.1) by
, it is a robust prediction. This has motivated us to conjecture that the first positive-parity excited charmed baryon
could be an admixture of
,
and
[
10].
It is worth mentioning that the Peking group [
22] has studied the strong decays of charmed baryons on the basis of the so-called
recombination model. For the
, the Peking group found that (i)
cannot be a radial excitation, as its decay into
is prohibited in the
model if
is the first radial excitation of
, and (ii) the states
,
and
are excluded, as they do not decay to
according to the
model. Moreover, the predicted ratios
are either too large or too small compared to experiment, for example,
Both symmetric states
and
, are thus ruled out. Hence, it appears that
dictates the inner structure of
.
† However, there are several problems with this assignment: (i) the quark model indicates a
state around 2910 MeV, which is close to the mass of
, whereas the mass of
is even higher [
17,
18], (ii)
can decay to an
F-wave
, and this has not been seen by BaBar and Belle, and (iii) the calculated width, 28.8 MeV, is too large compared to the measured one,
MeV. One may argue that the
model’s prediction can easily differ from the experimental measurement by a factor of 2–3 owing to its inherent uncertainties [
22].
Interestingly, the quantum numbers
for the
were correctly predicted on the basis of the diquark concept in Refs. [
24,
25] before the Belle experiment.
The highest state
was first discovered by BaBar in the
decay mode [
20] and confirmed by Belle in the decays
, which subsequently decay into
[
21]. Its spin-parity assignment is quite diverse. For example, it has been argued that
is the radial excitation of
with
, but the predicted mass is too large by the order of 40 MeV. Alternatively, it could be the first radial excitation of
(not
!) with
[
26]. The latter assignment has the advantage that the predicted mass is in better agreement with experiment. Because the mass of
is barely below the threshold of
, this observation motivated the authors of Ref. [
27] to suggest an exotic molecular state of
and
p with a binding energy of the order of 6 MeV and
for
. The quark potential model predicts a
state at 2900 MeV and a
state at 2910 MeV [
17,
18]. A similar result of 2906 MeV for
was also obtained in the relativistic quark model [
28].
2.1.2 states
The highest isotriplet charmed baryons,
, which decay to
, were first measured by Belle [
29] with widths of the order of 70 MeV. The possible quark states are
,
,
, and
. The states
and
are ruled out because their decays to
are prohibited in the heavy quark limit. Now the
baryon decays primarily into the
system in a
D-wave, whereas
decays into
in an
S-wave. Because heavy hadron chiral perturbation theory (HHChPT) implies a very broad
with a width of the order of 885 MeV (see Section 3.2 below), this
p-wave state is also excluded. Therefore,
are likely to be either
or
or a mixture of the two. In the quark-diquark model [
26], both of them have very close masses compatible with experiment. Given that for light strange baryons, the first orbital excitation of the
Σ also has the quantum numbers
(see Fig. 2), we will advocate a
state for
.
2.1.3 states
The states and form the doublet . Because the diquark transition is prohibited, cannot decay to . The dominant decay modes are for and for .
Many excited charmed baryon states
,
,
, and
have been seen at
B factories [
12,
30,
31]. Another state
, which is omitted from the PDG summary table, has been seen only by BaBar in the
mass projection of
[
33]. However, as we shall see below, it may form a sextet with
and
. The states
,
,
, and
could be the first positive-parity excitations of the
. The study of Regge phenomenology is very useful for the
JP assignment of charmed baryons [
26,
34]. The Regge analysis suggests that
for
and
for
[
26]. From Table 5 below, we shall see that
and
form a nice
antitriplet.
In the relativistic quark–diquark model [
26],
is a sextet
state. According to Table 2, possible candidates are
, and
. As pointed out in Ref. [
32], strong decays of these four states, which were studied in Ref. [
22] using the
model show that
does not decay to
and
, and has a width of 28 MeV, consistent with experiment. Therefore, the favored candidate for
is
, which has
.
The possible quark states for the
baryon in an antitriplet are
,
,
,
, and
(see Table 2). Because
is above the
threshold, the two-body mode
should exist although it has not been searched for in the
spectrum. Recall that the neutral
was observed recently by Belle in the
spectrum [
13]. According to the
model, the first four states are excluded as they do not decay into
[
22]. The only remaining possibility is
. This is the analog of
for
. Nevertheless, the identification of
with
encounters two potential problems: (i) its width is dominated by the
and
modes, which have not been seen experimentally, and (ii) the predicted width of the order of 47 MeV [
22] is too large compared to the measured one, which is of the order of 5.7 MeV.
2.1.4 states
Only two ground states have been observed thus far:
and
. The latter was seen by BaBar in the electromagnetic decay
[
39].
Molecular picture
Because
and
are barely below the
and
thresholds, respectively, it is tempting to conjecture an exotic molecular structure of
and
p for the former and a molecular state of
for the latter [
27,
40–
44]. Likewise,
could be a molecular state of
.
A coupled-channel calculation of the baryon–meson
system has been performed to look for an isospin–spin channel that is attractive enough to form a molecular state [
42,
45].
was found to be the most attractive one, followed by
. This suggests that
might be an
s-wave
molecular state with
and
an
s-wave
molecular state with
(see Fig. 3 of [
45]). Another possibility is a
molecular state with
for
and a
state with
for
. Because
has isospin 1, and moreover, as we have noted in passing,
will be too broad if it is assigned to
, we conclude that the second possibility is preferable (see also [
44]).
The possible spin-parity quantum numbers of the higher excited charmed baryon resonances that have been suggested in the literature are partially summarized in Table 4. Some of the predictions are already ruled out by experiment. For example, has , as seen by Belle. More experimental studies are certainly to pin down the quantum numbers.
Charmed baryon spectroscopy has been studied extensively in various models. It appears that the spectroscopy is well described by the heavy quark–light diquark picture elaborated by Ebert, Faustov, and Galkin (EFG) [
26] (see also Ref. [
35]). As noted in passing, the quantum numbers
of
were correctly predicted in a model based on the diquark idea before the Belle experiment [
24,
25]. Moreover, EFG have shown that all the available experimental data on heavy baryons fit the linear Regge trajectories nicely, namely, the trajectories in the
and
planes for orbitally and radially excited heavy baryons, respectively,
where
nr is the radial excitation quantum number,
α and
β are the slopes, and
and
are intercepts. The linearity, parallelism, and equidistance of the Regge trajectories were verified. The predictions of the spin-parity quantum numbers of charmed baryons and their masses in Ref. [
26] can be regarded as a theoretical benchmark (see Fig. 3).
Antitriplet and sextet states
The antitriplet and sextet states of charmed baryons are listed in Table 5. To date, the following states are established: the , , and states, (, , (, , and (, , respectively, and the and states, () and (), respectively. The mass difference in the antitriplet states clearly lies between 180 and 200 MeV. We note in passing that should carry the quantum numbers . From Table 5, we see that and form a nice antitriplet as the mass difference between and is consistent with that observed in other antitriplets. Likewise, the mass differences in the sextet (, predicted by the quark–diquark model are consistent with that measured in and sextets. Note that there is no sextet as the with these spin-parity quantum numbers will be too broad to be observed.
On the basis of the QCD sum rules, many charmed baryon multiplets classified according to
were recently studied in Ref. [
36]. Three sextets were proposed in this work:
for
and
for
. Note that
and
were treated as
p-wave baryons rather than the first positive-parity excitations, as we have discussed before. The results for the multiplet
led the authors of Ref. [
36] to suggest that there are two
,
, and
states with
and
. The mass splittings are
,
, and
MeV, respectively. The predicted mass of
is around
MeV. Using the central value of the predicted masses to label the states in the multiplet
(see Table I of Ref. [
36]), one obtains
for
and
for
. One can check that
MeV, and
is of the order of
MeV. Owing to the large theoretical uncertainties in the masses, it is not clear whether the QCD sum rule calculations are compatible with the mass differences measured in the
and
sextets. In any event, it will be interesting to test these two different model predictions for
and
sextets in the future.
2.2 Doubly charmed baryons
Evidence of doubly charmed baryon states has been reported by SELEX in
[
46]. Further observation of
was also announced by SELEX [
47]. However, none of the doubly charmed states discovered by SELEX has been confirmed by FOCUS [
48], BaBar [
49], Belle [
12] and LHCb [
50], although
events are produced in
B factories, for example, versus 1630
events observed at SELEX.
The doubly charmed baryons
with the quark contents
form an SU(3) triplet. They have been studied extensively using many different approaches: the quark model, light quark–heavy diquark model, QCD sum rules, and lattice simulation. The predicted doubly charmed baryon masses calculated in various models are tabulated in Refs. [
51,
52]. For recent QCD sum rule calculations, see e.g. [
53–
55]. Chiral corrections to the masses of doubly heavy baryons up to N
3LO were presented in Ref. [
56].
Figure 4 shows the results of recent lattice studies of doubly and triply charmed baryon spectra by different groups: RQCD [
57], HSC [
58], Brown
et al. [
59], ETMC [
60], ILGTI [
61], PACS-CS [
62], Durr
et al. [
63], Briceno
et al. [
64], Liu
et al. [
65], and Na
et al. [
66]. A new lattice calculation of
and
was available in Ref. [
67]. The various lattice results are consistent with each other and they fall into the ranges
and
Although lattice studies suggest that the mass of the low-lying
exceeds 3519 MeV, it is interesting to note that the authors of [
51] calculated the masses of doubly and triply charmed baryons on the basis of the Regge phenomenology and found
MeV, in good agreement with SELEX.
3 Strong decays
Owing to the rich mass spectrum and relatively narrow widths of the excited states, the charmed baryon system offers an excellent ground for testing the ideas and predictions of heavy quark symmetry and light flavor SU(3) symmetry. The pseudoscalar mesons involved in the strong decays of charmed baryons such as are soft. Therefore, heavy quark symmetry of the heavy quark and chiral symmetry of the light quarks will have interesting implications for the low-energy dynamics of heavy baryons interacting with Goldstone bosons.
The strong decays of charmed baryons are most conveniently described by the HHChPT, into which heavy quark symmetry and chiral symmetry are incorporated [
68–
70]. Heavy baryon chiral Lagrangians were first constructed in Ref. [
68] for strong decays of
s-wave charmed baryons and in Refs. [
9,
14] for
p-wave ones. Previous phenomenological studies of the strong decays of
p-wave charmed baryons based on HHChPT can be found in Refs. [
9,
10,
14,
71,
72]. The chiral Lagrangian involves two coupling constants,
g1 and
g2, for
P-wave transitions between
s-wave and
s-wave baryons [
68]; six couplings,
, for the
S-wave transitions between
s-wave and
p-wave baryons; and eight couplings;
, for the
D-wave transitions between
s-wave and
p-wave baryons [
9]. The general chiral Lagrangian for heavy baryons coupled to pseudoscalar mesons can be expressed compactly in terms of superfields. We will not write the relevant Lagrangians here; instead the reader is referred to Eqs. (3.1) and (3.3) of Ref. [
9]. The partial widths relevant for our purposes are [
9]
where MeV. The dependence on the pion momentum is proportional to , and for S-wave, P-wave and D-wave transitions, respectively. It is obvious that the couplings are dimensionless, whereas have canonical dimension .
3.1 Strong decays of s-wave charmed baryons
Because the strong decay
is kinematically prohibited, the coupling
g1 cannot be extracted directly from the strong decays of heavy baryons. In the framework of HHChPT, one can use some measurements as input to fix the coupling
g2, which, in turn, can be used to predict the rates of other strong decays. Among the strong decays
,
is the most well-measured. Hence, we shall use this mode to extract the coupling
g2. Using the 2006 data [
73]
, the coupling
g2 is extracted as
The predicted rates of other modes are shown in Table 6. The agreement between theory and experiment is clearly excellent, except that the predicted width for is slightly too large.
Using the new data from the 2014 Particle Data Group [
8] in conjunction with the new measurements of the
and
widths by Belle [
11], we obtain the new average
(see Table 3). Therefore, the coupling
g2 is reduced to
From Table 7 we see that the agreement between theory and experiment is further improved: The predicted
width is consistent with the first new measurement by Belle [
12], and the new calculated width for
is now in agreement with experiment. It is also clear that the
width is smaller than that of
by a factor of
, although they will become the same in the limit of heavy quark symmetry.
3.2 Strong decays of p-wave charmed baryons
Because
and
form the doublet
, it appears from Eq. (3.1) that the couplings
h2 and
h8 can in principle be extracted from
and from
, respectively. Likewise, the information on the couplings
and
can be inferred from the strong decays of
identified with
. Couplings other than
h2,
h8, and
can be related to each other via the quark model [
9].
Although the coupling
h2 can be inferred from the two-body decay
, this method is less accurate because this decay is kinematically barely allowed or even prohibited, depending on the mass of
. For the old mass measurement
MeV [
73],
and
are kinematically barely allowed. However, for the new measurement by the CDF,
MeV [
79], only the last mode is allowed. Moreover, the finite width effect of the intermediate resonant states could become important [
71].
We next turn to the three-body decays
of
and
to extract
h2 and
h8. As shown in Ref. [
10], the 2006 data
MeV [
73] and for the
mass lead to the resonant rate [
10]
as shown in Table 9. Assuming the pole contributions to
due to the intermediate states
and
, the resonant rate for the process
can be calculated in the framework of HHChPT [
9]. Numerically, we found
where
. It is clear that the limit on
gives an upper bound on
h8 of the order of
(in units of MeV
-1), whereas the decay width of
is entirely governed by the coupling
h2. Specifically, we have [
10]
It has been noted [
72] that the proximity of the
mass to the sum of the masses of its decay products will lead to an important threshold effect that will lower the
mass by
MeV compared to the observed mass. A more sophisticated treatment of the mass lineshape of
by the CDF yields
MeV [
79], which is 3.1 MeV smaller than the 2006 world average. Therefore, the strong decay
is very close to the threshold. With the new measurement of
, we have (in units of MeV) [
32]
By fitting to the measured
mass difference distributions and using
, CDF obtained
or
[
79]. This corresponds to a decay width
MeV [
79]. For the width of
, CDF observed a value consistent with zero and therefore calculated an upper limit of 0.97 MeV using a Bayesian approach. From the CDF measurements
MeV and
MeV, we obtain
Hence, the magnitude of the coupling h2 is greatly enhanced from 0.437 to 0.63 . Our h2 differs slightly from the value of 0.60 obtained by CDF. This is because the CDF used to calculate the mass dependence of , whereas we used .
The large difference between the values of the coupling h2 obtained in 2006 and 2015 is ascribed to the fact that the mass of is 3.1 MeV lower than the previous world average because of the threshold effect. To illustrate this, we consider the dependence of and on as depicted in Fig. 5. at MeV is clearly smaller than that at 308.9 MeV. This explains why h2 should become larger when becomes smaller.
The and baryons form the doublet . The and widths predicted using the coupling h2 obtained from Eq. (3.8) and assuming SU(3) flavor symmetry are shown in Table 9. The predicted two-body decay rates of and clearly exceed the current experimental limits because of the enhancement of h2 (see Table 9). Hence, there is a tension for the coupling h2, as its value extracted from from will imply and rates slightly above current limits. It is conceivable that SU(3) flavor symmetry breaking can help account for the discrepancy.
Some information on the coupling
can be inferred from the strong decays of
. From Eq. (3.1) and the quark model relation
from Ref. [
9], we obtain, for example,
MeV. Hence,
cannot be identified with
. Using the quark model relation
and the measured widths of
(Table 3), we obtain
The quark model relation then leads to
which improves the previous limit (3.8) by a factor of 3. The calculated partial widths of shown in Table 9 are consistent with experimental limits.
4 Lifetimes
4.1 Singly charmed baryons
Among singly charmed baryons, the antitriplet states
, and the
baryon in the sextet decay weakly. In 2006, the world averages of their lifetimes were [
73]
These results remain the same even in 2014 [
8]. As we shall see below, the lifetime hierarchy
is understood qualitatively but not quantitatively in the operator product expansion (OPE) approach.
On the basis of the OPE approach to analysis of inclusive weak decays, the inclusive rate of the charmed baryon is schematically represented by
where
is the relevant Cabibbo-Kobayashi-Maskawa matrix element. The
A0 term comes from the
c quark decay and is common to all charmed hadrons. There are no linear
corrections to the inclusive decay rate owing to the lack of gauge-invariant dimension-four operators [
81–
84], a consequence known as Luke’s theorem [
85]. Nonperturbative corrections start at order
and are model-independent. Spectator effects in inclusive decays due to the Pauli interference and
W-exchange contributions account for the
corrections and they have two noteworthy features: First, the estimate of the spectator effects is model dependent; the hadronic four-quark matrix elements are usually evaluated by assuming the factorization approximation for mesons and the quark model for baryons. Second, there is a two-body phase-space enhancement factor of
for spectator effects relative to the three-body phase space for heavy quark decay. This implies that spectator effects, which are of the order of
, are comparable to and even exceed the
terms.
The total width of the charmed baryon generally receives contributions from inclusive nonleptonic and semileptonic decays: . The nonleptonic contribution can be decomposed into
corresponding to the c-quark decay, W-exchange contribution, and destructive and constructive Pauli interference. The inclusive decay rate is known to be governed by the imaginary part of an effective nonlocal forward transition operator T. Therefore, corresponds to the imaginary part of Fig. 6(a) sandwiched between the same states. At the Cabibbo-allowed level, represents the decay rate of , and denotes the contribution from the W-exchange diagram . The interference () arises from destructive (constructive) interference between the u (s) quark produced in c-quark decay and the spectator u (s) quark in the charmed baryon . Note that the constructive Pauli interference is unique to the charmed baryon sector, as it does not occur in the bottom sector. From the quark content of the charmed baryons, it is clear that at the Cabibbo-allowed level, destructive interference occurs in and decays (Fig. 6(c)), whereas and can have constructive interference (Fig. 6(d)). Because contains two s quarks, it is natural to expect that . The W-exchange contribution (Fig. 6(b)) occurs only for and at the same Cabibbo-allowed level. In the heavy quark expansion approach, the above-mentioned spectator effects can be described in terms of the matrix elements of local four-quark operators.
The inclusive nonleptonic rates of charmed baryons in the valence quark approximation and in the limit
are expressed as [
86]:
where is the Cabibbo angle.
The results of a model calculation in Ref. [
86] are shown in Table 10. The lifetime pattern
clearly agrees with experiment. This lifetime hierarchy is qualitatively understandable. The baryon is the longest-lived among charmed baryons because of the smallness of W-exchange and partial cancellation between constructive and destructive Pauli interference, whereas is the shortest-lived owing to the presence of two s quarks in the , which greatly enhances the contribution of . Because is always positive, is negative, and the constructive interference has a larger magnitude than the destructive interference, this explains why . It is also clear from Table 10 that, although the qualitative feature of the lifetime pattern is comprehensive, the quantitative estimates of charmed baryon lifetimes and their ratios are still rather poor.
4.2 Doubly charmed baryons
The inclusive nonleptonic rates of doubly charmed baryons in the valence quark approximation and in the limit are expressed as
Because is positive and is negative, it is obvious that is longest-lived, whereas () is the shortest-lived if (). In general, we have
The predictions available in the literature are summarized in Table 11. Note that the lifetime of
was measured by SELEX to be
[
46].
Because the mass splitting between and and between and is less than 100 MeV (see also Eq. (2.6) for the lattice calculations),
it is clear that only electromagnetic decays are allowed for and .
5 Hadronic weak decays
5.1 Nonleptonic decays
In contrast to the significant advances made over the last 10 years or so in the study of hadronic weak decays in the bottom baryon sector, progress in the arena of charmed baryons, both theoretical and experimental, has been very slow.
In the naive factorization approach, the coefficients
a1 for the external
W-emission amplitude and
a2 for internal
W-emission are given by
and
, respectively. However, we have learned from charmed meson decays that the naive factorization approach never works for the decay rate of color-suppressed decay modes, although it usually works for color-allowed decays. Empirically, it was learned in the 1980s that if the Fierz-transformed terms characterized by
are dropped, the discrepancy between theory and experiment is greatly reduced [
91–
93]. This leads to the so-called large-
Nc approach to describing hadronic
D decays [
94]. Theoretically, explicit calculations based on QCD sum-rule analysis [
95–
97] indicate that the Fierz terms are indeed largely compensated by the nonfactorizable corrections.
As the
expansion method greatly reduces the discrepancy between theory and experiment for charmed meson decays, it is natural to ask if this scenario also works in the baryon sector. This issue can be settled by experimental measurement of the Cabibbo-suppressed mode
, which receives contributions only from the factorizable diagrams. As pointed out in Ref. [
98], the rate predicted by the large-
Nc approach is in good agreement with the measured value. In contrast, its decay rate predicted by the naive factorization approximation will be too small by a factor of 15. Therefore, the
approach also works for the factorizable amplitude of charmed baryon decays. This also implies that the inclusion of nonfactorizable contributions is inevitable and necessary. If nonfactorizable effects amount to a redefinition of the effective parameters
a1,
a2, and are universal (i.e., channel-independent) in charm decays, then we still have a new factorization scheme with the universal parameters
to be determined from experiment.
It is known that for heavy mesons, the nonfactorizable contributions will render color suppression of internal
-emission ineffective. However,
-exchange in baryon decays is not subject to color suppression even in the absence of nonfactorizable terms. A simple way to see this is to consider the large-
limit. Although the
-exchange diagram is down by a factor of
relative to the external
-emission one, this difference is compensated by the fact that the baryon contains
quarks in the limit of large
, thus allowing
different possibilities for
exchange between heavy and light quarks [
100]. That is, the pole contribution can be as important as the factorizable one. Experimental measurement of the decay modes
, which proceed only through the
-exchange contributions, indicates that
W-exchange indeed plays an essential role in charmed baryon decays.
Various theoretical approaches to weak decays of heavy baryons have been investigated, including the current algebra approach [
101–
117], factorization scheme, pole model technique [
98,
118–
122], relativistic quark model [
100,
123], and quark diagram scheme [
124,
125]. Model predictions of the branching fractions and decay asymmetries can be found in Tables VI–VII of [
1] for
decays, Table VIII for
decays and Table IX for
decays.
5.2 Discussion
5.2.1 decays
Experimentally, nearly all the branching fractions of the
are measured relative to the
mode. On the basis of ARGUS and CLEO data, the PDG has made a model-dependent determination of the absolute branching fraction,
[
8]. Recently, Belle reported a value of
[
126] from the reconstruction of
recoiling against the
production in
annihilation. Hence, the uncertainties are much smaller, and, most importantly, this measurement is model independent! More recently, BESIII has also measured this mode directly with the preliminary result
[
127]. Its precision is comparable to that of Belle’s result. Another approach is to exploit a particular decay,
, and its charge conjugate to measure
, also in a model independent manner [
128].
The branching fractions of the Cabibbo-allowed two-body decays of
are listed in Table 12. The central values of data taken from the PDG [
8] are scaled up by a factor of 1.37 because of the new measurement of
by Belle [
126]. BESIII recently measured 2-body, 3-body, and 4-body decay modes of
with significantly improved precision [
127]. For example, the result
obtained by BESIII has much better precision than the value of
quoted by the PDG [
8].
Many of the decay modes such as , , , and can proceed only through W-exchange. Experimental measurement of them implies the importance of W-exchange, which is not subject to color suppression in charmed baryon decays.
Some Cabibbo-suppressed modes such as
and
have been measured by Belle [
129] and BaBar [
130], respectively. Their branching fractions are of the order of
. The first measured Cabibbo-suppressed mode
, is of particular interest because it receives contributions only from the factorizable diagram and is expected to be color suppressed in the naive factorization approach. A calculation in Refs. [
131,
132] yields
From the experimental measurement
[
8]
†, it follows that
This is consistent with the result obtained by the approach, .
All the models except that of Ref. [
121] predict a positive decay asymmetry
α for the decay
(see Table VII of Ref. [
1]). Therefore, the measurement of
by CLEO [
133] is very surprising. If the negative sign of
α is confirmed in the future, this will imply that the
s-wave and
p-wave amplitudes for this decay have opposite signs, contrary to the model expectation. The implication of this has been discussed in detail in Refs. [
98,
99]. Because the error of the previous CLEO measurement is very large, it is crucial to more accurately measure the decay asymmetry for
.
5.2.2 decays
No absolute branching fractions have been measured. The branching ratios listed in Tables VI and VIII of Ref. [
1] are those relative to
. Several Cabibbo-suppressed decay modes such as
,
,
,
, and
have been observed [
8].
The Cabibbo-allowed decays
have been studied and are believed to be forbidden, as they do not receive factorizable and
pole contributions [
100,
120]. However, the
mode was seen earlier by FOCUS [
134], and this may indicate the importance of pole contributions beyond low-lying
intermediate states.
5.2.3 decays
No absolute branching fractions have been measured to date. However, there are several measurements of the ratios of branching fractions, for example [
8],
The decay modes
and
and
proceed only through
W-exchange. The measured branching ratio of
relative to
implies that the
W-exchange mechanism plays a significant role. The model of Körner and Krämer [
100] predicts
(see Table IX of Ref. [
1]), in agreement with experiment, but its prediction
is too small compared to the data.
5.2.4 decays
A unique feature of
decays is that the decay
proceeds only via external
W-emission, whereas
proceeds via the factorizable internal
W-emission diagram. Various model predictions of Cabibbo-allowed
are listed in Table IX of [
1] with the unknown parameters
a1 and
a2. From the decay
we learn that
. The hadronic weak decays of the
were recently studied in great detail in Ref. [
135], where most of the decay channels in
decays were found to proceed only through the
W-exchange diagram; moreover, the
W-exchange contributions dominate in the rest of the processes, with some exceptions. Observation of such decays would shed light on the mechanism of
W-exchange effects in these decay modes.
5.3 Charm-flavor-conserving nonleptonic decays
There is a special class of weak decays of charmed baryons that can be studied reliably, namely, heavy-flavor-conserving nonleptonic decays. Some examples are the singly Cabibbo-suppressed decays
and
. The idea is simple: In these decays, only the light quarks inside the heavy baryon will participate in weak interactions; that is, while the two light quarks undergo weak transitions, the heavy quark behaves as a “spectator”. As the emitted light mesons are soft, the
weak interactions among light quarks can be handled by the well known short-distance effective Hamiltonian. This special class of weak decays can usually be calculated more reliably than the conventional charmed baryon weak decays. The synthesis of the heavy quark and chiral symmetries provides a natural setting for investigating these reactions [
136]. The weak decays
with
were also studied in Refs. [
137–
139].
The combined symmetries of heavy and light quarks severely restrict the weak interactions allowed. In the symmetry limit, it is found that
and
nonleptonic weak transitions [
136] cannot occur. Symmetries alone permit three types of transitions:
,
and
transitions. However, in both the MIT bag and diquark models, only
transitions have nonzero amplitudes. The general amplitude for
is given by
where A and B are the S- and P-wave amplitudes, respectively. The S-wave amplitude can be evaluated using current algebra in terms of the parity-violating commutator term. For example, the S-wave amplitude of is given by
whereas the
P-wave amplitude arises from the ground-state baryon poles [
136]:
where
φ is the mixing angle of
, and
and
,
are their mass eigenstates. The matrix element
was evaluated in Ref. [
136] using two different models: the MIT bag model [
140,
141] and the diquark model.
The predicted rates are [
136]
and the corresponding branching fractions are
As stated above, the transition vanishes in the chiral limit. It receives a finite factorizable contribution as a result of a symmetry-breaking effect. At any rate, the predicted branching fractions for the charm-flavor-conserving decays and are of the order of and should be readily accessible in the near future.
5.4 Semileptonic decays
Exclusive semileptonic decays of charmed baryons: , and have been observed experimentally. Their rates depend on the form factors and (), which are defined as
These form factors have been evaluated using the non-relativistic quark model [
131,
132,
142–
144,
151], MIT bag model [
142,
143], relativistic quark model [
145,
149], light-front quark model [
146], and QCD sum rules [
147,
148,
150]. Experimentally, the only information available to date is the form-factor ratio measured in the semileptonic decay
. In the heavy quark limit, the six
form factors are reduced to two:
Assuming that the form factors exhibit dipole
q2 behavior, the ratio
is measured by CLEO to be [
152]
Various model predictions of the charmed baryon semileptonic decay rates and decay asymmetries are shown in Table 13. Dipole
q2 dependence of the form factors is assumed whenever the form factor momentum dependence is not available in the model. Four different sets of predictions for
, which are not listed in Table 13, were presented in the sum rule calculations of Ref. [
153]. The semileptonic decays of
were treated in Ref. [
154] within the framework of a constituent quark model. From Table 13, we see that the computed branching fractions of
, which fall in the range
are slightly smaller than the experimental values,
[
, of the PDG [
8]]. The branching fractions of
and
are predicted to lie in the ranges
and
, respectively, except that the QCD sum rule calculation in Ref. [
150] predicts a much larger rate for
. Experimentally, only the ratios of the branching fractions are available to date [
8]:
There have been active studies of semileptonic decays of doubly charmed baryons. The interested reader can consult [
155–
159] for further references.
Just as with the hadronic decays discussed in the last subsection, there are also heavy-flavor-conserving semileptonic processes, for example,
and
. In these decays only the light quarks inside the heavy baryon will participate in weak interactions, while the heavy quark behaves as a spectator. This topic was recently investigated in Ref. [
139]. Owing to the severe phase-space suppression, the branching fractions are of order
in the best cases, and typically
to
.
6 Electromagnetic and weak radiative decays
Although radiative decays are well measured in the charmed meson sector, e.g., and , only three of the radiative modes in the charmed baryon sector have been seen, namely, , and . This is understandable because MeV, and MeV. Hence, and are governed by the electromagnetic decays. However, it will be difficult to measure the rates of these decays because these states are too narrow to be experimentally resolvable. Nevertheless, we shall systematically study the two-body electromagnetic decays of charmed baryons and also weak radiative decays.
6.1 Electromagnetic decays
In the charmed baryon sector, the following two-body electromagnetic decays are of interest:
where we denote the spin baryons as and for the symmetric sextet 6 and antisymmetric antitriplet 3̄, respectively, and the spin baryon by .
An ideal theoretical framework for studying the above-mentioned electromagnetic decays is provided by the formalism in which the heavy quark symmetry and the chiral symmetry of light quarks are combined [
68–
70]. When supplemented by the nonrelativistic quark model, the formalism determines completely the low energy dynamics of heavy hadrons. The electromagnetic interactions of heavy hadrons consist of two distinct contributions: one from gauging electromagnetically the chirally invariant strong interaction Lagrangians for heavy mesons and baryons given in Refs. [
68–
70], and the other from the anomalous magnetic moment couplings of the heavy particles. The heavy quark symmetry reduces the number of free parameters needed to describe the magnetic couplings to the photon. There are two undetermined parameters for the ground-state heavy baryons. All of these parameters are related simply to the magnetic moments of the light quarks in the nonrelativistic quark model. However, the charmed quark is not particularly heavy (
GeV), and it carries a charge of
. Consequently, the contribution from its magnetic moment cannot be neglected. The chiral and electromagnetic gauge-invariant Lagrangian for heavy baryons can be found in Eqs. (3.8) and (3.9) of Ref. [
160], and are denoted by
and
, respectively.
The most general gauge-invariant Lagrangian Eq. (3.9) of Ref. [
160] for magnetic transitions of heavy baryons can be recast in terms of superfields [
161,
162]:
where , is the charge matrix for the light and s quarks, and is the charge of the heavy quark. In the above equation,
It follows that [
161,
162]
where
is the photon 4-momentum, and
is the polarization 4-vector. As stressed in Ref. [
160], SU(3) breaking effects due to light-quark mass differences can be incorporated by replacing the charge matrix
Q with
where and . To avoid any confusion with the current quark mass mq, we have used capital letters to denote the constituent quark masses. In the quark model, the coefficients a1 and a2 are simply related to the Dirac magnetic moments of light quarks,
whereas is connected to those of heavy quarks. Explicitly, is fixed by heavy quark symmetry as
Within the framework of HHChPT, the authors of [
163] proceeded to construct chiral Lagrangians at the level
and
and then calculated the electromagnetic decay amplitudes of charmed baryons up to
. It is not clear whether their
Lagrangian (see Eq. (12) of Ref. [
163]), characterized by the four couplings
and
f4, is equivalent to the first two terms of the
Lagrangian given by Eq. (6.2). The unknown couplings there were also estimated using the quark model.
The general amplitudes of electromagnetic decays are given by [
160]
The corresponding decay rates are [
160]
where
mi (
mf) is the mass of the parent (daughter) baryon. The coupling constants
can be calculated using the quark model for
a1,
a2, and
[
160,
164]:
The results calculated using the constituent quark masses,
MeV,
MeV,
MeV [
8], and
GeV, are summarized in the second column of Table 14. Some other model predictions are also listed there for comparison.
Radiative decays of
s-wave charmed baryons are considered in Ref. [
165] in the quark model, and the predictions are similar to ours. A similar procedure is followed in Ref. [
166] where the heavy quark symmetry is supplemented with light-diquark symmetries to calculate the widths of
and
. The authors of Ref. [
75] apply the relativistic quark model to predict various electromagnetic decays of charmed baryons. In addition to the magnetic dipole (M1) transition, the author of Ref. [
167] also considered and estimated the electric quadrupole (E2) amplitude for
arising from the chiral loop correction. The E2 contributions were analyzed in detail in Ref. [
168]. The E2 amplitudes appear at different higher orders for the three types of decays:
for
,
for
and
for
. Therefore, the E2 contribution to
is completely negligible. The electromagnetic decays were calculated in Refs. [
169–
171] using the QCD sum rule method, and they were studied within the framework of the modified bag model in Ref. [
172].
It is evident from Table 14 that the predictions in Refs. [
160,
164] and [
163] all based on HHChPT are quite different for the following three modes:
,
and
. Indeed, the results for the last two modes in Ref. [
163] are larger than all the other existing predictions by one order of magnitude! It is naively expected that all HHChPT approaches should agree with each other to the lowest order of chiral expansion provided that the coefficients are inferred from the nonrelativistic quark model. The lowest order predictions
keV and
keV obtained in Ref. [
163] are still very large. Note that a recent lattice calculation in Ref. [
174] yields
keV which is much smaller than the value of
keV predicted in Ref. [
163].
Chiral loop corrections to the M1 electromagnetic decays and to the strong decays of heavy baryons were computed at the one loop order in Refs. [
161,
162]. The leading chiral-loop effects we found are nonanalytic in the forms of
and
(or
and
, where
mq is the light quark mass). Some results are [
161,
162]
which should be compared with the corresponding quark-model results: 92, 20, and 0.4 keV, respectively (Table 14).
The electromagnetic decays
and
are of special interest. It was advocated in Ref. [
175] that a measurement of their branching fractions will allow us to determine one of the coupling constants in HHChPT, namely,
g1. They are forbidden at tree level in the SU(3) limit [see Eq. (6.10)]. In heavy baryon chiral perturbation theory, this radiative decay is induced via chiral loops where SU(3) symmetry is broken by the light current quark masses. By identifying the chiral loop contribution to
with the quark model prediction given in Eq. (6.10), it was found in Ref. [
164] that one of the two possible solutions agrees with the quark model expectation for
g1.
For the electromagnetic decays of
p-wave charmed baryons, the search for
and
has not yet succeeded. The interested reader is referred to Refs. [
14,
75,
166,
175–
180] for more details on the theoretical side. Some predictions are collected in Table 15, and they are more diversified than the
s-wave case. For electromagnetic decays of doubly charmed baryons, see, e.g., Refs. [
172,
181].
The electromagnetic decays considered to date do not critically test the heavy quark symmetry or chiral symmetry. The results follow simply from the quark model. There are examples in which both the heavy quark symmetry and chiral symmetry enter in a crucial way. These are the radiative decays of heavy baryons involving an emitted pion. Some examples that are kinematically allowed are
The contact interaction dictated by the Lagrangian
can be nicely tested by the decay
, whereas a test of the chiral structure of
is provided by the process
; see Ref. [
160] for the analysis.
6.2 Weak radiative decays
At the quark level, three different types of processes can contribute to the weak radiative decays of heavy hadrons, namely, single-, two-, and three-quark transitions [
183]. The single-quark transition mechanism arises from the so-called electromagnetic penguin diagram. Unfortunately, the penguin process
is highly suppressed, so it plays no role in radiative decays of charmed hadrons. There are two contributions from the two-quark transitions: one from the
W-exchange diagram accompanied by photon emission from the external quark, and the other from the same
W-exchange diagram but with a photon radiated from the
W boson. The latter is typically suppressed by a factor of
(where
k is the photon energy) compared to the former bremsstrahlung process [
182]. For charmed baryons, the Cabibbo-allowed decay modes via
(Fig. 7) or
are
Finally, the three-quark transition involving
W-exchange between two quarks and photon emission by the third quark is quite strongly suppressed because of the very small probability of finding three quarks that are adequately kinematically matched with the baryons [
183,
184].
The general amplitude of weak radiative decays of baryons is given by
where a and b are the parity-conserving and -violating amplitudes, respectively. The corresponding decay rate is
Nonpenguin weak radiative decays of charmed baryons, such as those in Eq. (6.13), are characterized by emission of a hard photon and the presence of a highly virtual intermediate quark between the electromagnetic and weak vertices. It was shown in Ref. [
185] that these features should make it possible to analyze these processes by perturbative QCD; that is, these processes can be described by an effective local and gauge invariant Lagrangian:
with
where , , and
For radiative decays of charmed baryons, one needs to evaluate the matrix element
. Because the quark-model wave functions most resemble the hadronic states in the frame where both baryons are static, the static MIT bag model was adopted in Ref. [
185] for the calculation. The predictions are
†A different analysis of the same decays in Ref. [
186] yielded the following results:
Evidently, these predictions (especially the decay asymmetry) differ greatly from those obtained in Ref. [
185].
Finally, it is worth remarking that, analogous to the heavy-flavor-conserving nonleptonic weak decays discussed in Section 6.3, there is a special class of weak radiative decays in which heavy flavor is conserved, for example, and . In these decays, weak radiative transitions arise from the light quark sector of the heavy baryon, whereas the heavy quark behaves as a spectator. However, the dynamics of these radiative decays is more complicated than that of their counterpart in nonleptonic weak decays, e.g., . In any event, it merits an investigation.
7 Conclusions
In this report, we began with a brief overview of the spectroscopy of charmed baryons and discussed their possible structure and spin-parity assignments in the quark model. For p-wave baryons, we have assigned to . As for first positive-parity excitations, with the help of the relativistic quark-diquark model and the model, we have identified with , with , and with , although the first and last assignments may encounter some potential problems.
It should be stressed that mass analysis alone is usually not sufficient to pin down the spin-parity quantum numbers of higher excited charmed baryon states, a study of their strong decays is necessary. For example, , , and for all have similar masses. Analysis of strong decays allows us to exclude the first two possibilities. It should be stressed that Regge phenomenology and the mass relations for antitriplet and sextet multiplets also provide very useful guidance for the spin-parity quantum numbers.
On the basis of various theoretical tools such as lattice QCD and the QCD sum rule method, there are many theoretical studies of charmed baryon spectroscopy, especially for doubly and triply charmed baryons. However, progress on the hadronic decays, radiative decays, and lifetimes has been very slow. Experimentally, nearly all the branching fractions of the are measured relative to the mode. The recent measurements by Belle and by BESIII are very encouraging. Moreover, BESIII recently measured 2-body, 3-body and 4-body nonleptonic decay modes of with significantly improved precision. It is conceivable that many new data emerging from LHCb, BESIII, and B factories in the near future and from the experiments at J-PARC and ANDA in the future can be used to test the underlying mechanism of hadronic weak decays.