Violent mergers can explain the inflated state of some of the fastest stars in the Galaxy

Basic Information

• Paper ID: 2510.12197
• Title: Violent mergers can explain the inflated state of some of the fastest stars in the Galaxy
• Authors: Aakash Bhat, Rüdiger Pakmor, Ken J. Shen, Evan B. Bauer, Abinaya Swaruba Rajamuthukumar
• Classification: astro-ph.SR (Stellar and Stellar Physics)
• Publication Date: October 15, 2025
• Journal: Astronomy & Astrophysics
• Paper Link: https://arxiv.org/abs/2510.12197

Abstract

This paper investigates the formation mechanism of a large population of hypervelocity stars (velocities between 1500-2500 km/s) observed in the Milky Way. These stars are believed to be surviving donor stars from Type Ia supernovae; however, surprisingly, they exhibit an inflated state with radii approximately one order of magnitude larger than expected for Roche lobe-filling donor stars. Through novel Arepo 3D hydrodynamic simulations, the authors study the violent merger and explosion of two massive CO white dwarfs and track the long-term evolution of a secondary star core remnant with mass 0.16 M⊙ using the MESA stellar evolution code, successfully explaining the observed characteristics of two hypervelocity stars, D6-1 and D6-3.

Research Background and Motivation

Problem Definition

1. Core Issue: Hypervelocity stars (1500-2500 km/s) observed in the Milky Way exhibit anomalously inflated states with radii ~0.1 R⊙, far exceeding the expected white dwarf radius of <0.01 R⊙
2. Observational Challenge: The existing D6 scenario (dynamically-driven double degenerate double detonation) can only partially explain the hottest and most compact escaping stars
3. Theoretical Deficiency: Among eight known candidate objects, nearly all are in inflated states, yet existing models predict inflation timescales (~1000 years) far shorter than kinematic ages (0.5-3 Myr)

Research Significance

• Cosmological Importance: These stars provide direct evidence of Type Ia supernova explosions, crucial for understanding white dwarf binary system evolution
• Stellar Physics: Offers a unique window into stellar structure and evolution under extreme conditions
• Observational Validation: Provides theoretical guidance for searching for escaping companions in supernova remnants

Limitations of Existing Approaches

1. D6 Scenario: Assumes Roche lobe-filling donor stars but cannot explain cooler hypervelocity white dwarfs (Teff < 9000K)
2. Merger-Disruption Scenario: While capable of producing higher velocities, still cannot fully match observed cooler stars
3. Timescale Problem: Existing models' inflation timescales are inconsistent with observed kinematic ages

Core Contributions

1. First 3D Violent Merger Simulation: Proposes a "violent merger" mechanism distinct from the D6 scenario, where the secondary star undergoes carbon ignition as it plunges toward the primary on the final orbit
2. New Velocity Mechanism: Demonstrates that the orbital velocity of merging white dwarfs can reach √(1+q) times the velocity at Roche lobe overflow, explaining higher escape velocities
3. Long Timescale Evolution Model: Tracks 100 Myr of stellar evolution using MESA code, for the first time matching the observed characteristics of D6-1 and D6-3
4. Thermal Timescale Analysis: Discovers that local thermal timescale rather than Kelvin-Helmholtz timescale determines the stellar inflation age, resolving the long-standing timescale problem
5. Supernova Remnant Predictions: Predicts the possible existence of dimmer, cooler escaping stars in supernova remnants, providing guidance for future observations

Methodology Details

Task Definition

Investigate the formation mechanism and long-term evolution of hypervelocity escaping stars produced by violent mergers in white dwarf binary systems, explaining the observed inflated states and physical characteristics.

Model Architecture

3D Hydrodynamic Simulation (Arepo)

1. Initial Conditions:
   • Primary star: 1.10 M⊙ CO white dwarf, 50/50 C/O core, 0.016 M⊙ 22Ne, 0.0003 M⊙ 4He shell
   • Secondary star: 0.70 M⊙ CO white dwarf, 0.32 M⊙ carbon, 0.36 M⊙ oxygen, 0.011 M⊙ 22Ne, 0.009 M⊙ helium shell
2. Evolution Process:
   • Avoids double detonation scenario (primary helium shell too small)
   • System contains almost no 4He when secondary approaches disruption
   • Secondary impact on primary surface directly triggers carbon ignition
   • Orbital velocity before explosion ~2700 km/s
3. Explosion Results:
   • Primary star completely exploded
   • Secondary star ~0.16 M⊙ remains bound, velocity ~2800 km/s
   • ~1.64 M⊙ material ejected

1D Stellar Evolution Simulation (MESA)

1. Initialization Method:
   • Uses entropy and composition profiles from Arepo final state
   • Averages over equipotential surfaces (rather than spherically symmetric averaging)
   • Defines four models with different mass cutoffs: 0.100, 0.120, 0.140, 0.154 M⊙
2. Physical Processes:
   • Nuclear burning, convection, thermohaline mixing
   • 56Ni decay heating (though content is minimal)
   • Rotation neglected (angular momentum mostly lost during disruption)

Technical Innovations

Improved Velocity Calculation

Traditional D6 scenario assumes explosion at Roche lobe overflow, but in violent mergers the secondary can approach the primary more closely. New orbital velocity estimate:

$v_{\text{circ}} = \sqrt{(1+q)} v_{\text{Roche}}$

where q is the mass ratio, explaining why lower-mass white dwarfs can also reach ~2000 km/s velocities.

Thermal Timescale Analysis

Local thermal timescale: $t_{\text{th}} = \frac{H^2}{D_{\text{th}}}$

where H is the pressure scale height and Dth is the thermal diffusion coefficient. This timescale (~10^5-10^6 years) matches observed kinematic ages, rather than the traditional Kelvin-Helmholtz timescale.

Experimental Setup

Numerical Simulation Parameters

1. Arepo Settings:
   • 3D magnetohydrodynamics code
   • Tracks evolution for 1000 seconds post-explosion
   • Includes detailed nuclear reaction networks and equations of state
2. MESA Settings:
   • Uses approx21_plus_co56 nuclear reaction network
   • FreeEOS (surface) + HELM + Skye (white dwarf core) equations of state
   • OPAL opacity tables
   • Evolution to 100 Myr

Observational Data Comparison

Uses eight known D6 candidate objects, focusing on:

• D6-1, D6-3: Cooler hypervelocity stars (Teff < 9000K)
• J1235, J1637: Newly discovered similar objects
• Observed values of surface gravity, effective temperature, and radius

Experimental Results

Main Results

Physical Parameter Matching

1. Velocity Range: Model-predicted 2800 km/s consistent with D6-1, D6-3 observations
2. Radius Evolution: 0.154 M⊙ model inflates to ~0.1 R⊙ on timescales of 10^5-10^6 years
3. Temperature Range: Evolution trajectory passes through Teff~4000-8000K region, matching cooler hypervelocity stars

HR Diagram and Kiel Diagram Analysis

• Evolution trajectories of lower-mass models (0.10-0.14 M⊙) pass through positions of D6-1 and D6-3
• First theoretical explanation for the existence of these cooler hypervelocity stars

Key Findings

Degeneracy Boundary

Models show stars near the degeneracy boundary in early evolution, with ~50% of stellar mass near the degeneracy limit (EF = kBT), significantly different from traditional white dwarfs.

Surface Contamination

Unlike previous D6 models, violent merger-produced escaping stars show minimal 56Ni surface contamination (only ~1% of surface mass fraction), better matching observations.

Timescale Explanation

Local thermal timescale (10^4-10^6 years) matches observed kinematic ages, resolving the long-standing timescale problem.

Related Work

Major Research Directions

1. D6 Scenario (Shen et al. 2018; Pakmor et al. 2022): Double degenerate double detonation via Roche lobe overflow
2. Merger-Disruption Scenario (Glanz et al. 2025): Merger disruption of HeCO white dwarfs
3. Observational Discoveries (El-Badry et al. 2023; Hollands et al. 2025): New hypervelocity star samples

Advantages of This Work

• First explanation of cooler hypervelocity stars (Teff < 9000K)
• Provides new velocity generation mechanism
• Timescale matches observations
• Predicts dim, cool escaping stars in supernova remnants

Conclusions and Discussion

Main Conclusions

1. Violent merger mechanism can produce higher-velocity and longer-inflation hypervelocity stars
2. Combined effects of mass loss and heating enable low-mass remnants to maintain inflated states over long periods
3. D6-1 and D6-3 likely originate from violent mergers rather than traditional D6 scenario
4. Supernova remnant searches should focus on dimmer, cooler candidates

Limitations

1. 3D to 1D Mapping: Spherical symmetry approximation may lose important 3D effects
2. Mass Cutoff Uncertainty: Boundary definition affects final stellar mass
3. Rotation Neglected: While reasonable, may affect long-term evolution
4. Limited Sample: Only explains part of observed sample

Future Directions

1. Parameter Space Exploration: Systems with different mass ratios and compositions
2. Extended Timescale Simulations: Extend 3D simulations to longer timescales
3. Observational Verification: Search for predicted dim escaping stars in supernova remnants
4. Statistical Analysis: Proportion of violent mergers in Type Ia supernovae

In-Depth Evaluation

Strengths

1. High Innovation: First proposal and verification of violent merger mechanism, filling theoretical gaps
2. Rigorous Methodology: Combines 3D hydrodynamics and 1D stellar evolution with sound technical approach
3. Problem-Oriented: Addresses observational phenomena unexplained by existing theory with clear scientific objectives
4. Credible Results: Multiple independent model verifications with clear physical picture
5. Predictive Value: Provides specific search strategies for future observations

Weaknesses

1. Computational Limitations: Limited 3D simulation time (1000 seconds) may miss long-term effects
2. Approximations: 1D spherical symmetry and mass cutoff introduce uncertainties
3. Sample Limitations: Primarily explains 2-3 objects with limited statistical significance
4. Mechanism Universality: Conditions and frequency of violent merger occurrence remain unclear

Impact

1. Theoretical Contribution: Provides new physical mechanism for hypervelocity star research
2. Observational Guidance: Changes search strategy for escaping stars in supernova remnants
3. Methodological Exemplar: 3D+1D combined approach has broader applicability
4. Field Advancement: Promotes research in white dwarf binary evolution and Type Ia supernova progenitors

Applicable Scenarios

1. Theoretical Research: White dwarf binary evolution, supernova explosion mechanisms
2. Observational Planning: Hypervelocity star searches, supernova remnant studies
3. Numerical Methods: Technical reference for multi-dimensional stellar evolution simulations
4. Educational Applications: Typical case study for extreme stellar physics

References

This paper cites important literature in the field, including:

• Shen et al. (2018a,b): Foundational work on D6 scenario
• Pakmor et al. (2022): 3D double detonation simulations
• Glanz et al. (2025): Merger-disruption mechanism
• El-Badry et al. (2023): Latest observational discoveries
• Hollands et al. (2025): Discovery and analysis of J1637

Overall Assessment: This is a high-quality astrophysics paper with significant contributions in theoretical innovation, numerical methodology, and observational interpretation. By introducing the violent merger mechanism, it successfully resolves long-standing observational puzzles in the field, providing a new physical picture for hypervelocity star formation and evolution. Despite some technical limitations, its scientific value and impact are substantial, representing important progress in the field.