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

Enunciado del Problema
Las nebulosas de viento de púlsar (PWNe) representan la población más grande de fuentes galácticas de rayos gamma de muy alta energía (VHE; $E > 100$ GeV) y sirven como laboratorios críticos para estudiar la aceleración de partículas y las interacciones entre púlsares y restos de supernova (SNR). Si bien instrumentos actuales como H.E.S.S., MAGIC y VERITAS han identificado firmemente numerosas PWNe, las predicciones realistas a nivel de población han carecido históricamente de un tratamiento detallado de la "fase de reverberación". Esta fase ocurre cuando la onda de choque inversa del SNR viaja de regreso hacia el centro de la explosión, colisionando con y comprimiendo la PWN. Este evento altera significativamente la dinámica de la nebulosa, amplifica los campos magnéticos, reenergiza las poblaciones de partículas y aumenta drásticamente las pérdidas radiativas. Los enfoques de modelado anteriores enfrentaron una dicotomía: los modelos simples de capa delgada de una sola zona son computacionalmente eficientes pero no logran capturar la compleja retroalimentación no lineal de la fase de reverberación, mientras que las simulaciones hidrodinámicas (HD) y magnetohidrodinámicas (MHD) multidimensionales son precisas pero prohibitivamente costosas para estudios estadísticos que involucran miles de fuentes durante $10^5$ años. En consecuencia, las predicciones para la población galáctica de PWNe en TeV, particularmente para los observatorios futuros, corren el riesgo de estar sesgadas si esta fase se descuida o se simplifica en exceso.

Metodología
Para abordar estas limitaciones, los autores emplean el marco híbrido TIDE+L, que combina la eficiencia computacional de un modelo dinámico de capa delgada (TIDE) para la fase de expansión libre con un tratamiento lagrangiano de la estructura del SNR durante la reverberación. Este enfoque permite la evolución autoconsistente del sistema PWN–SNR a través de todas las etapas, incluidos los complejos episodios de compresión de la reverberación, mientras se rastrea el espectro de partículas sujeto a pérdidas radiativas y adiabáticas.

El estudio construye una población sintética de aproximadamente 2.400 PWNe galácticas, de las cuales 1.600 fuentes se muestrean aleatoriamente en 1.000 realizaciones para garantizar la robustez estadística. La población se muestrea a partir de distribuciones de:

• Propiedades del púlsar: Periodo de rotación inicial (distribución gaussiana, truncada a 10 ms), campo magnético inicial (log-normal) e índice de frenado ($n=2.33$).
• Propiedades del SNR y del entorno: Energía de la explosión (gaussiana truncada), masa de eyección (derivada de la masa de la secuencia principal de edad cero mediante modelos de evolución estelar GENEVA) y densidad del medio interestelar (ISM) (distribución uniforme).
• Parámetros espectrales: Espectros de inyección de partículas, factores de contención y fracciones magnéticas.

Cada sistema se evoluciona hasta una edad máxima de $10^5$ años. Los flujos de rayos gamma resultantes se comparan luego con los umbrales de sensibilidad y la cobertura del cielo de instrumentos actuales y futuros: H.E.S.S., HAWC, el Observatorio del Telescopio Cherenkov (CTAO), LHAASO y el Observatorio Gamma de Gran Campo del Sur (SWGO). El análisis tiene en cuenta los efectos de extensión de la fuente, que degradan la sensibilidad para fuentes extendidas.

Contribuciones Clave

1. Modelado Realista de la Reverberación: El artículo introduce una síntesis de población estadísticamente robusta que incorpora explícitamente el tratamiento lagrangiano de la reverberación, superando las prescripciones analíticas simplificadas de los modelos anteriores de capa delgada.
2. Aplicación del Marco Híbrido: Demuestra la aplicación del marco TIDE+L a un estudio de población a gran escala, cerrando la brecha entre la física hidrodinámica detallada y los requisitos computacionales de la síntesis de poblaciones.
3. Perspectivas de Detección: El estudio proporciona las primeras estimaciones detalladas de la población detectable de PWNe galácticas para instrumentos de próxima generación, cuantificando específicamente el rendimiento de CTAO, SWGO y LHAASO.
4. Etapas Evolutivas: Categoriza la población en fases de expansión libre, compresión y post-compresión, analizando la detectabilidad y las características espectrales de cada una.

Resultados

• Detectabilidad: El modelo predice que CTAO detectará un orden de magnitud más de PWNe que las detectadas firmemente actualmente en el rango de TeV. Específicamente, se espera que CTAO detecte aproximadamente 172 fuentes (10.7% de la población), aumentando a 226 (14.1%) con una configuración mejorada "CTAO+". Esto es significativamente mayor que las ~52 fuentes predichas para H.E.S.S.
• Impacto de la Reverberación: Las comparaciones entre el marco TIDE+L y el modelo TIDE simplificado revelan que descuidar la dinámica realista de la reverberación conduce a una subestimación de las fuentes detectables. El modelo simplificado tiende a sobreestimar la compresión y la amplificación del campo magnético, lo que resulta en pérdidas de sincrotrón excesivas y una población de partículas de alta energía agotada, reduciendo así el flujo de rayos gamma predicho.
• Distribución Evolutiva: La mayoría de la población sintética (87.5%) se encuentra en la fase post-compresión. Sin embargo, las fuentes en esta fase suelen ser más antiguas, más grandes e intrínsecamente más tenues, lo que resulta en fracciones de detectabilidad más bajas en comparación con las de las fases de expansión libre o compresión.
• Evolución Espectral: Estudios de casos individuales ilustran que el factor de compresión (CF) gobierna la intensidad y duración de la amplificación de la luminosidad y el endurecimiento/ablandamiento espectral durante la reverberación. Las fuentes de alto CF exhiben cambios rápidos y agudos en la distribución de energía espectral (SED), mientras que las fuentes de bajo CF muestran una evolución más suave y retardada.
• Correlaciones: El estudio encuentra una correlación entre la luminosidad integrada de 1–10 TeV y la potencia de frenado del púlsar ($L_{1-10 \text{ TeV}} \propto \dot{E}^{1.15}$). Los autores señalan que recuperar la pendiente más suave observada en los datos de HGPS ($L \propto \dot{E}^{0.59}$) requiere aplicar cortes de selección específicos (por ejemplo, excluir fuentes con baja potencia de frenado o edades mayores), lo que sugiere que los sesgos de las encuestas y las limitaciones instrumentales influyen en las correlaciones observadas.
• Extensión de la Fuente: Las extensiones angulares simuladas son generalmente consistentes con los datos de HGPS, aunque el modelo predice una mayor fracción de fuentes con radios menores que el mínimo observado en HGPS. Esto se atribuye a la suposición de simetría esférica y al descuido de la mezcla y el movimiento propio del púlsar.

Importancia
El artículo afirma que el modelado realista de la fase de reverberación es esencial para predecir con precisión la población galáctica de PWNe en TeV. Al incorporar la dinámica compleja de las interacciones de la onda de choque inversa del SNR, el estudio proporciona una base teórica más fiable para interpretar los catálogos actuales de fuentes y anticipar los rendimientos de las próximas encuestas. Los resultados confirman que CTAO jugará un papel dominante en el próximo censo de poblaciones en TeV, resolviendo potencialmente la naturaleza de las fuentes no identificadas de alto flujo en HGPS. Además, el trabajo sugiere un vínculo evolutivo natural entre las PWNe de mediana edad y la aparición de halos extendidos en TeV, señalando que los sistemas más antiguos detectados por instrumentos futuros probablemente corresponderán a morfologías tipo halo en lugar de nebulosas compactas clásicas. Este estudio sirve como la primera parte de una serie, estableciendo la base estadística para la visibilidad, con un artículo posterior planificado para centrarse en los resultados de "supereficiencia" de la reverberación.