The Crab Nebula and its pulsar PSR B0531 + 21 is the space laboratory for the study of fundamental physical processes. From the beginning of the current century up to date, the giant pulses of the pulsar in the Crab Nebula are regularly observed on the Large Phased Array of the Pushchino Radio Astronomy Observatory of the Lebedev Physical Institute at a frequency 111 MHz. The connection between the scattering time scale of radio pulses and the dispersion measure of the pulsar, which was established earlier together with Jodrell Bank, is confirmed. The observed variations in the scattering of radio pulses and their partial correlation with the dispersion measure are explained by the eclipse of the pulsar by dense plasma clouds with fluctuations in the electron density significantly exceeding the corresponding fluctuations in the interstellar medium. The question of a possible connection between period failers (glitches), dispersion measure variations, radio pulses scattering and gamma-ray flares is discussed.
The Crab Nebula was first observed in 1731 by the English physicist and amateur astronomer John Bevis. It is included in the catalog of nebulae compiled by Charles Monsieur in 1758 as object M1. The Crab Nebula was so named by Lord Ross in 1844. Observing the nebula with high resolution, he noticed a structure similar to a crab claw. It is believed that it is a remnant of the hydrogen-rich type II supernova explosion, observed in China in 1054. In 1949, radio astronomers Bolton, Stanley, and Slee identified the Crab Nebula as a radio source. Buhler and Blandford [1] write: “Due to its high luminosity L ≈ 1.3 × 1038 erg∙s−1 (Hester 2008) and its proximity of ~2 kpc (Trimble 1973), the Crab can be studied in great detail. It is therefore one of our prime laboratories to study non-thermal processes in the Universe. Many discoveries have been made in the Crab and then been seen in other non-thermal sources (including active galactic nuclei, A gamma ray bursts and X-ray binaries) such as polarized synchrotron radiation or pulsed optical emission (Shklovsky 1953, Cocke et al. 1969). True to this legacy two remarkable discoveries have been made in the last years. Very high-energy (VHE & 100 GeV) gamma-ray emission has been detected from the pulsar, and high-energy (HE & 100 MeV) gamma-ray flares have been discovered from the nebula. These phenomena have not been observed in any other pulsar wind nebula to date”.
The Crab Nebula pulsar PSR B0531 + 21 (Crab Pulsar) was discovered by Staeilin and Reifenstein in 1968 as a radio source of giant pulses (GPs)1, rather than as a periodic source [2]. Fourier analysis did this later. The phenomena of scattering time scale τsc2 and the dispersion measures (DM) variations of the GPs in the Crab Pulsar were first recorded by Rankin and Counselman [3] from observations in Arecibo in the range of 74 - 430 MHz in 1969-1971. Analyzing the scattering variations, they concluded that these processes occur mainly within the Crab Nebula. Isaacman and Rankin [4] continued the measurements until 1974. They also found variations in the scattering time scale and the DM, but found no significant relationship between these phenomena.
Since the beginning of the current century and up to the present time, regular observations of GPs from the Crab Pulsar at 111 MHz have been carried out on the Large Phased Array of the Pushchino Radio Astronomy Observatory (LPA PRAO) of the Lebedev Physical Institute. This article systematizes the data obtained during the observation of the Crab Pulsar [5] [6] [7]. In particular, we show that the scattering of GPs of the Crab pulsar depends not only on the DM, but also on the degree of plasma inhomogeneity. The possible relationship between step changes in the rotation speed (“glitch”), the DM variations and scattering time scale is discussed. Results of the search for a connection between gamma-ray flares and disturbances in the pulsar magnetosphere are given.
2. Observations and Results
Systematic observations of the Crab Pulsar at a frequency of 111 MHz are carried out on LPA radio telescope from September 2002 up to date. The observations were initially carried out using a 128-channel spectrum analyzer with a channel bandwidth of 20 kHz, a total bandwidth 2560 kHz, and a readout interval 2.56 ms. Beginning from April 2006 the observations were conducted using 512-channel digital spectrum analyzer with a channel bandwidth of 5 kHz and readout interval 2.56 ms. The linear polarization of E.-W. Has been registered. GPs are analyzed using a special program that allows us to determine the scattering time scale pulsar radio pulses by simulating the passage of a pulse through a scattering medium. A theoretical scattering pulse profile was generated by successively convolving a gaussian profile, representing the intrinsic pulse, with: 1) an impulse scattering function of the form exp(−t/τsc); 2) a gaussian function representing the effects of interstellar dispersion across the filter bandpass; 3) an exponential function representing the post-detector RC time constant. The theoretical profile was least-squares fitted to the observed profile, with τsc as the adjustable parameter [8].
The scattering values for each GP were averaged over five days of observations. The mean scattering parameters were compared with dispersion measures derived from data obtained at Jodrell Bank Observatory (JBO)3. The monthly JBO data were interpolated to the comparison data.
Preliminary data on the correlation between the variations of scattering time scale τsc and the DM were obtained at the LPA PRAO by Kuzmin and Losovsky at a frequency of 111 MHz together with Jordan and Smith from the JBO in 2006-2007 [9]. However, due to the relatively short duration of this coincidence and the shift between the compared phenomena (the DM variations are 30 days late relative to the scattering variations), it was impossible to exclude the randomness of this event. It was necessary to continue the measurements to make sure that the correlation really took place.
In the course of observations in the period 2008-2013, a long non-stationary process was recorded, characterized by abnormally strong scattering time scale τsc reaching up to 115 ms (Figure 1). We can estimate the size of the massive plasma cloud that eclipsed the pulsar in 2008-2013. It is known that the DM (pc/cm3) = ne(cm−3 ) × D(pc), where ne—mean electron density and D—cloud size. The scattering value at a frequency of 111 MHz in January 2008 with DM = 56.763 was 20 ms. In December 2012, with the DM = 56.855, the scattering τsc was about 115 ms. Electron density of the filament ne, according to published data, is about 103 cm−3 [10]. To increase the DM by 0.092 pc/cm3, the cloud size “D” should be 0.092/103 = 0.92 × 10−4 pc = 2.84 × 1014 cm = 19 AU. During the observation period (2002-2020), there were several glitches4.
During the same period, a number of Gamma-ray flares > 100 Mev has been detected by Fermi Gamma-ray Space Telescope [11] [12] [13] and AGILE-X-ray and gamma-ray astronomical satellite [14] [15] [16] [17] (See Figure 1: dash-line segments).
Figure 2 shows the statistical relationship between the scattering time scale τsc (in milliseconds) and the conditional Crab Nebular dispersion measure, dm = (DM − 56.7) × 1000 in pc/cm3 of the Crab Pulsar at a frequency of 111 MHz in the period 2002-2020. The error of the scattering time scale measurement is estimated at 5%. The relationship between the scattering time scale and dispersion measure is given by expression
lg ( τ s c ) = − 2.51 ( 8 ) + 2.19 ( 4 ) × lg ( d m ) or τ s c = C × d m ( 2.19 ± 0.04 ) (1)
A comparison of the scattering time scale variations and the conditional dispersion measure dm = (DM − 56.7) × 1000 confirms their dependence (Figure 3). The correlation coefficient is 0.80 ± 0.05. The results obtained are in good agreement with measurements of the radio emission of the Crab Pulsar at a frequency of 610 MHz on the 42-foot radio telescope of the JBO and at a frequency of 1400 MHz on the Lovell radio Telescope. Based on the results of observations in the period from 2008 to 2012, the dependence between the τsc and the DM with a correlation coefficient of 0.56 ± 0.01 was obtained [18].
The scattering time scale τsc depends on: the electron density ne; the relative inhomogeneity of the medium Δne/ne; the dispersion measure DM, (or the distance between the source and the observer d); the frequency of the radio signal ν [10].
τ s c ~ ( d * Δ n e ) 2 / ν 4 = [ ( Δ n e / n e ) × D M ] 2 / ν 4 (2)
In an environment, where the parameter Δne/ne is almost constant, the scattering at a given frequency is proportional to the square of the DM:
τ s c ~ C × D M 2 (3)
This is what is observed for most pulsars [19]. The scattering variations in the first approximation are proportional to the variations of DM:
Δ τ s c ~ C 1 × Δ D M (4)
It follows from the Formula (2) that the scattering can be different with the same DM.
In the Crab Nebula, the parameter (Δne/ne) is subject to significant changes, so the variations of the τsc parameter depend not only on the dispersion measure, but also on the degree of plasma inhomogeneity (Δne/ne):
Δ τ s c ~ C 1 × Δ ( D M ) + C 2 × Δ ( Δ n e / n e ) (5)
The spectrum of electron density fluctuations
F ( q ) = C n 2 × q − α , (6)
where C n 2 (m−α−3) is the turbulence coefficient [20]; q-spatial frequency, q = 2π/l—wherel is the character scale of turbulence α-spectral index of the electron-density irregularity.
C n 2 = 2 × 10 − 3 ν β d − β / 2 Δ ν − ( β − 2 ) / 2 , (7)
where β—spectral index of pulse broadening β = 2 α / ( α − 2 ) .
Here ν is the frequency in GHz, Δν is the decorrelation bandwidth in MHz, d is in kpc.
Measurements of GPs from the Crab Pulsar over a wide range of frequencies from 40 to 2228 MHz show that the frequency dependence of scatter-broadening of the pulses is in agreement with a model with normal Gaussian distribution of the inhomogeneities (α = β = 4) [21], and
C n 2 = 2 × 10 − 3 ν 4 d − 2 Δ ν − 1 . (8)
The scattering τsc and the decorrelation bandwidth Δν are related by the ratio 2πτscΔν = 1. In particular, for the frequency ν = 111 MHz and τsc = 1/(2πΔν) turbulence level
C n 2 = 0.00198 τ s c d − 2 , (9)
whereτsc is in milliseconds.
The period from 2002 to 2007 with a small dispersion measure and small scattering variations is conditionally considered as quiet and the period from 2009 to 2013 with a large dispersion measure and significant scattering variations is considered as active.
We accepted mean scattering parameter τsc = 10 ms for the passive period and τsc = 100 ms for the active period. Error estimated as 20%.
According to Yao et al. [22] distance to Crab pulsar d = 1.311 kpc (Err = 14%)
In 2002-2007, during a quiet period C n 2 = 0.0115 ± 0.0038 m−7. The variations of DM measure ΔDM were 0.035, which at a radius of nebular 1.7 pc corresponds to Δne ~ 0.02 cm−3 and this is consistent with the estimate [18].
If, according to [23], the value of the dispersion measure of the Crab nebula can be taken as 4 pc/cm3, then ne ~ 2 cm−3 and Δne/ne ~ 0.01.
In 2009-2013, during the active period C n 2 = 0.115 ± 0.038 m−7. The variations of the DM were ΔDM = 0.105, which corresponds to the average value of the density variations Δne ~ 0.06 cm−3 and Δne/ne ~ 0.03.
In the interstellar medium in the direction of the Crab Pulsar [18] Δne = 2 × 10−6 cm−3.With ne ~ 0.03 cm−3 Δne/ne ~ 1.0 × 10−4.
Turbulence in the disturbed region higher than in quiet ones and much higher than in the interstellar medium.
3. Discussion
Pulsars gradually slow down as a result of the conversion of rotational energy into the energy of charged particles and electromagnetic waves. In the simple model of pulsar rotation frequency f, the first derivative of the rotation frequency df/dt and the second derivative of the rotation frequency d2f/dt2 are used. In this case power-low slowdown is the form df/dt = −k × fn. The value of exponent of f is known as braking index n. The value of the braking index
n = ( f × d 2 f / d t 2 ) / ( d f / d t ) 2 (10)
The braking index n is expected to be 3 in the case of braking by pure dipole radiation and 1 in the case of pulsar wind dominated torque. Deviations from value 3 for Crab pulsar (n ~ 2.5) indicates, than the part of the torque on the pulsar is due to the outflow of particles.
Jumping irregularities in the rotation of the pulsar (glitches) can be of two types: discrete glitches and slow glitches. Both types of glitches cause accelerated rotation of the pulsar, which occurs against the background of secular deceleration of the neutron star. With discrete glitches, the frequency increases suddenly, and then follows an exponential decrease in frequency to the previous value. Slow glitches are associated with slow frequency fluctuations.
The hypothesis about the connection of random irregularities of period with processes in the magnetosphere was firstly expressed Manchester and Taylor in the book “PULSARS” [10]. There are some models, describing the origin of glitches. One model treats glitches as star quakes, caused by the rearrangement of the flattened crust, which tends to become spherical, as the star’s rotation slows down. Another model considers a neutron star as a reservoir, filled with a superfluid, the mass of which, when the pulsar’s rotation slows down, transmits the angular momentum to the crust, which leads to a glitch [24].
The problem of the relationship between glitches, scattering changes, dispersion measures and gamma-ray bursts is discussed in detail Cadez et al. [25]. They devoted their work to the study of the pulsar in the Crab nebula, using radio data from the Jodrell Bank Observatory from 1998 to 2014 together with optical observations of the pulsar, using the ultra-fast photon counter, mounted on the Copernico telescope at the Astrophysical Observatory in Asiago (Italy) in 2008 and at the Southern European Telescope at La Silla observatory (Chile) in 2009. Analysis of the data showed, that jumps in the braking index are associated with large glitches, when the relative change in the rotational speed of the pulsar exceeds 10−8. A similar conclusion was obtained by Lyne et al. with the mean value n = 2.34 [26]. Different braking index values 2.45 and 2.57 are also reported between glitches by Kou and Tong [27].
The delay in the variations of the dispersion measure relative to the variations of the braking index is ~1010 days and is explained by the time of ionization of the nebula by the pulsar wind (Figure 4).
The model of Cadez et al. [25] suggests, that glitches and subsequent variations in the braking index are caused by instability in the magnetosphere, which changes the configuration of the magnetic field and the currents in the plasma, through which the pulsar interacts with the nebula. Perturbations in the pulsar’s magnetosphere can lead to glitches and jumps in the braking index, to the ejection of charged particles into the nebula, to the formation of filaments and to an
increase in the dispersion measure. Intermediary models have been discussed by Kou and Tong [27]. They consider that glitches are related to changes in interior of the neutron stars, but may lead to some effect in the outer magnetosphere. In turn, plasma instabilities in the nebula can disrupt the configuration of magnetic field lines, when fields of opposite polarity are pressed together, and cause gamma-ray flare, similar to flares in the solar corona. The same opinion says Buhler and Blangdford [1]: “The gamma-ray flares are therefore likely connected to explosive reconnection events triggered by current instabilities”. Striany et al. also support this idea [28], but they note: “However, evidence for magnetic field reconnection events in the Crab Nebula is elusive, and no optical or X-ray emission in coincidence with the gamma-ray flaring has been unambiguously detected to date (e.g., Weisskopf et al. 2012)”. According to our data (Figure 3) gamma-flares have tendency to concentrate during the period of enhanced disturbances in the Crab Nebula, that can be considered as support of this model. As for the correlation between of pulsar glitches and gamma-ray flares Buhler and Blandford [1] write: “The time scale of the recurrence of pulsar glitches is similar to the recurrence of the gamma-ray flares, however, there is no obvious correlation in time between these two events”.
4. Conclusions
The scattering of the Crab pulsar radio emission at a frequency of 111 MHz according to our measurements and the dispersion measure provided in Jodrell Bank, show secular variations they are well correlated (correlation coefficient 0.80 ± 0.05).
The observed changes in the pulsar pulse scattering and the dispersion measure and their partial correlation are explained by the passage of plasma clouds in front of the pulsar with a size of about 0.92 × 10−4 pc (3.0 × 1014 cm) with varying degrees of turbulence.
In 2002-2007, during the quiet period, in the Crab Nebula C n 2 = 0.0115 ± 0.0038 m−7 and Δne/ne ~ 0.01.
In 2008-2013, during the active period, in the Crab Nebula C n 2 = 0.115 ± 0.038 m−7 and Δne/ne ~ 0.03.
Turbulence in the disturbed region is higher than in quiet ones and much be higher than in the interstellar medium.
Various types of pulsar radio emission irregularities can be explained by physical processes occurring in the pulsar’s magnetosphere. Perturbations in the pulsar’s magnetosphere can lead to ejection of charged particles into the Nebula, lead to instability of the plasma, disrupt the configuration of the magnetic field lines and cause gamma-ray flares. According to our data gamma-flares have tendency to take place during the period of enhanced disturbances in the Crab Nebula that can be considered as support of this model. New observations are needed to support this hypothesis.
Funding
This work was supported by the research program of the Presidium of the Russian Academy of Sciences “Non-stationary phenomena in the objects of the Universe”
Acknowledgements
The author thanks the staff of the meter wave-length technical service laboratory of the Pushchino Radio Astronomy Observatory for the providing reliable observations. The author is very thankful to Managing Editor of IJAA Jane Gao for the fruitful cooperation and the anonymous reviewers for valuable notices.
Conflicts of Interest
The author declares no conflict of interest, regarding the publication of this paper.
Cite this paper
Losovsky, B.Ja. (2021) Pulsar B0531 + 21 in the Crab Nebula. International Journal of Astronomy and Astrophysics, 11, 470-480. https://doi.org/10.4236/ijaa.2021.114024
NOTESReferencesBuhler, R. and Blandford, R. (2014) The Surprising Crab Pulsar and Its Nebula: A Review. Reports on Progress in Physics, 77, Article ID: 066901. https://doi.org/10.1088/0034-4885/77/6/066901Staelin, D.H. and Reifenstein, E.C. (1968) Pulsating Radio Sources near the Crab Pulsar Nebula. Science, 162, 1481-1483. https://doi.org/10.1126/science.162.3861.1481Rankin, J.M. and Counselman, C.C. (1973) Pulsar NP 0532: Variability of Dispersion and Scattering. The Astrophysical Journal, 181, 875-889. https://doi.org/10.1086/152099Isaacman, R. and Rankin, J.M. (1977) The Crab Nebula Pulsar: Variability of Dispersion and Scattering. The Astrophysical Journal, 214, 214-232. https://doi.org/10.1086/155245Kuz’min, A.D., Belyatsky, Yu.A., Dumsky, D.V., Izvekova, V.A., Lapaev, K.A., Logvinenko, S.V., Losovsky, B.Ya. and Pugachev, V.D. (2012) Monitoring of the Radio Emission of the Crab Nebula Pulsar at Low Frequencies. Astronomy Reports, 55, 416-424. https://doi.org/10.1134/S1063772911040032Losovsii, B.Ya. (2017) The Dynamics and Structure of Plasma Inhomogeneities in the Crab Nebula. Astronomy Reports, 61, 187-192. https://doi.org/10.1134/S1063772917030052Losovsky, B.Ya., Dumsky, D.V. and Belyatsky, Yu.A. (2019) Relationship between Scattering and Dispersion Measure of the Crab Nebula Pulsar. Astronomy Reports, 63, 830-834. https://doi.org/10.1134/S1063772919090051Alurkar, S.K., Slee, O.B. and Bobra, A.D. (1986) Scattering of Pulsar Radiation and Electron Density Turbulence in the Interstellar Medium. Australian Journal of Physics, 39, 433-448. https://doi.org/10.1071/PH860433Kuzmin, A., Losovsky, B.Ya., Jordan, C.A. and Graham-Smith, F. (2008) Correlation of the Scattering and Dispersion Events in the Crab Nebula Pulsar. Astronomy & Astrophysics, 483, 13-14. https://doi.org/10.1051/0004-6361:20079211Manchester, R. (1977) Taylor J. Pulsars. Freeman, San Francisco; Moscow, 1980.Abdo, A.A., Ackermann, M., Ajello, M., et al. (2011) Gamma-Ray Flares from the Crab Nebula. Science, 331, 739-755.Ojha, R., Hays, E., Buehle, R. and Dutka, M. (2013) Fermi LAT Detection of a New Gamma-Ray Flare from the Crab Nebula Region. Atel, 4855.Buehler, R., Scargle, J.D., Blandford, R.D., et al. (2012) Gamma-Ray Activity in the Crab Nebula. The Exceptional Flare of 2011 April. The Astrophysical Journal, 749, 26. https://doi.org/10.1088/0004-637X/749/1/26Tavani, M., Bulgarelli, A., Vittorini, V., et al. (2011) Discovery of Powerful Gamma-Ray Flares from the Crab Nebula. Science, 331, 736-759.Striani, E., Tavani, M., Vittorini, V., et al. (2013) Variable Gamma-Ray Emission from the Crab Nebula: Short Flares and Long “Waves”. The Astrophysical Journal, 765, 52S. https://doi.org/10.1088/0004-637X/765/1/52Bulgarelli, A., Tavani, M., Parmiggiani, N., et al. (2018) AGILE Detection of Renewed Gamma-Ray Activity from Crab Nebula. Atel, 12148.Huang, X., Yuan, Q. and Fan, Y. (2021) Ray Flares from the Crab Nebula with Fermi-LAT. The Astrophysical Journal, 908, 65.McKee, J.W., Lyne, A.G., Stappers, B.W., Bassa, C.G. and Jordan, C.A. (2018) Temporal Variations in Scattering and Dispersion Measure in the Crab Pulsar and Their Effect on Timing Precision. Monthly Notices of the Royal Astronomical Society, 479, 4216-4224. https://doi.org/10.1093/mnras/sty1727Kuzmin, A.D. and Losovsky, B.Ya. (1999) Observations of the Mean Profiles of Radio Pulsar and the Structure of Their Emitting Regions at 102 MHz. Astronomy Reports, 43, 288-304.Cordes, J.M., Weisberg, J.M. and Boriakoff, V. (1985) Small-Scale Electron Density Turbulence in the Interstellar Medium. The Astrophysical Journal, 288, 221-247. https://doi.org/10.1086/162784Kuzmin, A.D., Kondrat’ev, V.I., Kostyuk, S.V., Losovsky, B.Ya., Popov, M.V., Soglasnov, V.F., D’Amico, N. and Montebugnoli, S. (2002) Frequency Dependence of the Scattering Pulse Broadening for Crab Pulsar. Astronomy Letters, 28, 251-256. https://doi.org/10.1134/1.1467260Yao, J.M., Manchester, R.N. and Wang, N. (2017) A New Electron-Density Model for Estimation of Pulsar and FRB Distances. The Astrophysical Journal, 835, 29-59. https://doi.org/10.3847/1538-4357/835/1/29Kuzmin, A.D., Losovsky, B.Ya., Logvinenko, S.V. and Litvinov, I.I. (2008) Deviation of the Arrival-Time Delay of Pulses from the Crab Pulsar from the Quadratic Frequency Relation. Astronomy Reports, 52, 910-916. https://doi.org/10.1134/S106377290811005XEspinoza, C.M., Lyne, A.G., Stappers, B.W. and Kramer, M. (2011) A Study of 315 Glitches in the Rotation of 102 Pulsars. Monthly Notices of the Royal Astronomical Society, 414, 1679-1704. https://doi.org/10.1111/j.1365-2966.2011.18503.xCadez, A., Zampiery, L., Barbiery, C., Calvani, M., Naletto, G., Barbiery, M. and Ponikwar, D. (2016) What Brakes the Crab Pulsar. Astronomy & Astrophysics, 587, A99. https://doi.org/10.1051/0004-6361/201526490Lyne, A.G., Jordan, C.A., Graham-Smith, F., Espinoza, C.M., Stappers, B.W. and Weltevrede, P. (2015) 45 Years of the Rotation of the Crab Pulsar. Monthly Notices of the Royal Astronomical Society, 446, 857-864. https://doi.org/10.1093/mnras/stu2118Kou, F.F. and Tong, H. (2015) Rotational Evolution of the Crab Pulsar in the Wind Braking Model. Monthly Notices of the Royal Astronomical Society, 450, 1990-1998. https://doi.org/10.1093/mnras/stv734Striani, E., Tavani, M., Vittorini, V., et al. (2013) Variable Gamma-Ray Emission from the Crab Nebula: Short Flares and Long “Waves”. The Astrophysical Journal, 765, 52-62. https://doi.org/10.1088/0004-637X/765/1/52