Collapse of Magnetized White Dwarfs as site of Heavy Element Formation and Kilonova Signal

This study presents the first end-to-end simulation demonstrating that the magnetized accretion-induced collapse of a white dwarf produces neutron-rich ejecta capable of robust $r$-process nucleosynthesis, generating a lanthanide-rich kilonova signal that matches the observed properties of AT 2023vfi/GRB 230307A without parameter tuning.

The Gist

The Big Picture: A Cosmic "Heavy Metal" Factory

Imagine a white dwarf star as a dense, dead ember of a star that has stopped burning. Usually, when these stars die, they just fade away or explode into a small, bright flash of light. But this paper asks a "What if?" question: What happens if that dead star is spinning like a top and has a magnetic field as strong as a giant magnet?

The authors found that under these specific conditions, the star doesn't just die quietly. It collapses into a neutron star and launches a massive, high-speed explosion that acts like a cosmic factory, forging the heaviest elements in the universe (like gold, platinum, and uranium) and lighting up the sky in a way we've never seen before.

The Cast of Characters

To understand how they figured this out, imagine a three-stage assembly line:

1. The Crash (The Simulation): They used a supercomputer to simulate the star collapsing. Think of this as a high-speed crash test, but instead of cars, it's a star imploding under its own gravity.
2. The Forge (The Chemistry): As the star explodes, they tracked the atoms being cooked up in the heat. This is like watching a chef mix ingredients in a pot, but the "ingredients" are subatomic particles, and the "heat" is so intense it creates new elements.
3. The Show (The Light): Finally, they calculated how the light from this explosion would look to an observer on Earth. This is like predicting how a fireworks display will look from different seats in the stadium.

The Plot Twist: Why This Time is Different

In the past, scientists thought that when a white dwarf collapses, the explosion is bathed in a flood of "neutrinos" (ghostly particles that pass through everything). These neutrinos usually act like a cosmic reset button, turning the heavy, neutron-rich material back into lighter, proton-rich stuff (like iron).

The Paper's Discovery:
Because this white dwarf was spinning fast and had a super-strong magnetic field, it acted like a high-pressure water hose.

• The Analogy: Imagine trying to wash a car with a garden hose (the old models). The water (neutrinos) soaks everything, changing the dirt. But in this new model, the magnetic field acts like a firehose nozzle. It shoots the material out so fast that the "neutrino rain" can't touch it before it's already far away.
• The Result: The material escapes while it is still "neutron-rich." This is the secret sauce needed to create the heavy elements (the "r-process") that make up gold and platinum.

The Evidence: Matching the Crime Scene

The researchers didn't just make up a theory; they compared their simulation to a real cosmic event that happened recently: GRB 230307A.

• The Event: A long gamma-ray burst (a flash of high-energy light) was detected, followed by a "kilonova" (a fading glow of light).
• The Match: The authors ran their simulation and looked at the light curve (how bright it gets and fades) from different angles.
  • If you look at the explosion from the side (equator), it looks one way.
  • If you look from the top (pole), it looks different.
• The "Aha!" Moment: When they looked at their simulation from the top (pole-on view), it matched the real observations of GRB 230307A perfectly, without them having to tweak any numbers. It was like finding a fingerprint that matched the suspect's hand exactly.

This suggests that GRB 230307A wasn't caused by two neutron stars smashing together (the usual suspect), but by a spinning, magnetized white dwarf collapsing.

The Aftermath: A Radio Silence

The paper also looked at whether we could hear this event with radio telescopes.

• The Analogy: Imagine the explosion is a loud party inside a thick, foggy room. At first, the "fog" (the expanding gas) is so thick that no radio waves can get out.
• The Finding: It takes about a month for the fog to clear enough for radio signals to escape. If there is a "magnetar" (a super-magnetic neutron star) left behind at the center, it might send out a radio burst (like a Fast Radio Burst), but we have to wait for the fog to lift before we can hear it.

Why Should We Care?

1. Where does Gold come from? This paper suggests that spinning, magnetized white dwarfs are a major factory for heavy elements in the universe, not just colliding neutron stars.
2. New Clues for Astronomers: It gives astronomers a new "recipe" to look for. If they see a gamma-ray burst followed by a specific type of infrared glow, they might know they are watching a white dwarf collapse, not a merger.
3. The "Perfect" Match: The fact that their computer model matched a real-life event so well without any "fudging" of the numbers is a huge win for physics. It proves that our understanding of how these magnetic stars work is getting very accurate.

In short: A spinning, magnetized dead star can collapse and shoot out a jet of material so fast that it creates gold and platinum, lighting up the sky in a way that perfectly matches a real event we saw last year. It's a new chapter in the story of how the universe makes its heavy stuff.

Technical Summary

1. Problem Statement

The Accretion-Induced Collapse (AIC) of white dwarfs (WDs) is a potential pathway for neutron star formation and transient events. However, previous theoretical models of unmagnetized AICs predicted:

• Proton-rich ejecta: Neutrino irradiation raises the electron fraction ($Y_e$) to $\sim 0.43–0.50$, preventing heavy $r$-process nucleosynthesis.
• Low ejecta mass: Typically $10^{-3}–10^{-2} M_\odot$.
• Observational signature: Fast, blue optical transients dominated by $^{56}$Ni, distinct from the red, lanthanide-rich kilonovae associated with neutron star mergers.

Recent observations of long-duration Gamma-Ray Bursts (GRBs) with kilonova-like emission (e.g., GRB 230307A/AT 2023vfi) challenge the traditional association of long GRBs with massive star collapse. The central problem addressed by this paper is whether magnetized, rapidly rotating WDs can produce neutron-rich outflows capable of heavy $r$-process nucleosynthesis and generate observable kilonova signals consistent with these events. A critical gap existed in connecting dynamical simulations to self-consistent electromagnetic signatures (light curves and spectra) without relying on parametric heating prescriptions.

2. Methodology

The authors present the first end-to-end computational pipeline connecting general-relativistic hydrodynamics to synthetic observables. The workflow consists of three coupled stages:

A. General-Relativistic Neutrino-Magnetohydrodynamics (GR$\nu$MHD)

• Code: Gmunu.
• Setup: 2D axisymmetric simulations of a $1.5 M_\odot$ rapidly rotating WD ($a_r = 0.75$, $\Omega \approx 5$ Hz) with strong initial magnetic fields ($B_{pol} = B_{tor} = 10^{12}$ G).
• Physics: Includes the LS220 equation of state, energy-integrated two-moment neutrino transport, and ideal MHD.
• Innovation: Unlike previous studies, this simulation evolves both hemispheres independently (no equatorial symmetry) over a large domain ($3 \times 10^5$ km) to capture full ejecta dynamics up to $t \sim 1$ s post-bounce.

B. Radiation Hydrodynamics with In-Situ Nucleosynthesis

• Code: kNECnn (an updated version of SNEC coupled with SkyNet).
• Process: 35 angular profiles extracted from the GR$\nu$MHD output are evolved as 1D Lagrangian problems for $\sim 30$ days.
• Nucleosynthesis: Uses the SkyNet network (7,836 isotopes, $\sim 140,000$ reactions) coupled in-situ. This tracks the full $r$-process evolution rather than using parametric heating.
• Heating: Calculates thermalized heating rates ($\dot{q}_{nucl}$) from $\gamma$-rays, $\beta$-decay, and fission fragments, accounting for spatially varying thermalization efficiencies.

C. Monte Carlo Radiative Transfer

• Code: Sedona.
• Process: Converts the hydrodynamic and nucleosynthetic outputs into synthetic light curves and spectra.
• Key Feature: The code was modified to accept spatially and temporally resolved heating rates derived directly from kNECnn, avoiding the uniform heating approximations used in previous kilonova studies.
• Opacity: Uses the Sobolev approximation with detailed atomic data (HULLAC/CMFGEN) for elements up to $Z=87$, capturing lanthanide/actinide opacities.

3. Key Contributions

1. First End-to-End AIC Pipeline: The study bridges the gap between 2D GR$\nu$MHD simulations and observable kilonova signals, moving beyond parametric assumptions to self-consistent nucleosynthesis and radiative transfer.
2. Magnetic Field as the Deciding Factor: The work demonstrates that strong magnetic fields ($10^{12}$ G) are essential for ejecting material on dynamical timescales (before neutrino irradiation can re-leptonize the flow), preserving low electron fractions ($Y_e \sim 0.24$) necessary for heavy $r$-process.
3. Spatially Resolved Heating: The implementation of non-uniform, composition-dependent heating rates in radiative transfer, revealing that heating varies by an order of magnitude across the ejecta.

4. Key Results

Dynamics and Ejecta Properties

• Mass and Energy: The simulation produces $\approx 0.18 M_\odot$ of ejecta with kinetic energy $\approx 8 \times 10^{50}$ erg.
• Composition: The ejecta are neutron-rich ($\langle Y_e \rangle \sim 0.24$). Approximately $0.14 M_\odot$ has $Y_e \leq 0.25$, favorable for strong $r$-process.
• Morphology: A magnetically dominated polar jet is surrounded by neutron-rich equatorial/mid-latitude outflows. The magnetic pressure dominates gas pressure ($\beta_{mag} < 1$) in the jet region.

Nucleosynthesis

• Heavy Element Production: The ejecta undergo robust $r$-process nucleosynthesis, producing elements up to and beyond the third peak ($A \gtrsim 195$).
• Abundance Pattern: The mass-weighted abundance pattern closely matches solar $r$-process residuals in the second peak, rare-earth region, third peak, and actinides.
• Dominant Elements: The four most abundant elements by mass are Xenon (Xe), Helium (He), Platinum (Pt), and Tellurium (Te), constituting $>50\%$ of the ejecta.
• Lanthanides: The lanthanide mass fraction is $X_{lan} \approx 6\%$, ensuring high opacity.

Kilonova Signal

• Light Curves: The emission is lanthanide-rich and dominated by near-infrared (NIR) wavelengths, peaking around 4 days. This contrasts sharply with the fast, blue transients predicted by unmagnetized models.
• Comparison to GRB 230307A: The synthetic light curves for polar viewing angles ($\theta_{obs} \lesssim 30^\circ$) show striking agreement with the photometric observations of AT 2023vfi (associated with GRB 230307A) in LSST and JWST bands.
  • No Parameter Tuning: The agreement is achieved using only the viewing angle as a free parameter; all ejecta properties are self-consistently derived from the simulation.
  • Self-Consistency: The near-polar angle required to match the kilonova is the same geometry required to observe the prompt $\gamma$-ray burst (jet axis), providing a robust physical consistency check.
• Spectra: Synthetic spectra show broad emission features in the $\sim 1.5–2.5 \mu$m range, consistent with the heavy $r$-process composition, though specific line identification is cautioned against due to theoretical atomic data limitations.

Radio Transparency

• The ejecta remain opaque to radio frequencies for $\sim 30$ days.
• Polar directions become transparent first ($\tau_{ff} < 1$ at $\nu \sim 5$ GHz by $t \approx 20$ days), suggesting that radio follow-up of AIC candidates should target $\sim 30$ days post-explosion.

5. Significance

• New Progenitor Channel: This work establishes magnetized AIC as a viable and compelling progenitor channel for long-duration GRBs with kilonova signatures, challenging the exclusive association of such events with neutron star mergers.
• Heavy Element Source: It identifies magnetized AICs as a significant, potentially dominant, source of heavy $r$-process elements (including actinides) in the universe, capable of reproducing the solar abundance pattern.
• Observational Diagnostics: The study provides specific predictions (NIR-dominated light curves, polar viewing geometry, delayed radio transparency) that can be used to distinguish AIC-powered kilonovae from those produced by binary neutron star mergers in future multi-messenger observations.
• Theoretical Advancement: By demonstrating that magnetic fields can overcome neutrino re-leptonization to produce neutron-rich outflows, the paper resolves a long-standing discrepancy between unmagnetized AIC models and the requirements for heavy element synthesis.