Population synthesis of Galactic middle-aged pulsar wind nebulae I. Detection prospects for current and future instruments

Declaração do Problema
As nebulosas de vento de pulsar (PWNe) representam a maior população de fontes galácticas de raios gama de energia muito alta (VHE; $E > 100$ GeV) e servem como laboratórios críticos para o estudo da aceleração de partículas e das interações entre pulsares e remanescentes de supernova (SNR). Embora instrumentos atuais como H.E.S.S., MAGIC e VERITAS tenham identificado firmemente inúmeras PWNe, previsões realistas ao nível da população historicamente careceram de um tratamento detalhado da "fase de reverberação". Esta fase ocorre quando o choque reverso do SNR viaja de volta em direção ao centro da explosão, colidindo com e comprimindo a PWN. Este evento altera significativamente a dinâmica da nebulosa, amplifica os campos magnéticos, reenergiza as populações de partículas e aumenta drasticamente as perdas radiativas. Abordagens de modelagem anteriores enfrentaram uma dicotomia: modelos simples de casca fina de uma zona são computacionalmente eficientes, mas falham em capturar o feedback complexo e não linear da fase de reverberação, enquanto simulações hidrodinâmicas (HD) e magnetohidrodinâmicas (MHD) multidimensionais são precisas, mas proibitivamente caras para estudos estatísticos envolvendo milhares de fontes ao longo de $10^5$ anos. Consequentemente, previsões para a população galáctica de PWNe em TeV, particularmente para observatórios futuros, correm o risco de ser enviesadas se esta fase for negligenciada ou simplificada em excesso.

Metodologia
Para abordar essas limitações, os autores empregam o framework híbrido TIDE+L, que combina a eficiência computacional de um modelo dinâmico de casca fina (TIDE) para a fase de expansão livre com um tratamento lagrangiano da estrutura do SNR durante a reverberação. Esta abordagem permite a evolução autoconsistente do sistema PWN–SNR em todos os estágios, incluindo os complexos episódios de compressão da reverberação, ao mesmo tempo que rastreia o espectro de partículas sujeito a perdas radiativas e adiabáticas.

O estudo constrói uma população sintética de aproximadamente 2.400 PWNe galácticas, das quais 1.600 fontes são amostradas aleatoriamente em 1.000 realizações para garantir robustez estatística. A população é amostrada a partir de distribuições de:

• Propriedades do pulsar: Período de rotação inicial (distribuição gaussiana, truncada em 10 ms), campo magnético inicial (log-normal) e índice de frenagem ($n=2.33$).
• Propriedades do SNR e ambientais: Energia da explosão (gaussiana truncada), massa do ejetado (derivada da massa da sequência principal de idade zero via modelos de evolução estelar GENEVA) e densidade do meio interestelar (MIS) (distribuição uniforme).
• Parâmetros espectrais: Espectros de injeção de partículas, fatores de contenção e frações magnéticas.

Cada sistema é evoluído até uma idade máxima de $10^5$ anos. Os fluxos de raios gama resultantes são então comparados com os limiares de sensibilidade e a cobertura do céu de instrumentos atuais e futuros: H.E.S.S., HAWC, o Observatório do Telescópio Cherenkov (CTAO), LHAASO e o Observatório Sul de Raios Gama de Grande Campo (SWGO). A análise leva em conta efeitos de extensão da fonte, que degradam a sensibilidade para fontes estendidas.

Principais Contribuições

1. Modelagem Realista de Reverberação: O artigo introduz uma síntese de população estatisticamente robusta que incorpora explicitamente o tratamento lagrangiano da reverberação, indo além das prescrições analíticas simplificadas de modelos anteriores de casca fina.
2. Aplicação do Framework Híbrido: Demonstra a aplicação do framework TIDE+L a um estudo de população em grande escala, preenchendo a lacuna entre a física hidrodinâmica detalhada e os requisitos computacionais da síntese de populações.
3. Perspectivas de Detecção: O estudo fornece as primeiras estimativas detalhadas da população detectável de PWNe galácticas para instrumentos de próxima geração, quantificando especificamente o rendimento do CTAO, SWGO e LHAASO.
4. Estágios Evolutivos: Categoriza a população em fases de expansão livre, compressão e pós-compressão, analisando a detectabilidade e as características espectrais de cada uma.

Resultados

• Detectabilidade: O modelo prevê que o CTAO detectará uma ordem de magnitude mais PWNe do que aquelas atualmente firmemente detectadas na faixa de TeV. Especificamente, espera-se que o CTAO detecte aproximadamente 172 fontes (10,7% da população), subindo para 226 (14,1%) com uma configuração atualizada "CTAO+". Isso é significativamente maior do que as ~52 fontes previstas para o H.E.S.S.
• Impacto da Reverberação: Comparações entre o framework TIDE+L e o modelo TIDE simplificado revelam que negligenciar a dinâmica realista da reverberação leva a uma subestimação das fontes detectáveis. O modelo simplificado tende a superestimar a compressão e a amplificação do campo magnético, resultando em perdas de sincrotron excessivas e em uma população de partículas de alta energia esgotada, reduzindo assim o fluxo de raios gama previsto.
• Distribuição Evolutiva: A maioria da população sintética (87,5%) está na fase pós-compressão. No entanto, fontes nesta fase são tipicamente mais antigas, maiores e intrinsecamente mais fracas, resultando em frações de detectabilidade mais baixas em comparação com aquelas nas fases de expansão livre ou compressão.
• Evolução Espectral: Estudos de caso individuais ilustram que o fator de compressão (CF) governa a força e a duração da amplificação da luminosidade e do endurecimento/amolecimento espectral durante a reverberação. Fontes de alto CF exibem mudanças rápidas e acentuadas na distribuição de energia espectral (SED), enquanto fontes de baixo CF mostram uma evolução mais suave e atrasada.
• Correlações: O estudo encontra uma correlação entre a luminosidade integrada de 1–10 TeV e a potência de frenagem do pulsar ($L_{1-10 \text{ TeV}} \propto \dot{E}^{1.15}$). Os autores observam que recuperar a inclinação mais rasa observada nos dados do HGPS ($L \propto \dot{E}^{0.59}$) requer a aplicação de cortes de seleção específicos (por exemplo, excluindo fontes com baixa potência de frenagem ou idades mais avançadas), sugerindo que vieses de levantamento e limitações instrumentais influenciam as correlações observadas.
• Extensão da Fonte: As extensões angulares simuladas são geralmente consistentes com os dados do HGPS, embora o modelo preveja uma fração maior de fontes com raios menores que o mínimo observado no HGPS. Isso é atribuído à suposição de simetria esférica e à negligência da mistura e do movimento próprio do pulsar.

Significância
O artigo afirma que a modelagem realista da fase de reverberação é essencial para prever com precisão a população galáctica de PWNe em TeV. Ao incorporar a dinâmica complexa das interações do choque reverso do SNR, o estudo fornece uma base teórica mais confiável para interpretar catálogos de fontes atuais e antecipar os rendimentos de levantamentos futuros. Os resultados confirmam que o CTAO desempenhará um papel dominante no censo populacional de TeV que se aproxima, potencialmente resolvendo a natureza de fontes não identificadas de alto fluxo no HGPS. Além disso, o trabalho sugere um link evolutivo natural entre PWNe de meia-idade e o surgimento de halos estendidos em TeV, observando que sistemas mais antigos detectados por instrumentos futuros provavelmente corresponderão a morfologias semelhantes a halos em vez de nebulosas compactas clássicas. Este estudo serve como a primeira parte de uma série, estabelecendo a base estatística para a visibilidade, com um artigo subsequente planejado para focar nos resultados de "super eficiência" da reverberação.