A new class of special functions arising in plasma linear susceptibility tensor calculations

This paper introduces and analyzes a new class of special functions derived from Bessel, Anger, and Weber functions to solve inhomogeneous Bessel differential equations, demonstrating how their recurrence relations enable a more efficient derivation of the linear susceptibility tensor in hot, magnetized plasmas by avoiding the slow convergence of traditional infinite Bessel series.

The Gist

Imagine you are trying to predict how a crowd of charged particles (like a hot, magnetized plasma) reacts when you shake the table they are sitting on. In physics, this "reaction" is called susceptibility. To calculate it, scientists usually have to solve a massive, messy puzzle involving waves and circles.

For decades, the standard way to solve this puzzle was like trying to count every single grain of sand on a beach to measure the beach's weight. The math involved adding up infinite lists of special numbers (called Bessel functions).

Here is the problem: When the particles are moving fast or the waves are long (a regime called "large gyro-radius"), this method of counting grains of sand becomes incredibly slow. It's like trying to fill a swimming pool with a teaspoon; you have to add millions of spoonfuls just to get a little water, and the computer takes forever to finish the calculation.

The "Magic Shortcut"

This paper introduces a new class of mathematical tools (special functions) that act like a magic shortcut or a pre-packed suitcase for this calculation.

Here is the breakdown of what the author, R. Ricci, did, using simple analogies:

1. The Old Way: The Infinite Staircase

Traditionally, to find the answer, physicists used a famous formula (Jacobi-Anger) that broke the problem down into an infinite staircase.

• The Analogy: Imagine you need to get to the top of a mountain. The old method says, "Take one step, then another, then another... forever."
• The Problem: If the mountain is huge (large particle radius), you have to take billions of steps. The computer gets tired and the calculation becomes impractical.

2. The New Way: The Helicopter

Ricci discovered that the function used in the old method (defined by Qin, Philips, and Davidson years ago) isn't just a random sum. It is actually a specific solution to a differential equation (a rule that describes how things change).

• The Analogy: Instead of walking up the infinite staircase step-by-step, Ricci realized you could just fly a helicopter straight to the top.
• He showed that this "helicopter function" (let's call it $G_\mu$) follows a very specific, simple set of rules (recurrence relations). You don't need to count the steps; you just need to know the starting point and the rule of the road.

3. The "Incomplete" Connection

The paper connects this new function to a family of older, well-known mathematical characters called Anger and Weber functions.

• The Analogy: Think of the old Bessel functions as "Complete" apples (whole fruits). The new functions are like "Incomplete" apples (slices). Ricci showed that by understanding how these slices fit together, you can reconstruct the whole fruit without ever having to slice it yourself.
• He proved that these new functions satisfy a specific "recipe" (Nielsen's requirement) that guarantees they are the correct solution to the plasma problem.

4. Why This Matters for Plasma Physics

In the real world, plasma is used in fusion energy (like the sun on Earth) and space weather prediction.

• The Old Problem: When calculating how plasma behaves in a magnetic field, the "infinite staircase" method would often crash computers or take days to run because the series converged too slowly.
• The New Solution: By using Ricci's new functions, scientists can calculate the plasma's reaction directly.
  • Instead of summing an infinite list of products (which is slow), the new method uses a finite, clean formula involving just a few terms.
  • It's like switching from a calculator that adds 1+1+1+1... to a calculator that just knows "4 x 100 = 400."

The Big Picture

The paper does two main things:

1. Mathematical Detective Work: It identifies exactly what these mysterious functions are. They aren't just random tools; they are solutions to a specific type of wave equation with a "kick" (an inhomogeneous term) on the side.
2. Practical Engineering: It provides a new, much faster way to calculate how hot plasma reacts to electromagnetic waves. This is crucial for designing better fusion reactors and understanding space storms.

In a nutshell: The author took a messy, slow, infinite math problem that plagued plasma physicists for years and found a "cheat code." Instead of grinding through billions of numbers, we can now use a sleek, elegant formula that gets the job done instantly, even in the most difficult conditions.

Technical Summary

1. Problem Statement

The calculation of the linear susceptibility tensor for a hot, magnetized plasma is a fundamental task in plasma kinetic theory. Traditionally, this involves solving the linearized Vlasov equation using the method of characteristics.

• The Bottleneck: Standard approaches rely on the Jacobi-Anger expansion to handle the phase factors associated with particle gyration. This results in final expressions containing infinite sums of products of Bessel functions of integer order.
• Numerical Difficulty: In regimes where the particle gyro-radius is large compared to the wavelength (large $z = k_\perp v_\perp / \omega_c$), these series converge extremely slowly. This necessitates summing a huge number of terms, making numerical evaluation computationally expensive and prone to errors.
• Previous Attempts: Qin, Phillips, and Davidson (2007) proposed a method to avoid these infinite sums by introducing a specific integral function. However, the mathematical nature of this function and its general properties were not fully explored, and the derivation remained somewhat ad-hoc.

2. Methodology

The author introduces and rigorously analyzes a new class of special functions, denoted as $G_\mu(z, \psi)$, defined by an improper integral (or equivalently a definite integral over one period):
$$G_\mu(z, \psi) = c_\mu \int_0^{2\pi} \frac{d\lambda}{2\pi} e^{-iz\sin(\lambda+\psi) - i\mu\lambda}$$
where $c_\mu = \pi e^{i\pi\mu} / \sin(\pi\mu)$, $z$ and $\psi$ are real variables, and $\mu$ is a complex parameter.

The paper employs the following analytical techniques:

1. Differential Equation Analysis: The author demonstrates that $G_\mu(z, \psi)$ is not just an integral but a specific solution to an inhomogeneous Bessel Ordinary Differential Equation (ODE).
2. Recurrence Relations: By exploiting the properties of the integrand, the author derives recurrence relations connecting $G_\mu$ with $G_{\mu \pm 1}$.
3. Connection to Classical Functions: The paper establishes explicit links between $G_\mu$ and known special functions, specifically Anger functions ($J_\mu$), Weber functions ($E_\mu$), and Incomplete Anger-Weber functions.
4. Nielsen's Functional Equations: The author frames the problem within the context of Nielsen's functional equations, showing that $G_\mu$ satisfies specific conditions required for solutions to inhomogeneous Bessel equations.
5. Application to Plasma Physics: The derived properties are applied to the linearized Vlasov equation to reformulate the calculation of the susceptibility tensor $\chi(k, \Omega)$.

3. Key Contributions

A. Mathematical Characterization of $G_\mu(z, \psi)$

• Inhomogeneous Bessel ODE: The paper proves that $G_\mu(z, \psi)$ satisfies:
  $$\mathcal{D}_{B, z, \mu} G_\mu(z, \psi) = \ell_\mu(z, \psi)$$
  where $\mathcal{D}_{B, z, \mu}$ is the Bessel differential operator and the source term is $\ell_\mu(z, \psi) = e^{-iz\sin\psi}(z\cos\psi - \mu)$.
• Series Expansion: An alternative representation is derived using the Jacobi-Anger formula:
  $$G_\mu(z, \psi) = \sum_n \frac{J_n(z)}{n+\mu} e^{-in\psi}$$
  This confirms the function's relation to the Fourier coefficients of the periodic solution.
• Relation to Incomplete Anger-Weber Functions: The author derives a compact closed-form expression relating $G_\mu$ to the incomplete Anger-Weber function $A_\mu(z, \psi)$:
  $$G_\mu(z, \psi) = e^{i\mu\psi} \left[ \frac{\pi}{\sin(\pi\mu)} J_\mu(z) - i\pi A_\mu(-z, \psi) \right]$$
• Recurrence Relations: New recurrence relations are established that allow shifting the index $\mu$ without generating infinite series, facilitating recursive numerical evaluation.

B. Reformulation of Plasma Susceptibility

• Avoidance of Infinite Sums: By utilizing the recurrence properties of $G_\mu$, the author derives the linear susceptibility tensor without resorting to the Jacobi-Anger expansion.
• Generalization: The derived expressions generalize previous results (e.g., Qin et al.) by handling arbitrary wave propagation angles ($\phi_k \neq 0$), whereas previous methods were often restricted to Stix coordinates ($\phi_k = 0$).
• Newberger's Sum Rule: The paper demonstrates that the integrals arising in the new formulation are equivalent to Newberger's sum rule (which sums products of integer-order Bessel functions into non-integer order products). However, the new method derives these results directly through the properties of $G_\mu$, bypassing the cumbersome algebraic summation of infinite series.

4. Key Results

The paper presents the explicit elements of the susceptibility tensor matrix $\chi(k, \Omega)$ in Cartesian coordinates (Eq. 87a–87i).

• Structure: The tensor elements are expressed as combinations of products of Bessel functions of non-integer order ($J_\mu$ and $J_{-\mu}$) and their derivatives, rather than infinite sums of integer-order Bessel functions.
• Efficiency: The resulting formulas involve finite combinations of special functions. This eliminates the convergence issues associated with large gyro-radii, as the evaluation of $J_\mu(z)$ for non-integer $\mu$ is numerically stable and efficient even for large arguments.
• Consistency: The results are shown to be mathematically equivalent to the traditional approach but derived through a more direct and algebraically simpler path.

5. Significance

• Numerical Robustness: The primary significance is the removal of the "slow convergence" bottleneck in plasma simulations. This allows for accurate and efficient calculation of linear wave properties in hot plasmas, particularly in regimes relevant to fusion research (e.g., ion cyclotron resonance heating) where large gyro-radii are common.
• Theoretical Unification: The work bridges a gap between plasma kinetic theory and the classical theory of special functions. It clarifies that the "trick" used by Qin et al. is actually a specific solution to a well-defined inhomogeneous Bessel problem, connecting it to the broader literature on Anger-Weber functions and Nielsen's conditions.
• Future Applications: The author notes that this systematic approach simplifies the algebraic complexity of the linear case and lays the groundwork for extending these methods to nonlinear plasma physics, where such calculations are currently prohibitively difficult.

In summary, the paper provides a rigorous mathematical foundation for a class of special functions that serve as a powerful computational tool, transforming the calculation of plasma susceptibility from a problem of summing slowly converging series into a problem of evaluating finite combinations of well-behaved special functions.