Understanding the formation and eruption of sigmoidal structure through data-driven modeling of magnetic evolution in solar active region 13500

This study utilizes data-driven magnetofrictional simulations to demonstrate that the November 28, 2023, eruption from solar active region 13500 was triggered by the progressive injection of magnetic helicity and energy, which formed a twisted flux rope that became torus-unstable when the ratio of current-carrying to total relative helicity reached a critical threshold of 0.30.

The Gist

The Solar Storm: A Tale of Twisted Rubber Bands and a Tipping Point

Imagine the Sun not as a ball of fire, but as a giant, chaotic ball of rubber bands. These aren't ordinary rubber bands; they are made of invisible magnetic force. Sometimes, these bands get tangled, twisted, and stretched tight, storing up a massive amount of energy—like a spring being wound up tighter and tighter.

This paper is a detective story about a specific event on November 28, 2023. On that day, a region on the Sun called Active Region 13500 (let's call it "The Storm Spot") suddenly snapped. It unleashed a massive explosion of plasma (a Coronal Mass Ejection, or CME) that hurtled toward Earth, causing a powerful geomagnetic storm.

The scientists wanted to answer two big questions:

1. How did the rubber bands get so twisted in the first place?
2. What was the exact moment they decided to snap?

To find the answers, they didn't just watch the Sun; they built a digital time machine (a computer simulation) to rewind the clock and watch the magnetic fields evolve day by day.

---

1. The Setup: The "S" Shape and the Tug-of-War

Before the explosion, the Storm Spot looked like a giant, glowing "S" (scientists call this a sigmoid). Think of this "S" shape as a sign that the magnetic rubber bands were getting dangerously twisted.

• The Inner Players: In the center of the spot, two magnetic poles (one positive, one negative) were stuck in a slow-motion tug-of-war. They were pulling in opposite directions, stretching the magnetic bands between them.
• The Outer Players: Meanwhile, the outer poles were drifting apart, like two people walking away from each other while holding a rope.

As these poles moved, they didn't just stretch the bands; they twisted them. Imagine taking a rubber band and twisting it with your fingers. The more you twist, the more energy you store.

2. The Simulation: Rewinding the Tape

The scientists used a special computer model called Magnetofriction (MF). Think of this model as a virtual laboratory where they could:

• Inject Energy: They fed the model real data from the Sun's surface, like a chef adding ingredients to a recipe.
• Watch the Evolution: They started the simulation 2.8 days before the explosion.

What happened in the simulation?

• Day 1-2: The magnetic field started as a simple, calm arch (like a rainbow).
• Day 3: As the "tug-of-war" continued, the arch got stretched and sheared.
• The Twist: Eventually, the arch twisted so much it turned into a Flux Rope—a giant, coiled spring of magnetic energy. This coiled spring looked exactly like the "S" shape seen in the real photos.

The simulation was so good that the digital "S" looked almost identical to the real one captured by telescopes. It was like watching a movie where the special effects matched reality perfectly.

3. The Breaking Point: When Does the Spring Snap?

The big mystery was: Why did it explode at that specific time?

The scientists looked at two main "pressure gauges":

A. The Twist Gauge (Kink Instability)

Imagine twisting a rubber band. If you twist it too much, it wants to kink and snap.

• The simulation showed that the core of the magnetic rope twisted about 1.4 times (1.4 turns) before it got unstable.
• This twisting made the rope rise slowly, like a balloon inflating and lifting off the ground. It rose from 50 kilometers high to 80 kilometers high.

B. The "Helicity Ratio" Gauge (The Tipping Point)

This is the most important discovery. The scientists invented a new way to measure the "twistiness" of the storm.

• They calculated a ratio: How much of the magnetic energy is actually doing the twisting vs. how much is just sitting there?
• The Rule: They found that when this ratio hit 0.3 (or 30%), the magnetic structure became unstable enough to trigger a full-blown eruption.
• The Analogy: Think of a seesaw. As long as the "twist" side is light, the seesaw stays balanced. But once the twist side gets heavy enough (reaching that 0.3 mark), the seesaw tips over, and the eruption happens.

4. The Result: A Perfect Match

The simulation predicted that the magnetic rope would rise slowly until it hit a "tipping point" (called Torus Instability).

• At the moment the real Sun exploded, the simulation showed the magnetic rope was at the exact height where the "overhead" magnetic field (the ceiling holding it down) became too weak to hold it.
• The "ceiling" broke, and the rope shot upward at 741 km/s (over 1,600 mph), heading straight for Earth.

Why Does This Matter?

This paper is a breakthrough because it proves that computer models can predict solar storms by watching how the magnetic "rubber bands" twist and stretch.

• The "Helicity Ratio" is the new crystal ball: If we can measure this ratio on the Sun, we might be able to predict exactly when a solar storm is about to happen, giving us better warning to protect our satellites and power grids.
• It's not just about energy: The study showed that even if the total energy goes down (because the sunspot was fading), the twist (helicity) was the real trigger.

In short: The Sun's magnetic field is like a giant, tangled slinky. This study showed us exactly how the slinky gets twisted, how high it bounces before it breaks, and the precise "tipping point" that sends a solar storm racing toward Earth.

Technical Summary

Here is a detailed technical summary of the paper "Understanding the formation and eruption of sigmoidal structure through data-driven modeling of magnetic evolution in solar active region 13500."

1. Problem Statement

Solar flares and Coronal Mass Ejections (CMEs) are driven by the storage and sudden release of magnetic energy in Active Regions (ARs). While the mechanisms for energy storage (flux emergence, shearing motions, sunspot rotation) are understood, predicting the exact onset of eruptions remains a challenge.

• Specific Case: The study focuses on AR 13500, which produced an M9.8 flare and a fast CME on November 28, 2023, triggering a severe geomagnetic storm.
• Scientific Gap: There is a need to understand the specific magnetic evolution leading to the formation of sigmoidal structures (S-shaped twisted fields) and to identify robust non-potential parameters (specifically magnetic helicity ratios) that serve as reliable thresholds for eruption onset, particularly in relation to the Torus Instability.

2. Methodology

The authors employed a data-driven Magnetofriction (MF) simulation to model the 3D magnetic evolution of AR 13500.

• Observational Data:
  • Source: Solar Dynamics Observatory (SDO) Helioseismic and Magnetic Imager (HMI) for photospheric vector magnetic fields.
  • Timeframe: Simulations started 2.8 days (approx. 67 hours) prior to the eruption (Nov 26–28, 2023).
  • Analysis: Vector velocities were derived using the DAVE4VM technique to calculate helicity flux injection ($dH/dt$) and energy injection ($dE/dt$).
• Simulation Model:
  • Code: PENCIL code with a custom driver module.
  • Physics: The model solves a simplified induction equation under the force-free assumption (low-$\beta$ plasma). It evolves the coronal magnetic field quasi-statically in response to photospheric boundary conditions.
  • Energy Constraint: To ensure the simulation injected sufficient energy and helicity to match observations, an ad-hoc parameter ($U$) was introduced to the electric field driver. Two runs were tested ($U=120$ m/s and $U=150$ m/s), with $U=150$ m/s providing the best match to observed Poynting flux.
• Key Metrics Calculated:
  • Magnetic Helicity ($H_V$), decomposed into current-carrying helicity ($H_J$) and mutual helicity ($H_{PJ}$).
  • Free Magnetic Energy ($E_f$).
  • Helicity Ratio: $H_J / H_V$.
  • Torus Instability Criterion: Decay index ($n$) of the overlying magnetic field.

3. Key Contributions

1. Reproduction of Sigmoidal Evolution: The study successfully simulated the transition of an AR from a potential field to a sheared arcade, and finally to a twisted Flux Rope (FR) that morphologically matches the observed S-shaped sigmoid.
2. Validation of Helicity Ratio as a Trigger: The research provides strong evidence linking the ratio of current-carrying helicity to total relative helicity ($H_J/H_V$) with the onset of torus instability.
3. Data-Driven Forecasting: It demonstrates that data-driven MF simulations can accurately reproduce the timing and topology of flux rope formation and slow-rise phases, offering a potential tool for assessing eruptive potential.

4. Key Results

A. Magnetic Evolution and Observations

• Flux Dynamics: The AR exhibited a decay phase with decreasing net magnetic flux, yet simultaneously injected significant magnetic helicity and energy via shearing motions of inner polarities (N2/P2) and proper motions of outer polarities (N1/P1).
• Helicity Injection: The helicity injection rate peaked at $13 \times 10^{37}$Mx$^2$/s at 18:00 UT on Nov 28, coinciding with the CME onset. The accumulated helicity reached$\approx 5 \times 10^{42}$Mx$^2$.
• Simulation Morphology: The simulated magnetic field lines formed a twisted flux rope by the 50th hour. Proxy emission maps (based on $J^2$) showed remarkable morphological agreement with SDO/AIA 304 Å and KSO H$\alpha$ observations of the sigmoid.

B. Flux Rope Dynamics

• Slow Rise: The simulated Flux Rope (FR) formed at $\approx 50$ hours and underwent a slow-rise motion, reaching a height of 80 Mm by the 67th hour (eruption time).
• Twist Evolution: The average twist in the FR core increased from 1.28 turns (at 55h) to 1.43 turns (at 65h), approaching the threshold for ideal kink instability.

C. Non-Potential Parameters and Instability Thresholds

• Free Energy: While total magnetic energy decreased (due to AR decay), free magnetic energy ($E_f$) increased, reaching 57% of the total energy in the $U=150$ m/s run, indicating high eruptive potential.
• Helicity Ratio ($H_J/H_V$):
  • At FR formation (50h): Ratio was 0.13.
  • At eruption onset (67h): Ratio increased to 0.30 (for $U=120$ m/s) and 0.37 (for $U=150$ m/s).
• Torus Instability: The decay index ($n$) of the overlying field exceeded the critical value of $n=1.5$ at heights consistent with the FR core location (approx. 30–34 Mm) at the time of eruption.

5. Significance and Conclusion

• Threshold Confirmation: The study supports the hypothesis that while $H_J/H_V \ge 0.1$ indicates an eruptive tendency, a threshold of $H_J/H_V \ge 0.3$ is associated with the onset of torus instability and the transition to a full eruption.
• Mechanism Insight: The results suggest a two-stage trigger: the accumulation of twist (kink instability) initiates a slow rise, which lifts the flux rope into a region where the overlying field decays rapidly enough (Torus instability) to cause a catastrophic eruption.
• Operational Utility: The data-driven MF approach proves to be a robust tool for reconstructing coronal magnetic fields and evaluating eruptive potential in real-time, particularly by monitoring the evolution of the helicity ratio.

Limitation: The MF model is quasi-static and cannot simulate the fast, dynamic phase of the eruption (acceleration) due to the imposed frictional coefficient, but it successfully captures the pre-eruptive buildup and the critical threshold conditions.