Are NICER and GW170817 constraints suggesting a compactified scenario for Neutron stars?

Resumo Técnico: As Restrições de NICER e GW170817 Sugerem um Cenário de Compactação para Estrelas de Nêutrons?

Enunciado do Problema
A equação de estado (EoS) da matéria em densidades que excedem a saturação nuclear permanece um problema fundamental não resolvido na física. Enquanto a teoria de campo efetivo quiral ($\chi$EFT) modela de forma confiável densidades subnucleares e a QCD perturbativa (pQCD) descreve regimes de alta densidade assintóticos, a região de densidade intermediária nos núcleos de estrelas de nêutrons (NS) é mal compreendida. Observações astrofísicas recentes do Neutron Star Interior Composition Explorer (NICER) e do evento de ondas gravitacionais GW170817 forneceram restrições precisas sobre massas, raios e deformabilidades de maré de NS. No entanto, permanece incerto se essas restrições favorecem uma classe específica de EoS (monotônica, não monotônica ou descontínua) ou se revelam um padrão estrutural universal na matéria densa. Especificamente, o artigo investiga se os dados atuais sugerem um cenário de "compactação" onde a gravidade naturalmente favorece configurações que maximizam a compactação estelar em vez de apenas maximizar a massa.

Metodologia
Os autores empregam uma estrutura abrangente de inferência Bayesiana para analisar três conjuntos distintos e agnósticos de classe de EoS:

1. Monotônica: Caracterizada por uma velocidade do som ($c_s^2$) que aumenta monotonicamente com a densidade (típica de EoS hadrônicas).
2. Não monotônica: Caracterizada por um pico em $c_s^2$ seguido por uma diminuição (associada a transições de fase suaves).
3. Descontínua: Caracterizada por um salto de densidade ($\Delta n$) em uma transição de fase (associada a transições de fase de primeira ordem fortes).

Construção da EoS:
A EoS é construída usando uma abordagem híbrida. A matéria de baixa densidade ($< 0.5 n_0$) é definida pela EoS tabulada BPS. A região de $0.5 n_0$ a $1.1 n_0$ é modelada via uma politrope com um índice adiabático $\Gamma$ restrito pelas bandas de incerteza da $\chi$EFT. Para densidades acima de $1.1 n_0$, o quadrado da velocidade do som $c_s^2$ é parametrizado como uma função linear por partes do potencial químico $\mu$ através de quatro segmentos randomizados, estendendo-se até o regime pQCD ($\mu \approx 2.6$ GeV).

Estrutura de Inferência:
Diferente de métodos anteriores de amostragem de rejeição que geram milhões de amostras de EoS e as filtram, este estudo utiliza o algoritmo de amostragem aninhada UltraNest para explorar diretamente a distribuição posterior. O modelo contém 10 parâmetros livres (potenciais químicos e valores do quadrado da velocidade do som para os segmentos, além de $\Gamma_{CET}$ e um parâmetro de salto de densidade $\delta n$ para a classe descontínua).
A análise incorpora:

• Dados do NICER: Distribuições posteriores conjuntas para massa e raio de PSR J0030+0451, PSR J0740+6620, PSR J0437−4715 e o recentemente relatado PSR J0614−3329.
• Dados do GW170817: Amostras posteriores de massas de componentes e deformabilidades de maré da fusão de estrelas de nêutrons binárias.
• Restrições de $\chi$EFT: Verossimilhanças gaussianas para pressão em baixas densidades.

O estudo compara as distribuições posteriores dos parâmetros da EoS e das propriedades estelares derivadas (sequências Massa-Raio) sob diferentes combinações de restrições observacionais.

Principais Resultados

1. Poder de Restrição: A análise Bayesiana confirma que os dados observacionais restringem efetivamente a região de baixa densidade da EoS (especificamente $\Gamma_{CET}$ e o primeiro segmento da velocidade do som $c_{s,1}^2$), que permanece robusto em todas as classes de EoS. Em contraste, os parâmetros de alta densidade ($c_{s,2}^2, c_{s,3}^2$) exibem distribuições posteriores mais amplas, indicando incerteza significativa remanescente.
2. Impacto de Novos Dados: A inclusão de PSR J0614−3329 e GW170817 desloca sistematicamente os parâmetros do potencial químico para valores mais altos e reduz os limites de massa máxima. No entanto, a evidência Bayesiana ($\ln Z$) não favorece fortemente uma classe de EoS sobre outra; as diferenças na evidência caem na faixa de "inconclusiva a fraca".
3. Massa Máxima vs. Compactação Máxima: Uma distinção crítica é traçada entre a EoS que produz a massa máxima possível ($M_{TOV}$) e a Eosa que produz a compactação máxima ($C = GM/Rc^2$).
   • A EoS de Massa Máxima é rígida desde densidades muito baixas, fornecendo suporte de pressão uniforme para sustentar as maiores massas. Essas configurações ocupam a borda de grande raio da posterior M-R.
   • A EoS de Compactação Máxima apresenta uma EoS de baixa-a-média densidade mais suave (minimizando o raio) seguida por um endurecimento pronunciado em altas densidades (suportando a massa).
4. Preferência Observacional: A EoS mais favorecida pelos dados combinados de NICER e GW170817 (maior $\ln Z$) espelha estruturalmente a EoS de Compactação Máxima, não a EoS de Massa Máxima. A EoS observacionalmente favorecida exibe relativa suavidade em densidades intermediárias, transitando para rigidez em altas densidades.
5. Transições de Fase: Para a classe descontínua, os dados preferem pequenos saltos de densidade ($\delta n \approx 0.09-0.10$ fm$^{-3}$) e transições de fase tardias (ocorrendo em densidades mais altas). Embora uma transição de fase de baixa densidade pudesse teoricamente produzir uma estrela altamente compacta com um grande núcleo de quarks, as previsões Bayesianas favorecem configurações onde a transição de fase ocorre mais tarde, confinando a matéria de quarks aos núcleos de estrelas suficientemente massivas.

Significância e Alegações
O artigo alega que as restrições astrofísicas atuais estão selecionando configurações de EoS que produzem estrelas compactas em vez de estrelas maximamente massivas. A convergência da EoS observacionalmente favorecida com a sequência de compactação máxima sugere uma propriedade fundamental da matéria densa: a EoS provavelmente sofre um endurecimento significativo em densidades intermediárias (potencialmente sinalizando o início de novos graus de liberdade) antes de tornar-se rígida novamente para suportar as estrelas de $\sim 2 M_\odot$ observadas.

Os autores concluem que esta estrutura "suave-para-rígida" é o que a gravidade naturalmente favorece para alcançar a compactação. Embora os dados não descartem definitivamente classes específicas de EoS, eles apontam para um padrão universal onde a equação de estado transita da suavidade para a rigidez, conduzindo as estrelas de nêutrons para seus limites de compactação. O artigo nota modestamente que novas restrições sobre pulsares de baixa massa (para definir melhor a EoS de baixa densidade) e pulsares de alta massa (para restringir a EoS de alta densidade) serão necessárias para refinar ainda mais estas conclusões e potencialmente descartar cenários específicos, como transições de fase precoces.