Spectroscopic survey of higher-lying states of B c meson family

In this work, we investigate the complete spectroscopy of the B c mesons, with a special focus on the consideration of the unquenched e ff ects. To account for such e ff ects, we employ the modified Godfrey-Isgur model and introduce a screening potential. The resulting mass spectrum of the concerned higher B c states is then presented, showing significant deviations after considering the unquenched e ff ects. This emphasizes the importance of considering the unquenched e ff ects when studying of the higher B c mesons. Furthermore, we determine the corresponding spatial wave functions of these B c mesons, which have practical applications in subsequent studies of their decays. These decays include two-body Okuba-Zweig-Iizuka allowed strong decays, dipion transitions between B c mesons, radiative decays, and some typical weak decays. With the ongoing high-luminosity upgrade of the Large Hadron Collider, we expect the discovery of additional B c states in the near future. The knowledge gained from the mass spectrum and the di ff erent decay modes will undoubtedly provide valuable insights for future experimental explorations of these higher B c mesons.


I. INTRODUCTION
The study of hadron spectroscopy offers a unique avenue to deepen our understanding of the non-perturbative nature of the strong interaction.Over the past few decades, significant progress has been made, both experimentally and theoretically, leading to extensive investigations of the exotic hadronic states, such as the charmonium-like states XYZ and the hidden-charm pentaquark states P c /P cs [1][2][3][4][5][6][7][8][9][10][11][12][13][14].Moreover, notable observations of light flavor hadrons by the BE-SIII Collaboration [15], as well as heavy flavor hadrons by the LHCb [16][17][18][19] and Belle collaborations [20][21][22], indicate that the construction of the conventional hadron family is an ongoing endeavor.The abundance of these phenomena in hadron spectroscopy underscores the field's position at the forefront of precision particle physics.
In contrast to the well-established charmonium and bottomonium families, the B c meson family remains relatively unexplored, with only a few low-lying B c states reported in experiments.Pivotal theoretical contributions were made by Refs.[23,24], which proposed the experimental detection of the B c mesons through the hadron colliders.In 1998, the CDF Collaboration reported the observation of a B c meson with the mass of M = (6.40 ± 0.39 ± 0.13) GeV, identified via the B ± → J/ψℓ ± ν decay [25].This measured mass is consistent with expectations for the ground state of the B c meson [24,[26][27][28][29].However, the full establishment of the B c meson family remains a formidable task, requiring extensive efforts to identify and explore its properties.
In 2014, the ATLAS Collaboration reported the observation of a structure consistent with the predicted B c (2S ) state, with a mass of (6842 ± 9) MeV [43].Additionally, the CMS and LHCb collaborations observed the excited B c (2 1 S 0 ) and B c (2 3 S 1 ) states in the B + c π + π − invariant mass spectrum, with masses determined as (6872.1 ± 2.2) MeV and (6841.2± 1.5) MeV, respectively [44,45].However, the current Particle Data Group (PDG) includes only two B c mesons, namely the B c (1S ) and B c (2S ) states [46].The limited experimental knowledge of the complete B c family, particularly the higher excited B c states, highlights the necessity for further studies in this area.
The study of the B c meson family plays a crucial role in advancing our understanding of the strong interaction and the non-perturbative regime of Quantum Chromodynamics (QCD).With the forthcoming high-luminosity upgrade of the Large Hadron Collider (LHC), there are promising opportunities to investigate the higher excited states of the B c meson, and it is anticipated that additional experimental data will become available in the near future.
Some theoretical studies have examined the B c spectrum using the quenched quark models, as discussed in references [47][48][49][50][51][52][53].Among these models, the Godfrey-Isgur (GI) model [47], proposed by Stephen Godfrey and Nathan Isgur in 1985, has demonstrated notable success in describing the spectra of low-lying hadrons.Additionally, other quenched quark models have exhibited remarkable achievements in predicting the spectra of low-lying mesons.However, it is now widely acknowledged that the unquenched effects play a significant role, as they can resolve the low-mass puzzles of several hadrons, such as the D s0 (2317) [54][55][56][57], D s1 (2460) [55,58], X(3872) [59][60][61], and Λ c (2940) [62].Therefore, when exploring the spectroscopy of higher excited states of the B c meson, it is imperative to consider the unquenched effects as well.
In this study, we employ the modified Godfrey-Isgur (MGI) model [63][64][65][66][67] to account for the unquenched effects.To reflect such effects, we utilize a screening potential introduced in Refs.[68,69].Furthermore, this method has been applied to the study of bottomonia [65] and charmonia family [66,67], and has yielded the predictions, providing valuable guidance for experimental investigations.It should be noted that the couple channel effect and the screening potential have a similar effect on the hadron masses, i.e. the higher excited state masses are reduced [60,70,71].In the following section, we provide a concise overview of the MGI model and present the mass spectrum of the B c family, along with a comparison to the quenched quark model, in order to elucidate the distinctions between the two models.Furthermore, we present the numerical spatial wave functions obtained from the MGI model, which serve as crucial inputs for investigating various properties.Specifically, we calculate the magnetic moments of the higher excited B c states to reveal their internal structures.This analysis can be valuable in differentiating between conventional B c states and exotic B c -like molecular tetraquark states with identical quantum numbers and similar masses.
To provide a comprehensive theoretical analysis, we also calculate the two-body Okubo-Zweig-Iizuka (OZI)-allowed strong decays, dipion decays, radiative decays, and typical weak decays by employing the numerical spatial wave functions of the relevant mesons.In the concrete calculations, the quark pair creation (QPC) model [72][73][74][75] is utilized to describe the behavior of the two-body OZI-allowed strong decays, while the electric dipole (E1) and magnetic dipole (M1) radiative transitions are analyzed by considering the radial decays of the higher excited B c states.For the investigations of dipion transitions, we adopt the quantum chromodynamics multipole expansion (QCDME) method.Additionally, we employ the covariant light-front quark model (CLFQM) to calculate a series of weak transition form factors and their corresponding weak decays.The spatial wave functions play a crucial role in these decay processes, and by utilizing the MGI model, we obtain reliable numerical results.As experimental data continue to accumulate, these decay processes can be further explored and potentially measured with higher precision.
This paper is organized as follows.In Section II, we present our analysis of the mass spectrum of B c mesons using the MGI model with the screening effects, and provide the associated numerical spatial wave functions.Section III is dedicated to the predictions of two-body OZI-allowed strong decays of the considered B c mesons by employing the QPC model.In Section IV, we investigate the dipion transitions between the B c states.Section V focuses on the predictions of the decay widths for the E1 and M1 radiative transitions and the magnetic moments of these higher B c mesons.Furthermore, in Section VI, we explore some concerned weak transition form factors and the corresponding weak decays for the B c meson.Finally, our findings are summarized in Section VII.

II. MASS SPECTRUM AND THE CORRESPONDING SPATIAL WAVE FUNCTIONS
In this study, we employ the MGI model to calculate the mass spectrum of the higher excited B c states.To accurately account for the screening effects, we incorporate a screening potential into our calculations.Furthermore, we derive the corresponding numerical spatial wave functions, which play a pivotal role in the investigations of the decay properties of the B c states.

A. The mass spectrum
To obtain the mass spectrum of the B c mesons, we apply the MGI model by introducing the screening potential [63-67, 76, 77].The involved Hamiltonian is where m 1 and m 2 denote the masses of the b and c quarks, respectively.V eff (p, r) = H conf + H hyp + H so is the effective potential representing the q ′ q interaction.In the non-relativistic limit, V eff (p, r) is transformed into the familiar nonrelativistic potential V eff (r), which can be written as where is the spin-independent potential, which includes the screening potential and the Coulomb-like potential.Here, α s (r) is a running coupling constant, and µ is a parameter which stands for the strength of the screening effects.The colour hyperfine interaction is given by where S 1(2) is the spin of quark or antiquark.And H so = H so(cm) + H so(tp) is the spin-orbit interaction, where is the color-magnetic term resulting from the one-gluon exchange, while denotes the Thomas precession term with the screening effects, where L is the orbital angular momentum between quark and antiquark.In addition, we need to smear the screened potential where is the smearing function, and σ is a smearing parameter.Then, we introduce the momentum dependent factors as follows , where ϵ i correspond to different types of the interactions, and Ṽi (r) are the effective potentials included in Eqs.(2.4)- (2.6).
For mesons composed of heavy quarks with equal masses, such as the charmonium and bottomonium states, the L-S coupling scheme is appropriate and the meson state can be labeled by the notation n 2S +1 L J .However, for mesons with constituents of different masses, such as the B c mesons, the spin-dependent terms in the Hamiltonian can mix the spinsinglet and spin-triplet states, and the resulting mixing states can be expressed as where θ represents the mixing angle.
To solve the Schrödinger equation [78,79], we use the simple harmonic oscillator (SHO) wave functions as a set of complete bases, i.e., where (2.12) Here, the radial wave functions are denoted by R nL (r, β) and R nL (p, β) in the coordinate and momentum spaces, respectively.Y L,M L (Ω r ) and Y L,M L Ω p stand for the spherical harmonic functions, and L L+  is a phenomenological parameter in the SHO wave function, and we set it to be β = 0.5 GeV in the calculation.
In Table I, we list the reported experimental data [46] of the masses of bottom-charmed, charmonium, bottomonium, charmed, charmed-strange, bottom, and bottom-strange mesons, which are quoted from the PDG [46] and are used to fit the parameters of the MGI model.Due to the limited availability of experimental data for the B c mesons, we also include experimental data from other mesons as auxiliary.By fitting the experimental data in Table I, we can obtain the mass spectrum of the B c mesons.Here, the accuracy of the fitting data is judged based on the χ 2 criterion as follows where m Th i , m Exp i , and m Er i are the theoretical value, the experimental value, and the error of the i-th data, respectively.Here, the errors m Er i = 1 MeV are the uniform values for all considered mesons.In the absence of experimental data on the masses of the B c mesons, we used the experimental data of other heavy flavor mesons listed in Table I to fit the parameters of the model.To make the mesons act in the same proportions in our fit, we chose a universe value of 1 MeV as the uncertainty.Thus, in this work we used it only as a mathematical tool to fit the model parameters, rather than a χ 2 fit in the usual sense.In this case, we have not given the χ 2 /d.o.f.value and have not considered the uncertainties of the parameters.
In Table II, we present the fitted parameters of the MGI model, and these parameters were determined by selecting the minimum χ 2 value.By utilizing the parameters in Table II, we can calculate the mass spectrum of the B c mesons, which are displayed in Table III.In addition to the mass spectrum, the Table III also includes the mixing angles of the P-wave, D-wave, and F-wave B c mesons.Furthermore, we provide a comparison with other theoretical studies [48,52,80,81].
In Fig. 1, we present a comparative analysis of the B c meson mass spectrum, considering both the screening effects and the results obtained by the GI model [52].Our findings reveal the significant impact of the screening effects on the mass spectrum of the higher excited B c states, as compared to the predictions based solely on the GI model [52].Several notable examples are highlighted below: (i) The B c (4 1 S 0 ) and B c (4 3 S 1 ) states exhibit mass reductions of approximately 76 MeV and 76 MeV, respectively, in our calculations compared to the GI model predictions [52].
(ii) The impact of the screening effects becomes even more pronounced in the B c (5 1 S 0 ) and B c (5 3 S 1 ) states, with a substantial mass difference of up to 120 MeV observed between our results and those obtained without considering the screening effects [52].This disparity is particularly evident in the higher excited B c states, as depicted in Fig. 1.
Hence, our findings underscore the necessity of accounting for the unquenched effects when investigating the spectroscopy of higher excited B c mesons.Future experimental studies with more precise data are expected to provide a valuable opportunity for testing and validating our predictions.
In the preceding subsection, we derived the mass spectrum of the B c mesons.Furthermore, we can extract the corresponding spatial wave functions, which serve as crucial inputs for investigating their decay properties.Consequently, in this subsection, we present these spatial wave functions to explore their characteristics.The figure 2 illustrates the S /P/D/Fwave spatial wave functions of the B c mesons.Notably, for the spatial wave functions of the S -wave B c mesons, noticeable node effects emerge when n exceeds 4, such as n = 5, and so on.These node effects can substantially impact certain physical results, which will be discussed in Section III.

III. TWO-BODY OZI-ALLOWED STRONG DECAYS
In this section, our focus lies on the study of two-body OZIallowed strong decays of the B c mesons by employing the QPC model in the concrete calculations.The QPC model, also known as the 3 P 0 model, was initially proposed by Micu in 1968 [72] and has since been further developed by the Orsay Group [82][83][84].Over time, the QPC model has been extensively utilized to investigate the two-body OZI-allowed strong decays of various hadrons [73][74][75][82][83][84].Here, we provide a brief introduction to the QPC model, with a specific focus on the transition matrix for the A → B + C process.
In the QPC model, the transition matrix for the A → B + C decay is defined as follows Here, M M J A M J B M J C represents the helicity amplitude, while P B and P C denote the momenta of the daughter mesons B and C, respectively, in the rest frame of the parent meson A. The states |A⟩, |B⟩, and |C⟩ correspond to the mock states associated with mesons A, B, and C, respectively.Furthermore, T represents the transition operator, which describes the creation of a quark-antiquark pair from the vacuum.In the nonrelativistic limit, the transition operator can be expressed as where the quark and antiquark are represented by indices 3 and 4, respectively.The state χ 34  1,−m corresponds to a spintriplet configuration, while ϕ 34 0 and ω 34 0 denote the SU(3) flavor and color singlets, respectively.The parameter γ is a dimensionless constant that characterizes the strength of quarkantiquark pair creation from the vacuum, and its value is determined through the fitting of experimental data.The term Y lm (p) = |p| l Y lm (p) represents the solid harmonic function.According to the Jacobi-Wick formula [85], the helicity amplitude can be transformed into the partial wave amplitude, which can be expressed as where L is the orbital angular momentum between final states B and C. The partial width of the A → B + C process can be given by where m A is the mass of the parent meson A.
In addition, the meson wave function is defined as the mock state [86], i.e., where χ A S M S , ϕ A , ω A , and Ψ A nLM L (p 1 , p 2 ) denote the spin, flavor, color, and spatial wave functions of meson A, respectively.In our concrete calculations, we utilize the numerical spatial wave functions obtained in Section II as inputs for these mesons.This approach helps to avoid the dependence on the β values by taking a SHO wave function.Furthermore, since experimental width data is lacking, we adopt the value of γ for q q as 0.4 √ 96π from Ref. [87], while the creation strength for s s satisfies γ s = γ/ √ 3.In Table IV, we present the two-body OZI-allowed strong decay widths and the corresponding branching ratios of the S -wave and P-wave B c states with n = 3, 4, 5 and n = 3, 4, respectively.These discussed B c (n 1 S 0 ) states with n = 3, 4, 5 can decay into the B * D channel.As we increase the principal quantum number n, such as n = 5, we observe a significant reduction in the partial width of their decays into the B * D channel.This reduction is primarily due to the node effects of the corresponding spatial wave functions, as illustrated in Fig. 2. For the B c (3 1 S 0 ) state, the B * D is the unique twobody OZI-allowed decay channel around 62 MeV.Thus, the B * D channel should be the promising channel to observe the B c (3 1 S 0 ) state.For the B c (4 1 S 0 ) state, it mainly decays into the BD * , B * D, and B * D * channels, and whose corresponding branching ratios are more than 99%.For the B c (5 1 S 0 ) state, it dominantly decays into the BD * , B * D, B * D * , and BD( 3 P 0 ) channels.
The node effects are more pronounced for the B c (n 3 S 1 ) mesons compared to the B c (n 1 S 0 ) mesons, with a notable example being the B c (4 3 S 1 ) state.Furthermore, for the B c (5 3 S 1 ) state, the decay widths of the B * s D * s and B * D(1 3 P 0 ) channels are 1.7 × 10 −3 MeV and 2.5 × 10 −3 MeV, respectively, reaching a order of magnitude of 10 −3 MeV.Thus, it is evident that the node effects in the B c (5 3 S 1 ) state are more significant than that in the B c (5 1 S 0 ) state.In the following, we summarize several main points: (i) The total width of the B c (3 3 S 1 ) state is about 82 MeV, and the main decay channels of the B c (3 3  For these 3P states of the B c meson, the partial widths of the B s D s and B * s D s decay channels are significantly suppressed.Specifically, for the B c (3 3 P 0 ) state, the B * s D s channel is kinematically forbidden.Similarly, for the 4P states of the B c meson, the decay channel of B * s D * s is also kinematically forbidden.Furthermore, for the B c (3 3 P 2 ) and B c (4P ′ 1 ) states, the decay widths of the B * D and B(1P 1 )D channels are 2.6 × 10 −3 MeV and 2.3 × 10 −3 MeV, respectively, reaching the order of magnitude of 10 −3 MeV.Our results show that the largest decay width of these discussed 3P-wave B c mesons is the B * D * channel, which has the estimated branching ratios of 63.1%, 54.9%, 49.0%, and 81.6% for the B c (3 3 P 0 ), B c (3P 1 ), B c (3P ′ 1 ), and B c (3 3 P 2 ) states, respectively.Thus, we suggest the future experiments to search for the 3P-wave B c mesons by the B * D * channel.For the 4P states of the B c meson, the largest decay width is also the B * D * channel.Compared with the 3P states of the B c meson, we find the total widths of the 4P states are more smaller than these 3P states of the B c meson.
In Table V, we show the two-body OZI-allowed strong decay widths and the corresponding branching ratios for the Dwave and F-wave B c states.For these 2D states of the B c meson, no decay channels are significantly suppressed, and the largest decay channels are also different.For B c (2 3 D 1 ) and B c (2D ′ 2 ) states, the largest decay width is also the BD * channel, whose branching ratios are 59.8% and 69.5%, respectively.However, the largest decay widths of the B c (2D 2 ) and B c (2 3 D 3 ) states are B * D and B * D * channels, respectively.For the B c (3D) states, the largest decay width is the B * D * channel, whose branching ratios are 58.4%,37.2%, 82.1%, and 64.4%For the B c (1 3 F 2 ) and B c (1 3 F 4 ) states, they can decay into the BD and B * D channels, and their dominant decay mode is the BD channel.In addition, it is worth noting that the 3 ) and B c (1F 3 ) states can decay into the B * D channel with the branching ratios 100%, but their decay widths exist obvious difference, which is similar to the decay behaviour of the B c (1 3 F 2 ) and B c (1 3 F 4 ) states.The B * D * is the dominant decay channel for the 2F states of the B c meson, and the branching ratios are 71.2%,48.0%, 58.6%, and 51.4% for the B c (2 3 F 2 ), B c (2F 3 ), B c (2F ′ 3 ), and B c (2 3 F 4 ) states, respectively.The B s D * s is a significantly suppressing channel for the 2F states of the B c meson, especially for the B c (2F 3 ) and B c (2 3 F 4 ) states, whose branching ratios are 2.9 × 10 −5 and 2.0 × 10 −5 , respectively.
In general, the information on the concerned two-body OZI-allowed strong decays presented of B c mesons in Tables IV-V can provide valuable guidance for further experimental searches for them.

IV. DIPION TRANSITIONS BETWEEN B c MESONS
In this section, we investigate the dipion transitions using the QCDME approach.The QCDME has been widely used to study the dipion or η hadronic transitions between the lowlying heavy quarkonium [65,[88][89][90][91][92][93][94].The idea of this approach is that the QZI-suppressed hadronic transition is represented by the parent heavy quarkonium first emitting a gluon to form an intermediate hybrid state, and then recombining itself to the daughter heavy quarkonium and the light meson(s) (like a η or a pair of pions) with emitting another gluon via the hadronization process.Following Ref. [90], one can calculate the decay width of the dipion transition as depicted by the following expression: where c 1 and c 2 are undetermined parameters, while l i( f ) and J i( f ) are the orbital and total angular momenta of meson A(B), respectively.Both mesons A and B possess an identical spin value denoted as s.Furthermore, the phase-space factors, denoted as G and H in this context, are defined as follows: respectively, with Here, m i and m f are the masses of the initial and final B c mesons, respectively, while m π is the mass of the pion.Additionally, the dynamical-associated part f l i f is written as where R i (r) and R f (r) are the radial wave functions of the parent and daughter B c mesons, respectively, while R kl (r) denotes the intermediate vibrational state.m kl is the mass of the intermediate vibrational state with the radial quantum number k and the orbital angular momentum l.
The spatial wave functions associated with the relevant B c mesons can be given by utilizing the MGI model as mentioned above.However, for the intermediate vibrational states, the appropriate approach involves resorting to the quark confining string (QCS) model, as established by Refs.[95][96][97].The potential of a hybrid state adopted in our calculations is [65] where V G (r) = −4α s (r)/(3r) is the one-gluon exchange potential and V S (r) = σ(r)r + c is the color confining potential with Here, we adopt {α 1 , α 2 , α 3 } = {0.25,0.15, 0.20} and [47,98].The effective vibrational potential V n (r) is chosen as [94,95] V In the calculations, we only consider the lowest string excitation, i.e., n = 1.And then, we choose α 1 = √ 1.5 [65].Before performing the concrete calculations, we should determine the parameters c 1 and c 2 in Eq. (4.1).Since there have no measurement of the dipion decays between different B c states, we would fit these parameters from the bottomomium segment, which is the same as in Ref. [99].Numerically, we take |c 1 | 2 = 61.8× 10 −6 and |c 2 | 2 = 1.93 × 10 −6 [90] in the calculations.
With the above preparations, we calculate the decay rates of dipion transitions between B c states.In Table VI, we present our results of the decay rates of dipion transitions We notice the discrepancies between the results from this work and Refs.[52,99].In fact, there are also notable differences between the results from the GI model [52] and Ref. [99].But, we still find some consistent conclusions for the decay rates, i.e., (i) For the B c (2S ) → B c (1S )π + π − transitions, although our results are about a half of those of Refs.[52,99], but our result is comparable to those of the other two theoretical groups.Thus, we come to the same conclusion about the decay rates for these transitions, which are sizable.(ii) For the B c (2P) → B c (1P)π + π − transitions, our results show that B c ( π + π − have significant decay rates and other transitions are suppressed.In general, this conclusion was supported by the GI model [52] and Ref. [99].(iii) Regarding the B c (1D) → B c (1S )π + π − transitions, it is noteworthy that our result is an order of magnitude smaller than the corresponding result from the GI model [52], but comparable to that of Ref. [99], taking B c (1 3 D 1 ) → B c (1 3 S 1 )π + π − and B c (1 3 D 3 ) → B c (1 3 S 1 )π + π − as examples.Obviously, the study of these dipion transitions of B c mesons should be further pursued both theoretically and experimentally.
In this work, we hope to provide the dipion mass distributions of some discussed dipion transitions of B c mesons, which are shown in Fig. 3.We want to emphasize that Here, the decay rates are given in units of keV.

Initial states Final states
This work GI [52] Ref. [99] the maxima of dipion mass distributions are all normalized to 1, since we only concern on their line shapes.Here, we not only select some decays of low-lying B c mesons far below the BD threshold, i.e., B c (2 3 S 1 ) → B c (1 3 S 1 )π + π − and B c (2 1 S 0 ) → B c (1 1 S 0 )π + π − , but also take some decays of B c mesons above (or near) the BD threshold, like At the present, this experimental information is still missing.Therefore, we strongly suggest our experimental colleagues to carry out the measurement of the dipion mass distributions of the B c dipion transitions, which is sensitive to reflect the difference in the result under the quenched and unquenched pictures.

V. RADIATIVE DECAYS AND MAGNETIC MOMENTS
In this section, we will delve into the electric dipole (E1) and magnetic dipole (M1) radiative decay behaviors and the magnetic moments exhibited by the relevant B c states.These observations hold significant implications for unraveling their intricate internal structures.
Here, we have normalized the maxima of the dipion distributions to 1.

A. The E1 transitions
To commence our explorations, we will delve into the E1 transition decays of the aforementioned B c states.The partial widths associated with the E1 radiative transition in the nonrelativistic quark model can be obtained using the following expression [123,124] with Additionally, we denote m c and m b as the masses of the charm and bottom quarks, respectively.The corresponding charges of the charm and bottom quarks are e c = 2/3 and e b = −1/3, respectively.The energy of the photon, denoted as ω, can be obtained through the conservation of energy and momentum, expressed as follows: where M i and M f are the masses of the parent and daughter mesons, respectively.
In Table VII, we present the widths of electric dipole decays for the discussed B c states and compare them with those obtained from other models.Upon scrutinizing the results in Table VII, we observe that the overall magnitude of our calculations remains below 0.1 MeV, which is several orders of magnitude smaller than the two-body OZI-allowed strong decays calculated earlier.Given the limited availability of experimental data, we can only discuss our findings in relation to other theoretical predictions.When comparing our results with those obtained from the GI model [52], we find substantial agreement in most cases, with only a few discrepancies.These variations could potentially arise from the introduction of the screening effects, such as the transitions B c (1D ′ 2 ) → B c (1P 1 )γ, B c (1D 2 ) → B c (1P ′ 1 )γ, and so on.Now let's delve into a detailed analysis of the obtained results presented in Table VII.Regarding the calculations for the initial state of 1P-wave, our findings exhibit consistency with the GI model [52], as well as with Refs.[49,51], and no instances of underestimation are observed.However, for the cases involving initial states of 1D-wave and 1F-wave, we have noticed a systematic tendency toward underestimation in the results.We notice that these results for the processes B c ( γ, are underestimated.Upon closer examination of the transition processes listed above, we can discern a fascinating pattern: the underestimation occurs consistently across similar transition processes when transitioning from a higher-order excited state to a lower-order excited state.Truly, this phenomenon is a remarkable finding.

B. The M1 transitions
We perform calculations and explore the M1 transition decays of the aforementioned B c states.The M1 transition is frequently employed to investigate the radiative transition behavior of quarkonia.Since the B c meson comprises an antibot-TABLE VII: Widths of the E1 transitions for the 1P-wave, 1D-wave, 1F-wave, 2S -wave, 2P-wave, and 3S -wave B c states compared with those from other theoretical work.Here, the width is in units of keV.[52] Ref. [51] Ref. [49] Initial states Final states This work GI [52] Ref. [51] Ref. [

Initial states Final states This work GI
under the nonrelativistic approximation, where and j 0 (x) is the spherical Bessel function, i.e., In Table VIII, we present the magnetic dipole decay widths of the aforementioned B c states.Upon examining the results, we observe that for a given initial state, the decay width increases as the principal quantum number n decreases.This finding aligns with the conclusions drawn by previous studies [52].Conversely, for a fixed final state, the decay width decreases as the initial state's principal quantum number n increases.This intriguing phenomenon warrants further investigations.We have observed two transitions that exhibit significant suppression, which are B c (4 3 S 1 ) → B c (4 1 S 0 )γ and B c (5 3 S 1 ) → B c (5 1 S 0 )γ, with corresponding values of 1.1 eV and 0.6 eV, respectively.Although experimental progress regarding the radiative transition behavior of B c mesons remains elusive, we expect that our discussions on this topic will offer theoretical insights for future experimental explorations of B c mesons and will prove valuable in this regard.

C. Magnetic moments
We delve into the magnetic moments of the discussed B c states, which serve as essential and significant physical observables of hadrons.Magnetic moments are fundamental quantities that have garnered considerable attention and sparked extensive discussions over the past decades, particularly in relation to the magnetic moments of the decuplet and octet baryons.Numerous theoretical models and approaches have been employed to investigate the magnetic moments of hadronic states, including the constituent quark model, the Bag model, the lattice QCD simulations, the chiral perturbation theory, the QCD sum rule, and others [12].
The magnetic moments of the B c states can be calculated by the following expectation values (5.8) In the case of the B c states, the spatial wave function fulfills the normalization condition, and the color wave function is unity due to the color confinement.Consequently, the magnetic moments of the B c states are determined by their flavor and spin-orbit wave functions.
In this study, we employ the constituent quark model to explore the magnetic moments of the discussed B c states.This approach is analogous to the investigations of the magnetic moments of the decuplet and octet baryons conducted in previous studies [126][127][128].Within the framework of the constituent quark model, the total magnetic moment of the B c state is composed of two distinct components: the spin magnetic moment and the orbital magnetic moment.The zcomponent of the spin magnetic moment operator, denoted as μspin z , and the orbital magnetic moment operator, denoted as μorbital z , can be expressed as   respectively.Here, µ i = q i /(2m i ) represents the magnetic moment of the i-th quark, where q i and m i denote the charge and mass of the i-th quark, respectively.Moreover, σiz cor-responds to the Pauli spin operator associated with the i-th quark.In the context of the B c meson, the subscripts c and b refer to the charm quark and the antibottom quark, respectively.Additionally, Lz represents the z-component of the orbital angular momentum operator between the c and b quarks.
In the case of the S -wave B c states, the contribution from the orbital magnetic moment is zero due to the absence of the orbital angular momentum (L z = 0).Therefore, we only need to consider the contribution from the spin magnetic moment for the S -wave B c states.On the other hand, for the orbital excited B c states, the total magnetic moments comprise both the spin magnetic moment and the orbital magnetic moment.The spin-orbit wave function | 2S +1 L J ⟩ for the orbital excited B c states can be expanded by coupling the orbital wave function Y L,m L with the spin wave function χ S ,m S , resulting in the following expression [168]: where C JM Lm L ,S m S is the Clebsch-Gordan coefficient.Here, we need to point out that the orbital magnetic moment and the spin magnetic moment can be obtained by sandwiching the orbital and spin magnetic moment operators between the relevant spin-orbit and flavor wave functions.
To illustrate the calculations of magnetic moments for the orbital excited B c states, let us consider the nP-wave B c states, which include the B c (n 3 P 0 ), B c (n 3 P 1 ), B c (n 1 P 1 ), and B c (n 3 P 2 ) states.In this case, due to the mixture of the spin-singlet and spin-triplet states, it is necessary to expand their spin-orbital wave functions | 2S +1 L J ⟩ as follows: , (5.12) Using the previously expanded spin-orbital wave functions, we can proceed to calculate the magnetic moments and the transition magnetic moments of the nP-wave B c states.The explicit expressions for these quantities are (5.13) where we introduce the notation µ L c b = m c µ¯b m c +m¯b + m b µ c m c +m¯b to simplify the magnetic moment calculations.It is important to note that the calculation method for the transition magnetic moments between the nP-wave B c states is similar to that of the magnetic moments of the nP-wave B c states, with the only difference being the wave functions of the initial and final states.For a more detailed calculation of the transition magnetic moments, please refer to Refs.[129-131, 140, 146, 148, 152, 154-159, 162, 166-168].
For the nP-wave B c states, it is important to consider the mixing of the B c (n 1 P 1 ) and B c (n 3 P 1 ) states, as discussed in Section II.Referring to Eq. (2.10), we can obtain the following expressions: (5.14) Here, θ nP represents the mixing angle between the B c (n 1 P 1 ) and B c (n 3 P 1 ) states, as provided in Table III.With these preparations, we can proceed to calculate the magnetic moments of the mixed states of the B c (n 1 P 1 ) and B c (n 3 P 1 ) states.The magnetic moments of these mixed states can be calculated using the following relations: (5.15) Therefore, the magnetic moments of the mixed states of the B c (n 1 P 1 ) and B c (n 3 P 1 ) states are determined not only by the magnetic moments and the transition magnetic moments of the B c (n 1 P 1 ) and B c (n 3 P 1 ) states but also by the mixing angles of the B c (n 1 P 1 ) and B c (n 3 P 1 ) states.Within the framework of the constituent quark model, the masses of the involved quarks play a crucial role in the study of their magnetic moment properties.In our calculations, we adopt the quark masses m c = 1.660GeV and m b = 4.730 GeV, which have been widely used to describe the hadronic magnetic moments quantitatively [133,167,[172][173][174][175][176].Table IX presents the magnetic moments of the S /P/D/F-wave B c mesons obtained from our calculations.Furthermore, we compare our results with those from other theoretical works and find that our findings are in good agreement with the theoretical predictions in Ref. [167].It is important to note that the investigations of the magnetic moments of B c states have not received much attention thus far.Therefore, we hope that our work will stimulate further theoretical and experimental efforts to explore the magnetic moments of the B c states.
As widely acknowledged, phenomena in the higher mass region are highly intricate.This domain encompasses various conventional B c states, which possess identical quantum numbers and similar masses, as well as predicted B c -like molecular states [178].Distinguishing between these states poses a critical challenge for both theoretical and experimental aspects.For instance, consider the conventional B c (2P ′ 1 ) state and the DB * molecular state with I(J P ) = 0(1 + ), which have closely aligned masses.However, it has been observed that their magnetic moments exhibit evident differences, i.e., µ B c (2P ′ 1 ) = 0.486µ N , µ DB * [0(1 + )] = 0.532µ N [167].There are many such examples.Consequently, investigating the magnetic moment properties provides a means to differentiate between states sharing identical quantum numbers and similar masses.
In theoretical aspects, the light-front quark model is a powerful phenomenological model to calculate the weak transition matrix elements of mesons and baryons decays.Specially, the B c weak decays have been studied by the CLFQM in Refs.[191][192][193][194][195][196][197][198].However, in the concrete calculations, the important input, i.e., the spatial wave function, is usually adopted as the Gaussian-like form, and may bring large uncertainty due to the phenomenological parameter β.In this work, we also employ the CLFQM to revisit the B c weak decays.In our calculations, we adopt the numerical spatial wave functions of the involved mesons, which benefit from the MGI model introduced in Section II, as opposed to SHO wave functions with a phenomenological parameter β used in previous works based on the CLFQM [191][192][193][194][195][196][197][198].This approach reduces the dependence on phenomenological parameters in determining the form factors.
In order to make a comprehensive discussion, we first discuss the B c → M weak transition form factors, where M represents a pseudoscalar meson (P), vector meson (V), scalar meson (S ), or axial meson (A).These transitions involve the quark-level transitions b → c(u) and c → s(d).Additionally, we employ these form factors as inputs to investigate various weak decays, including semileptonic decays and typical twobody nonleptonic decays.
Generally, the B c → M transitions induced by the V − A current can be expressed as [200] where p ′ and p ′′ are the momenta of the initial state meson B c and the final state meson P/V/A/S , respectively.Besides, we define P µ = p ′ µ + p ′′ µ and q µ = p ′ µ − p ′′ µ , while the convention ϵ 0123 = +1 is used.
In the frame of the CLFQM, the constituent (anti)quark inside a meson system are off-shell.The parent and daughter mesons have the four momenta P ′ = p ′ 1 + p 2 and P ′′ = p ′′ 1 + p 2 , where p ′(′′) 1 and p 2 are the four momenta of the quark and the antiquark, respectively.These momenta can be expressed in terms of the following internal variables (x i , ⃗ k ′ ⊥ ) (i = 1, 2): where they must satisfy the relation The B c → M weak transition form factors have been widely studied by CLFQM in Refs.[191][192][193][194][195][196][197].As derived in Refs.[192,195,197,[201][202][203], the B c → P weak transition form factors are where with Analogously, the B c → S weak transition form factors are written as [193,195,197,202,203] where The B c → V transition form factors can be evaluated by where g, f , a + , and a − are the scalar functions of q 2 .Their expressions in CLFQM are [191,192,195,197,[201][202][203] ) where (6.15) For the B c → 3 A( 1 A) transition, the form factors can be evaluated by the relations: where q 3 A( 1 A) , l 3 A( 1 A) , and c 3 A( 1 A) ± are functions of q 2 , with the concrete expressions in the CLFQM as [191,197,202] ⃗ k ′ ⊥ • ⃗ q ⊥ q 2 , (6.17) + Z 2 (3A (1)  2 − 2A (2)  4 − 1) + where ) and ω ′′ 1 A = 2.In the previous theoretical works [191][192][193][194][195][196][197], the phenomenological Gaussian-type wave functions are widely used.However, this treatment unavoidably results in the dependence of decay width on the parameter β, which is a parameter within the utilized SHO wave function.In this study, we capitalize on the knowledge acquired from the discussion on meson spectrum employing the MGI model in Section II to acquire the numerical spatial wave functions of the mesons under consideration.To accomplish this, we replace the form provided in Eq. (6.22) with a refined expression where the expansion coefficients c n represent the values of the corresponding eigenvectors, while l denotes the orbital angular momentum of the meson.By incorporating these modifications, we can effectively eliminate the associated uncertainties.To ensure proper normalization, the inclusion of the factor √ 4π is necessary, i.e., Besides, R nl (|p|) is the SHO wave function as (6.26)The parameter β = 0.5 GeV used in the above equation is consistent with Section II.In Eq. (6.26), we neglect the factor (−i) l since it does not affect the final results, but it does introduce a common factor of i to the P-wave final state weak transition form factors, which makes the form factors less concise.
Following the approach described in Refs.[201,202], we adopt the condition q + = 0.This implies that our form factor calculations are performed in the space-like region (q2 < 0), and therefore we need to extrapolate them to the time-like region (q 2 > 0).To perform the analytical continuations, we utilize the z-series parametrization [204], which has the form as [205] 1 where a 0 , a 1 , and a 2 are free parameters needed to be fitted in q 2 < 0 region, and z(q 2 ) is written as [204,206,207]  2 .m B c and m f are the masses of the B c meson and daughter meson, respectively.This parametrization are more convenient to reflect the character F B c P(S ) 1 (q 2 = 0) = F B c P(S ) 0 (q 2 = 0).Moreover, we have F (0) = a 0 .
To determine the values of the free parameters a 0 , a 1 , and a 2 2 , as stated in Eq. (6.27), we perform numerical calculations at 200 equally spaced points for each form factor, ranging from −20 GeV 2 to −0.1 GeV 2 , utilizing Eqs.(6.3)-(6.20).Subsequently, we fit the calculated points using Eq.(6.27).The fitted values of the free parameters, as well as F (0), F (q 2 max ), and the pole masses m pole , are compiled in Table X (XI, XII, XIII) for B c → P(V, S , A) transitions.Additionally, the q 2 dependence of the B c → P(V, S , A) transition form factors is illustrated in Fig. 4 (5, 6, 7).Particularly, in the B c → η c panel, we also present the LQCD's result [179].Our result is slightly larger than the LQCD in low q 2 region, but is anastomotic in large q 2 region.With the obtained weak transition form factors, we can future investigate the corresponding semileptonic decays.For the B c → Mℓν ℓ processes, the differential decay width can be obtained by where G F = 1.16637 × 10 −5 GeV −2 is the Fermi coupling constant, V CKM is the Cabibbo-Kobayashi-Maskawa (CKM) matrix element, λ(x, y, z) = x 2 + y 2 + z 2 − 2(xy + xz + yz) is the FIG. 4: The q 2 dependence of the weak transition form factors of the B c → B s , B, D, η c processes.In the B c → η c transition, the form factors F 1 and F 0 are also depicted as the results from LQCD, represented by a red circle and a blue square, respectively, as shown in Figure of Ref. [179].
If the discussed hadron has the near-threshold behaviour, the state should be a mixture of bc components and two meson components.The meson components should be important because the coupling is given by a strong mechanism, which gives the important decay widths, and they are enhanced by the denominators of the difference between the energy and the threshold energy position.Thus, the decays through the two meson components will be small and the results given in our work can be seen as an upper limit in most cases.

VII. SUMMARY
Although significant progress has been made in observing new hadronic states over the past two decades, the establishment of the B c meson family remains incomplete, with only B + c and B c (2S ) ± states listed in the PDG [46].With the upgrade being complete of the LHCb experiment in preparation for Run 3 and Run 4 of LHC, we have reason to believe this situation will be changed.Hence, the current work is timely and can provide valuable information for the experimental explorations of B c mesons.
According to the previous experience of the studies of hadron spectroscopy [60, 62, 63, 65-71, 210, 211], the importance of the unquenched effects has been realized step by step.The low mass puzzles involved in the X(3872), D s0 (2317), D s1 (2460), and Λ c (2940) can be well understood under this scenario.In this work, we provide a complete spectroscopy of the B c mesons under the unquenched picture.We first present the mass spectrum of the B c mesons, where the MGI model was applied to the concrete calculations, which can reflect the unquenched effects.The obtained mass spectrum of the B c mesons is valuable, but not sufficient for further experimental search for them.Therefore, in this work, we have investigated their various decay behaviors, including the two-body OZI-allowed strong decays, the dipion transitions betweeen B c mesons, the radiative decays, and some typical weak decays of B c (1 3 S 0 ).We must emphasize that the discussed decays of B c mesons are supported by the mass spectrum study, since we simultaneously obtain the information of the numerical spatial wave functions of these focused B c mesons associated with their masses, which are used as inputs for the calculations of the B c meson decays.This treatment avoids the parameter dependence in the decay studies.
With the accumulation of experimental data and enhancement of experimental capabilities, the investigations of hadron spectroscopy will enter a new stage.As an important part of the hadron family, the B c mesons become focal point, as our knowledge of the B c meson family is still infufficient.Facing this situation, we have reason to believe that it is full of challenges and opportunities.We expect that our experimental colleagues to seize this opportunity to continue to expand the realm of the observed hadrons.

FIG. 1 :
FIG. 1: (Color online) Mass spectrum of the B c mesons.Here, the red (left) lines and the blue (right) lines are our obtained results and the results from the GI model [52], respectively, while the short lines denote the thresholds of the B ( * ) (s) D ( * ) (s) channels.The masses of the mesons are given in units of MeV.

FIG. 2 :
FIG.2:The obtained spatial wave functions of the concerned B c mesons.Here, these B c mesons are classified by the notation n 2S +1 L J .
S 1 ) state are BD and B * D; (ii) The B c (4 3 S 1 ) state mainly decays into the B * D, BD * , and B * D * channels, and the BD also has the sizable contribution to the total width; (iii) The main decay modes of the B c (5 3 S 1 ) state include the B * D, BD * , B * D * , BD(1P 1 ), and BD(1P ′ 1 ) channels, while the dominant decay channel is the B * D * with branching ratio 70.1%.
for the B c (3 3 D 1 ), B c (3D 2 ), B c (3D ′ 2 ), and B c (3 3 D 3 ) states, respectively.For the 3D states of the B c meson, the decay widths of some channels are significantly suppressed.Such as the B s D * s channel for the B c (3 3 D 1 ) state and the B s D s channel for the B c (3 3 D 3 ) state, whose decay widths are 4.7×10 −3 MeV and 2.3 × 10 −4 MeV, respectively.

TABLE I :
The masses (in units of MeV) of the experimental values and our obtained results of the heavy flavor mesons in this work.

TABLE IV :
Partial widths and the corresponding branching ratios for these two-body OZI-allowed strong decays of the S -wave and P-wave B c mesons.Here, the decay widths of the discussed mesons are given in units of MeV.

TABLE V :
Partial widths and the corresponding branching ratios for these two-body OZI-allowed strong decays of the D-wave and F-wave B c states.Here, the decay width of the meson is given in units of MeV.

TABLE VI :
Decay rates of dipion transitions between B c states.

TABLE VIII :
Partial widths of the M1 transitions for the S -wave B c states compared with other theoretical works.Here, the width is in units of eV.

TABLE IX :
Magnetic moments of the S /P/D/F-wave B c mesons.Here, the magnetic moment of the hadron is in units of the nuclear magneton µ N = e/2m p .

TABLE X :
The form factors of the B c → B s , B, D, η c transitions in CLFQM.

TABLE XII :
The form factors of the B c → B s0 , B 0 , D * 0 , χ c0 transitions in CLFQM.

TABLE XIV :
The branching ratios of the B c → B s (B, D, η c )ℓν ℓ decays.Additionally, we compare our findings with the branching ratios reported in other theoretical studies.

TABLE XV :
The branching ratios and the Γ L /Γ T of the B c → B * s (B * , D * , J/ψ, ψ(2S ))ℓν ℓ decays.Additionally, we compare our findings with the branching ratios reported in other theoretical studies.