Proto-strange quark stars from density-dependent quark mass model

In this paper, we investigate the evolution of strange quark stars (SQS) from birth as proto-strange quark stars to maturity as stable SQSs at a zero temperature. We assume that self-bound free quarks form {entirely the compact} star and study its evolution {through a series of snapshots} using a density-dependent quark mass model. We consider $\beta$-equilibrated stellar matter at two major stages of the star's evolution: neutrino trapped regime and neutrino transparent regime during the deleptonization and cooling processes of the star. We fix the entropy density per baryon and the lepton fraction to investigate the nuclear equation of state (EoS), particle distribution, temperature profile inside the star, sound velocity, polytropic index, and the structure of the star. Our results show that stars with higher neutrino concentrations are slightly more massive than the neutrino-poor ones along the evolution lines of the SQS. We obtain EoSs in agreement with the conformal boundary set through sound velocity, and also the 2 M$_\odot$ mass constraint for NSs was satisfied at all stages of the star's evolution.


I. INTRODUCTION
Studying compact objects is an important fundamental research subject in contemporary nuclear physics, astrophysics, and cosmology.It plays a crucial part in understanding the early universe's nuclear structure, astrophysics, and state of matter.Compact objects enable us to explore the behavior of particles under exotic conditions such as higher baryon densities and temperatures in their core.The SQSs are compact objects whose existence was hypothesized for the first time by Ivanenko and Kurdgelaidze in 1965 [1,2].Witten later suggested that the fundamental state of hadronic matter could be a strange matter (SM) [3].Nowadays, it has become a common phenomenon to explore the inner core of NSs assuming that it is composed of deconfined quark matter taking into account the clear mass-radius window imposed by the data from the Neutron Star Interior Composition Explorer (NICER) [4][5][6] and the gravitational wave laser interferometer (Advanced LIGO, Virgo, and KAGRA) [7][8][9][10][11].Besides, the results found in [12] indicate that to satisfy the conformal bound for speed of sound, the massive NSs with M ≥ 2 M ⊙ are expected to have a deconfined quark core.At higher baryon densities (several times the nuclear saturation density n 0 = 0.152 fm −3 ) in the core of the compact objects, the particles are tightly packed overlapping each other such that they are capable of converting into free quarks.Therefore, two different forms of stars are expected: a hybrid NS (an NS with free quarks in its core) [13,14] or a pure SQS [15].Generally, quark stars (QS) are composed of high up (u) and down (d) quark asymmetry and can be formed entirely by deconfined SQM.A stable SQS can be formed from u, d and s (strange) quarks with leptons at β-equilibrium [3,[15][16][17][18].
Proto-neutron stars (PNS) or black holes are formed after the gravitational collapse in the core of a massive star after it has exhausted its fuel supply depending on the initial conditions, triggering a type II supernova explosion.The PNSs at the stage of collapse trap neutrinos and become neutrino-rich objects.The star then goes through several stages of evolution after the core bounce, including heat transfer, neutrino diffusion, deleptonization, and entropy variation until the stellar matter becomes neutrino-poor, catalyzes, and starts cooling to form cold NSs several years after its birth [19][20][21].A study of PNSs with exotic baryons in their core can be found in [22][23][24][25] while PNSs with free quarks in their core can be found in [26,27].There is still limited information regarding the transition from PNSs to proto-strange quark stars (PSS) during supernova explosion due to the difficulty involved in the deconfinement simulation process from hadronic matter to SQM [17,19,28].That notwithstanding, in this work we assume that the entire supernovae remnant (PNS) is composed of SQM.In this case, we do not expect any phase transitions since PSSs are assumed right from birth using thermodynamic conditions characteristic of PNSs evolution [29,30] to investigate their evolution.Since SQM has been hypothesized to be the true ground state of quantum chromodynamics (QCD) [3,16,31,32], hadronic matter can transition into QM, at this stage, there might be only QSs instead of NSs [33].It has been shown in Ref. [33][34][35][36][37][38][39][40] that NS merger events can also lead to the formation of PSSs, which may be bare or have a dynamically irrelevant crust of nucleonic matter, although there are arguments against the existence of a crust in SQS [41][42][43].
The observation of massive compact objects sets a strict constraint on the EoS of the stellar matter, discounting most of the conventional phenomenological quark models that are not capable of producing heavy quark stars with 2 M ⊙ threshold.Additionally, the observation of binary merger GW190814 event [44] by LIGO/Virgo Collaborations, involving a secondary compact object with mass determined to be around 2.50 − 2.67 M ⊙ at 90% confidence level has attracted extensive attention leading to further strict constraints.The observation of the first authoritative binary NS merger event with the emission of gravitational wave signal GW170817 by LIGO/Virgo Collaborations [7,10] reinforces these constraints.A recent data from NICER [4,5] has led to the measurement of a massive pulsar PSR J0740+6620 with mass 2.072 +0.067 −0.066 M ⊙ and radii 12.39 +1.30 −0.98 km at 68.30% confidence level [6] introducing a clear mass and radii window within which NSs can be explored.Besides, some compact stars are candidates for even more massive neutron stars, such as PSR J2215+5135 with a mass of 2.27 +0.17 −0.15 M ⊙ [45] and PSR J0952-0607 with a mass of 2.35 ± 0.17 M ⊙ [46].
In this work, we explore newly born proto-strange stars formed during the gravitational collapse of a massive star as a neutrino-rich object to the formation of a stable neutrino transparent SQS at a zero temperature (T = 0 MeV).We assume that PNSs are composed entirely of self-bound quarks in a deconfined state and study their evolution as PSSs.We employ a density-dependent quark mass model (DDQM) to calculate the EoS, particle abundances (Y i ), the structure of the star, temperature profile, polytropic index (γ), and speed of sound (c s ).We determine how the particle distribution inside the star during its evolution influences the softness or the stiffness of the EoS and the other NS characteristics during the early stages of the star's life.The DDQM involves two free model parameters D (which represents the linear confining properties of the QM) and C (which is related to the leading order perturbative behavior of the QM) that require fixing in the model framework.Different authors have adopted different approaches to determine suitable values of D and C that produce stable QSs with a maximum mass within the 2M ⊙ constraint [47][48][49][50].However, in this work, we use C = 0.7 (dimensionless) and √ D = 130.6 MeV motivated by the results in [15] obtained using Bayesian analysis.This choice of parameters produces QSs that satisfy the mass-radius constraints imposed by PSR J0740+6620 [6] and PSR J0030+0451 [51] pulsars.It also leads to a stable SQS at T = 0 MeV based on the Bodmer and Witten conjecture [3,16,21] which hypothesized that energy density (ε) per baryon density (n b ) of SQM at zero pressure (P = 0) should be smaller than the energy per nucleon (E/A) of 56 Fe, i.e., (ε/n b ) SQM ≤ 930 MeV.In the same vein, the two-flavor quark system should have ε/n b greater than the E/A of 56 Fe, (ε/n b ) 2QM > 930 MeV, otherwise protons and neutrons would decay into free u and d quarks.Also, we ensure that the fixed entropy SQSs determined within the model framework satisfy the 2 M ⊙ mass window required for NSs.Some other studies have also been done regarding PSSs using other quark models with different focuses on NS characteristics to describe the evolution of the SQSs in Ref. [52][53][54].
The paper is organized as follows: In Sec.II we present the general overview of the DDQM and elaborate on the zero temperature formalism in Subsec.II A and the finite temperature formalism in Subsec.II B. In Sec.III we discuss the properties of the SQM and the required conditions for calculating the EoS of the SQSs.We present our findings and discuss them in detail in Sec.IV and our concluding remarks in Sec.V.

II. THE DDQM MODEL
It is common knowledge that particle masses vary relative to the media in which they are measured.Such mediadependent masses are the so-called effective masses.In general, the effective masses vary, but the coupling constants also run in such media.Models in which the particle masses depend on chemical potential and/or temperature are referred to as the quasiparticle models [55,56].Models of these kinds have been extensively studied in the literature [47,48,52,[57][58][59].
The original intention of the DDQM model was to use the density-dependent quark mass model to study the nonperturbative characteristics of the QM [60,61].It was later extended to calculate the EoS [62][63][64][65], viscosity of SQM and dissipation of r-modes [66], diquark properties [67], and compact astrophysical objects [15,47,52] among others.The two main problems associated with the DDQM are the proper determination of the mass scaling and the thermodynamic consistency of density-dependent particle masses.These problems have widely been addressed in Refs.[47,49,65,68] and references therein.In this work, we consider the case in which strong interquark interactions are caused by density and temperature-dependent quark masses.We investigate the evolution of the PSSs that require both fixed entropy density per baryon and T = 0 MeV applications of the model, so we highlight each of these cases briefly below.

A. Zero temperature formalism
Considering a three-flavor quark ( u, d, and s) model at T = 0 MeV with a finite chemical potential, µ i .The Fermi momentum ν i in phase space is an important parameter, unlike the finite temperature treatment in the next subsection.The number density of the free quark system can be expressed as and the energy density where the summation index i runs over all the particles present (i.e., i = u, d, s, e(electron), and µ(muon); representing the quarks and the leptons present in the system at this stage) and g i is the degeneracy factor.Equations ( 1) and ( 2) are well-known expressions in literature.For free particles, ν i and µ i are related by the expressions It is general knowledge that quark masses depend on density and temperature.Ideally, the quark mass scaling is expected to be derived from the QCD theory, however, as it stands now, it is impossible due to the nonperturbative characteristics of the theory at low energies.In Ref. [68] a cubic root scaling law was reached at T = 0 MeV based on an in-medium chiral condensate and was extended to a finite temperature treatment in Ref. [48].Considering that the quarks interact with each other, we can introduce the interaction through density-dependent quark masses defined as where m i0 is the current quark mass and m I is the density dependent parameter representing quark interaction.We require that the relation for the number density from Eq. ( 1) and the energy density from Eq. ( 2) maintain their form after introducing Eq. ( 4), which is the original idea of Ref. [60].Such masses have been named equivalent mass to differentiate them from other mass concepts [69].Using the equivalent mass, both ε and n i are required to maintain their original form like the free particle system, however, in the quasiparticle model as the one under consideration, the n i remains unchanged but the ε changes.Such models are known in the literature as the quasiparticle models.A detailed derivation of the corresponding ε in the quasiparticle model using fundamental thermodynamic consistency can be found in Sec.IIB of Ref. [70] and in Sec.II of Ref. [65] they used general ensemble theory.We do not intend to repeat these derivations here, however, the main point to note is that the Fermi momentum is not directly connected to the chemical potential, on the contrary, the chemical potential is replaced by an effective chemical potential µ * i .Hence, the new form of ν i becomes ν i = µ * 2 i − m 2 i , while the real chemical potential µ i is related to the µ * i as then with µ I the interaction term.The thermodynamic potential Ω 0 is determined to be the pressure P is expressed as with and the energy density is given by We adopt the cubic root mass scaling law for the equivalent mass [68,71] expressed as where n b is the baryon density and D is a parameter associated with confinement, QCD string tension, chiral restoration density, and chiral condensate in a QCD vacuum.On the other hand, C corresponds to the strength of onegluon-exchange interaction, the first-order perturbative interaction, deconfinement phase transition, and it takes on positive values.

B. Finite temperature formalism
Since we are dealing with hot QM with masses dependent on temperature T , and baryon density n i of each particle i, it is convenient to work with state variables T , n i , and the volume V .We start with the thermodynamic free energy density F and derive the other quantities such as the energy density ε and pressure P through self-consistent thermodynamic treatment.We express the F in the same way as the free particle system with the free particle mass and chemical potential replaced with n b and T dependent particle mass m i → m i (T, n b ) and chemical potential where µ * i is the effective chemical potential of the ith particle and Ω 0 is the thermodynamic potential of the noninteracting particles which should be considered as an intermediate variable for now.The true thermodynamic potential Ω will be defined at the end of this section as a function of Ω 0 .Since the independent state variables defined above do not include the µ * i that appears in Ω 0 , we need to find a mechanism to link it to the independent variables.We create such a linkage through the particle number density n i i.e.
this expression has been used in previous thermodynamic treatments with µ * i replaced with the real chemical potential, µ i .Here, we use µ * i for thermodynamic consistency.Taking the derivative of Eq. ( 12) we have with and Substituting Eqs. ( 15) and ( 16) into Eq.( 14) and grouping the terms, we have Comparing the above expression with the thermodynamic relation we can identify where S is the entropy density of the particles, and The last term in Eq. ( 20) persists when the finite size effect of the system cannot be neglected whether the particle masses are fixed or not.However, it can be ignored in this study just like other previous studies of DDQM [47,48,65,71] considering an infinitely large system of QM like the one under consideration.In this case, the free energy density, Eq. ( 12), is independent of the volume of the system, thus the volume term in Eq. ( 20) does not appear.Now substituting Eqs. ( 12) and ( 21) into Eq.( 20) and considering an infinite QM system, the pressure becomes The energy density is determined by substituting Eqs. ( 12) and ( 19) into the Helmholtz free energy density relation, and the actual thermodynamic potential becomes Therefore, a given independent state variables T and n i , the µ * i can be determined by solving Eq. ( 13) and the other thermodynamic quantities such as S, µ i , P , and ε can also be calculated from Eqs. ( 19), ( 21), ( 22) and ( 23) respectively.
The net contribution of the Ω 0 for the particles in the system can be expressed in terms of particles and anti-particles present, with Ω + 0 and Ω − 0 representing the contributions of the particle and anti-particle respectively.Explicitly, where g i is the degeneracy factor and i = u, d, s, e, µ, and ν e (electron neutrino) represents the quarks and leptons present in the system.The net number density, comprising particle (n + i ) and antiparticle (n − i ) respectively.The derivatives that appear in S, µ * i , P, and, ε can also be determined as and where ϵ i = (p 2 + m 2 i ) is the single particle energy.The mass formula expressed in Eq. ( 11) was extended to the finite temperature case [48,49] to study the deconfinement phase transition, in the form where Λ = 280 MeV is the QCD scale parameter whose value depends on the renormalization scheme [72,73], C = 0.7 (dimensionless) and √ D = 130.6 MeV.We use Eqs.( 11) and ( 30) as the equivalent mass relation with current quark masses for u, d, and s represented by 2.16 MeV, 4.67 MeV and 93.4 MeV respectively as reported by the PDG [74].

III. PROPERTIES OF THE SQM
The SQM under investigation is composed of u, d, and s quarks with leptons.In the neutrino-trapped matter phase (corresponding to the first and second stages of the star's evolution in the snapshots, Fig. 1), we consider the three flavor quarks with electrons (e) and their corresponding neutrinos (ν e ), since muons (µ) only become relevant after the star has become neutrino-transparent in supernovae physics [26].In the neutrino transparent phase (the third and the last stages in the snapshots, Fig. 1) matter consists of the three flavor quarks, e, and µ.The quarks have degeneracy of g i = 6 (3 colors × 2 spins), and the e and µ with g i = 2 while neutrinos are characterized with g i = 1.The Fermi gas mixture then satisfies the β-equilibrium condition maintained by the weak interactions: d, s ↔ u + e + νe , s + u ↔ u + d, etc.In light of these interactions, the effective chemical potentials must satisfy the relations and neutrino transparent matter (32) with µ e = µ µ and µ νµ = µ νe , the subscripts of the chemical potentials µ represent the type of particle under consideration.These expressions also hold for the real chemical potentials of the quarks as well.
The charge neutrality condition is given by the expression the coefficients of the number densities (n) represent the electric charges of the particles and the subscripts represent the particle type, n µ = 0 for the neutrino-trapped matter.The baryon number is conserved through and the lepton fraction for the electron family (Y L,e ) for the neutrino-trapped matter is conserved by where Y L,e = Y e + Y νe and Y L,µ = Y µ + Y νµ = (n µ + n νµ )/n b ≈ 0 according to supernovae physics [26].Also, we calculate the particle fractions Y i of each constituent quark in the system through the expression where i runs over all the constituent quarks in the system.In Fig. 1 we show the snapshots of the particle abundances, Y i , in the evolution stages of the SQS from birth as PSS to maturity at T = 0 MeV.We used thermodynamic conditions that are characteristic of those of stellar collapse and core bounce [19].Generally, the initial condition before the eventual supernova explosion is denoted by t ∼ 0s in the literature.This stage is represented by the first panel (top-left) characterized by S/n b ∼ 1 and Y L,e ∼ 0.4 as the initial stage of the star's evolution, similar studies have been done in detail in [20,75] using different classes of EoSs.For instance, in [76] the authors choose S/n b = 0.8 and Y L,e = 0.35 to represent the initial condition of the star's evolution, and in [22,23,[25][26][27]77] the authors choose S/n b = 1 and Y L,e = 0.4 as their initial conditions, which is similar to the choice of this paper.Consequently, there is a closed set of thermodynamic conditions within which the evolution stages of PSSs can be chosen depending on the expected time interval under consideration.In the second panel (top-right), the neutrinos have already started to diffuse, leading to a deleptonization (Y L,e = 0.2) and an increase in entropy density per baryon (S/n b = 2).The top panels (ambient condition of core birth) show the stages in which the star is neutrino-rich at the early stages of the star's birth.The first snapshot from the left represents the Y i of the newly born star before the supernova explosion, t ∼ 0s, and the second panel is when the star is about 0.5 − 1s old after the core bounce.Two general scenarios occur immediately after the supernova explosion, from the first to the second stage (top panels) of Fig. 1: 1) If the supernova explosion is not strong enough to deleptonize the outer mantle of the PSS, the matter accretes again leading to the formation of a black hole.2) If the supernova explosion is strong enough to generate a significant pressure loss in the surrounding disk through deleptonization, the mantle collapses and accretion becomes less important and the proto-star continues its evolution.This work assumes the second scenario because if the first one occurs, the evolution of the PSS terminates at the second stage, and it will not continue to the next stage.
The bottom panels (ambient conditions after deleptonization) represent the Y i when the star is neutrino-poor.The first snapshot at this stage represents when the star is about 10 − 15s old, dominated by neutrino diffusion, deleptonization, and core heating.Here, the probability of the appearance of the strange quark increases relative to the previous stages.During the deleptonization and neutrino transparency period, the core continues to cool until the matter catalyzes when the entropy drops to S/n b = 0 at about 50s later, the core then continues cooling through thermal emissions reaching T = 0 MeV in about 100 years later [19].In between the second and the third stages (second-panel top and the first-panel bottom), when the star starts losing neutrinos rapidly, this can soften the EoS through deleptonization, this poses the possibility of black hole formation if the gravitational pressure is large enough.The evolution does not continue to the third stage if a black hole forms [19].After all the neutrinos have escaped from the core and the star is maximally heated (third stage in Fig. 1), the probability of s quark appearing in the stellar matter increases.This results in increasing the production of SQM, if the production of SQM is large enough, it will lead to the formation of a black hole [78,79].This phenomenon only occurs if there is accretion of matter at this stage.These processes have been described diagrammatically in [20].However, in this study, we assumed that the PSS evolved to its last stage of formation of cold-catalyzed SQS without forming a black hole along its evolution lines even though the possibility of black hole formation is not ruled out during the star's evolution.
Observing the snapshots across the top panels from left to right, we observe that the neutrino concentration reduces the relative abundance of u and d quarks and slightly suppresses the appearance of the s quarks to higher n b .Higher neutrino concentration means less u and d abundances and delayed appearance of the s particle.A similar observation is made in the evolution of PNSs where both nucleons and exotic baryons are present in the stellar matter.In these studies, higher neutrino concentrations were observed to delay the appearance of particles with strangeness [22,23,25,26].Also, the temperature rises in the core of the star as the entropy density per baryon increases and deleptonization of the stellar matter influences the appearance of the strange particle.The s quark appears earlier in a hotter matter than relatively cold ones, as expected.The strange quark shifts towards the lower baryon density region from the first, second, and third stages of the star's evolution which reflects the temperature rise in the core of the star along its evolution lines presented in Tab.I.In the last stage the appearance of the s quark delays to a relatively higher baryon density.In the first stage, s quark starts appearing at n b ∼ 0.66n 0 , in the second stage, s quark starts appearing at n b ∼ 0.53n 0 , in the third stage it starts appearing at n b ∼ 0.30n 0 and in the final stage it starts appearing at n b ∼ 0.59n 0 .Along the panels from top to bottom, left to right, the Y d and Y s increase while Y u decreases relatively, along the evolution lines of the star.It is worth mentioning that, in the lower panel, the muons do not appear in the Y i at T = 0 because they constitute less than 0.01% of the particle constituents.Also, the central baryon densities (n c ) are represented by vertical gray lines in Fig. 1.It shifts towards the higher baryon density regions along the evolution lines of the star in the neutrino-trapped matter (top panels).In the neutrino transparent matter (bottom panels) the n c increases as well along the evolution lines.Similar behavior is observed with the central energy densities (ε 0 ) shown in Tab.I. From Tab.I, we observe that the M bmax remains constant for the neutrino trapped regime and that of the neutrino transparent matter also remains constant whilst the other properties of the stellar matter changes through the evolution stages.In the same table, we fixed two different M b and determined the matter properties of these stars.The variation between the stellar radii for these stars along the evolution lines changes slightly but their gravitational masses change significantly.For a fixed M b , both the gravitational masses and the radii decrease throughout the evolution stages of the stars, and the central energy density and the temperature increase as the stars evolve, before reaching 'stability' at T = 0.It is worth mentioning here that the stiffness or otherwise of the EoS can be determined through the particle distribution in the stellar matter.It is known that the appearance of new degrees of freedom such as hyperons (strangeness-rich particles) in hadronic stars softens the EoS.A similar situation applies here, the higher the probability of the s quark appearing in the stellar matter the softer the corresponding EoS and vice versa.Consequently, the relative maximum masses of the PSSs can be predicted through the particle distribution since stiffer EoS means a higher maximum mass.
From previous works by [80][81][82] the isospin asymmetry (u − d quark asymmetry) was determined to be where n 3 = n d − n u is the isospin density and n b = (n d + n u )/3 is the baryon density for two flavor u − d QM.We can deduce that for pure neutron (proton) nucleon matter δ = 1(−1) based on their constituent quark structure.We observe along the evolution stages of the star that, the Y u decreases as Y d increases from the first stage of the star's evolution to maturity at the fourth stage.Hence, the stellar matter becomes more asymmetric as the star evolves.The u − d quark asymmetry decreases as n b increases towards the star's core.In the neutrino-trapped matter δ has a negative value at the higher density regions due to the higher abundance of u quark relative to the d quark beyond the n c .We can infer that the star is proton-rich at this point of its evolution.In Fig. 2 we show, in units of MeV fm −3 , the EoS for the SQM at different stages of evolution of the SQS, calculated from the equations presented in Subsecs.II A and II B and the required equilibrium conditions in Sec.III.We observe that the newly born star with higher neutrino concentration presents a stiffer EoS.As the star evolves through deleptonization its EoS starts to soften as the probability of s quark appearing in the stellar matter increases.The EoSs become stiffer with the increase of the lepton number, mainly at high densities.Hence, the EoS softens from the first to the fourth stage as the star evolves.This reflects the reduction of the maximum mass of the star as it evolves from the first to the last stage as shown on Tab.I. On the same table, we observe that fixed entropy stars are more massive than the colder ones, a similar observation was made in [52,54].From Fig. 1 we can deduce that lower lepton number concentration softens the EoS as we move along the panels.Additionally, we can deduce that higher isospin asymmetry corresponds to softer EoS.In Fig. 3, we determine the structure of the stars through their mass-radius (M-R) diagram.To find the M-R relation for the stars, we need to solve the Tolman-Oppenheimer-Volkoff (TOV) equations [83] assuming static spherically symmetric fluid.The equations are expressed as where r is the radial coordinate, M (r) is the gravitational mass, P (r) is the pressure, and ε(r) is the energy density, and we are using natural units (G = c = 1).Our analysis is motivated by the quasi-static approximation for proto-neutron star evolution in the spherical symmetric form, which results in the consideration that the star is in hydrostatic equilibrium.This eliminates the time derivatives of density, pressure, and the metric, leading to the use of average hydrodynamic effects that concentrate on the evolution of intensive thermodynamic quantities such as the Y L,l and the S/n b over the Kelvin-Helmholtz timescale [77,84], unlike the full hydrostatic simulation approach where matter is not necessarily in β-equilibrium.Since we have a static spherically symmetric spacetime with Schwarzschild gauge, the TOV equations are applicable (see e.g.[23-27, 85, 86]) because the system is completely in hydrostatic equilibrium throughout the star's evolution.In this case, the boundary conditions are P (R) = P surf , with P surf the surface pressure and R is the radius of the star.In the specific case of a cold-catalyzed NS, P surf = 0.The choice of P surf is arbitrary but has been found to have a significant effect on the characteristics of the mantle at the early stages of the star's evolution.On the other hand, it has been established that significantly low values of P surf have a negligible impact on the evolution of the internal structures of the star over long timescales [87,88].In Fig. 3, we show the M-R relation for the different values of the S/n b considered in this work.The precise position of the maximum mass M max and the corresponding R for each curve is shown by a dark point.We observe that the maximum mass of the PSSs decreases throughout its time evolution.The highest values of M max are achieved in the neutrino-trapped regime; in the first stage (S/n b = 1, Y L,e = 0.4) with M max = 2.20 M ⊙ and in the second stage (S/n b = 2, Y L,e = 0.2) with M max = 2.15 M ⊙ .After all the neutrinos have escaped from the core of the star, its maximum mass is further reduced; in the third stage (S/n b = 2, Y νe = 0) we have M max = 2.12 M ⊙ and the maximum mass continues to decrease until it reaches its smallest value, M max = 2.09 M ⊙ when the star cools and becomes a stable SQS at the fourth stage.From Fig. 3, we can deduce that over time, the SQSs shrink and become smaller, so their radii decrease along the stages of evolution, achieving the smallest values when the star catalyzes and contracts to become a cold SQS.Comparing this result to the one obtained in Fig. 1, we can deduce that lepton fraction and isospin asymmetry influence the mass of the star.A lower lepton fraction coupled with a high isospin asymmetry leads to a decrease in the maximum mass of the star and vice versa.As one always obtains, a stiffer EoS corresponds to a higher maximum mass along the evolution lines of the SQS.It is important to mention that, the structure of the PSSs and the cold SQS studied here have maximum mass within the 2 M ⊙ constraint required for NSs, such that, this model can accommodate the masses of PSR J0348+0432 (2.01 ± 0.04 M ⊙ ), PSR J0740+6620 (2.072 +0.067 −0.066 M ⊙ ), and PSR J1614-2230 (1.908 ± 0.016 M ⊙ ) [89,90], among others.Besides that, the radii determined here are well in the range of the one determined for the massive millisecond pulsar PSR J0740+6620 mentioned in the introduction.In Fig. 4, we show the temperature T variation in MeV as a function of n b /n 0 for different values of the S/n b .The core temperature, which corresponds to the central energy density of the stars where their maximum masses occur is shown with red stars on the plots with their corresponding n c , see details on Tab.I.In the first stage, when S/n b = 1, Y L,e = 0.4, the star is trapping neutrinos and heating up to cause an explosion, so it has a lower temperature profile.This stage has the lowest temperature profile compared to the other two stages studied that involve fixed entropies, with a core temperature of about 8.82 MeV at a central baryon density of n c = 3.95n 0 .In the second stage, a few seconds after the core bounce, the trapped neutrinos start to diffuse, heating the stellar matter thereby increasing the entropy density (S/n b = 2, Y L,e = 0.2) and causing deleptonization as the neutrinos escape from the core of the star.This process heats the star, increasing its core temperature to about 19.40 MeV at a central baryon density of n c = 4.34n 0 .From this stage, the neutrinos continue to escape the core of the star until they have completely escaped the stellar matter, at this stage (third stage) the stellar matter is maximally heated with S/n b = 2, Y νe = 0 which will start cooling through neutrino pair radiation and thermal emissions.This corresponds to a core temperature of 20.71 MeV at n c = 4.28n 0 .The second and the third stages of the star's evolution comprise the same S/n b = 2 with different neutrino concentrations but the temperature profiles of the stellar matter are significantly different.This is because the neutrino concentration suppresses the temperature of the stellar matter.Lastly, after becoming neutrino-transparent the star goes through catalyzation and continues to cool until it reaches T = 0 MeV several years later.We can also observe that the central energy densities of the stellar matter increase from the first to the second stage when neutrinos are trapped, and then slightly decrease in the third stage when the neutrinos have escaped from the core of the star and rise again when the star reaches T = 0 MeV.It is important to mention here that, the temperature profile for the proto-quark stars is significantly smaller [52,53,91] compared to its hadronic counterparts, PNSs [19,[22][23][24][25].As a result, observing the temperature profile of a compact object can enable us to determine its composition with reasonable certainty.
In Fig. 5, we show the pressure variations in two stars with M b = 1.85 M ⊙ (represented by dot-dash lines on the graph) and M b = 2.26 M ⊙ (represented by solid lines on the graph) as a function of the stellar radii.We observe that the star with a larger M b has a higher internal pressure and a higher M as shown on Tab.I. Additionally, in the neutrino transparent regime, when S/n b = 2 and Y νe = 0, and the stellar matter is maximally heated the pressure in the star drops due to expansion, and at T = 0 the star shrinks and its core pressure increases significantly.Figure 6 shows the temperature variations within the stars with M b = 1.85 M ⊙ (represented by dot-dash lines on the graph) and M b = 2.26 M ⊙ (represented by solid lines on the graph) as a function of the stellar radii.As expected, the temperature of the stellar matter increases steadily towards the tar's core.Also, higher S/n b corresponds to a higher temperature, and when all the neutrinos have escaped from the star's core the temperature reaches its maximum as shown on the graph.We can also notice that, despite having different temperatures in their interiors, the two stars tend to attain the same surface temperature at each stage.
In Fig. 7, we present the particle profile as a function of the stellar radii.We show the particle profiles of two stars with M b = 1.85 M ⊙ (represented by dot-dash lines on the graph) and M b = 2.26 M ⊙ (represented by solid lines on the graph) on the same graph.We observe that, the particle population of the d quarks and the e's decrease   slightly from the surface towards the core of stars, with the population of these particles for M b = 1.85 M ⊙ slightly higher than that for M b = 2.26 M ⊙ along the evolution lines.The u quark remains fairly constant in both stars for the same thermodynamic conditions along the evolution lines, however, their population reduces from the first to the second stage and increases from the third to the final stage when the star is neutrino transparent.The s quark and ν e population increase from the surface towards the core along the evolution lines and, where µ are present, their appearance is delayed and rises from the center towards the core.This signifies that the production of s, ν e , and µ are largely temperature dependent.Moreover, their population increases with the size of the star.The strongly interacting matter under quantum chromodynamics (QCD) theory shows different properties at both the perturbative and the nonperturbative regions of the theory.These properties can be distinguished through conformal symmetry arguments.The nonperturbative matter phase (confined quark matter) is not symmetric under the conformal symmetry transformation due to chiral symmetry breaking.In contrast, the perturbative matter phase (deconfined quark matter) is approximately symmetric under the conformal symmetry transformation.These properties can be quantitatively distinguished through the study of the sound velocity c s , that uses the EoS as input.It is known that c 2 s = 1/3 in exactly conformal matter while approaching the same value from below at high-density QM region, n b > 40 n 0 [92].In hadronic matter the c 2 s can rise up to about c 2 s ≳ 0.5, causality requires that c 2 s ≤ 1 and c 2 s > 0 for thermodynamic stability [12,93].In Fig. 8, we show the c 2 s in units of c 2 (constant speed of light) as a function of n b /n 0 for the stages of the star's evolution.We observe that c 2 s rises from the intermediate to the high baryon density regions as the star evolves and attains its lowest value at T = 0 MeV.Therefore, the c 2 s in hotter stars is higher than the ones in relatively colder stars.Also, the c 2 s monotonically increases with increasing n b /n 0 and approaches c 2 s = 1/3 from below at high-density regions, for all the stages considered, which is in agreement with the results obtained for conformal QCD matter [94][95][96].Thus, our results show that in the case The horizontal gray line represents the conformal limit where c 2 s = 1/3, below this limit matter is expected to behave in an asymptotically free manner, and above it, matter is expected to be in a confined state. of the evolution of SQSs within the framework of DDQM, it is possible to have compact objects with M max ≥ 2 M ⊙ formed by self-bound free quarks that conform with the conformal symmetry constraint for the speed of sound.The polytropic index as a function of n b /n 0 .The horizontal gray line, γ = 1, represents the conformal boundary.Below this limit, QM is considered to be in a deconfine state, and above it, QM is considered to be in a confined state.
The polytropic index γ is also used as a criterion for determining the stellar matter's inner composition.The γ is mathematically expressed as and uses the EoS as its input.For matter with exact conformal symmetry, matter without intrinsic scales, γ = 1, independent of any coupling strength.At this symmetry when matter is independent of any dimensionful parameter, energy density and pressure become proportional to each other which leads to γ = 1.This symmetry has also been proven to yield c 2 s = 1/3 as discussed above.On the other hand, chiral effective theory calculations and hadronic models predict γ ≈ 2.5 around and above nuclear saturation density [12,92,97].This is because at low to intermediate densities of QCD matter, the ground state does not possess the approximate chiral symmetry of the Lagrangian that governs it [98].The spontaneously broken symmetry results in the emergence of mass scales, like hadrons in the nuclear matter and scale-dependent interactions.In Fig. 9, we show γ as a function of n b /n 0 , the graph shows a steady decrease of γ with n b /n 0 , which slowly approaches the conformal limit at the high-density region when γ reaches a minimum value of 1.12.The γ ≈ 1.12 value is in good agreement with conformal matter classification considering that the limit set in Ref. [12] after examining a wide collection of EoSs is 1.75.The authors differentiated between hadronic and quarkionic matter phases with γ ≤ 1.75, which is the average between the γ calculated from perturbative QCD and chiral effective theory.

V. CONCLUSION
In this work, we investigate the evolution of SQSs from birth as PSSs to maturity when the star cools down to T = 0 MeV using a density-dependent quark mass model.The study assumes that the entire proto-star before the supernova explosion and after the core bounce is formed entirely by a deconfined SQM, so there are no phase transitions involved, PSSs are present from the beginning.We intend to extend the study to cover hybrid NSs in the future where PNS to PSS phase transitions become relevant [27].We considered four different stages of the star's evolution: Two stages when neutrinos are trapped in the core of the star and the other two stages when the neutrinos have escaped from the core.The neutrino-trapped stellar matter was investigated by fixing the entropy density per baryon (stage 1: S/n b = 1; stage 2: S/n b = 2) and the lepton fraction (stage 1: Y L,e = 0.4; stage 2: Y L,e = 0.2), we found that the EoS is stiffer when the entropy density per baryon of the stellar matter is low, and the lepton fraction is relatively high.In this phase, the stellar matter is less asymmetric, and the central energy and baryon densities increase from low entropy density per baryon to higher entropy density per baryon accompanied by an increase in core temperature from 8.82 to 19.40 MeV for stars with M max , see Tab.I.In the same table, we show the results for the core temperatures for two different stars with fixed M b , their values are significantly lower than the ones obtained for the stars with M max .In the neutrino-transparent phase, we consider the case where S/n b = 2, Y νe = 0 (stage 3) and T = 0 MeV (stage 4).In these phases, we found that the EoS for hot matter is stiffer than the one for cold SQM.The stellar matter becomes relatively more asymmetric compared to the previous stages and the central energy and baryon densities also increase from S/n b = 2 to T = 0 MeV.The star has a higher core temperature of T c = 20.71MeV for M max when all the neutrinos have escaped and S/n b = 2, for stars with a fixed M b the higher T c also occurs at this stage, these results have been elaborated on Tab.I.The stages of NS evolution from supernova explosion to maturity with the necessary thermodynamic conditions have been detailed in Refs.[19,75].
For the first time, we use DDQM to study the evolution of PSSs from birth, through deleptonization to the formation of cold-catalyzed SQS in a procedure similar to the one used in investigating PNSs in the literature [13,23,24].That notwithstanding, some progress has been made by several authors in describing the structure of the PSSs using different phenomenological quark models adopting different approaches.For instance, in [99], the authors explore strange quark matter and the structure of PSS using the same model as the one used in this work with a focus on conditions characteristic of the early formation of the star (e.g. they considered S/n b = 1 and 0.5 with Y L,e = 0.4).Their treatment did not take into account the deleptonization stages.In [86] the authors adopted the treatment of PSS similar to ours using the Polyakov chiral quark mean-field model considering conditions comparable to the birth of the star (S/n b = 1; Y L,e = 0.4) and neutrino transparent (S/n b = 2; Y ν,e = 0) stages.Their results qualitatively agree with ours in those instances.
Additionally, we investigated the particle abundances in the core of the star at different stages of the star's evolution and associated it with other properties of the SQM in Fig. 1.We observed that the isospin asymmetry of the stellar matter increases as the star evolves from the first to the fourth stage.Similar observations were made in Ref. [22] when they studied PNS with exotic baryons in its core.We also studied the EoS of the SQM and presented our findings in Fig. 2, we observed that the EoS softens as the star evolves from the first to the fourth stage which coincides with the direction of increasing isospin asymmetry.Some discussions on isospin asymmetry and its effect on the evolution of SQSs can be found in Refs.[52,54].We determined the structure of the star by calculating its M-R diagram in Fig. 3, we observed that both the mass and the radii decrease as deleptonization occurs and the star evolves from the first to the fourth stage.The PSSs down to the stable SQS satisfy the 2 M ⊙ constraint imposed by the observed pulsars PSR J0348+0432, PSR J0740+6620, and PSR J1614-2230 and the radii also satisfies the constraint imposed by the pulsar PSR J0740+6620 observed by NICER.We investigated the temperature profile in the star's core for the first three stages where entropy per baryon in the stellar matter was fixed and presented our outcome in Fig. 4. Generally, the temperature profile rises steadily with n b /n 0 .The first stage is when the star traps neutrinos, heats up and expands to cause an explosion, here, the stellar matter has the minimum temperature profile.The second stage is when the deleptonization starts with neutrinos escaping from the core, the stellar matter has the intermediate temperature profile.In the third stage when all the neutrinos have escaped from the core of the star, here, the stellar matter has its maximum temperature profile just before it starts cooling.This follows previous studies of PSSs [52,53,91] and PNSs [13,23,24].We also analyzed the pressure profiles of stars with a fixed M b in Fig. 5, where we observed that stars with higher baryonic masses have higher interior pressures.Also, the pressure inside these stars drops immediately after the neutrinos have escaped from their core, and then increases again at the last stage when the stars cool and shrink.The temperature profiles for the stars with fixed M b were also studied, and the results shown in Fig. 6, the temperatures for these stars are smaller at their surface and increase towards their center.Lastly, we also analyzed the particle fractions in the four stages of the star's evolution with the fixed M b 's in Fig. 7.We observed that the d quark and the e's decrease from the surface towards the star's core, the u quark remains fairly constant for both baryon masses and s quark, ν e and µ population increase from the surface towards the core in the direction of increasing temperature.Generally, some analysis for fixed baryon masses can be found in [100,101] and the references therein.
Furthermore, we investigated the behavior of the QM that are assumed to form the stars by studying the sound velocity and presented the result in Fig. 8.The sound velocity imposes a clear conformal boundary on the QM determined from the EoS, which helps to differentiate between the hadronic and quarkionic matter phases.The stars determined here produce c 2 s that lies below the c 2 s = 1/3 conformal line.This is evidence that the stars investigated are formed by self-bound free quark matter.Finally, we studied the polytropic index and presented our result in Fig. 9.The polytropic index also helps classify matter phases using conformal symmetry arguments similar to c s , here our results lie slightly above the conformal line, γ ≈ 1.12.However, this agrees with conformal matter classification following the new limit γ ≤ 1.75 set by the authors in [12].

FIG. 1 :
FIG. 1: In this diagram we show the snapshots of the particle distributions inside the star at various stages of its evolution.The upper panel shows the neutrino-trapped regime whilst the lower panel shows the neutrino-transparent regime of the star's evolution.The entropy densities S/n b are in units of Boltzmann constant k B .The vertical gray line in the Y i represents the central baryon densities, the position at which the maximum and radii occur.

0 FIG. 2 :
FIG.2:The nuclear EoS composed of free quarks and leptons in β-equilibrium at different stages of stellar evolution.

0 FIG. 3 :
FIG.3:The mass-radius relation for different stages in the evolution of SQS.The dark points on the curves represent the precise position of the star's maximum mass and the corresponding radii.

S/n b = 2 , 4 FIG. 4 :
FIG. 4:The temperature variation as a function of the baryon number density for neutrino trapped and neutrino transparent stellar matter is presented in the left panel.The red stars on the curves represent the position of the core temperature of the stellar matter.

0 FIG. 5 :FIG. 6 :
FIG.5: Pressure profile of a PSS star with a fixed baryon mass.We considered two stars with fixed baryon masses and calculated their corresponding pressure and radii.

FIG. 7 :
FIG. 7: Particle profile for two different stars.The dash-dot line represents M b = 1.85 M ⊙ and the solid line represents M b = 2.26 M ⊙ .

0 FIG. 8 :
FIG.8:The squared of sound speed as a function of n b /n 0 for both neutrino rich and neutrino poor stellar matter.The horizontal gray line represents the conformal limit where c 2 s = 1/3, below this limit matter is expected to behave in an asymptotically free manner, and above it, matter is expected to be in a confined state.

S/n b = 1 , 0 FIG. 9 :
FIG.9:The polytropic index as a function of n b /n 0 .The horizontal gray line, γ = 1, represents the conformal boundary.Below this limit, QM is considered to be in a deconfine state, and above it, QM is considered to be in a confined state.

TABLE I :
Stellar properties.Here, we show the maximum mass (M max ) and its corresponding radii (R), the central energy density (ε 0 ), the central baryon density n c , and the core temperature (T c ) for the evolution stages of the SQS considered.Additionally, we choose two fixed baryonic masses; M b = 1.85 M ⊙ and M b = 2.26 M ⊙ and calculate the other properties of the stellar matter for the four stages of the stellar evolution.