🔭 astrophysics

Kick bimodality of neutron stars and mode dependence of their parameters

This paper analyzes approximately 200 isolated radio pulsars to confirm a bimodal distribution of natal kick velocities, finding that while most observed differences in their properties stem from selection bias, low-velocity pulsars exhibit a significant overabundance of low magnetic fields compared to their high-velocity counterparts.

Original authors: Anton D. Lazarev, Sergei B. Popov

Published 2026-06-19

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Anton D. Lazarev, Sergei B. Popov

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant, cosmic dance floor. For a long time, astronomers have been watching a specific group of dancers: neutron stars. These are the incredibly dense, spinning remnants left behind after massive stars explode in supernovae.

One of the strangest things about these dancers is how they move. When a star explodes, the neutron star doesn't just stay put; it gets kicked off the dance floor at incredible speeds. Some jog slowly, while others sprint away at thousands of kilometers per hour.

For decades, scientists have debated: Do these kicks happen randomly, or is there a pattern? A popular theory suggests the kicks are bimodal, meaning there are two distinct "modes" or groups: a Low-Velocity Group (the joggers) and a High-Velocity Group (the sprinters).

This paper by Anton Lazarev and Sergei Popov is like a detective story. They took a sample of 202 "normal" radio pulsars (a type of neutron star) and tried to figure out which group each one belongs to. Then, they asked a crucial question: Do the joggers and the sprinters have different physical traits, or are they just the same dancers moving at different speeds?

Here is how they solved the mystery, explained simply:

1. The Time-Traveling Detective Work

To figure out how fast a neutron star was kicked, you can't just look at where it is now; you have to know where it started. But the stars don't have birth certificates.

The authors used a clever trick: Time Travel.

They took the current location and speed of each star.
They ran a computer simulation to rewind time, tracing the star's path backward through the Milky Way's gravity.
They stopped the clock at the moment the star was "born" (estimated by how fast the star is slowing down).
By comparing where the star was born to where it is now, they calculated the "kick" it received.

Because they didn't know the exact direction the star was moving toward or away from us (like trying to guess if a car is driving straight or slightly sideways from a distance), they had to run thousands of simulations for each star, guessing different angles. This gave them a probability: "This star is 80% likely to be a sprinter, 20% likely to be a jogger."

2. The Results: Who Belongs to Which Group?

After sorting the 202 stars, they found:

About 30% of the stars belong to the Low-Velocity (Jogger) Group.
About 70% belong to the High-Velocity (Sprinter) Group.

This matches previous theories, confirming that the "two-mode" kick idea is likely correct.

3. The Comparison: Are the Groups Different?

Now, the authors compared the two groups to see if they were fundamentally different or just the same things moving at different speeds. They looked at several "traits":

Distance and Age (The Selection Bias Trap):
The sprinters seemed younger and lived further away. But the authors realized this was a trick of the light (or rather, the telescope). It's hard to see slow-moving objects far away, and fast-moving objects leave the neighborhood quickly. So, the fact that sprinters looked younger and further away was mostly due to how we observe them, not because they are actually different types of stars.
Pulse Width (The Heartbeat):
Neutron stars pulse like lighthouses. The authors checked if the "width" of the flash was different between joggers and sprinters. Result: No difference. They pulse the same way. This suggests that, fundamentally, they are the same kind of object.
Magnetic Fields (The Big Discovery):
This is where it got interesting. Neutron stars have powerful magnetic fields. The authors found a clear difference:
- The Low-Velocity (Jogger) Group is full of stars with weaker magnetic fields.
- The High-Velocity (Sprinter) Group rarely has weak fields. In fact, among the stars with the weakest magnetic fields, none were found in the sprinter group.

4. Why Does This Matter?

The authors are honest: they don't know exactly why the joggers have weaker magnets.

Possibility A: It's a selection bias. Maybe weak magnets make stars harder to see if they are far away (where the sprinters live).
Possibility B: It's a physical difference. Maybe the way the supernova explosion happens determines both how hard the star is kicked and how strong its magnetic field ends up being.

If it's Possibility B, it's a huge clue. It would mean that the "kick" and the "magnetic field" are two sides of the same coin, born from the same violent explosion physics.

The Bottom Line

The paper concludes that while most differences between fast and slow neutron stars can be explained by how we look at them (observation bias), there is a real, significant difference in their magnetic fields. The slow-moving stars tend to have weaker magnetic fields.

This suggests that the "kick" a neutron star gets isn't just a random shove; it might be intimately linked to the internal physics of the star's magnetic field at the moment of its birth. It's like finding out that in a car race, the cars that accelerate slowly all happen to have smaller engines, hinting that the engine design determines both the speed and the acceleration style.

Problem Statement
Neutron stars (NSs) are observed to possess spatial velocities significantly exceeding those of their progenitor stars, a phenomenon attributed to natal kicks received during supernova explosions. While a bimodal distribution of these kick velocities (comprising low- and high-velocity modes) has been supported by theoretical modeling and population synthesis (e.g., Igoshev 2020), observational confirmation remains challenging. Deriving the natal kick distribution from observed 2D transverse velocities is non-trivial due to unknown radial velocity components, uncertainties in birthplaces, and kinematic evolution within the Galactic potential. Furthermore, existing studies have not thoroughly investigated whether specific physical parameters of radio pulsars correlate with these distinct kick velocity modes, particularly distinguishing between effects caused by selection bias and intrinsic physical differences.

Methodology
The authors analyzed a sample of 202 normal, isolated radio pulsars selected from the ATNF Pulsar Catalogue (v2.7.0). The selection criteria required robust parallax measurements or dispersion measure distances < 10 kpc, available proper motions, known spin parameters ( $P, \dot{P}$ ), and the exclusion of recycled pulsars, binaries, and exotic sources.

To classify these pulsars into low- or high-velocity modes, the authors employed a trajectory integration model:

Trajectory Integration: For each pulsar, 100 distance values were sampled (based on parallax or dispersion measure). Using Galactic coordinates and proper motions, 100 present-day positions were generated.
Velocity Reconstruction: Since radial velocities ( $v_r$ ) are unobserved, the authors assumed an isotropic velocity orientation relative to the Local Standard of Rest (LSR). They sampled the angle $\theta$ between the line of sight and the 3D velocity vector from a truncated uniform distribution (Uniform(0.05 $\pi$ , 0.95 $\pi$ )) to avoid extreme outliers, generating $10^3$ velocity vectors per distance. This resulted in $10^5$ position-velocity combinations per pulsar.
Backward Integration: Using the Galactic potential (McMillan 2017), trajectories were integrated backward in time up to $1.5 \times \tau_c$ (characteristic age). Birthplaces were identified either as intersections with the Galactic plane or at $t = -\tau_c$ if no intersection occurred.
Kick Velocity Calculation: At each potential birthplace, the kick velocity ( $v_{kick} = |\vec{v} - \vec{v}_{LSR}|$ ) was calculated.
Categorization: Using the bimodal kick distribution parameters from Igoshev (2020) (a weighted sum of two Maxwellians: $\sigma_1 = 56$ km s $^{-1}$ , $\sigma_2 = 336$ km s $^{-1}$ , with weights $w_1=0.2, w_2=0.8$ ), the authors calculated the probability of each pulsar belonging to the low- or high-velocity mode. Pulsars were categorized as "unambiguously" belonging to a mode if the probability exceeded 0.5 at all possible birthplaces; otherwise, they were categorized based on the majority of birthplaces.

Key Results

Mode Distribution: Approximately 30% (60 out of 202) of the sample belongs to the low-velocity mode, consistent with previous estimates for young pulsars. The remaining ~70% belong to the high-velocity mode.
Selection Bias Effects: As expected, high-velocity pulsars were found to be, on average, younger (shorter spin-down ages) and more distant than low-velocity pulsars. These differences were confirmed by Kolmogorov-Smirnov (K-S) tests ( $p < 0.01$ ) and are attributed to observational selection effects (difficulty measuring small velocities at large distances and high-velocity objects leaving the survey volume).
Pulse Width: No significant difference was found in pulse width ( $W_{10}$ and $W_{50}$ ) between the two modes ( $p \approx 0.5$ ), suggesting that this parameter is not subject to strong selection bias or intrinsic mode dependence.
Magnetic Field Correlation: A statistically significant difference ( $p \approx 0.02$ ) was detected in the surface magnetic field distributions. Pulsars in the low-velocity mode are overabundant with lower magnetic fields ( $B \lesssim 10^{12}$ G). Notably, no pulsars with $B \lesssim 10^{11}$ G were identified in the high-velocity mode. This discrepancy persists even when accounting for distance and age, though the authors note it may be partially influenced by magnetic field decay over time.

Significance and Claims
The paper claims to provide a detailed classification of normal radio pulsars into bimodal kick velocity groups and to analyze the resulting parameter distributions. The primary significance lies in the identification of a correlation between kick velocity modes and magnetic field strength. The authors state that while selection bias and magnetic field decay may contribute to the observed difference, the origin of the discrepancy is not fully clear. They suggest that if the bimodal kick structure arises from differences in supernova explosion physics, the magnetic field difference might share a similar origin, potentially linking NS kick mechanisms to the external dipolar magnetic field. However, the authors remain modest, acknowledging that detailed population synthesis modeling is required to confirm these origins and that selection effects cannot be entirely ruled out. The study does not propose new experimental applications but rather offers a statistical analysis of existing data to refine the understanding of NS kinematics and magnetic properties.

1. The Time-Traveling Detective Work

2. The Results: Who Belongs to Which Group?

3. The Comparison: Are the Groups Different?

4. Why Does This Matter?

The Bottom Line

More like this