pith. sign in

arxiv: 2606.20201 · v1 · pith:GLFQSKWJnew · submitted 2026-06-18 · 🌌 astro-ph.HE

Kick bimodality of neutron stars and mode dependence of their parameters

Pith reviewed 2026-06-26 16:11 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords neutron starspulsarsnatal kicksbimodal velocity distributionmagnetic fieldsradio pulsarskick modes
0
0 comments X

The pith

Neutron stars with magnetic fields below 10^11 G show only low-velocity natal kicks.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The analysis assigns roughly 200 isolated radio pulsars with measured spin and kinematic parameters to either the low- or high-velocity mode of a bimodal natal kick distribution, finding about 30 percent in the low-velocity mode. Differences in characteristic age, distance, and radio luminosity between the modes are attributed to selection effects. In contrast, parameters less affected by selection, such as pulse width, show no mode dependence. A statistically significant difference appears in magnetic field strength: pulsars with B ≲ 10^12 G are overabundant in the low-velocity mode, and no pulsars with B ≲ 10^11 G are found in the high-velocity mode.

Core claim

Analysis of observational data and theoretical modeling favors a bimodal distribution of the natal velocity kick of neutron stars. For ~200 normal isolated radio pulsars with well-determined spin and kinematic parameters, about 30% belong to the low-velocity mode. Lower field pulsars (B ≲ 10^12 G) are overabundant among objects from the low-velocity mode in comparison to the high-velocity one. Among pulsars with low field (≲ 10^11 G), no objects from the high-velocity mode of the kick distribution are identified.

What carries the argument

Bimodal natal kick velocity distribution used to classify individual pulsars into low- or high-velocity modes from their kinematic parameters.

Load-bearing premise

The bimodal kick distribution model accurately assigns pulsars to modes without substantial contamination from measurement uncertainties or selection biases in the kinematic parameters.

What would settle it

Identification of even one pulsar with magnetic field strength below 10^11 G that belongs to the high-velocity kick mode would falsify the reported absence.

read the original abstract

Analysis of observational data and theoretical modeling favors a bimodal distribution of the natal velocity kick of neutron stars. For $\sim200$ normal isolated radio pulsars with well-determined spin and kinematic parameters, we determine if they belong to the low- or high-velocity mode of the distribution. Our results demonstrate that about $30\%$ belong to the low-velocity mode. We then analyze the differences in the properties of the two sets of pulsars. For some parameters (characteristic ages, distances, and radio luminosities), we see a clear difference between the two modes. However, for these quantities, it can be easily attributed to selection bias. For those parameters that are not a subject of strong selection, such as pulse width, we do not observe any difference. Interestingly, we detect a significant difference in the magnetic field distribution between the two modes. Lower field pulsars ($B\lesssim 10^{12}$~G) are overabundant among objects from the low-velocity mode in comparison to the high-velocity one. Among pulsars with low field ($\lesssim 10^{11}$~G), we do not identify any objects from the high-velocity mode of the kick distribution. The origin of this discrepancy is not clear, and we discuss several possibilities.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper classifies ~200 normal isolated radio pulsars with well-determined spin and kinematic parameters into low- or high-velocity modes of a bimodal natal kick distribution, finding ~30% in the low-velocity mode. It reports differences in characteristic ages, distances, and radio luminosities (attributed to selection bias), no difference in pulse width, and a significant difference in magnetic field distributions, with low-B pulsars (B ≲ 10^12 G) overabundant in the low-velocity mode and no high-velocity objects identified at B ≲ 10^11 G.

Significance. If the per-pulsar mode assignments are robust and the B-field difference survives explicit bias tests, the result would suggest a physical connection between natal kick velocity and magnetic field strength, with implications for neutron star formation channels. The use of a sizable sample with kinematic constraints is a positive feature, though the absence of detailed classification and error-propagation methods limits immediate assessment of the central observational claim.

major comments (3)
  1. [Classification and sample selection] The procedure for assigning individual pulsars to low- or high-velocity kick modes from kinematic parameters is not described with sufficient detail (including the functional form of the bimodal model, handling of measurement uncertainties, or posterior probabilities), which is load-bearing for the claim that no high-velocity objects exist at B ≲ 10^11 G and for the reported overabundance at B ≲ 10^12 G.
  2. [Magnetic field analysis] No quantitative test or modeling of selection effects specific to the magnetic field distribution is presented, in contrast to the explicit attribution of age/distance/luminosity differences to bias; this omission directly affects the significance of the reported B-mode discrepancy.
  3. [Results on mode properties] The abstract and results do not report the impact of velocity measurement errors on mode contamination or the fraction of objects with ambiguous assignments, undermining in the clean separation used for the B-field comparison.
minor comments (2)
  1. [Sample description] Clarify the exact number of pulsars in each mode and any cuts applied beyond 'well-determined parameters' to allow reproducibility.
  2. [Parameter comparisons] The statement that pulse width shows no difference would benefit from a quantitative statistic (e.g., KS-test p-value) rather than qualitative description.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful and constructive review. The comments highlight areas where additional methodological detail and discussion of limitations will improve the manuscript. We address each major comment below and will incorporate revisions as noted.

read point-by-point responses
  1. Referee: [Classification and sample selection] The procedure for assigning individual pulsars to low- or high-velocity kick modes from kinematic parameters is not described with sufficient detail (including the functional form of the bimodal model, handling of measurement uncertainties, or posterior probabilities), which is load-bearing for the claim that no high-velocity objects exist at B ≲ 10^11 G and for the reported overabundance at B ≲ 10^12 G.

    Authors: We agree that the original description of the classification procedure lacked sufficient detail. In the revised manuscript we will expand Section 2 (or equivalent) to specify the functional form of the bimodal kick distribution (two-component Maxwellian with parameters taken from the literature), the Monte Carlo approach used to propagate kinematic uncertainties into velocity posteriors, and the explicit posterior probability threshold applied for mode assignment. We will also tabulate or report the number of objects with intermediate probabilities. revision: yes

  2. Referee: [Magnetic field analysis] No quantitative test or modeling of selection effects specific to the magnetic field distribution is presented, in contrast to the explicit attribution of age/distance/luminosity differences to bias; this omission directly affects the significance of the reported B-mode discrepancy.

    Authors: The referee correctly identifies an inconsistency in our treatment of selection effects. While we attributed differences in age, distance and luminosity to observational bias, we did not perform an analogous quantitative assessment for the magnetic-field distribution. In revision we will add a dedicated paragraph discussing possible B-dependent selection (e.g., radio luminosity–B correlations and survey sensitivity) and will note that the reported absence of high-velocity objects below 10^11 G remains an observational result whose robustness requires future bias modeling or larger samples. revision: yes

  3. Referee: [Results on mode properties] The abstract and results do not report the impact of velocity measurement errors on mode contamination or the fraction of objects with ambiguous assignments, undermining in the clean separation used for the B-field comparison.

    Authors: We will revise both the abstract and the results section to include quantitative statements on mode contamination arising from velocity errors and the fraction of pulsars whose posterior probabilities fall in an ambiguous range (e.g., 0.3–0.7). These statistics were computed during the analysis but omitted from the submitted text; their inclusion will directly address the concern about the robustness of the mode separation underlying the B-field comparison. revision: yes

Circularity Check

0 steps flagged

No significant circularity in classification or comparison

full rationale

The paper adopts a pre-existing bimodal kick model (from prior observational and theoretical work) to assign ~200 pulsars to low- or high-velocity modes using their measured kinematic parameters. It then performs a post-hoc comparison of magnetic-field distributions between the resulting groups. No equation or step reduces a claimed result to a fitted input by construction, no uniqueness theorem is imported from the authors' own prior work, and no ansatz is smuggled via self-citation. The B-field difference is an independent observational finding once the kinematic classification is fixed; the analysis is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review prevents enumeration of specific free parameters or axioms; the central claim rests on the prior bimodal kick model whose parameters are not detailed here.

pith-pipeline@v0.9.1-grok · 5754 in / 1171 out tokens · 29110 ms · 2026-06-26T16:11:55.146108+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

56 extracted references · 53 canonical work pages · 27 internal anchors

  1. [1]

    The ATNF Pulsar Catalogue

    The Australia Telescope National Facility Pulsar Catalogue. , keywords =. doi:10.1086/428488 , archivePrefix =. astro-ph/0412641 , primaryClass =

  2. [2]

    On the Nature of Pulsars. III. Analysis of Observations. , year = 1970, month = jun, volume =. doi:10.1086/150487 , adsurl =

  3. [3]

    , year = 1970, month = feb, volume =

    Possible Causes of the Secular Increase in Pulsar Periods. , year = 1970, month = feb, volume =

  4. [4]

    , keywords =

    MeerKAT Pulsar Timing Array parallaxes and proper motions. , keywords =. doi:10.1093/mnras/stae016 , archivePrefix =. 2401.06963 , primaryClass =

  5. [5]

    Birth and Evolution of Isolated Radio Pulsars

    Birth and Evolution of Isolated Radio Pulsars. , keywords =. doi:10.1086/501516 , archivePrefix =. astro-ph/0512585 , primaryClass =

  6. [6]

    The Velocity Distribution of Isolated Radio Pulsars

    The Velocity Distribution of Isolated Radio Pulsars. , keywords =. doi:10.1086/338805 , archivePrefix =. astro-ph/0106159 , primaryClass =

  7. [7]

    , keywords =

    VLBA Astrometry of the Galactic Double Neutron Stars PSR J0509+3801 and PSR J1930-1852: A Preliminary Transverse Velocity Distribution of Double Neutron Stars and its Implications. , keywords =. doi:10.3847/1538-4357/ad4883 , archivePrefix =. 2405.03914 , primaryClass =

  8. [8]

    The observed velocity distribution of young pulsars

    The observed velocity distribution of young pulsars. , keywords =. doi:10.1051/0004-6361/201731518 , archivePrefix =. 1708.08281 , primaryClass =

  9. [9]

    Modified pulsar current analysis: probing magnetic field evolution

    Modified pulsar current analysis: probing magnetic field evolution. , keywords =. doi:10.1093/mnras/stu1496 , archivePrefix =. 1407.6269 , primaryClass =

  10. [10]

    Pulsar Population Synthesis with Magnetorotational Evolution: Constraining the Decay of Magnetic field

    Pulsar Population Synthesis with Magnetorotational Evolution: Constraining the Decay of the Magnetic Field. , keywords =. doi:10.3847/1538-4357/ad5af8 , archivePrefix =. 2408.14786 , primaryClass =

  11. [11]

    Analysis of complete PSR

    The observed velocity distribution of young pulsars - II. Analysis of complete PSR. , keywords =. doi:10.1093/mnras/staa958 , archivePrefix =. 2002.01367 , primaryClass =

  12. [12]

    PSR J0357+3205: the tail of the turtle

    PSR J0357+3205: The Tail of the Turtle. , keywords =. doi:10.1088/0004-637X/765/1/36 , archivePrefix =. 1212.6664 , primaryClass =

  13. [13]

    The Evolution of Compact Binary Star Systems

    The Evolution of Compact Binary Star Systems. Living Reviews in Relativity , keywords =. doi:10.12942/lrr-2014-3 , archivePrefix =. 1403.4754 , primaryClass =

  14. [14]

    Can accreting isolated neutron stars be detected?

    Can accreting isolated neutron stars be detected?. arXiv e-prints , keywords =. doi:10.48550/arXiv.2512.10666 , archivePrefix =. 2512.10666 , primaryClass =

  15. [15]

    Galactic neutron stars. I. Space and velocity distributions in the disk and in the halo. , keywords =. doi:10.1051/0004-6361/200912222 , archivePrefix =. 0908.3182 , primaryClass =

  16. [16]

    , keywords =

    Interplay between neutrino kicks and hydrodynamic kicks of neutron stars and black holes. , keywords =. doi:10.1007/s10509-024-04343-1 , archivePrefix =. 2401.13817 , primaryClass =

  17. [17]

    Living Reviews in Computational Astrophysics , keywords =

    Hydrodynamics of core-collapse supernovae and their progenitors. Living Reviews in Computational Astrophysics , keywords =. doi:10.1007/s41115-020-0008-5 , archivePrefix =. 2006.05083 , primaryClass =

  18. [18]

    , keywords =

    Physical Correlations and Predictions Emerging from Modern Core-collapse Supernova Theory. , keywords =. doi:10.3847/2041-8213/ad319e , archivePrefix =. 2401.06840 , primaryClass =

  19. [19]

    , keywords =

    Combined analysis of neutron star natal kicks using proper motions and parallax measurements for radio pulsars and Be X-ray binaries. , keywords =. doi:10.1093/mnras/stab2734 , archivePrefix =. 2109.10362 , primaryClass =

  20. [20]

    , keywords =

    Natal kicks of compact objects. , keywords =. doi:10.1016/j.newar.2025.101734 , archivePrefix =. 2509.01430 , primaryClass =

  21. [21]

    , keywords =

    The Kick Velocity Distribution of Isolated Neutron Stars. , keywords =. doi:10.3847/2041-8213/adf286 , archivePrefix =. 2505.22102 , primaryClass =

  22. [22]

    Systemic velocities and merger locations of binary neutron stars

    The Galactic neutron star population - II. Systemic velocities and merger locations of binary neutron stars. , keywords =. doi:10.1093/mnras/stad3259 , archivePrefix =. 2310.14773 , primaryClass =

  23. [23]

    The mass distribution and gravitational potential of the Milky Way

    The mass distribution and gravitational potential of the Milky Way. , keywords =. doi:10.1093/mnras/stw2759 , archivePrefix =. 1608.00971 , primaryClass =

  24. [24]

    A New Electron Density Model for Estimation of Pulsar and FRB Distances

    A New Electron-density Model for Estimation of Pulsar and FRB Distances. , keywords =. doi:10.3847/1538-4357/835/1/29 , archivePrefix =. 1610.09448 , primaryClass =

  25. [25]

    Detection of the gravitational redshift in the orbit of the star S2 near the Galactic centre massive black hole

    Detection of the gravitational redshift in the orbit of the star S2 near the Galactic centre massive black hole. , keywords =. doi:10.1051/0004-6361/201833718 , archivePrefix =. 1807.09409 , primaryClass =

  26. [26]

    Vertical waves in the solar neighbourhood in Gaia DR2

    Vertical waves in the solar neighbourhood in Gaia DR2. , keywords =. doi:10.1093/mnras/sty2813 , archivePrefix =. 1809.03507 , primaryClass =

  27. [27]

    arXiv e-prints , keywords =

    Constraining the Pulsar 3D Velocity Distribution: The Impact of Spin-Velocity Alignment. arXiv e-prints , keywords =

  28. [28]

    Ages of radio pulsar: long-term magnetic field evolution

    Ages of radio pulsar: long-term magnetic field evolution. , keywords =. doi:10.1093/mnras/sty2945 , archivePrefix =. 1810.12922 , primaryClass =

  29. [29]

    Pulsar Spin-Velocity Alignment: Kinematic Ages, Birth Periods and Braking Indices

    Pulsar spin-velocity alignment: kinematic ages, birth periods and braking indices. , keywords =. doi:10.1093/mnras/stt047 , archivePrefix =. 1301.1265 , primaryClass =

  30. [30]

    A New HST Measurement of the Crab Pulsar Proper Motion

    A new HST measurement of the Crab Pulsar proper motion. , keywords =. doi:10.48550/arXiv.astro-ph/9901031 , archivePrefix =. astro-ph/9901031 , primaryClass =

  31. [31]

    The Vela Pulsar and its Synchrotron Nebula

    Vela Pulsar and Its Synchrotron Nebula. , keywords =. doi:10.1086/321533 , archivePrefix =. astro-ph/0007310 , primaryClass =

  32. [32]

    and Seljak, Uroš and Makarov, Alexey and Bahcall, Neta A

    Evidence for alignment of the rotation and velocity vectors in pulsars - II. Further data and emission heights. , keywords =. doi:10.1111/j.1365-2966.2007.12352.x , archivePrefix =. 0708.4251 , primaryClass =

  33. [33]

    Nature Astronomy , keywords =

    Evidence for three-dimensional spin-velocity alignment in a pulsar. Nature Astronomy , keywords =. doi:10.1038/s41550-021-01360-w , archivePrefix =. 2103.01839 , primaryClass =

  34. [34]

    , keywords =

    The Impact of Spin-kick Alignment on the Inferred Velocity Distribution of Isolated Pulsars. , keywords =. doi:10.3847/1538-4357/acb3c3 , archivePrefix =. 2210.12305 , primaryClass =

  35. [35]

    , keywords =

    Evidence for the spin-kick alignment of pulsars from the statistics of their magnetic inclinations. , keywords =. doi:10.1017/pasa.2025.10069 , archivePrefix =. 2412.12017 , primaryClass =

  36. [36]

    Initial spin periods of neutron stars in supernova remnants

    Initial spin periods of neutron stars in supernova remnants. , keywords =. doi:10.1007/s10509-012-1100-z , archivePrefix =. 1204.0632 , primaryClass =

  37. [37]

    Neutron stars initial spin period distribution

    Neutron star's initial spin period distribution. , keywords =. doi:10.1093/mnras/stt519 , archivePrefix =. 1303.5258 , primaryClass =

  38. [38]

    Post fall-back evolution of multipolar magnetic fields and radio pulsar activation

    Post-fall-back evolution of multipolar magnetic fields and radio pulsar activation. , keywords =. doi:10.1093/mnras/stw1902 , archivePrefix =. 1608.08806 , primaryClass =

  39. [39]

    Magnetar counterparts, kinematics and birth sites with HST and JWST

    Magnetar counterparts, kinematics and birth sites with HST and JWST. arXiv e-prints , keywords =. doi:10.48550/arXiv.2603.15750 , archivePrefix =. 2603.15750 , primaryClass =

  40. [40]

    , keywords =

    Neutrino Rocket Jet Model: An Explanation of High-velocity Pulsars and Their Spin-down Evolution. , keywords =. doi:10.3847/1538-4357/ac6cdd , archivePrefix =. 2206.02569 , primaryClass =

  41. [41]

    , keywords =

    On Neutron Star Natal Kicks in High-mass X-Ray Binaries: Insights from Population Synthesis. , keywords =. doi:10.3847/1538-4357/adc56f , archivePrefix =. 2504.19672 , primaryClass =

  42. [42]

    , keywords =

    A kinematically constrained kick distribution for isolated neutron stars. , keywords =. doi:10.1051/0004-6361/202554349 , archivePrefix =. 2503.01429 , primaryClass =

  43. [43]

    Constraining the Braking Indices of Magnetars

    Constraining the braking indices of magnetars. , keywords =. doi:10.1093/mnras/stv2465 , archivePrefix =. 1505.07013 , primaryClass =

  44. [44]

    The Dipole Magnetic Field and Spin-down Evolutions of The High Braking Index Pulsar PSR J1640-4631

    The Dipole Magnetic Field and Spin-down Evolutions of the High Braking Index Pulsar PSR J1640-4631. , keywords =. doi:10.3847/1538-4357/aa8f49 , archivePrefix =. 1709.03459 , primaryClass =

  45. [45]

    Universe , keywords =

    Estimation of Electrical Conductivity and Magnetization Parameter of Neutron Star Crusts and Applied to the High-Braking-Index Pulsar PSR J1640-4631. Universe , keywords =. doi:10.3390/universe6050063 , adsurl =

  46. [46]

    , keywords =

    The rocket effect mechanism in neutron stars in supernova remnants. , keywords =. doi:10.1093/mnras/stad1344 , archivePrefix =. 2305.00025 , primaryClass =

  47. [47]

    , keywords =

    Three dimensional magnetorotational core-collapse supernova explosions of a 39 solar mass progenitor star. , keywords =. doi:10.1093/mnras/stad1292 , archivePrefix =. 2212.00200 , primaryClass =

  48. [48]

    , keywords =

    Magnetorotational neutron star kicks. , keywords =. doi:10.1103/PhysRevD.110.083025 , archivePrefix =. 2410.09521 , primaryClass =

  49. [49]

    , keywords =

    The Origin of Pulsar Velocities and the Velocity--Magnetic Moment Correlation. , keywords =. doi:10.1086/167647 , adsurl =

  50. [50]

    The observational evidence pertinent to possible kick mechanisms in neutron stars

    The observational evidence pertinent to possible kick mechanisms in neutron stars. , keywords =. doi:10.48550/arXiv.astro-ph/9910103 , archivePrefix =. astro-ph/9910103 , primaryClass =

  51. [51]

    , keywords =

    Possible Early Linear Acceleration of Proto-neutron Stars via Asymmetric Neutrino Emission in Core-collapse Supernovae. , keywords =. doi:10.3847/2041-8213/ab30ca , archivePrefix =. 1907.04863 , primaryClass =

  52. [52]

    Core-Collapse Supernovae and Neutron Star Kicks

    Core-collapse Supernovae and Neutron Star Kicks. Radio Pulsars , year = 2003, editor =. doi:10.48550/arXiv.astro-ph/0212140 , archivePrefix =. astro-ph/0212140 , primaryClass =

  53. [53]

    arXiv e-prints , keywords =

    A Theory for Neutron Star and Black Hole Kicks and Induced Spins. arXiv e-prints , keywords =. doi:10.48550/arXiv.2311.12109 , archivePrefix =. 2311.12109 , primaryClass =

  54. [54]

    Hydrodynamical Neutron Star Kicks in Three Dimensions

    Hydrodynamical Neutron Star Kicks in Three Dimensions. , keywords =. doi:10.1088/2041-8205/725/1/L106 , archivePrefix =. 1010.0167 , primaryClass =

  55. [55]

    , keywords =

    Evolutionary Deformation toward the Elastic Limit by a Magnetic Field Confined in the Neutron-star Crust. , keywords =. doi:10.3847/1538-4357/acc06b , archivePrefix =. 2303.02312 , primaryClass =

  56. [56]

    arXiv e-prints , keywords =

    Effects of Rotation on the Gravitational Tug-Boat Mechanism for Neutron-Star Kicks and Implications for Spin-Kick Alignment. arXiv e-prints , keywords =