Eccentricities of millisecond pulsars with intermediate-mass progenitors

This paper presents an analytical model demonstrating that millisecond pulsars with intermediate-mass progenitors ($3\text{--}5\,{\rm M}_\odot$) form via stable Roche-lobe overflow, resulting in orbital eccentricities comparable to those of lower-mass systems due to a weak dependence on envelope mass, while suggesting that more massive white dwarfs likely originate from a distinct unstable channel.

The Gist

The Big Picture: Cosmic Dance Partners

Imagine the universe as a giant ballroom. In this ballroom, there are two very special dancers:

1. The Millisecond Pulsar (MSP): A super-dense, spinning neutron star that acts like a cosmic metronome, ticking with perfect precision.
2. The White Dwarf (WD): The burnt-out, cooling core of a dead star.

Usually, these two dance in a perfect circle. But sometimes, the dance is slightly elliptical (oval-shaped). The paper asks a simple question: Why do some of these pairs have oval orbits, and how did they get that way?

Specifically, the authors are trying to figure out how Millisecond Pulsars end up dancing with Carbon-Oxygen White Dwarfs (the heavier, "meatier" type of dead star) without their orbits becoming a mess.

The Old Story vs. The New Discovery

For a long time, astronomers knew how Pulsars got partners that were lightweight (Helium White Dwarfs).

• The Old Story: A low-mass star (like our Sun) swells up, spills its outer layers onto the Pulsar, and then shrinks away. The "spilling" process is smooth and stable. As the star shrinks, the Pulsar's gravity acts like a brake, smoothing out the orbit into a perfect circle. However, tiny, random bumps in the star's surface (like waves on a pond) leave a tiny bit of "wiggle" in the orbit. This explains the slight oval shape we see in lightweight pairs.

The Mystery:
Astronomers noticed that heavy White Dwarfs (Carbon-Oxygen) also dance with Pulsars, and surprisingly, they have the same amount of "wiggle" (eccentricity) as the lightweight ones.

• The Problem: Heavy stars usually behave differently. When they try to spill their layers, it's often a violent explosion (like a car crash) rather than a smooth pour. This should create a very messy, highly oval orbit. So, why are these heavy pairs so calm?

The Solution: The "Goldilocks" Zone

The authors (Hagai Bareli and Sivan Ginzburg) propose a new formation story for these heavy pairs. They found a "Goldilocks" scenario where things are just right.

The Analogy: The Overfilled Bathtub
Imagine a star as a bathtub filling with water.

• Low-mass stars are small tubs. When they overflow, the water spills slowly and gently.
• Very massive stars are huge tubs. If they overflow, the water rushes out so fast it causes a flood (a "Common Envelope" phase), crashing the two stars together.
• Intermediate-mass stars (The Goldilocks): These are medium-sized tubs. The authors suggest that if these stars start overflowing early in their life (before they get too big and convective), they can spill their water stably.

The Twist:
Usually, when a star stops spilling its water (detaches from the Pulsar), it's because it has run out of fuel.

• Light stars run out of fuel when their outer skin gets too thin.
• Medium stars in this study run out of fuel when their core suddenly lights up a new engine (Helium ignition). This causes the star to shrink rapidly, pulling away from the Pulsar just in time to avoid a crash.

The "Wiggle" Explanation

So, why do these medium-heavy stars have the same "wiggle" as the light ones?

The authors used a clever mathematical trick. They realized that the "wiggle" (eccentricity) depends on how much "skin" (envelope) the star has left when it stops spilling.

• The Intuition: You might think a medium star has much more skin left than a light star, so the wiggle should be huge.
• The Reality: The math shows that the wiggle is very stubborn. It only grows very slowly as the skin gets thicker (specifically, it scales with the 6th root of the mass).
• The Result: Even though the medium star has about 10 times more "skin" left over than the light star, the resulting wiggle is only about 1.5 times bigger. Given the natural chaos of the universe, that difference is negligible. They end up with almost the exact same oval shape.

The Three Groups of Dancers

The paper sorts all these cosmic couples into three groups based on their weight and how they formed:

1. The Lightweights (Helium WDs): Formed by low-mass stars spilling gently. They have a specific, predictable oval shape.
2. The Medium-Heavies (Carbon-Oxygen WDs < 0.6 Solar Masses): Formed by the "Goldilocks" intermediate stars spilling early and smoothly. Surprise! They have the same oval shape as the lightweights. This confirms they formed the same way.
3. The Super-Heavies (Carbon-Oxygen/ONe WDs > 0.6 Solar Masses): These are the troublemakers. They likely formed through a violent crash (Common Envelope). Their orbits are much more oval, and we don't have a perfect theory for them yet.

The Takeaway

This paper is a detective story. The authors solved the mystery of why heavy dead stars orbiting Pulsars look so calm. They proved that these stars didn't crash into their partners; instead, they performed a delicate, stable dance early in their lives, leaving behind a "wiggle" that is surprisingly similar to their lighter cousins.

They also built a simple formula (a "recipe") that predicts exactly how long the dance takes (orbital period) based on how heavy the star is, matching real-world observations perfectly.

In short: Nature found a way for medium-sized stars to be gentle giants, creating a perfect cosmic dance that looks just like the one performed by the smallest stars.

Technical Summary

1. Problem Statement

Millisecond pulsars (MSPs) with white dwarf (WD) companions exhibit exceptionally low orbital eccentricities ($e \sim 10^{-7}$ to $10^{-3}$). The formation mechanism for MSPs with low-mass Helium (He) WDs is well-established: they form via stable Roche-lobe overflow (RLO) from low-mass ($M \lesssim 2 M_\odot$) red giant branch (RGB) progenitors. The residual eccentricity is explained by the convective fluctuation–dissipation theory (Phinney 1992), where stochastic fluctuations in the giant's convective envelope excite the orbit, balanced by tidal dissipation.

However, the formation of MSPs with Carbon-Oxygen (CO) WDs (typically $M_{wd} \gtrsim 0.5 M_\odot$) remains less understood. While some CO WDs have similar eccentricities to He WDs, previous theoretical attempts to explain them (e.g., Cohen et al. 2024) relied on Case C overflow (AGB stars burning He in a shell), which typically requires a subsequent common envelope (CE) phase to shrink the orbit. This creates a discrepancy: the observed short orbital periods of many CO WDs do not match the long periods expected from AGB RLO, and the eccentricity predictions from Case C models were inconsistent with observations (predicting $e \sim 3 \times 10^{-3}$, failing to explain the lower end of the observed range).

The central question is: Can the stable RLO channel of intermediate-mass progenitors ($3 M_\odot \lesssim M \lesssim 5 M_\odot$) explain the masses, orbital periods, and eccentricities of observed MSP-CO WD binaries without invoking a common envelope phase?

2. Methodology

The authors employ a combination of numerical stellar evolution and analytical modeling:

• Numerical Simulations: They use the MESA (Modules for Experiments in Stellar Astrophysics) code to simulate binary evolution. They model intermediate-mass progenitors ($3-5 M_\odot$) transferring mass stably to a neutron star ($M_{ns} \approx 1.8 M_\odot$) during the main sequence (Case A) or early RGB (Case B).
• Analytical Modeling: They derive a simple analytical framework to relate the final orbital period ($P$) to the WD mass ($m_{wd}$) and progenitor mass ($M$). This involves adapting the radius-mass relations for non-degenerate cores (intermediate-mass) compared to degenerate cores (low-mass).
• Eccentricity Calculation: They apply Phinney's (1992) theory to calculate the final eccentricity. The theory posits that eccentricity is set by energy equipartition between the orbital epicyclic motion and convective eddies at the moment of Roche-lobe detachment.
  • Key variable: The envelope mass ($m_e$) remaining at detachment.
  • Scaling relation derived: $e \propto P (m_e/m_{wd})^{1/6}$.

3. Key Contributions and Physical Insights

A. Distinct Detachment Mechanisms

The paper highlights a fundamental difference in how low-mass and intermediate-mass stars detach from their Roche lobes:

• Low-mass ($M \lesssim 2 M_\odot$): Detach when the envelope becomes too light to support the H-burning shell (degenerate core). The remaining envelope mass is very small ($m_e \sim 10^{-3} M_\odot$).
• Intermediate-mass ($3-5 M_\odot$): Possess non-degenerate He cores. They detach when the core reaches the Schönberg–Chandrasekhar limit and subsequently ignites Helium. This ignition causes a rapid contraction of the envelope. Consequently, the envelope mass at detachment is significantly higher ($m_e \sim 0.08 M_\odot$), roughly an order of magnitude larger than in the low-mass case.

B. Analytical Model for Orbital Periods

The authors derive a new scaling relation for the final orbital period of CO WDs formed via stable Case A/B overflow:
$$P \approx 20 \left[ \frac{R_0^3}{G(m_{wd} - m_e)} \right]^{1/2} \max \left[ 1, \left( \frac{m_{wd} - m_e}{\beta M} \right)^{3\alpha/2} \right]$$
Where $\beta$ is the Schönberg–Chandrasekhar limit and $\alpha > 4.5$ is a steep power law index for the radius expansion. This model successfully reproduces numerical tracks and explains why intermediate-mass CO WDs occupy a specific region in the $P-m_{wd}$ plane (typically $P \sim 10-100$ days).

C. The Eccentricity Result (The Core Finding)

Despite the envelope mass ($m_e$) being $\sim 10\times$ larger for intermediate-mass progenitors, the authors find that the eccentricity remains largely unchanged compared to low-mass He WDs.

• Reasoning: The eccentricity scales as $e \propto m_e^{1/6}$. Because the exponent is so small ($1/6$), a factor of 10 increase in envelope mass only increases the eccentricity by a factor of $10^{1/6} \approx 1.5$.
• Result: The predicted eccentricity for intermediate-mass CO WDs follows the same linear scaling with period as He WDs:
  $$e \approx 2 \times 10^{-5} \left( \frac{P}{10 \text{ d}} \right)$$
  This contrasts sharply with previous Case C models (Cohen et al. 2024) which predicted much higher eccentricities due to much longer initial periods ($P \sim 10^3$ days).

4. Results

• Reproduction of Observations: The stable Case A/B channel successfully reproduces the observed masses and orbital periods of MSP-CO WD binaries with $m_{wd} \lesssim 0.6 M_\odot$.
• Eccentricity Match: The theoretical prediction for eccentricity ($e \sim 10^{-5}$ for $P \sim 10$ days) aligns perfectly with the observed population of intermediate-mass CO WDs. This confirms that these systems do not require a common envelope phase to explain their current state.
• Bimodal Population: The observed MSP-WD population appears bimodal:
  1. Intermediate-mass CO WDs ($m_{wd} \lesssim 0.6 M_\odot$): Formed via stable Case A/B overflow. Their eccentricities follow the Phinney (1992) scaling.
  2. Massive CO/ONe WDs ($m_{wd} \gtrsim 0.6 M_\odot$): Likely formed via unstable Case C overflow followed by common envelope inspiral. These systems exhibit a different (and currently unexplained) eccentricity distribution.

5. Significance

• Theoretical Unification: This work firmly establishes the stable Case A/B overflow of intermediate-mass stars as the primary formation channel for the majority of observed MSP-CO WD binaries. It unifies the eccentricity explanation for both He and CO WD companions under the single framework of convective fluctuation–dissipation theory.
• Resolution of Previous Discrepancies: It resolves the tension in previous models (Cohen et al. 2024) which struggled to explain low eccentricities in short-period systems without invoking complex, uncertain common envelope physics.
• Predictive Power: The model provides a parameter-free theoretical prediction for the eccentricity-period relation ($e \propto P$) for intermediate-mass CO WDs, matching observations without free parameters.
• Future Directions: The paper identifies that the most massive WDs ($m_{wd} \gtrsim 0.6 M_\odot$) remain an outlier, likely requiring a different formation history involving unstable mass transfer and common envelope evolution, the physics of which remains a challenge for future study.