Glitching Pulsars: Unraveling the Interactions of General Relativistic and Quantum Fields in the Strong Field Regimes

In this article we modify our previous model for the mechanisms underlying the glitch phenomena in pulsars. Accordingly, pulsars are born with embryonic cores that are made of purely incompressible superconducting gluon-quark superfluid (henceforth SuSu-cores). As the ambient medium cools and spins down due to emission of magnetic dipole radiation, the mass and size of SuSu-cores must grow discretely with time, in accordance with the Onsager-Feynmann analysis of superfluidity. Here we argue that the spacetime embedding glitching pulsars is dynamical and of bimetric nature: inside SuSu-cores the spacetime must be flat, whereas the surrounding region, where the matter is compressible and dissipative, the spacetime is Schwarzschild. It is argued here that the topological change of spacetime is derived by the strong nuclear force, whose operating length scales are found to increase with time to reach ( ) 1  cm at the end of the luminous lifetimes of pulsars. The here-presented model is in line with the recent radio- and gravitational wave observations of pulsars and merger of neutron stars.

α , is generally larger than half (see [1] [2] [3] and the references therein). Recalling that the average density of matter in the interiors of UCSs is larger than the nuclear one, 0 ρ , then these objects are well-suited for probing the interaction of general relativity with quantum fields in the strong field regime (Figure 1).
In fact pulsars and NSs have been successfully used for probing general relativity, such as the emission of gravitational waves (GWs) in binary pulsars and direct detection GWs through merger of neutrons stars and black holes (see [4] and [5] for detailed reviews). However, the internal structure of UCSs will continue to be a mystery and remain a challenge for astrophysicists not for a short time (see the review in [6] and the references therein). On the other hand, the glitch-phenomena observed to associate pulsars [7] [8] are events that may be considered to be connected to the dynamics of matter inside their interiors (see [6] [7] [9] [10] and the references therein). However, based on observations and theoretical arguments, one may try to constrain the nature of matter inside normal UCSs as follows: 1) The observed continuous spin-down of UCSs results from the combined loss of magnetic and rotational energies (see [6] and the references therein).
These are estimated to fall in the range: [10 37 , 10 40 ] erg/s, depending on the nature and age of the concerned UCSs and specifically whether these concerned ones are magnetars, pulsars or neutron stars [9] [11]. This would imply complete exhaustion of the stored removable energies 1 inside pulsars within about ten million years, provided that the heat conductivity operates on length scales that are larger than the nuclear ones. Consequently, very old and isolated UCSs, including those formed through the collapse of the first generation of stars, must be dark by now and therefore excellent candidates for stellar black hole. This argument would explain why neither UCSs nor black holes have never been observed in the mass range 2 5 2) The rest of thermal, kinetic and magnetic energies that are left from the collapse of the progenitors of UCSs are transported outwards into the outer shells and subsequently liberated away. Due to the relatively small thermal energy content in UCSs and in the absence of nuclear energy generation, the Tolman-Oppenheimer-Volkoff equation (TOV-equation) may still accept a positive gradient of the thermal energy, turning the central core to the coldest region inside UCSs, thereby facilitating a phase transition of the compressible dissipative nuclear matter into incompressible superfluid (see [12] for further details on the properties of incompressible fluid flows). Such a phase transition may still occur even when the matter's temperature is still beyond several million degrees [11] (Figure 2). 1 i.e. energies other than rest energy.  = , it is energetically favourable for neutrons to start deconfine the quarks and merge together to form an ocean of gluon-quark superfluid, though the resulting global gluon cloud must be much more energetic.
3) Theoretical studies show that in the regime of nuclear density and beyond, almost all EOSs used for modelling the state of matter in the interiors of UCSs tend to converge to the limiting case: P ε = (see [9] and the references therein).
However, this state corresponds to pure incompressible nuclear fluids [2]. Due to causality and stability reasons, incompressible nuclear fluids should be free of all other removable energies, such as thermal, kinetic and magnetic energies.
Hence the core is expected to end up as a condensate with zero-entropy, whose constituents communicate with each other with the speed of light and therefore the momentum exchange between them must saturate around a universal maximum. In this case, the coupling constant in the context of asymptotic freedom in QCD would converge to its universal minimum value, where quarks inside merged-neutrons are allowed to oscillate locally [13] [14]. Such a zero-entropy superfluid 2 cannot accept spatial variations or stratification due to gravity and 2 Vanishing energy loss due to acceleration/deceleration should be excluded.  This state may be maintainable, if involved quarks are attached to a spatially fixed lattice, where they   may behave collectively as a single macroscopic quantum entity. Indeed, although both entropy and energy states of colliding protons at the RHIC and LHC differ significantly from those at the centers of UCSs, these experiments have shown that the properties of the resulting fluids mimic those of perfect fluids (see [15] [16] for further details). 4) Based on astronomical observations, there appears to be a gab in the mass spectrum of compact objects: neither black holes nor NSs have never been observed in the mass range 2 5 Moreover, intensive astronomical observations of the newly detected NS-merger GW170817 [1] failed to unambiguously classify the nature of the remnant object, though theoretically it must be a stellar black hole. These both facts are in line with our scenario which predicts that massive pulsars should evolve on the cosmological time scale toward forming massive dark energy objects that become indistinguishable from stellar black holes. A schematic description of the scenario is shown in Figure 1

The Model and the Governing Equations
For modelling the internal structure of pulsars, the general relativistic field equations should be solved: where , G T µν µν are the Einstein and the Stress-Energy tensors, respectively, and κ is a coefficient (see [17] [18] for further details). The rotational, magnetic, thermal and other energies in UCSs are generally about three orders of magnitude smaller than the rest energy and therefore they may be safely neglected. In this case, Schwarzschild spacetime is most suited for modelling their internal structures. Hence the corresponding metric reads: , , , , s r r G P m r correspond to the radius, Schwarzschild radius, gravitational constant, pressure and the dynamical mass, respectively (see Sec. (4) in [2] for further details).
Recalling that GR is incapable of modelling gravitationally bound incompressible matter, then the spacetime embedding pulsars may be decomposed into two separate domains: a flat spacetime that embeds the SuSu-core and a surrounding Schwartzschild spacetime that embeds the ambient media as mentioned above. The decomposition of the domain is motivated by relativistic cau-Journal of Modern Physics sality which prohibits fluid stratification or spatial variation of the density of purely incompressible fluids; hence why the spacetime embedding the core must have zero-curvature, i.e. purely flat.
Based thereupon the line element in the core is described by the metric in spherical polar coordinates ( ) , , , t r θ ϕ : The physical properties of the SuSu-core are set to affect the structure of the ambient medium by allowing On the other hand, as the matter in the rotating core is in incompressible superfluid state, it must contain a discrete array of vortices. Their total number, N, the mass and size growth of the core are determined through the Onsager-Feynmann equation for modelling quantized circulations: correspond to rotational velocity, line-element along the circular path, the reduced Planck constant, the mass of the superfluid particles and the enclosed number of vortices, respectively [19] [20] Inside the core, the quantum fluid is assumed to have zero-entropy and has reached the critical supranuclear density, where bag E ∆ is the bag energy enhancement needed to confine the quarks inside the super-baryon and N denotes the total number of merged neutrons to forming the super-baryon, n E is the rest energy of a single neutron. In this case the energy and pressure of the incompressible gluon-quark superfluid inside the core read: ε is the energy density of the baryonic matter and are the energy density and pressure of the scalar field, which are assumed to be identical to bag E ∆ . As the matter inside the cores of pulsars is set to be incompressible and stationary, it is reasonable to assume constant φ = . In this case, φ  and φ ∇ must vanish and ( ) V φ can be considered as the energy density required for deconfining the quarks inside the super-baryon, i.e., the SuSu-core. Here we set Thus the total local pressure of incompressible gluon-quark superfluid inside SuSu-cores vanishes completely. Exterior to the core, the matter is said to be baryonic, dissipative, compressible and embedded in a Schwarzschild spacetime.

Initial and Boundary Conditions
The present scenario is based on a previous model for the origin of glitches in pulsars. Accordingly, pulsars and young neutron stars are expected to undergo billions of glitch events during their luminous lifetimes before ending as ultra-compact and invisible dark energy objects.
In the present study, we select several epochs in the lifetimes of pulsars and the corresponding internal structures are calculated subject to the following conditions:  The pulsar has the initial mass of 1.33  , which is made of a purely baryonic compressible matter.
The total density inside SuSu-cores amounts to: The solutions here are based on adapting the parameters  and γ of the EOS in such a manner so that the initial baryonic mass of 1.33  the pulsar remain preserved.
3) The ultimate final phase: Here the radius and mass of the SuSu-core are roughly equal to the critical values: Here the initial pulsar must have metamorphosed entirely into an invisible Su-Su-object.

Results
For enhancing the spatial accuracy of the calculations, an explicit adaptive mesh refinement (EAMR) has been developed, in which the aspect ratio, max min d d where the pressure vanishes is vitally important and therefore the resolution should be maximally refined ( Figure 3). Unlike AMR, solution methods based on EAMR generally converge much faster than their AMR-counterparts, hence contributing to the global efficiency and robustness of the numerical solution procedure (see [21] and the references therein).
Based on the previous study [22], five episodes in the cosmological evolution of pulsars have been selected (see Figure 4). For each episode, the TOV equation, modified to include dark energy input at the background of the here-presented bimetric scenario of spacetime, is solved. Figure 3. The distribution of the finite volume cells as function of radius. Here the explicit adaptive mesh refinement (EAMR) has been employed to increase the accuracy at both the interface between the core and the ambient medium as well as to outer radius of the object. Note that as the SuSu-core becomes more massive, the curvature of spacetime embedding the ambient medium is enhanced, which, in turn, compresses the ambient medium even more, thereby reducing the effective radius of the pulsar.
This would explain the mass-radius relations displayed in Figure 7 and Figure 8.
Here the corresponding Schwarzschild radius increases and propagates outwards to finally meet the decreasing radius of the contacting pulsar. The overlapping of both radii is expected to occur at the end of the pulsar's luminous lifetime (see Figure 4 and Profile 5 of Figure 5).
As purely incompressible superfluids in flat spacetime have zero-spatial variations ( Figure 9), then the gravitational potential inside SuSu-cores should attain a sequence of constant values. Their magnitudes correlate with the mass and size of the SuSu-core. Consequently, as the SuSu-core becomes more massive, the gravitational redshift of the pulsar increases to finally reach very large values at the end of the pulsar's luminous lifetime (Profile "6" in Figure 10). Here the Schwarzschild radius becomes almost equal to the effective radius of the object, enforcing the object to sink deeply in spacetime to become invisible.
Note that the profiles shown in Figure 9 and Figure 10 are not just schematic representations, but have been obtained using direct numerical computations. Figure 5. The distribution of the baryonic pressure inside an evolving pulsar after five selected glitch events at different epochs. Inside the core, in addition to the rest energy of the baryonic matter, there is an energy enhancement due to the scalar field, which is equivalent to the additional energy required for re-confining the sea of quarks. The transition between both regions is not smooth as the matter inside the core evolves quantum mechanically, whereas the ambient medium obeys the normal laws of continuum.  The total enclosed mass of the object versus radius after selected glitch events. As the mass of the SuSu-core increases, its compactness of the whole object increases as well to finally reach the critical Schwarzschild mass. Figure 8. The mass-radius relation of an evolving pulsar after selected glitch events. As the core becomes more massive, the Schwarzschild radius grows and converges asymptotically to the effective radius of the entire contracting object.
A. A. Hujeirat, R. Samtaney Figure 9. The radial distribution of the gravitational potential shown for different epochs in the lifetime of the pulsar, which is equivalent to the radial projection of spacetime inside the SuSu-core governed by Minkowski spacetime, as well as in the surrounding region governed by the Schwarzschild metric. Obviously, as the object become more massive, it sinks deeper in spacetime, becomes gravitational redshifted to finally ends as an invisible dark energy object. Figure 10. The radial distribution of the gravitational redshift (Z) of the pulsar at different evolutionary epochs. Profile "1" corresponds to the low "Z" immediately after the pulsar was born, whereas Profile "6" depicts the divergent limit of "Z" at the surface of the object when it becomes completely invisible.

Summary and Discussion
The here-presented pulsar model is analogous to the terrestrial rotating superfluid Helium in a container. While the SuSu-matter inside the core corresponds to rotating superfluid Helium, the ambient dissipative medium would represent the solid container.
However, the pulsar model here is much more sophisticated, as quantum and general relativistic effects are taken into account. In the following, we discuss the main features of the model. • The formation of incompressible gluon-quark superfluids is associated with energy enhancement of the gluon field at the surface of the core, which is needed for globally re-confining the enclosed quarks, thereby effectively enhancing the effective mass of the core. • The phase transition of the quantum fluid in the BL from compressible-dissipative into incompressible superfluid states is associated with changes of the spacetime topology from a curved spacetime into a purely flat one. The latter reaction is due to causality, which prohibits incompressible superfluids to be embedded by a curved spacetime. Indeed, the topology change of spacetime is associated with the emission of gravitational waves.
However, due to the tiny little volume of the BL, the detection of GW-emission during glitch-events would be much below the sensitivity of today GW-detectors.
• Both the phase transition of matter and the topological change of spacetime in the boundary layer can be provoked by the strong nuclear force (SNF) only, which appears to develop into a macroscopic force as the pulsars ages.
In fact, as the core is made of SuSu-matter, the core is equivalent to a super-baryon, in which the enclosed quarks are shielded by a gluon cloud [23]. As it will be shown later, the SNF transmitted by this cloud is found to correlate nicely with the size and mass of the core. The neutrons in the BL become increasingly locked to the core, thereby gradually adapting the physical conditions of the matter inside the core.
Following the scenario of [22], the occurrence of glitches in pulsars follows a well defined sequence, whose elements are { } R Ω be respectively the radius and angular frequency of the SuSu-core at the verge of the glitch event number "n". The rotational energies of the SuSu-core at time n t and 1 n t + read: 5  , BL n n R δ δΩ here denote the width of the boundary layer between the Su-Su-core and the ambient dissipative medium and the absolute difference between the rotational frequencies of the core before and after the glitch event, respectively. Based thereon, the width of the BL during different evolutionary epochs of the pulsar may be estimated as follows: where age τ stands for the age of the pulsar, provided the object is perfectly isolated and shielded from whatsoever external effects. In calculating the ratio at 0 age τ = , the initial values  Figure 8 of [22]. For determining n SB R we use the current glitch-observation of the Vela and Crab and try to extrapolate them to other epochs. Based on Equations (13,14,15) in [22], the correlation of the inertia of the ambient dissipative media, AM I , of the Crab and Vela pulsars reads: 2

10
, whereas the requirement that the ejected rotational energy off the core of the Vela pulsar in the form of vortices to be observationally noticeable, implies that  An urgent question here is: What is the nature of force in the boundary layer that is capable of triggering glitches?
Here we argue that the SNF is the deriving force that is capable of changing instantly both the physical properties of the quantum fluid as well as the topology of spacetime in the BL.
-Following [22], the core mimics a super-baryon (SB), in which the enclosed quarks are shielded by a gluon-cloud. The medium inside the SB interacts with overlaying supranuclear dense neutrons via vector mesons: the messengers of the SNF. In turn, the SNF locks the neutrons gradually to the core and enforces them to rotate rigidly with the core.
-The SNF generally operates effectively on nuclear length scales, i.e. when ( ) However, the inertia of a core having the neutron-size would be approximately 100 orders of magnitude smaller than that of the ambient medium. In this case, its rotational and thermal decoupling from the ambient medium would be almost impossible. On the other hand, in order for a SuSu-core to survive and affect the dynamics of the ambient medium, its inertia must have a minimum critical value. Here we recall that the initial frequency of a newly born pulsar is approximately 1400/s [24]. In order for a SuSu-core to be able to eject one single vortex, its radius must be larger than 10 3 cm, hence its minimum mass: correlates nicely with the size of the enclosed SuSu-core and it may even become of macroscopic size at the end of pulsars luminous lifetime.
-The next related question reads: What keeps the core dynamically stable against the weight of the overlying massive shells, or equivalently to be stable against further compression by the surrounding curved spacetime? As the gravitational potential inside the core is constant, then the medium has zero stratification due to gravity. In this case the only opposing force against gravitational collapse is the pressure due to the uncertainty principle as well as the pressure jump across the interface between SuSu-cores and the ambient compressible dissipative medium. The medium in the BL is still fermionic and the governing pressure is due to the Pauli exclusion principle. Moreover, similar to the bag energy confining the quarks inside individual baryons, the dynamics across the surface of the SuSu-core may be dominated by an extraordinary strong surface tension that is capable of confining the ocean of quarks quarks as well as opposing further compression.
-As the neutron's density in the BL approaches the critical value, crit ρ , and as the chemical potential becoming closer to the upper-limit, cr µ , the EOS should converge to 2 0 P a n ε → = , where n is the number density.
This limiting EOS corresponds to purely incompressible fluids and therefore the topology of the embedding spacetime ought to change into purely flat spacetime. In this case, two instantaneous reactions in the BL are expected to occur (see Figure 11 detailed description): Figure 11. The three-stage glitch scenario. In (a) the time-development of both rotational frequencies of the SuSu-core (blue-color) and of the ambient medium (red-color) are shown, using GR-numerical calculations. In (b) the rigid-body rotation and increasing size of the SuSu-core for four successive glitch-events are shown. The boundary layer here has the width R δ , area S δ and contains N δ vortices. In the lower panel we show the three-stage scenario of the glitch phenomena in pulsars. In (c1) the SuSu-core and boundary at the verge of a glitch event are shown. Once the SNF has locked the neutrons inside the BL, to the core and enforced them to adopt the same thermo-and hydrodynamical properties of the matter inside the core, then the spacetime embedding the BL undergoes a topological change into a flat spacetime and merges with that of the core (c2). However, once the resulting difference ∆Ω between n SB Ω of the core and the ambient medium n AM Ω has surpassed a critical value, then the core must undergo a transition into the next lower energy state by expelling a certain number of vortices (c3). In turn, the ambient medium absorbs the vortices and re-distribute them viscously, thereby giving rise to the prompt spin-up observed in gliching pulsars.
1) As the supranuclear dense matter in the BL becomes increasingly colder and its thermodynamical state and dynamics become indifferent from those of the matter inside the core, then it merges instantly with the core, thereby increasing its size and effective mass.
2) However, when where we have optimistically assumed that the viscous length scale is of order  [25].
• Based on the here-presented scenario, the mechanisms underlying the glitch phenomena in pulsars are due to dramatic changes of the physical conditions of the matter as well as of the topology of spacetime in the BL. Here the system operates on two different time scales: a relatively long times scale, δ Ω = Ω − Ω becomes increasingly larger. The above process provokes an instantaneous reaction of the SuSu-core: similar to superfluid Helium in a container, once n δ Ω has surpassed the critical value n cr δ Ω , then the SuS-core reacts promptly: the BL then merges with the SuS-core, thereby triggering the core to make a transition into the next lower and quantum-mechanically permitted energy state by expelling a certain number of vortices. quantum τ may be predicted to be order • By the end of luminous lifetime of pulsars, the discrete repetition of mergers of the BL with SuSu-core should have metamorphosed the entire pulsars into a SuSu-object, thereby doubling its initial mass and reducing its radius to roughly coincide with the corresponding event horizon. For remote observers the object becomes practically invisible and therefore indistinguishable from stellar black holes. The present evolutionary scenario may explain nicely, why neither black holes nor neutron stars have ever been observed in the mass range as well as why the NS-merger GW170817 [4] did not collapse to form a stellar BH. If nature indeed forbids the formation of BHs in this mass regime, then our understanding of the big bang, the nature of dark matter and dark energy in cosmology should be revised (see [3] [26] for further details).