Disk-Regulated Mass Transfer Between Rotating Non-Degenerate Stars: Insights from Be and sdOB Binaries

Resumen Técnico: Transferencia de Masa Regulada por Disco entre Estrellas No Degeneradas en Rotación

Planteamiento del Problema
La transferencia de masa entre estrellas no degeneradas es un proceso crítico pero aún poco comprendido en la evolución binaria. Una incertidumbre central es la eficiencia de acreción —la fracción de la masa transferida que es retenida por el acreedor—. Las prescripciones estándar de síntesis de poblaciones binarias (BPS) y las simulaciones detalladas de evolución estelar (por ejemplo, utilizando MESA) emplean típicamente una prescripción de acreción "limitada por rotación". En este marco, la acreción se restringe una vez que el acreedor se aproxima a su tasa de rotación crítica de ruptura, lo que a menudo conduce a eficiencias de acreción muy bajas.

Sin embargo, las restricciones observacionales de las binarias de estrellas Be, particularmente los sistemas Be + enana subdwarf O/B (sdOB), sugieren eficiencias de acreción significativamente más altas (a menudo $\gtrsim 30\%$ a $>50\%$). Los modelos existentes limitados por rotación subestiman sistemáticamente las masas de las estrellas Be en estos sistemas y no logran reproducir el espacio de parámetros observado de las binarias Be+sdOB. Además, estudios previos indican que ni las prescripciones limitadas por rotación ni las limitadas térmicamente pueden, por sí solas, satisfacer los requisitos de eficiencia derivados de las poblaciones observadas de Be+sdOB.

Metodología
Los autores proponen una prescripción de acreción de masa de "acoplamiento disco-estrella" motivada físicamente para abordar estas discrepancias. El mecanismo físico central postula que un disco de acreción o de decréción regula el momento angular transportado al acreedor, permitiendo que la acreción de masa continúe incluso cuando la estrella rota cerca de su límite crítico.

1. Regulación del Momento Angular: Basándose en las soluciones de disco delgado en estado estacionario de Popham & Narayan (1991), los autores modelan el momento angular específico acrecido por la estrella ($j_{acc}$) como una función de la tasa de rotación estelar ($\Omega/\Omega_{crit}$).

   • A tasas de rotación bajas, $j_{acc}$ es aproximadamente el valor Kepleriano en la superficie estelar.
   • A medida que la estrella se aproxima a una tasa de rotación umbral (adoptada como $\Omega/\Omega_{crit} \approx 0.9$), el esfuerzo de cizalladura del disco permite que el momento angular sea transportado hacia afuera de la estrella hacia el disco. Consecuentemente, $j_{acc}$ disminuye drásticamente, pudiendo volverse negativo, lo que permite una acreción de masa continua sin que la estrella exceda su rotación crítica.
   • Los autores implementan una forma funcional suavizada para $j_{acc}$ para asegurar la estabilidad numérica en las simulaciones 1D, evitando discontinuidades que surgirían de un umbral estricto.

2. Leyes de Conservación y Cálculo de Eficiencia: El modelo asume que el exceso de momento angular extraído de la estrella es transportado a través del disco y eliminado mediante la pérdida de masa en el borde exterior del disco ($R_{out}$). Al aplicar la conservación de la masa y del momento angular, la eficiencia de acreción ($\eta = \dot{M}_{acc}/\dot{M}_{tr}$) se calcula de forma autoconsistente. La eficiencia depende de la relación entre el momento angular específico en el borde del disco y el del radio de acreción.

3. Marco de Simulación: Los autores implementaron esta prescripción dentro del marco POSYDON, utilizando MESA (Modules for Experiments in Stellar Astrophysics) para evolucionar modelos binarios.

   • Deshabilitaron el esquema de viento implícito utilizado para prevenir la rotación supercrítica, dependiendo en su lugar de vientos explícitos aumentados por rotación y de la reducción en la tasa de acreción en el régimen supercrítico para restaurar la rotación subcrítica.
   • Se simuló una cuadrícula de modelos binarios con metalicidad solar, variando las masas primarias ($2-9 M_\odot$), las razones de masa ($0.4-0.95$) y los periodos orbitales.
   • El estudio también investigó el impacto de los tratamientos de sobreimpulso (overshooting) del núcleo (comparando la configuración estándar de POSYDON con una prescripción de sobreimpulso reducido) sobre las masas finales de los componentes.

Contribuciones Clave y Resultados

• Estimaciones Analíticas: Los autores derivaron expresiones analíticas para la eficiencia de acreción. Demuestran que, en el régimen de acreción por disco, la eficiencia depende de la rotación estelar, la razón de masa y la relación entre el radio estelar y el radio de la lóbulos de Roche. En el régimen de impacto directo, depende de la rotación y la razón de masa. Se encontró que las eficiencias de acreción ponderadas por masa global están cerca de los valores esperados cerca de la tasa de rotación umbral, donde el momento angular específico acrecido declina bruscamente.
• Comparación de la Evolución Binaria:
  • Modelo Limitado por Rotación: En este modelo estándar, los acreedores alcanzan la rotación crítica tras acrecer solo una pequeña fracción de la masa transferida ($\sim 0.2 M_\odot$ en un caso de prueba), lo que conduce a eficiencias bajas ($<10\%$ para la mayoría de los sistemas). Esto resulta en estrellas Be que son sistemáticamente menos masivas de lo observado.
  • Modelo de Acoplamiento Disco-Estrella: Este modelo permite que los acreedores acrezcan significativamente más masa (por ejemplo, $\sim 1.6 M_\odot$ en la fase rápida del caso de prueba) con eficiencias que alcanzan el $\sim 40-60\%$. El disco extrae el momento angular, permitiendo que la estrella se expanda y acrecione mientras permanece cerca del umbral de rotación crítica.
• Reproducción de Observaciones:
  • Al aplicarse a sistemas Be+sdOB, el modelo limitado por rotación falla en producir estrellas Be con masas consistentes con las observaciones, particularmente para sistemas de órbita amplia (Caso B).
  • El modelo de acoplamiento disco-estrella produce estrellas Be con masas que cubren una parte significativa del rango observado.
  • Los autores encontraron que las masas finales de los componentes son sensibles al sobreimpulso (overshooting) asumido. Un modelo que combina la prescripción de acoplamiento disco-estrella con un sobreimpulso reducido produce masas de componentes para sistemas Be+sdOB que están en mayor acuerdo con las observaciones que los modelos estándar de alto sobreimpulso.
• Periodos Orbitales: Aunque el modelo de acoplamiento disco-estrella aumenta la retención de masa, los periodos orbitales resultantes son comparables o ligeramente más largos que en el modelo limitado por rotación. Esto se debe a que la transferencia de masa más conservativa reduce la masa total perdida del sistema, lo que conduce a razones de masa más extremas que aceleran la expansión orbital, contrarrestando parcialmente la pérdida de momento angular desde el borde del disco.

Significancia y Reivindicaciones
El artículo afirma que la discrepancia entre las propiedades observadas de las binarias Be+sdOB y los resultados actuales de las simulaciones puede resolverse adoptando un modelo de transferencia de masa regulado por disco. Los autores argumentan que:

1. Las prescripciones de acreción limitadas por rotación subestiman sistemáticamente las masas de las estrellas Be en sistemas Be+sdOB.
2. Un mecanismo de acoplamiento disco-estrella, que permite la regulación del transporte de momento angular, permite eficiencias de acreción más altas consistentes con las restricciones observacionales (específicamente el requisito de eficiencias $>50\%$ en muchos sistemas).
3. La combinación de esta nueva prescripción de acreción con un sobreimpulso del núcleo reducido proporciona una explicación más físicamente consistente para las masas y periodos orbitales observados de las binarias Be+sdOB.

Los autores reconocen que la eficiencia real depende de factores inciertos como la estructura del disco, los campos magnéticos y las condiciones específicas de formación de disco y pérdida de masa (por ejemplo, vía puntos L2 o puntos calientes). Sin embargo, concluyen que el modelo de acoplamiento disco-estrella ofrece un mecanismo físico necesario para explicar la existencia de estrellas Be masivas en binarias de órbita amplia donde el frenado por marea (tidal spin-down) es ineficiente.