On the impact of the carbon fusion rate over the properties of superbursts -- Numerical simulations of superbursts with MESA

Resumen Técnico: Sobre el impacto de la tasa de fusión de carbono en las propiedades de los superestallidos

Planteamiento del Problema
Los superestallidos son explosiones termonucleares energéticas que ocurren en las cortezas de estrellas de neutrones dentro de Binarias de Rayos X de Baja Masa (LMXB), desencadenadas por la combustión inestable de carbono a temperaturas $T \le 10^9$ K. Aunque las simulaciones multizona previas (p. ej., Keek & Heger 2011, 2012) han modelado con éxito los tiempos de recurrencia y las curvas de luz de los superestallidos, estas dependían de la tasa clásica de fusión de carbono de Caughlan & Fowler (1988, CF88). Investigaciones recientes en física nuclear han cuestionado esta tasa, sugiriendo que, a las temperaturas relevantes ($T \le 10^9$ K), la tasa real podría ser significativamente diferente. Específicamente, los datos experimentales de la colaboración LUNA (Tumino et al. 2018) sugieren que la tasa podría incrementarse por un factor de $\sim 10^3$, mientras que los modelos de impedimento (Monpribat et al. 2022) sugieren que podría reducirse por un factor similar. El impacto de estas variaciones sustanciales en los observables de los superestallidos —tales como los tiempos de recurrencia, los tiempos de decaimiento, las profundidades de ignición y los rendimientos de nucleosíntesis— permanece inexplorado en el contexto de los códigos modernos de evolución estelar.

Metodología
Los autores emplearon el código de evolución estelar público MESA v.24.08.1 para simular la evolución de las envolturas de estrellas de neutrones sometidas a acreción de masa. Las simulaciones utilizaron una red de reacciones nucleares reducida, Approx_A56 (que cubre nucleidos hasta $A=56$), con ejecuciones selectas utilizando Approx_A64 para verificar los efectos de truncamiento.

Los parámetros y configuraciones clave incluyeron:

• Tasas de Acreción: Se probaron dos tasas fiduciales: una "moderada" ($\dot{M}_{moderate} \approx 5.88 \times 10^{-9} M_\odot \text{yr}^{-1}$) y una "alta" ($\dot{M}_{high} \approx 2 \times 10^{-8} M_\odot \text{yr}^{-1}$).
• Composición: El material acrecionado se modeló como 20% de $^{12}\text{C}$ y 80% de $^{56}\text{Fe}$, consistente con los requisitos de la literatura para desencadenar superestallidos. La envoltura inicial era de $^{56}\text{Fe}$ puro.
• Calentamiento de la Base: Se aplicó una luminosidad de base fija ($L_b$) en la parte inferior de la envoltura para asegurar la ignición inestable del carbono, calibrada para coincidir con las simulaciones previas de KEPLER.
• Tasas de Fusión de Carbono: Se compararon cuatro tasas de reacción distintas:
  1. CF88: La tasa clásica (modelo fiducial).
  2. LUNA: Basada en las mediciones de LUNA, prediciendo una tasa $\sim 10^3$ veces mayor que CF88 a $T \sim 2\text{--}3 \times 10^8$ K.
  3. HIN: Basada en el modelo de impedimento (hindrance), prediciendo una tasa $\sim 10^3$ veces menor que CF88.
  4. HIN-RES: El modelo HIN incluyendo una resonancia de baja energía, produciendo resultados intermedios entre HIN y CF88.
• Tratamientos Físicos: Las simulaciones incluyeron opacidades radiativas (prescripción de Wolf), Teoría del Largo de Mezcla para la convección y condiciones de contorno de presión superficial fijas. Se omitieron las correcciones de relatividad general.

Resultados Clave
El estudio cuantifica cómo las variaciones en la tasa de fusión de carbono alteran las propiedades de los superestallidos:

1. Escalas de Tiempo y Profundidad de Ignición:

   • Tasa Incrementada (LUNA): Conduce a tiempos de recurrencia y de decaimiento más cortos, y desencadena la ignición a una profundidad de columna más superficial ($y_{ignition}$). Para $\dot{M}_{moderate}$, el tiempo de recurrencia cae de $\sim 241$ días (CF88) a $\sim 100$ días.
   • Tasa Reducida (HIN): Resulta en tiempos de recurrencia y de decaimiento más largos, con la ignición ocurriendo a profundidades de columna mayores. Para $\dot{M}_{moderate}$, el tiempo de recurrencia aumenta a $\sim 364$ días.
   • Tasa Intermedia (HIN-RES): Produce resultados casi indistinguibles del modelo fiducial CF88.

2. Temperatura y Curvas de Luz:

   • La temperatura máxima alcanzada durante la explosión es sensible a la tasa. La tasa HIN produce las temperaturas pico más altas ($\sim 5 \times 10^9$ K), mientras que LUNA produce picos más bajos ($\sim 3 \times 10^9$ K).
   • Los perfiles de decaimiento de la luminosidad se desplazan en consecuencia. Las simulaciones de LUNA muestran un decaimiento más rápido, mientras que las de HIN muestran un decaimiento más lento.
   • En $\dot{M}_{high}$, la luminosidad pico es mayormente insensible a la elección de la tasa, pero los tiempos de decaimiento e intervalos de recurrencia aún varían significamente.

3. Nucleosíntesis y Composición de las Cenizas:

   • La distribución de las cenizas exhibe una estructura característica de doble pico (alrededor de $^{28}\text{Si}$ y $^{56}\text{Fe}$) en todos los modelos.
   • Sensibilidad a la Tasa: Para $A \le 44$, la abundancia de nucleidos-$\alpha$ es altamente sensible a la tasa. La tasa LUNA (mayor temperatura/eficiencia) produce significativamente más nucleidos-$\alpha$ en comparación con CF88/HIN-RES, mientras que la tasa HIN produce menos.
   • Efectos de Alta Tasa de Acreción: En $\dot{M}_{high}$, hay un incremento notable en la síntesis de núcleos con $28 \le A \le 50$ y una acumulación específica de $\sim 0.1\%$ de $^{16}\text{O}$, una característica menos prominente en tasas moderadas.
   • Fotodesintegración: Las simulaciones reemplazando el $^{56}\text{Fe}$ inerte por $^{72}\text{Ge}$ confirmaron que la distribución de cenizas del pico del hierro está impulsada principalmente por la fotodesintegración del $^{56}\text{Fe}$ más que por la fusión de carbono por sí sola. Eliminar los canales de fotodesintegración detuvo la nucleosíntesis más allá de $^{40}\text{Ca}$.

4. Degeneración con el Calentamiento de la Base:

   • El estudio demuestra que cambiar la tasa de fusión de carbono tiene un efecto en las escalas de tiempo del superestallido y en la profundidad de ignición comparable a cambiar la luminosidad de base ($L_b$) de la envoltura. Por ejemplo, reducir la tasa (HIN) e incrementar $L_b$ produce tiempos de recurrencia similares. Esto sugiere una degeneración donde el ajuste de la curva de luz por sí solo no puede restringir de forma única tanto la tasa de fusión como el calentamiento de la base.

Significado y Conclusiones
El artículo concluye que la tasa de fusión de carbono a $T \le 10^9$ K es un parámetro crítico que gobierna la física de los superestallidos. Los autores afirman que:

• Las incertidumbres en la tasa de fusión de carbono pueden explicar las variaciones observadas en los tiempos de recurrencia y decaimiento sin necesidad de invocar cambios en la tasa de acreción de masa o en el calentamiento de la base.
• La composición "canónica" de las cenizas derivada de las simulaciones previas de KEPLER (usando CF88) es robusta respecto a la estructura de doble pico, pero sensible a la tasa de reacción específica en cuanto a las abundancias de nucleidos-$\alpha$.
• La degeneración entre la tasa de fusión de carbono y el calentamiento de la base implica que el ajuste de una única curva de luz de superestallido no puede determinar de forma inequívoca la física nuclear subyacente o el estado térmico de la corteza de la estrella de neutrones.
• Los resultados validan el uso de MESA para el modelado de superestallidos, mostrando consistencia con los resultados previos de KEPLER al tiempo que ofrece un código público alternativo para explorar las incertidumbres de la física nuclear.

Los autores señalan modestamente que, aunque sus hallazgos sugieren que la tasa de fusión de carbono podría aliviar la necesidad de un alto calentamiento de la base en los modelos (si la tasa LUNA es la correcta), la incertidumbre en las tasas de acreción de masa para fuentes específicas (como 4U 1636-53) impide conclusiones definitivas sobre cuál tasa es la físicamente realizada. Además, aclaran que su tratamiento de la combustión de carbono como puramente termonuclear está justificado por las densidades y temperaturas alcanzadas, siendo distinto de los regímenes piocore nucleares.