Qué revelan los isótopos de calcio y estaño sobre la materia invisible en las estrellas de neutrones

Una estrella tan densa que ningún laboratorio puede reproducir

Una estrella de neutrones concentra más masa que el Sol en una esfera de apenas decenas de kilómetros de diámetro. Esa densidad supera con creces lo que cualquier laboratorio puede generar, por lo que los astrónomos la deducen a partir de las propias estrellas. Una de ellas, el pulsar J0740+6620, el más masivo conocido hasta la fecha, aporta ahora un valor de 2,08 ± 0,07 masas solares a dicha estimación. Dado que una estrella de ese peso no puede colapsar en un agujero negro, la materia más densa en su interior debe ser lo suficientemente rígida para soportar la carga. Otras mediciones, como cuatro determinaciones del radio con NICER y una señal de fusión, impiden que esa rigidez sea arbitraria.

Una calibración bayesiana convierte esas observaciones en distribuciones de probabilidad, en lugar de ofrecer una única respuesta fija sobre la materia nuclear. La masa del pulsar actúa como un límite superior a lo blando que puede ser la región de alta densidad, mientras que los radios y la deformación por fusión sirven de control para el espacio que dicha estrella puede ocupar. Juntas, estas restricciones reducen la familia de modelos aceptables, aunque no la cierran en una sola curva. Cada modelo superviviente conserva su incertidumbre cuantificada, de modo que el resultado es un rango de comportamientos posibles para estrellas densas con probabilidades explícitas.

Esa forma probabilística importa porque el mismo cálculo debe mantenerse conectado a la materia nuclear ordinaria cerca de la densidad de saturación. Una suposición libre sobre la alta densidad podría sostener estrellas masivas mientras contradice lo que sabemos de los núcleos en el laboratorio. El protocolo, por tanto, transmite la información de baja densidad hacia adelante, para que el núcleo denso y las capas externas formen parte de una única descripción física. De este modo, el pulsar más masivo no añade simplemente un número: ancla una región de materia inaccesible de otro modo a evidencia que podemos medir realmente.

Los núcleos de calcio y estaño con capas abiertas ajustan el modelo nuclear

La calibración nuclear basada en probabilidades incorpora ahora núcleos de calcio y estaño con capas abiertas para evaluar cómo maneja las partículas no emparejadas. Estos núcleos son clave porque sus masas y tamaños medidos dependen del apareamiento y de la variación de la energía nuclear con la densidad. El ajuste añade cuatro isótopos de calcio (42, 44, 46 y 50) y tres de estaño (112, 120 y 124). Sus energías de enlace y radios de carga proporcionan al modelo un conjunto más amplio de puntos de referencia a baja densidad que los núcleos de capas cerradas por sí solos. La señal de apareamiento de neutrones del estaño-120 ayuda, además, a fijar la fuerza de apareamiento añadida en el cálculo. Este término es importante porque los núcleos con capas abiertas pueden tener neutrones sin emparejar, y sus masas revelan con qué intensidad se unen. La inclusión de estos isótopos ha llevado la energía de enlace a densidad nuclear ordinaria a valores más bajos. También ha desplazado ligeramente hacia arriba los parámetros de la energía de simetría. Esta energía representa el coste adicional de tener más neutrones que protones en un núcleo. Su pendiente describe cómo aumenta ese penalización por exceso de neutrones al aumentar la densidad, y la pendiente inferida se mantiene baja frente a otras estimaciones. No obstante, los valores superiores a 50 MeV ya no quedan fuertemente excluidos una vez incluidos los datos de capas abiertas en el ajuste. El resultado es, por tanto, una distribución de probabilidad sobre posibles parámetros nucleares, no una respuesta única para la materia densa. Esa dispersión importa porque las energías de enlace tiran hacia una región mientras que los observables ricos en neutrones lo hacen desde otra. Así, las mediciones de calcio y estaño restringen algunas partes del modelo sin borrar la tensión entre ellas.

La corteza gana iones más pesados sin que cambie su resistencia mecánica

Una inferencia estadística actualizada, basada en mediciones de núcleos atómicos, desplaza el rango de probabilidad de la carga iónica en la corteza de las estrellas de neutrones hacia valores más altos. Este cambio refleja un mayor valor del término de energía de simetría a baja densidad, que cuantifica el coste energético de una mezcla desigual de neutrones y protones. Por tanto, los iones de la corteza tienden a cargas más pesadas no porque la red cristalina se reconstruya desde cero, sino porque aumenta la fracción de protones permitida. La misma inferencia deja prácticamente inalteradas las propiedades mecánicas globales de la corteza, incluida la presión y el módulo de cizalladura.

El tamaño de la celda que determina cómo se sitúan esos iones en el sólido también permanece cerca de su valor anterior. La densidad central, la densidad del gas y el radio no muestran cambios significativos al añadir los nuevos parámetros nucleares. Como los iones siguen incrustados en una red con casi el mismo espaciado, la resistencia de la corteza a la cizalladura no aumenta junto con la mayor carga. Incluso las correcciones por tamaño finito aplicadas a los iones dejan ese módulo de cizalladura esencialmente donde estaba.

El efecto es, por tanto, local en cuanto a composición y no global en cuanto a estructura. La ecuación de estado y la relación masa-radio permanecen cercanas al análisis anterior de estrellas de neutrones y siguen ajustarse a las mediciones de radio de NICER. Esa separación importa, porque una corteza puede albergar iones diferentes sin obligar a que el tamaño total o el perfil de presión de la estrella cambien.

Energía de simetría más blanda a baja densidad y comportamiento más rígido a alta densidad

La calibración bayesiana del modelo Skyrme arroja una distribución de probabilidad para la pendiente de la energía de simetría que resulta relativamente baja cerca de la densidad nuclear ordinaria. Dicho modelo es, en esencia, una receta matemática para calcular la energía nuclear. Esta calibración se apoya en datos de núcleos de calcio y estaño con capas abiertas, otros observables de núcleos finitos, cálculos directos de materia neutra pura, mediciones del tamaño mediante rayos X realizadas por NICER y los límites mareales derivados de las ondas gravitacionales. Su valor reside en que las correlaciones se conservan a lo largo del proceso, sin reiniciarse en cada paso. Gracias al resumen estadístico compacto, otros grupos pueden emplear esas incertidumbres vinculadas directamente, en lugar de suponer que cada parámetro nuclear varía de forma independiente.

Esta inferencia transforma un problema complejo de múltiples parámetros en un punto de partida común. Los isótopos de calcio y estaño con capas abiertas aportan energías de enlace, radios de carga y una restricción sobre el hueco de emparejamiento, lo que impide que el resultado dependa de una sola familia nuclear o de una masa mal descrita. A continuación, la etapa correspondiente a las estrellas de neutrones vincula esos cambios nucleares con la corteza y el núcleo estelar, manteniendo al mismo tiempo la curva global de masa-radio compatible con las mediciones del radio realizadas por NICER.

La limitación principal es que los datos de núcleos finitos restringen el modelo con mayor firmeza cerca de la densidad nuclear ordinaria, mientras que las densidades del núcleo se abordan mediante una extensión flexible a alta densidad. Dentro de esta calibración, el comportamiento permitido a alta densidad tiende a rigidizarse por encima de la densidad nuclear ordinaria si el modelo ha de seguir siendo compatible con las mediciones del radio de NICER y con una masa de púlsar de unas 2,08 masas solares, utilizada como restricción de masa máxima. Los límites mareales de las ondas gravitacionales acotan entonces hasta qué punto puede rigidizarse esa región de alta densidad. Además, el resumen gaussiano compacto no captura la asimetría ni la multimodalidad, especialmente en el caso de un parámetro superficial cuya distribución de probabilidad es amplia y de forma irregular. Por ello, otros grupos que utilicen dicho resumen deben tener presente que no representa la totalidad de la incertidumbre.

Dado que la inferencia sigue siendo probabilística, las futuras mediciones del radio y los efectos mareales pondrán a prueba cuánto debe rigidizarse la materia nuclear a alta densidad mientras se mantiene relativamente blanda la energía de simetría a baja densidad.