g Ω } n must be identical to the observed values of the Crab and to the Vela pulsars, respectively.

Table 1. A list of the main observational data of the Crab and Vela pulsars relevant for the present study (see [2] - [7] and the references therein).

Figure 1. A schematic description of the bimetric spacetime inside glitching pulsars: The incompressible superfluid core is embedded in a Minkowski spacetime whereas the ambient media are enclosed in Schwarzschild spacetime.

l { Δ Ω g Ω } n n α < , i.e. the sequence must converge to a finite value. Moreover, our test calculations have shown that for t t V e l a , t ( { Δ Ω g Ω } n ) < 0 as otherwise the magnetic field would fail to spin-down the crust and therefore to surpass Δ Ω c r required for triggering a prompt spin-down of the core into the next lower energy state.

Indeed, one possible sequence which fulfills the above-mentioned constraints, though it might not be unique, is shown in Figure 2.

2) The initial conditions used here are Ω 0 = Ω ( t = 0 ) = 1440 Hz (see [10] and the references therein). Here both the core and the ambient medium are set to initially rotate with the same frequency, i.e. Ω c n = Ω a m n ( t = 0 ) = Ω 0 . The initial total mass of the pulsars and that of the core are taken to be: M c = ϵ 0 M and M 0 = M ( t = 0 ) = ϵ 1 M , respectively, where ϵ 0 and ϵ 1 are parameters whose values are determined through a global iteration procedure. The initial magnetic field strength is taken to be B ( t = 0 ) = 10 13 Gauss.

Figure 2. The elements of the quantum sequence Δ Ω g n / Ω are shown as function of time in year units. Both the core and the ambient medium are set to rotate with the same frequency 1440s−1 initially. As the pulsar cools down, Δ Ω g n / Ω starts increasing to reach 4 × 10−9 after approximately 1000 yr, (which corresponds to the Crab phase/blue-star) and 8.15 × 10−7 after 11,000 yr (which corresponds to Vela phase/red-star). Δ Ω g n / Ω here is measured in units of Ω = 200 / s .

3) The elements of the sequence { Ω a m } n are obtained through the energy balance equation:

d d t ( 1 Ω a m 2 ) = α E M B 2 I a m , (1)

where I a m is the inertia of the ambient compressible dissipative medium, which, due to the increase of the SuSu-core, must decrease on the cosmic time and α E M = 4.9 × 10 4 is a non-dimensional constant.

Thus, the time-evolution of the core’s rotational frequency proceeds as follows: for a given Ω c n , the ambient medium is set to decrease its frequency continuously with time through the emission magnetic dipole radiation. This implies that the difference Δ Ω c a m should increase with time until Δ Ω c a m has surpassed the critical value Δ Ω c r , In this case three events are expected to occur promptly:

l The rigid-body rotating core changes abruptly its rotational state from E c n , r o t = 1 2 I c n ( Ω c n ) 2 into the next the quantum-mechanically permitted lower energy state: E c n + 1, r o t = 1 2 I c n + 1 ( Ω c n + 1 ) 2 . This process is associated with ejection of a certain number of vortices into the boundary layer (BL) between the core and the overlying dissipative medium.

l The ejected vortices by the core are then absorbed by the differentially rotating dissipative medium and re-distributed viscously. Hence the medium in the BL would experience the prompt spin-up: Ω a m n Ω a m n + Δ Ω g n , where Δ Ω g n is deduced from the sequence { Δ Ω g Ω } n .

l The radius of the core is set to increases as dictated by the Onsager-feynmann equation:

V d l = 2 π m N , (2)

where V , l , , m , N denote the velocity vector, the vector of line-element, the reduced Planck constant, mass of the superfluid particle pair and the number of vortices, respectively (see [11] for further details). Imposing zero-torque condition on the incompressible SuSu-core, i.e., d d t ( Ω I ) c = 0 . We then obtain the following recursive relation:

( Ω S ) c n + 1 = ( Ω S ) c n S n + 1 = ( Ω c n Ω c n + 1 ) S n (3)

where S c n π ( R c n ) 2 and R c n correspond to the cross-sectional area of the SuSu-core and to the corresponding radius, respectively. The increase in the dimension of the core implies that the matter in the geometrically thin boundary layer between the SuSu-core and the ambient medium should undergo a crossover phase transition into an incompressible superfluid, whose total energy density saturates around the critical value ρ c r 6 ρ 0 (see [12] and the references therin). The growth of the core proceeds on the cosmic time scale and ends when the pulsar has metamorphosed entirely into a maximally compact invisible dark energy object and therefore becomes observationally indistinguishable from a stellar black hole.

2. Solution Procedure & Results

The set of equations consists of the TOV equation for modeling the compressible dissipative matter in the shell overlaying the incompressible gluon-quark superfluid core, whereas the latter is set to obey the zero-torque condition and to dynamically evolve according to the Onsager-Feymann equation (for further details see Sec. 2 and Eq. 10 in [8]).

The global iteration loop is designed here to find the optimal values of the parameters: α 0 , α 1 , the elements of the sequence Δ Ω c r n and the decay rate of the magnetic field. These values should fulfill the initial and final conditions, the currently observed values of the time passages between two successive glitch-events Δ t g both of the Crab and the Vela pulsars, the current observed values of their magnetic fields masses.

Indeed, our intensive computations reveal that optimal fitting may be achieved for M c ( t = 0 ) 0.029 M , a sequence of Δ Ω g n , whose elements are shown in Figure 2. In Figure 3, Figure 4 the optimal values of Δ Ω c r n and Δ Ω g n are shown versus cosmic time, whereas in Figure 5 we show the time-development of the rotational frequencies of the core, the ambient medium and of the magnetic field during the first 10 to 100 years after the birth of the pulsar. The long-term evolution of the magnetic field and the growing mass of the core and of the entire object are shown in Figure 6 and Figure 7. Here the mass of the pulsar’s core grows with time to reach 0.15 M after 1000 years and reaches 0.55 M after 11,000 years; hence reproducing the exact total masses of 1.4 M for the Crab and 1.8 M for the Vela pulsars as revealed from observations. The relative ratio of inertia of both cores reads I C r a b / I V e l a 3 / 20 .

Due to the incompressible, superfluid and supreconducting character of the core, the evolution of the magnetic field is solely connected to the dynamics of the ambient compressible and dissipative matter in the shell as well as to its dimensions (see [13] for further details on the physical aspects of compressibility of fluid flows). As the mass and dimension of the core grow with time, the surrounding shell must shrink. In this case, conservation of the magnetic flux should strengthen the magnetic field intensity. This interplay between the loss of magnetic energy due to loss of rotational energy and enhancement by conservation of magnetic flux in combination with dynamo action and other mechanisms,

Figure 3. The elements of the sequence Δ Ω c r n versus cosmic time. Each element corresponds to the critical difference between the rotational of frequency of the core and that of the ambient medium, beyond which the core undergoes a prompt spin-down to the next lower energy state.

Figure 4. The elements of the sequence Δ t g versus cosmic time. Each element corresponds to the time passage between two successive glitch events. The actual values that correspond to the Crab and Vela pulsars are shown in blue and red stars. These time passages increase dramatically with time due to the decay of both magnetic field intensity and rotational energy.

Figure 5. The elements of the sequences Ω c n , Ω a m n and B a m n that corresponds to the rotational frequencies of the core, the ambient medium and of the magnetic field during the first 10 years. During a passage of time between two successive glitches, the core rotates rigidly with a constant frequency Ω c n = Ω c o r e n (dotted line), whereas the ambient medium spin-down in a continuous manner (solid line). During the glitch, the core spin-down abruptly, triggering a prompt spin-up of the ambient medium in the boundary layer between the rigid-body rotating core and the overlying differentially rotating medium. The enhanced spin-up of the ambient medium in combination with the decreasing volume enclosing this matter gives rise to magnetic field B a m which evolves in a similar discrete manner (dash-dot).

Figure 6. The cosmic evolution of the magnetic field-B, of a newly born pulsar in units of 1013 G. The superimposed blue and red stars correspond to the current B of the Crab and Vela pulsars.

Figure 7. The mass growth of a newly born pulsar having initially M 0 = 1.25 M and an embryonic SuSu-core of 0.029 M . After approximately 1000 yr the pulsar recover the mass of the Crab ( M 0 = 1.4 M ) and 10,000 yr later M 0 = 1.8 M that corresponds to the mass of the Vela. At the end of the luminous life time, which lasts for approximately 10 million years, the pulsar enters the dark phase with a total mass of 2.5 M , which corresponds to a maximally compact invisible dark energy object. The plotted mass here is in units of M 0 .

may clarify the very weak decay of magnetic field as pulsars evolve from the Crab to the Vela phase. Mathematically, let the magnetic energy in a shell of a newly born pulsar be:

E M = B 2 8 π d v ~ B 2 6 ( R 3 R c 3 ) , (4)

where R denotes the pulsar’s radius. Assuming E M to roughly decay as the rotational energy E Ω , then we obtain:

B = α Ω Ω a m 2 M a m 1 / 2 R 2 R c 2 R 3 R c 3 , (5)

where α Ω is constant coefficient.

On the other hand, dynamo action in combination with magnetic flux conservation and other enhancement mechanisms would contribute positively to the magnetic field, that can, for simplicity absorbed in the term: B + = α B / ( R 2 R c 2 ) . The coefficient α B is set to ensure that the magnetic field remains in the very sub-equipartition regime. Hence the interplay between magnetic loss and enhancement would yield an effective magnetic field that evolves according to:

B t o t = α B R 2 R c 2 α Ω Ω a m 2 M a m 1 / 2 R 2 R c 2 R 3 R c 3 . (6)

Consequently, our model predicts that the decreasing volume of the shell enclosing the ambient medium in combination with dynamo action in the boundary layer could potentially be the mechanism that keeps the decay of magnetic fields in pulsars extremely weak.

In fact our model predicts the glitch activity of a newly born pulsar, which evolves into a Crab phase, followed by a Vela phase and finally by an invisible phase, to be approximately two orders of magnitude larger than it was estimated by other models (see Figure 8 to be compared to [1] [5]). According to our

Figure 8. The glitch activity of a newly born pulsar versus cosmic time. In the very early times, the pulsar underwent millions of glitches, though the total ejected rotational energy was relatively very low. These activities start to be significant as the pulsar ages and become maximally effective between 100 and 60,000 years, followed by a decreasing phase, during which time-passages between successive glitches become increasingly longer.

model pulsar may undergo millions or up to billions of glitches during their luminous life time with passages of time between two successive glitch events that range from nanoseconds in the very early time up to hundreds or even thousands of years toward the end of their luminous life times (see Figure 4). The vast difference in the evolution of glitch activity between the two approaches here may be attributed to the strong non-uniformity of time-duration between glitch events.

Moreover, the model also predicts the occurrence of under- and overshootings that have been observed to accompany the glitch events in the Vela pulsar (see [14] and the references therein). In the case of the Vela, when the core expels certain number of vortices and moves to the next lower energy state, the enhanced rotational energy of the matter in the BL amplifies the magnetic field strength. Due to the non-locality of magnetic fields1, this enhancement is communicated to the crust via Alfven waves, V A , whereas the excess of rotational energy is communicated via shear viscosity with an effective propagational velocity V v i s . As these two speeds are generally different with V A > V v i s in most cases, the time-delay in the arrival of communication enforces the crust to react differently. Specifically, the arrival of magnetic enhancement prior to the rotational one leaves the crust subject to an enhanced magnetic braking and therefore to a stronger reduction of its rotational frequency (see the top panel of Figure 9).

Indeed, in the case of Vela, the propagational speed of Alfven waves may be esitmated to be of order V A ~ B / ρ 10 8 cm / s . Hence the enhanced MFs in the BL would be communicate to the crust within δ τ M F = Δ R / V A = ( R R c ) / V A 10 2 s . On the other hand, supplying the crust with rotational energy would proceed on the viscous time scale, which is estimated to be: δ τ v i s = ( Δ R ) 2 / ν v i s (see [15] and the references therein). Under normal astrophysical conditions we may safely assume that V v i s / V A = α 2 1 and that the length scale, l v i s , over which viscous interaction occurs, l v i s is much smaller than the width of the shell Δ R , or equevalently l v i s = α 2 Δ R , with α 2 = 1 . Therefore δ τ V i s = Δ R / ( α 1 α 2 V A ) 10 2 / ( α 1 α 2 ) 1 s , where we reasonably set α 1 = α 2 = 0.1 . Consequently, the observed undershooting most likely results due to the time delay of the arrival of communications via magnetic fields and viscous torque, which amounts to δ τ v i s / δ τ M F = 1 / ( α 1 α 2 ) 100 . On the other hand, the observed overshooting can be attributed to the case in which the viscous front transporting rotational energy from the BL outwards has reached the crust. As δ τ M F δ τ v i s , at the end of δ τ v i s , the magnetic field intensity in the BL should have returned to values comparable or even lower than prior to the glitch event.

*I.e. In the absence of magnetic monopoles.

Moreover, the observed order in which undershooting followed by overshooting is an indication for a time-delay in the arrival of communication resulting from V A > V v i s and from the significant difference of the locations of the BL and the crust. This order is expected to reverse if V A < V v i s .

Figure 9. Under and overshooting of Ω a m shortly before and immediately after a glitch event (top Figure). At a given instant of time, t, the rotational frequency of the matter in the BL, Ω a m , differs from that of the crust Ω c r u s t . Due to their different locations (see lower panel), the response of the crust to the dynamical changes of matter in the BL depends strongly on the speed of communication via magnetic fields (e.g. Alfven waves) and shear viscosity, which, under most astrophysical conditions, are considered to be different.

In fact, the under- and overshooting here may indicate that MFs are insensitive to the momentary rotational frequency of the crust, but rather to the activity and dynamics of the matter in the BL.

Extending this analysis to both the Crab and Vela pulsars, the relative time-delays is expected to be: δ τ v i s C r a b / δ τ v i s V e l a ~ ( Δ R C r a b / Δ R V e l a ) 2 3.4 or equivalently, the undershooting in the case of the Crab is expected to last 3.4 sec comapred to one second in the Vela case.

Finally, although the physics is entirely different, the situation here is strikingly similar to action of the solar dynamo, which is considered to be located in the so-called tachcline between the rigid-body rotating core and the overlying convection zone [16].


The calculations have been carried out using the computer cluster of the IWR, University of Heidelberg. RS acknowledges the use of KAUS baseline research funds

Cite this paper
Hujeirat, A. and Samtaney, R. (2020) How Massive Are the Superfluid Cores in the Crab and Vela Pulsars and Why Their Glitch-Events Are Accompanied with under and Overshootings?. Journal of Modern Physics, 11, 395-406. doi: 10.4236/jmp.2020.113025.

[1]   Roy, J., Yashwant Gupta, Y. and Lewandowski, W. (2012) MNRAS, 424, 2213-2221.

[2]   Yu, M., Manchester, R.N., Hobbs, G., et al. (2013) MNRAS, 429, 688-724.

[3]   Yuan, J., Kou, F. and Wang, N. (2019) AIP Conference Proceedings, 2127, Article ID: 020004.

[4]   Eya, I.O. and Urama, J.O. (2014) International Journal of Astrophysics and Space Science, 2, 16.

[5]   Espinoza, C.M., Lyne, A.G., Stappers, B.W. and Kramer, C. (2011) MNRAS, 414, 1679.

[6]   Fuentes, J.R., Espinoza, C.M., Reisenegger, A., et al. (2017) A&A, 608, A131.

[7]   Ozel, F. and Freire, P. (2016) Annual Review of Astronomy and Astrophysics, 54, 401.

[8]   Hujeirat, A.A. (2018) Journal of Modern Physics, 9, 4.

[9]   Hujeirat, A.A. and Samtaney, R. (2019) Journal of Modern Physics, 9, 4.

[10]   Haensel, P., Lasota, J.P. and Zdunik, J.L. (1999) A&A, 344, 151.

[11]   Hujeirat, A.A. (2018) Journal of Modern Physics, 9, 51-69.

[12]   Hujeirat, A.A. (2018) Journal of Modern Physics, 9, 4.

[13]   Hujeirat, A.A. and Thielemann, F.-K. (2009) MNRAS, 400, 903.

[14]   Ashton, G., Lasky, P.D., Graber, V. and Palfreyman, J. (2019) Rotational Evolution of the Vela Pulsar during the 2016 Glitch.

[15]   Hujeirat, A.A. and Thielemann, F.-K. (2009) A&A, 496, 609.

[16]   Camenzind, M. (2007) Compact Objects in Astrophysics. Springer, Heidelberg.