Knudsen number as a non-thermal parameter: possible origin of skewness in space plasma distributions

This paper theoretically investigates the origin of skewness in solar wind electron distributions by introducing a Skew-Kappa model and a Krook-like collision term into the Boltzmann equation, ultimately deriving a relationship where the skewness parameter is proportional to the effective Knudsen number ($\delta \sim K_N$).

The Gist

The Big Picture: Why Space Plasma is Weird

Imagine the solar wind not as a smooth, calm river, but as a chaotic crowd of people (electrons) running through a giant stadium. In a perfect, calm world (what physicists call "thermal equilibrium"), everyone would be running at roughly the same speed, forming a neat, bell-shaped curve.

But in space, things are messy. The electrons don't run in a neat bell curve. Instead, they have two weird features:

1. The Long Tail: A few super-fast runners are way ahead of the pack, stretching the curve out.
2. The Lopsidedness (Skewness): The crowd isn't just running fast; they are leaning heavily to one side. More electrons are running in one direction (away from the Sun) than the other, making the distribution look like a lopsided hill rather than a symmetrical mountain.

This paper asks: Why is the crowd lopsided?

The Main Discovery: The "Knudsen Number" is the Secret Sauce

The authors, Gallo-Méndez, Viñas, and Moya, propose a new way to understand this lopsidedness. They suggest that the "tilt" of the electron crowd is directly caused by how crowded the stadium is and how steep the "slopes" are that the electrons are running down.

They introduce a concept called the Knudsen number.

• The Analogy: Imagine the electrons are hikers.
  • If the hikers are in a dense forest where they bump into trees constantly (high collision), they move in a straight, predictable line.
  • If the hikers are in a wide-open desert where they rarely bump into anything (low collision), they can drift wildly.
  • The Knudsen number measures how "open" the desert is versus how "dense" the forest is.

The paper finds a direct link: The more "open" the space is (higher Knudsen number), the more the electron crowd leans over (higher skewness).

How They Figured It Out

The scientists used a mathematical tool called the Boltzmann Equation. Think of this equation as a giant rulebook that predicts how a crowd of particles moves and interacts.

1. The "Krook" Rule: To make the math work, they added a specific rule to the equation called a "Krook-like term." Imagine this as a "reset button." It represents the rare moments when electrons bump into each other and try to straighten out their paths.
2. The Skew-Kappa Distribution: They assumed the electrons follow a specific shape called a "Skew-Kappa" distribution. This is a fancy mathematical shape that allows for both the "long tail" of fast electrons and the "lopsided" tilt.
3. The Calculation: They ran the numbers to see what happens when you combine the "reset button" (collisions) with the "slopes" (changes in temperature and density).

The Result: A Simple Formula

After doing the heavy math, they found a surprisingly simple relationship. The amount the electron crowd leans (the skewness parameter, $\delta$) is proportional to the Knudsen number.

• In plain English: The steeper the temperature or density changes in the solar wind, and the fewer collisions the electrons have, the more the distribution tilts.
• The "Aha!" Moment: They found that this tilt is essentially a measurement of how far the plasma is from being calm and balanced. It's like a "stress gauge" for the solar wind.

Why This Matters (According to the Paper)

The paper doesn't claim this will fix satellites or predict space weather storms immediately. Instead, it provides a theoretical bridge:

1. Connecting the Dots: It connects three things that scientists have studied separately:
   • The shape of the electron crowd (Skewness).
   • How heat moves through space (Heat Flux).
   • How often particles bump into each other (Collisionality/Knudsen number).
2. Validating Observations: Spacecraft (like the WIND mission) have already seen these lopsided electron crowds. This paper explains why they exist using the laws of physics, rather than just saying "it happens."
3. A New Tool: It suggests that if we measure how lopsided the electrons are, we can actually calculate the "Knudsen number" and understand the underlying physics of the solar wind without needing to measure every single collision.

Summary

Think of the solar wind electrons as a crowd of runners. Usually, they run in a straight line. But because the "ground" they run on (temperature and density) changes, and because they rarely bump into each other, the whole crowd starts to lean.

This paper proves that the angle of that lean is directly determined by how "bumpy" the path is and how "empty" the space is. They turned a complex, lopsided shape into a simple number (the Knudsen number) that tells us exactly how the solar wind is behaving.

Technical Summary

1. Problem Statement

Space plasmas, particularly the solar wind electron population, frequently exhibit Velocity Distribution Functions (VDFs) that deviate significantly from thermal equilibrium (Maxwellian). These distributions often feature:

• Suprathermal tails: High-energy particles better described by Kappa distributions.
• Asymmetry (Skewness): A field-aligned asymmetry where the distribution is skewed relative to the mean velocity, often associated with the "strahl" electron population.

While the Skew-Kappa distribution (specifically the Core-Strahlo or CS model) has been successfully used to empirically fit observational data (e.g., from the WIND spacecraft), the physical origin of the skewness parameter ($\delta$) remains theoretically underexplored. Specifically, there is a lack of a rigorous derivation linking the macroscopic transport properties of the plasma (gradients, collisions) to the emergence and maintenance of this asymmetry.

2. Methodology

The authors employ a kinetic theory approach based on the Boltzmann Transport Equation (BTE) to derive a theoretical relationship between the skewness parameter and plasma macro-dynamics.

• Governing Equation: They start with the BTE for electrons, incorporating a Krook-like collisional operator on the right-hand side to represent collisional relaxation.
  $$\frac{\partial f_e}{\partial t} + \vec{v} \cdot \nabla f_e + \frac{q_e}{m_e}\left(\vec{E} + \frac{\vec{v}}{c} \times \vec{B}\right) \cdot \nabla_v f_e = \left(\frac{\partial f_e}{\partial t}\right)_{coll}$$
• Ansatz: They assume a steady-state solution in the form of a Skew-Kappa distribution ($f_e^{\kappa\delta}$), which includes a skewness parameter $\delta_e$ and a kappa index $\kappa$.
• Perturbation Analysis:
  • They assume small asymmetry ($\delta_e \ll 1$) and perform a second-order Taylor expansion of the distribution function with respect to $\delta_e$.
  • They introduce a velocity-dependent collision frequency $\nu_e(v)$. Unlike standard constant-frequency models, they propose a model that converges to a finite value at low velocities and follows a $v^{-3}$ power law at high velocities (consistent with Helander & Sigmar theory), specifically tailored to the Skew-Kappa form.
• Derivation:
  • The BTE is reduced to a balance between streaming (left-hand side) and collisional relaxation (right-hand side).
  • The equation is integrated over velocity space to derive a closed-form expression for $\delta_e$ in terms of macroscopic gradients (density $n$ and parallel thermal velocity $W_{\parallel, e}$).

3. Key Contributions

• Theoretical Derivation of Skewness: The paper provides the first analytical derivation linking the skewness parameter $\delta_e$ directly to macroscopic plasma gradients and collisional effects, moving beyond empirical fitting.
• Novel Collision Frequency Model: The authors propose a specific form for the collision frequency (Eq. 10/11) that is consistent with the Skew-Kappa distribution. This model resolves singularities at low velocities while preserving the correct high-velocity asymptotic behavior required for Coulomb collisions.
• Identification of the Knudsen Number Connection: The study establishes a direct proportionality between the skewness parameter and the effective Knudsen number ($K_N$), unifying kinetic asymmetry with transport theory.
• Decoupling of $\delta$ and $\kappa$: The analysis explains why observational data often shows that the skewness parameter $\delta$ is largely independent of the kappa index $\kappa$. The derived formula shows $\delta$ depends on gradients, while $\kappa$ only appears in a prefactor that saturates rapidly for large $\kappa$.

4. Key Results

The central result is the derived expression for the skewness parameter $\delta_e$:

$$\delta_e = \Psi_\alpha(\kappa) \left( \frac{1}{2} K_T + K_n \right)$$

Where:

• $\Psi_\alpha(\kappa)$ is a function of the kappa index that saturates quickly (approaching $\approx 2.5$) for $\kappa > 4$, explaining the weak dependence of skewness on the spectral index.
• $K_T = \lambda_{fp} / L_T$ is the Temperature Knudsen number (ratio of mean free path to temperature gradient scale).
• $K_n = \lambda_{fp} / L_n$ is the Density Knudsen number (ratio of mean free path to density gradient scale).
• $\lambda_{fp} = W_{\parallel, e} / \nu_{e,0}$ is the electron mean free path.

Implications of the Results:

1. Linear Relationship: Skewness is linearly proportional to the Knudsen numbers. This confirms that asymmetry arises from non-equilibrium transport driven by spatial gradients.
2. Consistency with Heat Flux: Since electron heat flux ($q_e$) is known to be proportional to the Knudsen number (Spitzer & Harm theory), and previous studies linked heat flux to skewness, this work provides the theoretical bridge: Skewness $\sim$ Heat Flux $\sim$ Knudsen Number.
3. Reynolds Number Connection: The authors discuss a potential link to turbulence, suggesting that the small-skewness limit is consistent with high Reynolds number flows where turbulent structures constrain particle free paths.

5. Significance

• Unification of Kinetic and Fluid Descriptions: The paper bridges the gap between kinetic descriptions (non-Maxwellian VDFs) and fluid transport parameters (Knudsen numbers, heat flux). It demonstrates that skewness is not just a fitting parameter but a measurable proxy for non-equilibrium transport.
• Validation of the CS Model: The results provide a physical justification for the Core-Strahlo (CS) model, showing that a single skewness parameter can robustly characterize asymmetry because it is driven by macroscopic gradients rather than the specific shape of the high-energy tail ($\kappa$).
• Observational Predictions: The study offers a testable prediction: the skewness parameter $\delta_e$ should scale with the effective Knudsen number derived from independent measurements of temperature and density gradients in the solar wind. This can be tested using data from missions like Parker Solar Probe, Helios, and WIND.
• Understanding Solar Wind Dynamics: By linking asymmetry to collisional age and gradients, the work enhances the understanding of how solar wind electrons evolve from the Sun to 1 AU, specifically regarding the scattering of strahl electrons into the halo population.

In conclusion, this paper establishes that skewness in space plasma distributions is a direct consequence of collisional transport in the presence of macroscopic gradients, quantified effectively by the Knudsen number.