_{1}

^{*}

The observed recurrence of glitches in pulsars and neutron stars carries rich information about the evolution of their internal structures. In this article, I show that the glitch-events observed in pulsars are exact quantum signatures for their metamorphosis into dark super-baryons (SBs), whose interiors are made of purely incompressible superconducting gluon-quark superfluids. Here the quantum nuclear shell model is adopted to describe the permitted energy levels of the SB, which are assumed to be identical to the discrete spinning rates
Ω_{SB} that SBs are allowed to rotate with. Accordingly, a glitch-event corresponds to a prompt spin-down of the superconducting SB from one energy level to the next, thereby expelling a certain number of vortices, which in turn spins up the ambient medium. The process is provoked mainly by the negative torque of the ambient dissipative nuclear fluid and by a universal scalar field
∅ at the background of a supranuclear dense matter. As dictated by the Onsager-Feynman equation, the prompt spin-down must be associated with increase of the dimensions of the embryonic SB to finally convert the entire pulsar into SB-Objects on the scale of Gyrs. Based on our calculations, a Vela-like pulsar should display billions of glitches during its lifetime, before it metamorphoses entirely into a maximally compact SB-object and disappears from our observational windows. The present model predicts the mass of SBs and
Δ
Ω/Ω in young pulsars to be relatively lower than their older counterparts.

Besides other puzzling activities, the glitch-phenomena observed in pulsars and young neutron stars are considered to be the key for decoding their internal structures. The energies associated with their prompt spin-up are vast and in most of the cases beyond 10^{40} erg/event [

Although the state of matter in this exotic density regime is neither clear nor verifiable, a transition into a quark phase under normal conditions is considered to be theoretically unlikely [

In fact, the rest energy of pentaquarks detected by the LHC was found to lay between 4.38 - 4.45 GeV [

Recalling that the bulk mass of a baryon mass originates mainly from the gluon-cloud embedding the quarks, then the number of gluon-channels connecting 10^{57} quarks inside a stellar-size super-baryon would increase dramatically (

indefinitely and it must be upper-bounded by μ 0 , beyond which the fluid acquires a constant μ and therefore a constant “ n ⌣ ”; practically becoming a purely incompressible nuclear fluid [

E = a 0 n ⌣ 2 → darkenergy E = a q s f × n ⌣ , (1)

where a , a q s f are constant coefficients.

In a pervious article [

In

where neutrons, under the effect of the enhanced gravitational field, start merging together to form an SB at the center of the pulsar (the central sphere in top ^{1} and yields an enhancement of the effective mass of the SB. In the present study, increasing of the number of flux-tubes between quarks is equivalent to injection of dark energy, which accounted for through modifying the EOS in the following manner: E t o t = E b + E ϕ , where ϕ corresponds to the scalar field ϕ .

^{1}The gluon cloud can viewed as a collection of flux-tubes connecting quarks, whose number must increases super-linearly relative to the number of quarks.

Similar to laminar flows over solid static bodies, a thin boundary layer (−BL) between the dissipative ambient medium and the rigid-body rotating SB is predicated to form. The dynamics of the fluid inside the viscous BL, which initially was rotating differentially (due to shear viscosity), would start to be increasingly influenced by the rest of the strong nuclear force in combination with compression-enhancement due to the injection of dark energy. These forces would quickly start to dominate over viscous stresses to subsequently provoke its transition into the superfluid phase. The matter inside the BL would lock to the SB and starts rotating coherently with it.

Our model is based on the assumption that the SB behaves as a single quantum entity, whose permitted rotational frequencies form a sequence of discrete values { Ω S B n } .

Once the rotational frequency of the ambient medium has decreased and hit a suitable element of { Ω S B n } , then the SB undergos a sudden spin-down and starts rotating with this lower value, thereby ejecting a certain number of vortices. In turn, these vortices are absorbed then by the ambient dissipative medium and provoke its prompt spin up, which mimics a single glitch event out of many that have been observed to associate the long-term evolution of pulsars and young neutron stars. Each time the SB undergoes a spin-down, its size and inertia must increase respectively to finally metamorphose the entire pulsar into a stellar-size SB. At that moment, the object should have attained its maximum compactness and would turn invisible completely.

Assume the Vela pulsar to have been born with the initial angular velocity Ω 0 = 1540 / sec and magnetic field B 0 = 10 13 Gauss. If the entire interior of the pulsar were made of a superfluid and its total rotational energy were stored in a network of vortex lines, then we may use Onsager-Feynman equation,

∮ v ⋅ d l = 2 π ℏ m N V to compute the total number of vortex line threading the

entire object. Here N V , v , d l , ℏ , m denote the total number of vortex lines, velocity field, vector of line-element, the reduced Planck constant and the mass of the superfluid particle pair, respectively. Assuming the number density of vortices across the entire object to be constant, then the enclosed number of vortex lines inside a sphere of radius “r” would be N V ( r ) = const × r 2 , which mimics rigid body rotation with Ω = const .

However, this picture oversimplifies the problem and leaves numerous questions open, in particular those related to the origin, dynamics and timing of glitches as well as to the physics and EOSs of supranuclear dense superfluids.

Similar to normal stars, the cores of pulsars are expected to initially consist of normal dissipative fluids. Under the effect of gravity and exotic fields, neutrons at the very central region would start merging together to form an embryonic SB. The inertia of the SB relative to the entire star behaves like I S B / I ⋆ ~ ( R S B / R ⋆ ) 5 . Thus, unless the SB is of a considerable macroscopic size, its dynamics would be completely insignificant.

Inside the rotating SB, the number density of the mini-vortices is constant and therefore the SB can be safely considered as a rigid body rotator (see

Consider a terrestrial ^{4}Helium-superfluid in a container. It was experimentally verified that the superfluid would not change its rotational frequency, unless its rotational frequency has hit a term, which the superfluid is allowed to spin with in accord with the laws of quantum mechanics [

4 Heliumsuperuid ↔ superuidinsidetheSB Solidcontainer ↔ AmbientdissipativemediumsurroundingtheSB (2)

Here the emission of magnetic dipole radiation enforces the ambient medium to spin-down continuously. While spinning down, the ambient medium certainly would hit an element of the Ω c -sequence, whose elements are the values that the SB is allowed to rotate with. In order to lower its spin, the SB must expels instantly a certain number of vortices. However, the ambient dissipative medium would then absorb the expelled vortices and spin-up promptly as illustrated in

To follow the time-evolution of the interactions between the BL, SB and the ambient medium, the Onsager-Feynman equation (OFE) is employed. In the finite and flat-spacetime the time-derivative of the OFE reads:

( S Ω ) n + 1 − ( S Ω ) n Δ t = h 2 m ( N n + 1 − N n Δ t ) ,

where Δ t ( ≐ t n + 1 − t n ) is a time-interval between two successive glitch events. Applying the equation to the step function of Ω c (see

( S Ω ) n + 1 − ( S Ω ) n = [ S Ω ] = ( h / 2 m ) δ N V . (3)

This equation states that the matter in the concerned region is allowed to change its state from ( S Ω ) n into ( S Ω ) n + 1 , if a certain number of vortices, δ N , is extracted out of the system.

Consider now the matter inside the BL (see

Right at the Ω -discontinuity (e.g. t = t n + 1 in

・ As the matter both in the BL and in the SB rotates rigidly with the same Ω n , then the SB-BL system is rotationally behaving as a single quantum entity. In this case the corresponding OFE reads:

S n + 1 Ω n = ( h / 2 m ) N V * ,

where N * = ( S n + 1 / S n ) N n ( > N n ) is the total number of vortices enclosed in the SB-BL system (see sphere (ii) in

・ The matter inside the modified SB, namely inside S n + 1 , starts to interact directly with the surrounding viscous matter, which rotates with the lower frequency Ω n + 1 ; an element of Ω -sequence which the SB is allowed to quantum-mechanically rotate with. Hence the SB spin-down promptly from Ω n into Ω n + 1 , thereby ejecting a certain number of vortices that amounts to:

δ N V = N V * − N n + 1 = [ N v ] = ( 2 m / h ) S n + 1 ( Ω n − Ω n + 1 ) .

The loss of rotational energy here is assumed to be the mechanism that triggers the phase transition of the superfluid neutrons inside the LB into an incompressible gluon-quark superfluid state, where the EOS changes from E = a 1 n ⌣ 2 into E = a 2 n ⌣ and where a 1 , a 2 are constants and n ⌣ is the number density of matter. The expulsion of this energy would certainly affect the dynamical stability of the SB and in principle could provoke a dynamical collapse of the SB. This collapse may still be escaped if the SB enlarges its dimensions in such a manner that its total rotational energy is conserved, i.e.:

E r o t n = 1 2 I S B n ( Ω n ) 2 = 1 2 I S B n + 1 ( Ω n + 1 ) 2 = E r o t n + 1

Recalling that the matter inside the SB has the uniform constant density ρ = ρ c r , then the following relation can be obtained:

S n + 1 = S n ( Ω n Ω n + 1 ) Γ , (4)

where Γ is a constant. Depending on the density of dark energy, Γ may vary from Γ = 4 / 5 in the case of pure baryonic matter up to Γ ≈ 1 in dark energy dominated case.

On the other hand, as Ω n / Ω n + 1 = 1 + ϵ , where ϵ ≪ 1 , I use, for simplicity, the limiting case: S n + 1 = S n ( Ω n / Ω n + 1 ) to predict the size of the SB shortly after each glitch-event.

To conclude: when applying the OFE to Ω -discontinuity at t = t n + 1 , then the equation may be split into the following two parts:

[ S Ω ] = ( h / 2 m ) δ N V ⇒ { i . [ S Ω ] = 0 ii . S [ Ω ] = ( h / 2 m ) [ N V ] (5)

Hence, each time the pulsar undergo a glitch event, the dimensions of the SB should increase and similarly its gravitational mass and inertia. The ambient medium however would experience an opposite development: it would become increasingly smaller, dynamically insignificant to finally vanish completely.

Since the discrete energy levels of an SB embedded in a strongly curved spacetime are completely unknown, I use pulsar’s observation to guide the search for a reliable Ω -sequence, whose elements are likely to be the values that SBs are allowed to rotate with.

In the present study, instead of { Ω c n } , I use for simplicity the sequence { α c n } , whose elements are defined as follows:

α c n = Ω c n / Ω c n + 1 = 1 + ( Δ Ω / Ω ) n ,

where the terms { ( Δ Ω / Ω ) n } n = 1 ∞ are observables that can be used to constraint the α c n -sequence. Moreover, in searching for the best suited sequence, the two limiting cases of the cross section of the SB should be taken into account:

i . S c ( t = 0 ) = S 0 ii . S c ( t → ∞ ) = S ∞ = S ⋆ , (6)

where S 0 , S ∞ denotes the initial and final cross-sections of the dynamically evolving SB, respectively.

Noting that Ω c n / Ω c n + 1 is a slowly varying function with time (e.g. on the scales of years), the following two relation can be obtained:

S c n = ( α c 0 × α c 1 × α c 2 × ⋯ × α c n ) S 0 = ( ∏ i = 0 n α c i ) S c 0 , and Ω c n = Ω 0 / ( α c 0 × α c 1 × α c 2 × ⋯ × α c n ) = Ω 0 / ( ∏ i = 0 n α c i ) . (7)

Using theoretical considerations in combination with glitch observation of pulsars, in particular those of the well-studied Crab and Vela pulsars, a list of constraints may be constructed that can be used to determine the elements of the α c n -sequence:

1) The observed spin-down of pulsars implies that α c n > 1 .

2) The Vela pulsar is roughly 10 times older than the Crab pulsar

3) α c n → 1 + 2 × 10 − 8 as n → N ˜ Crab and α c n → 1 + 2 × 10 − 6 as α c n → 1 + 2 × 10 − 6 , where N ˜ [ * ] denotes the total number of glitches that the pulsar has experienced so far.

Accordingly, N ˜ Crab must be much smaller than N ˜ V e l a .

4) α c n → 1 as n → 0 : This is a consequence of the basic assumption that SBs are initially of nucleon-sizes and rotating uniformly with the ambient media.

5) ∏ i = 0 N ˜ ( Crab ) α c i = Ω c 0 / Ω c N ˜ = 7.7 and ∏ i = 0 N ˜ ( Vela ) α c i = Ω c 0 / Ω c N ˜ = 22

6) ∏ i = 0 n ≫ N ˜ α c i = S c ∞ / S c 0

7) When applied to the Crab and the Vela pulsars, the lapse of time between two successive glitch-events at the present time should be of order two years, i.e.

δ t n = t n + 1 − t n ≈ 2 years . (8)

Constraint number 1 (henceforth C1) means that pulsars are observed to spin down, i.e. d Ω / d t < 0 and therefore Ω ( t − ϵ ) > Ω ( t + ϵ ) , which is equivalent to Ω n > Ω n + 1 . C2 means that the element of the sequence should yield the ages of two pulsars correctly. Ages of objects are obtained by summing over the intervals between successive glitches up to the present time, e.g.:

τ t o t N ( Crab ) = ∑ n = 1 N ( Crab ) δ τ c n ≈ 10 3 yrs ,

and

τ t o t N ( Vela ) = ∑ n = 1 N ( Vela ) δ τ c n ≈ 10 4 yrs ,

where δ τ c n is calculated using Equation (11). C3 means that the present observations of the Crab and Vela pulsars reveal Δ Ω / Ω equal to 2 × 10 − 8 for the Crab and 2 × 10 − 6 for the Vela. Hence the elements of α c -sequence and the related parameters should be carefully constructed, so that the time intervals between two successive glitches at the present time should be equal to two years approximately and that their metamorphosis processes should be completed within the first hundred million years after their birth. C4 implies that at time = 0, i.e. when the pulsar was born, the size of the SB was comparable to those of nucleons in atomic nuclei. C5 means that the multiplications of the first N ( Crab ) elements of the sequence { α c n } should yield 7.7 in the Crab case, whereas the multiplications of the first N ( Vela ) elements should yield 22 in the Vela case. Here I rely on the works by [

Note that since α c n / α c n − 1 > 1 , then the corresponding sequence must not converge and therefore the above-mentioned 6th-constraint cannot be fulfill; specifically that the final cross section of the SB should be upper-bounded by that of the pulsar, i.e.

S n → n → ∞ S ⋆ ∩ S ∞ ≤ S ⋆ .

As it will be shown in the next section, this problem can be circumvented by

allowing the dipole magnetic field to asymptotically decay to zero at the end of the pulsar’s lifetime. In that case, the time between two successive glitch-events would stretch to infinity (see e.g. Equation (10) and Figures 5-7).

On the other hand, constructing a convergent sequence to describe the evolution of glitches in pulsars is not only physically inconsistent, but it is also contradictory to the observations of the Crab and Vela pulsars. One possible improvised sequence that may fulfill most of the above-mentioned constraints may read as follows:

α c n = ( 1 + ϵ 1 n ϵ 2 ) ϵ 3 n 1 + ˚ C 4 , for n : 0 → ∞ , (9)

where ϵ 1 = 10 − 12 , ϵ 2 = 5 × 10 − 2 , ϵ 3 = 1 + 10 − 10 and ϵ 4 = 2 × 10 − 2 .

Indeed, ϵ i , i = 1 → 4 have been carefully chosen and optimized so to agree with glitch-observation in pulsars, and in particular to those of the Crab and Vela pulsars.

Using the elements of { α c n } , a sequence for the cross sections of the SB can be constructed, { S c n } , from which the sequence { R c n } is then calculated to determines the radii of the permitted energy levels.

Assuming glitch-events in pulsars to evolve according to a well-defined “magic” sequence, then the here-presented α-sequence may be considered as one possible approximation, though additional modification may still be required.

The set of equations to be solved here reads as follows:

i . d ( I Ω ) n m / d t = − β B 2 Ω n m 3 + δ ( Ω n m − Ω c ) Δ I ˙ S B Ω S B ii . [ S Ω ] = 0 & S [ Ω ] = ( h / 2 m ) [ N V ] iii . M S B = 2 π ∫ 0 r c ( ρ c r + E ϕ ) r 2 d r iv . M n m = M ⋆ − M S B b v . M S B → M S c h as n → ∞ , (10)

I n m , Ω n m , B in Equation (i/10) denote the inertia, angular frequency and the magnetic field of the normal dissipative medium surrounding the SB. This equation describes the conservation of the total angular momentum of the ambient medium subject to the two forces:

・ The first term on the RHS of Equation (i/10) stands for the magnetic torque, which acts to brake the rotation of the ambient medium through emission of magnetic dipole radiation.

・ The second term on the RHS of Equation (i/10) represents the rotational energy gained by absorbing the vortices ejected by the SB during a glitch event.

The delta function on the RHS of Equation (i/10) is defined as follows:

δ ( Ω − Ω c ) = { 1 if Ω = Ω c 0 otherwise .

Equation (ii/10) describes the prompt reaction of the SB as clarified in the previous section. Here [ N V ] = δ N V ∝ Δ Ω S B represents the prompt loss of rotational energy from the SB-BL system during a glitch-event. Once the cross-section of the SB has increased from S n to S n + 1 , its total effective mass is calculated from Equation (iii/10), where I use E ϕ ≐ α 0 r 2 + β 0 to calculate the density of the dark energy. Here β 0 is the canonical energy scale characterizing the effective coupling of quarks (see [

The algorithmic structure of the numerical solution procedure runs as follows:

1) Construct the sequence { α c n } , which fulfills the list of conditions listed in Equation (8).

2) Calculate the elements of the sequences that corresponds to { Ω S B n } , { S S B n } , { R S B n } and { I S B n } .

3) Calculate the elements of the sequence that correspond to the total mass of the SB, i.e. { M S B n } t o t , including dark energy as well as its total inertia { I S B n } t o t .

4) Calculate the elements of the sequences that corresponds to the mass { M n m n } and inertia { I n m n } of the ambient dissipative medium.

5) Calculate the elements of the sequence { δ τ c n } using a physically reasonable formula that relates the magnetic field decay to the inertia { I a m n } , taking into account that both the magnetic field B and I n m should decay asymptotically go to zero as t → ∞ (see Equation (11)).

6) Compute the continuous decay of Ω n m of the ambient medium as dictated by the emission of magnetic dipole radiation (Equation i/10).

In solving the set of Equations (10), 1540/s is used as an initial value for Ω c ( t = 0 ) , which is in the range proposed by [

Note that within the time interval between two successive glitch events: [ t n , t n + 1 ] , both the inertia I S B and angular velocity Ω S B are constant, whilst Ω n m of the ambient medium decreases continuously as dictated by the loss of magnetic dipole radiation (see Ω n m in

δ τ c n = 1 2 β I B 2 { ( 1 Ω n + 1 c ) 2 − ( 1 Ω n − 1 c ) 2 } ≐ 1 2 β I B 2 Δ ¯ ¯ . (11)

Based on test calculations, no reasonable values for δ τ c n , that could agree with observations were found, unless I n m / B 2 are set to decrease with time.

Actually the numerical solution procedure is intrinsically sequential. Once the elements of the sequence { α c n } , for n = 0 → ∞ are fixed, the determination of the other sequences are straightforward. Since the rotational frequency and inertia of the SB are constant during the time period between two successive glitch-events, then the equation describing the evolution of torque interaction between the SB and the ambient medium, i.e. Equation (i/9), can be integrated straightforwardly.

Based on the numerical solution of the equations listed in Equation (10), I show in

The numerical solutions of the set of Equation (10) on the cosmological time-scale are shown in Figures 6-8. Here the long-term process of the metamorphosis is obvious: while the total energy of the SB is growing, the mass of the ambient medium is diminishing to finally disappear after roughly 5 × 10 6 yrs, or equivalently, after the object has suffered of roughly forty million glitches. At this time the object attains its maximum compactness, i.e. α S → 1 , it sinks deeply in spacetime as the metric coefficient, g r r , in the vicinity of the surface becomes extraordinary large. This implies that the redshift would be very large and therefore the object disappears from all possible observational windows. In the absence of external forces, thermal, magnetic and rotational energies of pulsars must decrease during their lifetimes. Assuming newly born pulsars to initially rotate with Ω n m ( t = 0 ) = 1540 / s and that their initial magnetic fields are B ( t = 0 ) = 10 13 Gauss, then the present calculations predict that these energies would be vanishingly small when the object becomes several million

years old (see

Finally, given that the initial mass and inertia of SBs are comparable to those of nucleons in atomic nuclei and therefore negligibly small compared to those of the ambient medium, then the time interval between two succussive glitches must be respectively short, though it becomes increasingly longer as the SBs age. Obviously, after one million years, the decreasing-rate of the inertia of the ambient medium surpasses that of the magnetic field, thereby shortening the time interval between successive glitches considerably (see Equation (10) and

In this article, a new scenario for explaining the origin and dynamics of the glitch-phenomena observed in pulsars and neutron stars has been presented. Accordingly, birth of pulsars most likely is associated with the formation of embryonic SBs, whose cores are made of supranuclear superconducting gluon-quark superfluid of uniform density ρ = 3 × ρ 0 . Following a well-defined discrete quantum scheme in combination with a scalar field, the mass and inertia of an SB should grow with time to finally metamorphose the entire pulsar into an invisible super baryon. The SB-pulsar system is strikingly similar to a terrestrial rotating ^{4}Helium inside a container. If both components are set to initially rotate with the same frequency and subsequently increasing/decreasing the frequency of the container in a continuous manner, then the enclosed helium-superfluid would increase/decrease its frequency discretely and promptly, as dictated by the Onsager-Feynman equation.

The analogy to the SB-pulsars system is obvious: The SB corresponds to ^{4}Helium, whereas ambient dissipative medium to the solid container.

However, the situation in the SB-pulsar is still much more difficult as it requires exact determination of the permitted quantum energy levels of a gravitationally bounded super-baryon, which is completely unknown.

In the here-proposed scenario, the nuclear shell model is adopted to describe these energy levels, which are assumed to coincide with the angular frequency levels. A glitch event here would correspond to a transition from one energy level to the next, or equivalently, from one spin frequency to the next lower one. These transitions are provoked mainly by the two mechanisms:

1) Loss of rotational energy of the ambient medium through the magnetic torque interaction of pulsar with their surrounding;

2) Enhancement of the effective mass-energy of the SB by the scalar field. Once the entire pulsar has metamorphosed into an SB, the mass of the object would double, provided the pulsar has the compactness parameter α s = 1 / 2 initially. Hence scalar fields act as compressors to the nuclear matter of the ambient media and enforce their transition into incompressible superconducting gluon-quark superfluids.

In fact, based on the discussion of Section (2), the solution procedure here can be used to estimate the relative width of the BL as follows:

δ R B L n R n = Ω n Ω n + 1 − 1 ≈ Δ Ω Ω ≈ { 10 − 8 Crab 10 − 6 Vela (12)

where R n is the radius of the SB at time t n . Obviously, as Ω 0 / Ω 1 ≈ 1 , but much smaller than Ω N / Ω N + 1 for N ≫ 1 and that δ R B L ( n = 0 ) ≈ 1 fm, then δ R B L n is practically a measure for the length scale over which the residual of the strong nuclear force is communicated to the surrounding matter. As it is shown in

δ N ≥ 2 m h S S B ( Δ Ω ) = 5.32 × 10 − 4 S 0 , (13)

which is equivalent to require S 0 ≥ 1.88 × 10 3 cm 2 .

Unlike classical models of glitches, the expelled vortices here are subject to destructive frictional effects exerted by the surrounding dissipative medium and therefore they may diffuse out long before affecting the crust. It should be noted that even in the ideal case, vortex lines may interact with each other, reconnect and turn turbulent (see [

E V k i n ≤ ( h 4 m I Ω S ⋆ ) ≈ 10 31 erg . (14)

Thus the minimum number of vortices needed to be expelled by the SB in order to be observationally noticeable would be N ≥ I Ω Ω ˙ / E V k i n ≈ 10 7 . In the case of the Vela pulsar this implies that the radius of the enclosed SB must be larger than 10 − 2 R ⋆ , i.e. of a macroscopic length scale.

In the Crab case, the radius of the SB compared to that in the Vela can be inferred from the time-interval between two successive glitch-events as described in Equation (11). Inserting the present values of Ω , Ω ˙ of both pulsars, then we obtain:

δ τ Crab ≈ 1 2 β I B 2 ( 10 − 7 ) and δ τ Vela ≈ 1 2 β I B 2 ( 10 − 5 ) . (15)

Recalling that the magnetic fields and the time-interval between two successive glitch-events in both system are of the same order, then we conclude that the inertia of the ambient dissipative medium in both pulsars may be related to each other as follows:

I V e l a ≈ 10 − 2 I C r a b . (16)

This implies that the radius of the SB inside the Vela pulsar is at least 5/2 larger than the one inside the Crab, which is a reasonable outcome, when taking into account the difference in mass and age of both pulsars.

It should be noted however that if the size of the embryonic SB is relatively small, then the number of ejected vortices would be respectively small and therefore may diffuse out quickly without causing an observationally significant spin-up of the crust.

Finally, while much more efforts are needed still to refine and further constrain the universal “magic” sequence describing the evolution of glitch-events in pulsars and neutron stars, I do think that the present scenario is promising as it is capable of comfortably accommodating observational data to provide answers related to the quantum origin and dynamics of glitches in pulsars. Specifically, glitch-observations appear to support the formation and evolution of embryonic SBs at the centers of pulsars, whose mass-growth follows a well-defined quantum sequence. The existence of such embryonic objects are conditioned to the occurrence of a phase transition form dissipative nuclear fluid into incompressible superconducting gluon-quark superfluid, which in turn requires a universal scalar field ϕ as a catalyst.

As quarks make at most 2% of baryon mass, we may simply view the ϕ -field as a universal mechanism for enhancing the density of communication channels between quarks inside SBs, which in turn enhances its effective mass.

Finally, similar to quarks in nucleons of atomic nuclei, it was argues that stellar-size SBs cannot exist in free-space [

Hujeirat, A.A. (2018) Glitches: The Exact Quantum Signatures of Pulsars Metamorphosis. Journal of Modern Physics, 9, 554-572. https://doi.org/10.4236/jmp.2018.94038