Rotation and Vibration in Tetraquarks

A novel approach is introduced for obtaining precise solutions of the pairing Hamiltonian for tetraquarks, which utilizes an algebraic technique in infinite dimensions. The parameters involved in the transition phase are calibrated based on potential tetraquark candidates derived from phenomenology. Our investigation shows that the rotation and vibration transitional theory delivers a reasonable agreement with other works for heavy tetraquarks compared to other methods. To illustrate the concept, we compute the spectra of several tetraquarks, namely charm, bottom, bottom–charm and open charm and bottom systems, and contrast them with those of other particles.


I. INTRODUCTION
Scientists have detected a new particle, dubbed X(2900), by analyzing all the data collected so far by the LHCb experiment at CERN's Large Hadron Collider [1,2].This experiment is renowned for discovering exotic quark combinations, which help scientists study the strong force, one of the four fundamental forces in the universe.LHCb has identified several tetraquarks, made up of four quarks (or two quarks and two antiquarks), including the latest discovery of an entirely new type of tetraquark with a mass of 2.9 GeV/c 2 , which has only one charm quark.While scientists predicted this particle's existence in 1964, it is the first observed instance of a tetraquark with only one charm quark.Quarks cannot exist independently; they form composite particles, such as mesons (a quark and an antiquark) or baryons (three quarks or three antiquarks), like the proton.The LHCb detector located at the LHC focuses on studying B mesons, which are composed of a bottom or an anti-bottom quark.These mesons quickly decay into lighter particles shortly after being produced in proton-proton collisions at the LHC.Tetraquarks are believed to be pairs of distinct mesons that are temporarily bound together like a "molecule," according to some theoretical models, while others view them as a single cohesive unit of four particles.Identifying and measuring the properties of new kinds of tetraquarks, such as their quantum spin and parity, will provide a better understanding of these strange inhabitants of the subatomic realm.The recently discovered particle, called X(2900), contains an anticharm, an up, a down, and an antistrange quark (cud s]) and is considered the first open-charm tetraquark, as all previous tetraquark-like states observed by LHCb had a charm-anticharm pair, resulting in a net-zero "charm flavour."[1,2] Pairing interactions between fermionic or bosonic systems are common in many physical contexts such as Bose-Einstein Condensation and Superfluidity, airing correlations in nuclei: from microscopic to macroscopic models, high-temperature superconductors, [3][4][5][6][7][8].One example of the application of algebraic methods in hadron physics is the use of such interactions.[9][10][11][12][13][14].We establish explicit extensions of duality relations that relate the Hamiltonians and basis classification schemes associated with number-conserving unitary and number-nonconserving quasispin algebras for four-level pairing interactions.The Hamiltonian of the model can be defined using a linear combination of first-and second-order Casimir operators when one-and two-body interactions are present.The four-level pairing model describes a finite system that undergoes a second-order quantum phase transition between the rotation and vibration limits.Recently, we utilized the interacting boson approximation proposed by Arima and Iachello [15,16] to calculate wave functions in an interacting sl many-body boson system [17].It is important to note that, in general, the building blocks of the boson system are associated with both s and l bosons for single and quadrupole angular momentum.The bosonic pairing systems exhibit similarities in their Lie algebraic properties, but the differences are significant in terms of the irreducible representations (irreps) that the eigenstates transform under, which play a critical role in defining the system's spectroscopy.Finite pairing systems can be described by two complementary algebraic formulations: (1) a unitary algebra consisting of bilinear products of a creation and annihilation operator, and (2) a quasispin algebra that uses creation and annihilation operators for time-reversed pairs of particles [18][19][20][21][22][23][24].
Tetraquarks are exotic hadrons composed of four quarks that can include two quarks and two antiquarks or four quarks of the same flavor.Despite being first proposed in the 1960s, their existence was only confirmed in 2013 by the Large Hadron Collider experiments.In recent years, the study of tetraquarks has gained increasing interest due to their unique properties and potential implications in particle physics [13].In this context, we propose to apply an algebraic framework to investigate the properties of heavy tetraquarks [QQ][ Q Q].Our approach is based on the S U(1, 1) algebraic technique [19,25,26] and extends the sl boson system.We will derive a new solvable model for hadron physics that takes into account the vector quark pairing strengths and examine the mass spectra of tetraquarks.
In recent years, there has been a growing interest in exploring the properties of fully-heavy tetraquarks.Theoretically, several models have been proposed to describe these states, including the diquark-antidiquark model, the chromomagnetic interaction model [27].On the experimental side, various searches have been performed to identify fully-heavy tetraquarks in high-energy experiments.For instance, the LHCb collaboration searched for deeply bound bb bb tetraquark states, but no significant excess was found in the µ + µ − Υ(1S ) invariant-mass distribution [28].However, the CMS experiment reported a potential candidate of a fully bottom tetraquark T 4b = [bb][ bb ] around 18-19 GeV [29].Moreover, the LHCb collaboration has recently reported the observation of a narrow peak and a broad structure in the J/ψ-pair invariant mass spectrum, which could originate from hadron states consisting of four charm quarks [30].These experimental results provide valuable information for further theoretical investigations of fully-heavy tetraquarks.

II. THEORETICAL METHOD
In the context of describing a tetraquark system, diquark clusters play an important role.It is suggested that a tetraquark system, denoted by T = Q 1 Q 2 Q3 Q4 , consists of two point-like diquarks.To consider multi-level pairing in this context, we extend the interacting boson model using algebraic solutions of an sl-boson system [17].The dynamical symmetry group in this case is generated by s and l operators, where l represents the configuration of the multiquark states.In the Vibron Model, scalar s-bosons with spin and parity l π = 0 + and vector l-bosons with spin and parity l π = 1 − represent elementary spatial excitations.The generators in the finite-dimensional S U(1, 1) algebra satisfy the following commutation relations.
[S 0 (l), S ± (l)] = ±S ± (l) , We can use the S U(1, 1) algebra to describe the rotation and vibration transitional Hamiltonian of the It is worth mentioning that the quasi-spin algebras have been extensively discussed in previous studies, such as Refs.[17,23].
Taking into account the generators of the S U l (1, 1)-algebra for tetraquarks given by Eqs.(1a) and(1b), we can express the relevant quantities as linear combinations of these generators.
The infinite dimensional S U(1, 1) Lie algebra is defined by The real-valued control parameters c Q and c Q play a crucial role in determining the properties of tetraquarks.Specifically, l 1 and l 2 correspond to the first and second tetraquarks, respectively, while l3 and l4 correspond to the third and fourth tetraquarks.
Additionally, the integer n can take on values of 1, 2, 3, and so on.
To ensure that the fully-heavy tetraquarks satisfy the correct properties, we impose the condition S − (l)|lw = 0 on the lowest weight state.The state |lw can be defined as follows: where Hence, we have The system shows vibrational and rotational transitions due to continuous variations of the pairing strengths, c l , in the closed interval [0, 1].The all-heavy tetraquark pairing model undergoes a quantum phase transition.The vibration limit is reached when In our analysis, we obtained diverse values for the control parameters, c Q i and c Qi , in the interval [0, 1] with i = 1, . . ., 4, between the two limits.
The Hamiltonian of the heavy tetraquark pairing model is expressed in terms of the Casimir operators Ĉ2 using branching chains.The first two terms of the Hamiltonian, S + 0 S − 0 and S 0 1 , are associated with the S U(1, 1) algebra, while the remaining terms are constant in terms of the Casimir operators.In the duality relation for tetraquarks, the irreducible representations simplify the quasi-spin algebra chains (4a) and (4b), and the labels for the chains are related via the duality relations.The Hamiltonian for the heavy tetraquark pairing model is derived by utilizing the generators of the S U(1, 1) algebra.However, the pairing models of multi-level are also characterized by an overlaid U(n 1 + n 2 + . ..) algebraic structure with this branching: So, we can define the Hamiltonian with where g, α, β, γ 1 , γ 2 , γ 3 , and γ are real-valued parameters.
To find the non-zero energy eigenstates with k-pairs, we exploit a Fourier Laurent expansion of the eigenstates of Hamiltonians which contain dependences on several quantities in terms of unknown c-number parameters x i , and thus eigenvectors of the Hamiltonian for excitations can be written as and The coefficients x i are determined through the following set of equations In the pursuit of finding exact solutions for a spin-spin interaction system, Gaudin utilized a similar structure [56] as an ansatz, which has now been verified as a consistent operator form in constructing the Bethe ansatz wavefunction for the present tetraquark system.To obtain the energy spectra, the Bethe ansatz equation (BAE), a non-linear equation, is employed for a k-pair excitation.The quantum number k-pair excitation pertains to the overall number of bosons N and is linked to seniority numbers, specifically the quantum number ν l of S O(2l + 1).As per equation ( 4), the allowed seniority numbers for a fixed ν l include n l = ν l , ν l + 2, ν l + 4, and so on.This information is well-established in the field.Our approach to calculating the masses of heavy tetraquarks follows the procedure outlined in Ref. [19].To account for the bosonic nature of the excitations (vibrations and rotations), we use the totally symmetric representation (7) and define the boson number as the total number of vibrational states in the representation [N].
Pairing in tetraquarks is an interesting phenomenon that affects their rotational and vibrational behavior.The quantum phase transition occurs between the vibrational and rotational limits in the fully-heavy tetraquark pairing model, and the quark (antiquark) configuration can undergo vibrations and rotations described by the quantum numbers ν Q i , ν Qi , and J.While we will not consider bending and twisting in this analysis due to their higher mass requirements, we must account for the internal degrees of freedom of quarks and antiquarks.To address this complication, we apply the method of pairing strengths, following Refs.[13,17].This scheme illustrates the stringlike configuration of the tetraquark and the vibration-rotation pattern we aim to identify.Within the two-quark configuration, we must follow the operator Q with Q, but since we are dealing with tetraquarks, we could also have combinations of QQ and Q Q.To determine the appropriate rotation-vibration pattern, we need to ensure that the pairing number We can optimize the control parameters to find the exact symmetry of vibration and rotation that gives us the desired pairing number.By understanding the pairing behavior in tetraquarks, we can better understand their physical properties and potentially make predictions for future experiments.In summary, our work builds on previous research in the field, but we introduce new ideas related to pairing in tetraquarks and optimize control parameters to identify the appropriate symmetry of vibration and rotation.

III. RESULTS
The determination of tetraquark mass in the diquark-anti-diquark pairing model requires solving the eigenvalue problem of Eq. (7).However, in addition to the spins of diquark and antidiquark clusters, the J PC quantum numbers that define a tetraquark state also include the total spin, spatial inversion symmetry, and charge conjugation of the system.Recent studies have shown that the J PC quantum numbers of a Q 1 Q 2 Q3 Q4 system can be 0 ++ , 1 +− , and 2 ++ , as discussed in Ref. [62].These quantum labels are essential for characterizing the properties of the tetraquark system.The total spin of the tetraquark is determined by the combination of the spins of diquark and antidiquark clusters, and it affects the tetraquark's stability and decay properties.Spatial inversion symmetry is related to the tetraquark's mirror image, and it determines whether the system is symmetric or asymmetric with respect to spatial inversion.Charge conjugation, on the other hand, is related to the transformation of particles to their corresponding antiparticles and is a fundamental symmetry of the strong interaction.Understanding the impact of total spin, spatial inversion symmetry, and charge conjugation on tetraquark states is crucial for predicting their properties and behavior.By considering these quantum numbers, we can gain insights into the tetraquark's internal structure and its interactions with other particles.This knowledge is essential for advancing our understanding of the strong interaction and the behavior of exotic hadrons.For scalar, vector and tensor systems, we have: 1. Two states for the scalar system: 2. Three states for the vector system: Charge conjugation is a fundamental symmetry of the strong interaction that transforms particles into their corresponding antiparticles.This symmetry leads to different configurations in which |A and |B can interchange, while |C remains odd.
In the J P = 1 + configuration, we have one C-even and two C-odd states.This arrangement plays a crucial role in determining the properties and behavior of the system.Understanding the implications of these configurations is essential for predicting the tetraquark's stability and decay properties.
When considering tetraquarks, it is essential to choose appropriate values for the spin of the quark-antiquark pairs.In particular, the selection of spin states can impact the overall properties and behavior of the system.
In the case of a tetraquark composed of Q 1 , Q 2 , Q3 , and Q4 , the appropriate spin states depend on the charge conjugation of the system.Specifically, when C = +, the only allowed state is one where Q 1 Q3 has a spin of S Q 1 Q3 = 1.
3. One state for the tensor system: where this state has also S Q 1 Q3 = 1.

A. The charm system
The pairing tetraquark model considers two phases, rigid and non-rigid, which correspond to rotation and vibration symmetries, respectively.While both phases are idealized situations, they must coexist in reality, resulting in the emergence of vibrational-rotational modes in the transitional region.The parameters in this region are known as the phase parameters, where c Q i = 1 with i = 1, . . ., 4 corresponds to the rotational mode and c Q i = 0 corresponds to the vibrational mode.The mass spectrum of the pairing tetraquark model can be calculated with fixed phase parameters, and the transitional spectra from one phase to another can be obtained by adjusting the phase parameters within the closed interval [0, 1].
Our numerical values for the coefficients are c Q 1 = 0.92, c Q 2 = 1, and c Q3 = c Q4 = 0.These values yield the following mass values: for the T 4c tetraquark system.As shown in the (Fig. 1) , we overall calculate the mass for T 4c tetraquark system.

B. The bottom system
The scenario presented here bears some resemblance to the earlier case.However, this time, according to the transitional theory, the extraction phase coefficients must be computed with regard to the meson-meson thresholds, such as η b (1S )η b (1S ) and Υ(1S )Υ(1S ) for J PC = 0 ++ , η b (1S )Υ(1S ) for J PC = 1 +− , and Υ(1S )Υ(1S ) for J PC = 2 ++ .Our numerical values for the coefficients are c Q 1 = 0.97, c Q 2 = 1, c Q3 = 1, and c Q4 = 0.These values yield the following mass values: for the T 4b tetraquark system.As shown in the (Fig. 2) , we overall calculate the mass for T 4b tetraquark system.
FIG. 1: The predicted mass spectrum of the T 4c tetraquarks.All spectroscopies are in GeV.
FIG. 2: The predicted mass spectrum of the T 4b tetraquarks.All spectroscopies are in GeV.

C. The bottom-charm system
This study also considers the T 2bc = [bc][ bc] tetraquark structure, which combines c quarks with b quarks.Here, the [bc] diquark spin may be either 0 or 1, allowing for the possibility of all states analyzed in the previous section.Once again, the best method for extracting the control parameters in T 2bc tetraquarks, based on the transitional theory, is to utilize the corresponding meson-meson families.This method yields the following values: c Q 1 = c Q 2 = 1, and c Q3 = c Q4 = 0.The computed masses can be classified into the following categories: (i) The J PC = 0 ++ contains two scalar states with masses FIG. 3: The predicted mass spectrum of the T 2bc tetraquarks.All spectroscopies are in GeV.
(ii) The J PC = 1 +− contains two states with masses (iii) The J PC = 1 ++ contains one state with mass (iv) The J PC = 2 ++ contains one state with mass As shown in the figure (Fig. 3) , we overall calculate the mass for T 2bc tetraquark system.

D. The open charm and bottom system
In the most recent research, a thorough investigation was conducted on open charm (OC) and bottom (OB) tetraquarks comprising bottom and charm quarks, including cq q q, cq s q, cs s q, and cs s s for charm, and bq q q, bq s q, bs s q, and bs s s for bottom.The outcomes for the masses of these tetraquraks are summarized in Figs. 1, and the available experimental data is compared in Table .1. Based on the method used in this study, the resulting values are demonstrated in the Caption of Fig. 1.Our numerical values are 0 for open charms cq q q, cq s q, cs s q, and cs s s, respectively.In contrast, numerical values are 37 and c Q4 = 0 for open bottom bq q q, bq s q, bs s q, and bs s s, respectively.
As shown in the (Fig. 4) , finally we calculate the mass for open charm and bottom tetraquark systems.

IV. DISCUSSION
The calculation process involves a fixed set of Hamiltonian parameters while allowing the phase parameters to fluctuate during the transition.In Ref. [14], the authors demonstrated that the boson number's quantum value could be obtained by taking the N → ∞ limit.It was found that taking N to be a large number was sufficient to account for all known and unknown states up to the maximum value of the quantum number of the angular momentum and other relevant quantum numbers for the applications.In this current research, we are utilizing the same approach as in Ref. [14], setting N = 100 to ensure that all states up to the maximum quantum number are considered.
The Hamiltonian's pattern is comparable to the O(4) restriction proposed in mesons, where the control parameter is set to 1.Our study indicates that the control parameters c Q3 and c Q4 cannot be considered as 1 in the presence of heavy antiquarks, except for T 4b tetraquarks.This is due to the fact that the T 4b tetraquark's mass is two to three times more substantial than that of T 2bc and T 4c tetraquarks.For heavy mass tetraquarks, pairing strength plays a significant role, as evidenced by the larger c Q 1 value in the T 2bc case than in the T 4c case.The same reasoning applies to the open and bottom tetraquark system.The Hamiltonian's parameters for the discussed structures are presented in the Figures' captions.During the transition phase, we set α to be 1.5.As the pairing model's vibrational-rotational transition is a second-order quantum phase transition, the masses of the wave functions in the tetraquark's vibrational model are smooth with respect to parameter changes.This enables us to determine them in the transition region.
In Table I, we present the difference between the calculated masses of the tetraquarks and the threshold for meson-pairing.The values of ∆ represent the difference between the tetraquark mass M tetra and its lowest meson-meson threshold E th .If ∆ is negative, it implies that the tetraquark state lies below the fall-apart decay threshold and should therefore be stable.On the other hand, a state with a small positive ∆ could be observed as a resonance due to suppression by phase space.The states with high positive ∆ values are considered broad and difficult to detect in experimental analyses.
Our investigation indicates that slightly deviating the control parameter c Q 1 from 1 is more suitable for determining the tetraquark masses, especially for the extensive T 2bc families.Additionally, for T 4b = [bb][ bb ] states, the dominant contribution comes from the pairing of c Q3 and c Q4 quarks, indicating that phase parameters for Q3 and Q4 quarks become significant in computing tetraquark masses at high energy, around 18 − 19, GeV.On the other hand, at low energy, there is a competition between Q 1 and Q 2 .
The energy spectra of the fully-heavy tetraquarks under study, where c Q i values are in the range of 0.9 − 1.0, can be attributed to a rotational phase, based on the definition mentioned above.It is worth noting that a change of ±15% in all coefficients results in a maximum variation of 30%, 23%, and 17% in the masses of the T 4c , T 4b , and T 2bc tetraquark systems, respectively.However, it should be emphasized that the masses of the remaining tetraquark states undergo lesser modifications.
Our study on the vibrational and rotational transitions in open and bottom tetraquarks suggests that a slight deviation of the control parameter from the vibration limit is more appropriate for determining the tetraquark masses, particularly for the open bottom families.This finding emphasizes the importance of choosing the correct control parameter for accurate mass spectroscopy.
Furthermore, we observed that in open bottom tetraquark states, the dominant contribution comes from the pairing of the third quarks, indicating that the phase parameters for these quarks play a crucial role in computing tetraquark masses at high energies around 5-6 GeV.In contrast, for open charm tetraquarks, there is a competition between the first and second quarks.These observations suggest that the pairing of quarks and the interplay of their phase parameters can significantly impact the Our findings provide insights into the properties of fully-heavy tetraquarks and highlight the importance of considering the multiquark dynamics for a more comprehensive understanding of hadron spectroscopy.Further studies on other types of multiquarks and the inclusion of other degrees of freedom, such as pentaquarks and hexaquarks, may lead to a better understanding of the underlying physics and shed light on the nature of hadronic matter.
The comparison between our results obtained from the pairing model and the predictions of previous theoretical calculations are presented in Tables II and IV.Our findings show that the pairing model provides reasonable agreement with the other works, implying that it has the potential to play a crucial role in predicting fully-heavy tetraquark mesons.However, future work could aim to further improve the understanding of multiquark dynamics by including the large-N limit of the pure pairing Hamiltonian.

V. SUMMARY
The aim of this study was to investigate the mass spectra of tetraquarks in the transition region between vibration and rotation using the algebraic framework.To achieve this, a solvable extended transitional Hamiltonian based on S U(1, 1) algebra is proposed, which can describe both partial high energy states and quantum phase transition.The extracted mass spectra of various tetraquarks were in agreement with previous research and other theoretical approaches.However, it is important to consider other degrees of freedom, such as penta or hexa quarks, in future studies.Furthermore, the solvable technique introduced in this work could potentially be applied to diagonalize more complex multiquark systems.This approach is currently being applied to investigate other types of multiquarks in the following manuscript.

FIG. 4 :
FIG. 4: The predicted mass spectrum of the open charm and bottom tetraquarks.All spectroscopies are in GeV.

TABLE I :
Masses of fully-heavy tetraquark systems as computed within the theoretical framework presented herein.The meson-meson threshold is E th , and ∆ = M − E th represents the energy distance of the tetraquark with respected its lowest meson-pair threshold.The notation s and a indicates scalar and axial-vector diquarks.