Time-dependent cosmic-ray escape from wind bubbles: hard spectra formation

Resumen Técnico: Escape de rayos cósmicos dependiente del tiempo de las burbujas de viento: formación de espectros duros

Planteamiento del Problema
Las burbujas impulsadas por viento, generadas por fuentes tales como cúmulos estelares masivos jóvenes (YMSCs), estrellas de Wolf–Rayet y objetos compactos con acreción super-Eddington, son candidatas prometedoras para la aceleración de rayos cósmicos (CRs) hasta energías de PeV (PeVatrons). Si bien la aceleración difusiva en choques (DSA) en el choque de terminación del viento interno es ampliamente aceptada como el principal mecanismo de aceleración, el transporte dependiente del tiempo de estas partículas desde el choque de terminación hacia el choque frontal (el límite de escape hacia el medio interestelar) permanece poco explorado. Estudios previos a menudo dependieron de descripciones estacionarias o estimaciones de escalas de tiempo de transporte en lugar de resolver la ecuación de transporte completa y dependiente del tiempo. Este trabajo aborda la brecha en la comprensión de cómo la expansión dinámica de la burbuja y la evolución temporal de la distancia entre el choque de terminación y el choque frontal impactan en el espectro de energía de los rayos cósmicos que escapan.

Metodología
Los autores modelan la burbuja de viento como un sistema con simetría esférica donde tanto el choque de terminación (TS) como el choque frontal (FS) se expanden con velocidades dependientes del tiempo ($R_{ts} \propto t^{2/5}$ y $R_{fs} \propto t^{3/5}$ durante la fase de deceleración). El transporte de los CRs está gobernado por una ecuación de transporte esférica unidimensional dependiente del tiempo (ecuación de Fokker-Planck), que incluye la advección por el flujo del viento sacudido y la difusión.

Para resolver esta ecuación, los autores la transforman en un conjunto de ecuaciones diferenciales estocásticas (SDEs) utilizando el léma de Itô. Estas SDEs se integran utilizando una versión modificada del marco de trabajo de código abierto CRPropa3.2. Las características metodológicas clave incluyen:

• Inyección: Las partículas se inyectan en el TS con un espectro de ley de potencia pre-acelerado ($dN/dE \propto E^{-2}$) y un corte de alta energía dependiente del tiempo determinado por el límite de Hillas y las escalas de tiempo de aceleración.
• Advección: Se asume un perfil de velocidad radial $u(r) \propto r^{-2}$ en la región aguas abajo, normalizado a la velocidad post-choque.
• Difusión: El estudio investiga tres modelos de turbulencia: Bohm ($\delta=1$), Kraichnan ($\delta=3/2$) y Kolmogorov ($\delta=5/3$). El coeficiente de difusión depende de la fuerza del campo magnético local y del espectro de potencia de la turbulencia.
• Campo Magnético: El campo magnético se estima basándose en una fracción de la presión de choque. El estudio prueba varias dependencias radiales del campo magnético aguas abajo ($B \propto r^{-\alpha}$ con $\alpha = 0, -1, -2$) y niveles de magnetización ($\epsilon_B$).
• Parámetros de Simulación: Las simulaciones corren hasta los 5 Myr, cubriendo la fase impulsada por el viento antes de que las supernovas dominen el presupuesto de energía.

Resultados Clave
El tratamiento dependiente del tiempo revela desviaciones significativas de la expectativa estándar de los espectros de escape:

1. Endurecimiento Espectral: En configuraciones físicas específicas, los espectros de CRs que escapan pueden ser más duros que el espectro inyectado de $E^{-2}$. Este endurecimiento es más pronunciado para modelos de turbulencia con una dependencia de energía más fuerte (por ejemplo, difusión de Bohm) y dependencias de campo magnético más débiles.
2. Supresión de Baja Energía: Un hallazgo crítico es la supresión significativa de partículas de baja energía en el flujo de escape. Debido a que el choque frontal (límite de escape) se aleja de la fuente con el tiempo, las partículas de baja energía con coeficientes de difusión pequeños no pueden alcanzar el límite de escape y permanecen atrapadas en la región aguas abajo.
   • Para la difusión de Bohm, las partículas de baja energía ($E < 1$ TeV) están casi completamente depletadas del flujo de escape.
   • Para la difusión de Kraichnan, se observa un endurecimiento moderado ($\Delta \gamma \approx 0.1-0.2$) y una supresión parcial de la baja energía.
   • Para la difusión de Kolmogorov, el depletamiento es mínimo y el espectro de escape se asemeja estrechamente al espectro inyectado, aunque la energía máxima sigue siendo menor ($\sim 20$ TeV).
3. Dependencia de Parámetros: El grado de depletamiento de baja energía y el endurecimiento espectral dependen del modelo de turbulencia, la dependencia radial del campo magnético (una dependencia más débil conduce a espectros más duros) y el nivel de magnetización.
4. Implicaciones Observacionales: Al propagarse hacia la Tierra (simulado para una fuente a 3 kpc), la supresión de baja energía persiste, aunque el endurecimiento es parcialmente suavizado por la difusión galáctica en escalas de tiempo largas.

Significancia y Reivindicaciones
El artículo afirma que la naturaleza dependiente del tiempo de la expansión de la burbuja de viento es un factor crucial en la conformación del espectro de CR observado, independientemente de los procesos de pérdida de energía. Los autores enfatizan que:

• Formación de Espectros Duros: Las burbujas de viento pueden producir naturalmente espectros de CR duros (más duros que $E^{-2}$) sin requerir física de aceleración no estándar, siempre que el transporte sea dependiente del tiempo.
• Corte de Baja Energía: El "atrapamiento" de partículas de baja energía dentro de la burbuja es una predicción robusta de este modelo dependiente del tiempo. Esto podría conducir a características observables en la radiación multi-mensajero (rayos gamma y neutrinos) desde el interior de la burbuja y afectar la gramática acumulada de los CRs.
• Distinción de los SNR: Los espectros duros resultantes y los cortes de baja energía difieren de los típicamente esperados de los remanentes de supernova (SNR), lo que potencialmente ayuda a explicar el espectro de CR alrededor de la "rodilla" (knee) o observaciones específicas de rayos gamma de cúmulos estelares (por ejemplo, Westerlund 1).

Los autores concluyen que, si bien su modelo se centra en la fase impulsada por el viento, estas partículas de baja energía atrapadas pueden ser re-aceleradas posteriormente por choques de supernovas o escapar cuando el choque frontal cesa, contribuyendo potencialmente a la población de CR de baja energía con una forma espectral modificada. El trabajo futuro incorporará interacciones hadrónicas para predecir las emisiones de rayos gamma y neutrinos dependientes del tiempo.