Numerical simulations of waves and turbulence in coronal loops: observables and spectra

This study uses numerical simulations of a coronal loop to demonstrate that the upcoming Multi-slit Solar Explorer (MUSE) mission can detect signatures of phase mixing and turbulent cascades through synthesized high-resolution spectroscopic observations, specifically showing that intensity power spectra at MUSE resolution accurately reflect the underlying density turbulence spectrum.

The Gist

Imagine the Sun's outer atmosphere, the corona, as a giant, glowing forest of magnetic "trees" called coronal loops. For decades, scientists have been puzzled by a mystery: these loops are incredibly hot, but the energy source keeping them that way is hard to pin down. It's like trying to figure out how a campfire stays burning when you can't see the wood being added.

This paper is a computer simulation study that tries to solve that mystery by watching how "waves" and "turbulence" move inside these magnetic loops. The researchers are essentially building a digital twin of a solar loop to see if they can spot the heat-generating mechanisms before the next big space telescope, called MUSE, launches.

Here is the breakdown of their experiment using simple analogies:

1. The Setup: A Twisted Garden Hose

The researchers created a virtual, cylindrical magnetic tube (the loop) filled with hot plasma (superheated gas).

• The Environment: The inside of the tube is denser (thicker) than the outside, creating a boundary layer.
• The Disturbance: They didn't just shake the tube; they injected two types of "wiggles" into it:
  1. The Torsional Wave: Imagine twisting a garden hose back and forth. This is a smooth, spiraling motion.
  2. The Turbulent Component: Imagine shaking the hose randomly and chaotically, like a stormy day.
• The Mix: They ran simulations with different ratios of these two wiggles, from mostly smooth twisting to mostly chaotic shaking.

2. The Process: Mixing and Breaking

As these waves travel, two main things happen that generate heat:

• Phase Mixing (The "Traffic Jam"): Because the inside of the loop is denser than the outside, the waves travel at different speeds. Imagine a line of runners where the ones on the inside lane run slower than the ones on the outside. Eventually, the line gets stretched and twisted into a mess. This stretching creates tiny, fine-scale ripples. In physics, these tiny ripples are where energy turns into heat.
• Turbulent Cascade (The "Domino Effect"): The chaotic shaking creates a cascade. Big, slow waves crash into each other and break down into smaller, faster waves, which break down into even tinier ones, until the energy is finally dissipated as heat.

The paper found that these two processes often work together. The "traffic jam" (phase mixing) helps create the conditions for the "domino effect" (turbulence) to happen faster, heating the plasma more efficiently than either could alone.

3. The Observation: The "MUSE" Camera

The researchers didn't just look at the invisible physics; they simulated what a future telescope, MUSE (Multi-slit Solar Explorer), would actually see. MUSE is like a super-powered camera that can take incredibly sharp pictures of the Sun's light and color.

They synthesized three specific "images" from their simulation:

• Brightness (Intensity): How bright the loop looks. They saw that as the waves move, the loop starts to look like it has thin, parallel threads or strands, rather than being a smooth cylinder.
• Color Shift (Doppler Velocity): This shows how fast the gas is moving toward or away from the camera. They saw distinct patterns of motion, especially near the edges of the loop where the "traffic jam" (phase mixing) is strongest.
• Blur (Non-thermal Broadening): This measures how "fuzzy" the light is due to random motion. They found this blur was strongest at the loop's boundaries, confirming that the chaotic mixing is happening there.

4. The Verdict: Can We See It?

The most important conclusion is about resolution.

• The researchers compared their "perfect" high-resolution simulation with a "blurred" version that mimics what MUSE will see.
• The Good News: Even with the "blur" of the telescope, MUSE will still be able to see the main patterns. It can detect the formation of those thin threads and the specific signatures of the waves and turbulence.
• The Data: They analyzed the "texture" of the images (using something called power spectra). They found that the texture of the brightness images (what MUSE sees) matches the texture of the actual density inside the loop. This means that by looking at the brightness patterns MUSE captures, scientists can actually infer how the density and energy are distributed inside the loop, even though they can't see inside it directly.

Summary

In short, this paper says: "We built a digital solar loop and shook it with waves and turbulence. We found that these motions create tiny, heat-generating ripples. We then simulated what the upcoming MUSE telescope would see, and we are confident that MUSE is powerful enough to spot these patterns. If MUSE sees these specific 'threads' and 'blurs' in the Sun's light, it will confirm that waves and turbulence are indeed the engines heating the solar corona."

Technical Summary

Problem Statement
The heating mechanism of the solar corona and the complex dynamics resulting from plasma turbulence and wave propagation remain unresolved. While waves and oscillations in coronal magnetic structures are widely observed, the specific interplay between phase mixing and turbulent cascades in generating small scales—and the subsequent dissipation of energy—requires further investigation. A critical challenge is determining whether upcoming high-resolution spectroscopic missions, specifically the Multi-slit Solar Explorer (MUSE), can resolve these phenomena. Current observational capabilities often lack the spatial and temporal resolution to distinguish between the spectral signatures of density fluctuations, velocity fields, and the resulting emission line properties (intensity, Doppler shift, and non-thermal broadening) in the context of these dynamic processes.

Methodology
The authors employ direct numerical simulations (DNS) using the COHMPA (COmpressible Hall Magnetohydrodynamics simulator for Plasma Astrophysics) code to solve 3D compressible, viscous-resistive magnetohydrodynamics (MHD) equations. The model represents a coronal loop as a cylindrical, pressure-balanced magnetic flux tube with transverse inhomogeneities in density and magnetic field.

• Initial Conditions: The system is initialized with a superposition of two components: a torsional Alfvén wave (azimuthally polarized, non-compressive) and a turbulent transverse perturbation (random superposition of Fourier modes with a power-law spectrum and near-zero cross-helicity). A parameter $\alpha$ ($0 \le \alpha \le 1$) controls the relative weight of the turbulent component versus the wave component.
• Physical Parameters: The simulations utilize a low plasma beta ($\beta = 0.01$) typical of the corona. The domain aspect ratio is set to $\Lambda = 4$ to ensure that phase mixing and nonlinear cascades operate with comparable efficiencies, a regime where their synergistic effects are most pronounced.
• Synthetic Observables: To bridge the gap between simulation and observation, the authors use the FoMo code to synthesize spectral moments for the Fe IX 171 Å line. They calculate the emission intensity ($I_0$), Doppler shift ($I_1$), and non-thermal broadening ($I_2$). These maps are generated at both the full simulation resolution and a degraded resolution of $312 \text{ km} \times 312 \text{ km}$, matching the nominal spatial resolution of the MUSE instrument.
• Analysis: The study tracks the time evolution of dissipation rates and compares 1D power spectra of the synthetic observables ($I_0, I_1$) against the power spectra of the underlying physical fields (density $\rho$ and line-of-sight velocity $v_x$).

Key Contributions and Results

1. Synergistic Dissipation: The simulations confirm that phase mixing and turbulent cascades work synergistically. In regimes where phase mixing initially reduces spatial scales, it facilitates the activation of the nonlinear cascade. Consequently, the total dissipation time is shorter than that of either mechanism acting alone. The dissipation rate increases with the turbulent amplitude ($\alpha$) and the perturbation amplitude ($B_1$).
2. Structural Evolution:
   • For pure torsional waves ($\alpha=0$), the loop retains its cylindrical symmetry, and density remains unchanged.
   • When turbulence is present ($\alpha > 0$), the loop structure deforms, and the background density and temperature are advected, leading to filamentation. This deformation is more pronounced as $\alpha$ increases.
   • Small-scale fluctuations are generated primarily in the inhomogeneous boundary region of the loop. Phase mixing creates transverse small scales localized at the loop boundary, visible as thin stripes in the Doppler shift and non-thermal broadening maps.
3. Observability with MUSE:
   • The study demonstrates that MUSE's spatial resolution ($312 \text{ km}$) is sufficient to capture the large-scale patterns of phase mixing and the formation of longitudinal threads, even if the finest scales are not fully resolved.
   • The temporal resolution of MUSE ($1-4$ s) is finer than the simulation's unit time ($\sim 8.3$ s), suggesting that MUSE could track the detailed time evolution of these wave and turbulence dynamics.
4. Spectral Analysis:
   • Intensity ($I_0$): The power spectrum of the emission intensity $I_0$ closely matches the power spectrum of the density field $\rho$ at large spatial scales (300 km to 3000 km). The spectral indices are very similar, indicating that $I_0$ spectra can reliably infer the distribution of density fluctuations.
   • Doppler Shift ($I_1$): The power spectrum of the Doppler shift $I_1$ is systematically steeper than the spectrum of the line-of-sight velocity $v_x$. This steepening is attributed to the line-of-sight integration, which smooths fluctuations. Despite this difference, $I_1$ spectra still follow power laws and can provide information on velocity fluctuations, provided the spectral index difference is accounted for.

Significance
The paper establishes a forward-modeling framework to interpret future MUSE observations of coronal loops. It claims that the unprecedented resolution of MUSE will allow for the first time the computation of power spectra of coronal loop observables with sufficient fidelity to infer physical properties of the plasma. Specifically, the agreement between the spectral index of the intensity power spectrum and the density power spectrum suggests that $I_0$ maps can be used to diagnose the spectrum of density fluctuations inside a loop. Furthermore, the study validates that the combined effects of phase mixing and turbulence, which are difficult to isolate observationally, produce distinct signatures in synthetic observables that MUSE should be capable of detecting. The results provide a theoretical basis for using high-resolution spectroscopic data to constrain models of coronal heating and fluctuation dynamics.