Kick bimodality of neutron stars and mode dependence of their parameters
Pith reviewed 2026-06-26 16:11 UTC · model grok-4.3
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.
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.
Referee Report
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)
- [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.
- [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.
- [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)
- [Sample description] Clarify the exact number of pulsars in each mode and any cuts applied beyond 'well-determined parameters' to allow reproducibility.
- [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
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
-
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
-
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
-
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
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
Reference graph
Works this paper leans on
-
[1]
The Australia Telescope National Facility Pulsar Catalogue. , keywords =. doi:10.1086/428488 , archivePrefix =. astro-ph/0412641 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1086/428488
-
[2]
On the Nature of Pulsars. III. Analysis of Observations. , year = 1970, month = jun, volume =. doi:10.1086/150487 , adsurl =
-
[3]
, year = 1970, month = feb, volume =
Possible Causes of the Secular Increase in Pulsar Periods. , year = 1970, month = feb, volume =
1970
-
[4]
MeerKAT Pulsar Timing Array parallaxes and proper motions. , keywords =. doi:10.1093/mnras/stae016 , archivePrefix =. 2401.06963 , primaryClass =
-
[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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1086/501516
-
[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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1086/338805
-
[7]
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]
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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1051/0004-6361/201731518
-
[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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1093/mnras/stu1496
-
[10]
Pulsar Population Synthesis with Magnetorotational Evolution: Constraining the Decay of the Magnetic Field. , keywords =. doi:10.3847/1538-4357/ad5af8 , archivePrefix =. 2408.14786 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4357/ad5af8
-
[11]
The observed velocity distribution of young pulsars - II. Analysis of complete PSR. , keywords =. doi:10.1093/mnras/staa958 , archivePrefix =. 2002.01367 , primaryClass =
-
[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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/0004-637x/765/1/36
-
[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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.12942/lrr-2014-3 2014
-
[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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2512.10666
-
[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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1051/0004-6361/200912222
-
[16]
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]
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]
Physical Correlations and Predictions Emerging from Modern Core-collapse Supernova Theory. , keywords =. doi:10.3847/2041-8213/ad319e , archivePrefix =. 2401.06840 , primaryClass =
-
[19]
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]
Natal kicks of compact objects. , keywords =. doi:10.1016/j.newar.2025.101734 , archivePrefix =. 2509.01430 , primaryClass =
-
[21]
The Kick Velocity Distribution of Isolated Neutron Stars. , keywords =. doi:10.3847/2041-8213/adf286 , archivePrefix =. 2505.22102 , primaryClass =
-
[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]
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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1093/mnras/stw2759
-
[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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4357/835/1/29
-
[25]
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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1051/0004-6361/201833718
-
[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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1093/mnras/sty2813
-
[27]
arXiv e-prints , keywords =
Constraining the Pulsar 3D Velocity Distribution: The Impact of Spin-Velocity Alignment. arXiv e-prints , keywords =
-
[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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1093/mnras/sty2945
-
[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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1093/mnras/stt047
-
[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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.astro-ph/9901031
-
[31]
The Vela Pulsar and its Synchrotron Nebula
Vela Pulsar and Its Synchrotron Nebula. , keywords =. doi:10.1086/321533 , archivePrefix =. astro-ph/0007310 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1086/321533
-
[32]
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]
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]
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]
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]
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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1007/s10509-012-1100-z
-
[37]
Neutron stars initial spin period distribution
Neutron star's initial spin period distribution. , keywords =. doi:10.1093/mnras/stt519 , archivePrefix =. 1303.5258 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1093/mnras/stt519
-
[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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1093/mnras/stw1902
-
[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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2603.15750
-
[40]
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]
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]
A kinematically constrained kick distribution for isolated neutron stars. , keywords =. doi:10.1051/0004-6361/202554349 , archivePrefix =. 2503.01429 , primaryClass =
-
[43]
Constraining the Braking Indices of Magnetars
Constraining the braking indices of magnetars. , keywords =. doi:10.1093/mnras/stv2465 , archivePrefix =. 1505.07013 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1093/mnras/stv2465
-
[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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4357/aa8f49
-
[45]
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]
The rocket effect mechanism in neutron stars in supernova remnants. , keywords =. doi:10.1093/mnras/stad1344 , archivePrefix =. 2305.00025 , primaryClass =
-
[47]
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]
Magnetorotational neutron star kicks. , keywords =. doi:10.1103/PhysRevD.110.083025 , archivePrefix =. 2410.09521 , primaryClass =
-
[49]
The Origin of Pulsar Velocities and the Velocity--Magnetic Moment Correlation. , keywords =. doi:10.1086/167647 , adsurl =
-
[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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.astro-ph/9910103
-
[51]
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]
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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.astro-ph/0212140 2003
-
[53]
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]
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 =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/2041-8205/725/1/l106 2041
-
[55]
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]
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 =
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.