White dwarf planets in star clusters: gravitational scattering versus mass-loss effects

Resumen Técnico: Planetas de enanas blancas en cúmulos estelares: dispersión gravitacional frente a efectos de pérdida de masa

Planteamiento del problema
Las enanas blancas (WD, por sus siglas en inglés) sirven como laboratorios únicos para comprender la supervivencia y evolución a largo plazo de los sistemas planetarios. Si bien se sabe que la pérdida de masa post-secuencia principal altera las propiedades orbitales planetarias y agita los discos de escombros (provocando la contaminación atmosférica), la influencia del entorno de nacimiento estelar en estos sistemas permanece poco definida. Específicamente, no está claro cómo las interacciones dinámicas dentro de regiones densas de formación estelar (como la dispersión gravitacional y los sobrevuelos estelares) compiten con o modifican la evolución orbital impulsada por la pérdida de masa estelar. Estudios previos han modelado interacciones dinámicas en cúmulos o la evolución estelar de forma aislada, pero pocos han acoplado de manera consistente la evolución dinámica N-body con la evolución estelar desde la formación de la región de formación estelar hasta que las estrellas progenitoras se convierten en enanas blancas.

Metodología
Los autores realizaron simulaciones N-body de regiones de formación estelar utilizando el integrador kira dentro del entorno Starlab. Las simulaciones evolucionaron durante 1 Gyr, un tiempo escala suficiente para que las estrellas más masivas ($>2 M_\odot$) abandonen la secuencia principal y se conviertan en enanas blancas.

• Condiciones iniciales: Las simulaciones utilizaron $N_\star = 1000$ estrellas extraídas de una Función Inicial de Masa (IMF) de Maschberger (2013). A cada estrella con una masa inicial $M_\star \leq 2.5 M_\odot$ se le asignó un único planeta de masa de Júpiter. Los semiejes mayores iniciales ($a_p$) se muestrearon a partir de una distribución uniforme (0.1–50 au) o funciones delta en valores específicos (0.1, 1, 5 y 30 au). Las excentricidades e inclinaciones iniciales se fijaron en cero.
• Dinámica del cúmulo: Se probaron dos tipos de configuraciones iniciales de cúmulos:
  1. Distribuciones fractales: Generadas mediante un método de box-fractal con una dimensión fractal $D=1.6$ para representar regiones de formación estelar altamente subestructuradas y ligadas. Los radios iniciales se fijaron en 1 pc (alta densidad, $\sim 10^4 M_\odot \text{pc}^{-3}$) y 5 pc (menor densidad, $\sim 10^2 M_\odot \text{pc}^{-3}$).
  2. Esferas de Plummer: Distribuciones suaves y concentradas centralmente con un radio de semimasas de 0.1 pc.
• Evolución estelar: La evolución estelar y binaria se implementó utilizando tablas de consulta SeBa, sincronizadas con el paso de tiempo (timestep) de la simulación N-body. El código rastreó la pérdida de masa durante las fases de la rama gigante. Los planetas se consideraron engullidos si el radio estelar excedía el semieje mayor del planeta. Las simulaciones no incluyeron binarias estelares primordiales ni campos de marea galácticos externos.
• Conjuntos de simulación: Se ejecutaron diez realizaciones para cada conjunto de condiciones iniciales para evaluar la estocasticidad. El estudio se centró en el destino de los planetas que orbitan progenitores que evolucionan hacia enanas blancas (masas iniciales de 1–2.5 $M_\odot$).

Resultos clave

1. Dominio de la pérdida de masa: El principal motor de la evolución orbital para los planetas supervivientes es la pérdida de masa estelar, no la dispersión dinámica. Aunque los entornos densos destruyen o alteran significativamente las órbitas de hasta el 50% de los sistemas planetarios mediante la dispersión, el efecto de la pérdida de masa domina la dinámica para los sistemas restantes. Este comportamiento es en gran medida independiente de la densidad estelar inicial de la región de nacimiento.
2. Expansión del semieje mayor: Para los planetas que sobreviven al engullimiento y permanecen ligados, el semieje mayor aumenta significativamente debido a la pérdida de masa. En el límite adiabático, el semieje mayor final escala inversamente con la masa estelar final ($a_f \propto M_{inicial}/M_{final}$). Dado que los progenitores pierden $\sim$50–80% de su masa, las órbitas planetarias se expanden por factores de 2 a 5. Esta expansión crea una característica distintiva en la distribución del semieje mayor final (alcanzando un pico alrededor de 40–200 au para distribuciones iniciales de 0.1–50 au) que no está presente en la distribución de todos los planetas (incluyendo aquellos alrededor de estrellas de la secuencia principal) a los 100 Myr.
3. Cambios en la excentricidad: Aproximadamente el 40% de los planetas de enanas blancas experimentan un aumento en la excentricidad, en comparación con el 30% de todos los planetas. Crucialmente, estos cambios en la excentricidad son impulsados casi totalmente por encuentros dinámicos que ocurren antes de la fase de enana blanca. Una vez que la estrella se convierte en enana blanca, la pérdida de masa por sí sola induce cambios de excentricidad insignificantes para los planetas dentro de $\sim$100 au.
4. Captura y sistemas triples:
   • Planetas capturados: Aproximadamente el 10% de los sistemas planetarios de enanas blancas finales consisten en planetas capturados (planetas de masa libre capturados por una estrella). Estos sistemas típicamente exhiben semiejes mayores grandes ($>100$ au) y altas excentricidades ($e > 0.5$), asemejándose a una distribución térmica.
   • Sistemas triples: Las simulaciones produjeron sistemas triples de larga duración que contienen una enana blanca, un planeta y un tercer cuerpo (ya sea una estrella o un planeta capturado). Estas configuraciones son análogos potenciales a sistemas como PSR B1620-26 (AB) b.
5. Tasas de destrucción: En regiones densas, la tasa de destrucción de los sistemas planetarios es alta. Para los planetas que orbitan progenitores de enanas blancas, la fracción de planetas ligados supervivientes cae a $\sim$52% después de 1 Gyr, en comparación con el $\sim$71% para todas las estrellas. Esto se atribuye a las mayores tasas de interacción experimentadas por las estrellas más masivas.

Significado y conclusiones
El artículo concluye que el entorno dinámico del cúmulo de nacimiento tiene un impacto modesto en las propiedades orbitales de la mayoría de los planetas de enanas blancas supervivientes. Específicamente:

• Para $\sim$80% de los sistemas supervivientes, la configuración orbital actual está determinada por la pérdida de masa post-secuencia principal más que por la dispersión del cúmulo de nacimiento.
• Para el $\sim$20% restante, los encuentros dinámicos alteran significativamente el semieje mayor.
• La presencia de planetas de separación amplia y alta excentricidad alrededor de enanas blancas (por ejemplo, WD 0806-661 b) puede ser indicativa de un origen por captura dentro de un cúmulo, un proceso más eficiente en cúmulos que en el campo galáctico.
• Por el contrario, los planetas muy cercanos (por ejemplo, WD 1856+534 b) probablemente llegaron a sus órbitas actuales a través de mecanismos de dispersión (potencialmente involucrando otros planetas o estrellas) que ocurrieron temprano en la vida del cúmulo, seguidos de la truncación por marea, en lugar de la expansión por pérdida de masa por sí sola.

Los autores proporcionan una síntesis de población de planetas gigantes (1–100 au) alrededor de enanas blancas que emergen de cúmulos de nacimiento, ofreciendo restricciones para futuras observaciones de instalaciones como el Telescopio Espacial Nancy Grace Roman, Gaia y JWST. El estudio enfatiza que, si bien los entornos de nacimiento dan forma a la población a través de la destrucción y la captura, la expansión orbital fundamental de los sistemas supervivientes es una consecuencia robusta de la evolución estelar.