Dilution of accreted planetary matter in hot DA white dwarfs according to their mass

Planteamiento del Problema
Las observaciones de enanas blancas (WD, por sus siglas en inglés) revelan frecuentemente la presencia de elementos pesados contaminantes en sus atmósferas, lo cual se atribuye a la acreción de escombros planetarios. Las abundancias observadas representan un equilibrio entre el flujo de acreción y la dilución interna de estos elementos. Estudios estadísticos recientes (por ejemplo, Ould Rouis et al. 2024, en adelante OR24) han identificado una tendencia en la que las enanas blancas DA (ricas en hidrógeno) masivas exhiben tasas de contaminación significativamente menores en comparación con las enanas blancas DA de menor masa. Esta tendencia fue interpretada como evidencia de una menor ocurrencia de sistemas planetarios alrededor de las progenitoras de la secuencia principal masivas de estas estrellas. El problema central abordado es si los procesos de transporte interno, específicamente la interacción entre la difusión atómica y la convección termohalina, varían con la masa de la enana blanca de una manera que pueda explicar las tendencias de contaminación observadas sin invocar diferencias en la formación de sistemas planetarios.

Metodología
Los autores computaron la eficiencia de la difusión atómica y la convección termohalina utilizando modelos estáticos de enanas blancas DA con masas variables (0.6, 0.8 y 1.0 $M_\odot$), temperaturas efectivas (20,000, 25,000 y 30,000 K) y masas de la envoltura de hidrógeno. Estos modelos fueron generados utilizando el código STELUM, asumiendo composiciones del núcleo derivadas de cálculos evolutivos y un amortiguador mixto de C/O/He.

Para simular la acreción, los autores asumieron una masa total acrecionada de $10^{16}$ g (aproximando la masa de Fobos) inicialmente dispersa en una cáscara esférica. Probaron dos escenarios para asegurar la robustez:

1. Caso 1: Un perfil de decaimiento exponencial fijo para la materia acrecionada en relación con el radio estelar para todos los modelos.
2. Caso 2: Un perfil de decaimiento ajustado para asegurar gradientes de peso molecular medio ($\mu$) idénticos en todos los modelos, aislando el efecto de la densidad estelar.

Los coeficientes de difusión atómica se estimaron utilizando aproximaciones para elementos traza en un medio de hidrógeno completamente ionizado (Michaud et al. 2015). La convección termohalina se modeló utilizando la prescripción de Kippenhahn et al. (1980), que relaciona el coeficiente de mezcla con la difusividad térmica y la relación de densidad ($R_0$). Los autores verificaron que sus computaciones permanecieran dentro del rango de estabilidad donde esta aproximación es válida.

Resultados Clave

• Dominio de la Convección Termohalina: El estudio confirma que la convección termohalina es consistentemente más eficiente que la difusión atómica para diluir los elementos pesados acrecionados en enanas blancas DA calientes.
• Eficiencia Dependiente de la Masa: Contrario a la hipótesis de que el mezclado interno podría ser más eficiente en estrellas masivas (ocultando así la contaminación), los autores encontraron que la convección termohalina es menos eficiente en enanas blancas masivas. Esto se atribuye a las mayores densidades internas de las enanas blancas masivas (una consecuencia de la relación masa-radio), las cuales reducen la difusividad térmica ($\kappa_T \propto 1/\rho^2$) y, consecuentemente, disminuyen la eficiencia del mezclado termohalino.
• Dependencia de la Temperatura: La eficiencia del transporte termohalino aumenta con la temperatura efectiva debido a la dependencia $T^3$ de la difusividad térmica.
• Implicaciones para las Firmas de Contaminación: Debido a que la dilución es menos eficiente en enanas blancas masivas, una tasa de acreción dada debería resultar en firmas de contaminación superficial mayores en enanas blancas masivas en comparación con las de baja masa, asumiendo que la difusión atómica y el mezclado termohalino son los únicos procesos relevantes.

Significado y Conclusiones
El artículo concluye que la tendencia observada de menor contaminación en enanas blancas DA masivas no puede explicarse únicamente mediante los procesos de dilución interna (difusión atómica y convección termohalina). De hecho, la física de estos procesos de transporte sugiere que la contaminación debería ser más fácilmente detectable en enanas blancas masivas.

Por lo tanto, la tendencia observada por OR24 —donde las enanas blancas masivas parecen estar menos contaminadas— probablemente refleja una diferencia genuina en la ocurrencia o naturaleza de los sistemas planetarios alrededor de sus progenitoras de la secuencia principal, en lugar de ser un artefacto de la física estelar interna. Sin embargo, los autores mantienen una postura modesta, señalando que deben considerarse otros factores antes de extraer conclusiones definitivas:

• Mezclado Rotacional: Las enanas blancas masivas pueden rotar más rápido, lo que potencialmente induce un mezclado rotacional que podría contrarrestar la reducida eficiencia termohalina.
• Incertidumbres de la Progenitora: Vincular la masa de la enana blanca con la masa de su progenitora en la secuencia principal es complicado por las incertidumbres en la pérdida de masa y la posibilidad de fusiones binarias, las cuales constituyen una fracción significativa de las enanas blancas masivas.
• Limitaciones del Modelo: El presente estudio utiliza modelos estáticos y condiciones iniciales simplificadas para la acreción. Se requiere un modelado dependiente del tiempo de la acreción y la dilución para obtener conclusiones más sólidas.

En última instancia, el trabajo enfatiza que omitir la convección termohalina conduce a una subestimación de la dilución y de las tasas de acreción derivadas en enanas blancas ricas en H, pero incluso al contabilizar esto, la tendencia de la contaminación dependiente de la masa apunta hacia diferencias intrínsecas en la formación de sistemas planetarios.