arXiv:2105.02116v1 [physics.gen-ph] 5 Jan 2021 rbesses[] unu rusadqatmal- quantum and groups inte- as Quantum quantum such of gebras, theory [1]. the systems by grable led model is group quantum which particle com- discovered symmetries generalized recently these by a described describe consider is which to to is method systems effective plex complicated An symmetry in generalized used concepts. more require be we cannot instead which problems, well corre- physics are their in transformations and known symmetry groups as Lie algebras many- The sponding in crucial studies. very quantum is body structures internal their and tem in ewe atce ipie ycnrligtede- between of interaction the the use describe controlling to The interac- by parameters formation simplified complex particles 16]. more [12– between makes [15, physics tions particles mechanics hole group statistical black quantum in [11], physics and state 14], solid [6– 10], physics in [9, experimental as and such theoretical compli- 8], much of and areas phenomena cated non-linear in widespread complicated of systems. investigation quantum in quantum many-body of crucial is interactions are area theory application field particle whose leaving invariant, groups symmetry. algebras quantum group existing quantum the the that still Therefore, obtained have is algebras it deformed deformed are algebras Although particle the [2–5]. parameters two with versions generalized omto aaee hc eut sa or as de- son results some which by q algebras parameter Lie formation and groups Lie usual underlying ∗ † [email protected] [email protected] nesadn h neato aueo atce sys- particles of nature interaction the Understanding h s fdfre emosadbsn a become has bosons and fermions deformed of use The q dfre emo siltragba,adtheir and algebras, oscillator fermion -deformed ters ihdniy aebe bandb sn eeaiaino generalization a comp using the by from obtained and been have mass density, same high of radius Schwarzschild sponding atce fsas h odadhtQSasi temperatures in Q- hot and cold Q-S The maximum the stars. of giving particles parameters deformation the consider ydfrainprmtr in parameters deformation by n aiso h trices hntesa esco on a down, cool gets Moreove radius. the Km when 75 increase and star mass the e solar of TOV 13 radius of to and up solutions reach numerical can the Q-Stars from investigated been have h -trsrcue,as alda ryhlswt radiu a with holes gray as called also structures, Q-Star The N SU .INTRODUCTION I. q uesmer SS)algebras (SUSY) supersymmetry 2 = and q ( n p ,aeseicdfrain fthe of deformations specific are ), stecnttet ftesa.Bcueteitrcinbet interaction the Because star. the of constituents the as opc -tr ihqpDfre Fermions q,p-Deformed with Q-Stars Compact aetp nvriy et fPyisEngineering, Physics of Dept. University, Hacettepe ekn nvriy niern Faculty, Engineering University, Beykent q ,p q, etp,Cnaa Ankara-Turkey Cankaya, Beytepe, yzg,Sryr Istanbul-Turkey Sariyer, Ayazaga, dfre bo- -deformed frin nta facmlctditrcinLagrangian interaction complicated a of instead -fermions he Mecit Ahmet Dtd a ,2021) 6, May (Dated: meDil Emre olwn esn;te r h otgnrloscillators general most the are having they reasons; following aees[1 2,wihcno edsrbdb ohone both by pa- described deformed deformation be parameter cannot the which 32], adjusting [31, by rameters phases bosonic fermionic between transitions and phase the describe to enables in ydfrainprmtr 1,2–8.Teueof use The parameters 21–28]. deformation [15, two parameters interac- deformation complicated by the tions describe to have introduced been counterparts also fermion one-parameter deformed bosons, two-paramater deformed and of Similar applications [17–20]. those structures to composite particles resulting the their between and interactions effective the vestigate ( ( parameter one-parameter example, For interaction complicated Hamiltonians. determining of instead particles, opc tr.Tepen ol lob ecie by of described constituents an be the having also as algebras could type preons particle preons special The deformed the on of is stars. one compact stars, to compact related is of parti- neu- which deformed of of algebras application core cle direct the Another in the stars. interactions of tron fundamental force are 38] pairing gen- of [37, nuclear formalism mass investigation the many-body (BCS) dynamical and The Bardeen–Cooper–Schrieffer the quarks for as to Q-Star. eration properties related compact hadronic applications on some other study some our are there fermions, proper- laser real the fit [36]. pa- to ties deformation observed two been the have of rameters For values critical non-linear studies. some some phenomenological instance, explain in to especially flexibility phenomena fermions more of provide deformation ad- to parameter third is two a using As of [35]. vantage model fermion undeformed standard hte h eomdfrin a eteconstituents ask the to be us can lead fermions studies deformed These the 40]. whether [39, invariance group Ozta¸s ∗ nadto oteaoeapiain ftwo-parameter of applications above the to addition In ¨ f SU of † q dfre emoswt w parame- two with fermions -deformed ftotmsgetrta h corre- the than greater times two of s T utosadi sotie htthe that obtained is it and quations q/p a rsuebtenteinteracting the between pressure tar damc esrsaei reached. is state denser much a nd c eto trfml ihavery a with family star neutron act 0, = ,p q, ( ,i sfudta h oa mass total the that found is it r, n eomdbsn r osdrdt in- to considered are bosons deformed ) unu ru naine[9 0 n it and 30] [29, invariance group quantum ) T enprilsaecontrolled are particles ween 0MVand MeV 40 = SU q 2 emo oe 3,3]and 34] [33, model fermion (2) ,p q, T sipratdet the to due important is 0MeV 60 = SL we , q 2 quantum (2) q n two- and ) 2 of a family of compact stars. where T is temperature, ǫ is the single state energy and In astrophysics studies, Q-stars are described as the µ is chemical potential of q,p-fermions. Then, the total stars having a very high density and mass, and a ra- number of particles and energy for q,p-fermion system dius about 2 times greater than the are obtained from [16] of the corresponding mass. They are categorized in the ∞ most dense spectrum of the neutron stars with exotic N = η(ǫ,q,p)C√ǫdǫ, (3) constituents such as SUSY Q-balls [41–43]. Because of Z0 their very high density they are also called as the gray and holes, meaning very close to form a . The stud- ∞ ies for describing the internal structures of Q-Stars have U = ǫη(ǫ,q,p)C√ǫdǫ, (4) been attempted with fermions [44], and obtained very Z0 complicated interaction Lagrangians. This motivated us where C = (V/2π2)(8π2m)3/2. After solving these equa- to investigate whether the internal structure of these Q- µ/T Stars can be described by q,p-deformed fermions which tion according to the z = e > 1 low-temperature case, can control these complicated interactions by the defor- the total number of particles and energy can be obtained mation parameters. as Therefore, we first use the thermodynamical properties 2 V 2 3/2 3/2 3 q,p T of q,p-deformed fermions called as Fermionic Fibonacci N = (8π m) µ 1 I (5) 3π2 − 8 µ oscillators. In order to solve Tolman—Oppenheimer- "   # –Volkoff (TOV) equations for our star formed by q,p- and fermions, we obtain the pressure, total number density 2 and energy density functions for high densities, as neu- V 2 3/2 5/2 15 q,p T U = 2 (8π m) µ 1 I (6) tron stars, much larger than the nuclear density. Then, 5π " − 8 µ # we solve the TOV equations numerically to obtain the   mass-radius relations of our q,p star model. The results where the q,p-deformed integral Iq,p is given by are consistent with Q-Star configurations. Also, we reach ∞ ∂η(q,p) the consistent and similar results with the fermion Q- Iq,p = (ǫ µ)2 dǫ (7) − ∂ǫ Stars results. We choose the parameters giving maximum Z0 pressure for the related number density. The resulting Moreover, the pressure of the system is obtained from mass-radius relations are much beyond the the equation of state (EOS) U = (3/2)P V for deformed structure and are consistent with Q-Stars because the fermions, such as q,p-fermions, as stated in above studies allows a transi- tion from fermionic to bosonic phase and construction of 2 2 15 T 2 P = ε = (8π2m)3/2µ5/2 1 Iq,p (8) BCS superconducting phases as in quark stars by con- 3 15π2 − 8 µ trolling the deformation parameters. "   # where ε = U/V is the energy density. Moreover, we know from the thermodynamical relation that P = Ω where II. THERMODYNAMICS OF Q-STARS WITH Ω is grand potential whose minimum value corresponds− q,p-DEFORMED FERMIONS to the equilibrium state. In order to investigate the core of Q-Stars in terms of q,p-fermions we need to find the Two-parameter q,p-deformed fermions are described zero temperature thermodynamical functions of the sys- by the following algebra between the creation and anni- tem. In the T = 0 limit, Eq.(5) turns out to be hilation operators ψ† and ψ, respectively [31] (8π2m)3/2 † † 3/2 ψ ψ + qpψ ψ =0, i = j nT (0) = 2 µ0 (9) i j j i 6 3π q ψ ψ + ψ†ψ =0, i

3 and

2 5 T 2 P = n (0)µ 1 Iq,p . (12) 5 T 0 − 4 µ "  0  # 168

In the core of compact stars, such as neutron stars, the 166 density is generally taken to be at least three times larger 164

−3 ) than the nuclear density n0 =0.122 fm [45]. Here, we 3 162 investigate the behavior of the Q-Stars for densities up 160

to 10 times the nuclear density, and also for a particular P(MeV/fm 158 value of 6 nuclear density. Moreover, in 156 studies, such as quark stars, the mass of the constituents 154 is taken to be around 2 3 MeV [45, 46]. Because the − 6 Q-Stars are more compact objects than the quark stars, 4 6 we consider the mass of constituent deformed fermions 2 4 0 2 as 7 MeV. For the total number density, Eq.(9) becomes −2 0 −2 −4 at T =0 −4 −6 −6 3 3/2 q p (8π27) 2 µ n (0) = 0 1.3015 10−7 fm−3 T 3π2 × −5 3/2 −3 −3 FIG. 1. Behavior of pressure with q,p deformation parameters =5.7118 10 µ0 fm =0.122n fm (13) × at temperature T = 40 MeV where the unit conversion GeV3 = 1.3015 102 fm−3 × is used since natural units are used h = kB = 1. More- over, the unit of chemical potential is taken in MeV. The pressure can be obtained at T = 0 from Eq.(12)

P (0) = 2.2847 10−5µ5/2 MeV/fm3. (14) × 0 The value of µ0 in terms of the number factor n of nuclear density can be obtained from Eq.(13) as 175

2/3 170 µ0 = 165.86 n , (15) )

3 165 then the pressure becomes 160

P(MeV/fm 155 P (0) = 8.0944 n5/3 MeV/fm3. (16) 150 Then, we obtain the pressure (12) in terms of the number 145 n factor of nuclear density at finite temperatures by using 6 4 Eq.(15) 6 2 4 0 2 − − −2 0 5/3 5 q,p 4/3 2 −2 P =8.0944n 1 4.5439 10 I n T (17) −4 −4 − × −6 −6 h i q p For a particular n = 6 nuclear density of the Q-Star, the behavior of the pressure with q,p deformation pa- rameters is illustrated in Figs. 1 and 2 for temperatures FIG. 2. Behavior of pressure with q,p deformation parameters T = 40 MeV and T = 60 MeV, respectively. The maxi- at temperature T = 60 MeV mum pressure which is equivalently the minimum grand potential value as the equilibrium state is obtained at q = 0.5325 and p = 0.3195 for deformation parame- ters. For T = 40 MeV, P = 165.23 MeV/fm3, and for T = 60 MeV,P = 171.78 MeV/fm3 are obtained. T = 60 MeV, respectively. q = 0.5325 and p = 0.3195 We can also find the chemical potential (10) at finite values of deformation parameters giving the maximum temperatures by using Eq.(15) pressure yields the chemical potential, at T = 40 MeV, µ = 544.09 MeV, and at T = 60 MeV, µ = 539.63 MeV µ = 165.86 n2/3 1+9.0878 10−6Iq,pn−4/3T 2 . (18) are obtained. × h i The behavior of the chemical potential is represented Moreover, we can find the total number density in fi- in Figs. 3 and 4 for temperatures T = 40 MeV and nite temperatures by using unit conversions in Eq.(5) as 4

T = 40 MeV and T = 60 MeV. After a certain chemi- cal potential value about µ = 400 MeV, the temperature has no effect on the dependency between number density and chemical potential. In addition, the higher values of chemical potential yield increasing number densities 552 independent from the temperature. 550 Similarly, we can obtain the energy density and the pressure in finite temperatures in Eq.(6) and Eq.(8) as 548

(MeV) 3 µ 546 2 2 5/2 (8π 7) µ −7 −3 ε = 2 1.3015 10 MeV.fm 544 5π × × 15 T 2 542 1 Iq,p 6 − 8 µ 4 "   # 6 2 4 2 0 2 −5 5/2 15 q,p T −3 −2 0 =3.4269 10 µ 1 I MeV.fm −2 −4 8 µ −4 × " − # −6   q −6 p (20)

and FIG. 3. Behavior of chemical potential with q,p deformation parameters at temperature T = 40 MeV 15 T 2 P =2.2846 10−5µ5/2 1 Iq,p MeV.fm−3 × − 8 µ "   # (21) The behavior of energy density and the pressure with re- spect to the chemical potential for temperatures T = 0,

560 T = 40 MeV, T = 60 MeV are given in Figs. 6 and 7, respectively. In these figures, higher temperatures 555 gives higher pressure and energy density values. More- over, the energy density and pressure values show a dra- 550 matic increase after a chemical potential value around

(MeV) 545 µ µ = 300 MeV. In Figs. 5-7, we infer that the increas- ing values of chemical potential µ gives increasing values 540 of number density nT , energy density ε and pressure P values of the Q-, respectively. 535

6 Moreover, the behavior of pressure in Eq.(17) and 4 6 chemical potential in Eq.(18) with respect to the num- 2 4 0 2 ber density for temperatures T = 0, T = 40 MeV, −2 0 −2 −4 T = 60 MeV are represented in Figs. 8 and 9. According −4 −6 −6 to these figures, when the number density n increases, q p T pressure and chemical potential increases for the system. On the other hand, we see that there occurs a high par- FIG. 4. Behavior of chemical potential with q,p deformation ticle production in high temperature regimes, although parameters at temperature T = 60 MeV the chemical potential is zero, or even less than zero due to the temperature, in Figs. 5 or 9. follows

3 III. CONFIGURATIONS OF HOT Q-STARS 2 2 3/2 (8π 7) µ −7 −3 nT = 2 1.3015 10 fm 3π × × We investigate the possible configurations of formation 2 3 T and evolution of cold and hot Q-Stars in temperatures of 1 Iq,p " − 8 µ # T = 0, T = 40 MeV and T = 60 MeV. In the literature,   Q-Stars are expected to have a maximum radius about 2 3 T two times greater than the Schwarzschild radius meaning =5.7118 10−5µ3/2 1 Iq,p fm−3 (19) × " − 8 µ # it is very close to the mass-radius relations of the black   holes. However, they are not yet dense as much as to The dependence of total number density on chemical po- form a black hole, therefore the gray hole term is used tential is illustrated in Fig.5 for temperatures T = 0, to express Q-Stars which enables light to escape from 5

450

10 400

350

8 300 ) 3 250 )

0 6 (n T

n 200 P(MeV/fm

4 150

100

2 p=0.3195 p=0.3195 q=0.5325 50 q=0.5325 T=0 T=0 T=40 MeV 0 T=40 MeV 0 T=60 MeV T=60 MeV

0 100 200 300 400 500 600 700 800 0 100 200 300 400 500 600 700 800 µ(MeV) µ(MeV)

FIG. 5. Dependence of total number density on chemical FIG. 7. Dependence of pressure on chemical potential for potential for temperatures T = 0, T = 40 MeV, T = 60 MeV temperatures T = 0, T = 40 MeV, T = 60 MeV

700 400

600 350

500 300

250 400 ) ) 3 3

200 300 (MeV/fm ε P(MeV/fm 150

200 100

p=0.3195 100 p=0.3195 q=0.5325 50 q=0.5325

T=0 T=0 0 T=40 MeV 0 T=40 MeV T=60 MeV T=60 MeV

0 100 200 300 400 500 600 700 800 0 1 2 3 4 5 6 7 8 9 10 µ(MeV) n (n ) T 0

FIG. 6. Dependence of energy density on chemical potential FIG. 8. Behavior of pressure with number density for different for temperatures T = 0, T = 40 MeV, T = 60 MeV temperatures T = 0, T = 40 MeV, T = 60 MeV its gravitational field. The constituent of Q-Stars are mass m(r) in the spherical radius of r is given by expected to be SUSY particles, or Q-Balls with a large r particle density. Here we investigate the Q-Stars with m(r)= ε(r′)4πr′2dr′ (23) q,p-fermions constituent. Z0 The Tolman–Oppenheimer–Volkoff (TOV) equations In order to solve these equation, it is necessary to de- for spherically symmetric configurations stars are applied fine the boundary conditions at the center of the star. to Q-Stars with q,p-deformed fermions [47], The mass is considered as zero, the pressure and number dP (r) [ε(r)+ P (r)] m(r)+4πr3P (r) density nT are maximum at the center. When the star = (22) dr r [r 2m(r)] reaches the maximum possible radius its mass becomes −  maximum, however the pressure and number density where ε(r) and P (r) are the energy density and pressure reaches almost a zero value. According to the thermo- values at a distance r from the center of the star. The dynamical relations in previous section, we first choose 6

12 800

10 600

8 400 )

6

(MeV) 200 µ m(M

4 0

p=0.3195 q=0.5325 2 −200 T=0 T=0 T=40 MeV T=40 MeV T=60 MeV T=60 MeV 0 0 1 2 3 4 5 6 7 8 9 10 0 10 20 30 40 50 60 70 n (n ) r(Km) T 0

FIG. 9. Behavior of chemical potential with number density FIG. 10. Solution of TOV equation for accumulated mass for different temperatures T = 0, T = 40 MeV, T = 60 MeV inside the Q-Star with central density nT (r = 0) = 6n0 for temperatures T = 0, T = 40 MeV, T = 60 MeV an n = 6 nuclear density initial value to investigate the dependence of the mass and the central number density 6 T=0 as well as the pressure on the radius for different tem- T=40 MeV peratures T = 0, T = 40 MeV, T = 60 MeV defining a T=60 MeV configuration. We used the cut off value for the number 5 density nT =0.5n0 because the model below this cut off density gives unphysical low-density quasi-free deformed 4 fermions. We illustrate the results in Figs. 10-12.

In Fig. 10, it is clear that higher temperatures lead to ) 0

(n 3 the less massive and smaller size Q-Stars. We can easily T n see from Fig. 11 that temperature has no role on the me- chanical equilibrium of Q-Star core around the center of 2 the configuration because of its high density. The effect of temperature is significant on the outer peripheral re- gions and changes the total mass of the star for a certain 1 central number density, as in Figs. 10 and 13. In addition, the temperature leads to the increase of 0 0 10 20 30 40 50 60 70 pressure but this gives a quick fall off the pressure of Q- r(Km) Star (see Fig. 12) and yields a smaller M R relation. When the temperature increases in Fig. 11,− the number density nT decreases, and then a less denser and smaller FIG. 11. Solution of TOV equation for matter distribution size Q-Star exists. inside the Q-Star with central density nT (r = 0) = 6n0 for By using the findings of above configurations, we can temperatures T = 0, T = 40 MeV, T = 60 MeV obtain the total radius of the Q-Star core R for different number density values ranging to 10 central nuclear den- sity. Moreover, the total mass M(R) can also be obtained density value. This fact is also clear from the Fig. 14 that for these varying number density values up to 10n0. The the higher temperatures gives smaller Q-Star configura- stable configurations of Q-Stars with maximum mass at tions in total mass and total radius. This implies that the the maximum radius for different temperatures T = 0, temperature rises the pressure and lowers the density of T = 40 MeV, T = 60 MeV are given in Figs. 13 and the star, which yields a smaller mass and a smaller radius 14 such that the mass as a function of central number Q-Star configuration. Moreover, there exist some Q-Star density and of the rotal radius, respectively. configurations having two different radii values for a sin- Fig. 13 shows increasing temperatures gives Q-Star gle mass value between R = [60 75] Km values for T =0 configurations with smaller masses. When the star gets case, between R = [40 50] Km− values for T = 40 MeV cooler we obtain more massive Q-Stars for a same nuclear case and R = [30 40]− Km values for T = 60 MeV case. − 7

180 14 T=0 T=0 T=40 MeV T=40 MeV 160 T=60 MeV T=60 MeV 12

140

10 120 )

3 8

100 ) sun

80 M(M

P(MeV/fm 6

60 4

40

2 20

0 0 0 10 20 30 40 50 60 70 0 10 20 30 40 50 60 70 80 r(Km) R(Km)

FIG. 12. Solution of TOV equation for pressure inside the Q- FIG. 14. Stable configurations of Q-Stars for different tem- Star with central density nT (r = 0) = 6n0 for temperatures peratures T = 0, T = 40 MeV, T = 60 MeV with total mass T = 0, T = 40 MeV, T = 60 MeV as a function of maximum radius

tions of maximum pressure and equilibrium suitable for 14 the discussion of Q-Stars. The EOS and related ther- modynamical quantities have been derived for the q,p- 12 deformed fermion system as the constituents of the Q- Star structure. The deformed model parameters are cho- 10 sen such that the pressure becomes maximum, or equiv- alently the grand potential becomes minimum implying

) 8 the equilibrium state. It is obtained that the increase sun of temperature leads to a higher average pressure and M(M 6 chemical potential values in Figs. 1-4. The dependence of thermodynamical functions such as

4 total number density, energy density and the pressure on chemical potential are illustrated in Figs.5-7 for temper-

2 T=0 atures T = 0, T = 40 MeV and T = 60 MeV. We infer T=40 MeV from Fig. 5 that the particles can still be produced for T=60 MeV

0 zero, or even less than zero chemical potential values for 0 1 2 3 4 5 6 7 8 9 10 n (n ) hot Q-Star profiles. However, after a certain chemical T 0 potential value around µ = 400 MeV, the temperature has no effect on the dependency of particle number den- FIG. 13. Stable configurations of Q-Stars for different tem- sity. The behavior of energy density and the pressure is peratures T = 0, T = 40 MeV, T = 60 MeV with total mass such that the higher temperatures gives higher pressure as a function of central number density and energy density values with a dramatic increase after a chemical potential value of µ = 300 MeV. After analyzing the thermodynamical relations and Similar case is also valid for the Q-Stars with one same EOS of the Q-Star profiles with q,p-fermion constituents, radius and having different total masses at this radius for we investigate the evolution of star interior with TOV the same regions of three temperature values. equations. As a first sample profile, we consider a 6n0 nu- clear density central number density for cold and hot star configurations at T = 0 and T = 40 MeV, T = 60 MeV, IV. CONCLUSIONS respectively. We show the results in Figs. 10-12 according to which the higher temperatures lead to the less massive In this study, we have investigated the effect of inter- and smaller size Q-Stars. We also infer that the tempera- polating property of deformation parameters on the ther- ture has no role on the mechanical equilibrium of Q-Star modynamics of q,p-deformed fermions under the condi- core around the center of the configuration because of its 8 high density. On the other hand, the temperature has an R = [30 40] Km values for T = 60 MeV case. Similar effect on the outer regions and changes the total mass of case is also− valid for the Q-Stars with same radius values the star. Moreover, when the temperature the increases but having different total masses at this radius for three pressure of the star increases but this is led by a quick different temperature cases. fall off the pressure of Q-Star as in Fig. 12, and also by a smaller M R relation. In addition, the increase of tem- As a concluding remark, we could obtain the Q-Star perature leads− to the decrease in number density and a structures, also called as gray holes having a radius of less denser and smaller size Q-Star, as it is clear from the two times greater than the corresponding Schwarzschild Figs. 11 and 12 for the pressure and number density ver- radius by considering the constituents as a generaliza- sus radius, respectively. The pressure remains constant tions of q-deformed fermions by two deformation param- for radius values greater than 10 Km and the number eters q and p. Instead of defining the interactions be- density of the Q-Star decreases for higher temperatures, tween the fermions by complicated Lagrangian function, which leads to less denser star profiles for hotter star we choose to consider q,p-deformed fermions in order to configurations. control the interactions by deformation parameters such In order to investigate the stable configurations, we that the maximum Q-Star pressure, minimum grand po- also obtain the total mass and radius relations of the tential are obtained for the equilibrium configuration of Q-Star by changing the number density values up to 10 the star. The most important result of the study is the central nuclear density in the TOV equations. Figs. 13 existence of Q-Star configuration that can reach up to and 14 illustrate the dependence of total mass and the 13 solar mass and 75 Km radius for these q,p-fermion total radius on the central number density, respectively. constituents, which is consistent with the general Q-Star In Fig. 13, we see the effect of higher temperatures as giv- structure in the literature. Moreover, it is found that ing smaller mass values for Q-Star configuration. When when the star gets cool down from high temperatures the star gets cooler we obtain more massive Q-Stars for the total mass and radius of the star increases, and a a same nuclear density value. This fact is also obvi- more denser state is reached. ous from the Fig. 14 that the higher temperatures gives ACKNOWLEDGEMENTS smaller Q-Star configurations in total mass and total ra- dius. Moreover, there exist some Q-Star configurations having two different radii values for a single mass value The authors are grateful to Elif Dil, PhD at Hacettepe between R = [60 75] Km values for T = 0 case, be- University, Statistics Department, for writing the matlab tween R = [40 50]− Km values for T = 40 MeV case and codes for numerical calculations. −

[1] A. A. Altintas, Czechoslovak Journal of Physics 56, 1069 [16] A. Algin and A. Arikan, Journal of Statistical Mechanics (2006). 2017, 043105. [2] A. Algin and M. Arik, Modern Physics Letters A 20, 613 [17] P. N. S. A.M. Scarfone, Journal of Physics A: Mathemat- (2005). ical and General 41, 275211 (2008). [3] M. Jimbo, Letters in Mathematical Physics 11, 247– [18] P. N. S. A.M. Scarfone, Journal of Statistical Mechanics (1986). 2009, P02055. [4] L. T. N.Yu. Reshetikhin and L. D. Faddeev, Algebra i [19] A. Gavrilik and A. Rebesh, Modern Physics Letters B Analiz 1, 178– (1989). 26, 1150030 (2012). [5] S. Woronowicz, Communications in Mathematical [20] I. K. A.M. Gavrilik and Y. Mishchenko, Journal of Physics 111, 613 (1987). Physics A: Mathematical and General 44, 475303 (2011). [6] M. Arik and D. Coon, Journal of Mathematical Physics [21] A. Lavagno and P. N. Swamy, Phys. Rev. E 65, 036101 17, 524 (1976). (2002). [7] L. Biedenharn, Journal of Physics A: Mathematical and [22] A. Lavagno and P. N. Swamy, Physica A 305, 310 (2002). General 22, L873 (1989). [23] P. N. Swamy, International Journal of Modern Physics B [8] A. Macfarlane, Journal of Physics A: Mathematical and 20, 697 (2006). General 22, 4581 (1989). [24] P. N. Swamy, The European Physical Journal B 50, 291 [9] V. Spiridonov, Modern Physics Letters A 07, 1241 (2006). (1992). [25] P. N. Swamy, International Journal of Modern Physics B [10] W.-S. Chung, Physics Letters A 259, 437 (1999). 20, 2537 (2006). [11] F. B. A.A. Marinho and C. Chesman, Physica A 391, [26] A. Lavagno and P. N. Swamy, Foundations of Physics 40, 3424 (2012). 814 (2010). [12] A. Strominger, Phys. Rev. Lett. 71, 3397 (1993). [27] A. Algin and M. Senay, Phys. Rev. E 85, 041123 (2012). [13] E. Dil, Canadian Journal of Physics 93, 1274 (2015). [28] A. Gavrilik and Y. A. Mishchenko, Ukrainian Journal of [14] E. Dil, International Journal of Modern Physics A 32, Physics 58, 1171 (2013). 1750080 (2017). [29] M. A. et.al., Zeitschrift f¨ur Physik C Particles and Fields [15] A. A. A. Algin and E. Dil, Physica A 416, 499– (2014). 55, 89 (1992). 9

[30] A. Algin, Journal of Statistical Mechanics 2008, P10009. [40] R. Finkelstein, International Journal of Modern Physics [31] M. A. A. Algin and A. Arikan, The European Physical A 34, 1941009 (2019). Journal C 25, 487 (2002). [41] B. Lynn, Nuclear Physics B 321, 465 (1989). [32] M. A. A. Algin and A. Arikan, Phys. Rev. E 65, 026140 [42] B. L. S. Bahcall and S. Selipsky, Nuclear Physics B 331, (2002). 67 (1990). [33] M. Ubriaco, Modern Physics Letters A 11, 2325 (1996). [43] T. S. J.C. Miller and L. Nolan, Mon. Not. R. Astron. Soc. [34] M. Ubriaco, Phys. Rev. E 55, 291 (1997). 294, L25–L29 (1998). [35] K. Huang, Statistical Mechanics (Wiley, 1987). [44] B. L. S. Bahcall and S. Selipsky, Nuclear Physics B 325, [36] J. Katriel and A. Solomon, Phys. Rev. A 49, 5149 (1994). 606 (1989). [37] L. Tripodi and P. L. B. . . C.L. Lima, Physics Letters B [45] D. B. et.al., Nuclear Physics A 736, 203 (2004). 412, 7 (1997). [46] D. B. et.al., Phys. Rev. D 72, 065020 (2005). [38] K. S. et.al., Phys. Rev. Lett. 93, 152501 (2004). [47] J. Oppenheimer and G. Volkoff, Phys. Rev. 55, 377 [39] R. Finkelstein, International Journal of Modern Physics (1939). A 32, 1730001 (2017).