Bivariate incomplete-Bessel kernels for the first nonlinear Vlasov-Maxwell response

技术摘要：用于一阶非线性 Vlasov–Maxwell 响应的双变量不完全贝塞尔核

问题陈述
磁化等离子体中非线性波耦合的动力学理论传统上依赖于围绕均匀平衡态的 Vlasov–Maxwell 系统的弱非线性展开。标准方法将等离子体响应表示为双回旋谐波展开。虽然这种表示法是显式的，但在电磁波以任意角度相对于背景磁场传播时，其代数处理变得异常繁琐。在一阶非线性修正中，响应涉及两个耦合的谐波指数：一个源自输出模的外部特征积分，另一个源自作为源项的内部线性响应。这种结构导致了复杂的贝塞尔函数双重求和以及嵌套的共振分母，正如 Liu 和 Tripathi (1986) 及其后续工作中所见。

研究方法
本文提出了一种替代方案，避免了立即将拉莫尔相位（Larmor phase）展开为谐波级数。相反，作者保留了沿未扰动轨道的相位，并直接评估特征积分。

1. 线性基准： 作者回顾了近期的一种线性表述（Ricci 2026），其中特征积分通过不完全贝塞尔函数 $G_\mu(z, \psi)$ 进行评估。在悬数（susceptibility）张量投影之前，该函数取代了无穷次的回旋谐波求和。
2. 非线性扩展： 使用相同的特征法求解第一阶非线性 Vlasov 方程。非线性源项由一个模的洛伦兹力作用于另一个模的线性响应组成。通过将内部线性响应表示为这种不完全贝塞尔形式，外部特征积分作用于一个偏移的 $G_{\mu_q}$ 响应。
3. 核函数推导： 这一过程自然生成了一个双变量不完全贝塞尔函数，记作 $G^{(r)}_{\mu,\nu}(z, \psi; w, \chi)$。该函数作为一阶非线性问题的轨道解算子核（orbit-resolvent kernel），其中整数 $r$ 解释了由横向速度因子和回旋相位导数产生的有限谐波偏移。
4. 恢复经典结果： 作者利用这些双变量函数显式地推导了非线性分布函数。随后，他们证明了通过使用 Jacobi–Anger 公式和 Turkin 函数展开这些函数，可以恢复传统的双谐波 Liu–Tripathi 公式，从而验证了新表述的有效性。
5. 悬数投影： 最后，本文通过将双变量分布投影到速度矩来推导非线性悬数张量。这是通过“投影双变量核”（即双变量函数的角傅里叶投影）实现的，避免了返回完整的双谐波展开。

主要贡献与结果

• 定义双变量核： 本文引入了 $G^{(r)}_{\mu,\nu}$ 作为一阶非线性特征问题所产生的自然轨道积分。作者确立了该函数不仅是一个符号上的便利，更是通过迭代特征法产生的基本对象。
• 结构重组： 非线性分布函数使用作用于这些双变量核的微分算子进行紧凑表示。该表述将角度收缩（angular contractions）隔离到一组可重复使用的特殊函数块中。
• 数学恒等式： 文中推导了一套关于双变量函数的全面恒等式，包括叠加性质、谐波展开、递推关系以及微分方程（例如外部贝塞尔方程和内部 Turkin 递推）。
• 一致性验证： 本文提供了一个严谨的推导，表明双变量表述在数学上等价于经典的双谐波方法。通过展开双变量函数并应用 Graf 加法定理来收缩中间谐波指数，即可恢复传统公式。
• 显式悬数公式： 非线性悬数张量被表示为九个基本双变量块（由 $a, b \in \{-1, 0, 1\}$ 索引）的有限和，随后进行二维速度积分。角度依赖性被封装在投影核 $K^{(r,a)}_N$ 及其导数中。

意义与主张
本文声称，该表述的主要优势在于结构性和潜在的计算性：

• 结构性： 它识别出了作为一阶非线性回旋响应基础的特定特殊函数对象（双变量不完全贝塞尔核），将理论组织在特征流的几何结构之上，而非代数复杂的谐波求和之上。
• 计算性： 作者指出，对于数值实现，双变量表述提供了一个更具结构性的目标。与其处理具有偏移指数且具有长程振荡特性的贝塞尔函数双重求和，不如对投影核及其导数进行查表和复用。在直接谐波求和收敛缓慢或因抵消效应而导致数值敏感的机制下，这种灵活性可能具有优势。
• 可扩展性： 作者指出，这种方法可以递归地扩展到更高阶。预计 $n$ 阶分布可以用 $n$-元轨道核来表示，通过为每一次非线性相互作用增加一个特征时间积分来生成。

本文保持了审慎的态度，承认其数值优势并非普遍的计算成本降低，而是取决于特定的等离子体参数（拉莫尔半径、失谐量、阻尼）。文中明确指出，该构建适用于固定的特征摄动层级，并不涵盖会修改特征算子本身的非线性轨道修正（如捕获或庞德罗莫夫漂移）。