Next Article in Journal
Charged Particle Pseudorapidity Distributions Measured with the STAR EPD
Next Article in Special Issue
Galaxy Rotation Curve Fitting Using Machine Learning Tools
Previous Article in Journal
De Sitter-Invariant Black Holes
Previous Article in Special Issue
Neutron Star Binaries Produced by Binary-Driven Hypernovae, Their Mergers, and the Link between Long and Short GRBs
 
 
Correction published on 18 June 2024, see Universe 2024, 10(6), 265.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimates of the Surface Magnetic Field Strength of Radio Pulsars

by
Vitaliy Kim
1,2,*,
Adel Umirbayeva
1,3,† and
Yerlan Aimuratov
1,3,4,*,‡
1
Fesenkov Astrophysical Institute, Observatory 23, Almaty 050020, Kazakhstan
2
Pulkovo Observatory, Pulkovoskoe Shosse 65/1, 196140 Saint Petersburg, Russia
3
Department of Theoretical and Nuclear Physics, Faculty of Physics and Technology, Al-Farabi Kazakh National University, Al-Farabi Avenue 71, Almaty 050040, Kazakhstan
4
International Center for Relativistic Astrophysics Network, Piazza della Repubblica 10, 65122 Pescara, Italy
*
Authors to whom correspondence should be addressed.
Currently a Master student at Al-Farabi KazNU.
Currently a Postdoctoral associate at Al-Farabi KazNU.
Universe 2023, 9(7), 334; https://doi.org/10.3390/universe9070334
Submission received: 2 June 2023 / Revised: 12 July 2023 / Accepted: 13 July 2023 / Published: 14 July 2023 / Corrected: 18 June 2024
(This article belongs to the Special Issue Remo Ruffini Festschrift)

Abstract

:
We investigate the geometry of the magnetic field of rotation-powered pulsars. A new method for calculating an angle ( β ) between the spin and magnetic dipole axes of a neutron star (NS) in the ejector stage is considered within the frame of the magnetic dipole energy loss mechanism. We estimate the surface magnetic field strength ( B ns ) for a population of known neutron stars in the radio pulsar (ejector) stage. The evaluated B ns ( β ) may differ by an order of magnitude from the values without considering the angle β . It is shown that B ns ( β ) lies in the range 10 8 10 14 G for a known population of short and middle periodic radio pulsars.

1. Introduction

Radio pulsars are fast-spinning magnetized neutron stars (NS) demonstrating regular modulations (pulsations) of their radiation with a high stable period in the radio range. The axis of the magnetic field of the radio pulsar and its spin axis are not aligned, and the beam of radiation is emitted in a cone-shaped region (see Figure 1). Therefore, pulsar radiation is seen as pulses (beacon effect) by an external observer [1].
Radio pulsars are characterized by the rapid axial rotation (or spin) they have acquired due to the conservation of angular momentum during their formation. Their spin periods ( P s ) lie in a wide range: from 0.0014 s to 23.5 s [2,3], with the majority not exceeding a few seconds. Radio pulsars are usually divided into several groups depending on their spin period [4,5]:
  • Short-periodic pulsars, including millisecond pulsars ( P s < 0.1   s );
  • Middle-periodic pulsars ( 0.1   s < P s < 2   s );
  • Long-periodic pulsars ( P s > 2   s ).
According to [5], a pulsar wind mainly causes the spin-down process of long-periodic radio pulsars. However, for short and middle periodic radio pulsars, a primary mechanism of rotation energy loss ( E ˙ obs ) is believed to be magnetic dipole radiation (MDR).
The MDR mechanism of the energy loss was first considered in [6,7] for radio pulsars. It was shown that the magnetized NS could lose its rotational energy by MDR generation. This evolution stage of NS is also known as the “ejector” stage, and its energy loss ( E ˙ md ) for the generation of MDR expresses as:
E ˙ md = 2 3 μ ns 2 ω s 4 ( sin β ) 2 c 3 ,
where μ ns = B ns R ns 3 / 2 is the magnetic dipole moment of NS, R ns is the radius of NS, ω s = 2 π / P s is the angular rotation velocity, c is the speed of light, and β is the angle between spin and magnetic dipole axes; the value of β lies within 0–90 deg.
The expression for rotational energy loss is the following:
E ˙ obs = I ω s ω ˙ s = I 4 π 2 P ˙ s P s 3 ,
where I is a moment of inertia of NS, ω ˙ s = 2 π P ˙ s / P s 2 is a derivative of the angular rotation velocity, P ˙ s is a derivative of the spin period, i.e., rotational spin-up or spin-down.
Solving the system of Equations (1) and (2) by equating the losses, with μ ns and using canonical values for NS (see Section 2.2), one can derive an expression for B ns as
B ns sin β = 3.2× 10 19 P s P ˙ s G .
Equation (3) is a basic expression for estimating the magnetic field strength for rotation-powered pulsars. As seen from the equations above, E ˙ md significantly depends on angle β , reaching the maximal energy loss E ˙ md ( max ) with orthogonal axes ( β = 90 ). Indeed, in most cases, the magnetic field strength for radio pulsars is estimated by accepting the E ˙ md ( max ) case; however, angle β can vary widely from the maximum value. Thus, it is crucial to correctly estimate the angle between spin and magnetic dipole axes to evaluate B ns ( β ) for rotation-powered pulsars.
Various methods for estimating the β -parameter have been previously proposed in the literature [8,9,10,11]. Here, we offer a relatively simple method based on a geometric approach for calculating the angle between the spin and magnetic dipole axes of a neutron star (NS) in the ejector stage. Section 2.1 outlines a basic concept and geometry for β . On its basis, in Section 2.2, we evaluate the surface magnetic field strength B ns ( β ) for a population of known neutron stars from the ATNF pulsar catalog. We provide the main results with a discussion in Section 3 and Section 4.
The obtained data can help study properties and geometry of NS magnetic fields [12], study and model pulsar spin evolution [13], investigate stellar evolution in the late stages [14], etc.

2. Methods

2.1. Estimation of β -Parameter

According to [15], an equation of magnetic field lines, based on an assumption of a dipole magnetic field, in polar coordinates expressed as
r = r e ( cos ϕ ) 2 ,
where r e is an equatorial radius of the magnetic field, corresponding to NS magnetosphere radius (see Figure 2), ϕ is an angle measured from the magnetic equator r e towards the magnetic pole.
Using an Equation (4), it is possible to find an opening angle of emission cone γ assuming that the radio pulsar magnetosphere is limited by a light cylinder with radius R LC , i.e., r e = R LC = c × P s / 2 π [16]. At the base of the emission cone, the radius vector of the magnetic field line corresponds to the NS radius r b = R ns .
Solving Equation (4), we can find ϕ b corresponding to the angle between the magnetic equator and lateral surface of the emission cone (see Figure 2).
ϕ b = arccos R ns × 2 π c × P s 1 / 2 .
Using Equation (5), one can find an opening angle of the emission cone γ :
γ = 2 × ( 90 ϕ b ) .
At the next step, one can consider the dihedral angle (Figure 3) formed by two emission cone guides d and the rotation axis ω s . Then, the linear diameter a of the emission cone at a distance d will be expressed as follows (via the cosine theorem):
a 2 = 2 d 2 ( 1 cos γ ) a 2 = 2 ( d sin β ) 2 ( 1 cos ε ) .
Solution of the Equation (7) with respect to parameter sin β gives
sin β = 1 cos γ 1 cos ε .
The angle ε can be estimated by means of a pulse profile from observations of radio pulsars (see Figure 3 and Equation (9)). We suppose that a diameter a of the emission cone at a distance d covers the circle (dotted) around the larger cone formed by the magnetic axis μ rotating around the spin axis ω of the NS. The diameter of the cone a can be expressed in terms of the width of the observed pulse profile since the start and end times of the passage of the base of the emission cone through the observer correspond to the beginning (rise) and end (fall) of the curve in the pulse profile. We use data of the w 10 parameter from the pulsar catalog, which is a pulse width at 10% of the peak of intensity [17] supposing that w 10 approximately corresponds to the “travel” time (t) of the emission cone passing through an observer
ε = ω s t = 2 π P s t 2 π P s w 10 .

2.2. Data Selection and Evaluation of B n s ( β )

We use data from the ATNF pulsar catalog, available online1 to apply our method to interested pulsars. The catalog contains information on rotation-powered pulsars and counts for over 3000 objects, and it is the most extensive database that provides information about radio pulsars. The catalog is maintained by the Australia Telescope National Facility (ATNF). The catalog includes detailed information about pulsars, such as their positions, rotation periods, spin-down rates, dispersion measures, and other relevant parameters. The ATNF pulsar catalog was initially compiled using data from the Parkes radio telescope in Australia. Over time, data from other radio telescopes in Australia and worldwide were incorporated into the catalog. The catalog is regularly updated as new observations are made, and new pulsars are discovered; see [2] for more details.
The data selection was carried out according to the following criteria:
  • Objects with known spin period P s ;
  • Exclusion of objects with P s > 2 s ;
  • Known spin-down rate P ˙ s ;
  • Known w 10 parameter.
The sample resulted in 1468 objects from the ATNF pulsar catalog, which are NS with P s < 2 s and the above parameters presented.
Some parameters of NS, such as radius R ns , mass M ns , and moment of inertia I are in a narrow range close to canonical values; their variations should influence the evaluation insignificantly [1,15]. Therefore, in our calculations, we use canonical values of mass, radius, and moment of inertia as R ns 10 6 cm , M ns 1.4 M , I 10 45 g cm 2 , correspondingly.

3. Results

Table 1 (in full available in a machine-readable format) summarizes our calculations of β and B ns ( β ) parameters for the chosen population. The statistics on calculated β and B ns ( β ) is given in Table 2 and Table 3.
In the latter tables, we subdivided radio pulsars into three categories according to their spin periods to clarify and further underline the difference in evolutionary stages in Section 4.
For the population of middle-periodic pulsars ( 0.1   s < P s < 2   s ) counting to 1301 known objects, their values of B ns ( β ) lie within the range ∼ 10 10 10 14 G . Average values of B ns ( β ) are in good agreement with the canonical value of magnetic field 10 12 10 13 G for radio pulsars [12]. Their β lie in the wide range 0.61 47.51 deg , but for most cases does not exceed 10 deg, the median value of β for the population of middle-periodic pulsars corresponds to 6.98 deg.
Short-periodic pulsars ( P s < 0.1   s ), including 94 millisecond pulsars, altogether count to 167 known objects. Their B ns ( β ) values lie in the wide range ∼ 10 8 10 14 G , but when considering millisecond pulsars only, their B ns ( β ) cover ∼ 10 8 10 10 G range with 5.8× 10 8 G median value. Unlike the population of middle-periodic pulsars, millisecond objects have large values of β -parameter lying within the range 9.41 68.78 deg .
We built Figure 4 and Figure 5 on derived values of B ns ( β ) and β to show their general trend in relation to pulsars’ spin periods. We distinguish the above three groups by vertical lines on both plots. In Figure 4, we used two data sets as blue dots (1468 objects) for the calculated B ns ( β ) and gray dots (1579 objects) for B ns retrieved from the ATNF pulsar catalog with β fixed at 90 deg. Here, we sharply cut off data points with P s > 2 s, thus eliminating long-periodic pulsars. As mentioned in Section 1, the primary mechanism of their rotational energy loss is the generation of pulsar wind [4,5]. In such a case, the E ˙ value does not depend on the β . Within this approach, the magnetic field of the NS can be estimated by knowing the power of the ejected pulsar wind, which cannot be estimated directly from observations. Therefore, the estimation of the magnetic fields of long-periodic radio pulsars ( P s > 2   s ) is a model-dependent task and is beyond the scope of the current article.
Figure 5 shows the distribution of 1468 derived β -parameters according to their spin periods. Two distinct features can be noticed: (a) some dots pull in a chain showing the positive trend sequences positioned parallel to each other, and (b) dots distribution generally goes above some level, here marked as a solid black line. Both features are related to the w 10 -parameter, where the former points to the objects with similar values of w 10 . At the same time, the latter peculiarity indicates that all pulsars in our sample obey the condition w 10 P s / 2 , i.e., each pulse duration does not exceed half of the spin period.

4. Discussion

In our paper, we used the classical dipole model of the radio pulsar magnetosphere proposed by [16]. In this case, the magnetosphere of a neutron star has a dipole structure co-rotating with a pulsar. It is limited by the so-called light cylinder on which the linear velocity of the magnetic field lines reaches the speed of light. This model is canonical and relevant to this day [18].
Indeed, the width of the pulse profile of radio pulsars can vary depending on the frequency (wavelength) of the observed flux. However, significant deviations in the profile width are observed at lower frequencies (<200 MHz). According to [19], this phenomenon is present because the light cone becomes wider when observed at lower frequencies, thereby seeing areas further from the pulsar’s surface where the opening angle of the closed magnetic field lines is becoming broader. However, for higher frequencies (>200 MHz), this effect can be neglected [19]. The values of w 10 in the ATNF Pulsar Catalog are the average profile width in the range of frequencies between 400–2000 MHz.
We showed that estimates of the surface magnetic field strength ( B ns ) for a population of known neutron stars in the radio pulsar (ejector) stage should depend essentially on the angle β between spin and magnetic dipole axes of a neutron star. These estimates may differ by order of magnitude from those without considering the angle β (see Figure 4). The proposed method can be used when considering only rotation-powered pulsars with a MDR energy loss mechanism. This is not the case for the long-periodic pulsars with P s > 2 s; therefore, we sharply cut off such objects in Figure 4, although borderline transition cases may occur individually.
Within the framework of the proposed technique, it is not possible to estimate the evolution of β over time since these changes are associated with changes in the flow of currents in the core of NS and the interaction of the magnetosphere with the surrounding plasma [20]. Nevertheless, we can compare our results against angles obtained within the framework of other methods.
As was mentioned in Section 1 there are several approaches for estimating the β -parameter. They can be conditionally divided into two groups: geometric and polarimetric methods. The first is based on different geometric models for NS magnetic field and emission cone. Our method also belongs to the first group. The second is based on measuring the position angle of linear polarization from radio pulsars, which depends on β [21]. Interest in comparing β from these two approaches resulted in the following consideration.
In recent articles [5,11,22] an estimation of β -parameter was obtained within geometric method based on spherical trigonometry. A polar cap model was used by authors with an assumption that the line-of-sight passes through the center of the emission cone. The comparison between our data and data from [5,11,22] is shown in Figure 6 for matched 1242 and 246 radio pulsars and their statistics are given in Table 4. In most cases, the difference in estimates ( Δ β median) does not exceed 5 deg and is mainly caused by the difference in the methods (models) used. Negative and positive trends can be noticed correspondingly between our data and data by Ken’ko et al. 2023 [5] (Figure 6, left panel), with a vertical dotted line approximately marking the spin period where two methods give similar β estimation. Again, this is due to differences in geometric approaches for β estimation since the data themselves for both methods were taken from the same catalog (ATNF Pulsar Catalogue). No trends are seen between our data and the data by Nikitina et al. 2017 [22] (Figure 6, right panel), where blue dots are systematically positioned above red ones for matches pulsars. While these authors use analogous methods based on spherical trigonometry, their sample is relatively small, so the trends may not have enough data to manifest. Another reason could be in the data themselves, since in [22] the authors have used data from their observational facilities (Pushchino Radio Astronomy Observatory).
We further attempted to compare estimates obtained by polarimetric studies to determine the angle β performed for only a small part of the radio pulsar population (see Table 5). This method is based on measuring the position angle of linear polarization and is more reliable than geometric approaches. However, for some objects, when observed in different wavelength ranges (frequencies), it can give a significant scatter, especially for larger β . For example, as shown in [22] for PSR B1055-52 (aka J1057-5226) β -parameter estimation at 10 cm wavelength gives β 10 cm = 15 deg, but estimation at 20 cm wavelength gives β 20 cm = 24 deg; for PSR B1702-19 (aka J1705-1906) β 10 cm = 49 deg, β 20 cm = 70 deg, etc. The larger scatter in | β pol β geom | between geometric and polarimetric methods is mainly due to the assumption that the line-of-sight passes through the center of the base of the emission cone. Thus, the geometric estimates are the lower limits for the measured angle β [5].
As also seen in Figure 4, we obtained 110 objects (7.5% from 1468 pulsar sample) with estimated magnetic fields exceeding the so-called quantum critical threshold ∼ 4.4× 10 13 G [23]. These are blue dots over the solid red line, and all (except one short-periodic source PSR B0540-69) belong to the population of middle-periodic pulsars. The maximal value B ns ( β ) 7.56× 10 14 G refers to the pulsar PSR J1119-6127. According to [24] this radio pulsar demonstrates episodic SGR-like high-energy bursts reaching 2.8× 10 39 erg s 1 within 15– 150 keV range. The magnetic field of the NS derived from analysis of PSR J1119-6127 during its burst activity corresponds to B ns 10 14 G [24] that agrees with our estimate within an order of magnitude. The analysis of the rest of the high-B sub-sample can be interesting from the point of a possible relation between high-B radio pulsars and the population of isolated X-ray pulsars [25]: anomalous X-ray pulsars (AXP), soft gamma-repeaters (SGR), etc.
As seen in Figure 5 and Table 1, with increasing spin period P s , there is a tendency of the angle β to decrease. This agrees with the current view of the spin evolution of NS [20]: older neutron stars have lengthier spin periods and smaller values of β , excepting a millisecond pulsar population. According to [20] on the timescales 10 6 10 7 yr in the ejector stage a NS should align its magnetic and spin axes, i.e., the angle β tends to zero.
For the population of millisecond pulsars (MSPs), the evolution of the β -parameter may differ significantly from other radio-pulsar populations. The millisecond pulsars are neutron stars in close binary systems or descendants of close binary systems in the case of isolated MSPs, with a low-mass companion, where accretion flow from a normal companion recycled a NS to ultra-short spin periods [26]. Thus, MSPs are old neutron stars whose rotational evolution has gone all possible stages (ejectorpropelleraccretor) and then came back to the ejector stage through accretion recycling [26].
According to [27], the initial ejector stage for a neutron star in a binary system (with a normal star companion) lasts 10 5 10 6 yr , that is much shorter than in the case of an isolated NS and order of magnitude shorter than the timescale needed for aligning magnetic and spin axes of NS in ejector stage (see previous paragraph). Therefore, a NS in a binary system can move on to the following evolutionary stages (propeller and accretor) from the ejector stage with a β -parameter, which is significantly different from a zero value. Moreover, according to [28], the magnetic and spin axes of a neutron star in the stage of accretion tend to an orthogonal position, i.e., β -parameter increases to 90 deg on the timescale ∼ 10 5 yr . The maximal possible lifetime of a NS on the accretor stage in a low-mass binary system is comparable to the lifetime of its normal companion, ∼ ( 0.1 10 ) × 10 9 yr [29]. It exceeds the orthogonalization timescale by several orders of magnitude, sufficient to increase the β -parameter significantly. Thus, MSPs are old neutron stars that demonstrate large values of β -parameter compared to other types of radio pulsars in the ejector stage.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/2218-1997/9/7/334/s1, Table S1: Sample of radio pulsars with estimated β and B ns ( β ) . Other parameters (coordinates, P s , P ˙ s , w 10 ) are extracted from ATNF pulsar catalog [17].

Author Contributions

Conceptualization, writing—review and editing, project administration, V.K. and Y.A.; methodology, investigation, writing—original draft preparation, supervision, V.K.; validation, data curation, visualization, A.U.; formal analysis, resources, funding acquisition, Y.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research has been funded by the Science Committee of the Ministry of Education and Science of the Republic of Kazakhstan (Grant No. AP09258811).

Acknowledgments

Authors acknowledge Remo Ruffini and the organizers of the conference “Prof. Remo Ruffini Festschrift” for the invitation to give a talk and the support provided during the online meeting. We thank the Editors and Referees for their comments and suggestions to improve the presentation of the results.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
ATNFAustralia Telescope National Facility
AXPAnomalous X-ray pulsar
degDegree (unit)
GGauss (unit)
MDRMagnetic Dipole Radiation
msMillisecond (unit)
MSPMillisecond pulsar
NSNeutron star
SGRSoft gamma-repeater
sSecond (unit)
s/sSeconds per second (unit)
yrYear (unit)

Note

1

References

  1. Lyne, A.; Graham-Smith, F. Pulsar Astronomy; Cambridge University Press: Cambridge, UK, 2012. [Google Scholar]
  2. Manchester, R.N.; Hobbs, G.B.; Teoh, A.; Hobbs, M. The Australia Telescope National Facility Pulsar Catalogue. Astron. J. 2005, 129, 1993–2006. [Google Scholar] [CrossRef]
  3. Tan, C.M.; Bassa, C.G.; Cooper, S.; Dijkema, T.J.; Esposito, P.; Hessels, J.W.T.; Kondratiev, V.I.; Kramer, M.; Michilli, D.; Sanidas, S.; et al. LOFAR Discovery of a 23.5 s Radio Pulsar. Astrophys. J. 2018, 866, 54. [Google Scholar] [CrossRef]
  4. Malov, I.F.; Marozava, H.P. On Braking Mechanisms in Radio Pulsars. Astron. Rep. 2022, 66, 25–31. [Google Scholar] [CrossRef]
  5. Ken’ko, Z.V.; Malov, I.F. Evolution of the angles between magnetic moments and rotation axes in radio pulsars. Mon. Not. R. Astron. Soc. 2023, 522, 1826–1842. [Google Scholar] [CrossRef]
  6. Gunn, J.E.; Ostriker, J.P. Magnetic Dipole Radiation from Pulsars. Nature 1969, 221, 454–456. [Google Scholar] [CrossRef]
  7. Pacini, F. The Nature of Pulsar Radiation. Nature 1970, 226, 622–624. [Google Scholar] [CrossRef] [PubMed]
  8. Kuz’min, D.K.; Dagkesamanskaya, I.M. Evaluation of the angles of magnetic axis relative to rotation axis of the Pulsars. Astron. Lett. 1983, 9, 149–154. [Google Scholar]
  9. Lyne, A.G.; Manchester, R.N. The shape of pulsar radio beams. Mon. Not. R. Astron. Soc. 1988, 234, 477–508. [Google Scholar] [CrossRef]
  10. Malov, I.F. Angle Between the Magnetic Field and the Rotation Axis in Pulsars. Sov. Astron. 1990, 34, 189. [Google Scholar]
  11. Ken’ko, Z.V.; Malov, I.F. Comparison of the Angles between the Magnetic Moment and Rotation Axis for Two Groups of Radio Pulsars. Astron. Rep. 2022, 66, 669–692. [Google Scholar] [CrossRef]
  12. Igoshev, A.P.; Popov, S.B.; Hollerbach, R. Evolution of Neutron Star Magnetic Fields. Universe 2021, 7, 351. [Google Scholar] [CrossRef]
  13. Xie, Y.; Zhang, S.N. Rotational behaviors and magnetic field evolution of radio pulsars. Astron. Nachrichten 2014, 335, 775. [Google Scholar] [CrossRef]
  14. Portegies Zwart, S.F.; Van den Heuvel, E.P.J. The origin of single radio pulsars. New Astron. 1999, 4, 355–363. [Google Scholar] [CrossRef]
  15. Lipunov, V.M. Astrophysics of Neutron Stars; Springer: Berlin/Heidelberg, Germany; New York, NY, USA, 1992. [Google Scholar]
  16. Goldreich, P.; Julian, W.H. Pulsar Electrodynamics. Astrophys. J. 1969, 157, 869. [Google Scholar] [CrossRef]
  17. Manchester, R.N.; Hobbs, G.B.; Teoh, A.; Hobbs, M. VizieR Online Data Catalog: ATNF Pulsar Catalogue (Manchester+, 2005); VizieR Online Data Catalog: 2016; p. B/psr. Available online: https://ui.adsabs.harvard.edu/abs/2016yCat....102034M/abstract (accessed on 1 June 2023).
  18. Rahaman, S.M.; Mitra, D.; Melikidze, G.I.; Lakoba, T. Pulsar radio emission mechanism—II. On the origin of relativistic Langmuir solitons in pulsar plasma. Mon. Not. R. Astron. Soc. 2022, 516, 3715–3727. [Google Scholar] [CrossRef]
  19. Pilia, M.; Hessels, J.W.T.; Stappers, B.W.; Kondratiev, V.I.; Kramer, M.; van Leeuwen, J.; Weltevrede, P.; Lyne, A.G.; Zagkouris, K.; Hassall, T.E.; et al. Wide-band, low-frequency pulse profiles of 100 radio pulsars with LOFAR. Astron. Astrophys. 2016, 586, A92. [Google Scholar] [CrossRef]
  20. Philippov, A.; Tchekhovskoy, A.; Li, J.G. Time evolution of pulsar obliquity angle from 3D simulations of magnetospheres. Mon. Not. R. Astron. Soc. 2014, 441, 1879–1887. [Google Scholar] [CrossRef]
  21. Radhakrishnan, V.; Cooke, D.J. Magnetic Poles and the Polarization Structure of Pulsar Radiation. Astrophys. Lett. 1969, 3, 225. [Google Scholar]
  22. Nikitina, E.B.; Malov, I.F. Magnetic fields of radio pulsars. Astron. Rep. 2017, 61, 591–611. [Google Scholar] [CrossRef]
  23. Peng, Q.H.; Tong, H. The physics of strong magnetic fields in neutron stars. Mon. Not. R. Astron. Soc. 2007, 378, 159–162. [Google Scholar] [CrossRef]
  24. Göğüş, E.; Lin, L.; Kaneko, Y.; Kouveliotou, C.; Watts, A.L.; Chakraborty, M.; Alpar, M.A.; Huppenkothen, D.; Roberts, O.J.; Younes, G.; et al. Magnetar-like X-Ray Bursts from a Rotation-powered Pulsar, PSR J1119-6127. Astrophys. J. Lett. 2016, 829, L25. [Google Scholar] [CrossRef]
  25. Benli, O.; Ertan, Ü. On the evolution of high-B radio pulsars with measured braking indices. Mon. Not. R. Astron. Soc. 2017, 471, 2553–2557. [Google Scholar] [CrossRef]
  26. Lorimer, D.R. Binary and Millisecond Pulsars. Living Rev. Relativ. 2008, 11, 8. [Google Scholar] [CrossRef] [PubMed]
  27. Urpin, V.; Geppert, U.; Konenkov, D. Magnetic and spin evolution of neutron stars in close binaries. Mon. Not. R. Astron. Soc. 1998, 295, 907–920. [Google Scholar] [CrossRef]
  28. Biryukov, A.; Abolmasov, P. Magnetic angle evolution in accreting neutron stars. Mon. Not. R. Astron. Soc. 2021, 505, 1775–1786. [Google Scholar] [CrossRef]
  29. Kiziltan, B.; Thorsett, S.E. Millisecond Pulsar Ages: Implications of Binary Evolution and a Maximum Spin Limit. Astrophys. J. 2010, 715, 335–341. [Google Scholar] [CrossRef]
Figure 1. Scheme of rotating magnetized neutron star and its axes. In the figure denoted: μ is a magnetic dipole moment of NS, ω s is an angular rotation velocity, β is an angle between spin and magnetic dipole axes, r e is an equatorial radius of the magnetic field of NS, and R LC is a radius of the light cylinder of NS.
Figure 1. Scheme of rotating magnetized neutron star and its axes. In the figure denoted: μ is a magnetic dipole moment of NS, ω s is an angular rotation velocity, β is an angle between spin and magnetic dipole axes, r e is an equatorial radius of the magnetic field of NS, and R LC is a radius of the light cylinder of NS.
Universe 09 00334 g001
Figure 2. Scheme of a neutron star magnetosphere. In the figure denoted: R ns is a radius of the neutron star, r is a radius vector of the magnetosphere, ϕ is an angle measured from the magnetic equator r e towards the magnetic pole, with μ being an axis of a magnetic dipole moment of NS, γ is an opening angle of emission cone, ϕ b corresponds to the angle between the magnetic equator and lateral surface of the emission cone, and r e is an equatorial radius of the magnetic field of NS.
Figure 2. Scheme of a neutron star magnetosphere. In the figure denoted: R ns is a radius of the neutron star, r is a radius vector of the magnetosphere, ϕ is an angle measured from the magnetic equator r e towards the magnetic pole, with μ being an axis of a magnetic dipole moment of NS, γ is an opening angle of emission cone, ϕ b corresponds to the angle between the magnetic equator and lateral surface of the emission cone, and r e is an equatorial radius of the magnetic field of NS.
Universe 09 00334 g002
Figure 3. Scheme of the opening angle of the emission cone and the angle β between spin and magnetic dipole axes (left panel), and its view from the top (right panel). In the figure denoted: a is the linear diameter of the emission cone, and d is the distance. The rest of the symbols are identical to the ones in Figure 1 and Figure 2. Bottom panel shows the pulse profile with period ( P s ) and width of individual pulse at 10 % of maximal intensity ( w 10 ).
Figure 3. Scheme of the opening angle of the emission cone and the angle β between spin and magnetic dipole axes (left panel), and its view from the top (right panel). In the figure denoted: a is the linear diameter of the emission cone, and d is the distance. The rest of the symbols are identical to the ones in Figure 1 and Figure 2. Bottom panel shows the pulse profile with period ( P s ) and width of individual pulse at 10 % of maximal intensity ( w 10 ).
Universe 09 00334 g003
Figure 4. Two data sets: P s B ns with (blue dots, calculated) and without (grey dots, retrieved from ATNF) taking into account β parameter for radio pulsar population. The dotted vertical line corresponds to the borderline of P s = 0.01 s for the millisecond pulsar population. The dashed line corresponds to P s = 0.1 s borderline for the short-periodic pulsar population. The dash-dotted line corresponds to P s = 2 s borderline separating long-periodic pulsars. Blue dots are limited to P s < 2 s since we consider only rotation-powered pulsars with MDR mechanism of their energy loss (see Section 1). A horizontal solid red line corresponds to the quantum critical threshold B cr 4.4× 10 13 G (see Section 3).
Figure 4. Two data sets: P s B ns with (blue dots, calculated) and without (grey dots, retrieved from ATNF) taking into account β parameter for radio pulsar population. The dotted vertical line corresponds to the borderline of P s = 0.01 s for the millisecond pulsar population. The dashed line corresponds to P s = 0.1 s borderline for the short-periodic pulsar population. The dash-dotted line corresponds to P s = 2 s borderline separating long-periodic pulsars. Blue dots are limited to P s < 2 s since we consider only rotation-powered pulsars with MDR mechanism of their energy loss (see Section 1). A horizontal solid red line corresponds to the quantum critical threshold B cr 4.4× 10 13 G (see Section 3).
Universe 09 00334 g004
Figure 5. Distribution of β -parameter depending on P s for radio pulsar population with P s < 2 s . The solid black line corresponds to β = f ( P s ) with condition w 10 = P s / 2 . Other lines are identical to those in Figure 4. For all objects in the ATNF pulsar catalog, their observable w 10 -parameters do not exceed half of the spin period, i.e., w 10 P s / 2 .
Figure 5. Distribution of β -parameter depending on P s for radio pulsar population with P s < 2 s . The solid black line corresponds to β = f ( P s ) with condition w 10 = P s / 2 . Other lines are identical to those in Figure 4. For all objects in the ATNF pulsar catalog, their observable w 10 -parameters do not exceed half of the spin period, i.e., w 10 P s / 2 .
Universe 09 00334 g005
Figure 6. A comparison of β -parameter estimation for a radio pulsar population obtained by Ken’ko et al. (2023) [5] (green dots, 1242 objects, left panel) and Nikitina et al. (2017) [22] (light blue dots, 246 objects, right panel), and by our method (red dots, both panels). There are negative and positive trends between our data and data by Ken’ko et al. (2023) (left panel) correspondingly, with a vertical dotted line approximately marking the spin period where two methods give similar β estimations. No trends are seen between our data and the data by Nikitina et al. (2017) (right panel).
Figure 6. A comparison of β -parameter estimation for a radio pulsar population obtained by Ken’ko et al. (2023) [5] (green dots, 1242 objects, left panel) and Nikitina et al. (2017) [22] (light blue dots, 246 objects, right panel), and by our method (red dots, both panels). There are negative and positive trends between our data and data by Ken’ko et al. (2023) (left panel) correspondingly, with a vertical dotted line approximately marking the spin period where two methods give similar β estimations. No trends are seen between our data and the data by Nikitina et al. (2017) (right panel).
Universe 09 00334 g006
Table 1. Sample of radio pulsars with estimated β and B ns ( β ) . Other parameters (coordinates, P s , P ˙ s , w 10 ) are extracted from ATNF pulsar catalog [17]. The full table containing data on 1468 objects is available online in a machine-readable format (Supplementary Table S1).
Table 1. Sample of radio pulsars with estimated β and B ns ( β ) . Other parameters (coordinates, P s , P ˙ s , w 10 ) are extracted from ATNF pulsar catalog [17]. The full table containing data on 1468 objects is available online in a machine-readable format (Supplementary Table S1).
No.NameRADEC P s P ˙ s w 10 β B ns ( β )
PSRJ2000J2000(s)(s/s)(ms)(deg)( × 10 12 G)
1J0006+183400:06:04.8+18:34:59.00.692.10e-15195.01.2954.3
2B0011+4700:14:17.7+47:46:33.41.245.64e-16142.52.1123
3J0026+632000:26:50.5+63:20:00.80.321.51e-1648.03.223.94
4B0031-0700:34:08.8−07:21:53.40.944.08e-16120.02.2016.4
5J0038-250100:38:10.2−25:01:30.70.267.60e-1915.09.010.093
...........................
Table 2. Statistics for the β -parameter, an angle between spin and magnetic dipole axes, for different groups of radio pulsars.
Table 2. Statistics for the β -parameter, an angle between spin and magnetic dipole axes, for different groups of radio pulsars.
Puls. Group β Min β Max β Mean β Median σ Total
P s < 0.01 9.4168.7825.1520.9213.9394
0.01 < P s < 0.1 3.2261.8413.339.4610.5173
0.1 < P s < 2 0.6147.517.836.984.991301
All0.6168.789.217.437.651468
Table 3. Statistics for the surface magnetic field strength B ns ( β ) taking into account the β -parameter for different groups of radio pulsars.
Table 3. Statistics for the surface magnetic field strength B ns ( β ) taking into account the β -parameter for different groups of radio pulsars.
Puls. Group B ns min ( β ) B ns max ( β ) B ns mean ( β ) B ns med ( β ) σ Total
P s < 0.01 7.76 × 1071.19 × 10101.92 × 1095.75 × 1081.68 × 10994
0.01 < P s < 0.1 1.38 × 1097.62 × 10135.28 × 10122.96 × 10101.18 × 101373
0.1 < P s < 2 2.19 × 10107.56 × 10141.69 × 10138.33 × 10123.39 × 10131301
All7.76 × 1077.56 × 10141.53 × 10137.1 × 10123.24 × 10131468
Table 4. Statistics for the comparison of β obtained from geometric approaches from Ken’ko et al. (2023) [5] and Nikitina et al. (2017) [22] with our data. All pulsars have spin periods between 0.1 s and 2 s.
Table 4. Statistics for the comparison of β obtained from geometric approaches from Ken’ko et al. (2023) [5] and Nikitina et al. (2017) [22] with our data. All pulsars have spin periods between 0.1 s and 2 s.
Pulsar SampleSelection CriteriaPulsar Number | Δ β | Min | Δ β | Max | Δ β | Mean | Δ β | Median σ
Ken’ko et al. (2023) [5] P s , w 10 1242 0.007338 61.145636 5.825957 4.61572 5.013604
Nikitina et al. (2017) [22] P s , P ˙ s , w 10 246 0.037164 40.488991 5.544564 4.094941 5.092834
Table 5. Comparison of the obtained data from this work ( β geom ) with data obtained by polarization method ( β pol ) from Nikitina et al. (2017) [22].
Table 5. Comparison of the obtained data from this work ( β geom ) with data obtained by polarization method ( β pol ) from Nikitina et al. (2017) [22].
No.Name β geom β pol | β pol β geom |
PSR(deg)(deg)(deg)
1J0108-14313.93117.07
2B0656+145.181711.82
3J0905-512714.30227.70
4J1015-57192.2552.75
5B1055-526.62158.38
6J1349-61309.056454.95
7J1355-59256.49103.51
8B1509-583.07106.93
9B1702-1910.414938.59
10J1702-43105.73115.27
11J1723-365917.962810.04
12B1800-212.64129.36
13B1822-145.0782.93
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kim, V.; Umirbayeva, A.; Aimuratov, Y. Estimates of the Surface Magnetic Field Strength of Radio Pulsars. Universe 2023, 9, 334. https://doi.org/10.3390/universe9070334

AMA Style

Kim V, Umirbayeva A, Aimuratov Y. Estimates of the Surface Magnetic Field Strength of Radio Pulsars. Universe. 2023; 9(7):334. https://doi.org/10.3390/universe9070334

Chicago/Turabian Style

Kim, Vitaliy, Adel Umirbayeva, and Yerlan Aimuratov. 2023. "Estimates of the Surface Magnetic Field Strength of Radio Pulsars" Universe 9, no. 7: 334. https://doi.org/10.3390/universe9070334

APA Style

Kim, V., Umirbayeva, A., & Aimuratov, Y. (2023). Estimates of the Surface Magnetic Field Strength of Radio Pulsars. Universe, 9(7), 334. https://doi.org/10.3390/universe9070334

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop