Quark clusters, QCD vacuum and the cosmological 7Li, Dark Matter and Dark Energy problems

We propose a non-exotic electromagnetic solution (within the standard model of particle physics) to the cosmological 7Li problem based upon a narrow 2 MeV photo-emission line from the decay of light Glueballs (LGBs). These LGBs form within color superconducting, tens of Fermi in size, quark clusters (SQCs) in the radiation-dominated post-BBN epoch. The mono-chromatic line from the LGB ->gamma+gamma decay reduces Big-Bang nucleosynthesis (BBN) 7Be by 2/3 without affecting other abundances or CMB physics, provided the combined mass of the SQCs is greater than the total baryonic mass in the Universe. Following the LGB emission, the in-SQC Quantum-ChromoDynamics (QCD) vacuum becomes unstable and"leaks"(via quantum tunnelling) into the external space-time (trivial) vacuum inducing a decoupling of SQCs from hadrons. In seeking a solution to the 7Li problem, we uncovered a solution which also addresses the dark energy (DE) and dark matter (DM) problem making these critical problems intertwined in our model. Being colorless, charge neutral, optically thin and transparent to hadrons, SQCs interact only gravitationally making them a viable CDM candidate. The quantum tunnelling of the in-SQC QCD vacuum to the trivial vacuum offers an explanation of DE in our model and allows for a cosmology which evolves into a LambdaCDM universe at low redshift with a possible resolution of the Hubble tension. Our model distinguishes itself by proposing that the QCD vacuum within SQCs possesses the ability to tunnel into the exterior trivial vacuum, resulting in the generation of DE. This implies the possibility that DM and hadrons might represent distinct phases of quark matter within QCD, characterized by different vacuum properties. We discuss SQC formation in heavy-ion collision experiments at moderate temperatures and the possibility of detection of MeV photons from LGB ->gamma+gamma.


Introduction
The primordial abundances of the light elements produced in the first few minutes of the universe predicted by standard hot Big-Bang cosmology [1,2,3] are in excellent agreement with the abundances inferred from data (e.g.[4]).Big Bang Nucleosynthesis (BBN) starts when the deuteron (D) bottleneck is overcome at1 k B T ∼ 100 keV and terminates at k B T ∼ 30 keV (redshift z ∼ 4 × 10 8 ) due to electrostatic repulsion between nuclei (e.g.[5]); k B is the Boltzmann constant.Significant amounts of D, 3 H and 4 He build up followed by the production of much less-abundant elements such as 7 Be.With a half-life of ∼ 53 days, 7 Be decays into 7 Li via bound electron capture, with emission of a neutrino (e.g.[6]).This cannot occur, however, until recombination at z ∼ 1100 when 7 Be becomes singly ionized.The measured 7 Li abundance is ∼ 1/3 of what is expected from this process [7] defining the cosmological 7 Li problem (see [8] for a review).
The standard theory of electromagnetic cascades onto a photon background predicts a quasi-universal shape for the resulting non-thermal photon spectrum (e.g.[9]).In the case of non-thermal BBN, cosmological constraints using this quasi-universal shape make purely electromagnetic solutions to the 7 Li problem impossible (e.g.[10]) unless the injected photon energy falls below the pair-production threshold in which case the spectral shape is very different (see [11] and references therein).The effective pair production threshold is E γ,pair ≃ (m e c 2 ) 2 /22k B T ∼ 12 MeV × (keV/T) below which the double photon pair-creation process receives a Boltzmann suppression ( [10]; see also [12]); m e is the electron mass.If injected when BBN is over, sub-threshold photons act to post-process the abundances computed in the standard scenario.Photons injected with 1.59 MeV < E γ < 2.22 MeV can destroy 7 Be (suppressing 7 Li production) without affecting other BBN abundances, distorting the CMB background or injecting excess entropy [11,13].
In this paper, we sketch a model for the injection of sub-threshold photons (E γ < E γ,pair ) in the radiation-dominated post-BBN epoch, characterized by a universe temperature in the keV range.The condition 2 MeV < E γ,pair implies T CMB <∼ 6 keV, which corresponds to the post-BBN era when the universe is a few hours old and has a redshift below 6 keV/T CMB,0 ∼ 2.5 × 10 7 where T CMB,0 ∼ 2.73 keV.This model focuses on the nearly instantaneous emission (with a timescale of ∼ 10 −17 s) of a monochromatic 2 MeV line by quark clusters (QCs).These QCs are formed during the early universe's Quantum-ChromoDynamics (QCD) cross-over phase transition at a redshift of z QCD ≃ T QCD /T CMB,0 ∼ 6.7 × 10 11 , when the universe was a few microseconds old and had a temperature of T QCD ∼ 150 MeV.
Our model is grounded in the fundamental principles of QCD and relies on the properties of quark matter in a regime where quarks are deconfined and form Cooper pairs.Prior to delving into the details of our model, we provide a didactic introduction to the QCD phase diagram, the quark-gluon plasma (QGP), and other pertinent phases of quark matter in Appendix A (titled "Pedagogical Framework").Additionally, we offer a general overview of our model and the resulting interconnectedness of the 7 Li problem, Dark Energy (DE), and Dark Matter (DM) enigmas.
In our study, QCs consist of collections of up and down quarks.To ensure charge neutrality, the number density of down quarks (with a charge of −1/3) is twice that of up quarks (with a charge of +2/3).These QCs, which are assemblies of the quark-gluon plasma (QGP), have a radius of approximately 100 Fermi, rendering them transparent to photons.
Being composed of quarks, QCs are classified as baryons and are thus encompassed within the standard model of particle physics.They can be envisioned as macroscopic nucleons (not nuclei) carrying baryon number entirely through the quarks.Here and throughout the paper, baryonic matter includes all quarks which is more general than just protons and neutrons.The density of QCs is on the order of ∼ 10 39 cm −3 , which is approximately ten times the nuclear saturation density.Consequently, an approximately 100 Fermi QC possesses a baryon number A qc ∼ 10 6 .QCs are distinct entities from the significantly larger (A >> A qc ) cosmic strange-quark nuggets ( [14]) and Axion quark nuggets ( [15]), requiring different formation mechanisms (see bullet point #7 in Section 5).
In the post-BBN era, QCs transition into a color superconducting (CSC) phase, where their constituent quarks form Cooper pairs (see Appendix A) and give rise to what we refer to as Superconducting Quark Clusters (SQCs).For the purposes of this paper, we select the two-flavor color superconducting (2SC) phase as the reference phase, allowing us to convey the essence of our idea.Within the 2SC phase, a significant fraction of gluons (3/8 of them) do not interact with the paired quarks.Instead, they form their own condensate in the form of light glueballs (LGBs).Unlike the confined quark matter where GeV-scale glueballs are anticipated, LGBs within the 2SC phase exhibit a mass (M LGB ) in the MeV range.
Crucially, in the 2SC phase, it can be demonstrated that LGBs are electromagnetically unstable and decay into photons via the process LGB → γ + γ (by coupling with virtual quark loops which carry the electric charge).With a reasonable choice of QCD parameters, the 2SC phase accommodates LGBs with a mass of approximately M LGB c 2 ∼ 4 MeV.
During the post-BBN epoch, where kT CMB << M LGB c 2 , the decay of these MeV-scale LGBs generates a narrow, mono-chromatic ∼ 2 MeV line.This line emission presents a promising avenue for potentially resolving the enduring 7 Li problem.
The total mass of SQCs required to address the 7 Li problem is comparable to the observed value of Dark Matter (DM) in the Universe.Furthermore, the distinctive properties of SQCs -being colorless, cold, charge-neutral, optically thin, and decoupled from hadrons -position them as a plausible candidate for Cold Dark Matter (CDM).
We propose that the decay of LGBs within SQCs may lead to the destabilization of the in-SQC Quantum QCD vacuum.This destabilization, in turn, causes the in-SQC vacuum to undergo quantum tunnelling into the external (trivial) vacuum of space-time, a phenomenon we refer to as "leakage."Notably, this leakage process involves the transition to a vacuum state characterized by vanishing QCD (quark and gluon) condensates.
The tunnelling of the in-SQC QCD vacuum, behaving as Dark Energy (DE), occurs in the exterior space-time within our model.At low redshifts, the cosmological behavior aligns with the ΛCDM model.Therefore, in our model, the loss of gluonic content from SQCs due to LGB decay not only provides a potential solution to the 7 Li problem but also triggers the leakage process, leading to a cosmological scenario encompassing both Cold Dark Matter (CDM) in the form of SQCs and Dark Energy (DE) arising from the tunnelling of the in-SQC QCD vacuum.
The paper is structured as follows: In Section 2, we provide the derivation of the equations concerning the destruction of 7 Be through interaction with a 2 MeV monochromatic line.We explore the possibility of generating such a line through a lukewarm 2-flavor superconducting (2SC-like) quark phase, where a fraction (η G ) of the gluonic energy is converted into 2 MeV photons.The leakage phenomenon of the in-SQC QCD vacuum, involving a fraction ( f V ) of the SQC mass, is examined in Section 3. Here, we present the resulting cosmology and discuss how it may offer a plausible resolution to the Hubble tension, a discrepancy in the measurement of the expansion rate of the universe.In Section 4, we delve into the discussion of SQCs as a candidate for Cold Dark Matter (CDM).We analyze their properties and explore their viability as a constituent of the cosmic dark matter content.The limitations and predictions of our model are outlined in Section 5, offering insights into the boundaries and potential implications of our proposed framework.Finally, we conclude the paper in Section 6, summarizing the key findings and discussing the broader implications of our model.

Post-BBN 7 Be destruction
In our model, the CSC phase produces gluon condensation (i.e.M LGB c 2 ∼ 4 MeV mass LGBs) which decay to E 0 = M LGB c 2 /2 ∼ 2 MeV mono-energetic photons.To ensure that the energy of the 2 MeV line remains below the pair-creation threshold, we impose the condition 2 MeV < 12 MeV/T keV , which translates to T G < 6 keV.In other words, the QCs enter the color superconducting (CSC) phase and transform into SQCs during the radiation-dominated post-BBN era.Henceforth, quantities labeled with the letter "G" pertain to values associated with the LGB-decay/photon-burst event.This event occurs when the age of the universe exceeds approximately 10 hours, corresponding to a redshift The formation of LGBs takes place on a timescale comparable to that of hadronic processes when QCs enter the CSC phase.Subsequently, LGBs decay into photons almost instantaneously, with a timescale of τ LGB ∼ 10 −17 s, significantly shorter than the Hubble expansion timescale (as discussed in Section 2.1).
Subscripts "sG" and "eG" refer to the start and end, respectively, of this LGB-decay phase, denoted by t sG and t eG .The duration of this phase is approximately equal to the LGB decay timescale, i.e., t eG − t sG ∼ τ LGB .Consequently, we have z sG = z eG + δz, where δz is much smaller than z G .Effectively, this implies that z(t sG ) = z(t G ) = z(t eG ), indicating that the redshifts at the start, during, and end of the LGB-decay phase are approximately equal.
The photon production is a delta function, δ(t − t G ), in our model at time t = t G .As shown in Appendix B, the destruction rates for 7 Be nuclei due to such a sudden release of mono-energetic where Y is the abundance.Here, n B (t G ) is the universe's co-moving baryon number density at t G while n γ (E 0 , t G ) is the co-moving number density of photons from the LGB decay.The 7 Be photo-dissociation cross-section is σ Be (E 0 ) and σ CS (E 0 ) the Compton scattering cross-section.A reduction of 7 Be by 2/3 imposes that the RHS in Eq. (1) be unity (i.e.∼ −1.1).
Let us define η DM,sG as the initial total amount in mass of SQCs (the DM in our model) compared to baryons; i.e. before LGB decay.The corresponding SQC co-moving number density is with A sqc the cluster's baryon number.The emitted photons co-moving number density is then with ) the total number of E 0 photons emitted per SQC.Here, η G is the fraction of the gluonic energy per equivalent proton rest-mass converted to the mono-chromatic line at energy E 0 ; m p is the proton mass.Eq. ( 1) becomes where, in order to get the last expression, we made use of the fact that the SQC total mass after conversion of gluonic energy to photons is η DM,eG = (1 − η G )η DM,sG or η DM,sG = η DM,eG /(1 − η G ).
The 7 Be photo-dissociation cross-section, σ Be (E 0 ), is given by Eq. (III.8) in Ishida et al. [16].The ratio of σ Be (E 0 ) to Compton scattering cross-section, σ Be (E 0 )/σ CS (E 0 ), varies widely as shown in Figure 1.We see that for 2 MeV < E 0 < 2.2 MeV we have 3 × 10 −4 < σ Be /σ CS < 8 × 10 −4 .If the SQC total mass after conversion of gluonic energy to photons is the observed CDM amount, η DM,eG ∼ 5, then 5m p c2 /E 0 ∼ 2.4 × 10 3 .In this case, the 7 Li problem is solved if 0.65 which requires that on average ∼ 50% (i.e.η G ∼ 0.5) of the SQCs gluonic energy is converted to ∼ 2 MeV photons.As explained later in Section 5, the value of η G can be lower (less gluonic energy shed by SQCs) while still solving the 7 Li puzzle.This is because in addition to losing gluonic energy via LGB decay, SQCs will lose more during the leakage of its QCD vacuum so that if η DM,0 ∼ 5 in today's universe 2 then η DM,eG must have been higher due to leakage (see first bullet point in Section 5).The σ Be /σ CS ratio versus photon energy.The 7 Be photo-dissociation cross-section σ Be is from equation III.8 in [16].The Compton scattering cross-section, is given in Appendix (IV) in [10].
It is crucial to emphasize the requirement for a narrow band emission, preferably a mono-energetic photon source, in order to maximize the number of photo-dissociated 7 Be nuclei (as described by Eq. ( 1)).Without such a narrow emission, there would not be sufficient gluonic energy available from the SQCs -which serve as the Dark Matter (DM) in our model -to achieve the desired 2/3 reduction.In this context, the decay channel of LGBs presents an appealing mechanism for generating such a narrow emission (see Appendix C).The feasibility of this mechanism can be tested in ongoing experiments (as indicated in bullet point #5 in Section 5).

The mono-chromatic ∼ 2 MeV line
The QCD phase diagram is complex (e.g.[17,18,19]) and in principle one cannot exclude the existence of a CSC phase where gluon condensation would yield electromagnetically unstable X "particles" (here LGBs), as required in our model.Figure 2 shows a hypothetical phase diagram with the dashed line depicting a possible trajectory leading a quark cluster from its state at birth (in an unpaired phase) to the CSC phase.We require an unpaired phase at chemical potential µ < µ csc which bridges the hadronic phase and the CSC phase.I.e. the CSC phase is accessed, at low-T, from low density unpaired phase with a first-order line separating the two phases.Once in the CSC phase, a cluster becomes a SQC and produces LGBs in the radiation-dominated post-BBN epoch at keV temperatures.
In this paper, we will use the neutral 2SC phase as a reference CSC phase.It has the interesting property of converting a percentage (η G = 3/8) of its gluonic energy to LGBs (i.e.gluonic condensation) at low temperature with a subsequent decay to a mono-chromatic line via the LGB → γ + γ channel (see Appendix C for details).In the 2SC phase, an LGB with mass M LGB c 2 ∼ 4 MeV would decay to two E 0 = (M LGB c 2 /2) ∼ 2 MeV photons on timescales of τ LGB ∼ 10 −17 s.We set the CSC quark chemical potential at µ csc = 500 MeV with a corresponding number density n csc = µ3 csc /π 2 ∼ 10 39 cm −3 .Thus the density inside 3 an SQC is n sqc = n csc which is about ten times nuclear saturation density when it enters the CSC phase at time t G .
Once the cluster crosses into the CSC phase, the first order transition proceeds on hadronic timescales.During the transition, and because of latent heat released, a SQC gets heated to a temperature k B T sqc ∼ ∆ 2 csc /µ csc where ∆ csc is the CSC superconducting gap; here, ∆ csc and µ csc are in units of MeV (see Eq. A7).
LGBs cannot form at temperatures exceeding the melting temperature which is of the order of the LGB's restmass energy; k B T LGB,m ∼ M LGB c 2 .Thus we must ensure that T sqc < T LGB,m in addition ) depicts the cooling path of the quark cluster (formed at "1" with T QCD ∼ 150 MeV) traversing the unpaired phase and entering the CSC phase, at keV temperature, from the low density phase ("2").A first-order line separates the unpaired phase from the CSC phase releasing heat (the vertical "2 to 3" arrow).A SQC gets heated to a temperature not exceeding the LGB melting temperature k B T sqc,m ∼ M LGB c 2 (see Section 2.1).The ∼ 4 MeV LGBs decay to the 2 MeV narrow photon line via LGB → γ + γ. "1 ′ " shows the initial state of a NS core making its way to "2" via cooling and compression (see Section 5 and Appendix C).
to E 0 = E LGB /2 < 2.2 MeV where the 2.2 MeV upper limit keeps the line below the Deuteron photo-ionization threshold (see Section 1).In general for MeV is achieved for a reasonable range in µ as shown in Figure 3.The width of the ∼ 2 MeV line is given in Eq. (A8) in Appendix C, and for typical values, it is expected to be < 0.1 MeV.

The SQC size and baryon number
To first order, the SQC's photon mean-free-path is λ csc,γ = 1/(n csc σ qγ ) ∼ 10 2 fm/n csc,39 with σ qγ ∼ 10 −28 cm 2 the Thomson cross-section for photon scattering off up and down quarks.By setting a typical SQC radius to be R sqc,thin ∼ λ csc,γ , the SQC is optically thin to all photons including the 2 MeV photons.Being transparent to electromagnetic radiation of all energies, the SQC makes a good DM candidate (that is still within the standard model of particle physics) with a corresponding baryon number A sqc ∼ 10 6 .

SQCs and Dark Energy (DE)
Both hadrons and SQCs exist as states within QCD, but they exhibit different symmetries associated with their respective vacuum states.Hadrons, as discussed in [20], are characterized by confinement in their vacuum state, while SQCs are influenced by the formation and decay of LGBs to photons, leading to distinct symmetries.In Appendix C.1, we discuss the confinement of the SU c (2) phase within SQCs, where the vacuum state is speculated to undergo a tunnelling process, referred to as "leakage," into the external trivial vacuum.This tunnelling phenomenon is triggered by instabilities resulting from various mechanisms, including electromagnetic gluonic decay, such as the formation and subsequent decay of light gluonic bound states into photons in the CSC phase.
Although we chose to focus on the 2SC phase due to its simplicity, it is important to note that any CSC quark phase capable of converting gluonic energy into approximately 2 MeV photons and initiating a transition to the trivial vacuum could be equally effective in our proposed framework.This innovative concept offers an intriguing explanation for DE and establishes a cosmological framework that has the potential to address the Hubble tension while still yielding a universe consistent with the ΛCDM model at low redshifts.The total amount of confined SU c (2) QCD vacuum energy stored in the SQCs is sqc with ρ vac.QCD the density of the QCD condensates and V tot.sqc the total volume occupied by the SQCs which is constant in time (and assumed to not change on LGB decay)."Leakage" means that the density inside the SQC decreases from at t = t G (at the end of the LGB phase denoted by "eG") to at redshift z = 0 with the assumption that the leakage timescale is a fraction of the age of the universe t 0 ∼ 13.5 Gyrs (see below).The parameter f V = ρ vac.QCD /ρ csc is a measure of the contribution of the QCD condensates to the SQC rest-mass energy; naturally At the start of the LGB phase (denoted by "sG"), we have ρ sqc,sG = ρ csc .
We can get a rough estimate of the time t tun.(the e-folding time) it would take the in-SQC QCD vacuum to leak into the exterior space-time trivial QCD vacuum.A rigorous calculation would follow for example Coleman [21] and is beyond the scope of this paper.We consider instead a simple tunnelling problem across a square barrier with height given by the expectation value of the QCD vacuum and width L ∼ R sqc .The corresponding tunnelling probability is then P tun.∼ e −R sqc /δ where δ is the penetration depth (e.g.[22]).The tunnelling timescale is t tun.≃ (R sqc /c) × (1/P tun. ) giving us where x = R sqc /δ; the tunnelling timescale is highly sensitive to x.With δ of the order of a Fermi (which is not unrealistic), solutions with tunnelling timescales on the order of a billion years (t tun.∼ Gyr) require R sqc to be approximately 88 fm.If R sqc exceeds 90 fm, the tunnelling timescale t tun.becomes longer than the age of the universe, t 0 .Furthermore, the size of the SQC cannot be too small for leakage to occur on astrophysical timescales, ensuring that t G << t tun. .Therefore, an optimal range for the SQC radius is found to be 85 < R sqc (fm) < 88, which is consistent with our model.It is worth noting that in our model, optically thin SQCs with R sqc < 10 2 /n csc,39 fm are required to address the cosmological 7 Li problem.Furthermore, the size of the SQC cannot be too small for leakage to occur on astrophysical timescales, ensuring that t G << t tun. .Therefore, an optimal range for the SQC radius is found to be 85 < R sqc (fm) < 88, which is consistent with our model.It is worth noting that in our model, optically thin SQCs with R sqc < 10 2 /n csc,39 fm are required to address the cosmological 7 Li problem.
The time evolution of the density inside the SQC we can then write as ) .
The equation above incorporates the key parameters in our cosmology namely, ρ csc , η G , f V and t tun.which are all fundamentally related to QCD.The DM density in our model is the SQC density ρ sqc (t) averaged over the volume of the universe.Its time evolution, with V univ.(t) being the Hubble volume at time t, is (10) with the resulting cosmology analyzed in Appendix D.
Our model offers a resolution of the Hubble tension (see [23] for a recent review and references therein) and can be understood as a consequence of a CDM universe converting into a ΛCDM universe at z tun. .Figure 4 shows that H 0 ∼ 73 km s −1 Mpc −1 can be obtained for a range in η G and f V values with a leakage characteristic redshift 2 < z tun.< 6 (i.e. 1 < t tun.(Gyr) < 3.3).It is both noteworthy and encouraging that the timescale t tun., which is tightly constrained by the size of the SQC at approximately 100 Fermi, coincides with the timescale required to address the Hubble tension.This convergence underscores the interconnection of these two phenomena.Our cosmology yields a universe which is younger than the ΛCDM universe with an age of ∼ 13 Gyrs with η G and f V of the order of tens of percents each.It remains to be shown whether our cosmology is in agreement with other cosmological data and measurements which the flat ΛCDM model explains extremely well (see e.g.metrics and tests suggested in [24]).Furthermore, we caution that the details of how the QCD vacuum mixes with the space-time vacuum and how it evolves while preserving flatness remains to be understood.Nevertheless, being a vacuum the "DE" component in our model should obey an equation-of-state with parameter w = −1.

SQCs as cold dark matter (CDM)
Evidence for CDM is abundant but its nature remains unknown despite many theoretical investigations and dedicated experiments which have yet to detect any associated particle (e.g.[25,26,27,28,29,30]).In our model, the 7 Li problem is solvable if the total mass in SQCs is of the order of the measured CDM value.In addition, SQCs are colorless, cold, optically thin, electrically neutral and would make an ideal CDM candidate if they decouple from the strong force (or interact minimally with baryons) following leakage.η G : .10,f V : 0.10, z tun,73 : 5.35, age: 13.23 Gyr η G : .10,f V : 0.20, z tun,73 : 3.32, age: 13.09 Gyr η G : .10,f V : 0.  We have put forward the argument that SQCs that have undergone tunnelling would decouple from hadrons, resulting in their interaction being solely gravitational in nature.Consequently, these tunnelled SQCs would elude detection in current Dark Matter (DM) experiments.In the pre-LGB era, QCs should possess a sufficiently extended lifespan, allowing them to persist into the post-BBN era at z G without experiencing any growth in size due to interactions with hadrons.The rate at which non-tunnelled SQCs interact with hadrons depends on the fraction of SQCs that have survived until the present time.For instance, for a tunnelling timescale of t tun.∼ 1 Gyr, this survival fraction is e −t 0 /t tun.∼ 10 −6 , and it drops significantly to e −t 0 /t tun.∼ 10 −12 for t tun.∼ 0.5 Gyr.Below, we employ the term "baryon" specifically to refer to non-QC baryons (e.g.protons), allowing us to examine the interactions of QCs and SQCs with hadrons during these two distinct epochs.

The pre-LGB era
In the pre-LGB (i.e.pre-SQC at z > z G ∼ 2.5 × 10 7 ) epoch, QCs should interact with baryons.The QC-baryon reaction rate is σ qc−B n B v qc−B with n B (T) = η B n γ (T) the baryon's number density and n γ (T) ∼ 6.2 × 10 32 cm −3 × T 3 MeV ; η B ≃ 6.1 × 10 −10 is the baryonto-photon ratio.The QC-baryon interaction cross-section is σ qc−B while the QC-baryon relative velocity MeV is dominated by that of baryons; c is the speed of light and m p the proton mass.For the QC not to grow in size in the pre-LGB era (i.e. up to z G ) we impose MeV we have d(t/1 s) = −2T −3  MeV dT MeV so that 7.2 × 10 32 σ qc−B × T QCD,MeV T BBN,keV I.e.σ qc−B /M qc < 10 −12 cm 2 /g with M qc = A qc m p .The above ensures that the QC will neither grow in size in the pre-LGB era nor interact with baryons in the subsequent periods where baryons are more diffuse.We can express the effective cross-section in terms of the geometric one as σ qc−B = πR 2  qc × (τ crossing /τ conv.).Here, τ crossing = R qc /v sqc−B is the baryon QC crossing time and τ conv.the conversion timescale of a baryon to the CSC phase (i.e. for a baryon to be absorbed by the QC).The ratio τ crossing /τ B−qc,conv.accounts for the fraction of time the baryon spends inside the QC.
The fact that the conversion timescale τ conv.,B−qcexceeds the timescales associated with strong interactions provides support for the notion that the QC state is distinct from the hadronic phase.Alternatively, it suggests that the conversion process may involve two steps, where the baryon first undergoes deconfinement, possibly requiring an injection of energy, before transitioning into the unpaired quark phase characteristic of QCs.Another interpretation is that quarks within a hadron, where chiral symmetry is broken, perceive the QC as a distinct and separate phase of quark matter even prior it becomes an SQC.This intriguing aspect highlights the complexity and richness of the underlying physics involved in the conversion process and the distinct properties of QC and SQC phases.

The post-LGB era
From this point forward, when we refer to SQCs, we are referring to the clusters that did not undergo tunnelling.The SQCs that have undergone tunnelling, on the other hand, constitute the Dark Matter (DM) component within our proposed model.To impose a constraint on the SQC-baryon interaction cross-section σ SQC−B , we will utilize the measured lower limit on the proton lifetime (τ p > 3.6 × 10 33 years; [31]).This lower limit on the proton lifetime provides valuable information about the strength of the SQC-baryon interaction and allows us to further refine our model.
The SQC reaction rate per baryon is σ sqc−B n sqc v sqc−B and v sqc−B the SQC-baryon relative velocity.The SQC co-moving number density (not to be confused with the SQC density n csc ∼ 10 39 cm −3 ) is where ρ DM ∼ 0.36 GeV cm −3 is the local DM density [32] which is that of the tunnelled SQCs in our model.Imposing σ sqc−B n sqc v sqc−B < 1/τ p and using v sqc−B ∼ 240 km s −1 we arrive at σ sqc−B < 10 −48 cm 2 × A sqc × e t 0 /t tun. .
I.e.σ sqc−B /M sqc < (10 −24 cm 2 /g) × e −t 0 /t tun.with M sqc = A sqc m p .Consistent with expectations, the aforementioned analysis reveals that the efficiency of tunnelling, indicated by the presence of fewer remaining SQCs today, directly correlates with the permissible range for the SQC-baryon interaction cross-section σ sqc−B .Specifically, a more efficient tunnelling process allows for a higher allowable value of σ sqc−B .This relationship underscores the significance of the tunnelling mechanism in influencing the interaction dynamics between SQCs and baryons within our model.In terms of the baryon-to-SQC conversion timescale τ B−sqc,conv.(to be differentiated from the conversion of a baryon to an QC above), and with σ sqc−B = πR 2 sqc × (τ crossing /τ B−sqc,conv. ) we get τ B−sqc,conv.> 12.5 s n csc,39 × e −t 0 /t tun. .In section 3, we have established that a cosmological model capable of resolving the Hubble tension and simultaneously yielding the observed age of the universe necessitates a tunnelling timescale of approximately one gigayear (Gyr), corresponding to an exponential suppression factor of e −t 0 /t tun.∼ 10 −6 .Consequently, this implies that σ sqc−B < 10 −42 cm 2 × A sqc , or equivalently, that the conversion timescale τ B−sqc,conv.ex- ceeds 10 −5 seconds.
Intuitively, assuming σ sqc−B < σ qc−B would imply e −t 0 /t tun.> 10 −12 or t tun.> 0.4 Gyr which is also consistent wth our cosmology.We caution though that making direct comparisons between the QC and SQC states may not be justified.We should add that if SQC vacuum leakage is not a simple tunnelling process then it can deviate from a pure exponential process (e.g.[33]) which could affect our findings here.This is left as an avenue for future research.Furthermore, we can employ neutron stars (NSs) to impose even more stringent constraints on the SQC-baryon cross-section by studying the interactions between SQCs and baryonic matter within a NS.If σ sqc−B n NS R NS > A SQC , the SQC slows down inside the neutron star and is captured.For typical neutron star parameters, such as n NS ∼ 10 38 cm −3 and R NS = 10 6 cm, the condition for no capture implies that σ sqc−B < 10 −44 cm 2 × A sqc .Referring to Eq. ( 13), this condition is always satisfied if e −t 0 /t tun.> 10 −4 , indicating that neutron stars will not be converted through SQC capture.For a NS formed at time t NS , the mass of SQCs it accumulates over the Hubble time can be estimated as A sqc m p e −t/t tun.v sqc−NS dt × A sqc m p ∼ 10 11 g × t tun.,Gyr × (e −t NS /t tun.− e −t 0 /t tun.).
Here, we assume a SQC-neutron star relative velocity of v sqc−NS ∼ 10 7 cm s −1 and a representative local DM density of ρ DM ∼ 0.1 GeV cm −3 .Thus, even if a certain amount of SQC is captured, it is not evident whether it would be sufficient to convert a neutron star.
Due to the likelihood that the leakage timescale exceeds the typical timescale for the formation and subsequent collapse of structures during the early universe, we anticipate that primordial halos will be predominantly composed of non-leaked SQCs that formed after the Light Glueball (LGB) era.Baryonic matter will gravitationally accrete onto these halos, but as demonstrated earlier, the diffuse baryons will not interact with either leaked or non-leaked SQCs.However, we anticipate a critical change or "phase transition" in the growth of structures after the time t = t tun.This is because: (i) a significant fraction of the SQCs' mass is converted into Dark Energy (DE), and (ii) the SQCs would have decoupled from hadrons.Understanding the precise impact of these effects on subsequent structure growth requires more advanced investigations beyond the scope of this paper.Only state-of-the-art Cold Dark Matter (CDM) simulations, which incorporate the leakage of SQC mass as detailed in Appendix D, can provide insights into this aspect of our model.
To summarize, SQCs present a viable candidate for Cold Dark Matter (CDM) and may play a role in the formation of both mini-halos and larger structures in the universe, as indicated in previous studies [34,35].Once halos form and reach a state of virialization, SQCs are expected to interact most strongly with stars, particularly with matter at high densities such as neutron stars (NSs).However, if the tunnelling timescale t tun. is a small fraction of the age of the universe, as suggested earlier, then SQCs would have decoupled from hadrons by the time compact stars begin to form (as discussed in Section 3).In other words, even the densest stars, such as neutron stars, are unlikely to undergo conversion through capture of SQCs.However, if the cores of neutron stars can access the Color Superconducting Quark matter (CSC) phase through phenomena like cooling and mass accretion, it could have intriguing implications for astrophysics (as shown in equation ( 18)).This possibility opens up new avenues for research in understanding the behavior and properties of neutron stars within the framework of our model.

Discussion
Here we briefly discuss some limitations and distinctive features of our model and leave these as future investigations.

1.
The evolution of the SQC and of the DM densities: We first recall, from Eq. ( 9), that the density inside an SQC evolves from ρ csc to (1 − η G )ρ csc after the LGB phase.leakage further decreases the SQC density to an asymptotic value of (1 (which we associate with z = 0).The minimum possible value of an SQC would be given by the bare quarks masses or, ∼ 90% of their original mass.However, we find η G and η V to be of the order of tens of percents each (see Section 3), and an SQC would lose a fraction of its mass to DE over the age of the universe.Meaning that the DM content in today's universe is a fraction of that in the pre-BBN era according to our model.Averaging the total SQC mass over the Hubble volume gives us the DM density in our model.The ratio of total amount of DM to that of the baryonic matter evolves from a maximum value of η DM,sG before LGB decay to η DM,eG = (1 − η G )η DM,sG after the LGB decay to photons just before leakage starts.At full leakage, which we associate with z = 0, the ratio is Assuming η DM,0 ∼ 5 as measured today and extrapolating back to the LGB era, it means that the amount of DM content is larger than what we used in Section 2 where we set η DM,eG = η DM,0 ∼ 5 (i.e. when taking f V = 0).2.
The 7 Li problem revisited: The solution to the cosmological 7 Li problem presented in Section 2 did not take into account further loss of gluonic energy due to in-SQC vacuum leakage following LGB decay as discussed above.From Eq. ( 14), η DM,sG = 5/(1 − η G − f V ) and when plugged in Eq. ( 4), we get Eq. ( 5) becomes I.e. a smaller percentage of the gluonic energy of the SQCs converted to LGBs could resolve the 7 Li problem.For f V = 0 we recover the η G values we arrived at in Section 2 while for f V > 2/5 we get η G < 3/8 which is less than the maximum η G ∼ 3/8 expected from the 2SC-phase.
The ratio between the DE density and the DM energy density at z = 0 is (from Eq. ( 7)), To have ρ DE ∼ 2ρ DM as measured in today's universe requires f V ∼ 2 3 (1 − η G ) in which case Eq. ( 14) gives 0.18 < η G < 0.34 and 0.44 < f V < 0.55.On the other hand, the resolution of the Hubble tension suggests that a favored parameter combination lies around (η G , f V ) ∼ (0.1, 0.1), as illustrated in Figure 4.This discrepancy could potentially be mitigated by relaxing some of the approximations made in our cosmological model or by adopting a more realistic approach to the tunnelling process, moving beyond the simple assumption of an exponential decay.

3.
The CSC phase and Neutron stars: Figure 2 shows a suggested pathway, starting at point "1 ′ ", a NS core could take to enter the CSC phase since conversion following SQC capture is suppressed (see Section 4).A NS born with (or which acquires through evolution) a core in the unpaired phase could transition to the CSC phase by a sequence of cooling (to keV temperature via the URCA process; e.g.[36]) and compression (to µ csc = 500 MeV via mass accretion).Take a NS with a core making up a fraction η c of the total mass.The energy released from conversion of gluonic condensation (e.g.LGBs) to photons is Comparing this to the NS binding energy 3 5 GM 2 NS /R NS , we conclude that NSs with compactness parameter may be completely obliterated in the process; the NS mass and radius are in units of solar mass, M ⊙ , and 10 6 cm, respectively.NSs with higher compactness parameter would loose mass leaving behind a pure CSC core.In this latter case, the conversion to a CSC star puts a constraint on ρ csc due to the black hole limit 2GM NS /c 2 < R csc with R csc the radius of the CSC star.With NS,⊙ which is consistent with the µ csc ∼ 500 MeV (i.e. a 2SC-like phase) adopted in our model.Thus, if some NSs follow a path as suggested in Figure 2, the resulting photon fireball may have interesting implications to explosive astrophysics.

4.
SQC-hadron decoupling: Following leakage and loss of more gluonic energy, it is not unreasonable to assume that the quarks within the SQC become undressed and should in principle decouple (or at least experience some level of decoupling) from the strong interaction.They would still rely on gluons to remain bound while exhibiting minimal interaction with hadrons.We speculate that DM and hadrons may represent separate phases of quark matter within the framework of QCD, characterized by distinct vacuum properties.If our model is correct, it allows a unique connection between cosmology and the properties of the QCD vacuum in the CSC quark phase.An estimate of the parameter f V from cosmological observations may be an indication of the contribution of the QCD condensates in CSC to the SQC mass which may have implications, albeit an indirect one, to the mass of hadrons.The exact details of the SQC-Hadron decoupling remain to be worked out. 5.

SQCs and
LGBs in today's detectors: SQCs (the DM in our model; see Section 4) would interact only gravitationally and would thus evade detection in current DM experiments.Instead we propose that our model can be tested by experiments which can access the 2SC phase at temperatures below the LGB melting temperature; i.e. at < M LGB c 2 ∼ 4 MeV (see Appendix C).The Compressed Baryonic Matter experiment (CBM) at FAIR explores the QCD phase diagram in the region of high baryon densities (representative of neutron star densities) and moderate temperatures using highenergy nucleus-nucleus collisions (e.g.[37,38]).We note that the 2SC phase carries LGBs (i.e. can solve the cosmological 7 Li problem) at baryonic density as low as a few times that of nuclear matter (see Figure 3) and that LGBs form on strong interaction timescales and decay to MeV photons on timescales of ∼ 10 −17 s.Thus, experiments such as the CBM at FAIR could potentially be used to detect the photons from the LGB decay.
The width of the line can be found in Eq. (A8) in Appendix C, and for typical values, it is expected to be < 0.1 MeV.If this width were to be measured, it would provide a direct estimation of the 2SC dielectric constant, which could then be used to infer the ∆/µ ratio and consequently determine ∆ (as per Eq.A8).However, it remains unclear whether and how the CBM experiment could explore the temperature regime relevant to our study.More significantly, the 2 MeV line may be overshadowed by the anticipated high low-energy electromagnetic background in high-energy ion-ion collisions.Even if the background were to diminish as the QGP transitions into the 2SC phase and subsequently cools down to the LGB phase (T < 4 MeV), detecting the 2 MeV line would likely prove challenging.Calculating the strength of the LGB 2 MeV line is beyond the scope of this paper.

6.
The stability of the CSC phase: Our findings seem to hint at the standard neutral 2SC phase (adopted in our theoretical framework) as the unspecified CSC phase.However, the 2SC phase may be unstable at small superconducting gap (∆) values due to the mismatch in the up and down quarks chemical potential [39].It suggests that either the 2SC phase is stable in the regime of chemical potential (µ) and ∆ values we used or that another stable 2SC-like phase exits in nature and remains to be identified.A 2SCus phase, which has u-s pairs instead of u-d pairs, is also a candidate.This phase, however, faces a similar challenge in that the u and s Fermi momenta are split apart by the strange-quark mass (combined with electrical neutrality; e.g.[40]).The Color-Flavor-Locked phase may be one candidate if the strange-quark effective mass is small or alternatively a crystalline quark phase if strange-quark mass is heavier (e.g.[41]).However, these phases do not possess properties that allow for the conversion of gluonic energy to photons.Despite its limitations, the u-d 2SC phase remains as one of the candidate phases for dense quark matter (adopted in our model), pending the determination of its exact phase structure.
Another possibility are collective excitations (i.e.phonons) of SQCs which would decay directly to photons via phonon-to-photon conversion channels.Mono-chromatic photons of energy 1.59 MeV < E γ < 22.2 MeV would result if the SQC radius were 544.6 fm < R sqc ∼ hc/E γ ∼ 753.6 fm.Physically, the decaying particle is not an LGB but is instead a long-wavelength mode which is resonant with the quark cluster containing it.In this scenario, higher resonant modes should also be excited and these could affect the very sensitive D abundance which is not desired.Superconducting strings [42] is an interesting avenue to explore in this context.If such domain walls could form during the cosmic QCD phase, one could imagine a scenario where they would evolve into 2SC-dominated strings in the post-BBN era.It remains to be shown that LGBs as described in our model could form in this case.7.
Matter-Antimatter annihilation and SQC size: We hypothesize that each quark cluster is born with an anti-matter deficit of η B ≃ 6.1 × 10 −10 meaning that there is one extra baryon per η −1 B quark-antiquark pairs; η B is the baryon-to-photon ratio.After annihilation (on timescales of 1/n csc σ annih.∼ 10 −13 s with σ annih.∼ mbarn), a cluster has only baryons left in it; this assumes that annihilation does not destroy the cluster and instead it reduces it into a pure baryon cluster of radius R sqc,f ∼ R sqc,thin where "f" stands for final.Here, R sqc,thin ∼ 10 2 fm is the typical size of a SQC set by the photon mean-free-path in the 2SC phase.In other words, we claim that the maximum size of the "shrapnel" of the annihilated much bigger parent cluster to be of the order of R sqc,thin (see Section 2.2).In this case, a cluster's birth radius can be obtained from n csc R 3 sqc,thin /n csc R 3 sqc,b = η B which gives R sqc,b ∼ 10 5 fm/n csc,39 .Some constraints and implications to consider in the future include: (i) Annihilation should also yield pions.These would decay on weak-interaction timescales and if their mean-free-path turns out to be much smaller than that of photons they may affect the ∼ 10 2 fm SQCs; (ii) While cluster formation (when the universe has aged such that its temperature is in the tens of MeV) is followed rapidly by annihilation, we must avoid re-creation of matter-antimatter pairs.I.e.ensure that pair-creation timescales exceed the Hubble expansion timescale; (iii) In the framework we outline here, SQCs would require a formation mechanism which is different from that of the much larger A >> A qc cosmic strange-quark nuggets [14] (which require a first-order QCD phase transition) and Axion quark nuggets [15] (which require their co-existence with Axions).Additionally, the mechanism by which these nuggets can convert their stored gluonic energy into ∼ 2 MeV photons is a serious limitation.

Conclusion
We have proposed that a color superconducting (CSC) phase of lukewarm QCD matter could offer a non-exotic solution to the cosmological 7 Li, the CDM and the Dark Energy enigmas.The narrow 2 MeV photon line which destroys 7 Be in the radiationdominated post-BBN epoch we attribute to gluonic condensation (i.e.light glueballs or LGBs) in the CSC phase and its electro-magnetic decay modes (Section 2).The detailed properties of the CSC phase remain to be scrutinized although a neutral 2SC-like phase with a superconducting gap ∆ < 0.1µ is hinted at.
CDM, according to our model, consists of colorless, charge neutral, optically thin cosmic quark clusters in the CSC phase (SQCs) with R sqc ∼ 100 fm in size and baryon number A sqc ∼ 10 6 .They decouple from hadrons and interact only gravitationally thus evading detection in current DM experiments.If SQCs could be produced in experiments such as the Compressed Baryon Experiment at FAIR (see bullet point #5 in Section 5), they could be detected via the MeV photons from LGB → γ + γ decay giving support to our model.
The decoupling of SQCs from hadrons is due to leakage of the in-SQCs vacuum into the trivial vacuum of the exterior space-time which yields DE.As leakage proceeds, our cosmology gradually transitions from a non-DE to a DE (ΛCDM-like) universe at moderate redshift while allowing for a possible resolution of the Hubble tension (Section 3).It is crucial to recognize the significance of having optically thin SQCs with a size on the order of 100 Fermi in resolving the 7 Li problem.This requirement directly sets the timescale for the tunnelling process to be on the order of a Giga-year.The connection between the Fermi scales, which govern particle interactions, and the astrophysical scales suggests a fundamental connection between the 7 Li problem and Dark Energy (DE) unique to our model.
Our model does not introduce new physics to solve the 7 Li, DM and DE problems but instead makes use of still uncertain properties of QCD phases and its vacuum properties.We are not the first to discuss a connection between QCD vacuum and cosmology.It has been argued based on empirical properties of hadrons that confinement is a pre-requisite for retaining condensates inside hadrons, which then largely eliminates the problem of the smallness of the cosmological constant.On the other hand, our model, introduces the concept of tunnelling of the in-SQC vacuum into the exterior trivial vacuum, and the decoupling of SQCs from hadrons (see Section 3).This distinctive proposition raises the intriguing possibility that DM and hadrons could represent separate phases of quark matter within the framework of QCD, characterized by distinct vacuum properties that may turn out to have other useful physical applications.Cooper pairs due to the attractive interaction mediated by the strong force.The Cooper pairs in the 2SC phase are made of the red and green quarks only, and the blue quarks do not participate in pairing.Only five out of a total eight gluons become massive while the other three gluons remain massless and do not interact with the paired red and green quarks nor with the gapless blue quarks.These 3 gluons obey their own gluon-dynamics and form their own condensate made of bound states of pure gluons similar to how mesons are bound states of quarks and antiquarks.This gluon condensate is unstable electromagnetically and emits photons.The decay process of the electrically neutral condensate into photons occurs through its interaction with virtual quark loops, which carry the electric charge (see Appendix C).This electromagnetic decay plays a crucial role in resolving the cosmological Lithium-7 problem within our model, as explained in the following description.
The second key assumption in our model is that these QCs transition into a color superconducting state (CSC) in the post-BBN era, when the universe is a few hours old, at a redshift of z <∼ 2.5 × 10 7 .The 2-flavor color superconductivity (2SC) phase specifically involves two flavors of quarks, typically up and down quarks as seen in QCs.When QCs reach his phase and become superconducting (referred to as SQCs), a percentage of their gluons form Light Glueballs (LGBs).Glueballs are bound states of gluons, analogous to how mesons are bound states of quarks and antiquarks.
LGBs are considered to be relatively light, with a mass in the MeV range, compared to other possible glueball states with masses in the GeV range.LGBs are electromagnetically unstable, and by appropriately selecting 2SC parameters, they can decay into a monochromatic ∼ 2 MeV line, capable of selectively destroying primordial 7 Be while leaving other BBN products and CMB physics unaffected.
To resolve the 7 Li problem, approximately 20-30% of the gluonic energy in SQCs needs to be converted into a monochromatic 2 MeV line.Interestingly, the total mass of SQCs required for this solution corresponds to the observed DM in the universe.SQCs, being cold, colorless, charge-neutral, and optically thin to all photon energies, present an attractive DM candidate that is deeply intertwined with the enigma of 7 Li.We would like to emphasize the necessity of a narrow band emission (preferably a mono-energetic photon source) to optimize the number of photo-dissociated 7 Be nuclei (see Eq. ( 1)).Otherwise, there will not be enough gluonic energy from the SQCs (the DM in our model) for the 2/3 reduction.The LGB decay channel is an attractive mechanism for such a line (see Appendix C) and may be tested in current heavy-ion collision experiments (see bullet point #5 in Section 5) The nature of the QCD ground state continues to be a topic of active research and debate.There are many possible forms of the ground state wavefunction (in a many-body variational approach).Even when comparing the many ground state using free energies one can never be certain that the true ground state is found (see discussion in [44]).In principle, the vacuum state within SQCs after the formation and decay of LGBs can potentially exhibit different symmetries compared to the vacuum state of confined quark matter in hadrons; see Appendix C.1.This brings us to the third assumption in our model which has implications for DE.By assuming that the in-SQC QCD vacuum becomes unstable following LGB decay, we allow it to "leak" into the trivial vacuum of the surrounding space-time via quantum tunnelling (Appendix C.1).This "leakage" process behaves as a cosmological constant and occurs late in the evolution of the universe (at a redshift z tun. of a few units), leading to a transition to a ΛCDM-dominated era (see Figure 4).This "leakage" effectively separates the pre-ΛCDM and ΛCDM phases and offers a potential resolution to the Hubble tension (see Section 3).
Through our efforts to resolve the 7 Li problem, we have uncovered a solution that effectively tackles the challenges of DE and DM, thus revealing the intrinsic connection of these critical cosmological issues within our proposed model.For example, the in-SQC vacuum would be stable to decay if its gluonic sector were stable (as in hadrons).The cosmological 7 Be nuclei act as absorbers of LGB decay products (with the 2 MeV photons ideal in this scenario) and avoid affecting BBN and CMB properties.Without addressing the 7 Li, DM, and DE problems simultaneously, our model would not hold together.Despite its assumptions and limitations, it is important to reiterate that our framework operates within the standard model of particle physics but in an unexplored regime of QCD, emphasizing the need for further investigation and empirical validation.

Appendix B Post-BBN 7 Be destruction
In Poulin & Serpico [11] (and references therein), the photon-generating particle X decays exponentially on timescale τ X and forms at early time in the universe at t (<< τ X ).In our model, the X particle (i.e. the LGB) forms and decays to a narrow E 0 ∼ 2 MeV photons nearly instantaneously in the radiation-dominated post-BBN era at time t G when the quark cluster enters the CSC phase; here "G" stands for gluons as a reference to the symmetries and other (e.g.Poincaré) symmetries remain the same during the transition from deconfined-to-confined SU c (2), the vacuum looses the Z 2 -symmetry associated with the center of the SU c (2) group and "chooses" a specific direction or configuration.The vacuum of the confined SU c (2) is in principle different from that of the unconfined SU c (2) and it is not unreasonable to assume that it may be metastable.
Z 2 -symmetry breaking and false vacuum metastability and tunnelling into the true vacuum with consequences to cosmology have been discussed in the literature (e.g.[52] and references therein).In our case, the phenomenon of tunnelling or leaking could be attributed to two potential factors: (i) the zero-momentum loop at the origin of zero-point energy at zero temperature; (ii) the gluon condensate if we allow for a small temperature that is insignificant compared to any other scale in the problem.Which mechanism is exactly at play is currently uncertain within the scope of our model and what is presented here serves as preliminary outlines for potential future research directions.The main point is that the leakage suggests a mechanism to explain Λ-CDM cosmology in our model.and f V values (see examples in Figure 4) with a leakage characteristic redshift 2 < z tun.< 6 (i.e. 1 < t tun.(Gyr) < 4); see Section 3.

Figure 2 .
Figure 2. A possible QCD phase diagram; temperature versus quark chemical potential.The dashed curve ("1 to 2") depicts the cooling path of the quark cluster (formed at "1" with T QCD ∼ 150 MeV) traversing the unpaired phase and entering the CSC phase, at keV temperature, from the low density phase ("2").A first-order line separates the unpaired phase from the CSC phase releasing heat (the vertical "2 to 3" arrow).A SQC gets heated to a temperature not exceeding the LGB melting temperature k B T sqc,m ∼ M LGB c 2 (see Section 2.1).The ∼ 4 MeV LGBs decay to the 2 MeV narrow photon line via LGB → γ + γ. "1 ′ " shows the initial state of a NS core making its way to "2" via cooling and compression (see Section 5 and Appendix C).

Figure 4 .
Figure 4.The Hubble constant H 0 as a function of z tun.(the leakage characteristic redshift) in our model for different values of η G and f V .The H 0 = 73 km s −1 Mpc −1 value is shown as the horizontal line.The resulting age of the universe is shown for each case and is younger than the ΛCDM universe.

Funding:
R.O., D.L. and N.K. acknowledge the support of the Natural Sciences and Engineering Research Council of Canada (NSERC).P.J. is supported by a grant from the Natural Science Foundation PHY-2310003.Data Availability Statement: No new data were generated or analyzed in support of this research.

Figure A1 .
Figure A1.A simplified QCD phase diagram in the density-temperature plane; the density is represented by the quark chemical potential µ q .The hadronic phase exists at temperature T < T QCD ∼ 150 MeV and chemical potential µ q < µ h ∼ 350 MeV.The dashed line in the diagram represents unexplored phases that could potentially exist within the QGP phase before reaching the CSC phase.