Formation and rising phase of a flux rope through data-constrained simulations

This study utilizes a data-constrained, resistive magnetohydrodynamic simulation to demonstrate that a non-force-free initial magnetic field, characterized by a Lorentz force imbalance, can spontaneously trigger the formation and subsequent eruption of a flux rope in active region NOAA 12241 without requiring pre-existing structures or photospheric driving motions.

The Gist

Imagine the Sun's surface as a giant, chaotic trampoline made of invisible magnetic rubber bands. Sometimes, these bands get twisted, stretched, and tangled so tightly that they suddenly snap, launching a massive explosion of energy and matter into space. This is a solar eruption, and it can mess up our satellites, GPS, and power grids here on Earth.

Scientists have been trying to predict these explosions for years, but it's like trying to predict a thunderstorm by only looking at the ground. Usually, they try to guess what the "air" (the magnetic field) looks like high up in the sky based on what they see on the ground. But this guess often assumes the air is perfectly calm and balanced, which isn't true right before a storm.

This paper is about a new, more realistic way to simulate these storms. Here is the story of what they found, explained simply:

1. The Setup: A Tangled Mess, Not a Calm Lake

Most computer models start by assuming the Sun's magnetic field is in a perfect, peaceful balance (like a calm lake). But in reality, right before an explosion, the magnetic field is already wobbly and out of balance.

The researchers decided to start their simulation with a messy, unbalanced magnetic field. They took a real snapshot of a specific sunspot (called NOAA 12241) just minutes before a real explosion happened in 2014. Instead of smoothing out the wrinkles, they fed the "messy" data directly into their supercomputer.

The Analogy: Imagine you are trying to predict when a rubber band will snap.

• Old way: You assume the rubber band is sitting perfectly still on a table, then you slowly pull it until it snaps.
• This paper's way: You look at a rubber band that is already being twisted and pulled by invisible hands. You don't wait to pull it; you just let go and see what happens.

2. The Spark: The "Push" That Starts It All

Because they started with this messy, unbalanced field, there was an immediate "push" (called the Lorentz force) right from the very first second.

In the lower atmosphere of the Sun (the chromosphere), this push squeezed the gas together. Think of it like someone stomping on a soda can. The gas got squished, heated up instantly, and then exploded upward. This heat caused the cooler gas from the bottom to boil and shoot up into the hot upper atmosphere, a process called evaporation.

The Analogy: It's like a pressure cooker. The unbalanced magnetic force is the heat turning up the burner. The gas inside (plasma) gets squeezed, gets super hot, and turns into steam that rushes upward, filling a balloon.

3. The Balloon: A Magnetic Rope Takes Flight

As this hot, heavy gas rushed up, it got caught in the twisted magnetic bands. These bands wrapped around the gas, forming a giant, coiled structure called a flux rope. You can think of this as a giant, magnetic spring loaded with heavy, hot gas.

Once this "magnetic spring" was fully formed and loaded with gas, it didn't just sit there. The magnetic forces were so strong that they overpowered gravity. The spring started to rise.

The Analogy: Imagine a coiled spring (the flux rope) filled with water (the gas). If you squeeze the bottom of the spring, it shoots up. In this case, the "squeeze" was the magnetic imbalance, and the "spring" shot up into the sky at about 350 kilometers per second (that's nearly 800,000 miles per hour!).

4. The Detour: Why It Didn't Go Straight Up

You might expect the explosion to go straight up, like a rocket. But this one didn't. As it rose, it got pushed sideways.

Why? Because the magnetic "pressure" in the sky wasn't the same everywhere. Some areas had strong magnetic walls holding things back, while others had "open doors" with weak magnetic pressure. The rising rope naturally drifted toward the path of least resistance, like a balloon floating toward a gap in a crowd.

The Analogy: Imagine a hot air balloon rising in a canyon. If the canyon walls are high and tight on the left, but open and wide on the right, the balloon will drift to the right. The solar eruption did the same thing, drifting toward a region where the magnetic "walls" were weaker.

5. The Result: A Perfect Match

The researchers watched their simulation run for about 16 minutes. They saw the rope form, load up with gas, shoot up, and drift sideways. When they compared their computer movie to the actual video of the 2014 solar eruption taken by NASA satellites, it matched almost perfectly.

• Speed: The simulated rope moved at the same speed as the real one.
• Path: It drifted in the same direction.
• Cause: It proved that you don't need to manually "push" the Sun to make it erupt. If you just start with a realistic, messy magnetic field, the explosion happens all by itself.

Why This Matters

This study is a big deal because it shows that we don't need to guess how solar storms start. By using a more realistic starting point (acknowledging that the Sun is never perfectly calm), we can simulate eruptions that look exactly like the real thing.

It's like finally figuring out that you don't need to push a car to make it roll; if you just put it on a hill with the right slope, gravity does the rest. This helps scientists get better at predicting space weather, keeping our technology safe from the Sun's temper tantrums.

Technical Summary

1. Problem Statement

Solar eruptions (flares and Coronal Mass Ejections) are driven by the release of magnetic energy, yet the precise initiation mechanisms remain unclear due to the lack of direct magnetic field observations in the solar corona.

• Limitations of Current Models: Most data-constrained and data-driven simulations rely on Non-Linear Force-Free Field (NLFFF) extrapolations. These assume the coronal magnetic field is force-free (Lorentz force $\mathbf{J} \times \mathbf{B} = 0$), which is a valid approximation only in the low-$\beta$ corona. However, the photosphere and chromosphere are high-$\beta$ regions where the Lorentz force is non-negligible.
• The Gap: By forcing a force-free assumption at the lower boundary, models often neglect the inherent Lorentz force imbalances present in the real atmosphere. These imbalances may be crucial for triggering eruptions without requiring artificial external drivers (like shearing motions or flux emergence) to destabilize the system.
• Objective: The study aims to reproduce the solar eruption SOL2014-12-18T21:41 (associated with NOAA AR 12241) using a Non-Force-Free Field (NFFF) extrapolation to investigate if the inherent Lorentz force imbalance alone can trigger flux rope formation and eruption in a realistic, stratified atmosphere.

2. Methodology

The authors employed a fully resistive, compressible Magnetohydrodynamic (MHD) simulation framework.

• Initial Conditions (NFFF Extrapolation):
  • Instead of NLFFF, they used an NFFF extrapolation based on the principle of minimum energy dissipation rate.
  • The magnetic field is modeled as a superposition of two linear force-free fields and one potential field ($\mathbf{B} = \mathbf{B}_1 + \mathbf{B}_2 + \mathbf{B}_3$).
  • This approach preserves the non-zero Lorentz force observed in the photospheric vector magnetogram (SDO/HMI) taken minutes before the flare onset.
• Atmospheric Model:
  • The simulation domain includes a stratified atmosphere spanning from the photosphere to the corona.
  • It features a dense, cold lower layer (chromosphere-like), a narrow transition region, and a hot, extended corona.
  • Physics included: Thermal conduction, optically thin radiative cooling, and gravity.
• Numerical Setup:
  • Code: PLUTO (Mignone et al. 2007).
  • Domain: Extended southward to match the eruption direction, with a resolution of 1 arcsec per pixel.
  • Tracking: The GUITAR algorithm was used to identify and track the magnetic flux rope based on the twist number ($Tw$).

3. Key Contributions

• NFFF-Driven Eruption: Demonstrated that a flux rope can form and erupt self-consistently solely due to the initial Lorentz force imbalance derived from observational data, without assuming pre-existing flux ropes or applying artificial photospheric driving motions.
• Stratified Atmosphere Integration: Successfully coupled NFFF extrapolation with a realistic, compressible, stratified atmosphere, allowing for the study of mass loading (evaporation) and thermodynamic evolution during the eruption.
• Kinematic Analysis: Provided a detailed kinematic analysis of the flux rope's rise, deflection, and acceleration, comparing simulation results directly with observational data.

4. Results

The simulation successfully reproduced the eruption dynamics in two distinct phases:

A. Flux Rope Formation and Mass Loading (0 – ~6 min)

• Trigger: The non-zero Lorentz force in the initial state (specifically converging horizontal components near the polarity inversion line) induced a rapid converging flow.
• Compression & Heating: This flow caused localized compression in the upper transition region, increasing temperature by a factor of 8 (up to $4 \times 10^6$ K).
• Evaporation: Thermal conduction heated the plasma from above, triggering chromospheric evaporation. Dense, hot material filled the magnetic arcade, creating a mass-loaded flux rope.
• Reconnection: Opposite legs of the arcade underwent reconnection, increasing the twist and solidifying the flux rope structure.

B. Rise and Escape (~6 – 16 min)

• Rising Phase: The flux rope began to rise slowly, driven by a net upward force (combination of Lorentz force and gas pressure gradient).
• Deflection: As it ascended, the flux rope was deflected toward the southeast. This deflection was caused by the structure moving toward regions of low magnetic pressure (where magnetic reconnection at the leading edge reduced magnetic pressure and increased plasma $\beta$).
• Acceleration: After the initial slow rise, the flux rope entered a ballistic phase with constant acceleration of $424 \text{ m s}^{-2}$ (1.5 times solar gravity).
• Escape: The flux rope exited the simulation domain at a speed of approximately 350 km s$^{-1}$ after 16 minutes.
• Energy Evolution: The free magnetic energy initially decreased due to ohmic dissipation (triggering the eruption) and then slowly increased as plasma motions strengthened the arcades, before dropping as the rope left the domain.

5. Significance and Conclusions

• Validation of NFFF: The study confirms that NFFF extrapolations provide a more realistic initial state for eruption modeling than NLFFF, as they naturally include the Lorentz force imbalances necessary to trigger eruptions without external drivers.
• Comparison with Previous Work: Unlike previous incompressible/isothermal models of the same event (Prasad et al. 2023) where the eruption stalled, the inclusion of compressibility and a realistic atmosphere allowed the flux rope to escape the domain with high acceleration, matching observational speeds.
• Atmospheric Sensitivity: The authors noted that the choice of the transition region profile significantly impacts the dynamics. A smoother, higher transition region delayed the eruption and resulted in a slower, denser, and colder flux rope, highlighting the need for precise atmospheric modeling to reproduce synthetic emissions.
• Future Outlook: This framework establishes a robust numerical basis for studying flaring active regions. Future work will focus on generating synthetic emissions (to compare directly with SDO/AIA observations) and incorporating test-particle simulations to estimate non-thermal emissions.

In summary, this paper demonstrates that initial Lorentz force imbalances in a realistic, stratified atmosphere are sufficient to trigger the formation, mass loading, and rapid escape of a solar flux rope, offering a compelling alternative to models that rely on pre-existing structures or artificial driving.