Recognition: no theorem link
Evolution of Crab Pulsar: Magnetic Inclination Angle and Spin
Pith reviewed 2026-05-13 02:03 UTC · model grok-4.3
The pith
The routine pulsar evolution model with bulk viscosity best matches the Crab pulsar's magnetic inclination angle, spin period, and spin-down rate.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The routine evolution model simultaneously describes the spin-down caused by the magnetic dipole radiation (MDR) and gravitational wave radiation (GWR), damping of the free-body precession owing to the bulk viscosity, and GWR-induced quenching of the magnetic inclination angle χ. When applied to the Crab pulsar, the routine model best reproduces the magnetic inclination angle χ, the spin period P, and the spin period derivative Ṗ simultaneously, indicating the important role of bulk viscosity. The calculated magnetic inclination angle derivative χ̇ is (6.3×10^{-3} - 0.3) degree/century, in agreement with the observed tiny χ̇ = 0.62 degree/century.
What carries the argument
Bulk viscosity damping of free-body precession combined with gravitational wave radiation quenching of the magnetic inclination angle χ in the routine evolution model.
If this is right
- Shear viscosity has negligible impact on radio-pulsar evolution.
- R-modes are negligible under current observational limits on their amplitude.
- Electromagnetic torque and fallback disk accretion both suppress inclination angle growth but lead to poorer fits for the Crab pulsar.
- The routine model without these additions best explains the Crab's observed properties, underscoring bulk viscosity's role.
Where Pith is reading between the lines
- This suggests bulk viscosity may be essential for modeling the evolution of other young pulsars with similar characteristics.
- Agreement on the small inclination change rate could help constrain the internal temperature and composition of neutron stars.
- If accretion and torque are overestimated for Crab, revised models of fallback disks might be needed for young pulsars.
Load-bearing premise
The parameters of the routine model, including the strength of bulk viscosity, are correctly calibrated for the Crab pulsar without significant overfitting to its data.
What would settle it
A future measurement showing the Crab pulsar's magnetic inclination angle changing at a rate outside the range of 0.0063 to 0.3 degrees per century would contradict the model's prediction.
Figures
read the original abstract
The well-observed Crab pulsar helps one to uncover the underlying knowledge about pulsar evolution. The routine evolution model simultaneously describes the spin-down caused by the magnetic dipole radiation (MDR) and gravitational wave radiation (GWR), damping of the free-body precession owing to the bulk viscosity, and GWR-induced quenching of the magnetic inclination angle $\chi$. We explore the pulsar evolution based on this routine model supplemented with the effects of shear viscosity, r-mode, electromagnetic torque, and accretion, respectively, with the stellar thermal evolution as an important input. The impact of shear viscosity on radio-pulsar evolution is negligible, as it only slightly increases the magnetic inclination angle and promotes spin-down in magnetars. Under the observational limit for its saturation amplitude, the r-mode also turns out to be completely negligible. Yet, the electromagnetic torque (under certain conditions), along with the accretion based on our three-dimensional fallback disk accretion model, are all shown to suppress the growth of the magnetic inclination angle. When applied to the Crab pulsar, the routine model best reproduces the magnetic inclination angle $\chi$, the spin period $P$, and the spin period derivative $\Dot{P}$ simultaneously, indicating the important role of bulk viscosity. The inclusion of the electromagnetic torque and accretion works even worse, suggesting these two factors perhaps are overestimated for Crab pulsar. Intriguingly, the calculated magnetic inclination angle derivative $\Dot{\chi}$ is $(6.3\times 10^{-3} - 0.3)\, {\rm degree/century}$ with the routine model, also in agreement with the observed tiny $\Dot{\chi} = 0.62\, {\rm degree/century}$.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a 'routine' pulsar evolution model combining magnetic dipole radiation (MDR), gravitational wave radiation (GWR), bulk-viscosity damping of free-body precession, and GWR-induced quenching of the magnetic inclination angle χ. It supplements this baseline with shear viscosity, r-modes, electromagnetic torque, and a three-dimensional fallback-disk accretion model while incorporating stellar thermal evolution. When applied to the Crab pulsar, the routine model is reported to simultaneously reproduce the observed χ, spin period P, and period derivative Ṗ better than the extensions; the predicted χ̇ range (6.3×10^{-3}–0.3 deg/century) is stated to agree with the observed value of 0.62 deg/century, thereby highlighting the role of bulk viscosity.
Significance. If the two free parameters (r-mode saturation amplitude and bulk-viscosity coefficient) are fixed by independent microphysical or multi-pulsar constraints rather than optimized to the Crab data, the work would usefully isolate bulk viscosity as a dominant mechanism and supply a falsifiable prediction for χ̇. The simultaneous match to three observables is a positive feature. However, the significance is limited by the absence of explicit documentation that the viscosity strength was not adjusted post-hoc to enforce the observed χ while satisfying P and Ṗ.
major comments (2)
- [Abstract] Abstract: The central claim that the routine model 'best reproduces' Crab's χ, P, and Ṗ simultaneously (thereby demonstrating the importance of bulk viscosity) is load-bearing. The manuscript does not state whether the bulk-viscosity coefficient is determined from independent microphysical calculations or multi-pulsar fits, or whether it is varied to achieve the reported agreement. If the latter, the 'best reproduces' result and the derived χ̇ interval become consistent with the data by construction, and the statement that electromagnetic torque and accretion 'work even worse' does not isolate bulk viscosity as physically required.
- [Abstract] Abstract and § on the three-dimensional fallback disk model: The claim that accretion 'works even worse' rests on the new 3D fallback-disk model. No quantitative comparison is provided to existing 2D accretion models or to observational constraints on fallback-disk mass and viscosity for the Crab, leaving open the possibility that the suppression of χ growth is an artifact of the specific implementation rather than a robust physical effect.
minor comments (2)
- The notation for the magnetic inclination angle derivative (χ̇) and its units should be introduced consistently in the text and abstract to avoid ambiguity with the observed value.
- A table or figure summarizing the χ, P, Ṗ, and χ̇ residuals for the routine model versus each extension would improve clarity and allow direct assessment of the 'best reproduces' statement.
Simulated Author's Rebuttal
We thank the referee for the constructive comments on our manuscript. We address each of the major comments point by point below, providing clarifications and indicating where revisions will be made to strengthen the presentation.
read point-by-point responses
-
Referee: [Abstract] Abstract: The central claim that the routine model 'best reproduces' Crab's χ, P, and Ṗ simultaneously (thereby demonstrating the importance of bulk viscosity) is load-bearing. The manuscript does not state whether the bulk-viscosity coefficient is determined from independent microphysical calculations or multi-pulsar fits, or whether it is varied to achieve the reported agreement. If the latter, the 'best reproduces' result and the derived χ̇ interval become consistent with the data by construction, and the statement that electromagnetic torque and accretion 'work even worse' does not isolate bulk viscosity as physically required.
Authors: The bulk-viscosity coefficient in our routine model is taken from standard microphysical calculations for neutron star matter, as commonly used in the literature on pulsar precession damping (specific references are provided in the methods section of the manuscript). It was not varied or optimized to fit the Crab pulsar observations; the model parameters are fixed a priori, and the good match to the three observables (χ, P, Ṗ) is a result of the physics included. We agree that this should be stated more explicitly to avoid any ambiguity. In the revised manuscript, we will update the abstract and add a clarifying sentence in the model description section to explicitly note that the viscosity coefficient is determined from independent microphysical constraints rather than adjusted to the Crab data. This will also reinforce that the χ̇ range is a prediction, not a fit. revision: yes
-
Referee: [Abstract] Abstract and § on the three-dimensional fallback disk model: The claim that accretion 'works even worse' rests on the new 3D fallback-disk model. No quantitative comparison is provided to existing 2D accretion models or to observational constraints on fallback-disk mass and viscosity for the Crab, leaving open the possibility that the suppression of χ growth is an artifact of the specific implementation rather than a robust physical effect.
Authors: Our three-dimensional fallback disk accretion model is developed in this work to provide a more physically detailed treatment of the accretion process, including the geometry and dynamics not captured in simpler 2D approximations. The results show that accretion suppresses the growth of the magnetic inclination angle, leading to poorer agreement with the observed χ for the Crab compared to the routine model. We acknowledge that the manuscript does not include a direct quantitative comparison to prior 2D models or specific observational limits on the Crab's fallback disk parameters. To address this, we will revise the relevant section to include a discussion of how our 3D model differs from 2D implementations and reference available constraints from observations of young pulsars and supernova remnants to argue that the suppression effect is physically motivated rather than implementation-specific. revision: yes
Circularity Check
No significant circularity detected
full rationale
The paper defines a routine evolution model from first-principles components (MDR + GWR + bulk-viscosity damping of precession + GWR quenching of χ) with stellar thermal evolution supplied as an independent input. It then numerically integrates the coupled differential equations for P, Ṗ, and χ under this model and several extensions, comparing the outputs at the Crab's current age to the observed values. No equation or section in the supplied text shows a parameter (such as the bulk-viscosity coefficient) being adjusted to enforce agreement with Crab's χ while the same data are later presented as a 'prediction'; the reported χ̇ interval is obtained by propagating the same differential system forward rather than by algebraic rearrangement of the fit. Self-citations, if present, are not invoked to establish uniqueness or to close the derivation loop. The comparison to Crab data therefore functions as an external test rather than a self-referential identity.
Axiom & Free-Parameter Ledger
free parameters (2)
- r-mode saturation amplitude
- bulk viscosity coefficient
axioms (2)
- domain assumption Neutron star thermal evolution follows standard cooling curves that couple to viscosity and torque terms.
- domain assumption Gravitational wave radiation quenches magnetic inclination angle growth.
invented entities (1)
-
three-dimensional fallback disk accretion model
no independent evidence
Reference graph
Works this paper leans on
-
[1]
Abrikosov, A. A., & Khalatnikov, I. M. 1959, Rep. Prog. Phys., 22, 329, doi: 10.1088/0034-4885/22/1/310
-
[2]
Aloy, M. ., & Obergaulinger, M. 2020, MNRAS, 500, 4365, doi: 10.1093/mnras/staa3273
-
[3]
Alpar, M. A., & Sauls, J. A. 1988, ApJ, 327, 723, doi: 10.1086/166228
-
[4]
Anderson, R. H., Pethick, C. J., & Quader, K. F. 1987, Phys. Rev. B, 35, 1620, doi: 10.1103/PhysRevB.35.1620
-
[5]
A New class of unstable modes of rotating relativistic stars
Andersson, N. 1998, ApJ, 502, 708, doi: 10.1086/305919
-
[6]
Andersson, N., Comer, G., & Glampedakis, K. 2005, Nucl. Phys. A, 763, 212, doi: 10.1016/j.nuclphysa.2005.08.012
-
[7]
Archibald, R. F., Gotthelf, E. V., Ferdman, R. D., et al. 2016, ApJL, 819, L16, doi: 10.3847/2041-8205/819/1/L16
-
[8]
2024, JCAP, 2024, 067, doi: 10.1088/1475-7516/2024/05/067
Ardavan, H. 2024, JCAP, 2024, 067, doi: 10.1088/1475-7516/2024/05/067
-
[9]
Arons, J., & Barnard, J. J. 1986, ApJ, 302, 120, doi: 10.1086/163978
-
[10]
1993, ApJ, 403, 249, doi: 10.1086/172198
Arons, J., & Tavani, M. 1993, ApJ, 403, 249, doi: 10.1086/172198
-
[11]
Beskin, V. S., & Zheltoukhov, A. A. 2014, Physics Uspekhi, 57, 799, doi: 10.3367/UFNe.0184.201408e.0865
-
[12]
2021, MNRAS, 505, 1775, doi: 10.1093/mnras/stab1378
Biryukov, A., & Abolmasov, P. 2021, MNRAS, 505, 1775, doi: 10.1093/mnras/stab1378
-
[13]
Chen, W.-C., & Piekarewicz, J. 2015, Phys. Lett. B, 748, 284, doi: 10.1016/j.physletb.2015.07.020
-
[14]
2014, ApJL, 786, L13, doi: 10.1088/2041-8205/786/2/L13
Cheng, Q., & Yu, Y.-W. 2014, ApJL, 786, L13, doi: 10.1088/2041-8205/786/2/L13
-
[15]
Cheng, Q., Zhang, S.-N., Yu, Y.-W., & Zheng, X.-P. 2018, Phys. Rev. D, 97, 103012, doi: 10.1103/PhysRevD.97.103012
-
[16]
Chugunov, A. I. 2018, MNRAS, 482, 3045, doi: 10.1093/mnras/sty2867
-
[17]
Cutler, C. 2002, Phys. Rev. D, 66, 084025, doi: 10.1103/PhysRevD.66.084025
-
[18]
Cutler, C., & Jones, D. I. 2000, PhRvD, 63, 024002, doi: 10.1103/PhysRevD.63.024002
-
[19]
Dai, Z. G., & Liu, R.-Y. 2012, ApJ, 759, 58, doi: 10.1088/0004-637X/759/1/58 Dall’Osso, S., Shore, S. N., & Stella, L. 2009, MNRAS, 398, 1869, doi: 10.1111/j.1365-2966.2008.14054.x de Araujo, J. C., Coelho, J. G., & Costa, C. A. 2016, J. Cosmol. Astropart. Phys., 2016, 023, doi: 10.1088/1475-7516/2016/07/023
-
[20]
Dessart, L., Burrows, A., Livne, E., & Ott, C. D. 2007, ApJ, 669, 585, doi: 10.1086/521701
-
[21]
Dessart, L., Burrows, A., Ott, C. D., et al. 2006, ApJ, 644, 1063, doi: 10.1086/503626 D´ ıaz Teodori, M. A., Kajava, J. J. E., S´ anchez-Fern´ andez, C., et al. 2025, A&A, 695, A44, doi: 10.1051/0004-6361/202452243
-
[22]
Dong, J. M. 2021, MNRAS, 500, 1505, doi: 10.1093/mnras/staa3258
-
[23]
Dong, J. M., Lombardo, U., Zhang, H. F., & Zuo, W. 2016, ApJ, 817, 6, doi: 10.3847/0004-637X/817/1/6 Ek¸ si, K. Y. 2017, MNRAS, 469, 1974, doi: 10.1093/mnras/stx1035
-
[24]
Eki, K. Y., Anda, I. C., kntolu, S., et al. 2016, ApJ, 823, 34, doi: 10.3847/0004-637X/823/1/34
-
[25]
Fattoyev, F. J., Horowitz, C. J., Piekarewicz, J., & Reed, B. 2020, Phys. Rev. C, 102, 065805, doi: 10.1103/PhysRevC.102.065805
-
[26]
Faucher-Gigure, C.-A., & Kaspi, V. M. 2006, ApJ, 643, 332, doi: 10.1086/501516
-
[27]
Fonseca, E., Cromartie, H. T., Pennucci, T. T., et al. 2021, ApJL, 915, L12, doi: 10.3847/2041-8213/ac03b8
-
[28]
2013, ApJ, 775, 124, doi: 10.1088/0004-637X/775/2/124
Fu, L., & Li, X.-D. 2013, ApJ, 775, 124, doi: 10.1088/0004-637X/775/2/124
-
[29]
Gao, Z. F., Li, X. D., Wang, N., et al. 2016, MNRAS, 456, 55, doi: 10.1093/mnras/stv2465
-
[30]
2017, ApJ, 849, 19, doi: 10.3847/1538-4357/aa8f49
Gao, Z.-F., Wang, N., Shan, H., Li, X.-D., & Wang, W. 2017, ApJ, 849, 19, doi: 10.3847/1538-4357/aa8f49
-
[31]
Ghosh, P., & Lamb, F. K. 1979, ApJ, 232, 259, doi: 10.1086/157285
-
[32]
Ghosh, P. and Lamb, F. K. 1979, ApJ, 234, 296, doi: 10.1086/157498
-
[33]
1970, ApJL, 160, L11, doi: 10.1086/180513 14Liu et al
Goldreich, P. 1970, ApJL, 160, L11, doi: 10.1086/180513 14Liu et al
-
[34]
Haensel, P., Levenfish, K. P., & Yakovlev, D. G. 2000, A&A, 357, 1157, doi: 10.48550/arXiv.astro-ph/0004183
-
[35]
Harding, A. K., Stern, J. V., Dyks, J., & Frackowiak, M. 2008, ApJ, 680, 1378, doi: 10.1086/588037
-
[36]
Haskell, B., & Patruno, A. 2017, Phys. Rev. Lett., 119, 161103, doi: 10.1103/PhysRevLett.119.161103
-
[37]
Ho, W. C. G. 2015, MNRAS, 452, 845, doi: 10.1093/mnras/stv1339
- [38]
-
[39]
Jones, P. B. 1988, MNRAS, 233, 875, doi: 10.1093/mnras/233.4.875
-
[40]
Kou, F. F., & Tong, H. 2015, MNRAS, 450, 1990, doi: 10.1093/mnras/stv734
-
[41]
2014, EPJ Web Conf., 64, 01001, doi: 10.1051/epjconf/20136401001
Lai, Dong. 2014, EPJ Web Conf., 64, 01001, doi: 10.1051/epjconf/20136401001
-
[42]
Lander, S. K., & Jones, D. I. 2017, MNRAS, 467, 4343, doi: 10.1093/mnras/stx349
-
[43]
2018, MNRAS, 481, 4169, doi: https://doi.org/10.1093/mnras/sty2553
Lander, S K and Jones, D I. 2018, MNRAS, 481, 4169, doi: https://doi.org/10.1093/mnras/sty2553
-
[44]
2017, ApJ, 849, 47, doi: 10.3847/1538-4357/aa9074
Lei, W.-H., Zhang, B., Wu, X.-F., & Liang, E.-W. 2017, ApJ, 849, 47, doi: 10.3847/1538-4357/aa9074
-
[45]
2021, ApJ, 907, 87, doi: 10.3847/1538-4357/abcc70
Li, S.-Z., Yu, Y.-W., Gao, H., & Zhang, B. 2021, ApJ, 907, 87, doi: 10.3847/1538-4357/abcc70
-
[46]
Liu, B. S., & Li, X.-D. 2015, ApJ, 814, 75, doi: 10.1088/0004-637X/814/1/75
-
[47]
Liu, C. X., Feng, T. F., & Dong, J. M. 2024, MNRAS, 534, 1763, doi: 10.1093/mnras/stae2170
-
[48]
Lovelace, R. V. E., Romanova, M. M., & Bisnovatyi-Kogan, G. S. 1995, MNRAS, 275, 244, doi: 10.1093/mnras/275.2.244
-
[49]
Lovelace, R. V. E., Romanova, M. M., & Bisnovatyi-Kogan, G. S. 1999, ApJ, 514, 368, doi: 10.1086/306945
-
[50]
Lovelace, R. V. E., Sutton, J. M., & Craft, H. D. 1968, IAUC, 2113, 1
work page 1968
-
[51]
2013, Science, 342, 598, doi: 10.1126/science.1243254
Lyne, A., Graham-Smith, F., Weltevrede, P., et al. 2013, Science, 342, 598, doi: 10.1126/science.1243254
-
[52]
Lyne, A. G., Jordan, C. A., Graham-Smith, F., et al. 2014, MNRAS, 446, 857, doi: 10.1093/mnras/stu2118
-
[53]
Lyne, A. G., Pritchard, R. S., & Graham Smith, F. 1993, MNRAS, 265, 1003, doi: 10.1093/mnras/265.4.1003
-
[54]
2001, ApJ, 554, L63, doi: 10.1086/320927
Menou, K., Perna, R., & Hernquist, L. 2001, ApJ, 554, L63, doi: 10.1086/320927
-
[55]
Mestel, L., & Takhar, H. S. 1972, MNRAS, 156, 419, doi: 10.1093/mnras/156.4.419
-
[56]
D., Beniamini, P., & Giannios, D
Metzger, B. D., Beniamini, P., & Giannios, D. 2018, ApJ, 857, 95, doi: 10.3847/1538-4357/aab70c
-
[57]
2008, MNRAS, 390, 769, doi: 10.1111/j.1365-2966.2008.13778.x
Metzger, B. D., Quataert, E., & Thompson, T. A. 2008, MNRAS, 385, 1455, doi: 10.1111/j.1365-2966.2008.12923.x
- [58]
-
[59]
J., Lindblom, L., Cutler, C., et al
Owen, B. J., Lindblom, L., Cutler, C., et al. 1998, Phys. Rev. D, 58, 084020, doi: 10.1103/PhysRevD.58.084020
-
[60]
Page, D., Geppert, U., & Weber, F. 2006, Nucl. Phys. A, 777, 497, doi: https://doi.org/10.1016/j.nuclphysa.2005.09.019
-
[61]
Page, D., Lattimer, J. M., Prakash, M., & Steiner, A. W. 2004, Astrophys.J.Suppl.Ser., 155, 623, doi: 10.1086/424844
-
[62]
Philippov, A., Tchekhovskoy, A., & Li, J. G. 2014, MNRAS, 441, 1879, doi: 10.1093/mnras/stu591
-
[63]
Piro, A. L., & Ott, C. D. 2011, ApJ, 736, 108, doi: 10.1088/0004-637X/736/2/108
-
[64]
Raithel, C. A., ¨Ozel, F., & Psaltis, D. 2016, Phys. Rev. C, 93, 032801, doi: 10.1103/PhysRevC.93.032801
-
[65]
Rajbhandari, B., Owen, B. J., Caride, S., & Inta, R. 2021, Phys. Rev. D, 104, 122008, doi: 10.1103/PhysRevD.104.122008 ran Yang, H., & dong Li, X. 2023, ApJ, 945, 2, doi: 10.3847/1538-4357/acba09
-
[66]
2013, ApJ, 772, 1, doi: 10.1088/0004-637X/772/1/1 S´ a, P
Raskin, C., & Kasen, D. 2013, ApJ, 772, 1, doi: 10.1088/0004-637X/772/1/1 S´ a, P. M., & Tom´ e, B. 2005, Phys. Rev. D, 71, 044007, doi: 10.1103/PhysRevD.71.044007
-
[67]
2020, A&A, 642, A42, doi: 10.1051/0004-6361/202038193
Negreiros, Rodrigo. 2020, A&A, 642, A42, doi: 10.1051/0004-6361/202038193
- [68]
-
[69]
Shang, X., Wang, P., Zuo, W., & Dong, J. 2020, Phys. Lett. B, 811, 135963, doi: https://doi.org/10.1016/j.physletb.2020.135963
-
[70]
Shaw, B., Lyne, A. G., Stappers, B. W., et al. 2018, MNRAS, 478, 3832, doi: 10.1093/mnras/sty1294
-
[71]
1989, Nature, 342, 656, doi: 10.1038/342656a0
Shibazaki, N., Murakami, T., Shaham, J., & Nomoto, K. 1989, Nature, 342, 656, doi: 10.1038/342656a0
-
[72]
Shternin, P. S., & Yakovlev, D. G. 2008, Phys. Rev. D, 78, 063006, doi: 10.1103/PhysRevD.78.063006
-
[73]
2006, ApJ, 648, L51, doi: 10.1086/507518
Spitkovsky, A. 2006, ApJ, 648, L51, doi: 10.1086/507518
-
[74]
Staelin, D. H., & Reifenstein, III, E. C. 1968, Science, 162, 1481, doi: 10.1126/science.162.3861.1481
-
[75]
Taam, R. E., & van den Heuvel, E. P. J. 1986, ApJ, 305, 235, doi: 10.1086/164243
-
[76]
Tong, H., & Kou, F. F. 2017, ApJ, 837, 117, doi: 10.3847/1538-4357/aa60c6
-
[77]
2009, ApJ, 695, 1289, doi: 10.1088/0004-637X/695/2/1289
Johnston, S. 2009, ApJ, 695, 1289, doi: 10.1088/0004-637X/695/2/1289
-
[78]
2022, MNRAS, 513, 1365, doi: 10.1093/mnras/stac859
Xie, L., Wei, D.-M., Wang, Y., Li, L., & Jin, Z.-P. 2022, MNRAS, 513, 1365, doi: 10.1093/mnras/stac859
-
[79]
Xu, R. X., & Qiao, G. J. 2001, ApJ, 561, L85, doi: 10.1086/324381 Crab pulsar15
-
[80]
2001, Phys.Rep., 354, 1, doi: https://doi.org/10.1016/S0370-1573(00)00131-9
Yakovlev, D., Kaminker, A., Gnedin, O., & Haensel, P. 2001, Phys.Rep., 354, 1, doi: https://doi.org/10.1016/S0370-1573(00)00131-9
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.