Particle-in-Cell Simulation of the Parametric Decay Instability of Alfvén Waves with Absorbing Boundary Conditions

This paper presents fully kinetic one-dimensional simulations of Alfvén wave parametric decay instability using absorbing boundary conditions, revealing that at low plasma beta, nearly 92% of pump wave energy transfers to a backward-propagating Alfvén wave while the remainder heats electrons and ions only after the instability sufficiently develops, with heating rates approximately twice the linear growth rate.

The Gist

The Big Picture: A Cosmic Energy Transfer

Imagine a giant, invisible ocean made of charged particles (plasma) that fills space, stars, and fusion reactors. In this ocean, waves travel just like ripples on a pond. These are called Alfvén waves.

The scientists in this paper wanted to understand what happens when a big, powerful wave (the "pump") crashes into the plasma. Specifically, they were looking at a phenomenon called Parametric Decay Instability (PDI).

Think of PDI like a large, heavy drumstick hitting a drum. Instead of just making a sound, the energy from that single hit splits. The big wave breaks apart into two smaller things:

1. A smaller wave traveling in the opposite direction (like a reflection).
2. A "sound wave" traveling in the same direction (like a compression in the air).

The Experiment: A Controlled "Open Window"

Most previous studies on this topic were like studying a drum in a sealed, echoey room. The waves would bounce off the walls, hit the drum again, and create a confusing mess of energy that didn't look like the real world.

The researchers in this paper built a simulation with absorbing boundaries.

• The Analogy: Imagine the simulation room has walls made of special "black hole" foam. When a wave hits the wall, it disappears completely instead of bouncing back.
• Why it matters: This lets them see exactly how much energy is transferred to the particles (electrons and ions) without the "echoes" messing up the math. It's like listening to a single drum hit in a soundproof booth to hear exactly how the drumhead vibrates.

They also used a fully kinetic approach.

• The Analogy: Previous studies often treated the tiny electrons like a smooth, invisible fluid (like water). This study treated every single electron and ion like a distinct, bouncy ball. This is important because, in reality, these tiny balls can bounce around and heat up in ways a smooth fluid cannot.

The Results: Where Did the Energy Go?

The researchers pumped energy into the system and watched where it went. Here is the breakdown of the "energy pie":

• 92% went to the backward wave: The vast majority of the energy simply turned into the smaller wave traveling the other way. It was like the drumstick hitting the drum and mostly just sending a shockwave back up the stick.
• 6-7% went to the ions (heavy particles): The heavy particles (ions) got a little bit of heat.
• 1-2% went to the electrons (light particles): The tiny electrons got a very small amount of heat.

Key Finding: The heating didn't happen immediately. It was like a "slow burn." The instability had to grow strong enough first before the particles started getting hot. Once the instability kicked in, the particles heated up at a rate that was roughly twice as fast as the instability itself grew.

Why the Difference in Heating?

The paper explains why the heavy ions got more heat than the light electrons:

• The Ions: The "sound wave" created by the split became a bit "steep" (like a sharp cliff). The heavy ions crashed into this steep wave and got pushed, gaining energy.
• The Electrons: The electrons are so light and fast that they mostly just swam through the wave without getting caught. They didn't get "trapped" by the wave in the same way the ions did, so they stayed relatively cool.

The Takeaway

This study is a "baseline" test. It proves that if you look at a simple, one-dimensional line of plasma with realistic boundaries, you can accurately measure how energy splits between waves and particles.

The authors conclude that while this specific setup (a straight line) shows very little heating for electrons, it sets the stage for future, more complex 3D simulations. In those more realistic 3D worlds, they expect the electrons might get much hotter, which could change how we understand heating in fusion reactors and the solar wind.

In short: They built a perfect, echo-free digital lab to watch a big plasma wave break apart. They found that most of the energy just bounced back as a smaller wave, while a small fraction heated up the heavy particles, and a tiny fraction warmed the light ones.

Technical Summary

Technical Summary: Particle-in-Cell Simulation of the Parametric Decay Instability of Alfvén Waves with Absorbing Boundary Conditions

Problem Statement
The parametric decay instability (PDI) of Alfvén waves (AWs) is a fundamental mechanism for energy transfer, plasma heating, and turbulence generation in space, astrophysical, and fusion plasmas. While extensive theoretical and numerical studies exist, most prior simulations of Alfvén wave PDI have relied on hybrid models (kinetic ions, fluid electrons) and periodic boundary conditions. Periodic boundaries can lead to an overestimation of energy partitioning due to repeated wave-plasma interactions and may misrepresent instability signatures compared to open systems. Furthermore, the kinetic mechanisms governing the partitioning of pump-wave energy between ions and electrons, particularly in open systems, remain poorly understood. There is a specific need to determine the fraction of pump energy converted into internal energy and how it is distributed between species under realistic, non-periodic conditions.

Methodology
This study presents fully kinetic, one-dimensional (1D) Particle-in-Cell (PIC) simulations of the PDI of a left-hand circularly polarized (LHCP) Alfvén wave using the open-source SMILEI code. The simulations are conducted under conditions relevant to the Large Plasma Device (LAPD), specifically at low plasma beta ($\beta = 5 \times 10^{-4}$) and a normalized wave amplitude of $\delta B/B_0 = 0.01$.

Key methodological features include:

• Absorbing Boundary Conditions: Unlike previous studies using periodic domains, this work employs field masks at the domain boundaries to absorb outgoing waves. This prevents repeated wave-plasma interactions, offering a more realistic representation of an open system.
• Wave Injection: A pump AW is injected via an antenna module prescribing a current density profile, rather than simply prescribing field values, to better mimic laboratory conditions.
• Boundary Handling: While waves are absorbed at the boundaries, particles are confined by virtual reflecting walls placed within the physical region to prevent particle intervention from the mask zones, ensuring accurate measurement of heating.
• Parameters: The simulation uses a reduced speed of light and mass ratio ($m_i/m_e = 100$) to manage computational cost, with values verified as sufficient to capture the relevant physics. The domain is normalized to the ion inertial length ($d_i$).

Key Results
The simulations quantify the energy partitioning and heating dynamics as follows:

• Energy Partition: Approximately 92% of the pump wave energy is transferred to the backward-propagating "child" Alfvén wave. The remaining energy is partitioned between ions ($\sim 6-7\%$) and electrons ($\sim 1-2\%$).
• Heating Dynamics: Ion and electron heating occur only after the PDI has sufficiently developed. The heating is predominantly parallel to the ambient magnetic field.
• Growth Rates: The growth rates for ion and electron energy ($\gamma_i$ and $\gamma_e$) are found to be approximately twice the linear PDI growth rate ($\gamma \approx 0.0125\Omega_{ci}$, while $\gamma_{i,e} \approx 0.024\Omega_{ci}$). This quadratic relationship arises because particle energy depends on the square of the fluctuation amplitude, whereas the instability growth follows the amplitude linearly.
• Mechanisms:
  • Ions: Heating is attributed to ion Landau damping of the ion-acoustic wave and a small degree of front steepening in these waves. Phase-space analysis reveals a distinct "bump" in the ion distribution function near the ion-acoustic phase velocity.
  • Electrons: Electron heating is weaker and not attributed to Landau damping, as the ion-acoustic phase velocity lies within the bulk of the electron distribution. Instead, heating is likely driven by the steepening of the acoustic waves, which accelerates particles.
• Parameter Sensitivity: A parametric scan reveals that particle thermal energies and the PDI growth rate increase with pump amplitude ($\delta B/B_0$), decrease with plasma beta ($\beta$), and initially increase with the temperature ratio ($T_e/T_i$) before saturating. The numerical growth rates show qualitative agreement with theoretical predictions across these parameters.

Significance and Claims
The authors claim this work establishes the first fully kinetic study of Alfvén wave PDI utilizing absorbing wave boundary conditions. By demonstrating that the 1D limit ($k_\perp = 0$) yields specific energy partitioning and heating rates, the study provides a necessary baseline for future extensions to finite $k_\perp$ in higher dimensions. The paper posits that while electron heating is weak in this 1D configuration, future 2D/3D simulations may induce stronger electron heating (potentially via electron Landau damping of AWs), which could significantly alter PDI dynamics. The results offer critical insights into the energy transfer mechanisms in open plasma systems, relevant to understanding solar wind heating and fusion plasma transport.