Formation and Detectability of PSR J0653+0443-like Extremely Eccentric Binary Neutron Stars in the Milky Way

arXiv:2610.31583 · astro-ph.HE, astro-ph.SR · Submitted 2026-10-22 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: I'm Vera.

Jocelyn: And I'm Jocelyn, here with Subrahmanyan.

Vera: We're looking at the paper: "Formation and Detectability of PSR J0653+0443-like Extremely Eccentric Binary Neutron Stars in the Milky Way".

Jocelyn: Here's the gist of the paper. The gist Extremely eccentric BNS systems are rare among detectable Galactic BNSs, but their discovery provides unique insights into neutron star formation and evolution through compact binary population synthesis calculations incorporating self-consistent spin evolution and observational selection effects.

Vera: Model Constraint Using FAST GPPS BNS The authors constrain their compact binary population synthesis (CBPS) model by using two latest observational results: (1) the local BNS merger rate density inferred from the LVK O4 observing run, and (2) the orbital properties of Galactic BNS systems detectable by the FAST GPPS survey, including PSR J0653+0443 These observational constraints provide a direct diagnostic of the key parameters in the formation and evolution of BNS systems They utilize an updated version of CBPS code BSE to simulate the formation and evolution of BNS systems During the subsequent binary evolution, the common envelope (CE) phase plays a central role in tightening the orbit and enabling the formation of compact binaries In this work, they model the CE evolution using the widely adopted alpha-formalism (Webbink 1984) To simultaneously explain the low BNS merger rate density and the existence of short-orbital-period Galactic BNS systems, they follow a “median CE” description by introducing fHG, the fraction of Hertzsprung-gap (HG) donors that are forced to merge in the CE stage Among the models considered in this work (see Supplementary Table A2), they find that M1 with alpha = 5 and fHG = 0.8 still provides the best agreement with both observational constraints, consistent with their previous results in Chen et al.

Jocelyn: (2025) The predicted merger rate density R0 ∼ 149 Gpc−3yr−1 is consistent with current observational estimates Statistical Properties of Extremely Eccentric BNS The fraction of systems exceeding different eccentricity thresholds, i.e., e > 0.90, 0.95, and 0.99, within the synthetic FAST GPPS-detectable population is quantified Systems with e > 0.99 constitute a tiny subset of the overall population, with a fraction of 0.64%, confirming that objects such as PSR J0653+0443 occupy the extreme tail of the Galactic BNS distribution The occupation fraction fSNe (> e) of different SNe channels responsible for forming the second-born NS in BNS systems as a function of eccentricity threshold e is shown in Figure 2 A clear trend emerges that the fraction of CCSNe increases with eccentricity, while the contribution from ECSNe/USSNe decreases correspondingly This transition indicates that extremely eccentric systems like J0653+0443 are preferentially associated with CCSNe channels, likely reflecting the combined effects of their typically larger natal kicks or more substantial SN mass loss Spin Evolution and Observational Bias The authors investigate the spin evolution properties of extremely eccentric systems, finding that only ∼ 14.0% of the BNS systems with e > 0.99 have experienced stable Roche lobe overflow However, when restricting the pulsar spin period to P < 0.1s, as is the case for PSR J0653+0443, the recycled fraction increases to ∼ 40% for the same eccentricity range This naturally reconciles the recycled nature of PSR J0653+0443 with the statistical properties of the high-eccentricity population The formation time of systems with e > 0.99, defined as the elapsed time since the second supernova, is about ∼ 163+1168−156 Myr (uncertainties represent the 16%th and 84%th percentiles) Applying the same spin-period cut raises the median formation time to ∼ 1701+6452−1378 Myr, reflecting the observational bias that short-period pulsars remain detectable over longer timescales, thereby allowing older systems to enter the observable sample Future Detection Prospects for Extremely Eccentric BNS The probability of detecting at least one extremely eccentric system is given by P(N eGPPS ≥ 1) = 1 − exp(−λGPPS), as listed in the third row of Table 1 For systems with e > 0.99, this probability is P(N eGPPS ≥ 1) ≃ 26.7% Future radio facilities with substantially improved sensitivity and survey speed are expected to significantly enhance the detectability of rare BNS populations Adopting the SKA survey settings from Keane et al.

Vera: (2025) and under the same population model M1, they obtain illustrative estimates of λSKA ∼ 10, 5, and 1 for eccentricity thresholds of e > 0.90, 0.95, and 0.99, respectively In addition to radio surveys, future space-based GW detectors operating in the sub-mHz band are expected to provide an independent probe of extremely eccentric BNS populations Extremely eccentric BNS systems are rare among the GPPSdetectable population, with BNSs with orbital eccentricity e > 0.99 constituting only a small fraction (∼ 0.64%) of detectable systems in the FAST survey The formation of extremely eccentric BNSs is strongly linked to the second supernova event, as BNS systems with e ≥ 0.9 are preferentially produced through CCSNe, reflecting the increasing importance of large natal kicks in generating extreme orbital configurations The discovery of PSR J0653+0443 provides a valuable probe of the immediate post-supernova phase of binary evolution How it works The authors systematically investigate the population of extremely eccentric BNS systems as PSRJ0653+0443 (i.e., e ≳ 0.9) detectable by FAST GPPS survey, utilizing compact binary population synthesis (CBPS) with NS spin evolution They incorporate the detection selection effect of the FAST GPPS survey into the population synthesis framework This step allows the simulated population to be compared directly with the observed FAST BNS sample under realistic detection conditions They perform model constraint by exploring a grid of CBPS models spanning different values of alpha and fHG, which may strongly influence both the formation rate and orbital properties of BNS systems Appendix A: Constraint on CBPS Models The initial condition for the simulation on the formation and evolution of Galactic BNS systems are set as follows The initial mass spectrum for main-sequence (MS) binary stars are generated by assuming that, the primary mass M1 follows the initial mass function given by Kroupa (2001), i.e., p(M1) ∝ M−2.3 for M1 > 1M⊙ The common envelope (CE) stage in the binary evolution is approximated by the widely-used alpha-formalism (Webbink 1984), i.e., G(M1 − m1c)M1 / λS rL = αCE Gm1c M2 / 2a f− G M1M2 / 2ai They follow Chu et al.

Title and authors: Vera: Let’s start with the title and authors of this paper, "Formation and Detectability of PSR J0653+four hundred forty-three-like Extremely Eccentric Binary Neutron Stars in the Milky Way <ref:2610.31583#pg1,Formation and Detectability of PSR J0653+0443-like Extremely Eccentric Binary Neutron>." You can see how it immediately sets the stage by focusing on that specific high-eccentricity system.

Jocelyn: It’s interesting because it’s not just a theoretical exercise; it’s tied directly to a specific, highly eccentric object we know exists, PSR J0653+four hundred forty-three <ref:2610.31583#pg1>.

Subrahmanyan: The authors are using compact binary population synthesis with self-consistent neutron star spin evolution to model these extreme BNS systems detectable by the FAST GPPS survey.

Vera: They are performing a systematic investigation into the population of extremely eccentric BNS systems as PSR J0653+four hundred forty-three which is defined as having an eccentricity greater than zero point nine <ref:2610.31583#pg2,the population of extremely eccentric BNS systems as>.

Jocelyn: The implication there is that this specific system might be a strong constraint on the formation and evolution of BNSs, suggesting it probably formed via the isolated evolution of massive binary stars rather than in dense globular clusters or something similar.

Subrahmanyan: That’s a key point; they suggest observations of this extremely eccentric BNS system may place stringent constraints on how these systems evolve, which is vital for understanding the cosmic picture.

The paper's summary: Vera: So, what is the core summary of this work? Basically, they use the FAST GPPS survey data and the local merger rate density to constrain their compact binary population synthesis model by testing different parameter settings.

Jocelyn: They find that to simultaneously explain both the low BNS merger rate density and those short-orbital-period Galactic BNS systems, they have to use a "median CE" description involving a fraction fHG of Hertzsprung-gap donors that are forced into a common envelope phase.

Subrahmanyan: That fHG parameter is crucial because it directly impacts the orbital properties during the common envelope phase, which is where these binaries get their final tight orbits.

Vera: The paper tests different alpha values for the alpha-formalism used to model that CE evolution; they find that a specific combination, M1 with alpha equals five and fHG equals zero point eight, still gives them the best agreement with both observational constraints.

Jocelyn: So, the summary boils down to finding a single set of parameters—alpha and fHG—that works for both the merger rate density and the existence of these short-period systems.

Subrahmanyan: And what this means for us is that it provides a concrete prediction for the merger rate density, which they find to be around R zero about one hundred forty-nine per gigaparsec per year, consistent with current estimates.

The paper's improvements: Vera: The paper suggests some improvements to their existing model constraints. They look at the fraction of systems exceeding different eccentricity thresholds, like e > zero point nine zero, e > zero point nine five, and e > zero point nine nine <ref:2610.31583#pg2>.

Jocelyn: They found that systems with eccentricity greater than zero point nine nine are a tiny subset, only about a fraction of zero point six four percent of the overall population detectable by the FAST survey, which confirms that things like PSR J0653+four hundred forty-three are in the extreme tail <ref:2610.31583#pg3,like PSR J0653+0443 are>.

Subrahmanyan: They also looked at how these different SN channels, or SNe channels, contribute to forming the second neutron star in BNS systems as a function of eccentricity threshold e.

Vera: The trend they see is that the fraction of CCSNe increases with eccentricity, while the contribution from ECSNe and USSNe decreases correspondingly.

Jocelyn: That transition suggests that extremely eccentric systems like J0653+four hundred forty-three are preferentially associated with CCSNe channels, which likely reflects a combination of larger natal kicks or more substantial supernova mass loss <ref:2610.31583#pg1>.

Conclusion: Vera: So, to wrap up the paper, they’ve shown that extremely eccentric BNS systems are rare among the GPPS detectable population, with systems having orbital eccentricity greater than zero point nine nine constituting only about zero point six four percent of those detectable in the FAST survey <ref:2610.31583#pg3>.

Jocelyn: The discovery of PSR J0653+four hundred forty-three with its high orbital eccentricity around zero point nine nine two, is presented as a valuable probe into the immediate post-supernova phase of binary evolution and what that means for us <ref:2610.31583#pg1>.

Subrahmanyan: From my side, it’s clear that the formation of extremely eccentric BNSs is strongly linked to the second supernova event, because they are preferentially produced through CCSNe, which reflects how important those large natal kicks are in generating those extreme orbital configurations.

Vera: It all points to the fact that these high-eccentricity systems aren't just a statistical fluke; they have a physical mechanism tied to how the supernova happens.

Jocelyn: It’s also helpful that when they restricted the pulsar spin period to less than zero point one seconds, the recycled fraction increased to about forty percent for that eccentricity range, which helps reconcile PSR J0653+four hundred forty-three’s recycled nature with its statistical properties.

Subrahmanyan: So, this paper on "Formation and Detectability of PSR J0653+four hundred forty-three-like Extremely Eccentric Binary Neutron Stars in the Milky Way" gives us a concrete understanding of the formation physics linked to high eccentricity, while also giving us a chance to predict how rare these systems are and what kind of future surveys, like SKA, might find them <ref:2610.31583#pg1,Formation and Detectability of PSR J0653+0443-like Extremely Eccentric Binary Neutron>.

National Astronomical Observatories, Chinese Academy of Sciences · Institut für Theoretische Physik, Goethe University · School of Astronomy and Space Sciences, University of Chinese Academy of Sciences

astro-ph.HE, astro-ph.SR

Submitted: 2026-10-22

Updated: 2026-10-22

Comments: 4+7 pages, revised version submitted to MNRAS as a Letter

License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/

Importance score: 82/100

The gist: The gist Extremely eccentric BNS systems are rare among detectable Galactic BNSs, but their discovery provides unique insights into neutron star formation and evolution through compact binary

Key concepts

Compact Binary Population Synthesis (CBPS): This is a computational model used to simulate how binary systems like neutron stars form and evolve over cosmic time. The model tracks key physical processes, such as mass transfer and common envelope phases, to predict the resulting population of detectable binaries in the galaxy.
Common Envelope (CE) Phase: The CE phase is a critical stage in binary evolution where one star spirals into its companion's outer layers. This process dramatically shrinks the orbit of the pair, which is essential for creating tight, compact binaries like those observed today.
Observational Selection Effects: This refers to how the specific limitations of an observational survey (like FAST GPPS) bias which systems are detected. The model accounts for these biases by simulating only those systems that would be visible under realistic detection conditions, ensuring the results match real observations.

Summary

The gist Extremely eccentric BNS systems are rare among detectable Galactic BNSs, but their discovery provides unique insights into neutron star formation and evolution through compact binary population synthesis calculations incorporating self-consistent spin evolution and observational selection effects.

Model Constraint Using FAST GPPS BNS

The authors constrain their compact binary population synthesis (CBPS) model by using two latest observational results: (1) the local BNS merger rate density inferred from the LVK O4 observing run, and (2) the orbital properties of Galactic BNS systems detectable by the FAST GPPS survey, including PSR J0653+0443 These observational constraints provide a direct diagnostic of the key parameters in the formation and evolution of BNS systems They utilize an updated version of CBPS code BSE to simulate the formation and evolution of BNS systems During the subsequent binary evolution, the common envelope (CE) phase plays a central role in tightening the orbit and enabling the formation of compact binaries In this work, they model the CE evolution using the widely adopted alpha-formalism (Webbink 1984) To simultaneously explain the low BNS merger rate density and the existence of short-orbital-period Galactic BNS systems, they follow a “median CE” description by introducing fHG, the fraction of Hertzsprung-gap (HG) donors that are forced to merge in the CE stage Among the models considered in this work (see Supplementary Table A2), they find that M1 with alpha = 5 and fHG = 0.8 still provides the best agreement with both observational constraints, consistent with their previous results in Chen et al. (2025) The predicted merger rate density R0 ∼ 149 Gpc−3yr−1 is consistent with current observational estimates

Statistical Properties of Extremely Eccentric BNS

The fraction of systems exceeding different eccentricity thresholds, i.e., e > 0.90, 0.95, and 0.99, within the synthetic FAST GPPS-detectable population is quantified Systems with e > 0.99 constitute a tiny subset of the overall population, with a fraction of 0.64%, confirming that objects such as PSR J0653+0443 occupy the extreme tail of the Galactic BNS distribution The occupation fraction fSNe (> e) of different SNe channels responsible for forming the second-born NS in BNS systems as a function of eccentricity threshold e is shown in Figure 2 A clear trend emerges that the fraction of CCSNe increases with eccentricity, while the contribution from ECSNe/USSNe decreases correspondingly This transition indicates that extremely eccentric systems like J0653+0443 are preferentially associated with CCSNe channels, likely reflecting the combined effects of their typically larger natal kicks or more substantial SN mass loss

Spin Evolution and Observational Bias

The authors investigate the spin evolution properties of extremely eccentric systems, finding that only ∼ 14.0% of the BNS systems with e > 0.99 have experienced stable Roche lobe overflow However, when restricting the pulsar spin period to P < 0.1s, as is the case for PSR J0653+0443, the recycled fraction increases to ∼ 40% for the same eccentricity range This naturally reconciles the recycled nature of PSR J0653+0443 with the statistical properties of the high-eccentricity population The formation time of systems with e > 0.99, defined as the elapsed time since the second supernova, is about ∼ 163+1168−156 Myr (uncertainties represent the 16%th and 84%th percentiles) Applying the same spin-period cut raises the median formation time to ∼ 1701+6452−1378 Myr, reflecting the observational bias that short-period pulsars remain detectable over longer timescales, thereby allowing older systems to enter the observable sample

Future Detection Prospects for Extremely Eccentric BNS

The probability of detecting at least one extremely eccentric system is given by P(N eGPPS ≥ 1) = 1 − exp(−λGPPS), as listed in the third row of Table 1 For systems with e > 0.99, this probability is P(N eGPPS ≥ 1) ≃ 26.7% Future radio facilities with substantially improved sensitivity and survey speed are expected to significantly enhance the detectability of rare BNS populations Adopting the SKA survey settings from Keane et al. (2025) and under the same population model M1, they obtain illustrative estimates of λSKA ∼ 10, 5, and 1 for eccentricity thresholds of e > 0.90, 0.95, and 0.99, respectively In addition to radio surveys, future space-based GW detectors operating in the sub-mHz band are expected to provide an independent probe of extremely eccentric BNS populations

Conclusions and Discussions

Extremely eccentric BNS systems are rare among the GPPSdetectable population, with BNSs with orbital eccentricity e > 0.99 constituting only a small fraction (∼ 0.64%) of detectable systems in the FAST survey The formation of extremely eccentric BNSs is strongly linked to the second supernova event, as BNS systems with e ≥ 0.9 are preferentially produced through CCSNe, reflecting the increasing importance of large natal kicks in generating extreme orbital configurations The discovery of PSR J0653+0443 provides a valuable probe of the immediate post-supernova phase of binary evolution

How it works

The authors systematically investigate the population of extremely eccentric BNS systems as PSRJ0653+0443 (i.e., e ≳ 0.9) detectable by FAST GPPS survey, utilizing compact binary population synthesis (CBPS) with NS spin evolution They incorporate the detection selection effect of the FAST GPPS survey into the population synthesis framework This step allows the simulated population to be compared directly with the observed FAST BNS sample under realistic detection conditions They perform model constraint by exploring a grid of CBPS models spanning different values of alpha and fHG, which may strongly influence both the formation rate and orbital properties of BNS systems

Appendix A: Constraint on CBPS Models

The initial condition for the simulation on the formation and evolution of Galactic BNS systems are set as follows The initial mass spectrum for main-sequence (MS) binary stars are generated by assuming that, the primary mass M1 follows the initial mass function given by Kroupa (2001), i.e., p(M1) ∝ M−2.3 for M1 > 1M⊙ The common envelope (CE) stage in the binary evolution is approximated by the widely-used alpha-formalism (Webbink 1984), i.e., G(M1 − m1c)M1 / λS rL = αCE Gm1c M2 / 2a f− G M1M2 / 2ai They follow Chu et al. (2025) to simulate the spin evolution of −3 − logPorb [days] for the BNS system

Appendix B: Physical Properties of Extremely Eccentric BNS Systems

The top panel of Figure B2 shows the normalized distribution of kick velocity vk of the second SNe in the formation of BNS systems<ref:

Improvements for AI systems

Improved population synthesis models can now self-consistently incorporate neutron star spin evolution mechanisms, allowing AI to predict and model observable properties like pulsar spin periods based on binary formation channels such as core-collapse supernovae (CCSNe) versus electron-capture supernovae (ECSNe/USSNe).
The improved system can perform realistic observational selection effect modeling by explicitly incorporating the radio detectability of pulsars criterion, specifically identifying BNS systems as radio-loud if their final spin is above the death line.
The AI system can quantify the rarity of extremely eccentric BNS systems by calculating conditional probabilities, such as P(Ne GPPS >= 1 Ne GPPS >= 1), providing a statistically rigorous assessment of whether discoveries like PSR J0653+0443 are plausible within the standard binary evolution framework.
The system can predict the preferred supernova channel for forming extremely eccentric BNSs, as it shows that "systems with e > 0.9 are preferentially produced through CCSNe," linking orbital eccentricity directly to natal kick velocity distributions like those from CCSNe versus ECSNe/USSNe.
The improved AI can constrain the formation time of these extreme systems by showing how imposing a short spin-period cut (P < 0.1s) increases the median formation time to ∼ 1701+6452 Myr, naturally reconciling the recycled nature of systems like PSR J0653+0443 with their statistical properties.
The system can forecast future survey sensitivities by estimating expected detection numbers (lambda) for next-generation facilities like SKA, predicting that these surveys are expected to substantially increase these extremely eccentric systems.