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

Resumo Técnico: Núcleos de Bessel Incompletos Bivariados para a Primeira Resposta Vlasov–Maxwell Não Linear

Enunciado do Problema
A teoria cinética do acoplamento de ondas não lineares em plasmas magnetizados baseia-se tradicionalmente em uma expansão fracamente não linear do sistema Vlasov–Maxwell em torno de um equilíbrio homogêneo. A abordagem padrão expressa a resposta do plasma por meio de uma expansão de harmônicos ciclotrônicos dupla. Embora explícita, essa representação torna-se algebricamente complexa, particularmente quando ondas eletromagnéticas se propagam em ângulos arbitrários em relação ao campo magnético de fundo. Na primeira correção não linear, a resposta envolve dois índices harmônicos acoplados: um originado da integração da característica externa do modo de saída e outro da resposta linear interna que atua como fonte. Essa estrutura leva a somas duplas complexas de funções de Bessel e denominadores de ressonância aninhados, como visto nas formulações clássicas de Liu e Tripathi (1986) e trabalhos subsequentes.

Metodologia
Este artigo propõe uma formulação alternativa que evita a expansão imediata da fase de Larmor em uma série harmônica; em vez disso, os autores mantêm a fase ao longo da órbita não perturbada e avaliam as integrais características diretamente.

1. Linha de Base Linear: Os autores recordam uma formulação linear recente (Ricci 2026) onde a integral característica é avaliada em termos da função de Bessel incompleta $G_\mu(z, \psi)$. Esta função substitui a soma infinita de harmônicos ciclotrônicos antes que a susceptibilidade seja projetada.
2. Extensão Não Linear: A primeira equação de Vlasov não linear é resolvida utilizando o mesmo método característico. O termo de fonte não linear consiste na força de Lorentz de um modo atuando sobre a resposta linear de outro. Ao expressar a resposta linear interna em sua forma de Bessel incompleto, a integral característica externa atua sobre uma resposta $G_{\mu_q}$ deslocada.
3. Derivação do Núcleo: Este processo gera naturalmente uma função de Bessel incompleta bivariada, denotada como $G^{(r)}_{\mu,\nu}(z, \psi; w, \chi)$. Esta função serve como o núcleo do resolvente de órbita para o primeiro problema não linear, onde o inteiro $r$ contabiliza os deslocamentos harmônicos finitos gerados por fatores de velocidade transversais e derivadas de autofase.
4. Recuperação de Resultados Clássicos: Os autores derivam a função de distribuição não linear explicitamente em termos dessas funções bivariadas. Eles então demonstram que a expansão dessas funções usando a fórmula de Jacobi–Anger e a função de Turkin recupera a tradicional fórmula de dupla harmônica de Liu–Tripathi, validando assim a nova formulação.
5. Projeção da Susceptibilidade: Finalmente, o artigo deriva o tensor de susceptibilidade não linear projetando a distribuição bivariada sobre momentos de velocidade. Isso é alcançado através de "núcleos bivariados projetados", que são projeções de Fourier angulares das funções bivariadas, evitando a necessidade de retornar à expansão completa de dupla harmônica.

Principais Contribuições e Resultados

• Definição de Núcleos Bivariados: O artigo introduz $G^{(r)}_{\mu,\nu}$ como o núcleo de órbita natural gerado pelo primeiro problema característico não linear. Estabelece que esta função não é meramente uma conveniência notacional, mas o objeto fundamental resultante da iteração do método característico.
• Reorganização Estrutural: A função de distribuição não linear é expressa de forma compacta usando operadores diferenciais atuando sobre esses núcleos bivariados. A formulação isola as contrações angulares em um pequeno conjunto de blocos de funções especiais reutilizáveis.
• Identidades Matemáticas: Um conjunto abrangente de identidades para as funções bivariadas é derivado, incluindo propriedades de superposição, expansões harmônicas, relações de recorrência e equações diferenciais (ex: a equação de Bessel externa e a recorrência de Turkin interna).
• Verificação de Consistência: O artigo fornece uma derivação rigorosa mostrando que a formulação bivariada é matematicamente equivalente à abordagem clássica de dupla harmônica. A fórmula tradicional é recuperada expandindo a função bivariada e aplicando o teorema de adição de Graf para contrair os índices harmônicos intermediários.
• Fórmula Explícita de Susceptibilidade: O tensor de susceptibilidade não linear é apresentado como uma soma finita sobre nove blocos bivariados elementares (indexados por $a, b \in \{-1, 0, 1\}$), seguida por uma integral de velocidade bidimensional. A dependência angular é encapsulada nos núcleos projetados $K^{(r,a)}_N$ e suas derivadas.

Significância e Alegações
O artigo alega que a principal vantagem desta formulação é tanto estrutural quanto potencialmente computacional:

• Estrutural: Identifica o objeto de função especial específico (o núcleo de Bessel incompleto bivariado) que subjaz à primeira resposta ciclotrônica não linear, organizando a teoria em torno da geometria do fluxo característico em vez da complexidade algébrica das somas harmônicas.
• Computacional: Os autores sugerem que, para implementação numérica, a formulação bivariada oferece um alvo mais estruturado. Em vez de lidar com somas duplas de funções de Bessel longas e oscilatórias com índices deslocados, pode-se tabular e reutilizar os núcleos projetados e suas derivadas. Essa flexibilidade pode ser vantajosa em regimes onde a soma harmônica direta possui convergência lenta ou é numericamente delicada devido a cancelamentos.
• Escalabilidade: Os autores observam que esta abordagem se estende recursivamente para ordens superiores. Espera-se que a $n$-ésima distribuição possa ser expressa em termos de núcleos de órbita $n$-variantes, gerados pela adição de uma integral de tempo característica para cada interação não linear.

O artigo mantém uma postura modesta, reconhecendo que a vantagem numérica não é uma redução universal do custo computacional, mas depende de parâmetros específicos do plasma (raios de Larmor, desajustes, amortecimento). Afirma explicitamente que a construção se aplica à hierarquia perturbativa de característica fixa e não considera correções de órbita não lineares (como aprisionamento ou drifts ponderomotores) que modificariam o próprio operador característico.