Evidence for pulsars metamorphism and their possible connection to black holes and dark matter in cosmology

Pulsars and neutron stars are generally more massive than the Sun, whereas black holes have unlimited mass-spectrum, though the mass-gap between 2 - 5$M_\odot$, which applies for both classes, is evident and remains puzzling. Based on the solution of the TOV equation modified to include a universal scalar field $\phi$ at the background of supranuclear densities, we claim that pulsars must be born with embryonic super-baryons (SBs), that form through the merger of individual neutrons at their centers. The cores of SBs are made of purely incompressible superconducting gluon-quark superfluids (henceforth SuSu-fluids). Such quantum fluids have a uniform supranuclear density and governed by the critical EOSs for baryonic matter $P_b = \mathcal{E}_b$ and for $\phi-$induced dark energy $P_\phi= -\mathcal{E}_\phi.$ The incompressibility here ensures that particles communicate on the shortest possible timescale, superfluidity and superconductivity enforce SBs to spin-down promptly as dictated by the Onsager-Feynman equation, whereas their lowest energy state grants SBs lifetimes that are comparable to those of protons. The extraordinary long lifetimes suggests that conglomeration of SuSu-objects would evolve over several big bang events to possibly form dark matter halos that embed the galaxies in the observable universe. Having pulsars been converted into SuSu-objects, which is predicted to last for one Gyr or even shorter, then they become extraordinary compact and turn invisible. It turns out that recent observations on the quantum, stellar and cosmological scales remarkably support the present scenario.


SUPERFLUIDITY IN PULSARS
Pulsars and NSs are considered to be made of superfluids governed by triangular lattice of quantized vortices as prescribed by the Onsager-Feynman equation: v · dl = 2π m N. v, dl, , m here denote the velocity field, the vector of line-element, the reduced Planck constant and the mass of the superfluid particle pair, respectively. Accordingly, Crab pulsar should have approximately Nn = 8.6 × 10 17 neutron and Np ≈ 10 30 protonvortices (Fig. 1). Let the evolution of the number density of vortex lines, nv, obeys the following advectiondiffusion equation: where t, u f , νt denote the transport velocity at the cylindrical radius r = r f and dissipative coefficient in the local frame of reference, respectively. When νt = 0, then the radial component of u f in cylindrical coordinates reads: u max f ≈ −(Ω/Ω) r > 0. In the case of the Crab; this would imply that approximately 10 6 neutron vortices must be expulsed/annihilated each second, and therefore the object should switch off after 10 6 or 10 13 yr, depending on the underlying mechanism of heat transport (see Baym 1995;Link 2012, and the references therein). However both scenarios are contrary to observation, as numerous NSs has been found, which are older than 10 6 yrs, though non of them is older than 10 9 yrs. On the other hand, recent numerical calculations of superfluids reveal generation of large amplitude Kelvin waves that turn superfluids turbulent (see Baranghi 2008;Baggaley & Laurie 2014;Dix & Zieve 2014, and the references therein). It is therefore unlikely that trillions of Kilometer-long neutron and protons-vortices inside pulsars and NSs would behave differently. In this case, u f should be replaced by a mean turbulent velocity < u f > t with u max f being an upper limit 1 . As the number of vortex lines decreases with time due to emission of magnetic dipole radiation and therefore the separation between them increases non-linearly, it is reasonable to associate a timedependent turbulent length scale t(t), which covers the two limiting cases: t(t = 0) = 0 ≈ 10 −3 cm and t(t = ∞) = ∞ = R . This yields the geometrical mean < t >= √ 0 ∞ ≈ O(10) cm. Putting terms together and using νtur =< t >< u f > t to describe the effective turbulent viscosity, we obtain an upper limit for the global diffusion time scale: τ dif f = R 2 N S /νtur = O(10 9 ) yr. Similarly, a comparable time scale for the Ohmic diffusion in this turbulent medium can be constructed as well. This is in line with observations, which reveal that most isolated luminous NSs known are younger than 10 9 yr (see Espinoza et al. 2011, and the references therein). On the other hand, as stable degenerate NSs require the density gradient to be negative, then the very central region would be the first to be evacuated from Figure 1. A magnetized pulsar is born with an embryonic rigid-body rotating super-baryon (SB) at its center, which is made of incompressible gluon-quark superfluid. By interacting with the ambient medium, the SB expels certain number of vortex lines that are absorbed by the surrounding dissipative medium, thereby causing to promptly spin up: the exact quantum signatures of the glitches observed in pulsars and young neutron star systems.
vortex lines and all other removable energies that do not contribute significantly to the pressure. Let r f be the radius of the central region. The star cools and loses energy, r f would creep outwards with an average velocity:ṙ f ∼ R /τ dif f ≈ 10 −10 cm/s. As I show in the next sections, the nuclear matter inside r f would undergo a transition from compressible dissipative neutron fluid into an incompressible gluon-quark superfluid phase; the lowest possible energy state. Once r f = R − ( 1), then the object turns invisible.
In analogy with normal massive luminous stars, we expect massive NSs to also switch-off earlier than their less massive counterparts. The reason is that most models of EOSs for NSs predict a correlation of MNS ∝ αs, where αs( . = R S R ) denotes the compact-ness parameter, R , RS are star and the corresponding Schwarzschild radii (Haensel et al. 2007). In the extremal case however, when αs = 1, it is reasonable to expect that there will be no turbulence to dissipate. In this case the expression of for the global turbulent diffusion time scale should be modified as follows: In terms of NS-parameters, the turbulent viscosity reads: νtur = R 2 N S (Ω/Ω) ( 0/RN S ). Putting terms together, we obtain τ dif f = (1−αs) 2 (RNS/ 0)(Ω/Ω). For a Velatype pulsar with αs ≈ 2/3, the effect of compactness would shorten τ dif f by almost one order of magnitude.

THE INCOMPRESSIBILITY OF SUPRANUCLEAR DENSE FLUIDS
Modeling the internal structure of cold NSs while constraining their masses and radii to observations, would require their central densities to inevitably be much higher than the nuclear density −ρ0: a density regime in which all EOSs become rather uncertain and mostly acausal (see Hampel & Fischer 2012, and the references therein). On the other hand, at very high densities, almost all EOSs converge to the stiffest EOS: Camenzind 2007),where fluids become purely incompressible. The corresponding chemical potential here reads: whose solution is: Particles obeying this EOS communicate with the speed of the light. This implies that the number density must be upper-bounded by ncr, beyond which local thermodynamical quantities become constants, specifically µ = E φ = P local = const. In this case calculating the degenerate pressure from the local quantities becomes invalid as it yields a vanishing local pressure (P local = n 2 ∂ ∂n (E b /n) = 0). Even if the fluid were weakly-compressible only, ∇P local would be too smooth to hold the object gainst its own self-gravity. This becomes even blatant at the center, where the gradient of the pressure vanishes to meet the regularity condition. The usual adopted strategy to escape this pressure-deficiency is to enforce an unfounded inward-increase of E b as r → 0, resulting therefore in unreasonably large central densities.
Alternatively, one could argue that the strong gravitational field would enforce fusing of neutrons at the center of NSs and form a sea of gluon-quark fluid. However, a phase transition from normal nuclear matter into a gluon-quark plasma (GQP) under classical conditions was found to be unlikely (Baym & Chin 1976;Chapline & Nauenberg 1977). The GQPs in these studies were assumed to be locally conformal, compressible and governed by a non-local and constant bag pressure.
However, the situation may differ if there is a mechanism at the background of supranuclear densities, that among others, is capable of injecting energy into the system. This in turn would effectively increase the mass of the object, steepen the curvature of the embedding spacetime and compress the nuclear fluid up to the saturation level, beyond which it becomes purely incompressible. The underlying conjecture here is that compression would mainly affect the distance between individual neutrons, n, rather than the distance, q , between the quarks enclosed in individual neutrons.
In fact there are at least two observational facts that in favor of this argument: (i) Massless gluons mediate the strong force between quarks with the speed of light, whereas global compression of the central fluid-core proceeds almost in a quasi-stationary manner. Hence, as the quark contribution to baryon mass is negligibly small, compressibility of the gluon-quark plasma would violate the causality principle.
(ii) The quark potential Vq inside baryons increases with the radius to attains maximum at the boundaries -r b . Quark confinement however requires that Vq(r b ) ≈ 0.94 GeV. In this case the effective/reactive velocity at the boundary would read: < V >= ∆x/∆t = ∆x ∆E/h ≈ c, where we set ∆x = r b = 10 −13 cm. Consequently, the reaction time of GQPs inside neutrons to whatsoever external compressional effects is of order 10 −24 sec: the shortest possible time scale in the entire system. Moreover, recalling that the lifetime of neutrons in free space is extremely short compared to that of protons, then compression of GQPs inside neutrons most likely would provoke a runaway decay into protons and would give rise to neutrino-dominated electromagnetic eruptions. However, such events can be safely ruled out by observations. (iii) Recent RHIC and LHC experiments have revealed that the matter resulting from smashed protons behaves like fluids, whose constituents move collectively and relatively slowly. These observations are in line with the basic properties of incompressible fluids rather than randomly moving particles in plasmas. Moreover, the coupling constant in QCD is inversely proportional to the density: αS ∼ 1/ρ. Hence matter with ρ > ρ0 and T = 0 yields automatically αS 1. This means that the interaction strength is at low- , assuming the fluid to be incompressible (i.e., the internal energy density E in = const.) for a compactness parameter αs( . = R S /R ) near the critical value 8/9. Obviously, the solution diverges as the center is approached, where the matter becomes ultrabaric. The second profile (Pnum) corresponds to a numerical integration of the TOV equation using E in = const.) and assuming P (r = 0) = E in . est and therefore the quarks must be moving freely, practically unaffected by whatsoever external forces.
Consequently, as GQPs inside neutrons are incompressible, we expect nuclear fluid inside SBs to be indifferent. This implies that, when the distance between neutrons becomes critically small, e.g. n/ q = O(1), then the pions, π 0 , the carriers of the residual strong force between neutrons, become sufficiently energetic and overcome the repulsive barrier, where they go through a "gluonization" process, which subsequently enables them to mediate the strong force between quarks inside SBs efficiently (Fig. 3).
Unlike EOSs in compressible normal plasmas, classical EOSs in incompressible superfluids are nonlocal. In the latter case, constructing a communicator that merely depends on local exchange of information generally would not be sufficient for efficiently coupling different/remote parts of the fluid in a physically consistent manner. A relevant example is the solution of the TOV-equation for classical incompressible fluids (E = const.). In this case, the pressure depends, not only on the global compactness of the object, but it becomes even acausal whenever the global compactness is enhanced (see Fig. 2). This is similar to the case when solving the incompressible Navier-Stokes equations, where an additional Laplacian for describing the spatial variation of a nonlocal scalar field is constructed to generate a pseudo- Figure 3. A Feynman-like diagram of two interacting neutrons at the center of a pulsar (time runs from left-to-right). These two neutrons, that are made of two "d" and "u" quarks approach each other under the effect of the gravitational and the scalar fields. Once the distance between them falls below a certain critical value, the residual strong force undergo a "gluonization" procedure necessary for efficiently mediating the strong force between quarks. Here the two particles merge together and form a super-baryon. pressure (; actually a Lagrangian multiplier), which, again, does not respect causality (Hujeirat & Thielemann 2009).

THE ONSET OF SUPERFLUIDITY
In the presence of magnetic torque and in the absence of energy generation at the center of pulsars, all types of energies, including magnetic, rotational and rest of thermal energies must diffuse outwards. This energy loss enforces the degenerate core to cool down to T = 0 and to increase its density. In this case the corresponding de Broglie wave length, (λDe ∝ 1/T ), which coincides with the radius of the core, increases and the enclosed particles start occupying the same quantum state, hence entering the superfluid phase with vanishing entropy: (S = kBln(Ω) = 0). In fact this is the lowest energy state possible in the system, which is termed here L − State.
Indeed, recent RHIC and LHC-experiments reveal that GQPs inside smashed protons, which are generally at much higher energy states than the states considered here, were found to be nearly perfect with even smaller viscosity over entropy ratio (Shuryak 2017).
Consequently, we predict the interiors of SBs, which result from fusion of individual neutrons in the L − State, to be in a purely superfluid and superconducting phase.
The injected dark energy by φ may affect the central nuclear fluid in two way: it first strengthens the curvature of the embedding spacetime, enhancing compression of the central fluid, and secondly, it provokes the "gluonization" procedure: where ∆E φ , ∆EG, π 0 , ∆E SB g denote the local energy enhancement by φ and "G" ( . = gravitational force), the pion energy and the resulting excess of energy gained by the SB, respectively. This reaction is assumed to proceed directly and silently without destroying the superfluidity character of GQ-fluid inside SBs.
On the other hand, the rotational torque of the dissipative ambient medium exerted on the SBs would enforce them co-rotate uniformly. However, SBs are quantum identities and therefore must obey the laws of quantum mechanics. This means that their dimensions and energies can accept discrete values only. It turns out that SBs evolve in a discrete manner as prescribed by Onsager-Feymann equation in superfluids. SBs here must eject a certain number of vortex lines in order to increase their dimensions Hujeirat (2017). Following this scenario, the internal structure of pulsars and UCOs looks as follow: (a) Pulsars are born with SB-embryos at their centers, which are made of an incompressible superconducting gluon-quark superfluid, (b) a dissipative neutron fluid that surrounds the SB,and (c) a geometrically thin boundary layer inbetween(BL), where the residual of the strong nuclear force becomes dominant over the viscous forces. The viscosity of the neutron fluid here decreases strongly inwards and vanishes at the boundary of the SB. The BL here is practically the zone, where the neutron fluid is prepared to match the physical conditions governing GQ-fluids inside SBs prior to their merger with other neutrons. It turns out that these merger events are identical to the sudden glitch phenomena observed in pulsars and young neutron star systems Hujeirat (2017). Accordingly, once the rotational frequency of the ambient medium falls below a certain critical frequency Ωc, the SB then undergo a sudden spin-down to the next quantum-allowable frequency, thereby expelling a certain number of vortices, which are then absorbed by the ambient medium and causes the observed spin-up of crusts of pulsars (see Fig. 4). However, the excess of rotational energy is then viscously-redistributed into the entire ambient dissipative medium: a process that may last for weeks or even months. Indeed, very recent observations appear to confirm our scenario. Eya et al. (2017) found that the mean fractional moment of inertia in the glitching pulsars correlates weakly with the pulsar spin, implying therefore that glitches are provoked by a central core of different physical properties. Also Serim et al. (2017) discovered for the first time that the short X-ray outburst 25 days prior to the glitch in the accretion-powered pulsar SXP 1062 did not alter the spin-down of the source, which indicates that glitches are triggered by the core of the NS rather than by the outer shells.

CROSSOVER PHASE TRANSITION
For ρ > ρ0, short-range interactions between particles mediated by the exchange of vector mesons most likely would enhance the convergence of the EOSs towards P → E b ∼ n 2 b (see Haensel et al. 2007;Camenzind 2007, and the references therein). The chemical potential here: µ( . = (E b + p)/n b ) increases linearly with the number density n b . Matter with d µ/d n > 0 is classified here as H-State and depicted in red-color in Fig. (5). However In the absence of energy generation, degenerate matter cannot hold a correlation of the type µ ∼ n b indefinitely and it must terminate at a certain critical density ncr. This agrees with the two facts: (i) central densities in NSs, ρc, increase with their masses and (ii) ρc must be upper-bounded by ρc 12, 5 × ρ0 in order to fit the observed mass function (see Lattimer 2011, and the references therein).
On the other hand, in an ever expanding universe, the eternal-state of matter should be the one at which the internal energy attains a global minimum in spacetime (zero-temperature, zero-entropy and where Gibbs energy per baryon is lowest, i.e. in the L-State). Taking into account that µ(r)e V(r) = const. inside massive NSs together with the a posteriori results e V 1 (see Fig. 9), we conclude that µ ∼ E b /n b = const. Under these conditions the gradient of the local pressure P local vanishes and a non-local pressure PNL is necessary to avoid self-collapse of the object into a BH. 2 The slop of the chemical potential dµ/dn between the H and L-states in the n-space could attain positive, negative or even discontinuous values. However the case with dµ/dn > 0 should be excluded, as it implies that the eternal state of matter would be more energetic than the H-State, which is a contradiction by construction. Similarly, the case dµ/dn < 0 is forbidden as it would violate energy conservation (; dµ/dn < 0 ⇔ dP/dn < 0 ⇔ adding more particles yields a smaller pressure). Moreover, let us re-write the TOV equation in terms of µ: Obviously, as µ > 0, a negative dµ/dn would destabilize the hydrostatic equilibrium, unless external sources are included, e.g. bag energy and/or external fields. Therefore, although a first order phase transition may not be completely excluded, a crossover phase transition into an incompressible superfluid phase with µ = µ(n = ncr) = const. would be more likely.
Here, µ and P on both sides of the transition front are equal and, with the help of an external field, both (Etot/n) + and (Etot/n) − across the front can be made continuous (Fig. 2)).
In the present study, the crossover phase transition corresponds to "silent" mergers of baryons forming SBs. Inside SBs, the fluids have µ = n = const. and T = 0 and therefore vanishing entropy. To overcome the repulsive barrier between individual neutrons and provoke their merger, external forces are required. In the present study, we assume that, in addition to compression by self-gravity, there is a universal scalar field, φ, at the background of supranuclear dense matter, which provide the energy required for forming an eternal stable GQ-cloud inside SBs. In this case, the dark energy potential associated with φ should have similar effects as the generalized quarkpotentials inside individual neutrons, namely: where a φ , b φ , c EM φ are constant coefficients (see Sakumichi & Suganuma 2015, and the references therein). Without loss of generality, we may set c EM φ = 0 and γ = 2. As the scalar field is universal, its special and temporal variations are assumed to vanish. Therefore, 2 Classically incompressible fluids with E = const. have negative local pressure. Therefore an acausal non-local pressure is usually used for stabilizing the fluid configurations modelled by the incompressible Navier-Stokes equations. Figure 5. A schematic description of the chemical potential µ, µ ( . = ∂µ/∂n) versus the number density n at the centers of NSs. At a critical central density, ncr, the universal scalar field φ is set to provoke a crossover phase transition from compressible and dissipative neutron fluid into incompressible GQ-superfluid. The transition front creeps from inside-to-outside to reach the surface of the object on the scale of one Gyr.
the expression E φ = 1 2φ 2 + V (φ) + 1 6 (∇φ) 2 reduces to E φ = V (φ) = V φ . By requiring the chemical potential at r = 0 to be upper-limited by the deconfinement energy of quarks inside a single neutron, then a Gibbs function of the following type can be constructed: Substituting E b = a0n 2 and V φ (r) in Eq. (5) at r = 0, then f (n) reduces to: The Gibbs function here may accept several minima at nmin = (b φ /a∞) 1/2 , though f (n = nmin) doesn't necessary vanish. However f (n) > 0 and f (n) < 0 should be excluded, as they are energetically not suited for a smooth crossover phase transitions.
On the other hand, by varying a∞ and b φ , a set of minima could be found. One way to constrain b φ is to relate it to the canonical energy scale characterizing the effective coupling of quarks, i.e. b φ = 0.221 (see Bethke 2007, and the references therein). Indeed, as shown in Fig (3), f (n) attains a zero-minimum at n ≈ 3 n0 for a∞ = 1.0.

THE GOVERNING EQUATIONS AND THE SOLUTION PROCEDURE
The investigation here is based on numerical solving the TOV equation modified to include scalar fields (φ). The modified stress energy tensor reads: The superscripts "0" and "φ" correspond to baryonic and scalar field tensors: Uµ here is the 4-velocity, the subindices µ, ν run from 0 to 3 and gµν is a background metric of the form: gµν = e 2V dt 2 − e 2λ dr 2 − r 2 dθ 2 − r 2 sin 2 θd 2 ϕ 2 , (10) where V, λ are functions of the radius.
Assuming the configuration to be in hydrostatic equilibrium, then the GR field equations, Gµν = −8πG Tµν reduce into the generalized TOV equations: where m(r) = 4π Er 2 dr is the total enclosed mass: here denotes the interaction potential of the scalar field with the baryonic matter, i.e., the rate at which dark energy is injected into the system andφ is the time-derivative of φ.
Our reference object is a NS with 1.44 M , a radius R = 2 × RS, where RS is the Schwarzschild radius. The universal scalar field, φ is assumed to be spatially and temporarily constant, whereas V φ is taken from Eq.(4). The non-dimensional value of a φ is chosen to fulfill the a posteriori requirement: R = RS + , where 0 1. In most of the cases considered here, b φ is set to be identical to the canonical energy scale at which momentum transfer between quarks saturates, i.e., b φ = 0.221 GeV. The fluid in the post transition phase is governed by the EOS: P 0 = PL = E 0 = ρcrc 2 = const. For a given central density, the solution procedure adopted here is based on integrating the TOVequations for the pressure, enclosed mass and for the pseudo-gravitational potential from inside-to-outside, using either the first order Euler or the fourth order Runge-Kutte integration methods. For verification purposes the TOV equation has been solved using a polytropic EOS as well as for incompressible case (E b = const.), using both integration methods (Fig. 7). (8) and (9)/top we show the profiles of the negative pressure, the enclosed effective mass (i.e. the energy due to baryon and energy enhancement by the scalar field φ) and the inverse of the metric coefficient grr for different values of the dark energy coefficient a φ (see Eq. 4). Obviously, assuming the object to have αS = 1/2 initially, its final compactness must then increase with increasing a φ , i.e. with increasing the rate of dark energy injection. Using the limiting values a φ = 0.78 and b φ = 0.22, the scalar field is capable of injecting the energy required to deeply sink the object into spacetime and turn it invisible, where α cr S = 1 − , hence attaining its maximum possible value. Note that must be extremely small, but still finite in order to prevent the collapse of the object into a BH. These values of a φ and b φ yield a final effective mass that is twice as the initial baryonic mass.

Similarly, in Figures
Once the object has completely metamorphosed into a stellar-size SB and its compactness attained the critical limit α cr S , then the spatial variation of the metric coefficient grr across its surface becomes nearly singular (Fig. 9), which implies that all sorts of surface radiation will be extraordinary redshifted, hence the object becomes observationally indifferent from its BH-counterpart. On the other hand, the metric coef-ficient gtt, which is a monotonic function of the gravitational potential, appears to be fairly flat inside the object, but it undergoes a dramatic change across the surface to finally attains its weak-field values outside the object.
In addition, we have solved the TOV equation using two different EOSs: an EOS that corresponds to an incompressible GQ-superfluid inside a gradually growing SB together with a polytropic EOS for the ambient medium at different evolutionary epochs (Fig. 10). While the scalar field here is set to inject energy into the SB and therefore increasing its effective mass, the surrounding medium appears to stably and comfortably adjust to the new condition of the SB.
In Fig. (11) we display the Mass of fullydeveloped SuSu-objects and their progenitors versus ρcr superimposed on the locations of NSs as reported by Lattimer (2011). Most remarkable here is that Hulse-Taylor type pulsars are able to form SBs at their centers at much lower central density than usually required for a phase transition into quark fluids, whilst still end up twice as massive as their initial mass (Fig.  11).

SUMMARY & DISCUSSIONS
The here-presented model is motivated by the following issues: • Why neither NSs nor BHs have ever been observed in the mass-range 2 -5 M .
• Most sophisticated EOS used to model the internal structure of NSs are based on central densities that are far beyond the nuclear density: an unknown density regime with great uncertainty.
• What is the origin of glitches observed in pulsars and young neutron stars and wether these carry information that may disclose the internal structure of their cores?
• Could massive NSs end as maximally compact dark objects, i.e. as BH-candidates?
• How does the state of matter in pulsars and NSs evolve on the cosmological time scale and whether they have hidden connection to dark matter and dark energy in cosmology?
In this respect a scenario has been presented, which can be summarized as follows: (i) Pulsars are born with embryos at their centers (here termed super-baryons -SBs). The interiors of SBs are made of GQ-superfluids and governed by the EOS P = E = a∞ n 2 : a purely incompressible fluid state. As a consequence, there is a universal maximum density ρ uni q (= O(ρcr)), where the momentum transfer between quarks saturates, the coupling constant αasym attains its universal minimum, where quarks are moving freely. Recalling that the spatial variation of the coefficient grr of the Schwarzschild metric on the nuclear length scales is negligibly small (dgrr/dl 10 −19 ), the GQ-superfields may not accept stratification by gravitational fields.
(ii) In the presence of a universal scalar field φ at the background of supranuclear densities, the injected dark energy by φ is capable of provoking a phase transition from compressible dissipative neutron fluids into incompressible GQ-superfluids. The effect of both the gravitational and scalar fields is to mainly enhance and convert the residual strong force between nucleons into a strong force that holds quarks together. This action is termed here as a "gluonization" procedure, which is equivalent to energy injection into the system, thereby maximally enhancing the effective mass of the object and turns it invisible. It turns out that using an interaction potential of the type V φ = a φ r 2 + b φ appear to be most appropriate for maximizing the compactness of the object without significantly changing its dimensions (Fig. 5). (iii) In a recent study, I have shown that the glitch phenomena observed to associate the evolution of pulsars and NSs are in perfect-agreement with the formation scenario of embryonic SBs at their centers. (iv) Once the entire object has metamorphosed into a stellar-size SB, the spacetime in its interior would be fairly flat, but become exceedingly curved across its surface (Fig. 5). Hence SBs are practically trapped in spacetime, extraordinary redshifted and therefore completely invisible. (v) According to the here-presented scenario, all visible pulsars and NSs must contain SBs. The gravitational significance of these SBs depends strongly on their evolutionary phase and in particular on their ages and initial compactness. Accordingly, pulsars and young NSs should be less massive than old ones, and the very old NSs should turn invisible by now.
To quantify the mass-enhancement by φ, let M b be the mass of the NS at its birth and M φ being the mass enhancement due to φ. Requiring R > RS, then the following inequality holds: where αDE . = M φ /M b . Etot, ρ15, M1.44 denote the total energy, the density in units of 10 15 g/cc and the baryonic mass of the NS in units of 1.44 M , respectively. Assuming NSs to be born with Etot = E b , then by interacting with φ, they become more massive and more compact to finally reach R = RS + at the end of their luminous phase ( 1), which would last for approximately 10 9 yr or even less, depending on their initial compactness Fig. 4 and Fig. 7).
Similar to atomic nuclei, I conjecture that the enormous surface stress confining the sea of GQ-superfluid inside stellar-size SBs render their surfaces impenetrable by external low energy particles, hence maintaining the eternal stability and invisibility of SuSu-objects. This implies that the core of a fully-developed SB would be shielded by a protecting repulsive barrier. However, due to the deep gravitational wells of SBs, such incidents would be practically observable, though they are ruled out by observations completely. Even if there were no repulsive barriers, we expect SBs to still be stable against mass-enhancement from external sources. Let a certain amount of baryonic matter, δM b , be added to the object via accretion from external sources. Then the relative increase of R compared to RS scales as: δR δR S ρcr ρnew , whereρnew is the average density of the newly settled matter. Unlessρnew ρcr, which is forbidden under normal astrophysical conditions, the SuSu-object would react stably. However, in the case of super-Eddington accretion or merger, the newly settled matter must first decelerate, compressed and subsequently becomes virially hot, giving rise therefore toρnew ρcr. On the other hand, such events would lower the confinement stress at the surface and would turn the quantum jump of the energy density at R , which falls abruptly from approximately E ≈ 10 36 erg/cc at R down to zero outside it, into an extraordinary steep pressure gradient in the continuum. While such actions would smooth the strong curvature of spacetime across R , they would enable SuSu-objects to eject quark matter into space with ultra-relativistic speeds, which is forbidden. Nonetheless, even if this would occur instantly, then the corresponding time scale τ d would be of order Λj/c, where Λj is the jump width in centimeters. Relating Λj to the average spacing between two arbitrary particles (∼ n −1/3 ), this yields τ d ≈ 10 −24 s, which is many orders of magnitude shorter than any known thermal relaxation time scale between arbitrary luminous particles.
Although electromagnetic activities and jets have never been observed in dark matter halos, they are typical events for systems containing black holes. Recalling that supermasive GBECs are dynamically unstable (Hujeirat 2012), our results here address the following two possibilities: • If the onset of φ-baryon interaction indeed occurs at ncr, then the majority of the first generation of stars must have collapsed into UCOs and subsequently became SuSu-objects, rather than collapsing into stellar BHs with M 5 × M . These objects are expected to conglomerate into dark matter halos over several big bang events and to embed the galaxies in the observable universe. This conclusion is in line with recent observations of NASA's Spitzer Space Telescope, which reveal the existence of primordial galaxies, such as GN-z11, whose age might be even bigger than that of our universe (?).
• The passivity of DM to electromagnetic radiation may indicate that the collective effect of the repulsive forces of a cluster of SuSu-objects is repulsive on smaller length scales and attractive on the larger ones. Hence approaching luminous matter will be forced to deviate from face-to-face collisions with the cluster of SuSu-objects, though n-body and SPH-numerical calculations are needed here to verify this argument. Figure 7. For test purposes the profiles of the pressure and enclosed mass of a NS that have been obtained by solving the TOV equation, using a polytropic EOS (P = Kρ γ ). In the lower figure the the numerical errors for two different numerical integration methods are shown: the first order Euler and the fourth order Runga-Kutta methods using 10 4 grid points compared to a reference solution that has been obtained using 10 7 grid points. While the accuracy of Ruga-Kutte method is highly superior over Eulers, the effect of errors on the solutions presented here can be safely neglected. Figure 8. The radial distributions of the baryonic pressure (P bar ) and negative pressure (P φ ) inside an incompressible GQ-superfluid core (left). The enclosed mass of the baryonic matter and the gradual mass-enhancement due to dark energy is shown for different values of a φ (right).

Figure 9.
In the top panel we show the radial distributions of the metric coefficients grr and gtt inside a normal NS (; P L = Kρ γ and P φ = 0) and inside an SuSu-object (; P local = const. and P φ = −V φ ). Obviously, normal models of NSs have larger radii and considerably less compact than their SuSu-counterparts, which can be inferred from the very limited spacial variations of grr and gtt. In the lower panel, the compactness of a typical SuSu-object, expressed in terms of -1/grr is shown for different values of a φ . The object turns invisible if V φ is calculated with a φ = 0.78 and b φ = 0.22. Figure 10. The profiles of the total energy density Etot, the pressure P φ induced by φ and the combined pressure Ptot versus radius are shown for different evolutionary epochs τ 1 < τ 2 < τ 3 . Inside r f : P L = 0 and P φ = −V φ , whereas outside r f : P L = Kρ γ and P φ = 0. In each epoch, the object has an SB-core overlayed by a shell of normal compressible matter obeying a polytropic EOS. Obviously, the object appears to stably and comfortably adjust itself to the mass-redistribution inside r f , where the matter is in incompressible GQ-superfluid state. Figure 11. Upper mass limit of SuSu-objects versus critical density ncr (in units of n 0 ) is shown. The φ−baryon interaction is set to occur at ncr, which in turn provokes the phase transition into the incompressible GQsuperfluid state. The most probable mass-regime for SuSuobjects is marked here in blue colour. Accordingly, the progenitor of a SuSu-object with 3.36 M should be a pulsar/NS of 1.68 M , provided it has an initial compactness α S = 1/(1 + α DE ) = 1/2 and ncr = 3 n 0 . Similarly, a Hulse-Taylor type pulsar would end as an SuSu-object of 2.91 M , if its initial compactness was α S = 1/2 and if ncr = 4 n 0 . On the other hand, moderate and massive NSs with initial compactness α S 2/3, i.e., α DE 1/3, need less dark energy to become invisible SuSu-objects, though an unreasonably high ncr is required for triggering φ−matter interaction. NSs falling in this category are to be compared with the colored small cycles and triangles, which show the approximate locations of various NS-models as depicted in Fig. (4) of Lattimer & Prakash (2011). Figure 12. A schematic description of an invisible SuSuobject trapped in spacetime and surrounded by a repulsive barrier that protects/confines the enclosed sea of GQsuperfluid. The spacetime inside SuSu-objects is fairly flat, but extraordinary curved and nearly singular across their surfaces.