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

Technische Samenvatting: Bivariaat Onvolledige Bessel-kernen voor de Eerste Nietlineaire Vlasov–Maxwell Respons

Probleemstelling
De kinetische theorie van nietlineaire golfkoppeling in gemagnetiseerde plasma's steunt traditioneel op een zwak nietlineaire expansie van het Vlasov–Maxwell-systeem rond een homogeen evenwicht. De standaardbenadering drukt de plasmarespons uit via een dubbele cyclotron-harmonische expansie. Hoewel expliciet, wordt deze representatie algebraïsch omslachtig, met name wanneer elektromagnetische golven onder willekeurige hoeken ten opzichte van het achtergrondmagnetisch veld voortplanten. In de eerste nietlineaire correctie omvat de respons twee gekoppelde harmonische indices: één afkomstig van de externe karakteristieke integratie van de output-modus en een andere van de interne lineaire respons die als bron fungeert. Deze structuur leidt tot complexe dubbele sommen van Bessel-functies en geneste resonantiedenominatoren, zoals te zien in klassieke formuleringen door Liu en Tripathi (1986) en daaropvolgende werken.

Methodologie
Dit artikel stelt een alternatieve formulering voor die de directe expansie van de Larmor-fase in een harmonische reeks vermijdt. In plaats daarvan behouden de auteurs de fase langs de ongebaande baan en evalueren zij de karakteristieke integralen direct.

1. Lineaire Basis: De auteurs herinneren aan een recente lineaire formulering (Ricci 2026) waarbij de karakteristieke integraal wordt geëvalueerd in termen van de onvolledige Bessel-functie $G_\mu(z, \psi)$. Deze functie vervangt de oneindige som van cyclotron-harmonischen voordat de susceptibiliteitstensor wordt geprojecteerd.
2. Nietlineaire Uitbreiding: De eerste nietlineaire Vlasov-vergelijking wordt opgelost met dezelfde karakteristieke methode. De nietlineaire bronterm bestaat uit de Lorentz-kracht van één modus werkend op de lineaire respons van een andere. Door de interne lineaire respons in zijn onvolledige Bessel-vorm uit te drukken, werkt de externe karakteristieke integraal op een verschoven $G_{\mu_q}$-respons.
3. Kernafleiding: Dit proces genereert van nature een bivariaat onvolledige Bessel-functie, aangeduid als $G^{(r)}_{\mu,\nu}(z, \psi; w, \chi)$. Deze functie dient als de baan-resolvent kern voor het eerste nietlineaire probleem, waarbij de integer $r$ rekening houdt met eindige harmonische verschuivingen gegenereerd door transversale snelheidsfactoren en gyrofase-afgeleiden.
4. Herstel van Klassieke Resultaten: De auteurs leiden de nietlineaire distributiefunctie expliciet af in termen van deze bivariate functies. Ze demonstreren vervolgens dat het expanderen van deze functies met behulp van de Jacobi–Anger-formule en Turkin's functie de traditionele dubbel-harmonische Liu–Tripathi-formule terugbrengt, waarmee de nieuwe formulering wordt gevalideerd.
5. Susceptibiliteitsprojectie: Ten slotte leiden de auteurs de nietlineaire susceptibiliteitstensor af door de bivariate distributie te projecteren op snelheidsmomenten. Dit wordt bereikt door "geprojecteerde bivariate kernen", wat hoekige Fourier-projecties van de bivariate functies zijn, waardoor het noodzakelijk is niet terug te keren naar de volledige dubbel-harmonische expansie.

Belangrijkste Bijdragen en Resultaten

• Definitie van Bivariate Kernen: Het artikel introduceert $G^{(r)}_{\mu,\nu}$ als de natuurlijke baan-integraal die voortvloeit uit het eerste nietlineaire karakteristieke probleem. Het stelt vast dat deze functie niet louter een notatieve gemak is, maar het fundamentele object dat voortkomt uit de iteratie van de karakteristieke methode.
• Structurele Herorganisatie: De nietlineaire distributiefunctie wordt compact uitgedrukt met behulp van differentiaaloperatoren die op deze bivariate kernen inwerken. De formulering isoleert de hoekige contracties in een kleine set herbruikbare speciale-functieblokken.
• Wiskundige Identiteiten: Een uitgebreide set identiteiten voor de bivariate functies is afgeleid, inclusioneel superpositie-eigenschappen, harmonische expansies, recursie-relaties en differentiaalvergelijkingen (bijv. de externe Bessel-vergelijking en de interne Turkin-recursie).
• Consistentie-verificatie: Het artikel biedt een rigoureuze afleiding die aantoont dat de bivariate formulering wiskundig equivalent is aan de klassieke dubbel-harmonische benadering. De traditionele formule wordt teruggewonnen door de bivariate functie te expanderen en de Graf-additiestelling toe te passen om de tussenliggende harmonische indices te contrageren.
• Expliciete Susceptibiliteitsformule: De nietlineaire susceptibiliteitstensor wordt gepresenteerd als een eindige som over negen elementaire bivariate blokken (geïndexeerd door $a, b \in \{-1, 0, 1\}$), gevolgd door een tweedimensionale snelheidsintegraal. De hoekafhankelijkheid wordt gevangen in geprojecteerde kernen $K^{(r,a)}_N$ en hun afgeleiden.

Betekenis en Claims
Het artikel claimt dat het primaire voordeel van deze formulering zowel structureel als potentieel computationeel is:

• Structureel: Het identificeert het specifieke speciale-functie object (de bivariate onvolledige Bessel-kern) dat ten grondslag ligt aan de eerste nietlineaire cyclotron-respons, waarbij de theorie wordt georganiseerd rond de geometrie van de karakteristieke stroom in plaats van de algebraïsche complexiteit van harmonische sommen.
• Computationeel: De auteurs suggereren dat voor numerieke implementatie de bivariate formulering een meer gestructureerd doel biedt. In plaats van lange, oscillerende dubbele sommen van Bessel-functies met verschoven indices te hanteren, kan men de geprojecteerde kernen en hun afgeleiden tabelleren en hergebruiken. Deze flexibiliteit kan voordelig zijn in regimes waar directe harmonische sommatie traag convergeert of numeriek delicaat is vanwege annuleringen.
• Schaalbaarheid: De auteurs merken op dat deze benadering recursief uitbreidbaar is naar hogere orden. De $n$-de orde distributie wordt naar verwachting uitdrukbaar zijn in termen van $n$-variante baan-kernen, gegenereerd door voor elke nietlineaire interactie één karakteristieke integraal toe te voegen.

Het artikel behoudt een bescheiden standpunt door te erkennen dat het numerieke voordeel geen universele reductie in computationele kosten is, maar afhangt van specifieke plasmaparameters (Larmor-radii, detuning, demping). Er wordt expliciet vermeld dat de constructie van toepassing is op de vaste-karakteristiek perturbatiehiërarchie en geen rekening houdt met nietlineaire baan-correcties (zoals trapping of ponderomote driften) die de karakteristieke operator zelf zouden modificeren.