🔭 astrophysics

Hybrid Simulations of Supersonic Shear Flows: II) Cosmic Ray Viscosity

This study uses 2D hybrid simulations to demonstrate that cosmic rays act as long-range messengers introducing a viscosity-like effect in supersonic shear flows, thereby enhancing momentum transfer and influencing energy partitioning and particle acceleration even when their energy density is not dominant.

Original authors: Naixin Liang, Damiano Caprioli

Published 2026-04-15

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Naixin Liang, Damiano Caprioli

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a vast, invisible ocean in space. Unlike the water in our oceans, this "space ocean" is made of super-hot gas (plasma) where particles rarely bump into each other directly. Instead of colliding like billiard balls, they interact through magnetic fields and electric forces.

This paper is about what happens when two layers of this space ocean slide past each other at high speeds. Think of it like a river flowing next to a slower-moving stream, or the wind blowing over a calm lake. This sliding motion is called shear flow.

Here is the story of the paper, broken down into simple concepts:

1. The Problem: The Great Space Slip

When these two layers of space gas slide past each other, they are unstable. It's like trying to balance a stack of cards; eventually, they wobble and collapse. In space, this wobble turns into a chaotic mess of turbulence, creating tiny shockwaves and magnetic storms.

Usually, scientists thought this chaos was just caused by the gas itself. But this paper asks: What if there are "ghosts" in the machine?

2. The Ghosts: Cosmic Rays

Those "ghosts" are Cosmic Rays (CRs). These are high-energy particles (like protons) zooming through space at near-light speed. They are rare, but they are energetic and have a special superpower: they have huge "gyroradii."

The Analogy: Imagine the regular gas particles are like ants marching in a tight line. They can only talk to the ant right next to them.
The Cosmic Rays: These are like giant, fast-moving helicopters flying over the ants. Because they fly so high and fast, a single helicopter can see and touch ants that are miles apart.

3. The Discovery: Cosmic Ray "Viscosity"

The main discovery of this paper is that these "helicopters" (Cosmic Rays) act like a sticky glue or a viscous fluid between the sliding layers.

How it works: Because the Cosmic Rays can fly from the fast-moving layer to the slow-moving layer, they carry momentum with them. They effectively "hand off" speed from the fast layer to the slow layer.
The Result: This speeds up the mixing process. Instead of the layers sliding smoothly for a long time, the Cosmic Rays grab them and force them to mix and stop sliding much faster. The authors call this "Cosmic Ray Viscosity."

4. The Rules of the Game

The researchers ran computer simulations to see how this works under different conditions. They found two golden rules:

Rule #1: Size Matters. For the Cosmic Rays to act as glue, their "flight path" (gyroradius) must be smaller than the size of the sliding layers. If they are too big (like a helicopter flying so high it misses the ants entirely), they can't help mix the layers.
Rule #2: More Ghosts, Faster Mixing. The more Cosmic Rays you have, the faster the layers mix. Even if they aren't the most energetic part of the system, just having a few of them makes the whole process happen much quicker.

5. The Energy Trade-Off

When the layers mix, the energy has to go somewhere. It turns into heat, magnetic storms, or new high-speed particles.

The Accelerator: The Cosmic Rays act like a catalyst. They help the regular gas particles get a "boost" and turn into high-energy particles much faster than they would on their own.
The Catch: There is a trade-off. If the Cosmic Rays are too efficient at mixing the layers, they might stop the sliding motion so quickly that there isn't enough time left to accelerate the particles to their maximum possible speed. It's like a race car driver who hits the brakes too early to avoid a crash; they get to the finish line safely, but they didn't go as fast as they could have.

6. Why This Matters

This isn't just about computer games. This helps us understand real things in the universe:

Black Holes and Jets: The material shooting out of black holes often has these sliding layers.
Galaxy Clusters: The gas between galaxies moves in huge shear flows.
Star Formation: The clouds of gas that birth stars are full of turbulence.

The Big Takeaway:
In the universe, where particles rarely bump into each other, Cosmic Rays act as the invisible hands that mix everything together. They create a kind of "friction" that speeds up the chaos, heats up the gas, and accelerates particles to incredible speeds. Without them, the universe would be a much calmer, slower, and less energetic place.

1. Problem Statement

The paper addresses the dissipation of collisionless, nonrelativistic shear flows in astrophysical environments (e.g., jet boundaries, accretion disks, interstellar medium). While previous work (Paper I) established how shear flows dissipate into heat, magnetic fields, and spontaneously generated nonthermal particles, this study focuses on the dynamical role of pre-existing nonthermal particles (Cosmic Rays, or CRs).

Specifically, the authors investigate:

How CRs with large gyroradii act as "long-range messengers" to mediate momentum exchange between shear layers, effectively inducing a Cosmic Ray (CR) viscosity.
How this CR-induced viscosity alters the dissipation timescales of the shear flow.
The impact of CRs on the partitioning of free kinetic energy among thermal heating, CR reacceleration, and magnetic field amplification.
The efficiency of particle acceleration in the presence of CR seeds compared to spontaneous generation.

2. Methodology

The authors utilize 2D hybrid Particle-in-Cell (PIC) simulations using the dHybridR code.

Physics Model: Ions are treated as kinetic macro-particles evolving under the relativistic Lorentz force, while electrons are modeled as a massless, charge-neutralizing fluid (adiabatic equation of state).
Initial Conditions:
- Flow: A "Kolmogorov flow" (sinusoidal velocity profile $U_x(y) \propto \sin(2\pi y/L)$ ) is initialized, unstable to the Kelvin–Helmholtz (KH) instability.
- CR Seeds: A population of CRs is introduced with specific number densities ( $n_{CR}$ ) and isotropic momenta ( $p_{iso}$ ). These CRs are initially boosted with the local bulk flow.
- Magnetic Field: A uniform background field ( $B_0$ ) is aligned with the flow direction.
Parameters:
- Simulations cover a range of Alfvénic Mach numbers ( $M_A$ ) from subsonic (0.5) to highly supersonic (40).
- CR parameters are varied to test different energy densities ( $\eta_{CR}$ ) and gyroradii relative to the shear scale length ( $L$ ).
- Box size is typically $L = 200 d_i$ (ion inertial lengths), with a specific "Run L" using $L=2000 d_i$ to test scale effects.
Metrics:
- Dissipation Timescales: Defined by the decay of the momentum flux fraction $\Delta(t)$ . Key metrics include onset time ( $\tau_{90}$ ), halving time ( $\tau_{50}$ ), and viscous timescale ( $\tau_\nu$ ).
- Energy Partitioning: Tracking the fraction of energy in thermal gas, CRs, and magnetic fields.
- Acceleration: Monitoring the development of nonthermal tails and the weighted maximum energy ( $\bar{E}_{max}$ ).

3. Key Contributions

Quantification of CR Viscosity: The paper provides kinetic evidence that CRs induce an effective viscosity even when they are not energetically dominant, provided their gyroradii are smaller than the shear scale.
Trade-off Mechanism: It identifies a critical trade-off between acceleration rate and shear lifetime. While CRs accelerate the dissipation of the shear (viscosity), this rapid dissipation can truncate the time available for particle acceleration, potentially limiting the maximum energy achievable.
Kinetic vs. Fluid Dynamics: The study demonstrates that fluid/MHD approaches cannot capture the partitioning of energy into nonthermal particles or the specific microphysical momentum exchange rates driven by wave-particle interactions in collisionless plasmas.
Hillas Limit Constraints: The simulations show that in undriven, decaying shear flows, the maximum particle energy is often limited by the rapid depletion of the shear before the Hillas limit is reached, unless CR seeds are optimized.

4. Key Results

A. CR Viscosity and Dissipation Timescales

Enhanced Dissipation: Increasing CR number density ( $n_{CR}$ ) or energy density ( $\eta_{CR}$ ) significantly reduces the shear halving time ( $\tau_{50}$ ).
Gyroradius Threshold: The viscosity effect is most efficient when the CR gyroradius ( $r_g$ ) is smaller than or comparable to the shear scale ( $L$ ). If $r_g \gg L$ , CRs decouple from the shear, and the viscosity effect diminishes.
Thresholds: Significant acceleration of dissipation occurs when $n_{CR} \gtrsim 1\%$ and $\eta_{CR} \lesssim 0.1$ .

B. Particle Acceleration

Thermal Ion Acceleration: Pre-existing CR seeds significantly speed up the acceleration of thermal ions. In benchmark runs, the maximum energy of thermal ions was 4 times higher with CR seeds than without.
Spectral Evolution: The presence of CRs leads to faster development of nonthermal power-law tails (approx. $E^{-2.2}$ ).
The "Too Fast" Paradox: If CR density is too high, the shear dissipates so quickly that the acceleration process is cut short. Consequently, runs with very high $n_{CR}$ sometimes show less efficient acceleration than moderate runs because the "fuel" (the shear) is exhausted before particles can reach high energies.
CR Reacceleration: Pre-existing CRs gain energy (peak shifts $\sim 2\times$ ), but they do not develop extended power-law tails in these undriven, finite-box simulations.

C. Energy Partitioning

Direction of Energy Flow: The flow of energy between thermal gas and CRs is bidirectional and depends on the initial energy dominance:
- If CRs are initially dominant ( $\epsilon_{CR} \gg \epsilon_{gas}$ ), energy flows from CRs to the thermal gas (heating the gas).
- If the thermal gas is dominant, energy flows from the shear to both CRs and the thermal gas.
Magnetic Amplification: Magnetic turbulence is invariably generated, reaching $\sim 1\%$ of the initial free energy. While small compared to equipartition, these fluctuations are crucial for coupling the species.
Mach Number Dependence:
- Subsonic/Transonic: Shear acts primarily as a trigger for fluctuations that redistribute energy from CRs to thermal gas.
- Supersonic ( $M_A \gtrsim 10$ ): The shear acts as a direct energy source for both thermal heating and nonthermal acceleration. The fraction of energy in nonthermal particles saturates at $\sim 40\%$ for high $M_A$ .

D. Maximum Energy Limits

The weighted maximum energy ( $\bar{E}_{max}$ ) generally scales with CR parameters but is capped by the Hillas limit ( $E_H \approx 40 E_0$ for the benchmark).
In many cases, $\bar{E}_{max}$ falls below the Hillas limit because the CR viscosity depletes the shear faster than particles can be accelerated to the theoretical maximum.

5. Significance

This work establishes that CR viscosity is an intrinsic property of collisionless astrophysical plasmas, not just a fluid approximation.

Astrophysical Implications: The findings suggest that in environments like jet boundaries, accretion disks, and the interstellar medium, pre-existing nonthermal particles can drastically alter the evolution of shear flows, accelerating their decay and modifying how kinetic energy is converted into heat and magnetic fields.
Modeling Impact: The results highlight the necessity of kinetic simulations for accurately modeling energy partitioning in these systems, as fluid models fail to capture the nonthermal acceleration channels and the specific feedback loops between CRs and shear dissipation.
Future Directions: The authors note that these results apply to decaying shears. In driven systems (continuous energy injection), they anticipate that power-law tails could extend closer to the Hillas limit and magnetic turbulence could reach equipartition, a topic reserved for future study.