Probing the Dispersion and Rotation Measure Contributions from Supernova Remnants in Fast Radio Burst Source Environments with 1D SNR Simulation

This paper employs 1D supernova remnant simulations to demonstrate that the secular decrease in dispersion measure observed in repeating fast radio bursts is primarily driven by the expansion of unshocked ejecta from young magnetars, while shock-amplified magnetic fields account for rotation measure variations, thereby supporting a young core-collapse supernova origin for these sources and highlighting the necessity of modeling local environments for accurate cosmological inferences.

The Gist

Here is an explanation of the paper, translated into everyday language with creative analogies.

The Big Picture: Listening to the Universe's "Static"

Imagine the universe is a giant radio station. Fast Radio Bursts (FRBs) are like incredibly loud, short bursts of static coming from very far away. Astronomers love these bursts because they act as cosmic messengers.

As these radio waves travel through space to reach Earth, they have to pass through clouds of gas and plasma (ionized gas). This is like driving through a thick fog. The thicker the fog, the more the radio signal gets "stretched out" or delayed. Astronomers measure this delay using something called Dispersion Measure (DM).

Usually, this delay is constant. But recently, astronomers noticed something weird with two specific FRBs (FRB 20190520B and FRB 20121102): their "fog" is getting thinner over time. The delay is decreasing. This suggests the source of the radio burst is inside a changing, evolving environment, likely a Supernova Remnant (SNR)—the expanding debris cloud left behind after a massive star explodes.

The Experiment: Simulating a Cosmic Explosion

The authors of this paper wanted to figure out exactly what kind of explosion creates this specific type of "thinning fog." They built a 1D (one-dimensional) computer simulation.

Think of their simulation like a cross-section of a balloon being inflated. Instead of a 3D balloon, they are looking at a slice through the center, tracking how the gas moves, heats up, and ionizes as it expands outward from the center.

They tested two main scenarios for the star that exploded:

1. The Solo Star (Single-Star): A massive star that lived and died alone, blowing off its outer layers slowly like a gentle breeze before exploding.
2. The Stripped Star (Binary-Stripped): A massive star in a partnership with a companion star. The companion star acted like a cosmic vacuum cleaner, sucking away most of the exploding star's hydrogen "skin" (envelope) before the explosion happened.

The Key Findings

Here is what they discovered, broken down into simple concepts:

1. The "Fog" is Mostly Unshocked Debris

You might think the most dramatic part of an explosion—the shockwave hitting the air—would be the main source of the radio delay.

• The Analogy: Imagine a firecracker going off in a room. You'd think the shockwave hitting the walls is the loudest part.
• The Reality: The authors found that the shocked region (the shockwave itself) contributes very little to the radio delay (less than 10 units of DM).
• The Real Culprit: The delay comes from the unshocked ejecta—the raw, expanding debris cloud behind the shockwave. It's like the smoke from the firecracker that hasn't hit the walls yet. This cloud is massive and dense, and as it expands, it gets thinner, causing the radio delay to drop over time.

2. The "Solo" vs. "Stripped" Star Difference

The type of star matters a lot.

• Solo Stars (SS): These leave behind a huge, heavy cloud of hydrogen-rich debris. This creates a thick fog (high DM).
• Stripped Stars (BS): Because their companion star stole their hydrogen skin, these explosions leave behind a much smaller, lighter cloud of helium-rich debris. This creates a thin fog (low DM).
• The Metaphor: If the Solo Star explosion is a giant snowplow pushing a massive pile of snow, the Stripped Star explosion is a bicycle pushing a small pile of leaves. The snowplow creates a much bigger obstacle for the radio waves.

3. The Magnetic "Compass" (Rotation Measure)

The paper also looked at Rotation Measure (RM), which tells us how twisted the magnetic field is in the gas.

• They found that only the Solo Star (11 solar masses) model could reproduce the strong magnetic twisting observed in FRB 20121102.
• The Stripped stars didn't create enough magnetic "twist" to match the observations.
• Why? The Solo star's explosion creates a denser shockwave, which amplifies the magnetic field more effectively, like squeezing a sponge to make the water (magnetic energy) flow faster.

4. When Can We See the Burst?

For a radio burst to escape a supernova, the surrounding gas must become transparent (clear) enough for the signal to pass through without being absorbed.

• The authors calculated that for most of their models, the "fog" clears up enough for the radio signal to escape within 10 to 70 years after the explosion.
• This explains why we see these bursts: they are "young" supernovae, but not too young. If they were too young, the gas would be too thick, and the signal would be trapped.

Why Does This Matter?

Astronomers use FRBs to measure the total amount of "normal matter" (baryons) in the universe. To do this, they have to subtract the "local fog" (the DM from the source's own environment) from the total delay.

• The Problem: If we don't understand the local fog, we might think the universe has more matter than it actually does, or we might get the distance wrong.
• The Solution: This paper shows that the local fog isn't a fixed number. It depends on whether the star was a "Solo" or "Stripped" type, and it changes rapidly over the first few decades of the explosion's life.

The Takeaway

This paper is like a forensic investigation of a cosmic crime scene. By simulating the explosion of different types of stars, the authors figured out that:

1. The "fog" causing the radio delay is mostly the raw debris, not the shockwave.
2. The specific "signature" of the fog (how fast it thins and how magnetic it is) points strongly to a young, solo star (around 11 times the mass of our Sun) as the culprit for the famous FRB 20121102.
3. To accurately map the universe using these radio bursts, we must account for the specific "personality" of the star that exploded.

In short: The universe is a radio, the stars are the speakers, and the gas is the static. To hear the music clearly, we need to know exactly what kind of speaker (star) is playing and how the static (gas) is changing.

Technical Summary

Here is a detailed technical summary of the paper "Probing the Dispersion Measure and Rotation Measure Contributions from Supernova Remnants in Fast Radio Burst Source Environments with 1D SNR Simulation."

1. Problem Statement

Fast Radio Bursts (FRBs) are transient radio signals whose Dispersion Measures (DM) serve as a probe for the ionized baryon content of the universe. However, a significant portion of the observed DM originates from the local environment of the FRB source ($DM_{source}$).

• The Challenge: Recent observations of repeating FRBs, specifically FRB 20190520B and FRB 20121102, reveal secular (long-term) changes in DM. FRB 20190520B shows a steady decline, while FRB 20121102 exhibits a two-stage trend (early mild rise followed by a late decline).
• The Hypothesis: These variations suggest a dense, dynamically evolving local environment, likely a young magnetar embedded in a Supernova Remnant (SNR).
• The Gap: Previous studies relied on simplified analytic models (e.g., self-similar solutions) that often assume idealized progenitor structures and constant density profiles. There is a lack of physically consistent, time-dependent simulations that account for realistic progenitor evolution (single-star vs. binary-stripped), detailed shock physics, and non-equilibrium ionization (NEI) to accurately predict $DM_{source}$ and its time derivatives ($dDM/dt$) and Rotation Measure (RM).

2. Methodology

The authors employ a forward-modeling approach using detailed 1D hydrodynamic simulations coupled with Non-Equilibrium Ionization (NEI) calculations.

• Progenitor Models:
  • Utilized stellar evolution models (based on MESA) for two evolutionary channels: Single-Star (SS) and Binary-Stripped (BS).
  • Simulated two representative Zero-Age Main Sequence (ZAMS) masses: 11 $M_\odot$ and 30 $M_\odot$.
  • BS models account for mass loss via Roche-lobe overflow (RLOF), resulting in stripped helium stars (Type Ib/c progenitors) with significantly lower ejecta masses compared to their SS counterparts (Type II progenitors).
• CSM Construction:
  • Instead of assuming a simple $r^{-2}$ wind profile, the Circumstellar Medium (CSM) was constructed by simulating the interaction of time-variable stellar winds with the interstellar medium (ISM), creating complex structures like shells and cavities.
• SNR Simulation:
  • Used a 1D, spherically symmetric HD+NEI code (based on the ChN framework).
  • Simulated the evolution from $t_{init} = 3$ years post-explosion up to ~500 years.
  • Included radiative cooling and detailed ionization states for 30 elements.
  • Note: Cosmic-ray acceleration and magnetic field evolution were disabled for the DM calculation but a simplified ram-pressure-based magnetic field model was used for RM estimation.
• Observables Calculation:
  • DM: Calculated by integrating electron density ($n_e$) along the line of sight, considering both shocked and unshocked regions.
  • Optical Depth ($\tau$): Calculated for electron scattering and free-free absorption to determine the "escape time" ($t_{esc}$) for GHz radio emission.
  • RM: Calculated using shocked-region profiles with magnetic fields parameterized as a fraction ($\epsilon_B$) of ram pressure.
  • Matching: The simulated $dDM/dt$ and RM slopes were matched against observational data for FRB 20190520B and FRB 20121102 to constrain the age of the source ($t_{occur}$).

3. Key Contributions

• Realistic Progenitor Comparison: First detailed comparison of SS vs. BS progenitor channels in the context of FRB DM evolution, highlighting how binary stripping drastically alters CSM density and ejecta mass.
• Unshocked Ejecta Dominance: Demonstrated that while the shocked shell contributes little to the total DM, the unshocked ejecta is the dominant source of time-varying DM, following a power-law decay ($DM \propto t^{-\alpha}$ with $\alpha \approx 1.8\text{--}1.9$).
• Ionization Sensitivity: Showed that the DM contribution is highly sensitive to the ionization fraction of the unshocked ejecta ($\chi_{e,unsh}$), which is treated as a parameter ranging from weakly ionized ($\sim 0.01$) to fully ionized ($1.0$).
• RM Constraints: Provided a rigorous test of whether SNR shock amplification alone can explain the high RM of FRB 20121102.

4. Key Results

A. Dispersion Measure (DM) Evolution

• Magnitude: The shocked region contributes only a limited DM ($\lesssim 10$ pc cm$^{-3}$). The dominant component is the unshocked ejecta.
• Time Evolution: The DM declines as $t^{-1.8}$ to $t^{-1.9}$ in the ejecta-dominated phase, transitioning to a slower decline as the swept-up CSM becomes dominant.
• SS vs. BS:
  • SS models yield significantly higher DM (tens to hundreds of pc cm$^{-3}$) due to larger ejecta masses.
  • BS models yield much lower DM because binary stripping removes the hydrogen envelope, reducing the total electron column.
  • Mass Dependence: The relationship is non-monotonic. For example, the 11 $M_\odot$ SS model can produce higher DM than the 30 $M_\odot$ SS model at early times due to differences in mean molecular weight ($\mu_e$) and expansion velocity.
• Matching Observations:
  • To match the observed $dDM/dt$ of FRB 20190520B ($\approx -12.4$ pc cm$^{-3}$ yr$^{-1}$), the FRB must occur at $t_{occur} \approx 8\text{--}27$ years post-explosion.
  • Inferred local SNR DM contributions range from ~45 to ~200 pc cm$^{-3}$, depending on the progenitor channel and ionization state.

B. Optical Depth and Escape

• Transparency: The total optical depth is dominated by free-free absorption.
• Escape Time ($t_{esc}$):
  • For SS models, $t_{esc} \approx 13\text{--}17$ years (fiducial ionization).
  • For BS models, the environment is nearly transparent from very early times ($t_{esc} \approx 0$) due to lower densities.
  • This implies that GHz FRB emission can escape the SNR environment within the first few decades, consistent with the young ages inferred for these sources.

C. Rotation Measure (RM)

• Shock-Only Framework: In a model where RM is generated only by shock-amplified magnetic fields in the shocked shell:
  • Only the 11 $M_\odot$ Single-Star model successfully reproduces the observed RM amplitude and decay slope of FRB 20121102.
  • BS models and the 30 $M_\odot$ SS model fail to generate sufficient RM because their lower shocked densities result in weaker ram-pressure amplification of magnetic fields.
• Implication: If the SNR shock is the sole RM source, FRB 20121102 likely originates from a low-mass single-star progenitor. However, the authors acknowledge that a hybrid model (SNR + Magnetar Wind Nebula) is likely required to explain the full complexity of FRB 20121102.

5. Significance and Conclusions

• Cosmological Impact: The local SNR contribution to DM is non-negligible (tens to hundreds of pc cm$^{-3}$). Ignoring this or modeling it incorrectly can bias cosmological inferences regarding the baryon content of the universe ($DM_{IGM}$).
• Progenitor Diagnostics: The secular evolution of DM and RM serves as a diagnostic tool to distinguish between single-star and binary-stripped progenitors. The low DM of FRB 20190520B (if matched to a BS model) would imply a very low ejecta mass, whereas the high DM of other sources suggests SS progenitors.
• Physical Consistency: The study highlights that robust FRB cosmology requires physically consistent local-environment modeling that accounts for progenitor history, CSM structure, and ionization physics, rather than relying on idealized analytic scalings.
• Future Directions: The authors note limitations in 1D symmetry and the lack of self-consistent magnetic field evolution (CR acceleration). Future work will extend to 3D simulations and include Magnetar Wind Nebula (MWN) contributions to better explain the two-stage DM evolution of FRB 20121102.