Bottom-charmed baryons in a nonrelativistic quark model

In this work, we study the low-lying mass spectra for bottom-charmed baryons in a nonrelativistic quark model by solving the three-body Schr¨odinger equation. The lowest Ξ bc , Ξ ′ bc , Ω bc , and Ω ′ bc states are predicted to be about 6979, 6953, 7109, and 7092 MeV, respectively. Also, the strong decays for the low-lying excited states are investigated. Our results indicate that some of λ − mode P − wave bottom-charmed baryons are relatively narrow, which can be searched for in future experiments. For the low-lying ρ − mode and ρ − λ hybrid states, their strong decays are highly suppressed and they can survive as extremely narrow states. Moreover, the mass spectra and strong decays for bottom-charmed baryons preserve the heavy quark symmetry well. We hope our calculations can provide helpful information for further experimental and theoretical researches.


I. INTRODUCTION
In the past years, with the development of the large-scale accelerator facilities, plenty of heavy baryons have been observed and significant progress has been achieved in experiments.These discoveries have triggered wide attentions of theorists, which leads the study on mass spectra and internal structures of heavy baryons to a hot topic in hadron physics [1,2].Understanding the nature of heavy baryons and searching for the missing heavy resonances can help us to establish and complete the hadron spectroscopy and provide a good platform to better investigate the heavy quark symmetry.Until now, most experimental observations in heavy baryon sector belong to the singly heavy baryons, while the doubly heavy baryons are rare.Due to the lack of experimental data, our understanding of the doubly heavy baryons is still scarce, and more experimental and theoretical efforts are encouraged and needed.
In 2002, the SELEX Collaboration reported an evidence of a doubly charmed baryon Ξ + cc that has a mass of 3519 ± 1 MeV in the Λ + c K − π + final state [3] and pD + K − decay mode [4].However, the existence of Ξ + cc (3519) was disfavored by the following FOCUS, BaBar, Belle and LHCb Collaborations [5][6][7][8], and the following theoretical works do not also support this discovery.In 2017, a highly significant structure Ξ ++ cc (3621) with a mass of 3621.40 ± 0.72 ± 0.27 ± 0.14 MeV was observed in the Λ + c K − π + π + mass spectrum by the LHCb Collaboration [9].Subsequently, the LHCb Collaboration measured its lifetime [10,11] and other decay modes [12][13][14][15].Meanwhile, the LHCb Collaboration paid lots of attentions to hunt for more doubly heavy baryons, however no more signal has been found so far [16][17][18][19][20][21].In particular, the Ξ bc and Ξ ′ bc state was placed great expectations to be discovered in the near future, where one charm quark in Ξ cc is replaced by a bottom quark.
In this work, we concentrate on the bottom-charmed family that is made up of a bottom quark b, a charm quark c and a light quark (u, d, or s).When the light quark belongs to up or down quark, the bottom-charmed baryon is named as Ξ bc or Ξ ′ bc ; When the light quark is a strange quark, the bottom-charmed baryon is denoted as Ω bc or Ω ′ bc .The study of bottom-charmed baryons does not only provide an opportunity for us to investigate the heavy quark symmetry and chiral dynamics simultaneously, but also supplies a unique platform about the conventional baryons with three non-identical quarks.Moreover, the excited states and their strong decay behaviors are essential for better understanding the spectroscopy for bottom-charmed baryons and helping experimentalists to hunt for more hadrons.Among various properties, the Okubo-Zweig-Iizuka-allowed (OZI-allowed) twobody decay processes are particularly interesting where one light meson is, especially the pion, emitted from the parent baryon.In this process, the meson couples to the light quark, and the heavy quark subsystem behaves simply as a spectator.Hence this provides a good platform to investigate dynamics of chiral symmetry at the single quark level.However, the study on strong decays for bottom-charmed baryon is very scarce [71], and then it is time to explore this topic systematically.
In a previous work, a nonrelativistic quark model was adopted to study the mass spectrum of heavy baryons with two identical quarks, and gained significant achievements [79].However, the case of all three quarks with different masses were not studied.Here, we employ the same quark model to bottom-charmed system by solving the threebody Schrödinger equation in order to get the mass spectrum consistently.Besides the masses, the realistic wave functions are obtained simultaneously, and can be used for strong decay calculations.The lowest Ξ bc , Ξ ′ bc , Ω bc , and Ω ′ bc states are predicted to be about 6979, 6953, 7109, and 7092 MeV, respectively.Our results indicate that some of λ−mode Ξ bc (1P), Ξ ′ bc (1P), Ω bc (1P), Ω ′ bc (1P) states are narrow, which have good potentials to be observed by future experiments.Also, the strong decays of the low-lying ρ−mode and ρ−λ hybrid states are highly suppressed and can be searched for in the processes with electromagnetic and weak interactions.
This paper is organized as follows.The formalism of potential model and pseudoscalar meson emissions is briefly introduced in Sec.II.We present the numerical results and discussions for the bottom-charmed baryons in Sec.III.A summary is given in the last section.

A. Potential Model
In order to calculate the spectrum of the low-lying bottomcharmed baryons, we adopt a nonrelativistic quark model, where the Hamiltonian can be expressed as with the kinetic energy and the effective potential Here T CM is the center-of-mass energy, r i j is the distance between the ith and jth quarks.The linear confinement potential V con f r i j and one-gluon-exchange potential V coul r i j are as follows where C is the overall zero-point-energy parameter.The spindependent interaction V S D r i j is the sum of spin-spin term V S S r i j , spin-orbit term V LS r i j , and tensor term V T en r i j The parameters are taken from the original work [79] and listed in the Table I.Also, the overall constant C is adjusted by fixing the mass of ground state Ξ cc (3621), which is suitable for present bottom-charmed baryons systems.The bottom-charmed baryon is a three-quark system, where the Jacobi coordinates can be introduced to eliminate the center-of-mass energy.As illustrated in Figure 1, the Jacobi coordinates are defined as where r i and m i denote the position vector and the mass of the ith quark, respectively.The ρ stands for the relative coordinate between bottom and charm quarks, the λ represents the relative coordinate between the light and heavy subsystems, and the R is mass center coordinate.
The bottom-charmed baryons can be divided into four types according to the total wave functions, which are denoted as Ξ bc Ξ ′ bc , Ω bc , and Ω ′ bc , respectively.These notations follow those of the singly heavy baryons Ξ c , Ξ ′ c , Ξ b , and Ξ ′ b , and FIG.1: A typical sketch of the bottom-charmed baryons the naming scheme can be referred to Refs.[1,31,36].In the present work, the orbital wave functions are expanded in terms of a set of Gaussian basis functions that forms an approximate complete set [80,81].Then, the explicit expression for spatial part can be written as with and where C n,N are the expansion coefficients.It is worth noting that the ρ and λ mode can be hardly separated in the strict sense for the systems with three non-identical masses.However, in view of the heavy quark symmetry, the orbital wave functions for bottom-charmed baryons can be separated approximately and we prefer to adopt the above trial wave functions to solve the three-body Schrödinger equation.The range parameters of Gaussian functions are given as where n and N are the number of Gaussian functions, and a and A are the ratio coefficients.According to the Rayleigh-Ritz variational principle, one can have where the H i j are the matrix elements in the total color-flavorspin-orbital bases, E stands for the eigenvalue, and C j are the relevant eigenvectors.Finally, the spectrum of bottomcharmed baryons can be obtained by solving the generalized eigenvalue problem.
In constructing the total wave functions, we adopt the j − j coupling scheme to investigate the bottom-charmed baryons, where the states are defined as The label l ρ is the orbital quantum number between the two heavy quarks, the l λ the orbital quantum number between the light and heavy subsystems, S ρ the spin quantum number of two heavy quarks, J ρ the total angular momentum of heavy quark subsystem, j the total angular momentum of light quark system that is usually known as the light quark spin, and J P the spin-parity for hadrons.More details of different coupling schemes and their relations can be found in our previous work [77].

B. Pseudoscalar meson emissions
Besides the mass spectrum, the strong decays can reflect the internal structures of hadrons more explicitly.In this subsection, the approach of strong decays for bottom-charmed baryons is introduced briefly.In the quark model, the pseudoscalar meson can couple to the light quark inside a bottomcharmed baryons through the Yukawa interaction, which is considered to contribute predominantly to one-meson emission decays 2. The axialvector coupling between the pseudoscalar meson and a light quark can be written as where the q stands for the quark field, g q A is the quark-axialvector coupling and f p is the decay constant.In present work, f π = 93 MeVand f K =111 MeV are adopted that have been widely used in quark model calculations [72,[82][83][84][85][86][87].

2: One meson emission for bottom-charmed baryons
The wave function for the Y bc (Y bc = Ξ bc , Ξ ′ bc , Ω bc , and Ω ′ bc ) baryon with mass M Y bc in the rest frame can be expressed in the momentum representation as Then, the decay amplitude for Y i bc (P i ) → Y f bc (P f )+M p (q) can be obtained by and the q λ is defined as The g f denotes the flavor matrix, M i is the mass of initial baryon, M f is the mass of final baryon, q = (ω, q) is the 4momentum of outgoing pseudoscalar meson, m i is the constituent quark mass with m 1 = m b , m 2 = m c , and m 3 = m u/d/s .After calculating the mass spectrum, one can get the wave functions for bottom-charmed baryons, which are applied to estimate the root mean square and the range parameters α λ .These effective values α λ are obtained by equating the root mean square radius of simple harmonic oscillator wave functions to that obtained in the nonrelativistic quark model, which has been widely used in the previous studies of strong decays [88][89][90][91][92].Then, the helicity amplitude A h can be derived from the transition operator and effective parameters α λ in the harmonic oscillator wave functions.To calculate the strong decay widths for the pseudoscalar meson emissions, one also need to take into account the phase space factor.Finally, the strong decays for bottom-charmed baryons can be obtained within the helicity bases straightforwardly,

III. RESULTS AND DISCUSSION
In this section, we first calculated the mass spectrum of bottom-charmed baryons, and then estimate the strong decays for λ−mode low-lying excited states.Also, we discuss the relativistic corrections for the Roper-like resonances, the mixture of λ−mode states with the same spin-parity, and the suppression of strong decays for ρ−mode and ρ − λ hybrid states.In analogy with singly heavy baryons [93], we adopt the symbols ∼, ∧, ∨, and ⌣ on the top of capital Ξ and Ω to denote the ρ−mode P−wave, ρ−mode D−wave, ρ − λ hybrid, and ρ−mode radially excited states, respectively.II, III, IV, and V. Also, the full mass spectra are plotted in Figure 3 and 4 for reference.

States
n ρ n λ l ρ l λ S ρ J ρ j J P α λ Mass Ξ bc (1S ) 0 0 0 0 1 1 1 2 + 500 7354 + 291 7732 2 ) 0 0 0 2 0 0 5 2 + 350 7609 We first refer to the ground states for bottom-charmed baryons.From the Table VI and Figure 5, The lowest Ξ bc , Ξ ′ bc , Ω bc , and Ω ′ bc states are predicted to be about 6979, 6953, 7109, and 7092 MeV, respectively.It can be seen that our results are consistent with some works [23,24,36], while differ with others about 50 ∼ 100 MeV [27,28,43].In the singly heavy baryons Ξ (′)  c/b , the lowest states are Ξ c/b rather than Ξ ′ c/b .However, for the bottom-charmed baryons, the lowest states are Ξ ′ bc and Ω ′ bc , which is totally different with singly heavy sector.This is due to the different masses in spin-spin term V S S r i j in the nonrelativistic quark model for singly and bc − Ω ′ bc are 44 and 33 MeV, respectively.These quite small mass gaps are caused by spin-spin interaction which is inversely proportional to quark masses.Because of this, the pion emission between the ground states for bottom-charmed baryons is prohibited and only decay mode is by electroweak processes, which is in contrast with the strange sector allowing the strong decay Ξ * → Ξπ.In particular, For the lowest Ξ (′)  bc and Ω (′) bc states, future experiments can search for them in via b → c weak transitions.
Unlike singly heavy baryons, the heavy quark subsystem is more easily excited owing to its larger reduced mass, and then the ρ−mode excited states are lower than the λ−mode ones.It can be seen that our calculated mass spectra for bottom-charmed baryons faithfully reflect this specific feature.Meanwhile, the fine structures of excited states are small and the spectra are highly degenerate, especially for the D−wave states.Therefore, the low-lying P−wave excitations are more likely to be recognized both theoretically and experimentally.Moreover, the mass spectra for bottom-charmed baryons show quite similar patterns as other heavy-light systems, such as conventional charmed or bottom mesons, which suggests that the approximate light flavor SU(3) symmetry and heavy super-flavor symmetry are preserved well.Indeed, analogous to the doubly bottom baryons [77], the bottom and charm quarks stay close to each other like a static color source, and the light quark is shared by these two heavy quarks.
Here, we can also discuss the theoretical uncertainties arising from parameters α ss and α so for doubly heavy baryons.We first vary α ss in the range of 1.10∼1.30,and find that the discrepancies between theoretical results and experimental data for the mass splittings of singly heavy baryons can reach up to 39 MeV.Even with this large variation, the uncertainties for doubly heavy baryons are about 7 MeV, which are small enough.The same procedure should also be done for α so , but the established excited heavy baryons relevant with this pa- [36] [28] [52] [55] [38] [24] [27] [26] Mass(MeV) rameter are few.Here, we vary this value in a wide range of 0 ∼ 0.154, and find the uncertainties for doubly heavy baryons are about 10 MeV, which suggests that our predictions are also stable against the parameter α so .
B. Strong decays for λ−mode Ξ bc and Ω bc states The strong decays for λ−mode Ξ bc and Ω bc states are calculated and listed in Table VII and XI.For the two Ξ bc (1P) states, they have the light quark spins j = 1/2 and 3/2 , and the decay widths for Ξ bc (1/2 − , 1/2) and Ξ bc (3/2 − , 3/2) states are about 95 and 21 MeV, respectively.Owing to the limited phase space, both of them can only decay into the Ξ bc π channel.Because of relevant partial waves of decaying channels, we have found a rather broad J P = 1/2 − state with S −wave decay and a narrow J P = 3/2 − state with D−wave decay, which are roughly consistent with previous work [71].The distinctions of predicted decay widths may arise from the different phase spaces, wave functions, and phenomenological models.
For the two Ω bc (1P) states, they lie below the Ξ bc K threshold, and then the OZI-allowed strong decay is forbidden.The dominating decay channels should be Ω bc π and Ω bc γ.This situation is quite similar to D * s0 (2317) and D s1 (2460) resonances, where the isospin breaking decay and radiative decay modes dominate.Actually, based on the heavy super-flavor symmetry [37], the two Ω bc (1P) states can be related to the charmed mesons and Ω bb (1P) states [77].More theoretical studies on this topic can help us to better understand the D s0 (2317) resonance.
From the Table XI, it can be seen that most of the 2S and 1D states are broad, which can hardly be observed in experiments.For the Ω bc (2S ) state, the total width is relatively small, and the dominant decay mode is Ξ bc K within leading terms of the nonrelativistic transition amplitude Eq. ( 21).
For the radially excited Ξ bc (2S ) and Ω bc (2S ) states, our calculated decay widths are relatively narrow, and the dominant decay modes are Ξ ′ * bc π and Ξ ′ * bc K with the nonrelativistic transition amplitude Eq. ( 21), respectively.These relatively narrow decays width for radially excited state are usually obtained within the nonrelativistic reduction of the axial-vector coupling between the pseudoscalar meson and light quark.This is due to the kind of selection rule from the structure of the transition operator in the leading order of nonrelativistic expansion, and consequently due to the orthogonality of the orbital wave functions between initial and final baryons.We will continue to discuss the relativistic corrections for Roper-like resonances in the following subsection.
For the Ξ ′ bc (1D) and Ω ′ bc (1D) states, they can be divided into j = 3/2 and j = 5/2 triplets.the j = 3/2 triplet are predicted to be relatively broad, while the j = 5/2 are narrow states.It can be seen that the approximate heavy quark symmetry is preserved well in present calculations.

D. Relativistic corrections of order 1/m 2 for Roper-like resonances
The relativistic corrections of order 1/m 2 for convetional baryons are investigated in Refs.[83,84].The authors found that these corrections of order 1/m 2 are significant for the radially excited states, that is Roper-like resonances, while the effects for P−wave and D−wave states are small enough.In the present work, we also investigate these relativistic corrections for six radially excited bottom-charmed baryons.
We only take into account the ground states in the final states for comparison.The results are listed in Table XI, and it can be seen that these relativistic corrections for Roperlike resonances are significant.This specific feature has been found in other Roper-like resonances in the literature.Owning to the lack of experimental information for bottom-charmed baryons, more theoretical and experimental efforts are needed for further exploration.

E. Mixing
In our calculation, we adopt the j − j coupling scheme and the basis in the heavy quark limit to study mass spectra and strong decays for bottom-charmed baryons.Due to the finite mass of heavy quark subsystem, the physical observed resonances may correspond to the superposition of theoretical states in the quark model.For instance, the mixing scheme for λ−mode Ξ ′ bc (1P) and Ω ′ bc (1P) states can be formulated as where θ is the mixing angle.Also, the D−wave excited states with the same spin-parity can mix with each other.
In the heavy quark limit, the mixing angle should be zero.Actually, the heavy quark subsystem including bottom and charm quarks are relatively heavy, and heavy quark symmetry should be approximately preserved.Then, the above mixing angles are expected to be small enough.We calculate the mixture for λ−mode Ξ ′ bc (1P) states as an illustration and list them in Table XII.It can be found that the mixing angles are tiny, and the mixing effects for spectroscopy can be neglected in bottom-charmed sector.Moreover, there may also exist other mixing scheme, which are believed to be even smaller in these bottom-charmed baryons.) F. Low-lying ρ−mode and ρ − λ hybrid states In addition to the presence of λ−mode excitations, the bottom-charmed baryons also have low-lying ρ−mode and ρ − λ hybrid states.For example, when l ρ = l λ = 1, there exists thirteen Ξbc and five Ξ′ bc states, which can be seen in Figure 6.
For these low-lying states, the light meson emissions are supposed to be dominating.However, under the spectator assumption for the two heavy quarks, the orbital wave functions of heavy quark subsystems between initial ρ−mode or ρ − λ hybrid states and final ground states are orthogonal, which results in the vanishing amplitudes and strong decay widths.More explicitly, the light meson emission occurs for λ−mode excitations, and is irrelevant to the ρ−mode variables.Hence, the orthogonality of the different ρ−mode wave functions leads to the vanishing matrix element, which is shown in Figure 7.That is to say, our calculated decay modes preserve the heavy diquark symmetry automatically, where the heavy quark subsystems with different quantum numbers cannot transit into each other.Thus, these states should be extremely narrow and the weak and radiative decays may become dominating, which can provide good opportunities to be searched by future experiments.

IV. SUMMARY
In this work, we have studied the low-lying mass spectra for bottom-charmed baryons in a nonrelativistic quark model by solving the three-body Schrödinger equation.With the obtained realistic wave functions, we get the root mean square radius and study the strong decays of bottom-charmed baryons.The lowest Ξ bc , Ξ ′ bc , Ω bc , and Ω ′ bc states are predicted to be about 6979, 6953, 7109, and 7092 MeV, respectively.Our results indicate that some of λ−mode Ξ bc (1P), Ξ ′ bc (1P), Ω bc (1P), Ω ′ bc (1P) states are relatively narrow, which have good potentials to be observed by future experiments.Also, the strong decays of the low-lying ρ−mode and ρ − λ hybrid states are highly suppressed and can be hunted for in the electroweak processes.
Given the heavy quark symmetry, the heavy quark subsystem in bottom-charmed baryons play a role as a heavy antiquark, and then the heavy super-flavor symmetry emerges.Indeed, our results about the mass spectra and strong decays for bottom-charmed baryons support this claim.With the development of the large-scale accelerator facilities, we expect that more theoretical and experimental efforts are involved to search for more doubly heavy baryons and better understand the heavy quark symmetry.

A
. Mass spectrum In the quark model, for the Ξ bc or Ω bc family, there should be one ground state, five ρ−mode P−wave states, two λ−mode P−wave states, two ρ−mode D−wave states, two λ−mode D−wave states, thirteen ρ − λ hybrid D−wave states, and two radially excited states; for the Ξ ′ bc or Ω ′ bc family should exist two ground states, two ρ−mode P−wave states, five λ−mode P−wave states, six ρ−mode D−wave states, six λ−mode D−wave states, five ρ − λ hybrid D−wave states, and four radially excited states.Within the potential model, we obtain the masses and perform them in Table

FIG. 3 :
FIG.3:The calculated mass spectra for Ξ bc and Ω bc families

C.
Strong decays for λ−mode Ξ ′ bc and Ω ′ bc states

FIG. 7 :
FIG. 7: Schematic diagrams for strong decays of doubly heavy baryons within different exited modes.

TABLE I :
The relevant parameters adopted in this work.

TABLE II :
The mass spectrum of Ξ bc in MeV.

TABLE III :
The mass spectrum of Ω bc in MeV.

TABLE IV :
The mass spectrum of Ξ ′ bc in MeV.

TABLE V :
The mass spectrum of Ω ′ bc in MeV.

TABLE VI :
Masses of ground states for Ξ (′) bc and Ω(′)bc baryons compared with different calculations.The units are in MeV.

TABLE VII :
The predicted strong decay widths of Ξ bc (1P) and Ω bc (1P) states in MeV.The • • • stands for the closed channel.

TABLE X :
The predicted strong decay widths of Ξ ′ bc (2S , 1D) and Ω ′ bc (2S , 1D) states in MeV.The • • • stands for the closed channel, and × denotes the forbidden channel due to quantum numbers.