Fast Radio Bursts from White Dwarf Binary Mergers: Isolated and Triple-Induced Channels

Basic Information

• Paper ID: 2511.18678
• Title: Fast Radio Bursts from White Dwarf Binary Mergers: Isolated and Triple-Induced Channels
• Authors: Cheyanne Shariat (Caltech), Claire S. Ye (CITA), Smadar Naoz (UCLA), Sanaea Rose (Northwestern)
• Classification: astro-ph.HE (High Energy Astrophysics), astro-ph.GA (Galactic Astrophysics), astro-ph.SR (Stellar Astrophysics)
• Submission Date: November 24, 2025
• Paper Link: https://arxiv.org/abs/2511.18678v1

Abstract

Fast Radio Bursts (FRBs) are detected in both young and old stellar populations, indicating multiple formation pathways. This paper investigates a delayed channel: neutron star formation through merger or accretion-induced collapse (AIC) of massive white dwarfs (WDs) that emit FRBs. Through simulations of realistic stellar populations containing both binary and triple systems, the study finds: (i) triple-star dynamics opens new merger channels inaccessible to isolated binaries, significantly enhancing total merger rates; (ii) triple-induced mergers broaden the delay-time distribution (DTD), producing long-delay (≳1-8 Gyr) events and short-delay (≲100 Myr) rapid mergers that are essentially independent of metallicity; (iii) long delays naturally explain FRBs in quiescent host galaxies and galactic halos in old environments; (iv) after convolution with cosmic star formation history, the binary channel tracks the star formation rate (z_peak ~ 2), while the triple channel peaks later (z_peak ~ 1), yielding a combined local source rate R₀ ≈ 2×10⁴ Gpc⁻³ yr⁻¹ consistent with observations; (v) applying the same framework to Type Ia supernovae, triple systems extend the DTD tail, doubling the Ia efficiency relative to binaries.

Research Background and Motivation

Core Questions

1. The FRB Origin Mystery: Since the first FRB discovery in 2007, the physical origin of fast radio bursts has remained unsolved. Although the detection of an FRB from the magnetar SGR 1935+2154 in the Milky Way established magnetars as a primary channel, the phenomenological diversity of FRBs (repeating/non-repeating, polarization properties, host galaxy demographics) strongly suggests multiple progenitor channels.
2. Binary Nature of Host Environments: Observations show FRBs appearing both in star-forming galaxies (supporting a "prompt" channel with young magnetars from core-collapse supernovae) and in quiescent elliptical galaxy outskirts and globular clusters (such as the FRB source in M81), the latter suggesting a "delayed" channel.
3. Potential of the White Dwarf Merger Channel: Mergers or accretion-induced collapse (AIC) of massive white dwarfs may produce highly magnetized neutron stars, representing a promising FRB engine. However, existing research primarily focuses on isolated binaries, overlooking the important role of triple systems.

Research Significance

1. Ubiquity of Triple Systems: Approximately 35% of ~2M⊙ stars and ~50% of ~8M⊙ stars (progenitors of massive white dwarfs) reside in hierarchical triple systems, making triples a natural and potentially dominant environment for FRB progenitor formation.
2. Uniqueness of Dynamical Mechanisms: The Kozai-Lidov (EKL) mechanism in hierarchical triples can excite large eccentricity oscillations in the inner binary, leading to tidal circularization, mass transfer, or stellar mergers, providing evolutionary pathways for the inner binary that isolated evolution cannot achieve.
3. Need for Observational Constraints: The volumetric FRB rate as a function of redshift can distinguish different formation scenarios, but current constraints are weak due to ambiguities in the DM-redshift relation, selection biases, and the repeater problem.

Limitations of Existing Approaches

1. Neglect of Triple Contributions: Previous studies (e.g., Kashiyama et al. 2013; Lu et al. 2022; Kremer et al. 2021) primarily discussed white dwarf binary mergers as FRB progenitors without systematically considering triple dynamics.
2. Simplified Initial Conditions: Although Rajamuthukumar et al. (2023) studied triple effects on Type Ia supernovae, they allowed all CO+He WD mergers to produce SNe Ia (inconsistent with explosion models), and simulation time cutoffs may miss long-term dynamics.
3. Lack of Unified Framework: Missing is a systematic study simultaneously simulating binaries and triples across different metallicities, incorporating realistic stellar population initial conditions.

Core Contributions

1. First Systematic Triple FRB Study: Presents the first systematic study of FRB candidate formation including both isolated binaries and hierarchical triples, employing realistic initial conditions based on Gaia observations.
2. Identification of Unique Triple Channels: Discovers three main FRB formation pathways opened by triple dynamics:
   • Eccentric RG+MS systems (93%)
   • Eccentric WD+RG/MS systems (6%)
   • Eccentric WD+WD direct collisions (~1%)
3. Extension of Delay-Time Distribution: Demonstrates that triple-induced mergers extend the DTD tail to ≳8 Gyr, essentially independent of metallicity, naturally explaining FRBs in old environments.
4. Redshift Evolution Predictions: Predicts that the binary channel tracks the star formation rate (z_peak ~ 2), while the triple channel peaks at z ~ 1, yielding a combined local source rate R₀ ≈ 2×10⁴ Gpc⁻³ yr⁻¹ consistent with observations.
5. Additional Results for Type Ia Supernovae: Applying the same framework, finds that triples double the SNe Ia efficiency relative to binaries (~2×10⁻³ SNe M⊙⁻¹), with rates and redshift evolution in good agreement with observations.

Methodology Details

Task Definition

Objective: Through evolution of realistic stellar populations (containing both binaries and triples), identify white dwarf collapse pathways capable of producing FRB candidates and predict their delay-time distributions, rates, and redshift evolution.

Inputs:

• Initial conditions for hierarchical triple systems (masses, orbital parameters, metallicity)
• Cosmic star formation history
• Metallicity evolution

Outputs:

• Formation channels for FRB candidates
• Delay-time distribution (DTD)
• Redshift-dependent volumetric source rate R(z)

Triple Population Initial Conditions

Based on empirical distributions derived from Gaia observations by Shariat et al. (2025a):

1. Inner Binary Parameters (M₁, q_in, P_in, e_in): Sampled from the covariant distribution of Moe & Di Stefano (2017)
2. Outer Orbital Parameters:
   • Period: Lognormal distribution
   • Eccentricity: Thermal distribution
   • Mass ratio q_out: Power law (slope -1.4)
3. Mutual Inclination i_mut: Isotropic distribution
4. Stability Constraints:
   • Octupole hierarchical criterion: ϵ = (a₁/a₂)(e₂/(1-e₂²)) < 0.1
   • Mardling & Aarseth stability criterion

Evolution Model

1. Triple Long-Term Evolution

• Solve hierarchical three-body long-term equations to octupole order approximation
• Include first-order post-Newtonian general relativistic precession (Ford et al. 2000)
• Tidal evolution employs equilibrium tide model (Hut 1980) with different prescriptions for convective and radiative zones

2. Single-Star Evolution

• Track evolution of all three stars using SSE code (Hurley et al. 2000)

3. Binary Interactions

• Remove binaries undergoing mass transfer from triple evolution
• Track subsequent evolution using COSMIC code (Breivik et al. 2020)
• Key settings:
  • qcflag=4: Critical mass ratio follows Belczynski et al. (2008)
  • don_lim=-2: Donor mass loss rate follows Claeys et al. (2014)

4. Evolution Timescales

• Evolve each system for random time U(0, 12.5) Gyr, simulating star formation history
• Simulate three metallicities: Z⊙, 0.1Z⊙, 0.01Z⊙

FRB Candidate Identification

Employ an agnostic approach, exploring three potential FRB formation pathways:

1. O/Ne+Any Channel: O/Ne white dwarf merger/AIC with any companion (super-Chandrasekhar collapse)
2. WD+WD Channel: Merger of two CO white dwarfs (any mass)
3. Massive WD+WD Channel: Massive WD (M > 0.85M⊙, O/Ne or CO) merger with another CO WD

Key Assumptions:

• Do not require Chandrasekhar mass (but 70-100% of systems have total mass >1.35M⊙)
• Flag final outcomes of evolutionary trajectories without presetting FRB emission mechanism
• Treat these mergers as candidate pathways to quantify intrinsic rates and DTD

Redshift Evolution Calculation

Convolve DTD with cosmic star formation rate density:

$\psi(z) = 0.015 \frac{(1+z)^{2.7}}{1+[(1+z)/2.9]^{5.6}} \, M_\odot \, \text{yr}^{-1} \, \text{Mpc}^{-3}$

Combined with cosmic mean metallicity evolution (Madau & Fragos 2017):

$\log_{10}(\bar{Z}/Z_\odot) = 0.153 - 0.074z^{1.34}$

Assume multiplicity fractions:

• Isolated binaries: f_binary = 34%
• Triples: f_triple = 11.6%

Local rate estimation: $R_0 \approx 10^4 \, \text{Gpc}^{-3} \text{yr}^{-1} \left(\frac{\Gamma}{10^{-3} \text{yr}^{-1}}\right)$

where Γ is the Galactic rate, dependent on star formation rate and channel efficiency f_channel.

Experimental Setup

Simulation Parameters

1. Population Size: Evolve large numbers of triple systems to obtain statistically significant samples
2. Metallicity: Three values: Z⊙, 0.1Z⊙, 0.01Z⊙
3. Evolution Time: Random time 0-12.5 Gyr for each system
4. Control Experiments: Re-evolve the same inner binaries as isolated binaries using COSMIC

Evaluation Metrics

1. Formation Efficiency: Fraction of stellar systems producing FRB candidates (f_channel)
2. Delay-Time Distribution (DTD): Distribution of time from star formation to FRB candidate formation
3. Volumetric Source Rate R(z): Redshift-dependent FRB candidate formation rate Gpc⁻³ yr⁻¹
4. Peak Redshift z_peak: Redshift corresponding to maximum source rate

Comparison Baselines

1. Isolated Binary Evolution: Binary populations with same initial conditions and evolution codes
2. Observed FRB Rates:
   • James et al. (2022): 8.7×10⁴ Gpc⁻³ yr⁻¹
   • Shin et al. (2023): 7.3×10⁴ Gpc⁻³ yr⁻¹
   • Meng & Deng (2025): 3.9×10⁴ Gpc⁻³ yr⁻¹
   • Chen et al. (2024): 1.1×10⁴ Gpc⁻³ yr⁻¹
   • Hashimoto et al. (2022): 1.5×10⁴ Gpc⁻³ yr⁻¹

Experimental Results

Formation Channel Analysis

Three Main Triple FRB Formation Pathways (Figure 1):

1. Eccentric RG+MS (93%):
   • First mass transfer when primary evolves to red giant
   • Mass transfer can be stable or unstable (latter causes common envelope)
   • Results in separated WD+MS binary
   • Subsequent WD+MS mass transfer can lead to AIC or tight WD+WD binary
2. Eccentric WD+RG/MS (6%):
   • Initial wide inner binary (~>10 au) evolves independently to WD
   • EKL dynamics excite extreme eccentricity
   • After tidal circularization begins WD+MS mass transfer
   • Leads to AIC or WD+WD merger
3. Eccentric WD+WD (~1%):
   • Initial wide WD+WD binary (a_in ≳ 10 au)
   • Octupole EKL oscillations drive e_in ~> 0.999
   • Periapsis drops to ~< 10⁻² au, triggering direct collision
   • Though small in fraction, significantly contributes to long-delay (~> 1 Gyr) mergers

Metallicity Effects:

• Solar metallicity: Triple systems produce more FRB candidates (Table 1: 4.3×10⁻³ vs 3.2×10⁻³ for WD+WD)
• Low metallicity (Z ≲ 0.1Z⊙): Triple merger rates decrease due to early tidal circularization causing extreme mass ratios (q ≪ 1), making mass transfer more unstable

Delay-Time Distribution (Figure 2)

WD+WD Channel (Left Column)

Solar Metallicity:

• Triple DTD significantly wider, tail extends to ≳8 Gyr
• Binary DTD concentrated at short delays (< 1 Gyr)
• More triple events overall (enhancement factor ~1.3)

Physical Mechanism:

• Triples: Initial wide orbit (P_in ~ 10⁴ days) → EKL excites eccentricity → Tidal circularization to P_in ~ 10 days → Mass transfer → WD+WD → Gravitational wave merger
• Multi-stage process naturally produces long delays

Low Metallicity (0.1-0.01Z⊙):

• Reduced triple merger numbers (but dynamical channels preserved)
• Eccentric WD+WD collisions dominate long-delay tail
• Binaries more efficient at low Z (wider orbital range avoiding early unstable mass transfer)

Massive WD+WD Channel (Middle Column)

• Similar trends to WD+WD but fewer total events
• Solar Z: Triple contribution 2.3×10⁻³ vs binary 1.6×10⁻³
• More prominent long-delay tail

O/Ne WD+Any Channel (Right Column)

• Both binary and triple DTDs biased toward short delays (< 1 Gyr)
• Modest tail extension to several Gyr (accretion when secondary evolves to giant)
• Weaker metallicity dependence

Redshift Evolution (Figure 3)

Per-Channel Results (Upper Panel)

Binary Channel:

• Broadly tracks cosmic star formation history
• WD+WD peak: z_peak ~ 1.9
• O/Ne WD+any peak: z_peak ~ 1.3
• Reflects short-delay characteristics

Triple Channel:

• Systematic shift to lower redshifts
• WD+WD peak: z_peak ~ 1.0
• Massive WD+WD peak: z_peak ~ 0.9
• Reflects long-delay characteristics

Combined Results (Lower Panel)

Local Rate (z=0):

• WD+WD: R₀ = 1.8×10⁴ Gpc⁻³ yr⁻¹
• O/Ne WD+any: R₀ = 1.3×10⁴ Gpc⁻³ yr⁻¹
• Massive WD+WD: R₀ = 8.9×10³ Gpc⁻³ yr⁻¹

Peak Rate and Redshift (Table 2):

• WD+WD combined: R_peak = 8.7×10⁴ Gpc⁻³ yr⁻¹ at z_peak = 1.62
• Consistent with observations (Shin et al. 2023: 7.3×10⁴; Hashimoto et al. 2022: 1.5×10⁴)

Key Predictions:

• If triples significantly contribute, should observe excess relative to binary predictions at z ~ 1
• Binary+triple combined peak lies between pure binary (z ~ 2) and pure triple (z ~ 1)

Type Ia Supernova Results (Appendix A, Figure 4)

"Physical" Model (Massive CO WD primary):

• Binary efficiency: 6×10⁻⁴ SNe M⊙⁻¹
• Triple efficiency: 1.2×10⁻³ SNe M⊙⁻¹
• Combined: ~2×10⁻³ SNe M⊙⁻¹ (consistent with observed 1.3×10⁻³)

Local Rate:

• R₀,Ia ≈ 3×10⁴ Gpc⁻³ yr⁻¹

Redshift Evolution:

• Good agreement with observational compilation from Maoz & Graur (2017)
• Triples extend DTD tail, doubling Ia efficiency relative to binaries

Related Work

FRB Progenitor Models

1. Prompt Channels:
   • Young magnetars (core-collapse supernovae): Metzger et al. 2017; Bochenek et al. 2020
   • Delays ≲100 Myr, track star formation
2. Delayed Channels:
   • AIC of massive WDs: Nomoto & Kondo 1991; Schwab 2021
   • Compact object mergers: Wang et al. 2016; Sridhar et al. 2021
   • White dwarf mergers: Kashiyama et al. 2013; Lu et al. 2022; Kremer et al. 2021

Triple Dynamics Research

1. EKL Mechanism: Kozai 1962; Lidov 1962; Naoz 2016
2. White Dwarf Binary Formation:
   • Accreting WD binaries: Knigge et al. 2022; Shariat et al. 2025b
   • Ultra-compact WD binaries: Toonen et al. 2016, 2020; Shariat et al. 2023
3. Compact Object Mergers:
   • Neutron stars/black holes: Antonini et al. 2014; Fragione & Loeb 2019; Stegmann et al. 2022

Triples in Type Ia Supernovae

1. Prior Work:
   • Thompson 2011; Katz & Dong 2012; Hamers et al. 2013
   • Toonen et al. 2018; Michaely 2021
   • Rajamuthukumar et al. 2023 (allowing all CO+He mergers, potentially overestimating)
2. Improvements in This Work:
   • Realistic initial conditions (Gaia-based)
   • Only massive primaries (consistent with explosion models)
   • No simulation time cutoff, capturing long-term dynamics

Conclusions and Discussion

Main Conclusions

1. Triples Open New Channels: Hierarchical triples provide evolutionary pathways through EKL mechanism inaccessible to isolated binaries, significantly enhancing FRB candidate merger rates at solar metallicity (~1.3× enhancement).
2. DTD Extension: Triple-induced mergers extend DTD tail to ≳8 Gyr, essentially independent of metallicity (via eccentric WD+WD collision channel), naturally explaining FRBs in quiescent galaxies and globular clusters.
3. Host Environment Diversity: Long delays cause triple-channel FRBs not to fully track star formation, appearing in old stellar populations, consistent with observed host diversity.
4. Redshift Evolution Predictions:
   • Binaries: z_peak ~ 2 (tracking star formation)
   • Triples: z_peak ~ 1 (delayed contribution)
   • Combined local rate R₀ ≈ 2×10⁴ Gpc⁻³ yr⁻¹ consistent with observations
5. Additional Type Ia Insights: Triples double SNe Ia efficiency, extend DTD, with rates and redshift evolution in good agreement with observations.

Limitations

1. Physical Uncertainties:
   • Equilibrium tide model may underestimate efficiency at large eccentricities (dynamical tides)
   • Circularization and contraction mechanisms in WD+WD final stages (Fuller & Lai 2012)
   • Uncertainties in common envelope physics
2. Unknown FRB Emission Mechanism:
   • No assumption of specific FRB emission mechanism
   • Reported rates are formation rates (lower limits), not accounting for repeaters
   • Whether CO+CO mergers truly produce FRBs remains uncertain (may be SNe Ia)
3. Complexity of Observational Comparison:
   • Observed FRB rates depend on luminosity function, DM-redshift relation assumptions
   • Different treatment of repeaters
   • Selection effects and completeness corrections
4. Higher-Order Multi-Body Effects:
   • Neglect hexadecapole ((a₁/a₂)⁵) and higher-order terms
   • Expected secondary impact on wide outer orbits and unequal mass systems

Future Directions

1. Observational Tests:
   • CHIME/FRB Catalog 2, CHORD, DSA-110, DSA-2000 will expand samples
   • Precise measurements of rates, luminosity functions, repeater fractions, redshift evolution
   • Arcsecond-scale localization to identify local environments and potential companions
2. Theoretical Improvements:
   • Include dynamical tide models
   • Explore different common envelope prescriptions
   • Study magnetic field generation mechanisms for FRB emission
3. Multi-Messenger Connections:
   • FRB-SNe Ia spatial coincidences could support WD merger picture
   • Combined with gravitational wave detection (LISA for WD binaries)
4. Environmental Constraints:
   • FRB searches in low-metallicity environments
   • Deep observations of z ~ 1 target systems
   • Statistical analysis of host galaxy ages and stellar population properties

In-Depth Evaluation

Strengths

1. Methodological Innovation:
   • First systematic combination of binary and triple evolution for FRB studies
   • Uses realistic initial conditions based on Gaia, enhancing credibility
   • Agnostic approach (not presetting FRB mechanism) strengthens robustness
2. Deep Physical Insights:
   • Clearly identifies three triple formation channels and their physical mechanisms
   • Reveals non-trivial metallicity effects on relative binary/triple contributions
   • EKL eccentricity excitation + tidal circularization two-stage process explains long delays
3. Rigorous Experimental Design:
   • Control experiments (same inner binaries as isolated binaries) effectively isolate triple effects
   • Systematic exploration across three metallicities
   • Self-consistent convolution with cosmic star formation and metallicity evolution
4. Testable Predictions:
   • Clear redshift evolution differences (binary z2 vs triple z1)
   • Quantitative local rate predictions (R₀ ≈ 2×10⁴ Gpc⁻³ yr⁻¹)
   • Expected host galaxy age and metallicity distributions
5. Cross-Disciplinary Contribution:
   • Same framework successfully applied to SNe Ia, validating methodology
   • Provides unified understanding for multiple astrophysical phenomena

Shortcomings

1. Simplified Physical Models:
   • Equilibrium tide may underestimate efficiency, especially at large eccentricities and WD+WD final stages
   • Common envelope evolution uncertainties (qcflag choice)
   • Neglects magnetic field evolution and angular momentum transfer details
2. FRB Mechanism Ambiguity:
   • Unclear which WD mergers truly produce FRBs (vs SNe Ia or others)
   • No consensus whether CO+CO mergers produce magnetized neutron stars
   • Reported rates are candidate formation rates, not actual FRB event rates
3. Observational Comparison Limitations:
   • Large variation in observed rate estimates (1.1-8.7×10⁴ Gpc⁻³ yr⁻¹)
   • Repeater impact on rates not addressed (may bias predictions high)
   • DM-redshift relation uncertainties insufficiently discussed
4. Statistical Uncertainties:
   • Simulation sample sizes and statistical errors not reported
   • Rare channel (eccentric WD+WD collision ~1%) statistical significance questionable
   • Lacks Monte Carlo uncertainty propagation
5. Initial Condition Assumptions:
   • Isotropic mutual inclination may be oversimplified
   • Neglects dynamical environment effects (e.g., clusters) on initial conditions

Impact

1. Field Contribution:
   • Establishes triples as key and previously overlooked FRB progenitor channel
   • Provides natural framework for explaining FRB host diversity
   • Shifts FRB progenitor models from "prompt vs delayed" dichotomy toward multi-channel picture
2. Practical Value:
   • Provides clear predictions for upcoming surveys (CHIME Catalog 2, DSA-2000)
   • Guides observational strategies (e.g., z~1 target selection, low-metallicity searches)
   • Host galaxy demographics serve as diagnostic tool for progenitor channels
3. Reproducibility:
   • Uses public codes (COSMIC, SSE)
   • Initial condition sampling tools publicly available (https://github.com/cheyanneshariat/gaia_triples)
   • Detailed method description and explicit parameter settings
4. Theoretical Significance:
   • Demonstrates universal importance of triple dynamics in multiple compact object formation
   • Provides transferable framework for other phenomena (SNe Ia, gravitational wave sources)
   • Challenges traditional stellar population synthesis models neglecting multiplicity

Applicable Scenarios

1. Direct Applications:
   • FRB progenitor channel identification and rate predictions
   • Host galaxy demographic interpretation
   • Survey strategy optimization
2. Extensible Domains:
   • Other compact object mergers (NS+NS, BH+BH)
   • X-ray binary formation
   • Gravitational wave source statistics
3. Observational Constraints:
   • Redshift distribution fitting
   • Host galaxy age-FRB correlations
   • Metallicity dependence
4. Theoretical Exploration:
   • Different common envelope prescription impacts
   • Magnetic field generation mechanisms
   • Systematic multiplicity effects on population evolution

Key References

1. FRB Observations:
   • Lorimer et al. 2007: First FRB discovery
   • CHIME/FRB Collaboration et al. 2020: Milky Way magnetar FRB
   • Bhardwaj et al. 2021b; Kirsten et al. 2022: M81 globular cluster FRB
2. Triple Dynamics:
   • Kozai 1962; Lidov 1962: EKL mechanism
   • Naoz 2016: Review
   • Shariat et al. 2025a: Gaia-based triple initial conditions
3. White Dwarf Mergers:
   • Nomoto & Kondo 1991: AIC producing neutron stars
   • Kashiyama et al. 2013: WD mergers as FRB progenitors
   • Schwab 2021: Magnetized WD formation
4. Stellar Population Synthesis:
   • Hurley et al. 2000, 2002: SSE code
   • Breivik et al. 2020: COSMIC code
   • Moe & Di Stefano 2017: Stellar multiplicity statistics
5. Cosmology:
   • Madau & Dickinson 2014: Cosmic star formation history
   • Madau & Fragos 2017: Metallicity evolution

Overall Assessment: This is a high-quality research work with significant impact. Through rigorous numerical simulations and physical analysis, it systematically reveals for the first time the critical role of triple dynamics in FRB progenitor formation. The methodology is innovative, predictions are explicit and testable, providing a powerful theoretical framework for understanding FRB diversity and upcoming observational surveys. Main limitations lie in FRB emission mechanism uncertainties and simplifications in certain physical processes, but these do not compromise the robustness of core conclusions. This work not only advances FRB origin research but also provides strong evidence for the universal importance of triples in multiple astrophysical phenomena.