Maximal mass of neutron stars constrained by neutron star observations

Mediante el empleo de un marco de ponderación bayesiano que integra observaciones multimensajero (incluyendo GW170817, NICER y candidatos a objetos compactos) con ecuaciones de estado híbridas, este estudio determina que la masa máxima de una estrella de neutrones está robustamente restringida a aproximadamente 2.2–2.3 masas solares, mientras que el radio correspondiente depende más fuertemente del modelo hadrónico subyacente, cayendo típicamente cerca de 12 km.

Lo esencial

Imagina que el universo está lleno de un material misterioso y superdenso llamado "materia de estrella de neutrones". Es tan pesado que una sola cucharadita pesaría tanto como una montaña. Durante mucho tiempo, los físicos han intentado descifrar las "reglas del juego" de este material, específicamente, qué tan rígido o compresible es. Este conjunto de reglas se denomina Ecuación de Estado (EOS).

La gran pregunta que este artículo intenta responder es: ¿Cuál es el peso absoluto máximo que puede alcanzar una estrella de neutrones antes de colapsar en un agujero negro?

Aquí está la historia de cómo los autores resolvieron este acertijo, explicada en términos sencillos:

1. Los Dos Puntos de Partida (Las Recetas)

Para descifrar las reglas, los científicos comenzaron con dos "recetas" diferentes sobre cómo se comporta esta materia densa a densidades más bajas. Piensa en estas como dos teorías distintas sobre cómo se mezclan los ingredientes:

• Receta A (SFHo): Una receta "más blanda", lo que significa que la materia es un poco más fácil de comprimir.
• Receta B (DD2): Una receta "más rígida", lo que significa que la materia resiste más la compresión.

Sabían que estas recetas funcionaban bien al "inicio" de la escala de densidades, pero no sabían qué ocurría en las densidades extremas y superaltas que se encuentran en el centro de una estrella de neutrones. Para llenar ese vacío, utilizaron un "puente" matemático para conectar sus recetas con lo que sabemos sobre la física de partículas en las energías más altas posibles.

2. El Trabajo de Detective (Usando Pistas Reales)

En lugar de simplemente adivinar, los autores actuaron como detectives. Tomaron sus dos recetas y las pusieron a prueba contra pistas del mundo real recopiladas por telescopios y detectores de ondas gravitacionales. Utilizaron un método estadístico especial (llamado ponderación bayesiana) para ver qué versiones de sus recetas sobrevivieron a la prueba.

Estas son las pistas que utilizaron:

• El "Gran Choque" (GW170817): Cuando dos estrellas de neutrones chocaron entre sí, enviaron ondas a través del espacio. La forma en que se comportaron estas ondas les dijo a los científicos qué tan "compresibles" eran las estrellas.
• La "Linterna" (NICER): Un telescopio espacial tomó fotografías de puntos calientes en estrellas de neutrones en rotación. Al medir qué tan grandes parecían las estrellas y qué tan pesadas eran, obtuvieron una relación directa entre tamaño y peso.
• El Candidato "Ligero" (HESS J1731–347): Un objeto muy pequeño y ligero que podría ser una estrella de neutrones.
• El Candidato "Pesado" (GW190814): Un objeto misterioso que es más pesado que la mayoría de las estrellas de neutrones pero más ligero que la mayoría de los agujeros negros. Los científicos se preguntaron: ¿Podría esto ser en realidad una estrella de neutrones superpesada?

3. Los Resultados: Lo Que las Pistas les Dijeron

Los científicos ejecutaron sus dos recetas a través de estas pistas y examinaron los resultados.

El Límite de Peso (Masa Máxima):

• La Sorpresa: No importó mucho qué receta inicial (Blanda o Rígida) utilizaron. Las pistas del mundo real fueron tan fuertes que obligaron a ambas recetas a coincidir en la misma respuesta.
• El Veredicto: Cuando utilizaron las pistas más fiables (el "Gran Choque" y la "Linterna"), el peso máximo que una estrella de neutrones puede soportar es de aproximadamente 2.2 a 2.3 veces la masa de nuestro Sol.
• El Giro del "Pesado": Si asumen que ese objeto pesado misterioso (GW190814) es una estrella de neutrones, el límite salta a aproximadamente 2.6 a 2.7 veces la masa del Sol. Sin embargo, esto crea un conflicto con las pistas de "compresibilidad" del Gran Choque, lo que hace que la situación sea complicada.

El Límite de Tamaño (Radio):

• La Diferencia: A diferencia del peso, el tamaño de la estrella sí dependía de qué receta inicial utilizaron.
• El Veredicto: La receta "Blanda" predijo un radio de aproximadamente 11.8 km, mientras que la receta "Rígida" predijo aproximadamente 12.4 km.
• El Punto Dulce: Cuando se combinan todas las mejores pistas, el tamaño más probable para estas estrellas es de alrededor de 12 kilómetros (más o menos 1 km).

4. El Panorama General

El artículo concluye que al observar los "puntos finales" (las estrellas más pesadas y grandes posibles) y utilizando una mezcla de datos astronómicos reales, podemos acotar las reglas de la materia más densa del universo.

• El Peso: El universo parece tener un "límite de velocidad" para lo pesada que puede ser una estrella de neutrones, situándose cómodamente alrededor de 2.2 a 2.3 masas solares. Esto coincide con la estrella de neutrones más pesada que hemos visto realmente hasta ahora.
• El Tamaño: Son aproximadamente del tamaño de una ciudad pequeña, con unos 12 km de diámetro.
• La Conclusión: Las observaciones del mundo real (las pistas) son mucho más poderosas que las suposiciones teóricas iniciales. No importa con qué teoría comiences, los datos de las propias estrellas fuerzan la respuesta a converger hacia los mismos números.

En resumen, el universo nos ha dado una respuesta muy clara: Las estrellas de neutrones pueden volverse increíblemente pesadas, pero hay un techo duro, y son sorprendentemente pequeñas para lo que pesan.

Resumen técnico

Aquí se presenta un resumen técnico detallado del artículo "Masa máxima de estrellas de neutrones restringida por observaciones de estrellas de neutrones" de Gábor Kasza y György Wolf.

1. Planteamiento del Problema

El desafío central en la física nuclear y de iones pesados es determinar la Ecuación de Estado (EOS) de la materia fuertemente interactiva a densidades y potenciales químicos intermedios. Mientras que la materia a baja densidad está restringida por experimentos nucleares y QCD en retículo, y la materia a alta densidad por QCD perturbativa (pQCD), la región relevante para los núcleos de estrellas de neutrones (varias veces la densidad de saturación nuclear, $\rho_0$) permanece poco comprendida.

El problema específico abordado es determinar la masa máxima ($M_{TOV}$) y el radio correspondiente ($R_{TOV}$) de las estrellas de neutrones. Estos puntos finales de las secuencias masa-radio ($M(R)$) son críticos para comprender la rigidez de la EOS a densidades supranucleares. Los autores buscan cuantificar cómo las observaciones multimensaje actuales restringen estos puntos finales y en qué medida los resultados dependen de la elección del modelo base hadrónico subyacente.

2. Metodología

A. Construcción de EOS Híbridas

Los autores construyen familias de EOS híbridas causales conectando materia hadrónica a baja densidad con materia de quarks a alta densidad:

• Régimen de Baja Densidad: Se utilizan dos modelos hadrónicos representativos como bases:
  • SFHo: Un modelo de campo medio relativista (EOS más blanda, incompresibilidad $K=245$ MeV).
  • DD2: Una EOS hadrónica más rígida.
• Régimen de Alta Densidad: Un Modelo Sigma Lineal extendido (eLSM) basado en la simetría quiral $U(3) \times U(3)$, incorporando quarks constituyentes y nonetos de mesones. Este modelo está restringido para coincidir con la termodinámica de QCD en retículo a bajo potencial químico y con los resultados de pQCD a alto potencial químico.
• Interpolación: Para unir ambos regímenes, se aplica una interpolación polinómica de quinto orden a la densidad de energía en función de la densidad bariónica. Esto asegura la consistencia termodinámica y la continuidad de la presión y la velocidad del sonido ($c_s$) a través de la región de transición, evitando transiciones de fase de primer orden grandes.
• Restricción Asintótica: Se requiere que la EOS de mesones-quarks se conecte suavemente con los resultados de pQCD a densidades asintóticamente altas ($\mu \approx 2.6$ GeV).

B. Marco Bayesiano

El estudio emplea un marco de ponderación bayesiana para analizar las distribuciones de probabilidad de los puntos finales de la secuencia $M(R)$.

1. Prior: Se genera un amplio conjunto de EOS híbridas causales que abarcan un rango de rigidez a densidades supranucleares.
2. Ponderación de Verosimilitud: El conjunto se pondera frente a datos observacionales para generar distribuciones de probabilidad posteriores para $M_{TOV}$ y $R_{TOV}$.

C. Restricciones Observacionales

Las siguientes restricciones se aplican secuencialmente y en combinación:

1. Restricción pQCD: Asegura que la EOS se aproxime al punto de referencia N3LO de pQCD ($\mu=2.6$ GeV, $\rho=6.47$ fm$^{-3}$, $p=3823$ MeV/fm$^3$) manteniendo la causalidad ($c_s < 1$).
2. GW170817: Restricción de deformabilidad de marea ($\tilde{\Lambda} < 720$) y una restricción combinada más estricta ("$\Lambda_{190}$": $70 < \Lambda_{1.4} < 580$ y $9.1 < R_{1.4} < 12.8$ km).
3. NICER: Mediciones de masa-radio para cinco púlsares (J0030+0451, J0740+6620, J0614–3329, J1231–1411, J0437–4715) con una precisión de $\sim 10\%$.
4. HESS J1731–347: Un candidato para una estrella de neutrones de masa muy baja ($M \approx 0.77 M_\odot$), tratado con cautela debido a incertidumbres en la interpretación.
5. GW190814: El objeto de "brecha de masa" ($M \approx 2.59 M_\odot$). Los autores prueban la hipótesis de que este objeto es una estrella de neutrones aplicando una restricción gaussiana centrada en $2.59 M_\odot$.

3. Contribuciones Clave

• Diagnóstico Unificado: Demuestra que las distribuciones de probabilidad de los puntos finales de la secuencia $M(R)$ sirven como un diagnóstico sensible y complementario para restringir el comportamiento de la EOS a alta densidad dentro de un marco multimensaje.
• Independencia vs. Dependencia de la Base: Desentraña sistemáticamente la influencia de los datos observacionales de la elección teórica de la base hadrónica (SFHo vs. DD2).
• Identificación de Tensión: Identifica una tensión específica entre la EOS rígida DD2, la hipótesis de la estrella de neutrones de "brecha de masa" y las restricciones estrictas de deformabilidad de marea ($\Lambda_{190}$).

4. Resultados Clave

Distribuciones de Masa Máxima ($M_{TOV}$)

• Dominio de las Observaciones: El pico de la distribución de $M_{TOV}$ está impulsado principalmente por restricciones observacionales más que por la elección de la EOS hadrónica base (SFHo o DD2).
• Progresión de Restricciones:
  • Solo pQCD: Produce un pico ligeramente superior a $1.5 M_\odot$.
  • GW170817/HESS: No hay un desplazamiento significativo en la posición del pico.
  • NICER: Desplaza el pico a 2.1–2.2 $M_\odot$ y afila la distribución.
  • Brecha de Masa (GW190814): Desplaza el pico significativamente a 2.6–2.7 $M_\odot$.
• Conclusión Robusta: Al aplicar las restricciones más robustas (pQCD + GW170817 + NICER), ambos modelos de EOS convergen en una masa máxima preferida de 2.2–2.3 $M_\odot$. Esto es consistente con la masa observada de PSR J0952–0607 ($2.35 \pm 0.17 M_\odot$).

Distribuciones de Radio ($R_{TOV}$)

• Dependencia de la Base: A diferencia de la masa, la distribución del radio muestra una fuerte dependencia de la EOS hadrónica subyacente.
  • SFHo+eLSM: Pico cerca de 11.8 km (con restricciones robustas).
  • DD2+eLSM: Pico cerca de 12.4 km (con restricciones robustas).
• Preferencia General: Las restricciones combinadas favorecen un rango de radio de $12 \pm 1$ km.
• Impacto de la Deformabilidad de Marea: Aplicar la restricción más estricta $\Lambda_{190}$ estrecha las distribuciones de radio y reduce los valores preferidos, ya que una menor deformabilidad de marea se correlaciona con radios más pequeños.

La Tensión de la "Brecha de Masa"

• Si el objeto GW190814 se interpreta como una estrella de neutrones, la EOS basada en DD2 no logra proporcionar una descripción consistente cuando se combina con la restricción de marea $\Lambda_{190}$. La rigidez requerida para soportar una estrella de $2.6 M_\odot$ en el modelo DD2 resulta en deformabilidades de marea demasiado grandes para satisfacer los límites de GW170817.
• La EOS basada en SFHo (más blanda) maneja esta combinación ligeramente mejor, pero aún enfrenta desafíos.

5. Significado

Este trabajo proporciona una evaluación estadística rigurosa de la masa máxima de las estrellas de neutrones, avanzando más allá de estimaciones de punto único hacia distribuciones de probabilidad.

1. Validación de Modelos Actuales: Los resultados apoyan la existencia de estrellas de neutrones hasta $\sim 2.2-2.3 M_\odot$, validando los modelos actuales de física nuclear frente a datos astrofísicos.
2. Restricción sobre Transiciones de Fase: El procedimiento de ajuste suave y las distribuciones resultantes sugieren que es poco probable que ocurran transiciones de fase de primer orden grandes en el rango de densidades sondeado por las observaciones actuales, o al menos están estrictamente restringidas.
3. Futuras Direcciones: El estudio destaca que la precisión futura en las mediciones de radio (por ejemplo, de NICER o futuras misiones de rayos X) y la naturaleza definitiva del objeto GW190814 (estrella de neutrones vs. agujero negro) son críticas para resolver las incertidumbres restantes en la EOS de alta densidad. La tensión entre la hipótesis de la brecha de masa y las restricciones de deformabilidad de marea sugiere que, si GW190814 es efectivamente una estrella de neutrones, la EOS debe ser más blanda de lo que predice el modelo DD2, o las restricciones de marea requieren una reevaluación.