Synchrotron Emission from Cooled Particle Distributions

Resumen Técnico: Emisión de Sincrotrón de Distribuciones de Partículas Enfriadas

Planteamiento del Problema
La emisión de sincrotrón de electrones relativistas es una herramienta de diagnóstico primaria para fuentes astrofísicas de alta energía, incluyendo los brotes de rayos gamma (GRB), los transitorios ópticos azules rápidos y luminosos (LFBOT) y los eventos de disrupción de marea (TDE) con chorros. En estos entornos, los electrones son acelerados en frentes de choque y posteriormente se enfrían mediante pérdidas radiativas (sincrotrón e Compton inverso) y expansión adiabática. Este enfriamiento desplaza a los electrones en el espacio de energía, alterando fundamentalmente la función de distribución de electrones y, consecuentemente, los coeficientes de emisión y absorción de sincrotrón.

Mientras que los modelos de "una sola zona" a menudo asumen un enfriamiento uniforme y distribuciones simplificadas, los modelos de "volumen completo" tratan la estructura tridimensional del choque y del fluido aguas abajo de manera autoconsistente. En los modelos de volumen completo, la distribución de electrones aguas abajo se vuelve altamente dependiente del espacio debido a las variaciones en las historias de enfriamiento. Los enfoques previos para modelar estos efectos locales han dependido de la integración numérica para calcular la función de distribución enfriada y los coeficientes de radiación asociados en cada punto detrás del choque. Esta dependencia de los métodos numéricos aumenta significativamente el costo computacional de simular estos sistemas, obstaculizando la eficiencia del modelado de volumen completo.

Metodología
Los autores extienden el tratamiento "local" de la emisión de sincrotrón derivando expresiones analíticas para la función de distribución de electrones y los coeficientes de emisión y absorción resultantes, contabilizando tanto el enfriamiento adiabático como el de sincrotrón.

1. Formalismo General: El artículo comienza derivando la evolución de una función de distribución de electrones inyectada impulsivamente arbitraria, $(\partial n / \partial \gamma)_{\text{inj}}$, a medida que es advectada aguas abajo. Asumiendo que el número total de electrones en un elemento de fluido se conserva, los autores relacionan la distribución en el tiempo $t$ con la distribución de inyección mediante la relación de densidades numéricas y el mapeo de los factores de Lorentz $\gamma_{\text{inj}} \to \gamma(t)$. La tasa de enfriamiento incluye términos para las pérdidas por sincrotrón y la expansión adiabática (trabajo de $PdV$), gobernados por escalas de tiempo características $t_B$ y $t_{\text{ad}}$.
2. Casos Específicos: Los autores especializan este formalismo general a dos distribuciones inyectadas:
   • Una distribución de ley de potencia: $(\partial n / \partial \gamma)_{\text{inj, pl}} \propto \gamma^{-p}$.
   • Una distribución de Maxwell-Boltzmann relativista (térmica).
     Se derivan expresiones analíticas para las funciones de distribución "enfriadas" para ambos casos. Estas expresiones incluyen factores de corte proporcionales a $(1 - \gamma/\gamma_\infty)$, donde $\gamma_\infty$ representa el factor de Lorentz que poseería un electrón inyectado con energía infinita en el tiempo $t$.
3. Derivación de Funciones de Ajuste: Los autores calculan la emisividad promedio respecto al ángulo de pitch ($j_\nu$) y el coeficiente de absorción ($\alpha_\nu$) integrando las distribuciones enfriadas contra el núcleo de sincrotrón. Debido a que estas integrales no pueden resolverse en forma cerrada para todos los regímenes de frecuencia, los autores derivan funciones de ajuste analíticas.
   • Ley de Potencia: Las funciones de ajuste para el caso de la ley de potencia se construyen uniendo los límites asintóticos: bajas frecuencias ($x_1 \ll 1$), frecuencias intermedias ($1 \le x_1 < \eta^2$) y altas frecuencias ($x_\infty \gg 1$). Los autores utilizan el método del descenso más pronunciado para el régimen de alta frecuencia y funciones hipergeométricas para el régimen de baja frecuencia. La transición entre estos regímenes se maneja mediante funciones sigmoides con coeficientes ajustados a integraciones numéricas para $2 < p \le 5$.
   • Térmica: De manera similar, se derivan funciones de ajuste para la distribución térmica mediante la interpolación entre los límites de baja frecuencia (que involucran funciones hipergeométricas confluentes) y los límites de alta frecuencia (usando una aproximación de descenso más pronunciado para el punto de silla de la integral).

Contribuciones Clave y Resultados

• Funciones de Distribución Analíticas: El artículo proporciona formas analíticas exactas para las distribuciones de electrones de ley de potencia enfriada (Eq. 11) y térmica (Eq. 12). Estas formas contabilizan explícitamente el "acumulamiento" (pile-up) de electrones cerca del factor de Lorentz máximo $\gamma_\infty$ debido al enfriamiento, una característica que a menudo se omite en los modelos de corte abrupto.
• Funciones de Ajuste Analíticas: La principal contribución es el conjunto de funciones de ajuste analíticas para los coeficientes de emisividad y absorción de sincrotrón para las distribuciones de ley de potencia y térmica.
  • Para las distribuciones de ley de potencia, las funciones de ajuste (Ecs. 33, 52) reproducen las integraciones numéricas con errores relativos típicos a nivel de porcentaje, con errores máximos generalmente por debajo del 15% (alcanzando su pico alrededor de $\eta \approx 1.14$).
  • Para las distribuciones térmicas, las funciones de ajuste (Ecs. 79, 80) logran errores relativos medios $\lesssim 10\%$, con errores pico de 10–25% que ocurren cerca de la frecuencia de giro donde se encuentran los asíntotas de baja y alta frecuencia.
• Comparación con Modelos de Corte Abrupto: Los autores demuestran que aproximar la distribución enfriada como una ley de potencia "recortada" (corte abrupto en $\gamma_\infty$) introduce errores significativos, particularmente en altas frecuencias y cuando el rango de enfriamiento $\eta = \gamma_\infty/\gamma_1$ es pequeño. Se muestra que el corte suave derivado de la ODE de enfriamiento es necesario para un modelado local preciso.
• Aplicación a los Afterglows de GRB: Los autores prueban las funciones de ajuste de la ley de potencia insertándolas en un marco numérico de volumen completo para los afterglows de GRB (Ferguson & Margalit 2026) y comparando los resultados con el espectro analítico de Granot & Sari (2002). La comparación muestra que las nuevas funciones de ajuste reproducen el espectro estándar de "enfriamiento lento" con errores relativos consistentemente por debajo del 10% (típicamente $\sim 1\%$), validando su uso como sustituto de las integrales numéricas computacionalmente costosas.

Significado y Reivindicaciones
El artículo afirma que estas funciones de ajuste analíticas ofrecen una alternativa computacionalmente eficiente y precisa a la integración numérica para calcular los coeficientes de radiación de sincrotrón en modelos de volumen completo. Al evitar la necesidad de resolver numéricamente las integrales de enfriamiento en cada punto espacial, estas funciones permiten el tratamiento consistente del enfriamiento en los espectros de sincrotrón sin el costo computacional prohibitivo de los métodos anteriores.

Los autores enfatizan que, mientras que los modelos de una sola zona a menudo aproximan el enfriamiento imponiendo un índice espectral más pronunciado, los modelos de volumen completo requieren el promedio de las contribuciones de elementos de fluido con distintos grados de enfriamiento. Las funciones de ajuste derivadas permiten este proceso de promediado, facilitando la calibración de los modelos globales de una sola zona frente a las propiedades de radiación locales más precisas. Los autores señalan que su formalismo es aplicable a cualquier población de electrones inyectada impulsivamente que se enfríe con el tiempo, aunque restringen explícitamente su análisis actual a regímenes donde los electrones permanecen relativistas ($\gamma \gg 1$) y el índice de la ley de potencia es $p > 2$. Reconocen que los efectos como el enfriamiento por Compton inverso pueden incluirse trivialmente, mientras que otros procesos como el enfriamiento de Coulomb o el calentamiento por autoabsorción requerirían modificaciones adicionales al formalismo.