Drifters on the edge of town: $λ$ Boötis stars in clusters

Resumen Técnico: Errantes en los bordes de la ciudad: estrellas $\lambda$ Boötis en cúmulos

Planteamiento del Problema
Las estrellas $\lambda$ Boötis ($\lambda$ Boo) son estrellas de tipo A con una composición química peculiar, caracterizadas por abundancias similares a las solares de elementos ligeros (C, N, O, S), pero con deficiencias significativas de elementos del grupo del hierro. La explicación teórica predominante postula que estas estrellas acrecen gas prístino y pobre en metales del Medio Interestelar (ISM). Sin embargo, este escenario enfrenta un desafío significativo: descubrimientos recientes de estrellas $\lambda$ Boo dentro o en las periferias de cúmulos estelares contradicen la expectativa de que la intensa radiación ionizante de estrellas intermedias y masivas (>5 M$_\odot$) en tales entornos impediría la acreción de gas prístino del ISM o eliminaría los discos circunestelares necesarios para el proceso. Además, existe una tensión entre las tasas de acreción requeridas para producir las abundancias de $\lambda$ Boo ($10^{-14}$ a $10^{-9}$ M$_\odot$ año$^{-1}$) y las tasas mucho menores inferidas de observaciones de enanas blancas (<$10^{-17}$ M$_\odot$ año$^{-1}$). El problema específico abordado es si las historias dinámicas de las estrellas tipo A en cúmulos en evolución les permiten escapar del núcleo del cúmulo de alta radiación, acrecer gas prístino y, posteriormente, regresar, reconciliando así la presencia de estrellas $\lambda$ Boo en entornos de cúmulos con la hipótesis de acreción del ISM.

Metodología
Los autores emplean simulaciones de N-cuerpos puramente gravitacionales para modelar la formación y evolución dinámica de regiones de formación estelar durante 100 Myr.

• Configuración de la Simulación: Se realizaron cuatro conjuntos de simulaciones con números totales de estrellas ($N_*$) de 150, 500, 1000 y 3000. Las condiciones iniciales incluyeron subestructura espacial y cinemática generada mediante el método fractal de caja (dimensión fractal $D=1.6$) y una relación virial de $\alpha_{vir} = 0.3$, lo que condujo al colapso inicial y a la formación de cúmulos concentrados centralmente. Las masas estelares se muestrearon de una Función de Masa Inicial (IMF) de Maschberger (2013) con un límite superior de 50 M$_\odot$. Se excluyeron las binarias primordiales.
• Análisis de Post-procesamiento: Las simulaciones fueron analizadas ex post facto para determinar:
  1. Potencial de Acreción: Se rastrearon las estrellas tipo A para identificar cuándo se movían más allá del radio de Jacobi (tidal) ($r_J$), definido como el límite donde podrían encontrar gas prístino del ISM. Cuando estaban fuera de $r_J$, se calcularon las tasas de acreción de Bondi-Hoyle-Lyttleton (BHL) basadas en la velocidad de la estrella respecto a una nube de ISM estacionaria ($\rho = 10$ cm$^{-3}$).
  2. Entorno de Radiación: Se calcularon los flujos de UV lejano (FUV) y EUV incidentes en cada estrella tipo A basándose en la distancia a las estrellas masivas (>5 M$_\odot$ y >20 M$_\odot$, respectivamente).
  3. Supervivencia del Disco: Utilizando la rejilla de modelos de fotoevaporación de FRIED y prescripciones analíticas, los autores modelaron la pérdida de masa de los discos circunestelares de las estrellas tipo A. Un disco se consideró "destruido" si la masa total perdida superaba la masa inicial del disco (asumida como el 10% de la masa estelar). Una estrella fue considerada candidata para la formación de $\lambda$ Boo solo si acreía material del ISM mientras retenía su disco.
  4. Restricciones Temporales: Los cálculos se repitieron con una restricción que impedía la acreción después de 40 Myr para considerar el potencial enriquecimiento del ISM por supernovas de Tipo Ia.

Resultados Clave

• Historias Dinámicas: En todas las realizaciones de las simulaciones, la mayoría (90–98%) de las estrellas tipo A viajaron más allá del radio de Jacobi en algún momento durante la evolución de 100 Myr. Esto ocurrió tanto por expulsión del núcleo del cúmulo como por "deriva", donde las estrellas se movieron fuera del radio de marea y posteriormente regresaron a las regiones internas.
• Tasas de Acreción: Las estrellas tipo A que se movieron más allá del radio de Jacobi alcanzaron tasas de acreción BHL que oscilaron entre $10^{-12}$ y $10^{-5}$ M$_\odot$ año$^{-1}$. Las estrellas expulsadas a altas velocidades tuvieron tasas menores ($10^{-12}–10^{-10}$ M$_\odot$ año$^{-1}$), mientras que aquellas en órbitas centradas en el cúmulo que derivaron hacia afuera y regresaron tuvieron tasas más altas ($10^{-9}–10^{-5}$ M$_\odot$ año$^{-1}$), solapándose con el rango requerido para explicar las abundancias de $\lambda$ Boo.
• Supervivencia del Disco vs. Radiación: La supervivencia de los discos circunestelares dependió fuertemente del número total de estrellas (y, por ende, del número de fuentes ionizantes) en el cúmulo.
  • En regiones de baja masa ($N_* = 150$), aproximadamente el 55% de las estrellas tipo A que se aventuraron más allá del radio de Jacobi retuvieron sus discos por >10 Myr.
  • En regiones de alta masa ($N_* = 3000$), el 0% de las estrellas tipo A retuvo sus discos debido a los intensos campos de FUV/EUV, lo que imposibilita efectivamente la formación de $\lambda$ Boo en estos entornos.
  • En regiones intermedias ($N_* = 500, 1000$), una fracción de las estrellas tipo A (8.8% y 4.8% respectivamente) retuvo sus discos mientras acreaba.
• Sistemas Binarios: El estudio encontró que los sistemas binarios anchos (>100 au) que contienen estrellas tipo A pueden formarse por captura durante la disolución de las regiones de formación estelar. Crucialmente, una estrella podría acrecer material prístino del ISM, desarrollar abundancias de $\lambda$ Boo y luego capturar una compañera (o viceversa). Este mecanismo permite la existencia de sistemas binarios donde los componentes tienen abundancias químicas significativamente diferentes. La distribución de separación simulada de tales binarias intersecta con las binarias anchas observadas que contienen estrellas $\lambda$ Boo (por ejemplo, HD 87304 + CD-33 6615B).

Significancia y Reivindicaciones
El artículo afirma demostrar que la presencia de estrellas $\lambda$ Boo en o cerca de cúmulos estelares no invalida inherentemente el escenario de acreción del ISM. Los autores argumentan que:

1. La Escapa Dinámica es Suficiente: Las estrellas tipo A no necesitan ser expulsadas permanentemente para acrecer gas prístino; pueden derivar más allá del radio de marea, acrecer y regresar al centro del cúmulo.
2. Umbrales de Radiación: La viabilidad del escenario de acreción está fuertemente correlacionada con la masa del cúmulo. Los cúmulos de baja masa (similares a Theia 139 y HSC 1640, donde se han encontrado recientemente estrellas $\lambda$ Boo) poseen campos de radiación lo suficientemente débiles como para permitir la supervivencia del disco, mientras que los cúmulos masivos no lo hacen.
3. Diagnóstico Binario: La observación de sistemas binarios con abundancias dispares entre sus componentes no invalida la hipótesis de acreción. Tales sistemas pueden formarse dinámicamente después de que el evento de acreción haya ocurrido.

Los autores mantienen la modestia respecto a la validez absoluta de las tasas de acreción, señalando que la tensión entre las tasas requeridas y los límites de las enanas blancas permanece sin resolver. Concluyen que, si bien sus modelos dinámicos respaldan la posibilidad de que las estrellas $\lambda$ Boo se formen en cúmulos, se requiere más investigación para validar las tasas de acreción específicas y la interacción física entre los discos que acrecen y los campos de radiación.