The Mass Function of Neutron Stars from Core-Collapse Supernova Simulations

En lo profundo del corazón de las estrellas masivas, se desarrolla un drama final y violento. Cuando una estrella mucho más pesada que nuestro Sol se queda sin combustible, su núcleo colapsa bajo su propio peso, desencadenando una explosión de supernova que puede brillar más que una galaxia entera. Tras el desenlace, el núcleo a menudo sobrevive como una estrella de neutrones, un objeto tan denso que una sola cucharadita de su material pesaría tanto como una montaña. Durante décadas, los astrónomos han intentado comprender el nacimiento de estos remanentes cósmicos. Han medido las masas de las estrellas de neutrones encontradas en sistemas binarios, donde dos estrellas orbitan entre sí, pero las razones detrás de la distribución específica de estos pesos han permanecido esquivas. ¿Existe un peso estándar para una estrella de neutrones recién nacida, o nacen con una gran variedad de masas? Y, ¿cómo determina el caótico proceso de la muerte de una estrella si se convierte en una estrella de neutrones o colapsa completamente en un agujero negro?

Un nuevo estudio realizado por un equipo de investigadores reúne las simulaciones computacionales más avanzadas de estas explosiones estelares con las observaciones más recientes de estrellas de neutrones para responder a estas preguntas. Al mapear la estructura interna de la estrella moribunda con el peso del remanente que deja atrás, el equipo ha construido una imagen teórica de cómo nacen las estrellas de neutrones. Su trabajo sugiere que la variedad de masas que vemos hoy no es solo el resultado de que las estrellas ganen peso más tarde en su vida, sino que está escrita en el mismísimo momento de su creación. Los hallazgos ofrecen una explicación coherente de por qué la mayoría de las estrellas de neutrones se agrupan alrededor de un peso específico, mientras que un grupo más pequeño se extiende hacia masas mucho más pesadas, y proporcionan una nueva estimación de la frecuencia con la que las estrellas masivas no llegan a convertirse en estrellas de neutrones en absoluto, colapsando en su lugar en agujeros negros.

Los investigadores comenzaron analizando la intrincada relación entre el tamaño inicial de una estrella masiva y la estructura de su núcleo justo antes de explotar. Utilizaron modelos computacionales tridimensionales sofisticados que simulan la física compleja de una supernova, rastreando cómo se comporta el núcleo de la estrella en los momentos finales de su vida. Estas simulaciones revelaron un patrón claro: la "compacidad" interna del núcleo de la estrella —una medida de qué tan densamente empaquetada está la materia en un punto específico en su interior— determina directamente la masa de la estrella de neutrones que se forma. Un núcleo menos compacto tiende a dejar tras de sí una estrella de neutrones más ligera, mientras que un núcleo más compacto da como resultado una más pesada. Esta conexión permitió al equipo tomar la distribución conocida de las estrellas masivas en el universo y predecir la distribución de las masas de las estrellas de neutrones en el momento de su nacimiento.

Cuando compararon esta distribución de nacimiento teórica con las masas reales de las estrellas de neutrones observadas en el cielo, apareció una discrepancia. Las simulaciones predecían que las estrellas de neutrones recién nacidas más comunes deberían ser ligeramente más ligeras que la masa promedio observada, con un segundo grupo de estrellas más pesadas. Sin embargo, los datos observados mostraron un pico en un peso ligeramente superior. Los investigadores se dieron cuenta de que esta diferencia podía explicarse mediante dos factores principales: las estrellas de neutrones ganando masa después de nacer, y la forma en que nacen afectando si permanecen en un sistema binario.

Las estrellas de neutrones en sistemas binarios a menudo extraen material de sus estrellas compañeras, un proceso llamado acreción, que añade peso a la estrella de neutrones con el tiempo. El equipo aplicó correcciones a los datos observados para estimar cuáles habrían sido las masas en el momento del nacimiento, eliminando el peso ganado de sus compañeras. Este ajuste desplazó el pico observado hacia abajo para alinearlo mucho más estrechamente con la predicción teórica. Además, el estudio destacó el papel de los "patadas" (kicks) recibidas por las estrellas de neutrones al nacer. Así como un cohete es impulsado hacia adelante por su escape, una estrella de neutrones recibe un empuje repentino en la dirección opta de la explosión. Las simulaciones sugieren que las estrellas de neutrones más pesadas reciben patadas más fuertes en promedio. Estas poderosas patadas pueden romper el vínculo gravitacional que mantiene unido al sistema binario, expulsando a las estrellas de neutrones más pesadas al vacío donde son más difíciles de encontrar. Esta eliminación selectiva significa que los sistemas binarios que observamos están sesgados hacia las estrellas de neutrones más ligeras que recibieron patadas más suaves y se mantuvieron en órbita.

Al combinar los efectos de la formación de agujeros negros, la acreción de masa y estas patadas de nacimiento, el equipo creó un modelo refinado que coincide notablemente bien con los datos observados. Su modelo sugiere que el universo produce un amplio espectro de masas de estrellas de neutrones al nacer, que van desde aproximadamente 1.2 hasta 2.0 veces la masa de nuestro Sol. La masa de nacimiento más común es de alrededor de 1.35 veces la masa del Sol, lo que explica el pico prominente visto en las observaciones. La cola más pesada de la distribución, que se extiende hasta las 2.0 masas solares, es un resultado natural de las estructuras variables de las estrellas progenitoras, más que solo el resultado de la acreción posterior.

El estudio también proporcionó una nueva estimación de la fracción de estrellas masivas que terminan sus vidas como agujeros negros en lugar de estrellas de neutrones. Al analizar qué tipos de estrellas en sus simulaciones no lograron producir una estrella de neutrones, los investigadores calcularon que aproximadamente el 21 por ciento de las estrellas masivas en el rango de masa relevante colapsan directamente en agujeros negros. Este número no es una regla simple basada en el tamaño inicial de la estrella; más bien, depende de la estructura compleja y no lineal del núcleo de la estrella, lo que significa que estrellas de tamaños similares pueden tener destinos muy diferentes.

Este trabajo representa un paso significativo hacia la conexión de la teoría de cómo mueren las estrellas con la realidad de lo que observamos. Va más allá de la vieja idea de que todas las estrellas de neutrones nacen con casi la misma masa, reemplazándola con una imagen dinámica donde la estructura inicial de la estrella dicta una amplia gama de resultados. El acuerdo entre las complejas simulaciones y los datos observados sugiere que nuestra comprensión de las explosiones de supernova está madurando. Aunque todavía existen incertidumbres, como la influencia precisa de la composición química de la estrella o los detalles de las interacciones binarias, el estudio demuestra que la teoría moderna de las supernovas puede ahora explicar la población de objetos compactos que vemos en el cielo. Abre una nueva era en la que el nacimiento de estos objetos extremos puede rastrearse hasta las condiciones específicas de las estrellas que los crearon, convirtiendo una colección de mediciones en una historia coherente de muerte y renacimiento estelar.