A Variational Framework for Guiding-Center Kinetics, Anisotropic Equilibria, and Quasisymmetry in Stellarators

Resumen Técnico: Un Marco Variacional para la Cinética de Centro de Guía, Equilibrios Anisotrópicos y Cuasisimetría en Estelaradores

Planteamiento del Problema
El artículo aborda una brecha conceptual fundamental en la teoría de los estelaradores respecto a la relación entre la cuasisimetría (QS), la anisotropía de presión y los invariantes cinéticos. En el confinamiento magnético toroidal, los estelaradores carecen de la simetría geométrica continua de los tokamaks, lo que conduce a un movimiento de centro de guía no integrable y a un transporte neoclásico incrementado. Si bien se ha establecido que la QS es una condición geométrica (donde la intensidad del campo magnético $B$ depende de un único ángulo helicoidal) que restaura la integrabilidad, su estatus físico en relación con la estructura cinética de los plasmas magnetizados sigue siendo incompletamente comprendido.

Los tratamientos existentes suelen introducir la QS como una restricción geométrica externa independiente de la ecuación de estado del plasma, mientras que la teoría de equilibrio típicamente asume una presión escalar. Sin embargo, los plasmas magnetizados sin colisiones son intrínsecamente anisotrópicos debido al rápido movimiento de giro, lo que resulta en un tensor de esfuerzo giotrópico con presiones paralelas ($p_\parallel$) y perpendiculares ($p_\perp$) distintas. El artículo identifica una inconsistencia crítica: se sabe que la suposición de presión escalar está sobredeterminada en geometrías 3D (el Conundrum de Garren–Boozer), pero la necesidad de la presión anisotrópica para los equilibrios 3D no ha sido derivada rigurosamente desde los primeros principios dentro del contexto de la QS. Los autores se preguntan si la QS, el equilibrio anisotrópico y los invariantes cinéticos son características independientes o facetas de una única estructura variacional.

Metodología
Los autores desarrollan un marco variacional euleriano unificado basado en la acción de Vlasov–Maxwell de centro de guía. La metodología procede a través de los siguientes pasos:

1. Formulación de la Acción: Partiendo del Lagrangiano de partícula única de centro de guía, los autores realizan un promedio de fase de giro para reducir el espacio de fase de seis a cinco dimensiones. Construyen un funcional de acción total que comprende la contribcción de la partícula de centro de guía y la contribución del campo electromagnético.
2. Variaciones Dinámicamente Accesibles: Las variaciones admisibles de la función de distribución se restringen a formas "dinámicamente accesibles" (transformaciones canónicas infinitesimales) que preservan la topología del espacio de fase y los invariantes de Casimir, de forma análoga a la restricción de Newcomb en la MHD ideal.
3. Derivación de Ecuaciones: Al aplicar variaciones restringidas a la acción, los autores derivan la ecuación de Vlasov de centro de guía, las ecuaciones de Maxwell y el balance de fuerzas macroscópico.
4. Aplicación del Teorema de Noether: Los autores aplican el teorema de Noether a la acción para derivar leyes de conservación. Específicamente, analizan la invariancia de traslación espacial para obtener la ley de conservación del momento y el balance de fuerzas macroscópicas.
5. Identificación de Simetría: La cuasisimetría se identifica no como una restricción externa, sino como una simetría continua (simetría de Noether) de la acción de centro de guía. Los autores investigan las condiciones bajo las cuales la acción permanece invariante bajo una transformación de coordenadas generada por un campo vectorial $\xi$.
6. Análisis de Estabilidad: La estabilidad cinética de los equilibrios resultantes se analiza utilizando el método de energía–Casimir, examinando la segunda variación del funcional de energía libre.

Contribuciones Clave y Resultados

• Estructura Variacional Unificada: El artículo demuestra que la teoría cinética de centro de guía, el balance de fuerzas macroscópicas con presión anisotrópica y las restricciones de QS no son independientes, sino que emergen como diferentes facetas de una única estructura variacional derivada de la acción de Vlasov–Maxwell de centro de guía.
• Derivación del Balance de Fuerzas Anisotrópico: Sin invocar cierres fenomenológicos, la variación de la acción produce la ecuación de balance de fuerzas macroscópicas $\mathbf{J} \times \mathbf{B}/c = \nabla \cdot \mathbf{\Pi}$, donde $\mathbf{\Pi}$ es el tensor de esfuerzo giotrópico. Esto recupera directamente la ecuación de balance de fuerzas de Chew–Goldberger–Low (CGL) a partir de los momentos cinéticos, estableciendo que la anisotropía de presión es una necesidad estructural derivada de la conservación del momento magnético $\mu$.
• La QS como Simetría de Noether: Los autores identifican rigurosamente la QS como una simetría continua de la acción de centro de guía. Cuando la intensidad del campo magnético satisface la condición de QS, la acción es invariante bajo una transformación específica, produciendo un momento generalizado conservado $P_\xi = (q/c)\mathbf{A}^* \cdot \boldsymbol{\xi}$. Esta cantidad conservada restaura la integrabilidad de las trayectorias de las partículas en las superficies magnéticas.
• Restricciones de Solvabilidad Geométrica: El requisito de que el segundo invariante adiabático (la órbita de deriva promediada por rebote) sea una función de flujo conduce a una restricción geométrica específica. Los autores derivan una condición de integrabilidad independiente de coordenadas:
  $$f_T \equiv \nabla\psi \cdot (\nabla B \times \nabla(B \cdot \nabla B)) = 0$$
  Esta condición, introducida previamente en formulaciones variacionales restringidas (por ejemplo, por Rodriguez y Bhattacharjee), se muestra aquí como una condición de solvabilidad natural para la existencia del invariante de Noether y la integrabilidad del sistema, sin la necesidad de multiplicadores de Lagrange externos.
• Resolución del Conundrum de la Presión Escalar: El marco sugiere que los modelos de presión escalar son restricciones estructuralmente inestables del estado cinético general. En geometrías 3D, la diferencia de presión anisotrópica $(p_\perp - p_\parallel)$ proporciona los grados de libertad necesarios para satisfacer las ecuaciones de equilibrio global, resolviendo así la naturaleza sobredeterminada del problema de la presión escalar.
• Criterio de Estabilidad Cinética: Utilizando el método de energía–Casimir, el artículo deriva una condición suficiente para la estabilidad de los equilibrios cuasisimétricos: $\partial f_{0s} / \partial H_{gc,s} < 0$. Esto generaliza la condición de monotonicidad para los equilibrios de Vlasov y sirve como un criterio de Newcomb cinético.

Significancia y Reivindicaciones
El artículo afirma resolver una brecha conceptual persistente al demostrar que la QS y la anisotropía de presión son manifestaciones duales de un único principio variacional cinético. Al derivar tanto las restricciones geométricas como el equilibrio anisotrópico de una acción unificada, los autores argumentan que la teoría del equilibrio anisotrópico es más fundamental que la MHD de presión escalar para describir sistemas cuasisimétricos.

El trabajo establece que las restricciones geométricas tradicionalmente impuestas "a mano" en el diseño de estelaradores emergen naturalmente como condiciones de solvabilidad para la existencia de invariantes de Noether. Esto implica que los estelaradores cuasisimétricos están caracterizados por corrientes paralelas suaves que evitan las singularidades encontradas en los modelos de presión escalar. Los autores concluyen que su marco proporciona un tratamiento unificado de la estabilidad cinética y de fluidos y ofrece una base para trabajos futuros en la construcción de resolvedores de equilibrio anisotrópico y la optimización de diseños de estelaradores basados en el grado de ruptura de simetría en la acción de centro de guía.