Synchrotron Emission from Cooled Particle Distributions

Resumo Técnico: Emissão de Síncrotron de Distribuições de Partículas Resfriadas

Enunciado do Problema
A emissão de síncrotron de elétrons relativísticos é uma ferramenta de diagnóstico primária para fontes astrofísicas de alta energia, incluindo surtos de raios gama (GRBs), transientes ópticos azuis rápidos luminosos (LFBOTs) e eventos de desrupção de maré (TDEs) com jatos. Nesses ambientes, os elétrons são acelerados em frentes de choque e, subsequentemente, resfriam-se via perdas radiativas (síncrotron e Compton inverso) e expansão adiabática. Esse resfriamento desloca os elétrons no espaço de energia, alterando fundamentalmente a função de distribuição de elétrons e, consequentemente, os coeficientes de emissão e absorção de síncrotron.

Embora modelos de "zona única" frequentemente assumam um resfriamento uniforme e distribuições simplificadas, modelos de "volume total" tratam a estrutura tridimensional do choque e do fluido a jusante de forma autoconsistente. Em modelos de volume total, a distribuição de elétrons a jusante torna-se altamente dependente do espaço devido a diferentes históricos de resfriamento. Abordagens anteriores para modelar esses efeitos locais basearam-se na integração numérica para computar a função de distribuição resfriada e os coeficientes de radiação associados em cada ponto atrás do choque. Essa dependência de métodos numéricos aumenta significativamente o custo computacional de simular esses sistemas, dificultando a eficiência da modelagem de volume total.

Metodologia
Os autores estendem o tratamento "local" da emissão de síncrotron derivando expressões analíticas para a função de distribuição de elétrons e os resultantes coeficientes de emissão e absorção, contabilizando tanto o resfriamento adiabático quanto o de síncrotron.

1. Formalismo Geral: O artigo inicia com a derivação da evolução de uma função de distribuição de elétrons injetada impulsivamente arbitrária, $(\partial n / \partial \gamma)_{\text{inj}}$, à medida que ela é advectada a jusante. Assumindo que o número total de elétrons em um elemento de fluido é conservado, os autores relacionam a distribuição no tempo $t$ à distribuição de injeção através da razão de densidades numéricas e do mapeamento dos fatores de Lorentz $\gamma_{\text{inj}} \to \gamma(t)$. A taxa de resfriamento inclui termos para perdas de síncrotron e expansão adiabática (trabalho de PdV), governados por escalas de tempo características $t_B$ e $t_{\text{ad}}$.
2. Casos Específicos: Os autores especializam este formalismo geral para duas distribuições injetadas:
   • Uma distribuição de lei de potência: $(\partial n / \partial \gamma)_{\text{inj, pl}} \propto \gamma^{-p}$.
   • Uma distribuição de Maxwell-Boltzmann relativística (térmica).
     Expressões analíticas para as funções de distribuição "resfriadas" são derivadas para ambos os casos. Essas expressões incluem fatores de corte proporcionais a $(1 - \gamma/\gamma_\infty)$, onde $\gamma_\infty$ representa o fator de Lorentz que um elétron injetado com energia infinita possuiria no tempo $t$.
3. Derivação de Funções de Ajuste: Os autores calculam a emissividade média pelo ângulo de pitch ($j_\nu$) e o coeficiente de absorção ($\alpha_\nu$) integrando as distribuições resfriadas contra o kernel de síncrotron. Como essas integrais não podem ser resolvidas em forma fechada para todas as regiões de frequência, os autores derivam funções de ajuste analíticas.
   • Lei de Potência: As funções de ajuste para o caso de lei de potência são construídas unindo limites assintóticos: baixas frequências ($x_1 \ll 1$), frequências intermediárias ($1 \le x_1 < \eta^2$) e altas frequências ($x_\infty \gg 1$). Os autores utilizam o método do declive mais íngreme para o regime de alta frequência e funções hipergeométricas para o regime de baixa frequência. A transição entre esses regimes é tratada por funções sigmoides com coeficientes ajustados a integrações numéricas para $2 < p \le 5$.
   • Térmica: Da mesma forma, funções de ajuste para a distribuição térmica são derivadas por interpolação entre limites de baixa frequência (envolvendo funções hipergeométricas confluentes) e limites de alta frequência (usando uma aproximação de declive mais íngreme para o ponto de sela da integral).

Principais Contribuições e Resultados

• Funções de Distribuição Analíticas: O artigo fornece formas analíticas exatas para as distribuições de elétrons de lei de potência resfriada (Eq. 11) e térmica (Eq. 12). Essas formas explicitamente contabilizam o "empilhamento" (pile-up) de elétrons próximo ao fator de Lorentz máximo $\gamma_\infty$ devido ao resfriamento, uma característica frequentemente negligenciada em modelos de corte abrupto (sharp cutoff).
• Funções de Ajuste Analíticas: A principal contribuição é o conjunto de funções de ajuste analíticas para os coeficientes de emissividade e absorção de síncrotron para distribuições de lei de potência e térmica.
  • Para distribuições de lei de potência, as funções de ajuste (Eqs. 33, 52) reproduzem integrações numéricas com erros relativos típicos na ordem de porcentagem, com erros máximos geralmente abaixo de 15% (atingindo o pico em $\eta \approx 1.14$).
  • Para distribuições térmicas, as funções de ajuste (Eqs. 79, 80) alcançam erros relativos médios $\lesssim 10\%$, com erros de pico de 10–25% ocorrendo perto da frequência de virada (turnover frequency), onde os assíntotas de baixa e alta frequência se encontram.
• Comparação com Modelos de Corte Abrupto: Os autores demonstram que aproximar a distribuição resfriada como uma lei de potência "decapitada" (corte abrupto em $\gamma_\infty$) introduz erros significativos, particularmente em altas frequências e quando o intervalo de resfriamento $\eta = \gamma_\infty/\gamma_1$ é pequeno. O corte suave derivado da ODE de resfriamento mostra-se necessário para uma modelagem local precisa.
• Aplicação a Afterglows de GRB: Os autores testam as funções de ajuste de lei de potência inserindo-as em um arcabouço numérico de volume total para afterglows de GRB (Ferguson & Margalit 2026) e comparando o resultado com o espectro analítico de Granot & Sari (2002). A comparação mostra que as novas funções de ajuste reproduzem o espectro padrão de "resfriamento lento" com erros relativos consistentemente abaixo de 10% (tipicamente $\sim 1\%$), validando seu uso como substituto para integrações numéricas computacionalmente caras.

Significância e Alegações
O artigo alega que estas funções de ajuste analíticas oferecem uma alternativa computacionalmente eficiente e precisa à integração numérica para calcular os coeficientes de radiação de síncrotron em modelos de volume total. Ao evitar a necessidade de resolver as integrais de resfriamento numericamente em cada ponto espacial, essas funções permitem o tratamento consistente do resfriamento nos espectros de síncrotron sem o custo computacional proibitivo dos métodos anteriores.

Os autores enfatizam que, embora modelos de zona única frequentemente aproximem o resfriamento impondo um índice espectral mais íngreme, modelos de volume total requerem a média das contribuições de elementos de fluido com diferentes graus de resfriamento. As funções de ajuste derivadas permitem esse processo de média, facilitando a calibração de modelos de zona única globais contra propriedades de radiação locais mais precisas. Os autores observam que seu formalismo é aplicável a qualquer população de elétrons injetada impulsivamente que resfrie ao longo do tempo, embora restrinjam explicitamente sua análise atual a regimes onde os elétrons permanecem relativísticos ($\gamma \gg 1$) e o índice de lei de potência é $p > 2$. Eles reconhecem que efeitos como o resfriamento por Compton inverso podem ser incluídos trivialmente, enquanto outros processos como o resfriamento de Coulomb ou o aquecimento por autoabsorção exigiriam modificações adicionais ao formalismo.