The age sequence of young clusters in Perseus: Estimating ages from mass distributions

Resumen Técnico: La Secuencia de Edad de los Cúmulos Jóvenes en Perseo

Planteamiento del Problema
Establecer edades robustas para cúmulos estelares jóvenes (menores de 10 Myr) es un desafío crítico en los estudios de formación estelar. Los métodos de datación tradicionales, como la gyrocronología, el ajuste de isócronas y las mediciones del límite de agotamiento de Litio, suelen fallar o volverse poco fiables para poblaciones de la pre-secuencia principal. Una complicación significativa al determinar la Función de Masa Inicial (IMF) de cúmulos jóvenes es que la relación masa-luminosidad es una función fuerte de la edad; por lo tanto, suele ser necesario asumir una edad para derivar una distribución de masas. Este artículo aborda el problema inverso: en lugar de asumir una edad para determinar la función de masa, los autores proponen determinar edades relativas asumiendo una función de masa universal en una región de formación estelar.

Metodología
El estudio se centra en el complejo de formación estelar de Perseo, analizando seis cúmulos jóvenes: NGC 1333, Autochthe, IC 348, Heleus, Mestor, Electryon, Cynurus y Alcaeus. La metodología procede en tres etapas principales:

1. Selección Homogénea de Miembros: Utilizando datos de Gaia DR3, los autores construyen nuevas muestras homogéneas de miembros para seis agrupaciones de cúmulos (tratando a NGC 1333 y Autochthe como uno solo, y a Electryon y Cynurus como otro). Los criterios de selección incluyen posición espacial, movimiento propio, paralaje y un corte de edad (<10 Myr) basado en la comparación de isócronas. Las muestras se cruzan con fotometría 2MASS. Las posibles fuentes binarias se identifican y excluyen utilizando el estadístico RUWE (RUWE > 1.4) para probar los efectos de contaminación.
2. Derivación de la Distribución de Masas: Para cada cúmulo, las masas se estiman comparando la fotometría desreddenizada de 2MASS (bandas $J$ y $K_s$) contra isócronas teóricas de Baraffe et al. (2015). La extinción se calcula basándose en el vector de enrojecimiento relativo a la isócrona. Para propagar las incertidumbres fotométricas, una rutina de Monte Carlo genera 1,000 distribuciones de masa por cúmulo para un rango de edades asumidas (1, 2, 3, 4, 5 y 8 Myr).
3. Estimación de Edad mediante Comparación Estadística: El núcleo del método consiste en comparar las distribuciones de masa de todos los pares de cúmulos a través de todas las combinaciones de edad utilizando la prueba de Kolmogorov-Smirnov (KS). El supuesto subyacente es que los cúmulos estelares dentro del mismo complejo de Perseo se forman con una función de masa universal. Por lo tanto, la combinación de edad correcta para cualquier par de cúmulos debería resultar en distribuciones de masa estadísticamente indistinguibles (es decir, un valor p de la prueba KS $\ge$ 0.05). Los autores iteran a través de los 15 pares únicos de cúmulos y las 36 combinaciones de edad para encontrar el conjunto de edades que maximiza la consistencia de las distribuciones de masa en toda la región.

Resultados Clave
Al identificar las combinaciones de edad que producen distribuciones de masa estadísticamente indistinguibles, los autores establecen una secuencia de edad relativa clara para los cúmulos de Perseo:

• NGC 1333 + Autochthe: 1 Myr (el más joven)
• IC 348: 2 Myr
• Heleus: 2–3 Myr
• Mestor: 3–4 Myr
• Electryon + Cynurus: 4–5 Myr
• Alcaeus: 5–8 Myr (el más viejo)

El estudio encuentra que, si bien las edades absolutas para los cúmulos intermedios (Heleus, Mestor, Electryon) tienen rangos superpuestos, la secuencia relativa es robusta. Los resultados son insensibles a las variaciones de distancia dentro de los márgenes de error y permanecen inalterados cuando se excluyen las fuentes binarias.

Validación y Consistencia
Los autores validan su metodología utilizando un conjunto de cuatro cúmulos artificiales con edades de entrada conocidas y una función de masa uniforme. El método recuperó con éxito la secuencia de edad correcta (C, D, B, A) y los rangos de edad apropiados, confirmando la validez del enfoque sin necesidad de ajustes finos.

Además, la secuencia de edad derivada se alinea con:

• Fracciones de Disco: La secuencia coincide con la tendencia de disminución de las fracciones de disco circunestelar con el aumento de la edad, según lo reportado por Pavlidou (2022).
• Estimaciones de la Literatura: La secuencia es ampliamente consistente con estudios previos de ajuste de isócronas (por ejemplo, Kounkel et al. 2022; Olivares et al. 2023), aunque las edades absolas derivadas aquí tienden a ser ligeramente menores para algunos cúmulos, una diferencia que los autores atribuyen a los criterios de selección de muestra y a los cortes de magnitud en trabajos previos.

Significado
El artículo afirma introducir y validar un nuevo enfoque complementario para estimar las edades relativas de poblaciones estelares jóvenes. Al aprovechar la universalidad de la función de masa local, este método evita la dependencia circular de asumir una edad para derivar una función de masa. Los autores sugieren que esta técnica sirve como una verificación útil para los métodos establecidos de ajuste de isócronas, particularmente para establecer secuencias evolutivas relativas en regiones de formación estelar donde la datación absoluta sigue siendo incierta. El estudio no propone nuevas aplicaciones más allá de este marco comparativo, sino que enfatiza la robustez de la secuencia de edad derivada para el complejo de Perseo.