Doubly-heavy tetraquark at finite temperature in a holographic model

In this paper, we employ gauge/gravity duality to investigate the string breaking and melting of doubly-heavy tetraquark that includes two heavy quarks and two light antiquarks in a holographic model at finite temperature. Firstly, we investigate four configurations of $\rm{QQ\bar{q}\bar{q}}$ in the confined phase and consider different separation distances of the heavy quarks at varying temperatures. At high temperature, $\rm{QQ\bar{q}\bar{q}}$ melts at certain distances and the confined quarks are released. As the temperature continues to increase, some configurations of doubly-heavy tetraquark can not exist. Furthermore, we investigate three decay modes of $\rm{QQ\bar{q}\bar{q}}$ and compare the potential energy of $\rm{QQ\bar{q}\bar{q}}$ with that of $\rm{QQq}$ at finite temperature .


I. INTRODUCTION
The gauge/gravity duality has been widely recognized as a fundamental feature of quantum gravity, and extensive research has been carried out in this field over the past two decades, leading to many important findings [1].Holographic QCD offers a novel approach to study and compute the properties of various physical phenomena in QCD.Studying heavy quarkonium is beneficial for understanding the properties of quark-gluon plasma(QGP), as well as validate our understanding of the interactions between hadrons and fundamental particles [2,3].In holographic QCD, we can study the interaction between quarks and antiquarks by placing a pair of them inside a bulk.Then, we can utilize the tools of string theory to calculate the interaction between the two quarks and their corresponding potential energy.
The string breaking phenomenon is a result of the nonperturbative effects of QCD.Up to now, lattice QCD has been effectively utilized to investigate this phenomenon, albeit limited to meson modes at zero temperature and zero chemical potential [4].As is widely acknowledged, the environment in which tetraquark are situated can often be complex.Therefore, introducing high-temperature factors into the model to study the behavior of four quark potentials at elevated temperatures can aid in better understanding the behavior of tetraquark under extreme conditions.In previous studies, it has been discovered that under sufficiently high temperatures, thermal excitations produce a plasma of quarks and gluons [5].Subsequently, they discussed the potential at finite temperatures [5][6][7].The deformed AdS 5 model [8][9][10] and Einstein-Maxwell-Dilaton model are employed to compute the quark-antiquark potential [11][12][13][14][15][16][17][18][19][20][21][22][23][24][25][26][27][28] in these researches.
Recently, at the Large Hadron Collider, the LHCb collaboration observed a hadron state that contains four quarks [29,30].This tetraquark contains two charm quarks, a ū and a d quark, with a mass of about 3875 MeV/c 2 .This finding has renewed interest in studying the tetraquark theory.It is worth noting that lattice gauge theory is still a fundamental tool for studying non-perturbative phenomena in QCD, but research results on the potential of QQqq are limited [31][32][33][34][35][36].
The tetraquark model used in this paper was proposed by Andreev [37].This model assumes that the heavy quarks are significantly heavier than the typical energy scale of the system, allowing us to treat them as static.The interaction potential between the quarks is determined by their relative separation.The main reason for studying this model is that its results on both quark-antiquark and tetraquark potentials are consistent with lattice calculations and QCD phenomenology [4,38].At zero temperature, the timelike Wilson loop is realized by a U-shaped macroscopic string for any interquark separation.At finite temperature, we have seen that the string configuration is either a pair of straight strings or U-shaped [5].Our technique for extracting the potential of QQqq is similar to the one used in lattice QCD.We extract the potential from the expectation value of the QQqq Wilson loop, W QQqq (R, T ).The QQqq Wilson loop consists of heavy quark paths and light quark propagators.The separation distance between heavy quark pairs and the potential energy relationship of the model were calculated and analyzed at finite temperature [39][40][41][42][43][44][45][46][47].
Subsequently, three decay modes of tetraquark were studied to determine the most probable decay mode that occurs at high temperatures and to compare it with the decay of three quarks.
Sec. II provides a brief review of the theoretical foundation of the model and establishes a framework for studying QQqq at finite temperature.Sec.III involves numerical solutions for the energy and separation distance between heavy quarks for these configurations at different temperatures.This is followed by a discussion of the results in Sec.IV.In this section, we will also discuss the melting of QQqq strings at finite temperatures by analyzing potential energy trends.Later in Sec.V, the three types of decay modes in the QQqq model will be discussed and compared with the decay modes in the QQq model.Finally, Sec.VI presents the summary and conclusion of this paper.

II. PRELIMINARIES
In this paper, we extend the study of the potential for doubly-heavy tetraquarks from zero temperature [37] to finite temperatures.To begin our discussion on the potential of QQqq at finite temperatures, let us first review the specific holographic model utilized in this paper.Following [37,48,49], the background metric at finite temperatures is given by: ab dω a dω b , such model is a deformation of the Euclidean AdS 5 space of radius R, with a deformation parameter s [50].In the five-dimensional compact space (sphere) characterized by the black-ening factor f with coordinates ω a and f (r), the Nambu-Goto action of a string is expressed as here, γ represents an induced metric on the string world-sheet with a Euclidean signature, while α ′ is a parameter associated with the string.For the AdS 5 space, we assume that the blackening factor f takes the form f (r) = 1 − r r h

4
, where f (0) = 1 at the boundary and f (r h ) = 0 at the horizon.The Hawking temperature, which is consistent with the temperature of the dual gauge theory, can be defined as T = 1 4π |∂ r f | r=r h .Then, we consider baryon vertices which are string junctions [51].According to the AdS/CFT correspondence, they are represented by a five-brane wrapped around the internal space X at the point where three strings intersect and appear to be joined together in five dimensions [52,53].Correspondingly, the antibaryon vertex is represented by an antibrane in the AdS/CFT correspondence.At leading order α ′ , the brane's action is (6) , where T 5 represents the brane tension and ξ i denotes the worldvolume coordinates.If we choose static specifications ξ 0 = t and ξ a = θ a , where θ a represents the coordinates on X, then the resulting action is as follows Here, T v is a dimensionless parameter defined by T v = T 5 Rvol(X), where vol(X) represents the volume of X, and it serves as a free parameter.Additionally, we have the same action for both baryon and antibaryon vertices S v = S v at finite temperature.
Finally, we consider the light quark located at the end of the string, a scalar field that is coupled to the boundary of the worldsheet via the open string tachyon S q = dτ eT [54].Here, the integral is over a world-sheet boundary parameterized by τ ,and e is a boundary metric.Assuming a constant background given by T(x, r) = T 0 and worldsheets with boundaries along lines in the t direction.Thus, the action can be written as Here, m =RT 0 represents the mass of a point particle at rest, with T 0 as its mass.Therefore, the given action describes the behavior of a point particle with mass T 0 at rest.It should be noted that the same action also describes the behavior of light antiquarks located at string endpoints, and hence S q = S q .In this model, the interaction of quarks is described by the string tension, which is consistent with Ref. [3].
The model parameters are selected as follows: g = R 2 2πα ′ , k = τv 3 g and n = m g .All our parameters are consistent with Refs.[4,37,48,49] based on the lattice QCD.The value of s is fixed from the slope of the Regge trajectory of ρ(n) mesons in the soft wall model with the geometry Eq. ( 1).This gives s = 0.45GeV 2 .Then, fitting the value of the string tension σ to its value in [39] gives g = 0.176.According to the gauge/string duality g is related to the 't Hooft coupling.Next, the parameter m is adjusted to reproduce the lattice result for the QQ [39], it gives m = 0.538.The parameter n is expressed in terms of the parameters of as n = m g [4].For fixed the value of k, one should keep in mind two things.First, the value of k can be adjusted to fit the lattice data for the three-quark potential, as is done in [38] for pure SU (3) gauge theory.Unfortunately, at the moment, there are no lattice data available for QCD with two light quarks.Second, the range of allowed values for k is limited to − e 3 15 to − 1 4 e 1 4 .We take k = − 1 4 e 1 4 simply because it yields an exact solution to equations (20).Finally, the parameters of this article are as follows: 4 , and c = 0.623 GeV in this model from Ref. [37] at zero temperature.No other extra parameters are introduced in this article.

III. THE CONNECTED STRING CONFIGURATION OF QQqq
A. Small L The configuration of QQqq is illustrated in Fig. 1.The total action for this system is given by the sum of the Nambu-Goto actions, as well as the actions associated with the vertices and antiquarks.
We choose the static gauge ξ 1 = t and ξ 2 = r for the Nambu-Goto actions, the boundary conditions for x(r) become x A static string configuration at small heavy-quark separation distance.The heavy quarks Q are located on the boundary, whereas the light antiquarks q, baryon vertex V, and antibaryon vertex V are situated in the bulk of the five-dimensional space at r q, r v and r v, respectively.r h represents the position of the black-hole horizon, while r w indicates the position of a dynamic wall in the confined phase.
Considering the boundary conditions, we get the total action Here ∂ r x = ∂x ∂r , T = T 0 dt and the straight strings are located at x = 0.For string (1) and (2), which correspond to the first term in (7), we can derive the equation of motion (EoM) for x(r) using the Euler-Lagrange equation.Thus, it is found that I is a constant.We have ∂ r x = cot α when r = r v , and I can be expressed as Then ∂ r x can be obtained: Using equation (10), we can obtain an expression for the separation distance L, Next, we calculate the potential energy of doubly-heavy tetraquark.The energy of string (1) can be got from the first item in the equation ( 7), This expression is not well-defined, because the integral diverges at r = 0.So we should subtract the divergent term here c is a normalization constant.We can fix the constant by fitting the lattice results.For the heavy quarkonium, we can take 2c to fit the potential of heavy quarkonium with lattice.
The choice of the normalization constant c for the energy of a single baryon configuration is equal to 3c. [55] String (2) is calculated in the same way as string (1), and thus String (3) is described by the second term in equation (7), and represents a straight string stretched between the baryon vertex and antibaryon vertex.The energy can be calculated as String (4) and string (5) are both described by the third term in equation (7), and represent straight strings.
Then, we can get the energy of QQq q for this configuration.
It can be observed from equation ( 16) that energy is a function of r v , α, and r h .To obtain the energy of this configuration, we will solve the position of the light antiquark.One crucial condition for equilibrium is that the net forces exerted on the vertices and antiquarks must cancel out.According to the model, the force (shown in Fig. 1) balance equation in the r direction at r q is given by: Here e i is the string tension [56].By varying the action with respect to r q, we get the force = gw (r q) (0, −1) on the antiquark.Hence, the equation (17) becomes r q (the position of the light antiquark) is only a function of r h .This equation gives us the position r q at a fixed temperature.Then, at r v, the force balance equation is Here f v is the force on the antibaryon vertex, and each force is determined by Then, the force balance equation changes At fixed temperature, r v can be determined via the above equation.At r v , the force balance equation is here the force on the baryon vertex is , and the string tensions at the r v are e 3 = gw (r v ) (0, 1), , . α is the angle shown in Fig. 1.Putting these forces into equations (21), the force balance equation becomes By taking the variation of the action with respect to r v , we can deduce the relationship between r v and α.

B. Slightly larger L
The second configuration is shown in Fig. 2, and we consider the second configuration as one in which a string (3) contracts to a single point.The total action is now expressed by We choose the same static gauge as before and the boundary conditions are The separation distance of the slightly larger L is still determined by equation (11), which is the same as that used for the small L. Compared to the first configuration, the string (3) in the second configuration contracts to a point.Therefore, we need to consider other string tensions that satisfy equations (13) − (15).Then, the energy of the slightly larger L can be expressed as x e q q f q f q FIG. 2. A static configuration is formed by a slightly larger separation distance of a heavy-quark pair.The baryon vertex and the antibaryon vertex are in the same position.The force acting on the point is indicated by the black arrow.
As before, we will now proceed with solving the force balance equation.The location of r q can be determined using equations (18).Then, the force balance equation at the point Each force is determined by Then, the equation (26) becomes

C. Intermediate L
The third configuration, as shown in Fig. 3, is characterized by the compression of points r q, r v and r v into a single location.The configuration of the total action is given by Choosing the static gauge in the Nambu-Goto actions as before, then the boundary conditions can be obtained Equation (11) gives us the separation distance in this configuration.To calculate the energies, we only need to consider the range from r = 0 to r = r v , as specified by equations (13).
Naturally, the energy of the configuration is The force balance equation at the point x Each force is determined by Here r q = r v = r v , then, the equation (31) becomes D. Large L Fig. 4 illustrates the fourth configuration.The total action for this configuration, denoted as intermediate L, is given by equation (28).We choose another static gauge ξ 1 = t and ξ 2 = x in the Nambu-Goto actions and the boundary conditions are r (1) Then, the total action becomes The action for string (1) is given by the first term in equation (34).The subsequent step is to compute the first integral.
I is a constant.At the r 0 and r v points, we can obtain, respectively Here ∂ x r can be obtained from equation (36),( 37) The large L configuration has a turning point at r 0 , so the distance between heavy quarks is calculated in two parts.Then, the separation distance can be expressed as Here r ′ denotes ∂ x r.By substituting equation (38) into equation (39), the separation distance can be obtained.The energy of string (1) was calculated in two parts, which is similar to the separation distance.Therefore, the energy of string (1) is given by Same as the case of small L, we also need to subtract the divergent term here.Thus, equation (40) can be expressed as String ( 2) is calculated in the same way as string (1).Therefore, the total energy of the configuration is The force balance equation is the same as that for intermediate L, and the expressions of force are also the same as intermediate L. Equation (37) represents the functional relationship between r v and α when the temperature is fixed.We can first solve equations (32) and (37), and then substitute the numerical values into equations ( 39) and ( 42) to obtain the solutions for L and E QQq q.For α, at small separate distance of heavy quarks, as r v increases, the distance between the heavy quarks increases, and the curvature of the string between heavy quarks increases, causing α to decrease.When the string reaches α = 0 and r v increases further, the string will become an M-shape, at which point α becomes negative.

A. T=0.08 GeV
At a low temperature of 0.08 GeV, the configurations of QQqq are confined.Below the black-hole horizon, there exists a dynamic wall that prevents the string from crossing.In this part, we investigate four configurations of QQqq at this temperature.
First, we will give a discussion about the dynamic wall at finite temperature.Different from the quark-antiquark pair, the configuration of QQqq shows a "M" shape at large separation distance.Thus, the maximum value of r 0 at infinite separation distance gives the position of dynamic wall.By calculating the position of the dynamic wall, we obtain r w = 1.40GeV −1 at T = 0 and r w = 1.45GeV −1 at T = 0.08 GeV as shown in Fig. 5.
The increase of temperature leads to an increase of the position of dynamic wall in the r direction.For small L, we use equation (18) to calculate the position of the antiquark.The result shows that at a temperature of 0.08 GeV, r q = 1.13 GeV −1 .Subsequently, the antibaryon vertex position can be calculated using equation (20), which yields r v = 0.445 GeV −1 .
Within the range of 0 < r < r v , we can calculate α using equation (22), as illustrated in Fig. 6.It can be observed that α exhibits a decreasing trend as r v increases.Then, the separation distance of a heavy-quark pair and its corresponding energy can then be obtained using equations (11) and (16), respectively.These results are presented in Fig. 10 and Fig. 11.
When r v increases to r v, it reaches a critical value and transitions to the second configuration.In the second configuration, we can still determine the position of antiquark using equation (18).Subsequently, we use equation (27) to obtain α, as depicted in Fig. 7. α exhibits an increasing trend as r v increases, reaching its peak when r v = r v. Next, we can obtain the separation distance and energy of the second configuration using equations (11) and (25), and it is observed that as r v increases, L and E also exhibit an upward trend.
When r v reaches r v, it attains its maximum value.Beyond this point, the configuration switches to the third configuration.
In the third configuration, r v overlaps with the r q point.Using equations (32), we can establish a functional relationship between r v and α, as shown in Fig. 8. Clearly, α exhibits a linear decrease with increasing r v until it reaches 0. E and L can be obtained using equations (11) and (30).In this configuration, E varies linearly with L. When r v exceeds r q, the configuration shifts into the fourth configuration.
Then, in the configuration with large L, as r v increases below the dynamic wall, strings (1) and (2) exhibit turning points.α continues to decrease as r v increases, as shown in Fig. 9.However, unlike before, α becomes negative.As shown in Fig. 9, the maximum value of r v is 1.40 GeV −1 , which corresponds to the position of the dynamic wall at r w ≈ r 0 ≈ 1.45 GeV −1 .We can obtain E and L using equations (39) and (42).
At T = 0.08 GeV, the separation distance and energy plots for the tetraquark configuration are illustrated in Fig. 10 and Fig. 11, respectively.As can be observed from the figure, there is a smooth connection of L for each configuration.As r v increases, L also approaches infinity.This implies that r v cannot be infinite, and when the maximum value is exceeded, the configuration will decay into another state.We will discuss this further in the upcoming section.The energy increases as the separation distance increases.Notably, the energy is dominated by Coulomb potential at small separation distances and linear potential at large distances.

B. T=0.115GeV
At a temperature of T = 0.115 GeV, the configurations of QQqq are deconfined.For a deconfined QQqq, the melting of QQqq can happen at a certain separation distance.If the separation distance of heavy quark-anqitquark pair is small enough, the QQqq will not melt even at high temperature.When increasing the distance of heavy quark-anqitquark pair, the color screening becomes important and the QQqq will melt.In summary, the dynamic wall disappears and QQqq will melt at a sufficiently far distance in this stage.
Firstly, similar to T = 0.08 GeV, we focus on the first configuration.The position of the antiquark can be determined using equation (18), which yields r q = 1.1677GeV −1 when T = 0.115 GeV.Next, we can calculate the position of the antibaryon vertex using equation (20), which yields r v = 0.4626 GeV −1 , and then use equation (22) to determine α.L and E can still be obtained from equations (11) and (16), respectively.We can then proceed to investigate the second and third configurations using a similar approach as the T = 0.08 GeV calculation.This will enable us to determine the separation distance and energy.As the value of r v continues to increase, the value of L also exhibits a tendency to increase.
However, as L approaches a maximum value, it will tend towards infinity, indicating that the quark is now free and large L configurations will no longer be possible as shown in (11) and (25) are employed to calculate the separation distance and energy.However, at this temperature, as r v keeps increasing, the configuration will eventually collapse, leading to the quark becoming a free state.Therefore, there exists a maximum value of L, beyond which L will tend towards infinity.Consequently, large values of L are not possible at this temperature.

D. T=0.15GeV
At this temperature, using the force balance equation, the position of the quark or antibaryon vertex is determined to be r v = 0.4726 GeV −1 .However, we have discovered that at this temperature there is no solution for the position of the antiquark, r q.Therefore, we have computed the separation distance and energy for the first configuration, which are displayed in Fig. 14.As r v surpasses a certain value, the QQqq configuration collapses, resulting in the quark becoming free.

E. Short summary
Based on the calculation results of the QQqq potential at the above mentioned four temperatures, it shows that both r q and r v become larger and gradually approach the position of the dynamic wall as the temperature increases.Furthermore, the melting of the QQqq configuration happens at a small distance as the temperature increases.Comparing the two lines for T = 0 and T = 0.11 GeV in Fig. 15, we observe that as temperature increases, the same separation distance L corresponds to a lower energy value.Besides, at small distances, QQqq exhibits Coulombic behavior, whereas at large distances, the behavior is linear at finite temperature [60][61][62].16.Three disconneccted configuration of QQqq.

A. QQqq
During the confinement phase, quarks are confined within hadrons.However, as the distance between heavy quarks increases, the strings connecting them eventually break [63].
In this context, we shall consider three decay modes for QQqq configurations.
Fig. 16 displays the configuration diagram for the three possible decay modes.We will proceed to analyze the energy of each disconnected configuration.Qq consists of a fundamental string and an antiquark, with the total action given by S Qq = S NG + S q .QQq comprises two strings, a vertex, and a light quark.The total action can be written as NG + S vert + S q.On the other hand, Qqqq consists of four strings, a vertex, a antibaryon vertex, a light quark and two antiquarks.The total action is S Qq q q = 4 i=1 S (i) NG + S v + S v + 2S q + S q .Furthermore, q q q consists of a antibaryon vertex and three antiquarks, with the total action S q q q = S v + 3S q.We employ the same static gauge as before and obtain the total action for each configuration by specifying appropriate boundary conditions.The detailed calculation process of the energy of each configuration is given in 0.2 0.4 0.6 0.8 1.0 FIG. 17. (a)QQqq and his three disconnected configuration energy at temperature T = 0 GeV.
(b)QQqq and his three disconneccted configuration energy at temperature T = 0.08 GeV.The black line is the energy of E QQq q, the blue line is the energy of 2E Qq , the orange line is the energy of E QQq + E q q q, and the green line is the energy of E Qqq q + E Qq .
S 1 = S Qq + S Qq S 2 = S QQq + S q q q S 3 = S Qq q q + S Qq . ( The results of the calculation are shown in Fig. 17.Here, we define the string-breaking distance as the intersection point of the two energies.As shown in Fig. 17(b), we can see that when L QQq q = 0.1870 fm, QQqq will decay into Qq + Qq, and when L QQq q = 1.3036 fm, QQqq will decay into Qqqq + Qq.However, at zero temperature (as shown in Fig. 17(a)), the former decay takes place at L QQq q = 0.1874 fm, while the latter occurs at L QQq q = 1.3147 fm.The first scenario is linked to the process of vertex annihilation, whereas the second pertains to the occurrence of string breaking through the production of light quark pairs.The QQqq −→ Qq + Qq is the most possible decay mode.The presence of small temperature will increase the string-breaking distance a little bit.There is always a energy difference for QQqq and QQq + qqq as shown in Fig. 17.Thus, the decay mode QQqq −→ QQq + qqq also will not happen.
0.2 0.4 0.6 0.8 1.0 In this part, we will discuss the difference between QQqq and QQq at a temperature of 0.08 GeV.The QQq configuration has also garnered widespread attention.Similar to the QQqq configuration, the decay process QQq −→ Qqq + Qq will occur in the QQq configuration at high temperatures.As shown in Fig. 18, we calculated the two configurations separately as well as the energy after their decay.Then, QQqq decays to Qq at E = 2.3806 GeV(L = 0.1980 fm), and QQq decays to Qqq + Qq at E = 3.0167 GeV(L = 1.2646 fm).At lower energies and smaller separation distances, the QQqq configuration will decay, while QQq is more stable.This difference may result from two distinct mechanisms: string breaking by light quarks for QQq and string junction annihilation for QQqq.
Here, we also consider another relation in [29,37], which is This relationship is derived from heavy quark-diquark symmetry, as illustrated in Fig. 19.Furthermore, similar to the case at zero temperature, we find that relationship occurs at a very small separation distance when T = 0.08 GeV.After L = 0.2662 fm, the potential energy of QQqq and QQq + Qqq − q Q continues to increase.So there's a slight difference in energy, specifically, the energy of QQq + Qqq − q Q will be slightly higher than that of QQqq.In addition, we find that the energy difference between before and after equation will occur after L = 0.2396 fm when it is at zero temperature, and after L = 0.2662 fm when it is at T = 0.08 GeV.The increase in temperature will increase the critical distance for the appearance of energy difference.

VI. SUMMARY AND CONCLUSION
This paper focuses on the melting and breaking of strings in QQqq at finite temperatures by using the five-dimensional effective string model, and compares it with QQq under same conditions.During the confinement phase, a dynamic wall which strings cannot penetrate forms below the black-hole horizon.As the temperature increases, the system enters the deconfinement phase.In this phase, the dynamic wall disappears, causing the QQqq string to melt and the quarks to become free.We investigate the energy of QQqq at which string melting occurs at four different temperatures.Subsequently, the study of three decay modes of QQqq configuration at high temperature is continued, and the temperature at which the different decay modes occur is calculated.Finally, we compare the decay model of the QQqq configuration at T = 0.08 GeV with that of the QQq and conclude that the QQq configuration is more stable under the same conditions.Ultimately, we aim to provide additional insights for experiments through our study of the effective string model in the future work.
FIG. 3. Static configuration at an intermediate separation distance of a heavy-quark pair.The baryon vertex, antiquark and the antibaryon vertex are in the same position.The force acting on the point is indicated by the black arrow.

FIG. 4 .
FIG. 4. Static configuration at a large separation distance of a heavy-quark pair.There are two turning points at strings (1) and (2) and they are symmetric about the Y-axis.The force acting on the point is indicated by the black arrow.

FIG. 5 .FIG. 6 .
FIG.5.r 0 as a function of r v in the larger L configuration, where the black line is r 0 at T = 0, while the blue line is r 0 at T = 0.08 GeV.r w indicates the position of a dynamic wall.

4 FIG. 7 . 4 FIG. 8 .
FIG.7.α as a function of r v in the slightly larger L configuration.The unit for r v is GeV −1 .

FIG. 9 . 8 LFIG. 10 .
FIG.9.α as a function of r v in the large L configuration.The unit for r v is GeV −1 .

5 EFIG. 11 .
FIG. 11.Energy E as a function of separation distance L at T = 0.08 GeV.The black line represents the configuration with small L, blue line represents slightly larger L, orange line represents intermediate L, and red line represents larger L. The unit of E is in GeV and that of r v is in GeV −1 .

Fig. 12 [ 12 FIG. 13 . 15 FIG. 14 .
FIG.13.Energy E as a function of separation distance L at T = 0.12 GeV.The black line represents the configuration with small L and blue line represents slightly larger L. The unit of E is in GeV and that of r v is in GeV −1 .

5 EFIG. 15 .
FIG.15.The blue dash line is the energy at zero temperature, while the black line is the energy at T = 0.11 GeV.

FIG. 18 .
FIG. 18.(a) The energy of the QQqq and two Qq.(b) The energy of the QQq and Qqq+Qq.

FIG. 19 .
FIG. 19.(a) The energy of QQqq and QQq + Qqq − q Q at T = 0 GeV.(b) The energy of QQqq and QQq + Qqq − q Q at T = 0.08 GeV.