Dust attenuation in galaxies at cosmic dawn from the FirstLight simulations

Título: Atenuación por polvo en galaxias en el amanecer cósmico a partir de las simulaciones FirstLight

1. El Problema

El estudio de la atenuación del polvo en galaxias del "amanecer cósmico" (redshifts $z \approx 6 - 8$) es un desafío debido a la falta de observaciones detalladas. Tradicionalmente, se asume que las curvas de atenuación derivadas de galaxias locales o de bajo redshift (como la ley empírica de Calzetti et al. 2000) son aplicables a estas épocas tempranas. Sin embargo, las condiciones del medio interestelar (ISM) en el universo temprano difieren significativamente:

• Baja metalicidad: El contenido de metales (y por ende de polvo) es menor.
• Geometrías complejas: La distribución del polvo y las estrellas puede ser muy diferente a la de galaxias maduras.
• Efectos del CMB: A altos redshifts, la Radiación de Fondo de Microondas (CMB) es más intensa y puede calentar el polvo, afectando la emisión infrarroja.

Existe incertidumbre sobre si las relaciones estándar entre la pendiente del espectro UV ($\beta_{UV}$) y el exceso infrarrojo (IRX) se mantienen en estas condiciones extremas.

2. Metodología

Los autores combinan simulaciones cosmológicas de alta resolución con códigos de transferencia radiativa (RT) para crear "observaciones sintéticas":

• Simulaciones FirstLight (FL): Se utilizan simulaciones de "zoom-in" cosmológicas que cubren un volumen completo de masas y volúmenes, resolviendo galaxias con resolución de parsecs. Se selecciona una muestra de 120 galaxias (60 a $z=6$ y 60 a $z=8$) con masas estelares que van desde $10^5$ hasta $10^{10} M_\odot$.
• Extracción de Datos:
  • Se extraen las distribuciones de partículas estelares y densidad de polvo.
  • La masa de polvo se calcula asumiendo una relación constante polvo-metal ($DMR = 0.4$) basada en la metalicidad del gas.
• Espectros Intrínsecos (BPASS): Se utiliza el modelo de síntesis espectral de poblaciones estelares binarias (BPASS) para generar los espectros UV/ópticos intrínsecos de las partículas estelares según su edad, masa y metalicidad. Para galaxias masivas con muchas partículas, se realiza un muestreo estadístico para reducir el costo computacional sin perder precisión.
• Simulación de Transferencia Radiativa (POLARIS):
  • Se utiliza el código POLARIS, basado en la técnica de Monte Carlo 3D.
  • Modelo de Polvo: Se asumen granos esféricos de grafito y silicato con una distribución de tamaños de Mathis et al. (1977).
  • Proceso: Se simula la propagación de fotones desde las estrellas a través del ISM, calculando la absorción, dispersión y reemisión térmica.
  • Se ejecutan dos tipos de simulaciones: una para la temperatura del polvo (equilibrio termodinámico local) y otra para el campo de radiación transmitido y emitido.

3. Contribuciones Clave

• Análisis sistemático de curvas de atenuación: Se calculan curvas de atenuación para una muestra diversa de galaxias en el universo temprano, vinculando explícitamente la forma de la curva con la masa estelar de la galaxia.
• Validación de modelos en altos $z$: Se pone a prueba la aplicabilidad de los modelos estándar (Calzetti, SMC) en condiciones de baja metalicidad y alta redshift.
• Estudio del efecto del CMB: Se cuantifica cómo el calentamiento adicional por el CMB en $z=8$ altera la relación IRX-$\beta_{UV}$ en comparación con $z=6$.
• Parametrización: Se propone un modelo de regresión múltiple que relaciona la atenuación en 1500 Å con la masa estelar y la magnitud absoluta.

4. Resultados Principales

• Dependencia de la Masa Estelar en las Curvas de Atenuación:

  • Galaxias Masivas ($M_{stars} \approx 10^9 M_\odot$): Sus curvas de atenuación siguen el modelo de Calzetti (pendientes más planas). Esto se debe a una mayor opacidad óptica y geometrías más complejas (agrupamiento).
  • Galaxias de Baja Masa ($M_{stars} \approx 10^7 - 10^8 M_\odot$): Presentan curvas mucho más empinadas, consistentes con el modelo de la Nube de Magallanes Pequeña (SMC). La baja metalicidad y la menor densidad de polvo contribuyen a este comportamiento.
  • La transición de curvas tipo SMC a tipo Calzetti ocurre a medida que la galaxia crece en masa y metalicidad.

• Pendientes UV y Ópticas:

  • Las pendientes en el UV y óptico se ajustan mejor al modelo Calzetti modificado.
  • Existe una ligera preferencia hacia modelos de ley de potencia para las galaxias con mayor atenuación y masa.
  • La pendiente UV ($S_{UV}$) y la óptica ($S_{opt}$) muestran una correlación clara con la masa estelar.

• Relación IRX - $\beta_{UV}$:

  • En $z=6$: La relación sigue la tendencia del modelo de Calzetti, pero con un desplazamiento hacia valores de $\beta_{UV}$ ligeramente más bajos (espectros más azules) debido a las bajas metalicidades.
  • En $z=8$: Se observa un desplazamiento hacia valores más altos de IRX (mayor emisión infrarroja) para un mismo $\beta_{UV}$. Esto se atribuye al calentamiento del polvo por el CMB, que es más intenso a mayor redshift, aumentando la luminosidad IR incluso si la atenuación UV es similar.
  • Las galaxias de baja masa tienden a tener valores de IRX y $\beta_{UV}$ bajos debido a su juventud y baja atenuación.

• Límite de Masa para la Atenuación:

  • Se identifica que las galaxias con masas estelares inferiores a $\approx 10^8 M_\odot$ son ópticamente delgadas y no absorben una cantidad significativa de luz UV ($A_{1500} < 0.1$).

5. Significado e Impacto

Este trabajo demuestra que no es válido aplicar una única ley de atenuación (como la de Calzetti) a todas las galaxias del amanecer cósmico. La elección del modelo de corrección de extinción debe depender estrictamente de la masa estelar y la metalicidad de la galaxia objetivo.

• Para observaciones futuras (JWST, ELT): Los resultados sugieren que las estimaciones de tasas de formación estelar basadas en luz UV en galaxias de baja masa podrían estar subestimadas si se asume una atenuación tipo Calzetti, ya que sus curvas son más empinadas (tipo SMC).
• Física del ISM: Resalta la importancia de la geometría polvo-estrella y los efectos del entorno cosmológico (CMB) en la interpretación de los datos observacionales de alta redshift.
• Validación de Simulaciones: Confirma que las simulaciones cosmológicas de alta resolución, cuando se acoplan con códigos de transferencia radiativa, pueden reproducir la diversidad observada y predecir desviaciones de los modelos empíricos locales.