Kinetic Simulations of Laser-Driven Compression and Heating of Magnetised Cryogenic Hydrogen Targets using PIConGPU

Este artículo presenta simulaciones totalmente cinéticas con PIConGPU que demuestran que la compresión impulsada por láser de objetivos de hidrógeno criogénico magnetizado produce un mecanismo dominante de aceleración iónica no térmica mediante capas dobles de separación de carga, el cual se mantiene robusto bajo campos magnéticos a escala de laboratorio pero es significativamente suprimido y alterado por campos a escala de kilotesla que magnetizan electrones calientes y extienden los tiempos de compresión.

Lo esencial

La Gran Imagen: Aplastando un Cubito de Hielo Minúsculo con Luz

Imagina que tienes un pequeño cilindro de hidrógeno congelado (como un cubito de hielo microscópico) y quieres aplastarlo en el centro para crear una presión extrema. Para lograrlo, los científicos están utilizando láseres potentes. Este artículo es una simulación por computadora que actúa como un "túnel de viento digital", prediciendo exactamente qué sucede cuando estos láseres golpean el hielo.

Los investigadores están probando dos tipos diferentes de pulsos láser:

1. El "Aplauso" (30 femtosegundos): Un estallido de energía súper rápido y agudo, como un martillo golpeando un clavo.
2. El "Empujón" (150 femtosegundos): Un empuje más largo y sostenido, como una mano que presiona firmemente hacia abajo sobre un resorte de manera lenta pero constante.

También probaron qué sucede si añaden un campo magnético gigante a la mezcla, actuando como una jaula invisible alrededor del hielo.

El Descubrimiento Principal: Dos Tipos de Partículas

Cuando los láseres golpean el hidrógeno, no solo lo calientan; crean un extraño "atascos de tráfico" de partículas. La simulación reveló que el hidrógeno se divide en dos grupos distintos, muy parecido a una multitud de personas reaccionando a un evento repentino:

1. Los "Velocistas" (Iones Rápidos): Un pequeño grupo de partículas recibe un fuerte golpe y se lanza hacia adentro a velocidades increíbles (millones de electronvoltios).
2. Los "Caminantes" (Iones Masivos): El resto de las partículas se mueven hacia adentro mucho más lentamente, como una multitud que avanza a paso lento.

La Analogía del "Espejo Mágico":
El artículo explica que los "Velocistas" no son empujados directamente por el láser. En cambio, el láser crea una pared móvil de carga eléctrica (un "frente de separación de cargas") que actúa como un espejo móvil.

• Cuando el láser golpea el hielo, empuja los electrones hacia afuera, dejando un vacío.
• Este vacío crea un campo eléctrico masivo (¡unos 3 billones de voltios por metro!).
• A medida que este "espejo" eléctrico se mueve hacia adentro, rebota los iones de hidrógeno positivos contra él.
• Al igual que una pelota de tenis que rebota en una raqueta que se mueve hacia ti, los iones ganan velocidad. El artículo encontró una regla simple: si el espejo se mueve a velocidad $v$, la pelota rebota de vuelta a velocidad $2v$.

La Diferencia entre el "Aplauso" y el "Empujón"

El tipo de pulso láser cambia cómo se comportan estos "Velocistas":

• El "Aplauso" (30 fs): Debido a que el láser es tan corto, el espejo eléctrico se mueve a una velocidad constante durante una fracción de segundo. Esto crea un grupo ordenado y uniforme de Velocistas, todos golpeando el centro con exactamente la misma velocidad. Es como una salva de flechas perfectamente sincronizada.
• El "Empujón" (150 fs): Debido a que el láser dura más tiempo, el espejo eléctrico sigue acelerando mientras se mueve. Esto significa que los Velocistas son lanzados a diferentes velocidades a lo largo del tiempo. Algunos son lentos, otros rápidos. Es como un chorro de agua donde la velocidad varía, creando un "barrido" de energías en lugar de un grupo único y agudo.

El Experimento del Campo Magnético: La Jaula Invisible

Los investigadores luego activaron un campo magnético para ver si atraparía las partículas y aplastaría el hielo con más fuerza. Probaron campos que van desde lo que podemos construir en un laboratorio (20 Tesla) hasta campos teóricos extremos (10,000 Tesla).

• El Campo a Escala de Laboratorio (20 T): Esto es como una brisa suave. Las partículas se mueven tan rápido y son tan energéticas que simplemente ignoran el campo magnético. Atraviesan el campo sin problemas. La simulación mostró cero cambio en los resultados.
• El Campo Extremo (1,000–10,000 T): Esto es como una jaula de acero. A estos niveles, el campo magnético es lo suficientemente fuerte como para atrapar a los electrones en movimiento rápido.
  • El Resultado: Cuando los electrones quedan atrapados, no pueden escapar para formar ese "espejo móvil" más. Sin el espejo, los "Velocistas" (los iones rápidos) desaparecen. El láser pierde su capacidad de patear los iones hacia adentro.
  • El Giro: Aunque los "Velocistas" se han ido, el campo magnético en realidad ayuda a que los "Caminantes" (los iones masivos) permanezcan comprimidos durante el doble de tiempo. Es como si la jaula magnética mantuviera la presión por más tiempo, permitiendo que la multitud en movimiento lento aplaste el centro de manera más efectiva antes de rebotar hacia afuera.

Un Efecto Secundario Sorprendente: El Efecto del Globo

Podrías pensar que una jaula magnética apretaría todo más fuerte. Sin embargo, la simulación mostró algo contra intuitivo: el borde exterior del objetivo de hidrógeno en realidad se expandió más cuando el campo magnético fue fuerte.

La Analogía: Imagina un globo. Si aprietas el medio, los extremos pueden hincharse. El campo magnético atrapa los electrones calientes, pero también cambia la forma en que empujan contra las capas externas del objetivo. En lugar de colapsar ordenadamente, la "piel" exterior del objetivo se hincha más hacia el espacio.

El "Truco Geométrico"

El artículo señala una forma inteligente de probar esto en el mundo real. Los campos de 10,000 Tesla utilizados en la simulación son imposibles de construir para un objetivo diminuto de 15 micras. Sin embargo, la física depende de la relación entre la trayectoria de la partícula y el tamaño del objetivo.

Los autores argumentan que si usas un objetivo mucho más grande (como un chorro de hidrógeno 1,000 veces más grande), no necesitarías 10,000 Tesla. Podrías usar un campo modesto de 10 Tesla (que es fácil de construir) y obtener exactamente el mismo efecto de atrapamiento magnético. Es como cómo un auto de juguete pequeño y un auto real pueden girar de la misma manera si ajustas la velocidad del volante en relación con su tamaño.

Resumen

• Los láseres crean una pared eléctrica móvil que rebota los iones hacia adentro.
• Los láseres cortos crean un grupo uniforme de iones rápidos; los láseres largos crean un grupo mixto.
• Los imanes débiles no hacen nada.
• Los imanes súper fuertes detienen los iones rápidos pero ayudan a que los iones lentos permanezcan comprimidos por más tiempo.
• Los imanes fuertes también hacen que el borde exterior del objetivo se hinche, en lugar de encogerse.
• Los objetivos grandes pueden experimentar estos efectos de "superimán" usando imanes normales, de tamaño de laboratorio.

Resumen técnico

Resumen Técnico: Simulaciones Cinéticas de Compresión y Calentamiento Impulsados por Láser de Blancos de Hidrógeno Criogénico Magnetizados utilizando PIConGPU

Enunciado del Problema
El artículo aborda los mecanismos fundamentales que gobiernan el acoplamiento láser-plasma, el transporte de electrones calientes y la aceleración de iones en el régimen de alta densidad de energía, centrándose específicamente en las interacciones con blancos de hidrógeno criogénico de densidad sólida. Los modelos hidrodinámicos de radiación estándar (rad-hid), que tratan al plasma como un fluido cuasi-neutral, no logran capturar estructuras electrostáticas intrínsecamente no cuasi-neutrales ni canales de partículas no térmicas que surgen en escalas de tiempo sub-picosegundo. Los autores buscan establecer una base numérica predictiva para las próximas campañas experimentales en el Helmholtz-Zentrum Dresden-Rossendorf (HZDR), contrastando el láser operativo DRACO ($\tau = 30$ fs) con la futura instalación PEnELOPE ($\tau = 150$ fs). Además, el estudio investiga la influencia de campos magnéticos axiales externos sobre estas dinámicas, buscando comprender los umbrales de magnetización y su impacto en la compresión y la aceleración de iones.

Metodología
Los autores emplearon simulaciones totalmente cinéticas, bidimensionales, con tres componentes de velocidad (2D3V) de Partículas en Celda (PIC) utilizando el código de código abierto PIConGPU.

• Blanco: Un cilindro de hidrógeno criogénico de densidad sólida ($R = 7.5$ µm, $n_0 = 5.24 \times 10^{28}$ m$^{-3}$).
• Configuración del Láser: Una geometría de impulso directo simétrica de tres haces con polarización lineal ($\lambda = 800$ nm). Se modelaron dos duraciones de pulso: $\tau = 30$ fs (impulsivo) y $\tau = 150$ fs (sostenido), con potenciales vectoriales normalizados de $a_0 = 12.7$ y $a_0 = 22.0$.
• Módulos de Física: Las simulaciones utilizaron el módulo FLYonPIC para cinética atómica no-LTE autoconsistente e ionización por campo/colisional. El enfoque totalmente cinético resolvió explícitamente las ecuaciones de Maxwell y el sistema Vlasov-Maxwell, evitando cierres de flujo de calor fluido.
• Escaneo de Campo Magnético: Se aplicó un campo magnético axial estático externo ($B_z$), escaneando desde 0 T hasta 10 kT. Aunque 20 T representa un campo inicial alcanzable en laboratorio, los campos a escala de kT se utilizaron como proxy para la física magnetizada en blancos más grandes mediante un argumento de equivalencia geométrica (manteniendo constante la relación entre el radio de Larmor del electrón y el radio del blanco).
• Diagnósticos: Las métricas clave incluyeron campos eléctricos espacio-temporales ($E_r$), espectros de energía de iones, relaciones de momento hacia el interior para definir el tiempo de compresión ($t_{comp}$) y relaciones de compresión locales en $r=1$ µm.

Contribuciones y Resultados Clave

1. Bifurcación Cinemática y la Ley de Reflexión $2v_{hb}$:
   Las simulaciones revelan una bifurcación robusta de la población de iones en un haz rápido (1–5 MeV) y un flujo de masa más lento (1–100 keV). Los autores demuestran que la población de iones rápidos se origina en la reflexión directa de una capa doble electrostática en movimiento y no cuasi-neutral (la frente de separación de cargas).

   • Una escala de reflexión simple, $E_k \approx 2v_{hb}$ (donde $v_{hb}$ es la velocidad instantánea de la frente de separación de cargas), predice cuantitativamente la energía de los iones rápidos.
   • Para el impulsor de 30 fs, la frente se mueve a una velocidad casi constante, produciendo una banda de iones rápidos monoenergética.
   • Para el impulsor de 150 fs, la frente acelera continuamente durante el ascenso del pulso, barrido la energía de los iones rápidos desde cientos de keV hasta varios MeV.
   • Este mecanismo se identifica como una característica genuinamente cinética, que queda fuera de las suposiciones de cierre de los modelos rad-hid que no pueden resolver la estructura coherente de la capa doble.

2. Umbrales y Efectos del Campo Magnético:

   • Régimen de 20 T: Los campos alcanzables en laboratorios actuales (20 T) tienen un impacto insignificante en los observables macroscópicos porque los radios de Larmor de los electrones calientes de MeV y los protones permanecen órdenes de magnitud mayores que el radio del blanco.
   • Régimen a Escala de kT: Los campos en el rango de kilotesla magnetizan progresivamente a la población de electrones calientes de MeV ($r_L < R$). Esto conduce a:
     • Apagado de la Aceleración de Iones: El mecanismo de separación de cargas impulsado por láser se suprime a medida que los electrones calientes se magnetizan y se fijan a las líneas de campo, impidiendo la formación de la fuerte capa doble electrostática dirigida hacia el interior. En consecuencia, la banda de iones rápidos de MeV desaparece de los espectros.
     • Compresión Extendida: El tiempo de compresión neta hacia el interior ($t_{comp}$) para el impulsor de 30 fs se más que duplica bajo una fuerte magnetización (por ejemplo, de ~0.52 ps a ~1.1 ps en $a_0=12.7$). El empuje hacia el interior sostenido permite que la onda de iones de masa llegue al núcleo, creando un segundo pico de densidad más grande que está ausente en el caso no magnetizado.
     • Expansión de la Envoltura: Contrariamente a las expectativas ingenuas de confinamiento, la envoltura externa del blanco se expande aún más con el aumento de $B_z$. Esto se atribuye a electrones calientes fijados magnéticamente que sostienen un campo ambipolar hacia el exterior en tiempos tardíos, expandiendo la envoltura iónica incluso mientras el núcleo se comprime.

3. Equivalencia Geométrica:
   Los autores argumentan que la física a escala de kT observada en el blanco de 15 µm es cinemáticamente equivalente a intensidades de campo moderadas (~10 T) aplicadas a chorros de hidrógeno criogénico mucho más grandes (por ejemplo, 7.5 mm de diámetro), basándose en la conservación de la relación $r_L/R$. Esto sugiere que la física reportada es accesible para experimentos a mayor escala, aunque no se mantiene la similitud exacta del plasma (escalar la densidad como $1/R^2$).

Significado y Afirmaciones
El artículo afirma dos hallazgos centrales:

1. Perspectiva Mecanicista: La aceleración de iones en la superficie frontal en este régimen está impulsada por una capa doble electrostática no cuasi-neutral cuya velocidad instantánea establece las energías de los iones rápidos mediante una ley de reflexión $2v_{hb}$. Este mecanismo es distinto de la Aceleración por Manto Normal al Blanco (TNSA) y requiere un tratamiento totalmente cinético para ser capturado.
2. Impacto de la Magnetización: Un campo magnético axial aplicado, una vez que reduce el radio de Larmor de los electrones calientes de MeV por debajo del radio del blanco, suprime coherentemente este canal de aceleración. Si bien esto suprime el haz de iones rápidos, extiende la fase de compresión neta hacia el interior y altera la dinámica de la envoltura del blanco.

El estudio enfatiza que estos hallazgos proporcionan una base predictiva para las campañas de DRACO y PEnELOPE. Los autores señalan limitaciones inherentes a la geometría 2D3V (supresión de corrientes de retorno longitudinales) y la falta de similitud exacta del plasma en su mapeo geométrico, advirtiendo que los valores absolutos de campo y las relaciones de compresión deben interpretarse como indicativos en lugar de directamente escalables a experimentos 3D sin una investigación adicional. Se propone trabajo futuro para incluir simulaciones 3D, estudios de campos auto-generados mediante incidencia inclinada y comparación directa con observables experimentales.